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)]
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) ]
def listOfCalculables(self, params) : obj = params["objects"] outList = calculables.zeroArgs() outList += calculables.fromCollections(calculables.Muon, [obj["muon"]]) outList += calculables.fromCollections(calculables.Electron, [obj["electron"]]) outList += calculables.fromCollections(calculables.Photon, [obj["photon"]]) outList += self.calcListOther(obj, params["triggerList"]) outList += self.calcListJet(obj, params["etRatherThanPt"], params["lowPtThreshold"], params["lowPtName"]) return outList
def listOfCalculables(self, params) : if params["vertexMode"] : assert params["photonId"] in ["photonIDTightFromTwikiPat"],"In vertexMode but requested %s"%params["photonId"] obj = params["objects"] _etRatherThanPt = params["etRatherThanPt"] return calculables.zeroArgs() +\ calculables.fromCollections(calculables.Jet,[obj["jet"]]) +\ calculables.fromCollections(calculables.Muon,[obj["muon"]]) +\ calculables.fromCollections(calculables.Electron,[obj["electron"]]) +\ calculables.fromCollections(calculables.Photon,[obj["photon"]]) +\ [ calculables.XClean.xcJet( obj["jet"], gamma = obj["photon"], gammaDR = 0.5, muon = obj["muon"], muonDR = 0.5, correctForMuons = not obj["muonsInJets"], electron = obj["electron"], electronDR = 0.5 ), calculables.Jet.Indices( obj["jet"], ptMin = params["thresholds"][3], etaMax = 3.0, flagName = params["jetId"]), calculables.Jet.Indices( obj["jet"], ptMin = params["lowPtThreshold"], etaMax = 3.0, flagName = params["jetId"], extraName = params["lowPtName"]), calculables.Jet.Indices( obj["jet"], ptMin = params["highPtThreshold"], etaMax = 3.0, flagName = params["jetId"], extraName = params["highPtName"]), calculables.Muon.Indices( obj["muon"], ptMin = 10, combinedRelIsoMax = 0.15), calculables.Electron.Indices( obj["electron"], ptMin = 10, simpleEleID = "95", useCombinedIso = True), calculables.Photon.Indices(obj["photon"], ptMin = 25, flagName = params["photonId"]), calculables.Gen.genIndices( pdgs = [22], label = "Status3Photon", status = [3]), calculables.Gen.genMinDeltaRPhotonOther( label = "Status3Photon"), calculables.Gen.genIndices( pdgs = [22], label = "Status1Photon", status = [1]), calculables.Gen.genIsolations(label = "Status1Photon", coneSize = 0.4), calculables.Gen.genPhotonCategory(label = "Status1Photon"), calculables.XClean.IndicesUnmatched(collection = obj["photon"], xcjets = obj["jet"], DR = 0.5), calculables.XClean.IndicesUnmatched(collection = obj["electron"], xcjets = obj["jet"], DR = 0.5) ] \ + [ calculables.Jet.SumP4(obj["jet"]), calculables.Jet.SumP4(obj["jet"], extraName = params["lowPtName"]), calculables.Jet.SumP4(obj["jet"], extraName = params["highPtName"]), calculables.Jet.DeltaPhiStar(obj["jet"], extraName = ""), calculables.Jet.DeltaPhiStar(obj["jet"], extraName = params["lowPtName"]), calculables.Jet.DeltaPseudoJet(obj["jet"], _etRatherThanPt), calculables.Jet.AlphaTWithPhoton1PtRatherThanMht(obj["jet"], photons = obj["photon"], etRatherThanPt = _etRatherThanPt), calculables.Jet.AlphaT(obj["jet"], _etRatherThanPt), calculables.Jet.AlphaTMet(obj["jet"], _etRatherThanPt, obj["met"]), calculables.Jet.MhtOverMet((obj["jet"][0], obj["jet"][1]+params["highPtName"]), met = obj["met"]), calculables.Jet.MhtOverMet((obj["jet"][0], obj["jet"][1]+params["highPtName"]), met = "%sPlus%s%s"%(obj["met"], obj["photon"][0], obj["photon"][1])), calculables.Other.metPlusParticles(met = obj["met"], particles = obj["photon"]), calculables.Other.minDeltaRToJet(obj["photon"], obj["jet"]), calculables.Other.SumP4(obj["photon"]), calculables.Vertex.ID(), calculables.Vertex.Indices(), calculables.Jet.deadEcalDR(obj["jet"], extraName = params["lowPtName"], minNXtals = 10), calculables.Other.lowestUnPrescaledTrigger(params["triggerList"]), ]
def listOfCalculables(self, pars) : outList = calculables.zeroArgs() outList += calculables.fromCollections(calculables.Muon, [pars["muon"]]) outList += calculables.fromCollections(calculables.Electron, [pars["electron"]]) outList +=[calculables.Muon.Indices( pars["muon"], ptMin = 10, combinedRelIsoMax = 0.15), calculables.Muon.IndicesAnyIsoIsoOrder(pars['muon'], "CombinedRelativeIso"), calculables.Muon.LeadingIsoAny(pars['muon'], ptMin = 18, iso = "CombinedRelativeIso"), calculables.Electron.Indices( pars["electron"], ptMin = 10, simpleEleID = "95", useCombinedIso = True), calculables.Other.PtSorted(pars['muon']) ] return outList
def listOfCalculables(self, pars) : obj = pars["objects"] lepton = obj[pars["lepton"]["name"]] outList = calculables.zeroArgs() outList += calculables.fromCollections(calculables.Muon, [obj["muon"]]) outList += calculables.fromCollections(calculables.Electron, [obj["electron"]]) outList += calculables.fromCollections(calculables.Photon, [obj["photon"]]) outList += calculables.fromCollections(calculables.Jet, [obj["jet"]]) outList += [ calculables.Jet.IndicesBtagged(obj["jet"],pars["bVar"]), calculables.Jet.Indices( obj["jet"], ptMin = 20, etaMax = 3.5, flagName = "JetIDloose"), calculables.Muon.Indices( obj["muon"], ptMin = 10, combinedRelIsoMax = 0.15), calculables.Muon.IndicesTriggering(obj["muon"]), calculables.Electron.Indices( obj["electron"], ptMin = 10, simpleEleID = "80", useCombinedIso = True), calculables.Photon.Indices( obj["photon"], ptMin = 25, flagName = "photonIDLooseFromTwikiPat"), calculables.XClean.IndicesUnmatched(collection = obj["photon"], xcjets = obj["jet"], DR = 0.5), calculables.XClean.IndicesUnmatched(collection = obj["electron"], xcjets = obj["jet"], DR = 0.5), calculables.XClean.xcJet(obj["jet"], applyResidualCorrectionsToData = False, gamma = obj["photon"], gammaDR = 0.5, electron = obj["electron"], electronDR = 0.5, muon = obj["muon"], muonDR = 0.5, correctForMuons = not obj["muonsInJets"]), calculables.XClean.SumP4(obj["jet"], obj["photon"], obj["electron"], obj["muon"]), calculables.XClean.SumPt(obj["jet"], obj["photon"], obj["electron"], obj["muon"]), calculables.Vertex.ID(), calculables.Vertex.Indices(), calculables.Other.lowestUnPrescaledTrigger(zip(*pars["lepton"]["triggers"])[0]), calculables.Top.mixedSumP4(transverse = obj["met"], longitudinal = obj["sumP4"]), calculables.Other.pt("mixedSumP4"), calculables.Top.SemileptonicTopIndex(lepton), calculables.Top.fitTopLeptonCharge(lepton), calculables.Top.TopReconstruction(lepton,obj["jet"],"mixedSumP4"), calculables.Other.Mt(lepton,"mixedSumP4", allowNonIso=True, isSumP4=True), calculables.Muon.IndicesAnyIsoIsoOrder(obj[pars["lepton"]["name"]], pars["lepton"]["isoVar"]), calculables.Other.PtSorted(obj['muon']), calculables.Other.Covariance(('met','PF')), calculables.Other.abbreviation( "TrkCountingHighEffBJetTags", "NTrkHiEff", fixes = calculables.Jet.xcStrip(obj['jet']) ), calculables.Other.abbreviation( "nVertexRatio", "nvr" ), calculables.Other.abbreviation('muonTriggerWeightPF','tw'), calculables.Jet.pt( obj['jet'], index = 0, Btagged = True ), calculables.Jet.absEta( obj['jet'], index = 3, Btagged = False) ] outList += calculables.fromCollections(calculables.Top,[('genTop',""),('fitTop',"")]) outList.append( calculables.Top.TopComboQQBBLikelihood(pars['objects']['jet'], pars['bVar'])) outList.append( calculables.Top.OtherJetsLikelihood(pars['objects']['jet'], pars['bVar'])) outList.append( calculables.Top.TopRatherThanWProbability(priorTop=0.5) ) return outList
def calcList(jet, met, photon, muon, electron, muonsInJets, jetIdFlag) : outList = [ calculables.XClean.xcJet(jet, applyResidualCorrectionsToData = False, gamma = photon, gammaDR = 0.5, muon = muon, muonDR = 0.5, correctForMuons = not muonsInJets, electron = electron, electronDR = 0.5), calculables.Jet.Indices( jet, ptMin = ptMin, etaMax = 3.0, flagName = jetIdFlag), calculables.Jet.Indices( jet, ptMin = lowPtThreshold, etaMax = 3.0, flagName = jetIdFlag, extraName = lowPtName), calculables.Jet.Indices( jet, ptMin = highPtThreshold, etaMax = 3.0, flagName = jetIdFlag, extraName = highPtName), calculables.Jet.SumP4(jet), calculables.Jet.SumP4(jet, extraName = lowPtName), calculables.Jet.SumP4(jet, extraName = highPtName), calculables.Jet.DeltaPhiStar(jet, extraName = lowPtName), calculables.Jet.DeltaPseudoJet(jet, etRatherThanPt), calculables.Jet.AlphaT(jet, etRatherThanPt), calculables.Jet.AlphaTMet(jet, etRatherThanPt, met), calculables.Jet.MhtOverMet((jet[0], jet[1]+highPtName), met = "%sPlus%s%s"%(obj["met"], obj["muon"][0], obj["muon"][1])), calculables.Jet.deadEcalDR(jet, extraName = lowPtName, minNXtals = 10), ] return outList+calculables.fromCollections(calculables.Jet, [jet])
def listOfCalculables(self, pars) : return ( calculables.zeroArgs() + calculables.fromCollections(calculables.Top,[('genTop',""),('fitTop',"")]) + [ calculables.Vertex.ID(), calculables.Vertex.Indices(), ] )
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
def listOfCalculables(self, pars) : outList = calculables.zeroArgs() outList += [ calculables.Vertex.ID(), calculables.Vertex.Indices(), ] outList += calculables.fromCollections(calculables.Top,[('genTop',""),('fitTop',"")]) return outList
def listOfCalculables(self,pars) : return (calculables.zeroArgs() + calculables.fromCollections(calculables.Muon,[pars["muon"]]) + [calculables.Muon.Indices( pars["muon"], ptMin = 10, combinedRelIsoMax = 0.15), calculables.Muon.TriggeringIndex(pars['muon'], ptMin = 10), calculables.Other.lowestUnPrescaledTrigger(zip(*pars["triggers"])[0]), calculables.Vertex.ID(), calculables.Vertex.Indices(), ])
def listOfCalculables(self,pars) : from calculables import Muon,Vertex muon = ("muon","PF") return ( calculables.zeroArgs() + calculables.fromCollections(calculables.Muon, [muon]) + [ Muon.Indices( muon, ptMin = 10, combinedRelIsoMax = 0.15), Vertex.ID(), Vertex.Indices(), ] )
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), ]
def listOfCalculables(self,pars) : return (calculables.zeroArgs() + calculables.fromCollections(calculables.Muon,[pars["muon"]]) + [calculables.Muon.Indices( pars["muon"], ptMin = 10, combinedRelIsoMax = 0.15), calculables.Muon.IndicesTriggering( pars['muon'] ), calculables.Vertex.ID(), calculables.Vertex.Indices(), calculables.Other.abbreviation( "nVertexRatio", "nvr" ), calculables.Other.abbreviation('muonTriggerWeightPF','tw') ])
def listOfCalculables(self, params) : outList = calculables.zeroArgs() for muon in nameList(params["recoAlgos"],"muon") : outList += calculables.fromCollections(calculables.Muon, [muon]) outList += [calculables.Muon.Indices(muon, ptMin = 10, combinedRelIsoMax = 0.15)] for obj in dict(params["recoAlgos"]).values() : outList += self.calcListJet(obj) return outList
def calcListJet(self, obj) : outList = [ calculables.XClean.xcJet(obj["jet"], gamma = None, gammaDR = 0.5, muon = obj["muon"], muonDR = 0.5, correctForMuons = not obj["muonsInJets"], electron = None, electronDR = 0.5), calculables.Jet.Indices( obj["jet"], obj["jetPtMin"], etaMax = 3.0, flagName = "JetIDloose"), ] return outList+calculables.fromCollections(calculables.Jet, [obj["jet"]])
def listOfCalculables(self,config) : jets = config["jets"] minJetPt = config["minJetPt"] listOfCalculables = calculables.zeroArgs() listOfCalculables += calculables.fromCollections(calculables.Jet,[jets]) listOfCalculables += [ calculables.Jet.Indices( jets, ptMin = minJetPt, etaMax = 3.0, flagName = "JetIDloose"), calculables.Jet.SumP4( jets), calculables.Jet.DeltaPhiStar( jets ), calculables.Jet.AlphaT ( jets, config["etRatherThanPt"]), calculables.Jet.DeltaPseudoJet( jets, config["etRatherThanPt"]), ] return listOfCalculables
def listOfCalculables(self,params) : return calculables.zeroArgs() +\ calculables.fromCollections("calculablesJet",[params["jets"]]) +\ [calculables.jetIndices( params["jets"], ptMin = 50, etaMax = 3.0, flagName = "JetIDloose")]
def listOfCalculables(self,pars) : return calculables.zeroArgs() +\ calculables.fromCollections(calculables.Electron,[electron]) +\ [calculables.Electron.Indices( electron, ptMin = 10, simpleEleID = "95", useCombinedIso = True)]
def listOfCalculables(self,params) : return calculables.zeroArgs() +\ calculables.fromCollections(calculables.jet,jetAlgoList) +\ [calculables.jet.Indices( jet, 30.0, etaMax = 3.0, flagName = "JetIDloose") for jet in jetAlgoList]
def listOfCalculables(self,params) : return calculables.zeroArgs() +\ calculables.fromCollections(calculables.Muon,[muon]) +\ [calculables.Muon.Indices( muon, ptMin = 10, combinedRelIsoMax = 0.50)]
def listOfCalculables(self,params) : return calculables.zeroArgs() +\ calculables.fromCollections("calculablesElectron",[params["electrons"]]) +\ [calculables.electronIndices( params["electrons"], ptMin = 20, simpleEleID = "80", useCombinedIso = params["useCombinedIso"])]
def listOfCalculables(self,params) : return calculables.zeroArgs() +\ calculables.fromCollections(calculables.Photon,[params["photon"]]) +\ [calculables.Photon.Indices(collection = params["photon"], ptMin = 100, flagName = "photonIDLooseFromTwikiPat")]