Exemplo n.º 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) ]
Exemplo n.º 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"]
     _dR = params["dR"]
     _vPtMin = params["vetoPtMin"]
     
     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 = _dR, muon = _muon, muonDR = _dR, electron = _electron, electronDR = _dR,
                                 correctForMuons = not params["objects"]["muonsInJets"]),
              calculables.jetIndices( _jet, _jetPtMin, etaMax = 3.0, flagName = params["jetId"]),
              calculables.muonIndices( _muon, ptMin = _vPtMin, combinedRelIsoMax = 0.15),
              calculables.electronIndices( _electron, ptMin = _vPtMin, simpleEleID = "95", useCombinedIso = True),
              calculables.photonIndicesPat(  ptMin = _vPtMin, flagName = "photonIDLoosePat"),
              calculables.indicesUnmatched(collection = _photon, xcjets = _jet, DR = _dR),
              calculables.indicesUnmatched(collection = _electron, xcjets = _jet, DR = _dR)
              ] \
              + [calculables.deltaPseudoJet(_jet, True),
                 calculables.alphaT(_jet, True),
                 calculables.jetSumP4(_jet,1.0)]
Exemplo n.º 3
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),
                     ]
Exemplo n.º 4
0
 def listOfCalculables(self,params) :
     _jet =  params["objects"]["jet"]
     _muon = params["objects"]["muon"]
     _photon = params["objects"]["photon"]
     
     calcs =  calculables.zeroArgs() +\
             calculables.fromCollections("calculablesJet",[_jet]) +\
             calculables.fromCollections("calculablesLepton",[_muon]) +\
             [ calculables.xcJet( _jet,  gamma = _photon, gammaDR = 0.5, 
                                  muon = _muon, muonDR = 0.5),
               calculables.jetIndices( _jet, ptMin = 20.0, etaMax = 3.0, flagName = params["jetId"]),
               calculables.photonIndicesPat(  ptMin = 20, flagName="photonIDLoosePat"),
               calculables.muonIndices( _muon, ptMin = 20.0, combinedRelIsoMax = 0.15),
               calculables.genIndices( label = "Z", pdgs = [23]),
               ]
     return calcs
Exemplo n.º 5
0
 def listOfCalculables(self,params) :
     return calculables.zeroArgs() +\
            calculables.fromCollections("calculablesJet",[params["jets"]]) +\
            [calculables.jetIndices( params["jets"], ptMin = 50, etaMax = 3.0, flagName = "JetIDloose")]