MT2Analyzer.MT2Analyzer Class Reference

Inheritance diagram for MT2Analyzer.MT2Analyzer:

Public Member Functions
def	__init__ (self, cfg_ana, cfg_comp, looperName)
def	beginLoop (self, setup)
def	computeMT2 (self, visaVec, visbVec, metVec)
def	declareHandles (self)
def	getMT2AKT (self, event, TMPobjects40jc, met, collectionPostFix, postFix)
def	getMT2Hemi (self, event, TMPobjects40jc, met, collectionPostFix, postFix)
def	makeMT2 (self, event)

Public Attributes
	jetPt
	met

Static Public Attributes
	groupingViaAKt
	groupingViaKt
	hemisphereViaAKt
	hemisphereViaKt
	mt2ViaKt_had
	get hemispheres via AntiKT -1 antikt, 1 kt, 0 CA More...
	mt2ViaKt_Xj_had
	objects
	get hemispheres via AntiKT -1 antikt, 1 kt, 0 CA More...

Detailed Description

Definition at line 33 of file MT2Analyzer.py.

Constructor & Destructor Documentation

__init__()

def MT2Analyzer.MT2Analyzer.__init__	(		self,
			cfg_ana,
			cfg_comp,
			looperName
	)

Definition at line 34 of file MT2Analyzer.py.

    def __init__(self, cfg_ana, cfg_comp, looperName ):
        super(MT2Analyzer,self).__init__(cfg_ana,cfg_comp,looperName) 
        self.jetPt = cfg_ana.jetPt
 

Member Function Documentation

beginLoop()

def MT2Analyzer.MT2Analyzer.beginLoop	(		self,
			setup
	)

Definition at line 44 of file MT2Analyzer.py.

    def beginLoop(self, setup):
        super(MT2Analyzer,self).beginLoop(setup)
        self.counters.addCounter('pairs')
        count = self.counters.counter('pairs')
        count.register('all events')
 

References CSCSPEvent.counters(), cms::cuda::AtomicPairCounter::Atomic2.counters, NMaxPerLumi.counters, analyzer.Analyzer.counters, HistogramManager.counters, counter.Counters.counters, and cscdqm::ChamberCounterKeyType.counters.

computeMT2()

def MT2Analyzer.MT2Analyzer.computeMT2	(		self,
			visaVec,
			visbVec,
			metVec
	)

Definition at line 50 of file MT2Analyzer.py.

    def computeMT2(self, visaVec, visbVec, metVec):
        
        import array
        import numpy
 
        metVector = array.array('d',[0.,metVec.px(), metVec.py()])
        visaVector = array.array('d',[0.,visaVec.px(), visaVec.py()])
        visbVector = array.array('d',[0.,visbVec.px(), visbVec.py()])
 
        davismt2.set_momenta(visaVector,visbVector,metVector);
        davismt2.set_mn(0);
 
        return davismt2.get_mt2()
 

Referenced by MT2Analyzer.MT2Analyzer.getMT2Hemi().

declareHandles()

def MT2Analyzer.MT2Analyzer.declareHandles

(

self

)

Definition at line 38 of file MT2Analyzer.py.

    def declareHandles(self):
        super(MT2Analyzer, self).declareHandles()
       #genJets                                                                                                                                                                     
        self.handles['genJets'] = AutoHandle( 'slimmedGenJets','std::vector<reco::GenJet>')
        self.handles['met'] = AutoHandle( self.cfg_ana.metCollection, 'std::vector<pat::MET>' )
 

References analyzer.Analyzer.cfg_ana, core.AutoHandle.AutoHandle.handles, core.Analyzer.Analyzer.handles, and HTTP.RequestManager.handles.

getMT2AKT()

def MT2Analyzer.MT2Analyzer.getMT2AKT	(		self,
			event,
			TMPobjects40jc,
			met,
			collectionPostFix,
			postFix
	)

