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]]))
