g4Ex_gap.MyTrackingActionD Class Reference

Inheritance diagram for g4Ex_gap.MyTrackingActionD:

Collaboration diagram for g4Ex_gap.MyTrackingActionD:

Public Member Functions
	PostUserTrackingAction (self, atrack)
	PreUserTrackingAction (self, atrack)
	myPrintout (self, atrack)

Detailed Description

Definition at line 309 of file g4Ex_gap.py.

Member Function Documentation

myPrintout()

g4Ex_gap.MyTrackingActionD.myPrintout	(		self,
			atrack
	)

Definition at line 364 of file g4Ex_gap.py.

 def myPrintout(self, atrack):
   part         = atrack.GetDynamicParticle()
   pid          = part.GetPDGcode()
   vx  = atrack.GetVertexPosition()
   pos = atrack.GetPosition()
   mom = atrack.GetMomentum()
   print('TA',pid,atrack.GetTotalEnergy()/GeV,ecut*GeV) 
   print('start tracking',atrack.GetDynamicParticle().GetPDGcode(),atrack.GetKineticEnergy()/GeV,vx.z/m,pos.z/m,mom.z/GeV)
 
# ------------------------------------------------------------------

PostUserTrackingAction()

g4Ex_gap.MyTrackingActionD.PostUserTrackingAction	(		self,
			atrack
	)

Definition at line 310 of file g4Ex_gap.py.

 def PostUserTrackingAction(self,atrack):
   if particleGun and atrack.GetTrackID() == 1:
     # record dead of 400GeV proton
     if atrack.GetTrackStatus() != atrack.GetTrackStatus().fAlive :
      part         = atrack.GetDynamicParticle()
      pid          = part.GetPDGcode()
      vx           = atrack.GetVertexPosition()
      mom  = atrack.GetMomentum()
      ekin = atrack.GetKineticEnergy()/GeV
      pos  = atrack.GetPosition()
      w = atrack.GetWeight()
      parentid = int(w)/100000-10000
      pythiaid = int(w)%100000-10000
      h['ntuple'].Fill(float(pid), float(mom.x/GeV),float(mom.y/GeV),float(mom.z/GeV),\
                   float(pos.x/m),float(pos.y/m),float(pos.z/m),\
                   float(vx.x/m),float(vx.y/m),float(vx.z/m),pythiaid,parentid)
       

PreUserTrackingAction()

g4Ex_gap.MyTrackingActionD.PreUserTrackingAction	(		self,
			atrack
	)

Definition at line 327 of file g4Ex_gap.py.

 def PreUserTrackingAction(self,atrack):
   global trackHistory
# need to be careful with energy cut, anti-protons and neutrons can always produce pions and kaons
   part         = atrack.GetDynamicParticle()
   pid          = part.GetPDGcode()
   # if debug: self.myPrintout(atrack)
   moid = atrack.GetParentID()
   trackHistory[atrack.GetTrackID()]=[pid,moid]
   tmoid = str(moid)
   if moid in trackHistory: 
         mo = pdg.GetParticle(trackHistory[moid][0])
         if mo : tmoid = mo.GetName()
   tid = str(pid)
   if pdg.GetParticle(pid): tid = pdg.GetParticle(pid).GetName()
   mom = atrack.GetMomentum()
   p = ROOT.TMath.Sqrt(mom.x*mom.x+mom.y*mom.y+mom.z*mom.z)
   if debug and abs(pid)==13: print('track',atrack.GetTrackID(),tid,tmoid,moid,atrack.GetKineticEnergy()/MeV,p/MeV)
   if pid==12:
      if moid in trackHistory: 
        gmoid = trackHistory[moid][1]
        if gmoid in trackHistory:
          tgmoid = str(gmoid)
          mo = pdg.GetParticle(trackHistory[gmoid][0])
          if mo : tgmoid = mo.GetName()
          print('   <--',gmoid,tgmoid)
          gmoid = trackHistory[gmoid][1]
          if gmoid in trackHistory:
           tgmoid = str(gmoid)
           mo = pdg.GetParticle(trackHistory[gmoid][0])
           if mo : tgmoid = mo.GetName()
           print('      <--',gmoid,tgmoid)
 
   qed          = pid in qedlist  # use cut only for photons, electrons, protons and neutrons
   if atrack.GetKineticEnergy()/GeV < ecut and (qed or allPart): 
      G4TrackingManager().SetStoreTrajectory(False) 
      atrack.SetTrackStatus(atrack.GetTrackStatus().fStopAndKill)
 

---

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

• muonShieldOptimization/g4Ex_gap.py