gen.LHEAnalyzer.LHEAnalyzer Class Reference

Inheritance diagram for gen.LHEAnalyzer.LHEAnalyzer:

Public Member Functions
def	__init__ (self, cfg_ana, cfg_comp, looperName)
def	beginLoop (self, setup)
def	declareHandles (self)
def	process (self, event)

Public Attributes
	lheh

Detailed Description

Definition at line 8 of file LHEAnalyzer.py.

Constructor & Destructor Documentation

def gen.LHEAnalyzer.LHEAnalyzer.__init__	(		self,
			cfg_ana,
			cfg_comp,
			looperName
	)

Definition at line 10 of file LHEAnalyzer.py.

    def __init__(self, cfg_ana, cfg_comp, looperName ):
        super(LHEAnalyzer,self).__init__(cfg_ana,cfg_comp,looperName)
        self.lheh=Handle('LHEEventProduct')

Member Function Documentation

def gen.LHEAnalyzer.LHEAnalyzer.beginLoop	(		self,
			setup
	)

Definition at line 18 of file LHEAnalyzer.py.

    def beginLoop(self, setup):
        super(LHEAnalyzer,self).beginLoop(setup)

def gen.LHEAnalyzer.LHEAnalyzer.declareHandles

(

self

)

Definition at line 14 of file LHEAnalyzer.py.

    def declareHandles(self):
        super(LHEAnalyzer, self).declareHandles()
#        self.mchandles['lhestuff'] = AutoHandle( 'externalLHEProducer','LHEEventProduct')

def gen.LHEAnalyzer.LHEAnalyzer.process	(		self,
			event
	)

Definition at line 21 of file LHEAnalyzer.py.

References funct.abs(), gen.LHEAnalyzer.LHEAnalyzer.lheh, SiStripPI.max, core.Analyzer.Analyzer.readCollections(), and mathSSE.sqrt().

    def process(self, event):

        # if not MC, nothing to do
        if not self.cfg_comp.isMC: 
            return True
        event.lheHT=0
        event.lheHTIncoming=0 #restrict HT computation to particles that have status<0 mothers
        event.lheNj=0
        event.lheNb=0
        event.lheNc=0
        event.lheNl=0
        event.lheNg=0
        event.lheV_pt = 0
        try:
          event.input.getByLabel( 'externalLHEProducer',self.lheh)
        except :
          return True
        if not  self.lheh.isValid() :
            return True
        self.readCollections( event.input )
        hepeup=self.lheh.product().hepeup()
        pup=hepeup.PUP
        l=None
        lBar=None
        nu=None
        nuBar=None 
        for i in range(0,len(pup)):
          id=hepeup.IDUP[i]
          status = hepeup.ISTUP[i]
          idabs=abs(id)

          mothIdx = max(hepeup.MOTHUP[i][0]-1,0) #first and last mother as pair; first entry has index 1 in LHE; incoming particles return motherindex 0
          mothIdxTwo = max(hepeup.MOTHUP[i][1]-1,0) 
          
          mothStatus  = hepeup.ISTUP[mothIdx] 
          mothStatusTwo  = hepeup.ISTUP[mothIdxTwo] 

          hasIncomingAsMother = mothStatus<0 or mothStatusTwo<0
          
          if status == 1 and ( ( idabs == 21 ) or (idabs > 0 and idabs < 7) ) : # gluons and quarks
              pt = sqrt( pup[i][0]**2 + pup[i][1]**2 ) # first entry is px, second py
              event.lheHT += pt
              if hasIncomingAsMother: event.lheHTIncoming += pt
              event.lheNj +=1
              if idabs==5:
                event.lheNb += 1
              if idabs==4:
                event.lheNc += 1
              if idabs in [1,2,3]:
                event.lheNl += 1
              if idabs==21:
                event.lheNg += 1
          if idabs in [12,14,16] :  
              if id > 0 :
                nu = i
              else :
                nuBar = i
          if idabs in [11,13,15] :  
              if id > 0 :
                l = i
              else :
                lBar = i
          v=None
          if l and lBar : #Z to LL
              v=(l,lBar)
          elif l and nuBar : #W 
              v=(l,nuBar)
          elif lBar and nu : #W 
              v=(nu,lBar)
          elif nu and nuBar : #Z to nn 
              v=(nu,nuBar)
          if v :
            event.lheV_pt = sqrt( (pup[v[0]][0]+pup[v[1]][0])**2 +  (pup[v[0]][1]+pup[v[1]][1])**2 )

        return True

setattr(LHEAnalyzer,"defaultConfig",

Member Data Documentation

gen.LHEAnalyzer.LHEAnalyzer.lheh

Definition at line 12 of file LHEAnalyzer.py.

Referenced by gen.LHEAnalyzer.LHEAnalyzer.process().