HTMHTAnalyzer Class Reference

#include <HTMHTAnalyzer.h>

Inheritance diagram for HTMHTAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &)
	Get the analysis. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	HTMHTAnalyzer (const edm::ParameterSet &)
	Constructor. More...
	~HTMHTAnalyzer () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Public Attributes
int	evtCounter

Private Attributes
double	_ptThreshold
std::string	_source
int	_trig_JetMB
int	_verbose
MonitorElement *	hHT
std::vector< std::string >	HLTPathsJetMBByName_
MonitorElement *	hMHT
MonitorElement *	hMHTPhi
MonitorElement *	hMHx
MonitorElement *	hMHy
MonitorElement *	hNevents
MonitorElement *	hNJets
MonitorElement *	jetME
std::string	metname
edm::ParameterSet	parameters
edm::InputTag	theJetCollectionForHTMHTLabel

Detailed Description

DQM monitoring source for HTMHT

Author

K. Hatakeyama, Rockefeller University

Definition at line 31 of file HTMHTAnalyzer.h.

Constructor & Destructor Documentation

HTMHTAnalyzer::HTMHTAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 29 of file HTMHTAnalyzer.cc.

                                                        {

  parameters = pSet;
  _ptThreshold = 30.;

}

HTMHTAnalyzer::~HTMHTAnalyzer ( )

override

Destructor.

Definition at line 37 of file HTMHTAnalyzer.cc.

{ }

Member Function Documentation

void HTMHTAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::TriggerResults &	triggerResults
	)

Get the analysis.

Definition at line 75 of file HTMHTAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), gather_cfg::cout, edm::Event::getByLabel(), mps_fire::i, edm::HandleBase::isValid(), LogTrace, metname, nanoDQM_cfi::MHT, njet, edm::HLTGlobalStatus::size(), edm::TriggerNames::triggerIndex(), and edm::Event::triggerNames().

                                                           {

  LogTrace(metname)<<"[HTMHTAnalyzer] Analyze HT & MHT";

  jetME->Fill(4);

  // ==========================================================  
  // Trigger information 
  //


  //
  //
  // Check how many HLT triggers are in triggerResults 
  int ntrigs = triggerResults.size();
  if (_verbose) std::cout << "ntrigs=" << ntrigs << std::endl;
  
  //
  //
  // If index=ntrigs, this HLT trigger doesn't exist in the HLT table for this data.
  const edm::TriggerNames & triggerNames = iEvent.triggerNames(triggerResults);
  
  //
  //
  // count number of requested Jet or MB HLT paths which have fired
  for (unsigned int i=0; i!=HLTPathsJetMBByName_.size(); i++) {
    unsigned int triggerIndex = triggerNames.triggerIndex(HLTPathsJetMBByName_[i]);
    if (triggerIndex<triggerResults.size()) {
      if (triggerResults.accept(triggerIndex)) {
        _trig_JetMB++;
      }
    }
  }
  // for empty input vectors (n==0), take all HLT triggers!
  if (HLTPathsJetMBByName_.empty()) _trig_JetMB=triggerResults.size()-1;

  if (_trig_JetMB==0) return;
   
  // ==========================================================

  // **** Get the Calo Jet container
  edm::Handle<reco::CaloJetCollection> jetcoll;

  // **** Get the SISCone Jet container
  iEvent.getByLabel(theJetCollectionForHTMHTLabel, jetcoll);

  if(!jetcoll.isValid()) return;

  // ==========================================================
  // Reconstructed HT & MHT Information

  int    njet=0;
  double MHx=0.;
  double MHy=0.;
  double MHT=0.;
  double MHTPhi=0.;
  double HT=0.;

  for (reco::CaloJetCollection::const_iterator calojet = jetcoll->begin(); calojet!=jetcoll->end(); ++calojet){
    if (calojet->pt()>_ptThreshold){
      njet++;
      MHx += -1.*calojet->px();
      MHy += -1.*calojet->py();
      HT  += calojet->pt();
    }
  }

  TVector2 *tv2 = new TVector2(MHx,MHy);
  MHT   =tv2->Mod();
  MHTPhi=tv2->Phi();

  //std::cout << "HTMHT " << MHT << std::endl;

  hNevents->Fill(1.);
  hNJets->Fill(njet);
  if (njet>0){
    hMHx->Fill(MHx);
    hMHy->Fill(MHy);
    hMHT->Fill(MHT);
    hMHTPhi->Fill(MHTPhi);
    hHT->Fill(HT);
  }

}

void HTMHTAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

override

Definition at line 41 of file HTMHTAnalyzer.cc.

References DQMStore::IBooker::book1D(), LogTrace, metname, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                             {
  evtCounter = 0;
  metname = "HTMHTAnalyzer";

  // PFMET information
  theJetCollectionForHTMHTLabel = parameters.getParameter<edm::InputTag>("JetCollectionForHTMHTLabel");
  _source                       = parameters.getParameter<std::string>("Source");

  LogTrace(metname)<<"[HTMHTAnalyzer] Parameters initialization";
  ibooker.setCurrentFolder("JetMET/MET/"+_source);

  HLTPathsJetMBByName_ = parameters.getParameter<std::vector<std::string > >("HLTPathsJetMB");

  // misc
  _verbose     = parameters.getParameter<int>("verbose");
  _ptThreshold = parameters.getParameter<double>("ptThreshold");

  jetME = ibooker.book1D("metReco", "metReco", 4, 1, 5);
  jetME->setBinLabel(4,"HTMHT",1);

  hNevents = ibooker.book1D("METTask_Nevents",   "METTask_Nevents",   1,0,1);
  hNJets   = ibooker.book1D("METTask_NJets",     "METTask_NJets",     100, 0, 100);
  hMHx     = ibooker.book1D("METTask_MHx",       "METTask_MHx",       500,-500,500);
  hMHy     = ibooker.book1D("METTask_MHy",       "METTask_MHy",       500,-500,500);
  hMHT     = ibooker.book1D("METTask_MHT",       "METTask_MHT",       500,0,1000);
  hMHTPhi  = ibooker.book1D("METTask_MhTPhi",    "METTask_MhTPhi",    80,-4,4);
  hHT      = ibooker.book1D("METTask_HT",        "METTask_HT",        500,0,2000);

}

Member Data Documentation

double HTMHTAnalyzer::_ptThreshold

private

Definition at line 66 of file HTMHTAnalyzer.h.

std::string HTMHTAnalyzer::_source

private

Definition at line 56 of file HTMHTAnalyzer.h.

int HTMHTAnalyzer::_trig_JetMB

private

Definition at line 63 of file HTMHTAnalyzer.h.

int HTMHTAnalyzer::_verbose

private

Definition at line 52 of file HTMHTAnalyzer.h.

Referenced by MassReplace.MassSearchReplaceAnyInputTagVisitor::doIt(), MassReplace.MassSearchReplaceParamVisitor::enter(), cmsPerfSuite.PerfSuite::mkCandleDir(), cmsPerfSuite.PerfSuite::optionParse(), cmsPerfSuite.PerfSuite::printFlush(), cmsPerfSuite.PerfSuite::runcmd(), and cmsPerfSuite.PerfSuite::runCmdSet().

int HTMHTAnalyzer::evtCounter

Definition at line 45 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hHT

private

Definition at line 80 of file HTMHTAnalyzer.h.

std::vector<std::string > HTMHTAnalyzer::HLTPathsJetMBByName_

private

Definition at line 61 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHT

private

Definition at line 77 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHTPhi

private

Definition at line 78 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHx

private

Definition at line 75 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hMHy

private

Definition at line 76 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hNevents

private

Definition at line 71 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::hNJets

private

Definition at line 73 of file HTMHTAnalyzer.h.

MonitorElement* HTMHTAnalyzer::jetME

private

Definition at line 69 of file HTMHTAnalyzer.h.

std::string HTMHTAnalyzer::metname

private

Definition at line 54 of file HTMHTAnalyzer.h.

edm::ParameterSet HTMHTAnalyzer::parameters

private

Definition at line 50 of file HTMHTAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::InputTag HTMHTAnalyzer::theJetCollectionForHTMHTLabel

private

Definition at line 58 of file HTMHTAnalyzer.h.