#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

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

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_

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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 107 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"))),
       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 ( )

inlineoverride

Definition at line 60 of file ECALMultifitAnalyzer_HI.cc.

References analyze(), and bookHistograms().

60 {}

Member Function Documentation

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

overrideprivate

Definition at line 125 of file ECALMultifitAnalyzer_HI.cc.

References edm::SortedCollection< T, SORT >::begin(), caloJetToken_, reco::deltaR(), runTauDisplay::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, MonitorElement::Fill(), relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByToken(), CaloGeometry::getGeometry(), metsig::jet, mDeltaRJetThreshold, mDeltaRPhotonThreshold, mRechitEnergyThreshold, mRecoJetPtThreshold, mRecoPhotonPtThreshold, phi, edm::ESHandle< T >::product(), RecHitCollection_EB_, RecHitCollection_EE_, and recoPhotonsCollection_.

Referenced by ~ECALMultifitAnalyzer_HI().

                                                                                          {
   using namespace edm;
 
   edm::ESHandle<CaloGeometry> geomH;
   iSetup.get<CaloGeometryRecord>().get(geomH);
   const CaloGeometry *geom = geomH.product();
 
   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());
       }
     }
   }
 }

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

overrideprivate

Definition at line 217 of file ECALMultifitAnalyzer_HI.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DEFINE_FWK_MODULE, 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, btvTracks_cfi::maxChi2, mRechitEnergyThreshold, mRecoJetPtThreshold, mRecoPhotonPtThreshold, and DQMStore::IBooker::setCurrentFolder().

Referenced by ~ECALMultifitAnalyzer_HI().

                                                                                                             {
   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;
 }