21 direction =
'incoming'
25 outputFile =
'comparePDF.root'
30 if direction ==
'incoming':
35 alpha = (slot - 0.5) * 360 / 16 - 90
36 ca = math.cos(math.radians(alpha))
37 sa = math.sin(math.radians(alpha))
38 x = xloc * ca - (yloc + 120) * sa
39 y = xloc * sa + (yloc + 120) * ca
44 set_log_level(LogLevel.WARNING)
50 eventinfosetter = register_module(
'EventInfoSetter')
51 eventinfosetter.param(
'evtNumList', [1000])
52 main.add_module(eventinfosetter)
55 histo = register_module(
'HistoManager')
56 histo.param(
'histoFileName', outputFile)
57 main.add_module(histo)
60 gearbox = register_module(
'Gearbox')
61 main.add_module(gearbox)
64 geometry = register_module(
'Geometry')
65 geometry.param(
'useDB',
False)
66 geometry.param(
'components', [
'MagneticField',
'TOP'])
67 main.add_module(geometry)
70 particlegun = register_module(
'ParticleGun')
71 particlegun.param(
'pdgCodes', [pdg])
72 particlegun.param(
'nTracks', 1)
73 particlegun.param(
'varyNTracks',
False)
74 particlegun.param(
'momentumGeneration',
'fixed')
75 particlegun.param(
'momentumParams', [p])
76 particlegun.param(
'thetaGeneration',
'fixed')
77 particlegun.param(
'thetaParams', [theta])
78 particlegun.param(
'phiGeneration',
'fixed')
79 particlegun.param(
'phiParams', [phi])
80 particlegun.param(
'vertexGeneration',
'fixed')
81 particlegun.param(
'xVertexParams', [x])
82 particlegun.param(
'yVertexParams', [y])
83 particlegun.param(
'zVertexParams', [z])
84 particlegun.param(
'independentVertices',
False)
85 main.add_module(particlegun)
88 simulation = register_module(
'FullSim')
89 main.add_module(simulation)
92 topdigi = register_module(
'TOPDigitizer')
93 topdigi.param(
'useWaveforms',
False)
94 topdigi.param(
'simulateTTS',
False)
95 topdigi.param(
'electronicJitter', 0.0)
96 topdigi.param(
'timeZeroJitter', 0.0)
97 main.add_module(topdigi)
100 trackmaker = register_module(
'TOPMCTrackMaker')
101 main.add_module(trackmaker)
104 toppdf = register_module(
'TOPPDFChecker')
105 main.add_module(toppdf)
108 progress = register_module(
'Progress')
109 main.add_module(progress)