#include <DQMHcalDiJetsAlCaReco.h>

Inheritance diagram for DQMHcalDiJetsAlCaReco:

Public Member Functions
	DQMHcalDiJetsAlCaReco (const edm::ParameterSet &)

	~DQMHcalDiJetsAlCaReco () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Attributes
bool	allowMissingInputs_

edm::EDGetTokenT< EcalRecHitCollection >	ec_

int	eventCounter_

std::string	fileName_
	Output file name if required. More...

std::string	folderName_
	DQM folder name. More...

edm::EDGetTokenT< HBHERecHitCollection >	hbhe_

edm::EDGetTokenT< HFRecHitCollection >	hf_

MonitorElement *	hiDistrEtThirdJet_

MonitorElement *	hiDistrProbeJetEnergy_

MonitorElement *	hiDistrProbeJetEta_

MonitorElement *	hiDistrProbeJetPhi_

MonitorElement *	hiDistrRecHitEnergyEBEE_

MonitorElement *	hiDistrRecHitEnergyHBHE_

MonitorElement *	hiDistrRecHitEnergyHF_

MonitorElement *	hiDistrRecHitEnergyHO_

MonitorElement *	hiDistrTagJetEnergy_

MonitorElement *	hiDistrTagJetEta_

MonitorElement *	hiDistrTagJetPhi_

edm::EDGetTokenT< HORecHitCollection >	ho_

edm::EDGetTokenT< reco::CaloJetCollection >	jets_
	object to monitor More...

bool	saveToFile_
	Write to file. More...

Detailed Description

Definition at line 26 of file DQMHcalDiJetsAlCaReco.h.

Constructor & Destructor Documentation

DQMHcalDiJetsAlCaReco::DQMHcalDiJetsAlCaReco ( const edm::ParameterSet & iConfig )

Definition at line 41 of file DQMHcalDiJetsAlCaReco.cc.

References ec_, fileName_, folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hbhe_, hf_, ho_, jets_, and saveToFile_.

                                                                            : eventCounter_(0) {
   //
   // Input from configurator file
   //
   folderName_ = iConfig.getUntrackedParameter<string>("FolderName", "ALCAStreamHcalDiJets");
 
   jets_ = consumes<CaloJetCollection>(iConfig.getParameter<edm::InputTag>("jetsInput"));
   ec_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ecInput"));
   hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
   ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
   hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInput"));
 
   saveToFile_ = iConfig.getUntrackedParameter<bool>("SaveToFile", false);
   fileName_ = iConfig.getUntrackedParameter<string>("FileName", "MonitorAlCaHcalDiJets.root");
 }

DQMHcalDiJetsAlCaReco::~DQMHcalDiJetsAlCaReco ( )

override

Definition at line 57 of file DQMHcalDiJetsAlCaReco.cc.

57 {}

Member Function Documentation

void DQMHcalDiJetsAlCaReco::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotected

Definition at line 114 of file DQMHcalDiJetsAlCaReco.cc.

References edm::SortedCollection< T, SORT >::begin(), RecoEcal_EventContent_cff::ec, ec_, edm::SortedCollection< T, SORT >::end(), reco::LeafCandidate::energy(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), eventCounter_, MonitorElement::Fill(), edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, hbhe_, photonIsolationHIProducer_cfi::hf, hf_, hiDistrEtThirdJet_, hiDistrProbeJetEnergy_, hiDistrProbeJetEta_, hiDistrProbeJetPhi_, hiDistrRecHitEnergyEBEE_, hiDistrRecHitEnergyHBHE_, hiDistrRecHitEnergyHF_, hiDistrRecHitEnergyHO_, hiDistrTagJetEnergy_, hiDistrTagJetEta_, hiDistrTagJetPhi_, photonIsolationHIProducer_cfi::ho, ho_, edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, jets_, LogDebug, and reco::LeafCandidate::phi().

                                                                                  {
   eventCounter_++;
 
   CaloJet jet1, jet2, jet3;
   Float_t etVetoJet;
 
   edm::Handle<CaloJetCollection> jets;
   iEvent.getByToken(jets_, jets);
 
   if (!jets.isValid()) {
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't getjet product!" << std::endl;
     return;
   }
 
   if (jets->size() > 1) {
     jet1 = (*jets)[0];
     jet2 = (*jets)[1];
     if (fabs(jet1.eta()) > fabs(jet2.eta())) {
       CaloJet jet = jet1;
       jet1 = jet2;
       jet2 = jet;
     }
     //     if(fabs(jet1.eta())>eta_1 || (fabs(jet2.eta())-jet_R) < eta_2){
     //     return;}
   } else {
     return;
   }
 
   hiDistrTagJetEnergy_->Fill(jet1.energy());
   hiDistrTagJetEta_->Fill(jet1.eta());
   hiDistrTagJetPhi_->Fill(jet1.phi());
 
   hiDistrProbeJetEnergy_->Fill(jet2.energy());
   hiDistrProbeJetEta_->Fill(jet2.eta());
   hiDistrProbeJetPhi_->Fill(jet2.phi());
 
   if (jets->size() > 2) {
     jet3 = (*jets)[2];
     etVetoJet = jet3.et();
     hiDistrEtThirdJet_->Fill(etVetoJet);
   } else {
     etVetoJet = 0.;
     hiDistrEtThirdJet_->Fill(etVetoJet);
   }
 
   Handle<EcalRecHitCollection> ec;
   iEvent.getByToken(ec_, ec);
 
   if (!ec.isValid()) {
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get ec product!" << std::endl;
     return;
   }
 
   for (EcalRecHitCollection::const_iterator ecItr = (*ec).begin(); ecItr != (*ec).end(); ++ecItr) {
     hiDistrRecHitEnergyEBEE_->Fill(ecItr->energy());
   }
 
   Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(hbhe_, hbhe);
 
   if (!hbhe.isValid()) {
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get hbhe product!" << std::endl;
     return;
   }
 
   for (HBHERecHitCollection::const_iterator hbheItr = hbhe->begin(); hbheItr != hbhe->end(); hbheItr++) {
     hiDistrRecHitEnergyHBHE_->Fill(hbheItr->energy());
   }
 
   Handle<HORecHitCollection> ho;
   iEvent.getByToken(ho_, ho);
 
   if (!ho.isValid()) {
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get ho product!" << std::endl;
     return;
   }
 
   for (HORecHitCollection::const_iterator hoItr = ho->begin(); hoItr != ho->end(); hoItr++) {
     hiDistrRecHitEnergyHO_->Fill(hoItr->energy());
   }
 
   Handle<HFRecHitCollection> hf;
   iEvent.getByToken(hf_, hf);
 
   if (!hf.isValid()) {
     LogDebug("") << "DQMHcalDiJetsAlCaReco: Error! can't get hf product!" << std::endl;
     return;
   }
 
   for (HFRecHitCollection::const_iterator hfItr = hf->begin(); hfItr != hf->end(); hfItr++) {
     hiDistrRecHitEnergyHF_->Fill(hfItr->energy());
   }
 
 }  // analyze

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

