HLTJetMETValidation Class Reference

#include <HLTJetMETValidation.h>

Inheritance diagram for HLTJetMETValidation:

Public Member Functions
	HLTJetMETValidation (const edm::ParameterSet &)
	~HLTJetMETValidation () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 54 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(), HLTBitAnalyser_cfi::hltresults, HLTriggerResults, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, hltTriggerMap, mps_fire::i, edm::HandleBase::isValid(), packedPFCandidateRefMixer_cfi::pf, PFJetAlgorithm, pfJetBenchmark_cfi::pfJets, trig_iter, triggerEventObject_, L1TEGammaOffline_cfi::triggerNames, 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  )

overrideprivatevirtual

Implements DQMEDAnalyzer.

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, dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), hltTrgJet, hltTrgMet, HLTConfigProvider::removeVersion(), dqm::dqmstoreimpl::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  )

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 36 of file HLTJetMETValidation.cc.

References hltConfig_, hltTrgJet, hltTrgJetLow, hltTrgMet, hltTrgMetLow, HLTConfigProvider::init(), dqmiolumiharvest::j, 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, edm::TriggerNames::triggerName(), and EgHLTOffTrigSelection_cfi::trigName.

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 94 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 96 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 97 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 95 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 98 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 100 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 101 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 99 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 90 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 92 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 93 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 91 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 107 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 109 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 110 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 108 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 81 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 83 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 84 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 82 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 85 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 87 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 88 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 86 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 77 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 79 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 80 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 78 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 103 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 105 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 106 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 104 of file HLTJetMETValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HLTJetMETValidation::_triggerResults

private

Definition at line 112 of file HLTJetMETValidation.h.

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

private

Definition at line 70 of file HLTJetMETValidation.h.

Referenced by analyze().

int HLTJetMETValidation::evtCnt

private

Definition at line 115 of file HLTJetMETValidation.h.

Referenced by analyze(), and HLTJetMETValidation().

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

private

Definition at line 69 of file HLTJetMETValidation.h.

Referenced by analyze().

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

private

Definition at line 71 of file HLTJetMETValidation.h.

Referenced by analyze().

HLTConfigProvider HLTJetMETValidation::hltConfig_

private

Definition at line 117 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

bool HLTJetMETValidation::HLTinit_

private

Definition at line 128 of file HLTJetMETValidation.h.

Referenced by getHLTResults().

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

private

Definition at line 72 of file HLTJetMETValidation.h.

Referenced by analyze().

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

private

Definition at line 118 of file HLTJetMETValidation.h.

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

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

private

Definition at line 119 of file HLTJetMETValidation.h.

Referenced by analyze(), and dqmBeginRun().

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

private

Definition at line 120 of file HLTJetMETValidation.h.

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

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

private

Definition at line 121 of file HLTJetMETValidation.h.

Referenced by analyze(), and dqmBeginRun().

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

private

Definition at line 125 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

std::vector<bool> HLTJetMETValidation::hlttrigs

private

Definition at line 124 of file HLTJetMETValidation.h.

std::string HLTJetMETValidation::patternJetTrg_

private

Definition at line 75 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

std::string HLTJetMETValidation::patternMetTrg_

private

Definition at line 75 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

std::string HLTJetMETValidation::patternMuTrg_

private

Definition at line 75 of file HLTJetMETValidation.h.

Referenced by dqmBeginRun().

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

private

Definition at line 68 of file HLTJetMETValidation.h.

Referenced by analyze().

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

private

Definition at line 126 of file HLTJetMETValidation.h.

Referenced by analyze(), and getHLTResults().

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

private

InputTag of TriggerEventWithRefs to analyze.

Definition at line 67 of file HLTJetMETValidation.h.

Referenced by analyze().

std::string HLTJetMETValidation::triggerTag_

private

Definition at line 75 of file HLTJetMETValidation.h.

Referenced by bookHistograms().

bool HLTJetMETValidation::writeFile_

private

Definition at line 130 of file HLTJetMETValidation.h.