Inheritance diagram for ParticleGun_FastCalo_ChargeFlip_Config.MyParticleSampler:

Collaboration diagram for ParticleGun_FastCalo_ChargeFlip_Config.MyParticleSampler:

Public Member Functions
	__init__ (self, energy, eta, pid, shift_z=0)
	shoot (self)

Public Attributes
	pid = pid
	shift_z = shift_z
	mom1 = PG.EEtaMPhiSampler(energy=energy,eta=eta,mass=mass)

Detailed Description

Definition at line 10 of file ParticleGun_FastCalo_ChargeFlip_Config.py.

Constructor & Destructor Documentation

◆ init()

ParticleGun_FastCalo_ChargeFlip_Config.MyParticleSampler.__init__	(	self,
		energy,
		eta,
		pid,
		shift_z = 0 )

Definition at line 11 of file ParticleGun_FastCalo_ChargeFlip_Config.py.

    def __init__(self,energy,eta,pid,shift_z=0):
        self.pid = pid
        self.shift_z = shift_z
        pdg_table = ROOT.TDatabasePDG.Instance() #Gives values in GeV
        mass = pdg_table.GetParticle(self.pid()).Mass()*1000.
        self.mom1 = PG.EEtaMPhiSampler(energy=energy,eta=eta,mass=mass)
 

Member Function Documentation

◆ shoot()

ParticleGun_FastCalo_ChargeFlip_Config.MyParticleSampler.shoot ( self )

Definition at line 18 of file ParticleGun_FastCalo_ChargeFlip_Config.py.

    def shoot(self):
        pid = self.pid()
        
        shift_z = self.shift_z
 
        mom = self.mom1.shoot()
        pos_temp = mom.Vect().Unit()
        
        # Would it hit the barrel, or the endcap?
        if abs(pos_temp.Z())/3550.<pos_temp.Perp()/1148.: # Hit the barrel!
            pos_temp *= 1148./pos_temp.Perp()
        else: # Hit the endcap!
            pos_temp *= 3550./abs(pos_temp.Z())
 
        # Shift position of vector in the Z direction
        pos_temp_2 = ROOT.TVector3()
        pos_temp_2.SetXYZ(pos_temp.X(), pos_temp.Y(), pos_temp.Z()+shift_z)
        pos_temp_2 *= 1. / pos_temp_2.Mag(); # reduce magnitude of vector
 
        # recalculate; Would it hit the barrel, or the endcap?
        if abs(pos_temp_2.Z())/3550.<pos_temp_2.Perp()/1148.:
            pos_temp_2 *= 1148./pos_temp_2.Perp()
        else:
            pos_temp_2 *= 3550./abs(pos_temp_2.Z())
 
        pos = ROOT.TLorentzVector(pos_temp_2.X(),pos_temp_2.Y(),pos_temp_2.Z(), pos_temp_2.Mag())
        
        #print "pid ",pid
        
        return [ PG.SampledParticle( pid , mom , pos ) ]
 