#include <DQMHcalDiJetsAlCaReco.h>

Inheritance diagram for DQMHcalDiJetsAlCaReco:

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

	~DQMHcalDiJetsAlCaReco () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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...

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 24 of file DQMHcalDiJetsAlCaReco.h.

Constructor & Destructor Documentation

◆ DQMHcalDiJetsAlCaReco()

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

Definition at line 41 of file DQMHcalDiJetsAlCaReco.cc.

                                                                            : 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");
 }

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

◆ ~DQMHcalDiJetsAlCaReco()

DQMHcalDiJetsAlCaReco::~DQMHcalDiJetsAlCaReco ( )

override

Definition at line 57 of file DQMHcalDiJetsAlCaReco.cc.

57 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 114 of file DQMHcalDiJetsAlCaReco.cc.

                                                                                  {
   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

References ec_, reco::LeafCandidate::energy(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), eventCounter_, dqm::impl::MonitorElement::Fill(), photonIsolationHIProducer_cfi::hbhe, hbhe_, photonIsolationHIProducer_cfi::hf, hf_, hiDistrEtThirdJet_, hiDistrProbeJetEnergy_, hiDistrProbeJetEta_, hiDistrProbeJetPhi_, hiDistrRecHitEnergyEBEE_, hiDistrRecHitEnergyHBHE_, hiDistrRecHitEnergyHF_, hiDistrRecHitEnergyHO_, hiDistrTagJetEnergy_, hiDistrTagJetEta_, hiDistrTagJetPhi_, photonIsolationHIProducer_cfi::ho, ho_, iEvent, edm::HandleBase::isValid(), metsig::jet, singleTopDQM_cfi::jets, jets_, LogDebug, and reco::LeafCandidate::phi().

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 60 of file DQMHcalDiJetsAlCaReco.cc.

                                                                     {
   // 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);
  
   //==================================================================================
 }