HLTJetMETValidation Class Reference

#include <HLTJetMETValidation.h>

Inheritance diagram for HLTJetMETValidation:

Public Member Functions
	HLTJetMETValidation (const edm::ParameterSet &)
	~HLTJetMETValidation () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &iRun, edm::EventSetup const &iSetup) override
void	dqmBeginRun (edm::Run const &iRun, edm::EventSetup const &iSetup) override
void	getHLTResults (const edm::TriggerResults &, const edm::TriggerNames &triggerNames)

Private Attributes
std::vector< MonitorElement * >	_meGenJetEta
std::vector< MonitorElement * >	_meGenJetEtaTrg
std::vector< MonitorElement * >	_meGenJetEtaTrgLow
std::vector< MonitorElement * >	_meGenJetEtaTrgMC
std::vector< MonitorElement * >	_meGenJetPhi
std::vector< MonitorElement * >	_meGenJetPhiTrg
std::vector< MonitorElement * >	_meGenJetPhiTrgLow
std::vector< MonitorElement * >	_meGenJetPhiTrgMC
std::vector< MonitorElement * >	_meGenJetPt
std::vector< MonitorElement * >	_meGenJetPtTrg
std::vector< MonitorElement * >	_meGenJetPtTrgLow
std::vector< MonitorElement * >	_meGenJetPtTrgMC
std::vector< MonitorElement * >	_meGenMET
std::vector< MonitorElement * >	_meGenMETTrg
std::vector< MonitorElement * >	_meGenMETTrgLow
std::vector< MonitorElement * >	_meGenMETTrgMC
std::vector< MonitorElement * >	_meHLTJetEta
std::vector< MonitorElement * >	_meHLTJetEtaTrg
std::vector< MonitorElement * >	_meHLTJetEtaTrgLow
std::vector< MonitorElement * >	_meHLTJetEtaTrgMC
std::vector< MonitorElement * >	_meHLTJetPhi
std::vector< MonitorElement * >	_meHLTJetPhiTrg
std::vector< MonitorElement * >	_meHLTJetPhiTrgLow
std::vector< MonitorElement * >	_meHLTJetPhiTrgMC
std::vector< MonitorElement * >	_meHLTJetPt
std::vector< MonitorElement * >	_meHLTJetPtTrg
std::vector< MonitorElement * >	_meHLTJetPtTrgLow
std::vector< MonitorElement * >	_meHLTJetPtTrgMC
std::vector< MonitorElement * >	_meHLTMET
std::vector< MonitorElement * >	_meHLTMETTrg
std::vector< MonitorElement * >	_meHLTMETTrgLow
std::vector< MonitorElement * >	_meHLTMETTrgMC
MonitorElement *	_triggerResults
edm::EDGetTokenT< reco::CaloMETCollection >	CaloMETColl
int	evtCnt
edm::EDGetTokenT< reco::GenJetCollection >	GenJetAlgorithm
edm::EDGetTokenT< reco::GenMETCollection >	GenMETColl
HLTConfigProvider	hltConfig_
bool	HLTinit_
edm::EDGetTokenT< edm::TriggerResults >	HLTriggerResults
std::vector< std::string >	hltTrgJet
std::vector< std::string >	hltTrgJetLow
std::vector< std::string >	hltTrgMet
std::vector< std::string >	hltTrgMetLow
std::map< std::string, bool >	hltTriggerMap
std::vector< bool >	hlttrigs
std::string	patternJetTrg_
std::string	patternMetTrg_
std::string	patternMuTrg_
edm::EDGetTokenT< reco::PFJetCollection >	PFJetAlgorithm
std::map< std::string, bool >::iterator	trig_iter
edm::EDGetTokenT< trigger::TriggerEventWithRefs >	triggerEventObject_
	InputTag of TriggerEventWithRefs to analyze. More...
std::string	triggerTag_
bool	writeFile_

Detailed Description

Definition at line 55 of file HLTJetMETValidation.h.

Constructor & Destructor Documentation

HLTJetMETValidation::HLTJetMETValidation ( const edm::ParameterSet &  ps )

explicit

Definition at line 13 of file HLTJetMETValidation.cc.

