#include <DQMOffline/ECALMultifitAnalyzer_HI/plugins/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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

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 127 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 &
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 222 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;
 
 }