Definition at line 64 of file MT2Analyzer.py.

    def getMT2AKT(self, event, TMPobjects40jc, met , collectionPostFix, postFix):
 

getMT2Hemi()

def MT2Analyzer.MT2Analyzer.getMT2Hemi	(		self,
			event,
			TMPobjects40jc,
			met,
			collectionPostFix,
			postFix
	)

Definition at line 92 of file MT2Analyzer.py.

    def getMT2Hemi(self, event, TMPobjects40jc, met, collectionPostFix, postFix):
 
        if len(TMPobjects40jc)>=2:
 
            pxvec  = ROOT.std.vector(float)()
            pyvec  = ROOT.std.vector(float)()
            pzvec  = ROOT.std.vector(float)()
            Evec  = ROOT.std.vector(float)()
            grouping  = ROOT.std.vector(int)()
            
            for jet in TMPobjects40jc:
                pxvec.push_back(jet.px())
                pyvec.push_back(jet.py())
                pzvec.push_back(jet.pz())
                Evec.push_back(jet.energy())
 
            hemisphere = Hemisphere(pxvec, pyvec, pzvec, Evec, 2, 3)
            grouping=hemisphere.getGrouping()
 
            pseudoJet1px = 0 
            pseudoJet1py = 0 
            pseudoJet1pz = 0
            pseudoJet1energy = 0 
            multPSJ1 = 0
 
            pseudoJet2px = 0 
            pseudoJet2py = 0 
            pseudoJet2pz = 0
            pseudoJet2energy = 0 
            multPSJ2 = 0
                
            for index in range(0, len(pxvec)):
                if(grouping[index]==1):
                    pseudoJet1px += pxvec[index]
                    pseudoJet1py += pyvec[index]
                    pseudoJet1pz += pzvec[index]
                    pseudoJet1energy += Evec[index]
                    multPSJ1 += 1
                if(grouping[index]==2):
                    pseudoJet2px += pxvec[index]
                    pseudoJet2py += pyvec[index]
                    pseudoJet2pz += pzvec[index]
                    pseudoJet2energy += Evec[index]                    
                    multPSJ2 += 1
 
            pseudoJet1pt2 = pseudoJet1px*pseudoJet1px + pseudoJet1py*pseudoJet1py
            pseudoJet2pt2 = pseudoJet2px*pseudoJet2px + pseudoJet2py*pseudoJet2py
 
            if pseudoJet1pt2 >= pseudoJet2pt2:
                setattr(event, "pseudoJet1"+collectionPostFix+postFix, ROOT.reco.Particle.LorentzVector( pseudoJet1px, pseudoJet1py, pseudoJet1pz, pseudoJet1energy ))
                setattr(event, "pseudoJet2"+collectionPostFix+postFix, ROOT.reco.Particle.LorentzVector( pseudoJet2px, pseudoJet2py, pseudoJet2pz, pseudoJet2energy ))
                setattr(event, "multPseudoJet1"+collectionPostFix+postFix, multPSJ1 )
                setattr(event, "multPseudoJet2"+collectionPostFix+postFix, multPSJ2 )
            else:
                setattr(event, "pseudoJet2"+collectionPostFix+postFix, ROOT.reco.Particle.LorentzVector( pseudoJet1px, pseudoJet1py, pseudoJet1pz, pseudoJet1energy ))
                setattr(event, "pseudoJet1"+collectionPostFix+postFix, ROOT.reco.Particle.LorentzVector( pseudoJet2px, pseudoJet2py, pseudoJet2pz, pseudoJet2energy ))
                setattr(event, "multPseudoJet1"+collectionPostFix+postFix, multPSJ2 )
                setattr(event, "multPseudoJet2"+collectionPostFix+postFix, multPSJ1 )
 
            setattr(event, "mt2"+collectionPostFix+postFix, self.computeMT2(getattr(event,'pseudoJet1'+collectionPostFix+postFix), getattr(event,'pseudoJet2'+collectionPostFix+postFix), met) )
            return self.computeMT2(getattr(event,'pseudoJet1'+collectionPostFix+postFix), getattr(event,'pseudoJet2'+collectionPostFix+postFix), met)
 
 