References evtCnt.

    : triggerEventObject_(
          consumes<TriggerEventWithRefs>(ps.getUntrackedParameter<edm::InputTag>("triggerEventObject"))),
      PFJetAlgorithm(consumes<PFJetCollection>(ps.getUntrackedParameter<edm::InputTag>("PFJetAlgorithm"))),
      GenJetAlgorithm(consumes<GenJetCollection>(ps.getUntrackedParameter<edm::InputTag>("GenJetAlgorithm"))),
      CaloMETColl(consumes<CaloMETCollection>(ps.getUntrackedParameter<edm::InputTag>("CaloMETCollection"))),
      GenMETColl(consumes<GenMETCollection>(ps.getUntrackedParameter<edm::InputTag>("GenMETCollection"))),
      HLTriggerResults(consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("HLTriggerResults"))),
      triggerTag_(ps.getUntrackedParameter<std::string>("DQMFolder", "SingleJet")),
      patternJetTrg_(ps.getUntrackedParameter<std::string>("PatternJetTrg", "")),
      patternMetTrg_(ps.getUntrackedParameter<std::string>("PatternMetTrg", "")),
      patternMuTrg_(ps.getUntrackedParameter<std::string>("PatternMuTrg", "")),
      HLTinit_(false) {
  evtCnt = 0;
}

HLTJetMETValidation::~HLTJetMETValidation ( )

override

Definition at line 29 of file HLTJetMETValidation.cc.

{}

Member Function Documentation

