# The figure plots column 1 and 2 

# A. critical \alpha:

                        #a_exp  a_crit
    H_alpha=  np.array([[0.474, 0.44],
                        [0.412, 0.50],
                        [0.416, 0.64],
                        [1.28 , 1.27]])

    M_alpha=  np.array([[0.392, 0.52],
                        [1.02 , 1.08],
                        [1.76 , 1.73]])

    L_alpha = np.array([[0.300, 0.43],
                        [0.387, 0.52],
                        [0.543, 0.60],
                        [0.486, 0.70],
                        [1.150, 1.30]])

    U_alpha = np.array([[0.325, 0.28],
                        [0.529, 0.40],
                        [0.695, 0.59],
                        [0.623, 0.61],
                        [0.883, 0.84],
                        [1.340, 1.31]])


# B. critical a/LTi:
# Here the last column (3) indicates whether the value
# was stopped at 4 times the nominal value 
# (0 = stopped at 4times, 1 = actual GK threshold)

                        #a_exp a_crit open_star(=0)
    M_RLTi =  np.array([[1.85, 4.2,   1],
                        [2.37, 5.7,   1],
                        [2.07, 8.3,   0],
                        [2.25, 9.0,   0]])

    L_RLTi =  np.array([[0.79, 2.40,  1],
                        [1.53, 2.70,  1],
                        [4.59, 5.50,  1],
                        [5.22, 10.4,  1],
                        [4.55, 18.0,  0]])

    U_RLTi =  np.array([[1.28, 3.10,  1],
                        [2.29, 2.75,  1],
                        [4.62, 5.56,  1],
                        [7.40, 30.0,  0],
                        [9.09, 36.0,  0]])