#include <DQMSourceEleCalib.h>

Inheritance diagram for DQMSourceEleCalib:

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

	~DQMSourceEleCalib () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

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

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

	~DQMEDAnalyzer () override=default

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

Private Member Functions
void	fillAroundBarrel (const EcalRecHitCollection *, int, int)
	fills local occupancy graphs More...

void	fillAroundEndcap (const EcalRecHitCollection *, int, int)

DetId	findMaxHit (const std::vector< std::pair< DetId, float >> &, const EcalRecHitCollection , const EcalRecHitCollection )
	find the MOX More...

Private Attributes
MonitorElement *	ElectronsNumber_
	Number of electrons. More...

MonitorElement *	ESCoP_
	ESCoP. More...

int	eventCounter_

std::string	fileName_
	Output file name if required. More...

std::string	folderName_
	DQM folder name. More...

MonitorElement *	HitsVsAssociatedHits_
	recHits over associated recHits More...

MonitorElement *	LocalOccupancyEB_

MonitorElement *	LocalOccupancyEE_

MonitorElement *	OccupancyEB_
	Occupancy. More...

MonitorElement *	OccupancyEEM_

MonitorElement *	OccupancyEEP_

unsigned int	prescaleFactor_
	Monitor every prescaleFactor_ events. More...

edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEB_
	object to monitor More...

edm::EDGetTokenT< EcalRecHitCollection >	productMonitoredEE_
	object to monitor More...

edm::EDGetTokenT< reco::GsfElectronCollection >	productMonitoredElectrons_
	electrons to monitor More...

MonitorElement *	recHitsPerElectron_
	Number of recHits per electron. More...

bool	saveToFile_
	Write to file. More...

Detailed Description

DQM Source for phi symmetry stream

Author
Stefano Argiro' Andrea Gozzelino - Universita� e INFN Torino

Definition at line 27 of file DQMSourceEleCalib.h.

Constructor & Destructor Documentation

DQMSourceEleCalib::DQMSourceEleCalib ( const edm::ParameterSet & ps )

Definition at line 37 of file DQMSourceEleCalib.cc.

References fileName_, folderName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), prescaleFactor_, productMonitoredEB_, productMonitoredEE_, productMonitoredElectrons_, and saveToFile_.

                                                               : eventCounter_(0) {
   folderName_ = ps.getUntrackedParameter<string>("FolderName", "ALCAStreamEcalSingleEle");
   productMonitoredEB_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("AlCaStreamEBTag"));
   productMonitoredEE_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("AlCaStreamEETag"));
 
   saveToFile_ = ps.getUntrackedParameter<bool>("SaveToFile", false);
   fileName_ = ps.getUntrackedParameter<string>("FileName", "MonitorAlCaEcalSingleEle.root");
   productMonitoredElectrons_ = consumes<reco::GsfElectronCollection>(ps.getParameter<InputTag>("electronCollection"));
   prescaleFactor_ = ps.getUntrackedParameter<unsigned int>("prescaleFactor", 1);
 }

DQMSourceEleCalib::~DQMSourceEleCalib ( )

override

Definition at line 48 of file DQMSourceEleCalib.cc.

48 {}

Member Function Documentation

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

overrideprotected

Definition at line 73 of file DQMSourceEleCalib.cc.

