SVDRecoDigitsQuickCheck Class Reference

Inheritance diagram for SVDRecoDigitsQuickCheck:

Public Member Functions
def	initialize (self)
def	event (self)
def	terminate (self)

Detailed Description

quick check of SVDRecoDigits

Definition at line 267 of file testNewReco.py.

Member Function Documentation

event()

def event

(

self

)

look at old and new reco digits

Definition at line 302 of file testNewReco.py.

    def event(self):
        '''look at old and new reco digits'''
 
        recodigitList = Belle2.PyStoreArray("SVDRecoDigits")
        recodigitListNew = Belle2.PyStoreArray("SVDNewRecoDigits")
 
        print("number of recodigits = " + str(recodigitList.getEntries()) +
              " (old) VS " + str(recodigitListNew.getEntries()) + " (new)")
 
        for d in recodigitList:
            self.time.Fill(d.getTime())
            self.charge.Fill(d.getCharge())
            isU = 0
            if(d.isUStrip()):
                isU = 0.5
            self.test[d.getSensorID().getLayerNumber() - 3].Fill(d.getSensorID().getLadderNumber(),
                                                                 d.getSensorID().getSensorNumber() + isU)
 
        for d in recodigitListNew:
            self.timeNew.Fill(d.getTime())
            self.chargeNew.Fill(d.getCharge())
            isU = 0
            if(d.isUStrip()):
                isU = 0.5
            self.testNew[d.getSensorID().getLayerNumber() - 3].Fill(d.getSensorID().getLadderNumber(),
                                                                    d.getSensorID().getSensorNumber() + isU)
 

initialize()

def initialize

(

self

)

define histograms

Definition at line 270 of file testNewReco.py.

    def initialize(self):
        '''define histograms'''
 
        
        self.test = []
        self.testNew = []
        self.time = TH1F("rd_time", "RecoDigit Time", 300, -100, 200)
        self.timeNew = TH1F("rdNew_time", "New RecoDigit Time", 300, -100, 200)
        self.charge = TH1F("rd_charge", "RecoDigit Charge", 300, 0, 100000)
        self.chargeNew = TH1F("rdNew_charge", "New RecoDigit Charge", 300, 0, 100000)
        
 
        geoCache = Belle2.VXD.GeoCache.getInstance()
 
        ladderN = 0
        sensorN = 0
 
        for layer in geoCache.getLayers(Belle2.VXD.SensorInfoBase.SVD):
            layerNumber = layer.getLayerNumber()
            ladderN = 0
            for ladder in geoCache.getLadders(layer):
                ladderN = ladderN + 1
                sensorN = 0
                for sensor in geoCache.getSensors(ladder):
                    sensorN = sensorN + 1
            self.test.append(TH2F("rd_layer" + str(layerNumber), "Layer " + str(layerNumber) +
                                  " Ladder VS Sensor.Side", ladderN, 0.5, ladderN + 0.5, 2 * sensorN, +0.75, sensorN + 0.75))
            self.testNew.append(TH2F("rdNew_layer" + str(layerNumber), "Layer " + str(layerNumber) +
                                     " Ladder VS Sensor.Side", ladderN, 0.5, ladderN + 0.5, 2 * sensorN, +0.75, sensorN + 0.75))
 
        print(self.test)
 

terminate()

def terminate

(

self

)

write

Definition at line 329 of file testNewReco.py.

    def terminate(self):
        '''write'''
 
        f = TFile("quicktestSVDRecoDigitOldDefault.root", "RECREATE")
        for hist in self.test:
            hist.GetXaxis().SetTitle("ladder #")
            hist.GetYaxis().SetTitle("sensor # + 0.5 is isU")
            hist.Write()
        for hist in self.testNew:
            hist.GetXaxis().SetTitle("ladder #")
            hist.GetYaxis().SetTitle("sensor # + 0.5 is isU")
            hist.Write()
 
        self.time.GetXaxis().SetTitle("recodigit time (ns)")
        self.time.Write()
        self.timeNew.GetXaxis().SetTitle("recodigit time (ns)")
        self.timeNew.Write()
 
        self.charge.GetXaxis().SetTitle("recodigit charge (ke-)")
        self.charge.Write()
        self.chargeNew.GetXaxis().SetTitle("recodigit charge (ke-)")
        self.chargeNew.Write()
 
        f.Close()
 
# b2conditions.prepend_globaltag("svd_test_svdRecoConfiguration")
 
 

---

The documentation for this class was generated from the following file:

• svd/examples/testNewReco.py