Ntuple Class Reference

Inheritance diagram for Ntuple:

Collaboration diagram for Ntuple:

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

Public Attributes
	file
	root file
	tree
	root tree
	data
	data structure to be written to tree
	reso
	histogram of t0 residuals
	pull
	histogram of t0 pulls

Detailed Description

makes ntuple for MC studies of TOPCosmicT0Finder performance 

Definition at line 58 of file testCosmicT0Finder.py.

Member Function Documentation

event()

def event

(

self

)

event processing: fill ntuple and histograms 

Definition at line 84 of file testCosmicT0Finder.py.

    def event(self):
        ''' event processing: fill ntuple and histograms '''
 
        mcParticles = Belle2.PyStoreArray('MCParticles')
        timeZeros = Belle2.PyStoreArray('TOPTimeZeros')
        for timeZero in timeZeros:
            extHit = timeZero.getRelated('ExtHits')
            if not extHit:
                B2WARNING('no related extHits')
                continue
            barHit = None
            for barhit in mcParticles[0].getRelationsWith('TOPBarHits'):
                if barhit.getModuleID() == extHit.getCopyID():
                    barHit = barhit
            if not barHit:
                B2WARNING('no corresponding TOPBarHit')
                continue
 
            moduleID = timeZero.getModuleID()
            if extHit.getCopyID() != moduleID:
                B2ERROR('moduleID in extHit differs: ' + str(extHit.getCopyID()) +
                        ' ' + str(moduleID))
            if barHit.getModuleID() != moduleID:
                B2ERROR('moduleID in barHit differs: ' + str(barHit.getModuleID()) +
                        ' ' + str(moduleID))
 
            self.data.slot = moduleID
            self.data.nfot = timeZero.getNumPhotons()
            self.data.pdg = barHit.getPDG()
            pos = barHit.getLocalPosition()
            self.data.x = pos.x()
            self.data.y = pos.y()
            self.data.z = pos.z()
            self.data.t = barHit.getTime()
            phiBar = math.pi / 2 - math.pi / 8 * (moduleID - 0.5)
            mom = barHit.getMomentum()
            mom.RotateZ(phiBar)
            self.data.p = mom.Mag()
            self.data.theta = math.degrees(mom.Theta())
            self.data.phi = math.degrees(mom.Phi())
            dr = extHit.getPosition() - barHit.getPosition()
            dr.RotateZ(phiBar)
            self.data.dx = dr.x()
            self.data.dy = dr.y()
            self.data.dz = dr.z()
            self.data.dt = extHit.getTOF() - self.data.t
            momExt = extHit.getMomentum()
            momExt.RotateZ(phiBar)
            self.data.dp = momExt.Mag() - self.data.p
            self.data.dtheta = math.degrees(momExt.Theta()) - self.data.theta
            self.data.dphi = math.degrees(momExt.Phi()) - self.data.phi
            self.data.t0 = timeZero.getTime()
            self.data.t0err = timeZero.getError()
            self.data.timeShift = timeShift
            self.data.dt0 = self.data.dt + self.data.t0 - self.data.timeShift
            self.file.cd()
            self.tree.Fill()
            self.reso.Fill(self.data.dt0)
            self.pull.Fill(self.data.dt0/self.data.t0err)
            if timeZero.getChi2().GetEntries() > 0:
                timeZero.getChi2().Write()
                timeZero.getPDF().Write()
                timeZero.getHits().Write()
 

initialize()

def initialize

(

self

)

initialize 

Definition at line 61 of file testCosmicT0Finder.py.

terminate()

def terminate

(

self

)

terminate: write and close root file 

Definition at line 148 of file testCosmicT0Finder.py.

---

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

• top/examples/testCosmicT0Finder.py