5 SVD Default PulseShape Calibration importer (MC).
6 Script to Import Calibrations into a local DB
12 from ROOT
import Belle2
13 from ROOT.Belle2
import SVDStripCalAmp
18 now = datetime.datetime.now()
40 gain_origami_U = 271.4
41 gain_origami_V = 218.8
48 peakTime_origami_U = 66
49 peakTime_origami_V = 52
55 '''default pulse shape calibrations importer'''
63 tmp_calAmp = SVDStripCalAmp()
65 tmp_calAmp.peakTime = 75
66 tmp_calAmp.pulseWidth = pulseWidth
68 tmp_calAmp,
"PulseShapeCalibrations_default_" + str(now.isoformat()) +
69 "_INFO:_peakTime=fromPhase3calibrations_pulseWidth=130_gain=fromPhase3calibrations")
73 for layer
in geoCache.getLayers(Belle2.VXD.SensorInfoBase.SVD):
74 layerNumber = layer.getLayerNumber()
75 for ladder
in geoCache.getLadders(layer):
76 ladderNumber = ladder.getLadderNumber()
77 for sensor
in geoCache.getSensors(ladder):
78 sensorNumber = sensor.getSensorNumber()
81 print(
"setting PulseShape for " +
82 str(layerNumber) +
"." + str(ladderNumber) +
"." + str(sensorNumber) +
"." + str(side))
86 peakTime = peakTime_L3_V
91 peakTime = peakTime_fwd_V
93 if sensorNumber == layerNumber - 1:
95 peakTime = peakTime_bkw_V
98 peakTime = peakTime_origami_V
102 peakTime = peakTime_L3_U
104 if sensorNumber == 1:
106 peakTime = peakTime_fwd_U
108 if sensorNumber == layerNumber - 1:
110 peakTime = peakTime_bkw_U
112 gain = gain_origami_U
113 peakTime = peakTime_origami_U
115 tmp_calAmp.gain = 1 / gain
116 tmp_calAmp.peakTime = peakTime
119 for strip
in range(0, Nstrips):
122 calAmp_payload.set(layerNumber, ladderNumber, sensorNumber, bool(side), strip, tmp_calAmp)
128 use_central_database(
"svd_onlySVDinGeoConfiguration")
129 use_local_database(
"localDB_defaultPulseShapeCalibrations/database.txt",
"localDB_defaultPulseShapeCalibrations")
134 eventinfosetter = register_module(
'EventInfoSetter')
135 eventinfosetter.param({
'evtNumList': [1],
'expList': 0,
'runList': 0})
136 main.add_module(eventinfosetter)
138 main.add_module(
"Gearbox")
139 main.add_module(
"Geometry", components=[
'SVD'])
144 progress = register_module(
'Progress')
145 main.add_module(progress)