References MT2Analyzer.MT2Analyzer.computeMT2(), and FastTimerService_cff.range.

makeMT2()

def MT2Analyzer.MT2Analyzer.makeMT2	(		self,
			event
	)

Definition at line 155 of file MT2Analyzer.py.

    def makeMT2(self, event):
#        print '==> INSIDE THE PRINT MT2'
#        print 'MET=',event.met.pt()
 
        self.met = ROOT.pat.MET(self.handles['met'].product()[0])
 
        import array
        import numpy
 
        objects40jc = [ j for j in event.cleanJets if j.pt() > 40 and abs(j.eta())<2.5 ]
        objectsXjc = [ j for j in event.cleanJets if j.pt() > self.jetPt and abs(j.eta())<2.5 ]
 
        setattr(event, "mt2ViaKt"+self.cfg_ana.collectionPostFix+"had", -999)
        setattr(event, "mt2ViaKt"+self.cfg_ana.collectionPostFix+"_Xj_had", -999)
        setattr(event, "pseudoViaKtJet1"+self.cfg_ana.collectionPostFix+"_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaKtJet2"+self.cfg_ana.collectionPostFix+"_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaKtJet1"+self.cfg_ana.collectionPostFix+"_Xj_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaKtJet2"+self.cfg_ana.collectionPostFix+"_Xj_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
 
        setattr(event, "mt2ViaAKt"+self.cfg_ana.collectionPostFix+"had", -999)
        setattr(event, "mt2ViaAKt"+self.cfg_ana.collectionPostFix+"_Xj_had", -999)
        setattr(event, "pseudoViaAKtJet1"+self.cfg_ana.collectionPostFix+"_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaAKtJet2"+self.cfg_ana.collectionPostFix+"_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaAKtJet1"+self.cfg_ana.collectionPostFix+"_Xj_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
        setattr(event, "pseudoViaAKtJet2"+self.cfg_ana.collectionPostFix+"_Xj_had", ROOT.reco.Particle.LorentzVector( 0, 0, 0, 0 ))
 

Member Data Documentation

groupingViaAKt

MT2Analyzer.MT2Analyzer.groupingViaAKt

static

Definition at line 84 of file MT2Analyzer.py.

groupingViaKt

MT2Analyzer.MT2Analyzer.groupingViaKt

static

Definition at line 74 of file MT2Analyzer.py.

hemisphereViaAKt

MT2Analyzer.MT2Analyzer.hemisphereViaAKt

static

Definition at line 83 of file MT2Analyzer.py.

hemisphereViaKt

MT2Analyzer.MT2Analyzer.hemisphereViaKt

static

Definition at line 73 of file MT2Analyzer.py.

jetPt

MT2Analyzer.MT2Analyzer.jetPt

Definition at line 36 of file MT2Analyzer.py.

met

MT2Analyzer.MT2Analyzer.met

Definition at line 159 of file MT2Analyzer.py.

Referenced by objects.METAnalyzer.METAnalyzer.applyDeltaMet().

mt2ViaKt_had

MT2Analyzer.MT2Analyzer.mt2ViaKt_had

static

get hemispheres via AntiKT -1 antikt, 1 kt, 0 CA

Definition at line 184 of file MT2Analyzer.py.

mt2ViaKt_Xj_had

MT2Analyzer.MT2Analyzer.mt2ViaKt_Xj_had

static

Definition at line 188 of file MT2Analyzer.py.

objects

MT2Analyzer.MT2Analyzer.objects

static

get hemispheres via AntiKT -1 antikt, 1 kt, 0 CA

Definition at line 69 of file MT2Analyzer.py.