SVDDefaultPulseShapeImporter.py

#!/usr/bin/env python
# -*- coding: utf-8 -*-
 
 
 
"""
SVD Default PulseShape Calibration importer (MC).
Script to Import Calibrations into a local DB
"""
import basf2 as b2
from ROOT import Belle2
from ROOT.Belle2 import SVDStripCalAmp
import datetime
 
 
now = datetime.datetime.now()
 
pulseWidth = 130
 
# gain here is actually ADUequivalent (e-/ADC)
# we set the INVERSE of it in the payload, so that the payload actually contains the gain
'''
# Phase 3 - old
gain_L3_U = 375
gain_L3_V = 375
gain_bkw_U = 375
gain_bkw_V = 375
gain_origami_U = 375
gain_origami_V = 375
gain_fwd_U = 375
gain_fwd_V = 375
'''
# Phase 3 - matching data
gain_L3_U = 282.5
gain_L3_V = 275.5
gain_bkw_U = 250.1
gain_bkw_V = 246.6
gain_origami_U = 271.4
gain_origami_V = 218.8
gain_fwd_U = 247.4
gain_fwd_V = 245.7
peakTime_L3_U = 67
peakTime_L3_V = 58
peakTime_bkw_U = 66
peakTime_bkw_V = 52
peakTime_origami_U = 66
peakTime_origami_V = 52
peakTime_fwd_U = 60
peakTime_fwd_V = 51
 
 
class defaultPulseShapeImporter(b2.Module):
    '''default pulse shape calibrations importer'''
 
    def beginRun(self):
        '''begin run'''
 
        iov = Belle2.IntervalOfValidity.always()
 
        # gain, peakTime,
        tmp_calAmp = SVDStripCalAmp()
        tmp_calAmp.gain = 275  # set after the loop on sensors
        tmp_calAmp.peakTime = 75  # set after the loop on sensors
        tmp_calAmp.pulseWidth = pulseWidth
        calAmp_payload = Belle2.SVDPulseShapeCalibrations.t_calAmp_payload(
            tmp_calAmp, "PulseShapeCalibrations_default_" + str(now.isoformat()) +
            "_INFO:_peakTime=fromPhase3calibrations_pulseWidth=130_gain=fromPhase3calibrations")
 
        geoCache = Belle2.VXD.GeoCache.getInstance()
 
        for layer in geoCache.getLayers(Belle2.VXD.SensorInfoBase.SVD):
            layerNumber = layer.getLayerNumber()
            for ladder in geoCache.getLadders(layer):
                ladderNumber = ladder.getLadderNumber()
                for sensor in geoCache.getSensors(ladder):
                    sensorNumber = sensor.getSensorNumber()
                    for side in (0, 1):
                        Nstrips = 768
                        print("setting PulseShape for " +
                              str(layerNumber) + "." + str(ladderNumber) + "." + str(sensorNumber) + "." + str(side))
                        if side == 0:  # V
                            if layerNumber == 3:  # L3 V
                                gain = gain_L3_V
                                peakTime = peakTime_L3_V
                            else:
                                Nstrips = 512
                                if sensorNumber == 1:  # FW V
                                    gain = gain_fwd_V
                                    peakTime = peakTime_fwd_V
                                else:  # BKW V
                                    if sensorNumber == layerNumber - 1:  # FW V
                                        gain = gain_bkw_V
                                        peakTime = peakTime_bkw_V
                                    else:  # BARREL V
                                        gain = gain_origami_V
                                        peakTime = peakTime_origami_V
                        if side == 1:  # U
                            if layerNumber == 3:  # L3 U
                                gain = gain_L3_U
                                peakTime = peakTime_L3_U
                            else:
                                if sensorNumber == 1:  # FW U
                                    gain = gain_fwd_U
                                    peakTime = peakTime_fwd_U
                                else:  # BKW U
                                    if sensorNumber == layerNumber - 1:  # FW U
                                        gain = gain_bkw_U
                                        peakTime = peakTime_bkw_U
                                    else:  # BARREL U
                                        gain = gain_origami_U
                                        peakTime = peakTime_origami_U
 
                        tmp_calAmp.gain = 1 / gain
                        tmp_calAmp.peakTime = peakTime
 
                        # print(str(Nstrips))
                        for strip in range(0, Nstrips):
                            # print("setting Gain for strip " + str(strip) + " to " + str(tmp_calAmp.gain))
 
                            calAmp_payload.set(layerNumber, ladderNumber, sensorNumber, bool(side), strip, tmp_calAmp)
 
        Belle2.Database.Instance().storeData(Belle2.SVDPulseShapeCalibrations.calAmp_name, calAmp_payload, iov)
 
 
b2.conditions.prepend_globaltag("svd_onlySVDinGeoConfiguration")
 
main = b2.create_path()
 
# Event info setter - execute single event
eventinfosetter = b2.register_module('EventInfoSetter')
eventinfosetter.param({'evtNumList': [1], 'expList': 0, 'runList': 0})
main.add_module(eventinfosetter)
 
main.add_module("Gearbox")
main.add_module("Geometry")
 
main.add_module(defaultPulseShapeImporter())
 
# Show progress of processing
progress = b2.register_module('Progress')
main.add_module(progress)
 
# Process events
b2.process(main)