10 from optparse
import OptionParser
15 set_log_level(LogLevel.INFO)
17 outroot =
'arichbtest.root'
19 parser = OptionParser()
20 parser.add_option(
'-r',
'--run', dest=
'runno', default=
'068',
22 parser.add_option(
'-p',
'--path', dest=
'path',
23 default=
'/afs/f9.ijs.si/data/belle/data/beamtests/2011/run',
26 parser.add_option(
'-y',
'--year', dest=
'year', default=
'2011',
27 help=
'beam test year')
29 parser.add_option(
'-t',
'--track-mask', dest=
'mask', default=
'0x5',
30 help=
'track mask 4 bits')
32 parser.add_option(
'-m',
'--avgagel', dest=
'avgagel', default=
'0',
33 help=
'average thc calculation based on the mean agel')
40 help=
'Output filename',
44 parser.add_option(
'-n',
'--neve', dest=
'neve', default=20000,
45 help=
'Number of events to process')
47 (options, args) = parser.parse_args()
49 fname = options.path +
'/run' + options.runno +
'.dat'
51 if not os.path.isfile(fname):
52 print(
'File does not exist:')
55 print(
'Appending .gz...')
56 if not os.path.isfile(fname):
57 print(
'File does not exist:')
61 if options.output == outroot:
62 outroot =
'run_' + options.runno +
'.root'
64 outroot = options.output
66 mask = int(options.mask, 16)
67 runno = int(options.runno)
69 print(
'TrackMask:' + options.mask)
72 averageagel = int(options.avgagel)
74 eventinfosetter = register_module(
'EventInfoSetter')
75 eventinfosetter.param(
'evtNumList', [int(options.neve)])
76 eventinfosetter.param(
'runList', [int(options.runno)])
77 eventinfosetter.param(
'expList', [1])
80 paramloader = register_module(
'Gearbox')
82 xmlgeometry =
'file://%s/arich/modules/arichBtest/data/%s/arichBtest%s.xml' \
83 % (os.getcwd(), options.year, options.year)
84 paramloader.param(
'fileName', xmlgeometry)
86 paramloader.param(
'fileName', xmlgeometry)
99 geobuilder = register_module(
'Geometry')
100 geobuilder.param(
'components', [
'ARICHBtest'])
102 btest = register_module(
'arichBtest')
103 btest.param(
'mwpcTrackMask', [mask])
105 btest.param(
'runList', [fname])
106 btest.param(
'outputFileName', outroot)
108 if options.year ==
'2013':
110 print(
'Beam momentum ' + str(momentum))
111 btest.param(
'beamMomentum', momentum)
114 g4sim = register_module(
'FullSim')
117 g4sim.param(
'RegisterOptics', 1)
121 g4sim.param(
'PhotonFraction', 0.3)
123 g4sim.param(
'EnableVisualization',
True)
127 g4sim.param(
'UICommandsAtIdle', [
'/vis/open VRML2FILE',
'/vis/drawVolume',
128 '/vis/scene/add/axes 0 0 0 100 mm'])
131 geosaver = register_module(
'ExportGeometry')
132 geosaver.param(
'Filename',
'Belle2Geo.root')
134 arichrec = register_module(
'ARICHReconstructor')
135 arichrec.param(
'inputTrackType', 1)
136 arichrec.param(
'beamtest', 3)
137 arichrec.param(
'trackPositionResolution', 0.0)
138 arichrec.param(
'trackAngleResolution', 0.0)
140 profile = register_module(
'Profile')
141 profile.param(
'outputFileName',
'profileusage.ps')
144 main.add_module(eventinfosetter)
145 main.add_module(paramloader)
147 main.add_module(geobuilder)
148 main.add_module(btest)
149 main.add_module(arichrec)
151 main.add_module(geosaver)
156 print(
'Event Statistics:')