13 from rawdata
import add_unpackers
14 from svd
import add_svd_reconstruction
15 from ROOT
import Belle2
20 print(
"usage: basf2", argvs[0],
"expNo runNo globalTag [outputFolder]")
23 expNo =
'e' +
'{:0=4d}'.format(int(argvs[1]))
24 runNo =
'r' +
'{:0=5d}'.format(int(argvs[2]))
28 indir =
'/hsm/belle2/bdata/Data/Raw/' + expNo +
'/' + runNo +
'/sub00'
29 files = sorted(glob.glob(indir +
'/*.root'))
31 basf2.B2ERROR(
'No files found in ' + indir)
37 outputFile = outdir +
'/beamBkgHitRates-' + expNo +
'-' + runNo +
'.root'
41 ''' select events according to given trigger types '''
44 ''' event processing '''
48 if not trg_summary.isValid():
49 basf2.B2ERROR(
'No TRGSummary available - event ignored')
52 for trgType
in trigTypes:
53 if trg_summary.getTimType() == trgType:
59 basf2.reset_database()
60 basf2.use_central_database(global_tag)
63 main = basf2.create_path()
64 emptypath = basf2.create_path()
67 main.add_module(
'RootInput', inputFileNames=files)
70 main.add_module(
'Gearbox')
73 main.add_module(
'Geometry')
76 main.add_module(
'TRGGDLUnpacker')
77 main.add_module(
'TRGGDLSummary')
80 main.add_module(
'Progress')
83 if len(trigTypes) > 0:
85 main.add_module(selector)
86 selector.if_false(emptypath)
89 add_unpackers(path=main,
90 components=[
'PXD',
'SVD',
'CDC',
'ECL',
'TOP',
'ARICH',
'KLM'])
93 main.add_module(
'ARICHFillHits')
94 main.add_module(
'TOPChannelMasker')
95 main.add_module(
'ActivatePXDGainCalibrator')
96 main.add_module(
'PXDClusterizer')
97 add_svd_reconstruction(main, applyMasking=
True)
99 'SVDZeroSuppressionEmulator',
101 ShaperDigits=
'SVDShaperDigitsUnmasked',
102 ShaperDigitsIN=
'SVDShaperDigitsZS5')
103 main.add_module(
'TFCDC_WireHitPreparer',
104 wirePosition=
'aligned',
105 flightTimeEstimation=
'outwards',
106 filter=
'cuts_from_DB')
110 'BeamBkgHitRateMonitor',
111 outputFileName=outputFile,
113 svdShaperDigitsName=
'SVDShaperDigitsZS5',
114 cdcTimeWindowLowerEdgeSmallCell=4720,
115 cdcTimeWindowUpperEdgeSmallCell=5020,
116 cdcTimeWindowLowerEdgeNormalCell=4720,
117 cdcTimeWindowUpperEdgeNormalCell=5020)
123 print(basf2.statistics)