Пример #1
0
    def listOfSteps(self,params) :
        _jet  = params["objects"]["jet"]
        _electron = params["objects"]["electron"]
        _muon = params["objects"]["muon"]
        _photon = params["objects"]["photon"]
        _met  = params["objects"]["met"]
        
        outList=[
            steps.progressPrinter(),
            steps.histogrammer("genpthat",200,0,1000,title=";#hat{p_{T}} (GeV);events / bin"),

            steps.preIdJetPtSelector(_jet,100.0,0),
            steps.preIdJetPtSelector(_jet, 80.0,1),
            steps.leadingUnCorrJetPtSelector( [_jet],100.0 ),
            steps.hltFilter("HLT_Jet50U"),
            steps.vertexRequirementFilter(),
            steps.techBitFilter([0],True),
            steps.physicsDeclared(),
            steps.monsterEventFilter(),
            steps.hbheNoiseFilter(),
            
            steps.multiplicityFilter("%sIndices%s"%_photon, nMax = 0),
            steps.multiplicityFilter("%sIndices%s"%_electron, nMax = 0),
            steps.multiplicityFilter("%sIndices%s"%_muon, nMax = 0),
            steps.multiplicityFilter("%sIndicesUnmatched%s"%_electron, nMax = 0),
            steps.multiplicityFilter("%sIndicesUnmatched%s"%_photon, nMax = 0),
            steps.multiplicityFilter("%sIndicesOther%s"%_muon, nMax = 0),
            steps.uniquelyMatchedNonisoMuons(_jet),
            steps.singleJetHistogrammer(_jet),
            steps.histogrammer(_met,100,0.0,500.0,title=";"+_met+" (GeV);events / bin", funcString = "lambda x: x.pt()"),
            steps.histogrammer("%sSumP4%s"%_jet,100,0.0,500.0,title=";MHT (GeV);events / bin", funcString = "lambda x: x.pt()"),

            steps.multiplicityFilter("%sIndices%s"%_jet, nMin = params["nJetsMinMax"][0], nMax = params["nJetsMinMax"][1]),
            steps.multiplicityFilter("%sIndicesOther%s"%_jet, nMax = 0),
            steps.jetPtSelector(_jet,100.0,0),
            steps.jetEtaSelector(_jet,2.5,0),
            steps.jetPtSelector(_jet,80.0,1),
            steps.singleJetHistogrammer(_jet),
            steps.histogrammer(_met,100,0.0,500.0,title=";"+_met+" (GeV);events / bin", funcString = "lambda x: x.pt()"),
            steps.histogrammer("%sSumP4%s"%_jet,100,0.0,500.0,title=";MHT (GeV);events / bin", funcString = "lambda x: x.pt()"),
            
            steps.variableGreaterFilter(350.0,"%sSumEt%s"%_jet, suffix = "GeV"),
            steps.singleJetHistogrammer(_jet),
            steps.histogrammer(_met,100,0.0,500.0,title=";"+_met+" (GeV);events / bin", funcString = "lambda x: x.pt()"),
            steps.histogrammer("%sSumP4%s"%_jet,100,0.0,500.0,title=";MHT (GeV);events / bin", funcString = "lambda x: x.pt()"),
            
            steps.variablePtGreaterFilter(140.0,"%sSumP4%s"%_jet,"GeV"),
            steps.variableGreaterFilter(0.55,"%sAlphaT%s"%_jet),
            steps.singleJetHistogrammer(_jet),
            steps.histogrammer(_met,100,0.0,500.0,title=";"+_met+" (GeV);events / bin", funcString = "lambda x: x.pt()"),
            steps.histogrammer("%sSumP4%s"%_jet,100,0.0,500.0,title=";MHT (GeV);events / bin", funcString = "lambda x: x.pt()"),
            
            ]
        return outList
Пример #2
0
 def listOfSteps(self,params) :
     return [
         steps.progressPrinter(),
         steps.jetPtSelector(params["jets"],100.0,0),
         steps.jetPtSelector(params["jets"],100.0,1),
         steps.jetEtaSelector(params["jets"],2.5,0),
         steps.lowestUnPrescaledTrigger(params["triggers"]),
         steps.vertexRequirementFilter(),
         steps.techBitFilter([0],True),
         steps.physicsDeclared(),
         steps.monsterEventFilter(),
         steps.hbheNoiseFilter(),
         
         steps.hltPrescaleHistogrammer(params["triggers"]),
         ]
Пример #3
0
    def listOfSteps(self,params) :
        _jet = params["objects"]["jet"]
        _muon = params["objects"]["muon"]
        _photon = params["objects"]["photon"]
        stepList = [
            steps.progressPrinter(),

            steps.jetPtSelector(_jet,100.0,0),
            steps.leadingUnCorrJetPtSelector( [_jet],100.0 ),
            steps.hltFilter("HLT_Jet50U"),
            steps.vertexRequirementFilter(5.0,24.0),
            steps.techBitFilter([0],True),
            steps.physicsDeclared(),
            steps.monsterEventFilter(10,0.25),
            steps.hbheNoiseFilter(),
        
            steps.multiplicityFilter("%sIndices%s"%_jet, nMin = params["nJetsMinMax"][0], nMax = params["nJetsMinMax"][1]),
            steps.multiplicityFilter("%sIndicesOther%s"%_jet, nMax = 0),
            steps.variableGreaterFilter(350.0,"%sSumPt%s"%_jet),
            #steps.multiplicityFilter("%sIndices%s"%_photon, nMax = 0),
            steps.multiplicityFilter("%sIndices%s"%_muon, nMax = 0),

            steps.alphaTetaDependence(_jet)

            # signal distributions for {g,Z}x{forward,mid,central}
            ]
        return stepList