#include <ECALMultifitAnalyzer_HI.cc>

Inheritance diagram for ECALMultifitAnalyzer_HI:

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

	~ECALMultifitAnalyzer_HI () 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

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

Private Attributes
edm::EDGetTokenT< reco::CaloJetCollection >	caloJetToken_

MonitorElement *	eb_chi2

MonitorElement *	eb_chi2_e5

MonitorElement *	eb_chi2_e5_eta

MonitorElement *	eb_chi2_eta

MonitorElement *	eb_chi2_jet30

MonitorElement *	eb_chi2_photon15

MonitorElement *	eb_errors

MonitorElement *	eb_errors_e5

MonitorElement *	eb_errors_e5_eta

MonitorElement *	eb_errors_eta

MonitorElement *	eb_errors_jet30

MonitorElement *	eb_errors_photon15

MonitorElement *	ee_chi2

MonitorElement *	ee_chi2_e5

MonitorElement *	ee_chi2_e5_eta

MonitorElement *	ee_chi2_eta

MonitorElement *	ee_chi2_jet30

MonitorElement *	ee_chi2_photon15

MonitorElement *	ee_errors

MonitorElement *	ee_errors_e5

MonitorElement *	ee_errors_e5_eta

MonitorElement *	ee_errors_eta

MonitorElement *	ee_errors_jet30

MonitorElement *	ee_errors_photon15

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomH

double	mDeltaRJetThreshold

double	mDeltaRPhotonThreshold

double	mRechitEnergyThreshold

double	mRecoJetPtThreshold

double	mRecoPhotonPtThreshold

edm::EDGetTokenT< EcalRecHitCollection >	RecHitCollection_EB_

edm::EDGetTokenT< EcalRecHitCollection >	RecHitCollection_EE_

edm::EDGetTokenT< std::vector< reco::Photon > >	recoPhotonsCollection_

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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQMOffline/ECALMultifitAnalyzer_HI/plugins/ECALMultifitAnalyzer_HI.cc Description: [one line class summary] Implementation: [Notes on implementation]

Definition at line 57 of file ECALMultifitAnalyzer_HI.cc.

Constructor & Destructor Documentation

◆ ECALMultifitAnalyzer_HI()

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

explicit

Definition at line 108 of file ECALMultifitAnalyzer_HI.cc.

     : recoPhotonsCollection_(consumes<std::vector<reco::Photon>>(iConfig.getParameter<edm::InputTag>("recoPhotonSrc"))),
       caloJetToken_(consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("recoJetSrc"))),
       RecHitCollection_EB_(consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("RecHitCollection_EB"))),
       RecHitCollection_EE_(consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("RecHitCollection_EE"))),
       geomH(esConsumes()),
       mRechitEnergyThreshold(iConfig.getParameter<double>("rechitEnergyThreshold")),
       mRecoPhotonPtThreshold(iConfig.getParameter<double>("recoPhotonPtThreshold")),
       mRecoJetPtThreshold(iConfig.getParameter<double>("recoJetPtThreshold")),
       mDeltaRPhotonThreshold(iConfig.getParameter<double>("deltaRPhotonThreshold")),
       mDeltaRJetThreshold(iConfig.getParameter<double>("deltaRJetThreshold")) {}

◆ ~ECALMultifitAnalyzer_HI()

ECALMultifitAnalyzer_HI::~ECALMultifitAnalyzer_HI ( )

inlineoverride

Definition at line 60 of file ECALMultifitAnalyzer_HI.cc.

60 {}

Member Function Documentation

◆ analyze()

