Inheritance diagram for ParticleGun_SamplingFraction.MyParticleSampler:

Collaboration diagram for ParticleGun_SamplingFraction.MyParticleSampler:

Public Member Functions
	__init__ (self, pid=11, momentum=50000., eta1=0., eta2=1.4, bec=0, radius=1500., z=3740.5)
	shoot (self)

Public Attributes
	pid = pid
	momentum = momentum
	eta1 = eta1
	eta2 = eta2
int	mass = pdg_table.GetParticle(self.pid()).Mass()*1000.
int	bec = bec
	radius = radius
	z = z

Detailed Description

Projective showers starting at entrance of calorimeter, flat in eta, constant energy

Definition at line 10 of file ParticleGun_SamplingFraction.py.

Constructor & Destructor Documentation

◆ init()

ParticleGun_SamplingFraction.MyParticleSampler.__init__	(	self,
		pid = 11,
		momentum = 50000.,
		eta1 = 0.,
		eta2 = 1.4,
		bec = 0,
		radius = 1500.,
		z = 3740.5 )

Definition at line 15 of file ParticleGun_SamplingFraction.py.

    def __init__(self,pid=11,momentum=50000.,eta1=0.,eta2=1.4,bec=0,radius=1500.,z=3740.5):
        self.pid = pid
        self.momentum = momentum
        self.eta1 = eta1
        self.eta2 = eta2
        pdg_table = ROOT.TDatabasePDG.Instance() #Gives values in GeV
        self.mass = pdg_table.GetParticle(self.pid()).Mass()*1000.
        self.bec=bec
        self.radius=radius
        self.z=z
 

Member Function Documentation

◆ shoot()

ParticleGun_SamplingFraction.MyParticleSampler.shoot ( self )

Definition at line 26 of file ParticleGun_SamplingFraction.py.

    def shoot(self):
       rtn=[]
       eta = random.uniform(self.eta1, self.eta2) 
       phi = random.uniform(0, math.tau)  # tau = 2 * pi
       v4 = ROOT.TLorentzVector()
       pt = self.momentum / math.cosh(eta)
       v4.SetPtEtaPhiM(pt, eta, phi, self.mass)
       if self.bec==0:
           radius= self.radius
           x=radius*math.cos(phi)
           y=radius*math.sin(phi)
           z=radius*math.sinh(eta)
       else:
           z=self.z
           radius=z/math.sinh(eta)
           x=radius*math.cos(phi)
           y=radius*math.sin(phi)
       t=math.sqrt(x*x+y*y+z*z)
       vp = ROOT.TLorentzVector(x,y,z,t)
       p = PG.SampledParticle(pid=self.pid(),mom=v4,pos=vp)
       #print "E,eta,phi,mass ",e,eta,phi,self.mass,"  position ",x,y,z," pid=",p.pid
       rtn.append(p)
       return rtn
 