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 )