#include <DQMSourceEleCalib.h>

Inheritance diagram for DQMSourceEleCalib:

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

	~DQMSourceEleCalib ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)

void	beginJob ()

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context)

void	beginRun (const edm::Run &r, const edm::EventSetup &c)

void	endJob ()

void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

void	endRun (const edm::Run &r, const edm::EventSetup &c)

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
DQMStore *	dbe_

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

DQM Source for phi symmetry stream

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

Definition at line 26 of file DQMSourceEleCalib.h.

Constructor & Destructor Documentation

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

Definition at line 40 of file DQMSourceEleCalib.cc.

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

                                                                 :
 eventCounter_(0)
 {
   dbe_ = Service<DQMStore>().operator->();
   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 ( )

Definition at line 56 of file DQMSourceEleCalib.cc.

57 {}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 94 of file DQMSourceEleCalib.cc.

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

                                              {
  
 //  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::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file DQMSourceEleCalib.cc.

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

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

void DQMSourceEleCalib::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	context
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 87 of file DQMSourceEleCalib.cc.

                                 {
   
 }

void DQMSourceEleCalib::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 82 of file DQMSourceEleCalib.cc.

                                                                            {
 
 }

void DQMSourceEleCalib::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file DQMSourceEleCalib.cc.

References dbe_, fileName_, DQMStore::save(), and saveToFile_.

                               {
   
   if (saveToFile_) {
      dbe_->save(fileName_);
   }
   
 }

void DQMSourceEleCalib::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	c
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 167 of file DQMSourceEleCalib.cc.

168 {

169 }

void DQMSourceEleCalib::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 171 of file DQMSourceEleCalib.cc.

                                                                      {
 
 }

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

private

fills local occupancy graphs

Definition at line 237 of file DQMSourceEleCalib.cc.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), eta(), MonitorElement::Fill(), EBDetId::ieta(), EBDetId::iphi(), LocalOccupancyEB_, phi, and reco::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 259 of file DQMSourceEleCalib.cc.

References edm::SortedCollection< T, SORT >::begin(), HTMLExport::elem(), edm::SortedCollection< T, SORT >::end(), MonitorElement::Fill(), h2_mapping_cfi::ics, EEDetId::ix(), EEDetId::iy(), LocalOccupancyEE_, and reco::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 185 of file DQMSourceEleCalib.cc.

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

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

Member Data Documentation

DQMStore* DQMSourceEleCalib::dbe_

private

Definition at line 62 of file DQMSourceEleCalib.h.

Referenced by beginJob(), DQMSourceEleCalib(), and endJob().

MonitorElement* DQMSourceEleCalib::ElectronsNumber_

private

Number of electrons.

Definition at line 68 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

MonitorElement* DQMSourceEleCalib::ESCoP_

private

ESCoP.

Definition at line 70 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

int DQMSourceEleCalib::eventCounter_

private

Definition at line 63 of file DQMSourceEleCalib.h.

Referenced by analyze().

std::string DQMSourceEleCalib::fileName_

private

Output file name if required.

Definition at line 99 of file DQMSourceEleCalib.h.

Referenced by DQMSourceEleCalib(), and endJob().

std::string DQMSourceEleCalib::folderName_

private

DQM folder name.

Definition at line 93 of file DQMSourceEleCalib.h.

Referenced by beginJob(), and DQMSourceEleCalib().

MonitorElement* DQMSourceEleCalib::HitsVsAssociatedHits_

private

recHits over associated recHits

Definition at line 79 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

MonitorElement* DQMSourceEleCalib::LocalOccupancyEB_

private

Definition at line 75 of file DQMSourceEleCalib.h.

Referenced by beginJob(), and fillAroundBarrel().

MonitorElement* DQMSourceEleCalib::LocalOccupancyEE_

private

Definition at line 76 of file DQMSourceEleCalib.h.

Referenced by beginJob(), and fillAroundEndcap().

MonitorElement* DQMSourceEleCalib::OccupancyEB_

private

Occupancy.

Definition at line 72 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

MonitorElement* DQMSourceEleCalib::OccupancyEEM_

private

Definition at line 74 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

MonitorElement* DQMSourceEleCalib::OccupancyEEP_

private

Definition at line 73 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

unsigned int DQMSourceEleCalib::prescaleFactor_

private

Monitor every prescaleFactor_ events.

Definition at line 90 of file DQMSourceEleCalib.h.

Referenced by DQMSourceEleCalib().

edm::EDGetTokenT<EcalRecHitCollection> DQMSourceEleCalib::productMonitoredEB_

private

object to monitor

Definition at line 82 of file DQMSourceEleCalib.h.

Referenced by analyze(), and DQMSourceEleCalib().

edm::EDGetTokenT<EcalRecHitCollection> DQMSourceEleCalib::productMonitoredEE_

private

object to monitor

Definition at line 85 of file DQMSourceEleCalib.h.

Referenced by analyze(), and DQMSourceEleCalib().

edm::EDGetTokenT<reco::GsfElectronCollection> DQMSourceEleCalib::productMonitoredElectrons_

private

electrons to monitor

Definition at line 87 of file DQMSourceEleCalib.h.

Referenced by analyze(), and DQMSourceEleCalib().

MonitorElement* DQMSourceEleCalib::recHitsPerElectron_

private

Number of recHits per electron.

Definition at line 66 of file DQMSourceEleCalib.h.

Referenced by analyze(), and beginJob().

bool DQMSourceEleCalib::saveToFile_

private

Write to file.

Definition at line 96 of file DQMSourceEleCalib.h.

Referenced by DQMSourceEleCalib(), and endJob().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation