24 from array
import array
27 from ROOT.Belle2
import SVDCoGCalibrationFunction
28 from ROOT.Belle2
import SVDCoGTimeCalibrations
32 import matplotlib.pyplot
as plt
38 branches = [
'SVDShaperDigits',
'SVDShaperDigitsFromTracks',
'EventT0']
40 trk_outputFile =
"TrackFilterControlNtuples_" + localdb +
".root"
47 "/group/belle2/dataprod/Data/release-03-02-02/DB00000635/proc00000009/\
48 e0008/4S/r01309/skim/hlt_bhabha/cdst/sub00/cdst.physics.0008.01309.HLT*"
51 if filename.rfind(".root") != -1:
52 inputFileList.append(filename)
54 with open(filename, 'r') as f:
55 inputFileList = [line.strip() for line in f]
59 use_central_database(
'data_reprocessing_proc9')
61 if '_CHECK' not in localdb:
62 use_central_database(
'svd_NOCoGCorrections')
63 if os.path.isdir(localdb):
64 print(
'WARNING, you are using CoG corrections in ' + localdb +
' exiting')
67 use_local_database(localdb +
"/database.txt", localdb, invertLogging=
True)
69 localdb = localdb.split(
'_')
70 if not os.path.isdir(localdb[0]):
71 print(
'WARNING, the localDB ' + localdb[0] +
' containing the corrections is missing, exiting.')
74 use_local_database(localdb[0] +
"/database.txt", localdb[0], invertLogging=
True)
76 if os.path.isdir(localdb):
77 print(
'WARNING, you are using CoG corrections in ' + localdb +
' exiting')
80 use_local_database(localdb +
"/database.txt", localdb, invertLogging=
True)
86 rootinput = register_module(
'RootInput')
87 rootinput.param(
'inputFileNames', inputFileList)
88 rootinput.param(
'branchNames', branches)
89 main.add_module(rootinput)
91 main.add_module(
"Gearbox")
92 main.add_module(
"Geometry", useDB=
True)
95 trkFlt = register_module(
'TrackFilter')
96 trkFlt.param(
'outputFileName', trk_outputFile)
97 trkFlt.param(
'outputINArrayName',
'SelectedTracks')
98 trkFlt.param(
'outputOUTArrayName',
'ExcludedTracks')
99 trkFlt.param(
'min_NumHitSVD', nSVD)
100 trkFlt.param(
'min_NumHitCDC', nCDC)
101 trkFlt.param(
'min_Pvalue', pVal)
109 add_svd_reconstruction(main)
111 for moda
in main.modules():
112 if moda.name() ==
'SVDCoGTimeEstimator':
113 moda.param(
"ShaperDigits",
'SVDShaperDigitsFromTracks')
114 moda.param(
"RecoDigits",
'SVDRecoDigitsFromTracks')
115 if moda.name() ==
'SVDSimpleClusterizer':
116 moda.param(
"Clusters",
'SVDClustersFromTracks')
117 moda.param(
"RecoDigits",
'SVDRecoDigitsFromTracks')
118 if moda.name() ==
'SVDSpacePointCreator':
119 moda.param(
"SVDClusters",
'SVDClustersFromTracks')
121 calib = SVDCoGTimeCalibrationImporterModule()
122 calib.set_localdb(localdb)
123 calib.notApplyCorrectForCDCLatency(
False)
124 main.add_module(calib)
127 progress = register_module(
'ProgressBar')
128 main.add_module(progress)