overrideprotected

Definition at line 60 of file DQMHcalDiJetsAlCaReco.cc.

References DQMStore::IBooker::book1D(), folderName_, hiDistrEtThirdJet_, hiDistrProbeJetEnergy_, hiDistrProbeJetEta_, hiDistrProbeJetPhi_, hiDistrRecHitEnergyEBEE_, hiDistrRecHitEnergyHBHE_, hiDistrRecHitEnergyHF_, hiDistrRecHitEnergyHO_, hiDistrTagJetEnergy_, hiDistrTagJetEta_, hiDistrTagJetPhi_, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

                                                                     {
   // create and cd into new folder
   ibooker.setCurrentFolder(folderName_);
 
   // book some histograms 1D
   hiDistrRecHitEnergyEBEE_ = ibooker.book1D("RecHitEnergyEBEE", "the number of hits inside jets", 100, 0, 800);
   hiDistrRecHitEnergyEBEE_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyEBEE_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHBHE_ = ibooker.book1D("RecHitEnergyHBHE", "the number of hits inside jets", 100, 0, 800);
   hiDistrRecHitEnergyHBHE_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHBHE_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHF_ = ibooker.book1D("RecHitEnergyHF", "the number of hits inside jets", 150, 0, 1500);
   hiDistrRecHitEnergyHF_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHF_->setAxisTitle("# rechits", 2);
 
   hiDistrRecHitEnergyHO_ = ibooker.book1D("RecHitEnergyHO", "the number of hits inside jets", 100, 0, 100);
   hiDistrRecHitEnergyHO_->setAxisTitle("E, GeV", 1);
   hiDistrRecHitEnergyHO_->setAxisTitle("# rechits", 2);
 
   hiDistrProbeJetEnergy_ = ibooker.book1D("ProbeJetEnergy", "the energy of probe jets", 250, 0, 2500);
   hiDistrProbeJetEnergy_->setAxisTitle("E, GeV", 1);
   hiDistrProbeJetEnergy_->setAxisTitle("# jets", 2);
 
   hiDistrProbeJetEta_ = ibooker.book1D("ProbeJetEta", "the number of probe jets", 100, -5., 5.);
   hiDistrProbeJetEta_->setAxisTitle("#eta", 1);
   hiDistrProbeJetEta_->setAxisTitle("# jets", 2);
 
   hiDistrProbeJetPhi_ = ibooker.book1D("ProbeJetPhi", "the number of probe jets", 50, -3.14, 3.14);
   hiDistrProbeJetPhi_->setAxisTitle("#phi", 1);
   hiDistrProbeJetPhi_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetEnergy_ = ibooker.book1D("TagJetEnergy", "the energy of tsg jets", 250, 0, 2500);
   hiDistrTagJetEnergy_->setAxisTitle("E, GeV", 1);
   hiDistrTagJetEnergy_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetEta_ = ibooker.book1D("TagJetEta", "the number of  tag jets", 100, -5., 5.);
   hiDistrTagJetEta_->setAxisTitle("#eta", 1);
   hiDistrTagJetEta_->setAxisTitle("# jets", 2);
 
   hiDistrTagJetPhi_ = ibooker.book1D("TagJetPhi", "the number of tag jets", 50, -3.14, 3.14);
   hiDistrTagJetPhi_->setAxisTitle("#phi", 1);
   hiDistrTagJetPhi_->setAxisTitle("# jets", 2);
 
   hiDistrEtThirdJet_ = ibooker.book1D("EtThirdJet", "Et of the third jet", 90, 0, 90);
 
   //==================================================================================
 }