References edm::SortedCollection< T, SORT >::begin(), DetId::det(), EcalBarrel, EcalEndcap, ElectronsNumber_, edm::SortedCollection< T, SORT >::end(), ESCoP_, eventCounter_, MonitorElement::Fill(), fillAroundBarrel(), fillAroundEndcap(), findMaxHit(), edm::Event::getByToken(), HitsVsAssociatedHits_, triggerObjects_cff::id, EBDetId::ieta(), EBDetId::iphi(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), Max(), OccupancyEB_, OccupancyEEM_, OccupancyEEP_, edm::Handle< T >::product(), productMonitoredEB_, productMonitoredEE_, productMonitoredElectrons_, recHitsPerElectron_, edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), and ecaldqm::zside().

                                                                              {
   //  if (eventCounter_% prescaleFactor_ ) return; //FIXME
   eventCounter_++;
   int numberOfHits = 0;
   int numberOfElectrons = 0;
   int numberOfAssociatedHits = 0;
   // reads the recHits
   edm::Handle<EcalRecHitCollection> rhEB;
   edm::Handle<EcalRecHitCollection> rhEE;
 
   iEvent.getByToken(productMonitoredEB_, rhEB);
   iEvent.getByToken(productMonitoredEE_, rhEE);
 
   EcalRecHitCollection::const_iterator itb;
 
   // reads the electrons
   edm::Handle<reco::GsfElectronCollection> pElectrons;
   iEvent.getByToken(productMonitoredElectrons_, pElectrons);
 
   if (pElectrons.isValid()) {
     ElectronsNumber_->Fill(pElectrons->size() + 0.1);
     numberOfElectrons = pElectrons->size();
     for (reco::GsfElectronCollection::const_iterator eleIt = pElectrons->begin(); eleIt != pElectrons->end(); ++eleIt) {
       ESCoP_->Fill(eleIt->eSuperClusterOverP());
       numberOfAssociatedHits += eleIt->superCluster()->size();
       DetId Max = findMaxHit(eleIt->superCluster()->hitsAndFractions(), rhEB.product(), rhEE.product());
       if (!Max.det())
         continue;
       if (Max.subdetId() == EcalBarrel) {
         EBDetId EBMax(Max);
         fillAroundBarrel(rhEB.product(), EBMax.ieta(), EBMax.iphi());
       }
       if (Max.subdetId() == EcalEndcap) {
         EEDetId EEMax(Max);
         fillAroundEndcap(rhEE.product(), EEMax.ix(), EEMax.iy());
       }
     }
   }  // is valid electron
 
   // fill EB histos
   if (rhEB.isValid()) {
     numberOfHits += rhEB->size();
     for (itb = rhEB->begin(); itb != rhEB->end(); ++itb) {
       EBDetId id(itb->id());
       OccupancyEB_->Fill(id.iphi(), id.ieta());
     }  // Eb rechits
   }    // is Valid
   if (rhEE.isValid()) {
     numberOfHits += rhEE->size();
     for (itb = rhEE->begin(); itb != rhEE->end(); ++itb) {
       EEDetId id(itb->id());
       if (id.zside() > 0) {
         OccupancyEEP_->Fill(id.ix(), id.iy());
       }  // zside>0
       else if (id.zside() < 0) {
         OccupancyEEM_->Fill(id.ix(), id.iy());
       }  // zside<0
 
     }  // EE reChit
   }    // is Valid
   if (numberOfElectrons)
     recHitsPerElectron_->Fill((double)numberOfHits / ((double)numberOfElectrons));
   if (numberOfHits)
     HitsVsAssociatedHits_->Fill((double)numberOfAssociatedHits / ((double)numberOfHits));
 }  // end of the analyzer

void DQMSourceEleCalib::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	irun,
		edm::EventSetup const &	isetup
	)

overrideprotected

Definition at line 51 of file DQMSourceEleCalib.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), ElectronsNumber_, ESCoP_, folderName_, HitsVsAssociatedHits_, LocalOccupancyEB_, LocalOccupancyEE_, OccupancyEB_, OccupancyEEM_, OccupancyEEP_, recHitsPerElectron_, and DQMStore::IBooker::setCurrentFolder().

                                                                     {
   // create and cd into new folder
   ibooker.setCurrentFolder(folderName_);
 
   recHitsPerElectron_ = ibooker.book1D("recHitsPerElectron_", "recHitPerElectron", 200, 0, 200);
   ElectronsNumber_ = ibooker.book1D("ElectronsNumber_", "electrons in the event", 40, 0, 40);
   ESCoP_ = ibooker.book1D("ESCoP", "ESCoP", 50, 0, 5);
 
   OccupancyEB_ = ibooker.book2D("OccupancyEB_", "OccupancyEB", 360, 1, 361, 171, -85, 86);
   OccupancyEEP_ = ibooker.book2D("OccupancyEEP_", "Occupancy EE Plus", 100, 1, 101, 100, 1, 101);
   OccupancyEEM_ = ibooker.book2D("OccupancyEEM_", "Occupancy EE Minus", 100, 1, 101, 100, 1, 101);
   HitsVsAssociatedHits_ = ibooker.book1D("HitsVsAssociatedHits_", "HitsVsAssociatedHits", 100, 0, 5);
   LocalOccupancyEB_ = ibooker.book2D("LocalOccupancyEB_", "Local occupancy Barrel", 9, -4, 5, 9, -4, 5);
   LocalOccupancyEE_ = ibooker.book2D("LocalOccupancyEE_", "Local occupancy Endcap", 9, -4, 5, 9, -4, 5);
 }