void ECALMultifitAnalyzer_HI::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 125 of file ECALMultifitAnalyzer_HI.cc.

                                                                                          {
   using namespace edm;
  
   const CaloGeometry *geom = &iSetup.getData(geomH);
  
   Handle<std::vector<reco::Photon>> recoPhotonsHandle;
   iEvent.getByToken(recoPhotonsCollection_, recoPhotonsHandle);
  
   Handle<reco::CaloJetCollection> recoJetHandle;
   iEvent.getByToken(caloJetToken_, recoJetHandle);
  
   Handle<EcalRecHitCollection> ebHandle;
   iEvent.getByToken(RecHitCollection_EB_, ebHandle);
  
   Handle<EcalRecHitCollection> eeHandle;
   iEvent.getByToken(RecHitCollection_EE_, eeHandle);
  
   for (EcalRecHitCollection::const_iterator hit = ebHandle->begin(); hit != ebHandle->end(); ++hit) {
     eb_chi2->Fill(hit->chi2());
     eb_errors->Fill(hit->energyError());
     double eta = geom->getGeometry(hit->detid())->getPosition().eta();
     double phi = geom->getGeometry(hit->detid())->getPosition().phi();
     eb_chi2_eta->Fill(eta, hit->chi2());
     eb_errors_eta->Fill(eta, hit->energyError());
     if (hit->energy() > mRechitEnergyThreshold) {
       eb_chi2_e5->Fill(hit->chi2());
       eb_errors_e5->Fill(hit->energyError());
       eb_chi2_e5_eta->Fill(eta, hit->chi2());
       eb_errors_e5_eta->Fill(eta, hit->energyError());
     }
  
     for (std::vector<reco::Photon>::const_iterator pho = recoPhotonsHandle->begin(); pho != recoPhotonsHandle->end();
          ++pho) {
       if (pho->et() < mRecoPhotonPtThreshold)
         continue;
       double dr = reco::deltaR(eta, phi, pho->eta(), pho->phi());
       if (dr < mDeltaRPhotonThreshold) {
         eb_chi2_photon15->Fill(hit->chi2());
         eb_errors_photon15->Fill(hit->energyError());
       }
     }
     for (std::vector<reco::CaloJet>::const_iterator jet = recoJetHandle->begin(); jet != recoJetHandle->end(); ++jet) {
       if (jet->pt() < mRecoJetPtThreshold)
         continue;
       double dr = reco::deltaR(eta, phi, jet->eta(), jet->phi());
       if (dr < mDeltaRJetThreshold) {
         eb_chi2_jet30->Fill(hit->chi2());
         eb_errors_jet30->Fill(hit->energyError());
       }
     }
   }
  
   for (EcalRecHitCollection::const_iterator hit = eeHandle->begin(); hit != eeHandle->end(); ++hit) {
     ee_chi2->Fill(hit->chi2());
     ee_errors->Fill(hit->energyError());
     double eta = geom->getGeometry(hit->detid())->getPosition().eta();
     double phi = geom->getGeometry(hit->detid())->getPosition().phi();
     ee_chi2_eta->Fill(eta, hit->chi2());
     ee_errors_eta->Fill(eta, hit->energyError());
     if (hit->energy() > mRechitEnergyThreshold) {
       ee_chi2_e5->Fill(hit->chi2());
       ee_errors_e5->Fill(hit->energyError());
       ee_chi2_e5_eta->Fill(eta, hit->chi2());
       ee_errors_e5_eta->Fill(eta, hit->energyError());
     }
  
     for (std::vector<reco::Photon>::const_iterator pho = recoPhotonsHandle->begin(); pho != recoPhotonsHandle->end();
          ++pho) {
       if (pho->et() < mRecoPhotonPtThreshold)
         continue;
       double dr = reco::deltaR(eta, phi, pho->eta(), pho->phi());
       if (dr < mDeltaRPhotonThreshold) {
         ee_chi2_photon15->Fill(hit->chi2());
         ee_errors_photon15->Fill(hit->energyError());
       }
     }
     for (std::vector<reco::CaloJet>::const_iterator jet = recoJetHandle->begin(); jet != recoJetHandle->end(); ++jet) {
       if (jet->pt() < mRecoJetPtThreshold)
         continue;
       double dr = reco::deltaR(eta, phi, jet->eta(), jet->phi());
       if (dr < mDeltaRJetThreshold) {
         ee_chi2_jet30->Fill(hit->chi2());
         ee_errors_jet30->Fill(hit->energyError());
       }
     }
   }
 }

