g4Ex_gap.MyGeneratorAction Class Reference

Inheritance diagram for g4Ex_gap.MyGeneratorAction:

Collaboration diagram for g4Ex_gap.MyGeneratorAction:

Public Member Functions
	GeneratePrimaries (self, anEvent)

Static Public Attributes
	pos = G4ThreeVector(0*cm, 0*cm, -world_r/2.)

Detailed Description

Definition at line 191 of file g4Ex_gap.py.

Member Function Documentation

GeneratePrimaries()

g4Ex_gap.MyGeneratorAction.GeneratePrimaries	(		self,
			anEvent
	)

Definition at line 195 of file g4Ex_gap.py.

 def GeneratePrimaries(self,anEvent):
  global debug,nevTot,particleGun,fullTungsten
  t_0 = time.time()
  npart = 0
  if particleGun:
# to evaluate interaction position
   ztarget   = G4ThreeVector(0*cm, 0*cm,  -50.*m)
   vertex    = G4PrimaryVertex(ztarget,0.)
   G4particle = G4PrimaryParticle( 2212 )
   G4particle.Set4Momentum( 0,0,400*GeV,400.0011*GeV )
# store mother ID
   w      =  2212 + 0 * 100000
   G4particle.SetWeight(w)
   vertex.SetPrimary( G4particle )
   npart += 1
  else:
   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
  # create new primaries and set them to the vertex
    particles = myPythia.GetListOfParticles()
    if fullTungsten: z = rnr.Exp(10*cm) -50.*m # tungsten interaction length
    else: z = rnr.Exp(16*cm) -50.*m # Mo/H20/W average interaction length
    ztarget   = G4ThreeVector(0*cm, 0*cm,   z)
    vertex    = G4PrimaryVertex(ztarget,0.)
    v = ROOT.TLorentzVector()
    for p in particles:
       if p.GetStatusCode()!=1 : continue
       pid = p.GetPdgCode()
       mkey   =  p.GetMother(0)+1  
       mother =  myPythia.GetParticle(mkey)
       if JpsiMainly and abs(pid) != 13  : continue
       if tauOnly    and abs(pid) != 16 : continue
       if JpsiMainly and mother.GetPdgCode() != 443 : continue
       if tauOnly  : 
         mmkey   =  mother.GetMother(0)+1  
         mmother =  myPythia.GetParticle(mmkey)
         if abs(mother.GetPdgCode()) != 431 and abs(mmother.GetPdgCode()) != 431 : 
          myPythia.EventListing()
          continue
       qed = pid in qedlist  # use cut only for photons, leptons/neutrinos, protons and neutrons
       p.Momentum(v)
       Ekin = (v.E()-v.M())
       if Ekin*GeV < ecut and (qed or allPart): continue
       G4particle = G4PrimaryParticle( pid )
       G4particle.Set4Momentum( v.Px()*GeV,v.Py()*GeV,v.Pz()*GeV,v.E()*GeV )
# store mother ID
       curPid =  p.GetPdgCode()      + 10000  # make it positive
       moPid  =  mother.GetPdgCode() + 10000
       w      =  curPid + moPid * 100000
       G4particle.SetWeight(w)
       vertex.SetPrimary( G4particle )
       npart += 1
   # if debug: myPythia.EventListing()
  anEvent.AddPrimaryVertex( vertex )
  if debug and not particleGun: print('new event at ',ztarget.z/m)
  myTimer['geant4_conv']+=time.time()-t_0

Member Data Documentation

pos

g4Ex_gap.MyGeneratorAction.pos = G4ThreeVector(0*cm, 0*cm, -world_r/2.)

static

Definition at line 193 of file g4Ex_gap.py.

---

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

• muonShieldOptimization/g4Ex_gap.py