#include <DQMSourceEleCalib.h>

Inheritance diagram for DQMSourceEleCalib:

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

	~DQMSourceEleCalib () 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 () noexcept(false) 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)

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

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...

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM Source for phi symmetry stream

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

Definition at line 30 of file DQMSourceEleCalib.h.

Constructor & Destructor Documentation

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

Definition at line 40 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 55 of file DQMSourceEleCalib.cc.

56 {}

Member Function Documentation

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

overrideprotected

Definition at line 84 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
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 60 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 214 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_, phi, and mathSSE::return().

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 236 of file DQMSourceEleCalib.cc.

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

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