ParticleGun_egammaET.egammaETSampler Class Reference

Inheritance diagram for ParticleGun_egammaET.egammaETSampler:

Collaboration diagram for ParticleGun_egammaET.egammaETSampler:

Public Member Functions
	__init__ (self, pid, eta=[-2.5, 2.5], phi=[0, PG.TWOPI], mu1=0.5, kT1=0.1, mu2=200, kT2=20, y0=0.005, PtMin=0, PtMax=3e3, nBins=None)
	pt (self)

Public Attributes
	m = PG.MASSES[abs(pid)]
	eta = eta
	phi = phi
	hist = PG.TH1(etSpectrumFullRange)

Detailed Description

Definition at line 18 of file ParticleGun_egammaET.py.

Constructor & Destructor Documentation

__init__()

ParticleGun_egammaET.egammaETSampler.__init__	(		self,
			pid,
			eta = [-2.5, 2.5],
			phi = [0, PG.TWOPI],
			mu1 = 0.5,
			kT1 = 0.1,
			mu2 = 200,
			kT2 = 20,
			y0 = 0.005,
			PtMin = 0,
			PtMax = 3e3,
			nBins = None )

Parameters for the MVA-shaped spectrum : higher density in the < 100 GeV range
PtMin = 0 # minimum Pt
PtMax = 3000 # maximum Pt (3 TeV)
nBins # number of bins (one every 100 MeV by default)
mu1 = 0.5       # mu1,kT1 : smooth but steep ramp-up from 0 to 1 GeV (requested by TauCP)
kT1 = 0.1
mu2 = 200       # mu2,kT2 : smooth, slow ramp-down in the 100-300 GeV range
kT2 = 20
y0 = 0.005      # y0 : baseline for low-density at high ET up to PtMax

Definition at line 20 of file ParticleGun_egammaET.py.

    mu1 = 0.5, kT1 = 0.1, mu2 = 200, kT2 = 20, y0 = 0.005, PtMin = 0 , PtMax = 3e3, nBins=None):
    """ 
    Parameters for the MVA-shaped spectrum : higher density in the < 100 GeV range
    PtMin = 0 # minimum Pt
    PtMax = 3000 # maximum Pt (3 TeV)
    nBins # number of bins (one every 100 MeV by default)
    mu1 = 0.5       # mu1,kT1 : smooth but steep ramp-up from 0 to 1 GeV (requested by TauCP)
    kT1 = 0.1
    mu2 = 200       # mu2,kT2 : smooth, slow ramp-down in the 100-300 GeV range
    kT2 = 20
    y0 = 0.005      # y0 : baseline for low-density at high ET up to PtMax
    """
    self.m = PG.MASSES[abs(pid)]
    self.eta = eta
    self.phi = phi
    
    # Create and fill a very fine-grained histogram 
    from ROOT import TH1D
    etSpectrumFullRange = TH1D("ETSpectrumFullRange", 
                               "Reference ET spectrum for egamma MVA calib",
                               int(nBins or (PtMax - PtMin)*10), PtMin , PtMax)
    for i in xrange(etSpectrumFullRange.GetNbinsX()):
         x = etSpectrumFullRange.GetBinCenter(i+1)
         y1 = dbnFermiDirac(x,mu1,kT1)
         y2 = dbnFermiDirac(x,mu2,kT2)
         y = y0 - y1 + y2
         etSpectrumFullRange.SetBinContent(i+1,y)
    self.hist = PG.TH1(etSpectrumFullRange) #< wrap *after* populating
  

Member Function Documentation

pt()

ParticleGun_egammaET.egammaETSampler.pt

(

self

)

Definition at line 50 of file ParticleGun_egammaET.py.

  def pt(self):
    return self.hist.GetRandom() * GeV

Member Data Documentation

eta

ParticleGun_egammaET.egammaETSampler.eta = eta

Definition at line 34 of file ParticleGun_egammaET.py.

hist

ParticleGun_egammaET.egammaETSampler.hist = PG.TH1(etSpectrumFullRange)

Definition at line 48 of file ParticleGun_egammaET.py.

m

ParticleGun_egammaET.egammaETSampler.m = PG.MASSES[abs(pid)]

Definition at line 33 of file ParticleGun_egammaET.py.

phi

ParticleGun_egammaET.egammaETSampler.phi = phi

Definition at line 35 of file ParticleGun_egammaET.py.

---

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

• ParticleGun_egammaET.py