Beispiel #1
0
 def listOfCalculables(self,params) :
     _jet = params["objects"]["jet"]
     _muon = params["objects"]["muon"]
     _electron = params["objects"]["electron"]
     _photon = params["objects"]["photon"]
     _jetPtMin = params["objects"]["jetPtMin"]
     _etRatherThanPt = params["etRatherThanPt"]
     
     return calculables.zeroArgs() +\
            calculables.fromCollections("calculablesJet",[_jet]) +\
            calculables.fromCollections("calculablesMuon",[_muon]) +\
            calculables.fromCollections("calculablesElectron",[_electron]) +\
            calculables.fromCollections("calculablesPhoton",[_photon]) +\
            [ #calculables.xcJet( _jet,  gamma = _photon, gammaDR = 0.5, muon = _muon, muonDR = 0.5, electron = _electron, electronDR = 0.5),
              calculables.jetIndices( _jet, _jetPtMin, etaMax = 3.0, flagName = params["jetId"]),
              calculables.muonIndices( _muon, ptMin = 20, combinedRelIsoMax = 0.15),
              calculables.electronIndices( _electron, ptMin = 20, simpleEleID = "95", useCombinedIso = True),
              calculables.photonIndicesPat(  ptMin = 20, flagName = "photonIDLoosePat"),
              #calculables.indicesUnmatched(collection = _photon, xcjets = _jet, DR = 0.5),
              #calculables.indicesUnmatched(collection = _electron, xcjets = _jet, DR = 0.5)
              ] \
              + [ calculables.jetSumP4(_jet, mcScaleFactor = params["mcMhtScaleFactor"]),
                  calculables.deltaPhiStar(_jet, ptMin = 50.0),
                  calculables.jetSumP4Low(_jet, ptMin = 30.0),
                  calculables.deltaPseudoJet(_jet, _etRatherThanPt),
                  calculables.alphaT(_jet, _etRatherThanPt) ]
Beispiel #2
0
    def listOfCalculables(self,params) :
        _jet = params["objects"]["jet"]
        _muon = params["objects"]["muon"]
        _electron = params["objects"]["electron"]
        _photon = params["objects"]["photon"]
        _jetPtMin = params["objects"]["jetPtMin"]
        _etRatherThanPt = params["etRatherThanPt"]
        _met = params["objects"]["met"]
        _correctForMuons = not params["objects"]["muonsInJets"]

        lowPtThreshold = 30.0
        lowPtName = "lowPt"
        
        return calculables.zeroArgs() +\
               calculables.fromCollections("calculablesJet",[_jet]) +\
               calculables.fromCollections("calculablesMuon",[_muon]) +\
               calculables.fromCollections("calculablesElectron",[_electron]) +\
               calculables.fromCollections("calculablesPhoton",[_photon]) +\
               [ calculables.xcJet(_jet,
                                   gamma = _photon,
                                   gammaDR = 0.5,
                                   muon = _muon,
                                   muonDR = 0.5,
                                   correctForMuons = _correctForMuons,
                                   electron = _electron,
                                   electronDR = 0.5),
                 calculables.jetIndices( _jet, _jetPtMin,      etaMax = 3.0, flagName = params["jetId"]),
                 calculables.jetIndices( _jet, lowPtThreshold, etaMax = 3.0, flagName = params["jetId"], extraName = lowPtName),
                 calculables.muonIndices( _muon, ptMin = 20, combinedRelIsoMax = 0.15),
                 calculables.electronIndices( _electron, ptMin = 20, simpleEleID = "95", useCombinedIso = True),
                 calculables.photonIndicesPat(  ptMin = 20, flagName = "photonIDLooseFromTwikiPat"),
                 calculables.indicesUnmatched(collection = _photon, xcjets = _jet, DR = 0.5),
                 calculables.indicesUnmatched(collection = _electron, xcjets = _jet, DR = 0.5)
                 ] \
                 + [ calculables.jetSumP4(_jet, mcScaleFactor = 1.0),
                     calculables.deltaPhiStar(_jet, ptMin = lowPtThreshold, extraName = lowPtName),
                     calculables.deltaPseudoJet(_jet, _etRatherThanPt),
                     calculables.alphaT(_jet, _etRatherThanPt),
                     calculables.alphaTMet(_jet, _etRatherThanPt, _met),
                    #calculables.mhtMinusMetOverMeff(_jet, _met, _etRatherThanPt),
                     calculables.mhtMinusMetOverMeff(_jet, "metP4PF", _etRatherThanPt),
                     calculables.vertexID(),
                     calculables.vertexIndices(sumPtMin=20),
                     ]