g4Ex.MyGeneratorAction Class Reference

Inheritance diagram for g4Ex.MyGeneratorAction:

Collaboration diagram for g4Ex.MyGeneratorAction:

Public Member Functions
	GeneratePrimaries (self, anEvent)

Detailed Description

Definition at line 98 of file g4Ex.py.

Member Function Documentation

GeneratePrimaries()

g4Ex.MyGeneratorAction.GeneratePrimaries	(		self,
			anEvent
	)

Definition at line 100 of file g4Ex.py.

 def GeneratePrimaries(self,anEvent):
   global debug,nevTot
   t_0 = time.time()
   npart = 0
   while npart == 0:
    myPythia.GenerateEvent()      
    nevTot+=1
    myTimer['pythia']+=time.time()-t_0
# pythia interaction happens at 0,0,0
    #x = rnr.Uniform(-3.,3.) 
    #y = rnr.Uniform(-3.,3.)
    # leave this smearing for later
    pos    = G4ThreeVector(0*cm, 0*cm, -50*m)
    vertex = G4PrimaryVertex(pos,0.)
  # create new primaries and set them to the vertex
    particles = myPythia.GetListOfParticles()
    for p in particles:
       if p.GetStatusCode()!=1 : continue
       pid = p.GetPdgCode()
       if tauOnly and abs(pid) != 16: continue
       if pid in notWanted : continue 
       G4particle = G4PrimaryParticle( pid )
       v = TLorentzVector()
       p.Momentum(v)
       if v.E()*GeV  < ecut : continue
       G4particle.Set4Momentum( v.Px()*GeV,v.Py()*GeV,v.Pz()*GeV,v.E()*GeV )
       vertex.SetPrimary( G4particle )
# store mother ID
       mkey   =  p.GetMother(0)+1  
       mother =  myPythia.GetParticle(mkey)
       curPid =  p.GetPdgCode()      + 10000  # make it positive
       moPid  =  mother.GetPdgCode() + 10000
       w      =  curPid + moPid * 100000
       G4particle.SetWeight(w)
       npart += 1
   if tauOnly and debug: myPythia.EventListing()
   anEvent.AddPrimaryVertex( vertex )
   myTimer['geant4_conv']+=time.time()-t_0

---

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

• muonShieldOptimization/g4Ex.py