VBF.VBF Class Reference

Inheritance diagram for VBF.VBF:

Public Member Functions
def	__init__
def	__str__
def	calcP4
def	findCentralJets

Public Attributes
	centralJets
	cjvPtCut
	deta
	dijetphi
	dijetpt
	dphi
	dphidijethiggs
	higgsp4
	jets
	leadJets
	met
	mjj
	mva
	self.ptvis = diLepton.pt() new VBF MVA, based on 4 variables More...
	otherJets
	ptvis
	vbfMvaCalc
	visjeteta

Detailed Description

Computes and holds VBF quantities

Definition at line 5 of file VBF.py.

Constructor & Destructor Documentation

def VBF.VBF.__init__	(		self,
			jets,
			diLepton,
			vbfMvaCalc,
			cjvPtCut
	)

jets: jets cleaned from the diLepton legs.
diLepton: the di-tau, for example. Necessary to compute input variables for MVA selection

Definition at line 7 of file VBF.py.

    def __init__(self, jets, diLepton, vbfMvaCalc, cjvPtCut):
        '''jets: jets cleaned from the diLepton legs.
        diLepton: the di-tau, for example. Necessary to compute input variables for MVA selection
        '''
        self.cjvPtCut = cjvPtCut
        self.vbfMvaCalc = vbfMvaCalc
        self.jets = jets
        # the MET is taken from the di-lepton, because it can depend on it
        # e.g. recoil corrections, mva met
        self.met = diLepton.met()
        self.leadJets = jets[:2] # taking first 2 jets (leading ones)
        self.otherJets = jets[2:]
        self.centralJets = self.findCentralJets( self.leadJets, self.otherJets )

        # delta eta
        self.deta = self.leadJets[0].eta() - self.leadJets[1].eta()

        # below, the variables for the MVA selection
        # delta phi
        self.dphi = deltaPhi(self.leadJets[0].phi(), self.leadJets[1].phi())
        dijetp4 = self.leadJets[0].p4() + self.leadJets[1].p4()
        # mass of the di-jet system
        self.mjj = dijetp4.M()
        # pt of di-jet system
        self.dijetpt = dijetp4.pt()
        # phi of di-jet system
        self.dijetphi = dijetp4.phi()
        # higgs momentum (defined as the di-lepton momentum + the met momentum)
        # don't access longitudinal quantities!
        self.higgsp4 = diLepton.p4() + self.met.p4()
        # delta phi between dijet system and higgs system
        self.dphidijethiggs = deltaPhi( self.dijetphi, self.higgsp4.phi() )
        # ? 
        visDiLepton = diLepton.leg1 ().p4 () + diLepton.leg2 ().p4 ()
        self.visjeteta = min (
            abs (self.leadJets[0].eta () - visDiLepton.eta ()), 
            abs (self.leadJets[1].eta () - visDiLepton.eta ()))
        # visible higgs pt = di-lepton pt
        self.ptvis = visDiLepton.pt()
        ## self.ptvis = diLepton.pt()
        # new VBF MVA, based on 4 variables
        if self.vbfMvaCalc is not None:
            self.mva = self.vbfMvaCalc.val( self.mjj,
                                            abs(self.deta),
                                            self.visjeteta,
                                            self.ptvis )
        else:
            self.mva = -99.

#  double mjj      , // the invariant mass of the two tag jets
#  double dEta     , // the pseudorapidity difference between the two tag jets
#  double dPhi     , // the phi difference between the two tag jets
#  double ditau_pt , // the vector sum of the pT of the tau + electron/muon + MET
#  double dijet_pt , // the vector sum of the pT of the two tag jets
#  double dPhi_hj  , // the phi difference between the di-tau vector and the di-jet vector
#  double C1       , // the pseudorapidity difference between the *visible* di-tau vector and the closest tag jet
#  double C2         // the *visible* pT of the di-tau