MakeMCParticles Class Reference

Inheritance diagram for MakeMCParticles:

Collaboration diagram for MakeMCParticles:

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

Detailed Description

 make MCParticles from reconstructed cosmic tracks 

Definition at line 24 of file makeMCParticlesFromGCRData.py.

Member Function Documentation

event()

def event

(

self

)

 event processing 

Definition at line 33 of file makeMCParticlesFromGCRData.py.

    def event(self):
        ''' event processing '''
 
        mcParticles = Belle2.PyStoreArray('MCParticles')
        tracks = Belle2.PyStoreArray('Tracks')
        absPDG = abs(Belle2.Const.muon.getPDGCode())
        mass = Belle2.Const.muon.getMass()
        for track in tracks:
            extHits = track.getRelationsWith('ExtHits')
            selExtHits = []
            for extHit in extHits:
                if extHit.getDetectorID() != Belle2.Const.TOP:
                    continue
                if abs(extHit.getPdgCode()) != absPDG:
                    continue
                if extHit.getPosition().Perp() < 123.5:
                    continue
                if extHit.getPosition() * extHit.getMomentum() > 0:
                    continue
                selExtHits.append(extHit)
            if len(selExtHits) == 0:
                continue
            sortedExtHits = sorted(selExtHits, key=lambda x: (x.getTOF()))
            extHit = sortedExtHits[0]
            pos = extHit.getPosition()
            mom = extHit.getMomentum()
            tof = extHit.getTOF()
            part = mcParticles.appendNew()
            part.setPDG(extHit.getPdgCode())
            part.setMass(mass)
            part.setProductionVertex(pos)
            part.setProductionTime(tof)
            part.setMomentum(mom)
            part.setEnergy(math.sqrt(mom.Mag2() + mass**2))
            part.setValidVertex(True)
            part.setStatus(Belle2.MCParticle.c_PrimaryParticle)
            part.addStatus(Belle2.MCParticle.c_StableInGenerator)
 
        if mcParticles.getEntries() > 0:
            self.return_value(1)
        else:
            self.return_value(0)
 
 
# Create path

initialize()

def initialize

(

self

)

 initialization 

Definition at line 27 of file makeMCParticlesFromGCRData.py.

---

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

• top/analysis/makeMCParticlesFromGCRData.py