void DQMSourceEleCalib::fillAroundBarrel	(	const EcalRecHitCollection *	recHits,
		int	eta,
		int	phi
	)

private

fills local occupancy graphs

Definition at line 179 of file DQMSourceEleCalib.cc.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, MonitorElement::Fill(), EBDetId::ieta(), EBDetId::iphi(), LocalOccupancyEB_, and phi.

Referenced by analyze().

                                                                                               {
   for (EcalRecHitCollection::const_iterator elem = recHits->begin(); elem != recHits->end(); ++elem) {
     EBDetId elementId = elem->id();
     LocalOccupancyEB_->Fill(elementId.ieta() - eta, elementId.iphi() - phi, elem->energy());
   }
   return;
 }

void DQMSourceEleCalib::fillAroundEndcap	(	const EcalRecHitCollection *	recHits,
		int	ics,
		int	ips
	)

private

Definition at line 189 of file DQMSourceEleCalib.cc.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), MonitorElement::Fill(), EEDetId::ix(), EEDetId::iy(), and LocalOccupancyEE_.

Referenced by analyze().

                                                                                               {
   for (EcalRecHitCollection::const_iterator elem = recHits->begin(); elem != recHits->end(); ++elem) {
     EEDetId elementId = elem->id();
     LocalOccupancyEE_->Fill(elementId.ix() - ics, elementId.iy() - ips, elem->energy());
   }
   return;
 }

DetId DQMSourceEleCalib::findMaxHit	(	const std::vector< std::pair< DetId, float >> &	v1,
		const EcalRecHitCollection *	EBhits,
		const EcalRecHitCollection *	EEhits
	)

private

find the MOX

Definition at line 143 of file DQMSourceEleCalib.cc.

References EcalBarrel, edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), and MTVHistoProducerAlgoForTrackerBlock_cfi::maxHit.

Referenced by analyze().

                                                                         {
   double currEnergy = 0.;
   DetId maxHit;
   for (std::vector<std::pair<DetId, float>>::const_iterator idsIt = v1.begin(); idsIt != v1.end(); ++idsIt) {
     if (idsIt->first.subdetId() == EcalBarrel) {
       EcalRecHitCollection::const_iterator itrechit;
       itrechit = EBhits->find((*idsIt).first);
       if (itrechit == EBhits->end()) {
         edm::LogInfo("reading") << "[findMaxHit] rechit not found! ";
         continue;
       }
       // FIXME: wnat to use the fraction i.e. .second??
       if (itrechit->energy() > currEnergy) {
         currEnergy = itrechit->energy();
         maxHit = (*idsIt).first;
       }
     } else {
       EcalRecHitCollection::const_iterator itrechit;
       itrechit = EEhits->find((*idsIt).first);
       if (itrechit == EEhits->end()) {
         edm::LogInfo("reading") << "[findMaxHit] rechit not found! ";
         continue;
       }
 
       // FIXME: wnat to use the fraction i.e. .second??
       if (itrechit->energy() > currEnergy) {
         currEnergy = itrechit->energy();
         maxHit = (*idsIt).first;
       }
     }
   }
   return maxHit;
 }