genutils Namespace Reference

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 common_cff.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

---

Generated for CMSSW Reference Manual by  1.8.11