void HLTJetMETValidation::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 206 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meHLTJetEta, _meHLTJetEtaTrg, _meHLTJetEtaTrgLow, _meHLTJetEtaTrgMC, _meHLTJetPhi, _meHLTJetPhiTrg, _meHLTJetPhiTrgLow, _meHLTJetPhiTrgMC, _meHLTJetPt, _meHLTJetPtTrg, _meHLTJetPtTrgLow, _meHLTJetPtTrgMC, _meHLTMET, _meHLTMETTrg, _meHLTMETTrgLow, _meHLTMETTrgMC, CaloMETColl, evtCnt, HcalObjRepresent::Fill(), relval_steps::gen(), GenJetAlgorithm, ttbarCategorization_cff::genJets, HepMCValidationHelper::genMet(), RazorAnalyzer::genmet, GenMETColl, edm::Event::getByToken(), getHLTResults(), HLTriggerResults, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, hltTriggerMap, mps_fire::i, edm::HandleBase::isValid(), packedPFCandidateRefMixer_cfi::pf, PFJetAlgorithm, pfJetBenchmark_cfi::pfJets, trig_iter, triggerEventObject_, and edm::Event::triggerNames().

                                                                                     {
  using namespace std;
  using namespace edm;
  using namespace reco;
  using namespace l1extra;
  using namespace trigger;

  evtCnt++;

  // get The triggerEvent
  Handle<TriggerEventWithRefs> trigEv;
  iEvent.getByToken(triggerEventObject_, trigEv);

  // get TriggerResults object
  bool gotHLT = true;
  std::vector<bool> myTrigJ;
  for (size_t it = 0; it < hltTrgJet.size(); it++)
    myTrigJ.push_back(false);
  std::vector<bool> myTrigJLow;
  for (size_t it = 0; it < hltTrgJetLow.size(); it++)
    myTrigJLow.push_back(false);
  std::vector<bool> myTrigM;
  for (size_t it = 0; it < hltTrgMet.size(); it++)
    myTrigM.push_back(false);
  std::vector<bool> myTrigMLow;
  for (size_t it = 0; it < hltTrgMetLow.size(); it++)
    myTrigMLow.push_back(false);

  Handle<TriggerResults> hltresults;
  iEvent.getByToken(HLTriggerResults, hltresults);
  if (!hltresults.isValid()) {
    // if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No
    // HLTRESULTS";
    gotHLT = false;
  }

  if (gotHLT) {
    const edm::TriggerNames &triggerNames = iEvent.triggerNames(*hltresults);
    getHLTResults(*hltresults, triggerNames);

    //---pick-up the jet trigger decisions
    for (size_t it = 0; it < hltTrgJet.size(); it++) {
      trig_iter = hltTriggerMap.find(hltTrgJet[it]);
      if (trig_iter == hltTriggerMap.end()) {
        // std::cout << "Could not find trigger path with name: " <<
        // _probefilter.label() << std::endl; if (evtCnt==1)
        // edm::LogWarning("HLTJetMETValidation") << "Could not find trigger
        // path with name: " << _probefilter.label();
      } else {
        myTrigJ[it] = trig_iter->second;
      }
      // std::cout<<hltTrgJet[it].c_str()<<" "<<myTrigJ[it]<<std::endl;
    }
    for (size_t it = 0; it < hltTrgJetLow.size(); it++) {
      trig_iter = hltTriggerMap.find(hltTrgJetLow[it]);
      if (trig_iter == hltTriggerMap.end()) {
        // std::cout << "Could not find trigger path with name: " <<
        // _probefilter.label() << std::endl; if (evtCnt==1)
        // edm::LogWarning("HLTJetMETValidation") << "Could not find trigger
        // path with name: " << _probefilter.label();
      } else {
        myTrigJLow[it] = trig_iter->second;
      }
      // std::cout<<hltTrgJetLow[it].c_str()<<" "<<myTrigJLow[it]<<std::endl;
    }
    //---pick-up the met trigger decisions
    for (size_t it = 0; it < hltTrgMet.size(); it++) {
      trig_iter = hltTriggerMap.find(hltTrgMet[it]);
      if (trig_iter == hltTriggerMap.end()) {
        // std::cout << "Could not find trigger path with name: " <<
        // _probefilter.label() << std::endl; if (evtCnt==1)
        // edm::LogWarning("HLTJetMETValidation") << "Could not find trigger
        // path with name: " << _probefilter.label();
      } else {
        myTrigM[it] = trig_iter->second;
      }
      // std::cout<<hltTrgMet[it].c_str()<<" "<<myTrigM[it]<<std::endl;
    }
    for (size_t it = 0; it < hltTrgMetLow.size(); it++) {
      trig_iter = hltTriggerMap.find(hltTrgMetLow[it]);
      if (trig_iter == hltTriggerMap.end()) {
        // std::cout << "Could not find trigger path with name: " <<
        // _probefilter.label() << std::endl; if (evtCnt==1)
        // edm::LogWarning("HLTJetMETValidation") << "Could not find trigger
        // path with name: " << _probefilter.label();
      } else {
        myTrigMLow[it] = trig_iter->second;
      }
      // std::cout<<hltTrgMetLow[it].c_str()<<" "<<myTrigMLow[it]<<std::endl;
    }
  }

  // --- Fill histos for PFJet paths ---
  // HLT jets namely hltAK4PFJets
  Handle<PFJetCollection> pfJets;
  iEvent.getByToken(PFJetAlgorithm, pfJets);
  double pfJetPt = -1.;
  double pfJetEta = -999.;
  double pfJetPhi = -999.;

  if (pfJets.isValid()) {
    // Loop over the PFJets and fill some histograms
    int jetInd = 0;
    for (PFJetCollection::const_iterator pf = pfJets->begin(); pf != pfJets->end(); ++pf) {
      // std::cout << "PF JET #" << jetInd << std::endl << pf->print() <<
      // std::endl;
      if (jetInd == 0) {
        pfJetPt = pf->pt();
        pfJetEta = pf->eta();
        pfJetPhi = pf->phi();
        for (size_t it = 0; it < hltTrgJet.size(); it++) {
          _meHLTJetPt[it]->Fill(pfJetPt);
          _meHLTJetEta[it]->Fill(pfJetEta);
          _meHLTJetPhi[it]->Fill(pfJetPhi);
          if (myTrigJ[it])
            _meHLTJetPtTrgMC[it]->Fill(pfJetPt);
          if (myTrigJ[it])
            _meHLTJetEtaTrgMC[it]->Fill(pfJetEta);
          if (myTrigJ[it])
            _meHLTJetPhiTrgMC[it]->Fill(pfJetPhi);
          if (myTrigJ[it] && myTrigJLow[it])
            _meHLTJetPtTrg[it]->Fill(pfJetPt);
          if (myTrigJ[it] && myTrigJLow[it])
            _meHLTJetEtaTrg[it]->Fill(pfJetEta);
          if (myTrigJ[it] && myTrigJLow[it])
            _meHLTJetPhiTrg[it]->Fill(pfJetPhi);
          if (myTrigJLow[it])
            _meHLTJetPtTrgLow[it]->Fill(pfJetPt);
          if (myTrigJLow[it])
            _meHLTJetEtaTrgLow[it]->Fill(pfJetEta);
          if (myTrigJLow[it])
            _meHLTJetPhiTrgLow[it]->Fill(pfJetPhi);
        }
        jetInd++;
      }
    }  // loop over pfjets
  } else {
    // std::cout << "  -- No PFJets" << std::endl;
    // if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No
    // PFJets";
  }

  // GenJets
  Handle<GenJetCollection> genJets;
  iEvent.getByToken(GenJetAlgorithm, genJets);
  double genJetPt = -1.;
  double genJetEta = -999.;
  double genJetPhi = -999.;

  if (genJets.isValid()) {
    // Loop over the GenJets and fill some histograms
    int jetInd = 0;
    for (GenJetCollection::const_iterator gen = genJets->begin(); gen != genJets->end(); ++gen) {
      if (jetInd == 0) {
        genJetPt = gen->pt();
        genJetEta = gen->eta();
        genJetPhi = gen->phi();
        for (size_t it = 0; it < hltTrgJet.size(); it++) {
          _meGenJetPt[it]->Fill(genJetPt);
          _meGenJetEta[it]->Fill(genJetEta);
          _meGenJetPhi[it]->Fill(genJetPhi);
          if (myTrigJ[it])
            _meGenJetPtTrgMC[it]->Fill(genJetPt);
          if (myTrigJ[it])
            _meGenJetEtaTrgMC[it]->Fill(genJetEta);
          if (myTrigJ[it])
            _meGenJetPhiTrgMC[it]->Fill(genJetPhi);
          if (myTrigJ[it] && myTrigJLow[it])
            _meGenJetPtTrg[it]->Fill(genJetPt);
          if (myTrigJ[it] && myTrigJLow[it])
            _meGenJetEtaTrg[it]->Fill(genJetEta);
          if (myTrigJ[it] && myTrigJLow[it])
            _meGenJetPhiTrg[it]->Fill(genJetPhi);
          if (myTrigJLow[it])
            _meGenJetPtTrgLow[it]->Fill(genJetPt);
          if (myTrigJLow[it])
            _meGenJetEtaTrgLow[it]->Fill(genJetEta);
          if (myTrigJLow[it])
            _meGenJetPhiTrgLow[it]->Fill(genJetPhi);
        }
        jetInd++;
      }
    }
  } else {
    // std::cout << "  -- No GenJets" << std::endl;
    // if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  -- No
    // GenJets";
  }

  // --- Fill histos for PFMET paths ---
  // HLT MET namely hltmet
  edm::Handle<CaloMETCollection> recmet;
  iEvent.getByToken(CaloMETColl, recmet);

  double calMet = -1;
  if (recmet.isValid()) {
    typedef CaloMETCollection::const_iterator cmiter;
    // std::cout << "Size of MET collection" <<  recmet.size() << std::endl;
    for (cmiter i = recmet->begin(); i != recmet->end(); i++) {
      calMet = i->pt();
      for (size_t it = 0; it < hltTrgMet.size(); it++) {
        _meHLTMET[it]->Fill(calMet);
        if (myTrigM.size() > it && myTrigM[it])
          _meHLTMETTrgMC[it]->Fill(calMet);
        if (myTrigM.size() > it && myTrigMLow.size() > it && myTrigM[it] && myTrigMLow[it])
          _meHLTMETTrg[it]->Fill(calMet);
        if (myTrigMLow.size() > it && myTrigMLow[it])
          _meHLTMETTrgLow[it]->Fill(calMet);
      }
    }
  } else {
    // std::cout << "  -- No MET Collection with name: " << CaloMETColl <<
    // std::endl; if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  --
    // No MET Collection with name: "<< CaloMETColl;
  }

  edm::Handle<GenMETCollection> genmet;
  iEvent.getByToken(GenMETColl, genmet);

  double genMet = -1;
  if (genmet.isValid()) {
    typedef GenMETCollection::const_iterator cmiter;
    for (cmiter i = genmet->begin(); i != genmet->end(); i++) {
      genMet = i->pt();
      for (size_t it = 0; it < hltTrgMet.size(); it++) {
        _meGenMET[it]->Fill(genMet);
        if (myTrigM.size() > it && myTrigM[it])
          _meGenMETTrgMC[it]->Fill(genMet);
        if (myTrigM.size() > it && myTrigMLow.size() > it && myTrigM[it] && myTrigMLow[it])
          _meGenMETTrg[it]->Fill(genMet);
        if (myTrigMLow.size() > it && myTrigMLow[it])
          _meGenMETTrgLow[it]->Fill(genMet);
      }
    }
  } else {
    // std::cout << "  -- No GenMET Collection with name: " << GenMETColl <<
    // std::endl; if (evtCnt==1) edm::LogWarning("HLTJetMETValidation") << "  --
    // No GenMET Collection with name: "<< GenMETColl;
  }
}

