def test9k(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) pat = mytof.Particle(13, -1., 0.1056583745, -.1, 0., .02, 0., 0.33, -1.5) tof = mytof.Tof(.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_9k.txt')
def test9c(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) pat = mytof.Particle(13, -1., 0.1056583745, .10 / 2**.5, .10 / 2**.5, .1, 0.0, 0.00, -1.5) tof = mytof.Tof(1.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_9c.txt')
def test8b(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) pat = mytof.Particle(130, 0., pdg2m(130), .5 / 2**.5, .5 / 2**.5, 1., 0.5, 0.25, 1.5) tof = mytof.Tof(1.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_8.txt')
def test6b(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) pat = mytof.Particle(-13, 1., 0.1056583745, .3 / 2**.5, .3 / 2**.5, -1., -.90 / 2**.5, .90 / 2**.5, 1.5) tof = mytof.Tof(1.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_6b.txt')
def test4(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) pat = mytof.Particle(13, -1., 0.1056583745, .5 / 2**.5, .5 / 2**.5, 1., 0.5, 0.25, 1.75) tof = mytof.Tof(.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_4.txt')
def test6c(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) phi = math.pi / 180 * 88. pat = mytof.Particle(-13, 1., 0.1056583745, .3 * math.cos(phi), .3 * math.sin(phi), -.1, -.99, 0.001, 1.5) tof = mytof.Tof(1.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_6c.txt')
def test1(): os.system('rm -rf tof_save.txt') field = mytof.Field(1., 1., 2.) # def __init__(self, pdgid, charge, m, px, py, pz, vx=0., vy=0., vz=0.): pat = mytof.Particle(13, -1., 0.1056583745, .5 / 2**.5, .5 / 2**.5, 0., 0., 0., 0.) # def __init__(self, r_in, r_out, z, sigma, eff): tof = mytof.Tof(.5, 2., 3., 25.e-9, .999) print tof.get_time_v2(pat, field) os.system('cp -f tof_save.txt tof_test_1.txt')
def beforeFill_pid(n0, pdg0, q0, en0, px0, py0, pz0, tof0, mu0, tof, mu, tofdet, mu0det, field): for i in xrange(n0[0]): ptcl = mytof.Particle(pdg0[i], q0[i], mytof.MASS[pdg0[i]], px0[i], py0[i], pz0[i]) tof0[i] = tofdet.get_time( ptcl, field ) mu0[i] = mu0det.get_time( ptcl, field )