r0 = 2.818e-15 alpha = 662 * 1.602e-19 / (9.109e-31 * 3e8 ** 2) theta = arange(0, 161) * pi / 180 cross_section = ( r0 ** 2 * (1 / (1 + alpha * (1 - cos(theta)))) ** 2 * (1 + cos(theta) ** 2) / 2 * (1 + (alpha ** 2 * (1 - cos(theta)) ** 2) / ((1 + cos(theta) ** 2) * (1 + alpha * (1 - cos(theta))))) ) # plot(theta, cross_section) # show() R = 190e-3 l = log(2) / (30.17 * yrs(1, 0, 0)) activity = 122.1e6 * exp(-l * (yrs(2013, 1, 2) - yrs(1979, 3, 22))) print "activity is %d" % activity flux = activity / (4 * pi * R ** 2) print "flux is %d" % flux fileset = [ ("AX20.ASC", (3350, 3850)), ("AX40.ASC", (3350, 3800)), ("AX60.ASC", (3400, 3750)), ("AX80.ASC", (3300, 3750)), ("AX90.ASC", (3300, 3900)), ("AX140.ASC", (3300, 3750)), ] dataset = []
rcParams['text.usetex'] = True rcParams['text.latex.unicode'] = True rcParams['font.family'] = 'serif' datapath = "data/effectivity/" imgpath = "img/" # For activities we have # n(t) = n_0 exp(- l t) <=> # n(0) = n_0 => # n(t) = n(0) exp(- l (t-0)) # DO NOT CHANGE THE DATES!!! gamma_R = 0.155 # Distance between preparate (on electron detector) # and gamma detector l_cs = log(2) / (30.17 * yrs(1, 1, 0)) l_am = log(2) / (432.2 * yrs(1, 1, 0)) cesium_activity = 404e3 * exp( -l_cs * (yrs(2013, 1, 5) - yrs(1979, 12, 01))) americium_activity = 401e3 * exp( -l_am * (yrs(2013, 1, 5) - yrs(1986, 02, 01))) gamma_area = pi * (50.8e-3) ** 2 / 4 data = [ ("ABGG", mins(1, 13)), ("CC", mins(15, 26)), ("CA", mins(22, 48))] dataset = []
def eff_error(E): return sqrt((E * scales[1, 0]) ** 2 + (scales[1, 1]) ** 2) def cross_section(theta): return r0 ** 2 * (1 / (1 + alpha * (1 - cos(theta)))) ** 2 * \ (1 + cos(theta) ** 2) / 2 * (1 + (alpha ** 2 * (1 - cos(theta)) ** 2) \ / ((1 + cos(theta) ** 2) * \ (1 + alpha * \ (1 - cos(theta))))) # Flux phi R = 0.19 + 0.0045 # Distance from preparate to electron detector R_error = 0.005 l = log(2) / (30.17 * yrs(1, 1, 0)) l_error = l / (30.17 * yrs(1, 1, 0)) * 24 * 3600 * 0.5 n0 = 122100e3 n0_error = 50e3 t = yrs(2013, 1, 3) t_error = 24 * 3600 * 0.5 t0 = yrs(1979, 3, 22) t0_error = 24 * 3600 * 0.5 activity = n0 * exp( -l * (t - t0)) print "activity is %d" % activity phi = activity * 0.851 / (4 * pi * R ** 2) phi_error = 0.851 / (4 * pi) * sqrt( (2 * activity / R ** 3 * R_error) ** 2 + (activity / n0 / R ** 2 * n0_error) ** 2 + (activity * (t - t0) / R ** 2 * l_error) ** 2 +