void HLTJetMETValidation::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 133 of file HLTJetMETValidation.cc.

References _meGenJetEta, _meGenJetEtaTrg, _meGenJetEtaTrgLow, _meGenJetEtaTrgMC, _meGenJetPhi, _meGenJetPhiTrg, _meGenJetPhiTrgLow, _meGenJetPhiTrgMC, _meGenJetPt, _meGenJetPtTrg, _meGenJetPtTrgLow, _meGenJetPtTrgMC, _meGenMET, _meGenMETTrg, _meGenMETTrgLow, _meGenMETTrgMC, _meHLTJetEta, _meHLTJetEtaTrg, _meHLTJetEtaTrgLow, _meHLTJetEtaTrgMC, _meHLTJetPhi, _meHLTJetPhiTrg, _meHLTJetPhiTrgLow, _meHLTJetPhiTrgMC, _meHLTJetPt, _meHLTJetPtTrg, _meHLTJetPtTrgLow, _meHLTJetPtTrgMC, _meHLTMET, _meHLTMETTrg, _meHLTMETTrgLow, _meHLTMETTrgMC, DQMStore::IBooker::book1D(), hltTrgJet, hltTrgMet, HLTConfigProvider::removeVersion(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and triggerTag_.

                                                                      {
  //----define DQM folders and elements
  for (size_t it = 0; it < hltTrgJet.size(); it++) {
    // std::cout<<hltTrgJet[it].c_str()<<"
    // "<<hltTrgJetLow[it].c_str()<<std::endl;
    std::string trgPathName = HLTConfigProvider::removeVersion(triggerTag_ + hltTrgJet[it]);
    // std::cout << "str = " << triggerTag_+hltTrgJet[it].c_str() << std::endl;
    // std::cout << "trgPathName = " << trgPathName << std::endl;
    iBooker.setCurrentFolder(trgPathName);
    _meHLTJetPt.push_back(iBooker.book1D("_meHLTJetPt", "Single HLT Jet Pt", 100, 0, 500));
    _meHLTJetPtTrgMC.push_back(iBooker.book1D("_meHLTJetPtTrgMC", "Single HLT Jet Pt - HLT Triggered", 100, 0, 500));
    _meHLTJetPtTrg.push_back(iBooker.book1D("_meHLTJetPtTrg", "Single HLT Jet Pt - HLT Triggered", 100, 0, 500));
    _meHLTJetPtTrgLow.push_back(
        iBooker.book1D("_meHLTJetPtTrgLow", "Single HLT Jet Pt - HLT Triggered Low", 100, 0, 500));

    _meHLTJetEta.push_back(iBooker.book1D("_meHLTJetEta", "Single HLT Jet Eta", 100, -10, 10));
    _meHLTJetEtaTrgMC.push_back(
        iBooker.book1D("_meHLTJetEtaTrgMC", "Single HLT Jet Eta - HLT Triggered", 100, -10, 10));
    _meHLTJetEtaTrg.push_back(iBooker.book1D("_meHLTJetEtaTrg", "Single HLT Jet Eta - HLT Triggered", 100, -10, 10));
    _meHLTJetEtaTrgLow.push_back(
        iBooker.book1D("_meHLTJetEtaTrgLow", "Single HLT Jet Eta - HLT Triggered Low", 100, -10, 10));

    _meHLTJetPhi.push_back(iBooker.book1D("_meHLTJetPhi", "Single HLT Jet Phi", 100, -4., 4.));
    _meHLTJetPhiTrgMC.push_back(
        iBooker.book1D("_meHLTJetPhiTrgMC", "Single HLT Jet Phi - HLT Triggered", 100, -4., 4.));
    _meHLTJetPhiTrg.push_back(iBooker.book1D("_meHLTJetPhiTrg", "Single HLT Jet Phi - HLT Triggered", 100, -4., 4.));
    _meHLTJetPhiTrgLow.push_back(
        iBooker.book1D("_meHLTJetPhiTrgLow", "Single HLT Jet Phi - HLT Triggered Low", 100, -4., 4.));

    _meGenJetPt.push_back(iBooker.book1D("_meGenJetPt", "Single Generated Jet Pt", 100, 0, 500));
    _meGenJetPtTrgMC.push_back(
        iBooker.book1D("_meGenJetPtTrgMC", "Single Generated Jet Pt - HLT Triggered", 100, 0, 500));
    _meGenJetPtTrg.push_back(iBooker.book1D("_meGenJetPtTrg", "Single Generated Jet Pt - HLT Triggered", 100, 0, 500));
    _meGenJetPtTrgLow.push_back(
        iBooker.book1D("_meGenJetPtTrgLow", "Single Generated Jet Pt - HLT Triggered Low", 100, 0, 500));

    _meGenJetEta.push_back(iBooker.book1D("_meGenJetEta", "Single Generated Jet Eta", 100, -10, 10));
    _meGenJetEtaTrgMC.push_back(
        iBooker.book1D("_meGenJetEtaTrgMC", "Single Generated Jet Eta - HLT Triggered", 100, -10, 10));
    _meGenJetEtaTrg.push_back(
        iBooker.book1D("_meGenJetEtaTrg", "Single Generated Jet Eta - HLT Triggered", 100, -10, 10));
    _meGenJetEtaTrgLow.push_back(
        iBooker.book1D("_meGenJetEtaTrgLow", "Single Generated Jet Eta - HLT Triggered Low", 100, -10, 10));

    _meGenJetPhi.push_back(iBooker.book1D("_meGenJetPhi", "Single Generated Jet Phi", 100, -4., 4.));
    _meGenJetPhiTrgMC.push_back(
        iBooker.book1D("_meGenJetPhiTrgMC", "Single Generated Jet Phi - HLT Triggered", 100, -4., 4.));
    _meGenJetPhiTrg.push_back(
        iBooker.book1D("_meGenJetPhiTrg", "Single Generated Jet Phi - HLT Triggered", 100, -4., 4.));
    _meGenJetPhiTrgLow.push_back(
        iBooker.book1D("_meGenJetPhiTrgLow", "Single Generated Jet Phi - HLT Triggered Low", 100, -4., 4.));
  }
  for (size_t it = 0; it < hltTrgMet.size(); it++) {
    // std::cout<<hltTrgMet[it].c_str()<<"
    // "<<hltTrgMetLow[it].c_str()<<std::endl;
    std::string trgPathName = HLTConfigProvider::removeVersion(triggerTag_ + hltTrgMet[it]);
    iBooker.setCurrentFolder(trgPathName);
    _meHLTMET.push_back(iBooker.book1D("_meHLTMET", "HLT Missing ET", 100, 0, 500));
    _meHLTMETTrgMC.push_back(iBooker.book1D("_meHLTMETTrgMC", "HLT Missing ET - HLT Triggered", 100, 0, 500));
    _meHLTMETTrg.push_back(iBooker.book1D("_meHLTMETTrg", "HLT Missing ET - HLT Triggered", 100, 0, 500));
    _meHLTMETTrgLow.push_back(iBooker.book1D("_meHLTMETTrgLow", "HLT Missing ET - HLT Triggered Low", 100, 0, 500));

    _meGenMET.push_back(iBooker.book1D("_meGenMET", "Generated Missing ET", 100, 0, 500));
    _meGenMETTrgMC.push_back(iBooker.book1D("_meGenMETTrgMC", "Generated Missing ET - HLT Triggered", 100, 0, 500));
    _meGenMETTrg.push_back(iBooker.book1D("_meGenMETTrg", "Generated Missing ET - HLT Triggered", 100, 0, 500));
    _meGenMETTrgLow.push_back(
        iBooker.book1D("_meGenMETTrgLow", "Generated Missing ET - HLT Triggered Low", 100, 0, 500));
  }
}

void HLTJetMETValidation::dqmBeginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 36 of file HLTJetMETValidation.cc.

References hltConfig_, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, HLTConfigProvider::init(), patternJetTrg_, patternMetTrg_, patternMuTrg_, AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerNames().

                                                                                     {
  bool foundMuTrg = false;
  std::string trgMuNm;
  bool changedConfig = true;

  //--define search patterns
  TPRegexp patternJet(patternJetTrg_);
  TPRegexp patternMet(patternMetTrg_);
  TPRegexp patternMu(patternMuTrg_);

  if (!hltConfig_.init(iRun, iSetup, "HLT", changedConfig)) {
    edm::LogError("HLTJetMETValidation") << "Initialization of HLTConfigProvider failed!!";
    return;
  }

  std::vector<std::string> validTriggerNames = hltConfig_.triggerNames();
  for (size_t j = 0; j < validTriggerNames.size(); j++) {
    //---find the muon path
    if (TString(validTriggerNames[j]).Contains(patternMu)) {
      // std::cout <<validTriggerNames[j].c_str()<<std::endl;
      if (!foundMuTrg)
        trgMuNm = validTriggerNames[j];
      foundMuTrg = true;
    }
    //---find the jet paths
    if (TString(validTriggerNames[j]).Contains(patternJet)) {
      hltTrgJet.push_back(validTriggerNames[j]);
    }
    //---find the met paths
    if (TString(validTriggerNames[j]).Contains(patternMet)) {
      hltTrgMet.push_back(validTriggerNames[j]);
    }
  }

  //----set the denominator paths
  for (size_t it = 0; it < hltTrgJet.size(); it++) {
    if (it == 0 && foundMuTrg)
      hltTrgJetLow.push_back(trgMuNm);  //--lowest threshold uses muon path
    if (it == 0 && !foundMuTrg)
      hltTrgJetLow.push_back(hltTrgJet[it]);  //---if no muon then itself
    if (it != 0)
      hltTrgJetLow.push_back(hltTrgJet[it - 1]);
    // std::cout<<hltTrgJet[it].c_str()<<"
    // "<<hltTrgJetLow[it].c_str()<<std::endl;
  }
  int itm(0), itpm(0), itmh(0), itpmh(0);
  for (size_t it = 0; it < hltTrgMet.size(); it++) {
    if (TString(hltTrgMet[it]).Contains("PF")) {
      if (TString(hltTrgMet[it]).Contains("MHT")) {
        if (0 == itpmh) {
          if (foundMuTrg)
            hltTrgMetLow.push_back(trgMuNm);
          else
            hltTrgMetLow.push_back(hltTrgMet[it]);
        } else
          hltTrgMetLow.push_back(hltTrgMet[it - 1]);
        itpmh++;
      }
      if (TString(hltTrgMet[it]).Contains("MET")) {
        if (0 == itpm) {
          if (foundMuTrg)
            hltTrgMetLow.push_back(trgMuNm);
          else
            hltTrgMetLow.push_back(hltTrgMet[it]);
        } else
          hltTrgMetLow.push_back(hltTrgMet[it - 1]);
        itpm++;
      }
    } else {
      if (TString(hltTrgMet[it]).Contains("MHT")) {
        if (0 == itmh) {
          if (foundMuTrg)
            hltTrgMetLow.push_back(trgMuNm);
          else
            hltTrgMetLow.push_back(hltTrgMet[it]);
        } else
          hltTrgMetLow.push_back(hltTrgMet[it - 1]);
        itmh++;
      }
      if (TString(hltTrgMet[it]).Contains("MET")) {
        if (0 == itm) {
          if (foundMuTrg)
            hltTrgMetLow.push_back(trgMuNm);
          else
            hltTrgMetLow.push_back(hltTrgMet[it]);
        } else
          hltTrgMetLow.push_back(hltTrgMet[it - 1]);
        itm++;
      }
    }
    // std::cout<<hltTrgMet[it].c_str()<<"
    // "<<hltTrgMetLow[it].c_str()<<std::endl;
  }
}

void HLTJetMETValidation::getHLTResults ( const edm::TriggerResults &  hltresults, const edm::TriggerNames &  triggerNames  )

private

Definition at line 447 of file HLTJetMETValidation.cc.

References accept(), edm::HLTGlobalStatus::accept(), HLTinit_, hltTriggerMap, edm::HLTGlobalStatus::size(), AlCaHLTBitMon_QueryRunRegistry::string, trig_iter, and edm::TriggerNames::triggerName().

Referenced by analyze().

                                                                                                                {
  int ntrigs = hltresults.size();
  if (!HLTinit_) {
    HLTinit_ = true;

    for (int itrig = 0; itrig != ntrigs; ++itrig) {
      // std::cout << "trigger " << itrig << ": " << trigName << std::endl;
    }
  }

  for (int itrig = 0; itrig != ntrigs; ++itrig) {
    std::string trigName = triggerNames.triggerName(itrig);
    bool accept = hltresults.accept(itrig);

    // if (accept) _triggerResults->Fill(float(itrig));

    // fill the trigger map
    typedef std::map<std::string, bool>::value_type valType;
    trig_iter = hltTriggerMap.find(trigName);
    if (trig_iter == hltTriggerMap.end())
      hltTriggerMap.insert(valType(trigName, accept));
    else
      trig_iter->second = accept;
  }
}

Member Data Documentation

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetEta

private

Definition at line 95 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetEtaTrg

private

Definition at line 97 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetEtaTrgLow

private

Definition at line 98 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetEtaTrgMC

private

Definition at line 96 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPhi

private

Definition at line 99 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPhiTrg

private

Definition at line 101 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPhiTrgLow

private

Definition at line 102 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPhiTrgMC

private

Definition at line 100 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPt

private

Definition at line 91 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPtTrg

private

Definition at line 93 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPtTrgLow

private

Definition at line 94 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenJetPtTrgMC

private

Definition at line 92 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenMET

private

Definition at line 108 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenMETTrg

private

Definition at line 110 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenMETTrgLow

private

Definition at line 111 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meGenMETTrgMC

private

Definition at line 109 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetEta

private

Definition at line 82 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetEtaTrg

private

Definition at line 84 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetEtaTrgLow

private

Definition at line 85 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetEtaTrgMC

private

Definition at line 83 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPhi

private

Definition at line 86 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPhiTrg

private

Definition at line 88 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPhiTrgLow

private

Definition at line 89 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPhiTrgMC

private

Definition at line 87 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPt

private

Definition at line 78 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPtTrg

private

Definition at line 80 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPtTrgLow

private

Definition at line 81 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTJetPtTrgMC

private

Definition at line 79 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTMET

private

Definition at line 104 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTMETTrg

private

Definition at line 106 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTMETTrgLow

private

Definition at line 107 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> HLTJetMETValidation::_meHLTMETTrgMC

private

Definition at line 105 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HLTJetMETValidation::_triggerResults

private

Definition at line 113 of file HLTJetMETValidation.h.

edm::EDGetTokenT<reco::CaloMETCollection> HLTJetMETValidation::CaloMETColl

private

Definition at line 71 of file HLTJetMETValidation.h.

Referenced by analyze().

int HLTJetMETValidation::evtCnt

private

Definition at line 116 of file HLTJetMETValidation.h.

Referenced by analyze(), and HLTJetMETValidation().

edm::EDGetTokenT<reco::GenJetCollection> HLTJetMETValidation::GenJetAlgorithm

private

Definition at line 70 of file HLTJetMETValidation.h.

Referenced by analyze().

edm::EDGetTokenT<reco::GenMETCollection> HLTJetMETValidation::GenMETColl

private

Definition at line 72 of file HLTJetMETValidation.h.

Referenced by analyze().

HLTConfigProvider HLTJetMETValidation::hltConfig_

private

Definition at line 118 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

bool HLTJetMETValidation::HLTinit_

private

Definition at line 129 of file HLTJetMETValidation.h.

Referenced by getHLTResults().

edm::EDGetTokenT<edm::TriggerResults> HLTJetMETValidation::HLTriggerResults

private

Definition at line 73 of file HLTJetMETValidation.h.

Referenced by analyze().

std::vector<std::string> HLTJetMETValidation::hltTrgJet

private

Definition at line 119 of file HLTJetMETValidation.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgJetLow

private

Definition at line 120 of file HLTJetMETValidation.h.

Referenced by analyze(), and dqmBeginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgMet

private

Definition at line 121 of file HLTJetMETValidation.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

std::vector<std::string> HLTJetMETValidation::hltTrgMetLow

private

Definition at line 122 of file HLTJetMETValidation.h.

Referenced by analyze(), and dqmBeginRun().

std::map<std::string, bool> HLTJetMETValidation::hltTriggerMap

private

Definition at line 126 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

std::vector<bool> HLTJetMETValidation::hlttrigs

private

Definition at line 125 of file HLTJetMETValidation.h.

std::string HLTJetMETValidation::patternJetTrg_

private

Definition at line 76 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

std::string HLTJetMETValidation::patternMetTrg_

private

Definition at line 76 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

std::string HLTJetMETValidation::patternMuTrg_

private

Definition at line 76 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

edm::EDGetTokenT<reco::PFJetCollection> HLTJetMETValidation::PFJetAlgorithm

private

Definition at line 69 of file HLTJetMETValidation.h.

Referenced by analyze().

std::map<std::string, bool>::iterator HLTJetMETValidation::trig_iter

private

Definition at line 127 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

edm::EDGetTokenT<trigger::TriggerEventWithRefs> HLTJetMETValidation::triggerEventObject_

private

InputTag of TriggerEventWithRefs to analyze.

Definition at line 68 of file HLTJetMETValidation.h.

Referenced by analyze().

std::string HLTJetMETValidation::triggerTag_

private

Definition at line 76 of file HLTJetMETValidation.h.

Referenced by bookHistograms().

bool HLTJetMETValidation::writeFile_

private

Definition at line 131 of file HLTJetMETValidation.h.