make MCParticles from reconstructed cosmic tracks
Definition at line 16 of file makeMCParticlesFromGCRData.py.
◆ event()
event processing
Definition at line 25 of file makeMCParticlesFromGCRData.py.
26 ''' event processing '''
30 absPDG = abs(Belle2.Const.muon.getPDGCode())
31 mass = Belle2.Const.muon.getMass()
33 extHits = track.getRelationsWith(
'ExtHits')
35 for extHit
in extHits:
36 if extHit.getDetectorID() != Belle2.Const.TOP:
38 if abs(extHit.getPdgCode()) != absPDG:
40 if extHit.getPosition().Perp() < 123.5:
42 if extHit.getPosition() * extHit.getMomentum() > 0:
44 selExtHits.append(extHit)
45 if len(selExtHits) == 0:
47 sortedExtHits = sorted(selExtHits, key=
lambda x: (x.getTOF()))
48 extHit = sortedExtHits[0]
49 pos = extHit.getPosition()
50 mom = extHit.getMomentum()
52 part = mcParticles.appendNew()
53 part.setPDG(extHit.getPdgCode())
55 part.setProductionVertex(pos)
56 part.setProductionTime(tof)
58 part.setEnergy(math.sqrt(mom.Mag2() + mass**2))
59 part.setValidVertex(
True)
60 part.setStatus(Belle2.MCParticle.c_PrimaryParticle)
61 part.addStatus(Belle2.MCParticle.c_StableInGenerator)
63 if mcParticles.getEntries() > 0:
◆ initialize()
The documentation for this class was generated from the following file: