Inheritance diagram for ParticleGun_SamplingFraction.MyParticleSampler:

Collaboration diagram for ParticleGun_SamplingFraction.MyParticleSampler:

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

def	shoot (self)

Public Attributes
	pid

	momentum

	eta1

	eta2

	mass

	bec

	radius

	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()

def 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()

def 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
  