SVDClustersDebug Class Reference

Inheritance diagram for SVDClustersDebug:

Collaboration diagram for SVDClustersDebug:

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

Public Attributes
	fb2u

Detailed Description

class to debug clusterizer

Definition at line 40 of file SVDClustersDebug.py.

Member Function Documentation

event()

def event

(

self

)

event

Definition at line 51 of file SVDClustersDebug.py.

    def event(self):
        '''event'''
 
        evt = Belle2.PyStoreObj("EventMetaData")
 
        clustersZS3 = Belle2.PyStoreArray(cluster3)
        clustersZS5 = Belle2.PyStoreArray(cluster5)
 
        nCl3 = 0
        nCl5 = 0
 
        for d in clustersZS5:
            if(str(d.getSensorID()) == "4.5.2"):
                if(d.isUCluster()):
                    nCl5 = nCl5 + 1
 
        for d in clustersZS3:
            if(str(d.getSensorID()) == "4.5.2"):
                if(d.isUCluster()):
                    nCl3 = nCl3 + 1
 
        self.fb2u.Fill(nCl3, nCl5)
 
        noisePayload = Belle2.SVDNoiseCalibrations()
#        print(type(noisePayload))
 
        if nCl3 != nCl5:
            print("different number of clusters in event " + str(evt.getEvent()) +
                  "! ZS3 = " + str(nCl3) + ", ZS5 = " + str(nCl5) + "    ")
 
            # for d in clustersZS5:
            #    if(str(d.getSensorID()) == "4.5.2") :
            #        if(d.isUCluster()):
            #            print(" ZS5: seed charge = "+str(d.getSeedCharge())+",\
            # size = "+str(d.getSize())+", position = "+str(d.getPosition())+", charge = "+str(d.getCharge())+"   ")
 
            # for d in clustersZS3:
            #    if(str(d.getSensorID()) == "4.5.2") :
            #        if(d.isUCluster()):
            #            print(" ZS3: seed charge = "+str(d.getSeedCharge())+",\
            # size = "+str(d.getSize())+", position = "+str(d.getPosition())+", charge = "+str(d.getCharge())+"   ")
 
            for d3 in clustersZS3:
                if(str(d3.getSensorID()) == "4.5.2"):
                    if(d3.isUCluster()):
                        found = False
                        for d5 in clustersZS5:
                            if(str(d5.getSensorID()) == "4.5.2"):
                                if(d5.isUCluster()):
                                    if(d5.getSeedCharge() == d3.getSeedCharge()):
                                        found = True
                        if not found:
                            print(" ONLY ZS3: seed charge = " + str(d3.getSeedCharge()) + ", size = " + str(d3.getSize()) +
                                  ", position = " + str(d3.getPosition()) + ", charge = " + str(d3.getCharge()) + "   ")
                            reco = d3.getRelationsTo(reco3)
#                            print("test size cluster: "+str(reco.size()))
                            for r in reco:
                                charge = r.getCharge()
                                noise = noisePayload.getNoiseInElectrons(r.getSensorID(), True, r.getCellID())
                                print(" SNR = " + str(charge / noise) + ", charge = " + str(charge) + ", noise = " + str(noise))
 

initialize()

def initialize

(

self

)

initialize

Definition at line 43 of file SVDClustersDebug.py.

terminate()

def terminate

(

self

)

terminate

Definition at line 112 of file SVDClustersDebug.py.

---

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

• svd/scripts/SVDClustersDebug.py