Functions
def	allDaughters (particle, daughters, rank)

def	bosonToX (particles, bosonType, xType)

def	findStatus1Leptons (particle)

def	isNotFromHadronicShower (l)

def	isNotHadronicId (pdgId, includeSMLeptons=True)

def	isPromptLepton (lepton, beforeFSR, includeMotherless=True, includeTauDecays=False)

def	lastGenCopy (gp)

def	realGenDaughters (gp, excludeRadiation=True)

def	realGenMothers (gp)

Function Documentation

def genutils.allDaughters	(	particle,
		daughters,
		rank
	)

Fills daughters with all the daughters of particle.
Recursive function.

Definition at line 19 of file genutils.py.

Referenced by bosonToX().

 def allDaughters(particle, daughters, rank ):
     '''Fills daughters with all the daughters of particle.
     Recursive function.'''
     rank += 1 
     for i in range( particle.numberOfDaughters() ):
         dau = GenParticle(particle.daughter(i))
         dau.rank = rank
         daughters.append( dau )
         daughters = allDaughters( dau, daughters, rank )
     return daughters
 
 

def genutils.bosonToX	(	particles,
		bosonType,
		xType
	)

Definition at line 31 of file genutils.py.

References funct.abs(), allDaughters(), and ALCARECOTkAlBeamHalo_cff.filter.

 def bosonToX(particles, bosonType, xType):
     bosons = filter(lambda x: x.status()==3 and x.pdgId()==bosonType, particles)
     daughters = []
     if len(bosons)==0:
         return [], False
     boson = bosons[0]
     daus = []
     allDaughters( boson, daus, 0)
     xDaus = filter(lambda x: x.status()==3 and abs(x.pdgId())==xType, daus)
     # print printOut(xDaus)
     return xDaus, True 
 

def genutils.findStatus1Leptons ( particle )

Returns status 1 e and mu among the particle daughters

Definition at line 4 of file genutils.py.

References funct.abs().

 def findStatus1Leptons(particle):
     '''Returns status 1 e and mu among the particle daughters'''
     leptons = []
     for i in range( particle.numberOfDaughters() ):
         dau = particle.daughter(i)
         if dau.status() == 1:
             if abs(dau.pdgId())==11 or abs(dau.pdgId())==13:
                 leptons.append( dau )
             else:
                 continue
         else:
             leptons = findStatus1Leptons( dau, leptons )
     return leptons
 
 

def genutils.isNotFromHadronicShower ( l )

Definition at line 65 of file genutils.py.

References funct.abs().

Referenced by gen.GeneratorAnalyzer.GeneratorAnalyzer.makeMCInfo().

 def isNotFromHadronicShower(l):
     for x in xrange(l.numberOfMothers()):
         mom = l.mother(x)
         if mom.status() > 2: return True
         id = abs(mom.pdgId())
         if id > 1000000: return True
         if id > 100: return False
         if id <   6: return False
         if id == 21: return False
         if id in [11,12,13,14,15,16]: 
             if l.status() > 2: return True
             return isNotFromHadronicShower(mom)
         if id >= 22 and id <= 39: return True
     return True
 

def genutils.isNotHadronicId	(	pdgId,
		includeSMLeptons = `True`
	)

Definition at line 43 of file genutils.py.

References funct.abs().

Referenced by isPromptLepton().

 def isNotHadronicId(pdgId,includeSMLeptons=True):
     if abs(pdgId) in [11,12,13,14,15,16]:
         return includeSMLeptons
     i = (abs(pdgId) % 1000)
     return i > 10 and i != 21 and i < 100
 

def genutils.isPromptLepton	(	lepton,
		beforeFSR,
		includeMotherless = `True`,
		includeTauDecays = `False`
	)

Definition at line 49 of file genutils.py.

References funct.abs(), and isNotHadronicId().

Referenced by objects.LeptonAnalyzer.LeptonAnalyzer.matchAnyLeptons().

 def isPromptLepton(lepton, beforeFSR, includeMotherless=True, includeTauDecays=False):
     if abs(lepton.pdgId()) not in [11,13,15]:
         return False
     if lepton.numberOfMothers() == 0:
         return includeMotherless;
     mom = lepton.mother()
     if mom.pdgId() == lepton.pdgId():
         if beforeFSR: return False
         return isPromptLepton(mom, beforeFSR, includeMotherless, includeTauDecays)
     elif abs(mom.pdgId()) == 15:
         if not includeTauDecays: return False
         return isPromptLepton(mom, beforeFSR, includeMotherless, includeTauDecays)
     else:
         return isNotHadronicId(mom.pdgId(), includeSMLeptons=False)
 
 

def genutils.lastGenCopy ( gp )

Definition at line 111 of file genutils.py.

References cosmictrackingParticleSelector_cfi.pdgId.

 def lastGenCopy(gp):
     me = gp.pdgId();
     for i in xrange(gp.numberOfDaughters()):
         if gp.daughter(i).pdgId() == me:
             return False
     return True
 
 

def genutils.realGenDaughters	(	gp,
		excludeRadiation = `True`
	)

Get the daughters of a particle, going through radiative X -> X' + a
   decays, either including or excluding the radiation among the daughters
   e.g. for  
              X -> X' + a, X' -> b c 
       realGenDaughters(X, excludeRadiation=True)  = { b, c }
       realGenDaughters(X, excludeRadiation=False) = { a, b, c }

Definition at line 80 of file genutils.py.

Referenced by gen.GeneratorAnalyzer.GeneratorAnalyzer.makeMCInfo().

 def realGenDaughters(gp,excludeRadiation=True):
     """Get the daughters of a particle, going through radiative X -> X' + a
        decays, either including or excluding the radiation among the daughters
        e.g. for  
                   X -> X' + a, X' -> b c 
            realGenDaughters(X, excludeRadiation=True)  = { b, c }
            realGenDaughters(X, excludeRadiation=False) = { a, b, c }"""
     ret = []
     for i in xrange(gp.numberOfDaughters()):
         dau = gp.daughter(i)
         if dau.pdgId() == gp.pdgId():
             if excludeRadiation:
                 return realGenDaughters(dau)
             else:
                 ret += realGenDaughters(dau)
         else:
             ret.append(dau)
     return ret
 

def genutils.realGenMothers ( gp )

Get the mothers of a particle X going through intermediate X -> X' chains.
   e.g. if Y -> X, X -> X' realGenMothers(X') = Y

Definition at line 99 of file genutils.py.

Referenced by gen.GeneratorAnalyzer.GeneratorAnalyzer.makeMCInfo().

 def realGenMothers(gp):
     """Get the mothers of a particle X going through intermediate X -> X' chains.
        e.g. if Y -> X, X -> X' realGenMothers(X') = Y"""
     ret = []
     for i in xrange(gp.numberOfMothers()):
         mom = gp.mother(i)
         if mom.pdgId() == gp.pdgId():
             ret += realGenMothers(mom)
         else:
             ret.append(mom)
     return ret
 

Functions

Function Documentation