References edm::SortedCollection< T, SORT >::begin(), caloJetToken_, reco::deltaR(), flavorHistoryFilter_cfi::dr, eb_chi2, eb_chi2_e5, eb_chi2_e5_eta, eb_chi2_eta, eb_chi2_jet30, eb_chi2_photon15, eb_errors, eb_errors_e5, eb_errors_e5_eta, eb_errors_eta, eb_errors_jet30, eb_errors_photon15, ee_chi2, ee_chi2_e5, ee_chi2_e5_eta, ee_chi2_eta, ee_chi2_jet30, ee_chi2_photon15, ee_errors, ee_errors_e5, ee_errors_e5_eta, ee_errors_eta, ee_errors_jet30, ee_errors_photon15, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, dqm::impl::MonitorElement::Fill(), relativeConstraints::geom, geomH, edm::EventSetup::getData(), iEvent, metsig::jet, mDeltaRJetThreshold, mDeltaRPhotonThreshold, mRechitEnergyThreshold, mRecoJetPtThreshold, mRecoPhotonPtThreshold, phi, RecHitCollection_EB_, RecHitCollection_EE_, and recoPhotonsCollection_.

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 215 of file ECALMultifitAnalyzer_HI.cc.

                                                                                                             {
   iBooker.setCurrentFolder("EcalCalibration/HIN_MultiFitAnalyzer");
  
   eb_chi2 = nullptr;
   eb_chi2_eta = nullptr;
   eb_chi2_e5 = nullptr;
   eb_chi2_e5_eta = nullptr;
   eb_errors = nullptr;
   eb_errors_eta = nullptr;
   eb_errors_e5 = nullptr;
   eb_errors_e5_eta = nullptr;
   eb_chi2_photon15 = nullptr;
   eb_errors_photon15 = nullptr;
   eb_chi2_jet30 = nullptr;
   eb_errors_jet30 = nullptr;
  
   ee_chi2 = nullptr;
   ee_chi2_eta = nullptr;
   ee_chi2_e5 = nullptr;
   ee_chi2_e5_eta = nullptr;
   ee_errors = nullptr;
   ee_errors_eta = nullptr;
   ee_errors_e5 = nullptr;
   ee_errors_e5_eta = nullptr;
   ee_chi2_photon15 = nullptr;
   ee_errors_photon15 = nullptr;
   ee_chi2_jet30 = nullptr;
   ee_errors_jet30 = nullptr;
  
   const int nBins = 500;
   const float maxChi2 = 70;
   const float maxError = 0.5;
  
   TH2F *hProfile_Chi2 = new TH2F("hProfile_Chi2", "", nBins, -5, 5, nBins, 0, maxChi2);
   TH2F *hProfile_Err = new TH2F("hProfile_Err", "", nBins, -5, 5, nBins, 0, maxError);
  
   eb_chi2 = iBooker.book1D("rechit_eb_chi2", "Rechit eb_chi2;chi2 fit value;", nBins, 0, maxChi2);
   eb_chi2_eta = iBooker.book2D("rechit_eb_eta_Vs_mean_Chi2", hProfile_Chi2);
   eb_chi2_e5 = iBooker.book1D(Form("rechit_eb_chi2_e%2.0f", mRechitEnergyThreshold),
                               Form("Rechit eb_chi2, e>%2.0fGeV;chi2 fit value;", mRechitEnergyThreshold),
                               nBins,
                               0,
                               maxChi2);
   eb_chi2_e5_eta = iBooker.book2D(Form("rechit_eb_chi2_e%2.0f_eta", mRechitEnergyThreshold), hProfile_Chi2);
   eb_errors = iBooker.book1D("rechit_eb_errors", "Rechit eb_errors;error on the energy;", nBins, 0, maxError);
   eb_errors_eta = iBooker.book2D("rechit_eb_errors_eta", hProfile_Err);
   eb_errors_e5 = iBooker.book1D(Form("rechit_eb_errors_e%2.0f", mRechitEnergyThreshold),
                                 "Rechit eb_errors, e>5GeV;error on the energy;",
                                 nBins,
                                 0,
                                 maxError);
   eb_errors_e5_eta = iBooker.book2D(Form("rechit_eb_errors_e%2.0f_eta", mRechitEnergyThreshold), hProfile_Err);
   eb_chi2_photon15 = iBooker.book1D(Form("rechit_eb_chi2_photon%2.0f", mRecoPhotonPtThreshold),
                                     "Rechit eb_chi2 near photons;chi2 fit value;",
                                     nBins,
                                     0,
                                     maxChi2);
   eb_errors_photon15 = iBooker.book1D(Form("rechit_eb_errors_photon%2.0f", mRecoPhotonPtThreshold),
                                       "Rechit eb_errors near photons;error on the energy;",
                                       nBins,
                                       0,
                                       maxError);
   eb_chi2_jet30 = iBooker.book1D(Form("rechit_eb_chi2_jet%2.0f", mRecoJetPtThreshold),
                                  "Rechit eb_chi2 near jets;chi2 fit value;",
                                  nBins,
                                  0,
                                  maxChi2);
   eb_errors_jet30 = iBooker.book1D(Form("rechit_eb_errors_jet%2.0f", mRecoJetPtThreshold),
                                    "Rechit eb_errors near jets;error on the energy;",
                                    nBins,
                                    0,
                                    maxError);
  
   ee_chi2 = iBooker.book1D("rechit_ee_chi2", "Rechit ee_chi2;chi2 fit value;", nBins, 0, maxChi2);
   ee_chi2_eta = iBooker.book2D("rechit_ee_chi2_eta", hProfile_Chi2);
   ee_chi2_e5 = iBooker.book1D(Form("rechit_ee_chi2_e%2.0f", mRechitEnergyThreshold),
                               "Rechit ee_chi2, e>5GeV;chi2 fit value;",
                               nBins,
                               0,
                               maxChi2);
   ee_chi2_e5_eta = iBooker.book2D(Form("rechit_ee_chi2_e%2.0f_eta", mRechitEnergyThreshold), hProfile_Chi2);
   ee_errors = iBooker.book1D("rechit_ee_errors", "Rechit ee_errors;error on the energy;", nBins, 0, maxError);
   ee_errors_eta = iBooker.book2D("rechit_ee_errors_eta", hProfile_Err);
   ee_errors_e5 = iBooker.book1D(Form("rechit_ee_errors_e%2.0f", mRechitEnergyThreshold),
                                 "Rechit ee_errors, e>5GeV;error on the energy;",
                                 nBins,
                                 0,
                                 maxError);
   ee_errors_e5_eta = iBooker.book2D(Form("rechit_ee_errors_e%2.0f_eta", mRechitEnergyThreshold), hProfile_Err);
   ee_chi2_photon15 = iBooker.book1D(Form("rechit_ee_chi2_photon%2.0f", mRecoPhotonPtThreshold),
                                     "Rechit ee_chi2 near photons;chi2 fit value;",
                                     nBins,
                                     0,
                                     maxChi2);
   ee_errors_photon15 = iBooker.book1D(Form("rechit_ee_errors_photon%2.0f", mRecoPhotonPtThreshold),
                                       "Rechit ee_errors near photons;error on the energy;",
                                       nBins,
                                       0,
                                       maxError);
   ee_chi2_jet30 = iBooker.book1D(Form("rechit_ee_chi2_jet%2.0f", mRecoJetPtThreshold),
                                  "Rechit ee_chi2 near jets;chi2 fit value;",
                                  nBins,
                                  0,
                                  maxChi2);
   ee_errors_jet30 = iBooker.book1D(Form("rechit_ee_errors_jet%2.0f", mRecoJetPtThreshold),
                                    "Rechit ee_errors near jets;error on the energy;",
                                    nBins,
                                    0,
                                    maxError);
  
   delete hProfile_Chi2;
   delete hProfile_Err;
 }