#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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);
  
   //==================================================================================
 }