def calc_sig_b(b, gaussmeter_file): sig_gauss = np.genfromtxt(gaussmeter_file) sig_b = np.zeros(len(b)) for j in range(len(b)): if (abs(b[j]) < sig_gauss[0, 0]): sig_b[j] = quad_add([sig_gauss[0, 1], sig_gauss[0, 2]]) if (abs(b[j]) > sig_gauss[0, 0] and abs(b[j]) < sig_gauss[1, 0]): sig_b[j] = quad_add([sig_gauss[1, 1], sig_gauss[1, 2]]) if (abs(b[j]) > sig_gauss[1, 0]): sig_b[j] = quad_add([sig_gauss[2, 1], sig_gauss[2, 2]]) return sig_b
def calc_sig_b(b,gaussmeter_file): sig_gauss = np.genfromtxt(gaussmeter_file) sig_b = np.zeros(len(b)) for j in range(len(b)): if (abs(b[j]) < sig_gauss[0,0]): sig_b[j] = quad_add([sig_gauss[0,1],sig_gauss[0,2]]) if (abs(b[j]) > sig_gauss[0,0] and abs(b[j]) < sig_gauss[1,0]): sig_b[j] = quad_add([sig_gauss[1,1],sig_gauss[1,2]]) if (abs(b[j]) > sig_gauss[1,0]): sig_b[j] = quad_add([sig_gauss[2,1],sig_gauss[2,2]]) return sig_b
def calc_sig_b_dev(b, gaussmeter_file): sig_gauss = np.genfromtxt(gaussmeter_file) sig_b = np.zeros(len(b)) for j in range(len(b)): for k in range(len(sig_gauss[:, 0])): if (abs(b[j]) < sig_gauss[k, 0]): sig_b[j] = quad_add([sig_gauss[k, 1], sig_gauss[k, 2]]) break return sig_b
def calc_sig_b_dev(b,gaussmeter_file): sig_gauss = np.genfromtxt(gaussmeter_file) sig_b = np.zeros(len(b)) for j in range(len(b)): for k in range(len(sig_gauss[:,0])): if (abs(b[j]) < sig_gauss[k,0]): sig_b[j] = quad_add([sig_gauss[k,1], sig_gauss[k,2]]) break return sig_b
def i_stats(i, multimeter_file): # sig_mult = np.genfromtxt(multimeter_file) sig_mult=[0, 0, 0] return (np.mean(i), quad_add([np.std(i), sig_mult[2]]))
def i_stats(i, multimeter_file): # sig_mult = np.genfromtxt(multimeter_file) sig_mult = [0, 0, 0] return (np.mean(i), quad_add([np.std(i), sig_mult[2]]))