#include <EcalBarrelRecHitsValidation.h>

Inheritance diagram for EcalBarrelRecHitsValidation:

Public Member Functions
	EcalBarrelRecHitsValidation (const edm::ParameterSet &ps)
	Constructor. More...

	~EcalBarrelRecHitsValidation () override
	Destructor. More...

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 &i, edm::Run const &, edm::EventSetup const &) override
	Analyze. More...

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 Attributes
edm::EDGetTokenT< EBDigiCollection >	EBdigiCollection_token_

edm::EDGetTokenT< EBUncalibratedRecHitCollection >	EBuncalibrechitCollection_token_

MonitorElement *	meEBUncalibRecHitAmplMap_ [36]

MonitorElement *	meEBUncalibRecHitMaxSampleRatio_

MonitorElement *	meEBUncalibRecHitMaxSampleRatioGt100adc_

MonitorElement *	meEBUncalibRecHitPedMap_ [36]

MonitorElement *	meEBUncalibRecHitsAmpFullMap_

MonitorElement *	meEBUncalibRecHitsAmplitude_

MonitorElement *	meEBUncalibRecHitsAmplitudeGt100adc_

MonitorElement *	meEBUncalibRecHitsChi2_

MonitorElement *	meEBUncalibRecHitsChi2Gt100adc_

MonitorElement *	meEBUncalibRecHitsJitter_

MonitorElement *	meEBUncalibRecHitsJitterGt100adc_

MonitorElement *	meEBUncalibRecHitsOccupancy_

MonitorElement *	meEBUncalibRecHitsOccupancyGt100adc_

MonitorElement *	meEBUncalibRecHitsPedestal_

MonitorElement *	meEBUncalibRecHitsPedestalGt100adc_

MonitorElement *	meEBUncalibRecHitsPedFullMap_

bool	verbose_

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

Definition at line 38 of file EcalBarrelRecHitsValidation.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 16 of file EcalBarrelRecHitsValidation.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and mps_fire::i.

                                                                               {
   
   // ---------------------- 
   EBdigiCollection_token_          = consumes< EBDigiCollection > (ps.getParameter<edm::InputTag>("EBdigiCollection") );
   EBuncalibrechitCollection_token_ = consumes< EBUncalibratedRecHitCollection > (ps.getParameter<edm::InputTag>("EBuncalibrechitCollection") );
   
   // ---------------------- 
   // verbosity switch 
   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
     
   // ----------------------   
   meEBUncalibRecHitsOccupancy_             = nullptr;
   meEBUncalibRecHitsAmplitude_             = nullptr;
   meEBUncalibRecHitsPedestal_              = nullptr;
   meEBUncalibRecHitsJitter_                = nullptr;
   meEBUncalibRecHitsChi2_                  = nullptr;
   meEBUncalibRecHitMaxSampleRatio_         = nullptr;
   meEBUncalibRecHitsOccupancyGt100adc_     = nullptr;
   meEBUncalibRecHitsAmplitudeGt100adc_     = nullptr;
   meEBUncalibRecHitsPedestalGt100adc_      = nullptr;
   meEBUncalibRecHitsJitterGt100adc_        = nullptr;
   meEBUncalibRecHitsChi2Gt100adc_          = nullptr;
   meEBUncalibRecHitMaxSampleRatioGt100adc_ = nullptr;
   meEBUncalibRecHitsAmpFullMap_            = nullptr;
   meEBUncalibRecHitsPedFullMap_            = nullptr;
   for (int i=0; i<36 ; i++) 
     {
       meEBUncalibRecHitAmplMap_[i] = nullptr;
       meEBUncalibRecHitPedMap_[i]  = nullptr;
     }
 }

EcalBarrelRecHitsValidation::~EcalBarrelRecHitsValidation ( )

override

Destructor.

Definition at line 48 of file EcalBarrelRecHitsValidation.cc.

                                                          {   
 
 }

Member Function Documentation

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

overrideprotected

Definition at line 109 of file EcalBarrelRecHitsValidation.cc.

References EcalMGPASample::adc(), edm::SortedCollection< T, SORT >::begin(), photonPostprocessing_cfi::eMax, EcalCondObjectContainer< T >::end(), edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), HcalObjRepresent::Fill(), EcalCondObjectContainer< T >::find(), edm::DataFrameContainer::find(), edm::EventSetup::get(), edm::Event::getByToken(), EcalCondObjectContainer< T >::getMap(), EBDetId::ic(), EBDetId::ieta(), EBDetId::iphi(), ecalpyutils::ism(), EBDetId::ism(), edm::HandleBase::isValid(), LogDebug, edm::Handle< T >::product(), edm::ESHandle< T >::product(), and simplePhotonAnalyzer_cfi::sample.

                                                                             {
 
   
   const EBUncalibratedRecHitCollection *EBUncalibRecHit = nullptr;
   Handle< EBUncalibratedRecHitCollection > EcalUncalibRecHitEB;
   e.getByToken( EBuncalibrechitCollection_token_, EcalUncalibRecHitEB);
   if (EcalUncalibRecHitEB.isValid()) {
     EBUncalibRecHit = EcalUncalibRecHitEB.product();
   } else {
     return;
   }
 
   bool skipDigis = false;
   const EBDigiCollection *EBDigi = nullptr;
   Handle< EBDigiCollection > EcalDigiEB;
   e.getByToken( EBdigiCollection_token_, EcalDigiEB);
   if (EcalDigiEB.isValid()) {
     EBDigi = EcalDigiEB.product();    
   } else {
     skipDigis = true;
   }
 
   edm::ESHandle<EcalPedestals> ecalPeds; 
   c.get<EcalPedestalsRcd>().get(ecalPeds);
 
   // ---------------------- 
   // loop over UncalibRecHits
   for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EBUncalibRecHit->begin(); uncalibRecHit != EBUncalibRecHit->end() ; ++uncalibRecHit)
     {
       EBDetId EBid = EBDetId(uncalibRecHit->id());
       
       // general checks
       if (meEBUncalibRecHitsOccupancy_)  meEBUncalibRecHitsOccupancy_  -> Fill(EBid.ieta(), EBid.iphi());
       if (meEBUncalibRecHitsAmplitude_)  meEBUncalibRecHitsAmplitude_  -> Fill(uncalibRecHit->amplitude());
       if (meEBUncalibRecHitsPedestal_)   meEBUncalibRecHitsPedestal_   -> Fill(uncalibRecHit->pedestal());
       if (meEBUncalibRecHitsJitter_)     meEBUncalibRecHitsJitter_     -> Fill(uncalibRecHit->jitter());
       if (meEBUncalibRecHitsChi2_)       meEBUncalibRecHitsChi2_       -> Fill(uncalibRecHit->chi2());
       if (meEBUncalibRecHitsAmpFullMap_) meEBUncalibRecHitsAmpFullMap_ -> Fill(EBid.ieta(), EBid.iphi(), uncalibRecHit->amplitude()); 
       if (meEBUncalibRecHitsPedFullMap_) meEBUncalibRecHitsPedFullMap_ -> Fill(EBid.ieta(), EBid.iphi(), uncalibRecHit->pedestal()); 
 
       // general checks, with threshold at 3.5 GeV = 100 ADC counts
       if ( uncalibRecHit->amplitude() > 100 ) 
     {
       if (meEBUncalibRecHitsOccupancyGt100adc_)  meEBUncalibRecHitsOccupancyGt100adc_  -> Fill(EBid.ieta(), EBid.iphi());
       if (meEBUncalibRecHitsAmplitudeGt100adc_)  meEBUncalibRecHitsAmplitudeGt100adc_  -> Fill(uncalibRecHit->amplitude());
       if (meEBUncalibRecHitsPedestalGt100adc_)   meEBUncalibRecHitsPedestalGt100adc_   -> Fill(uncalibRecHit->pedestal());
       if (meEBUncalibRecHitsJitterGt100adc_)     meEBUncalibRecHitsJitterGt100adc_     -> Fill(uncalibRecHit->jitter());
       if (meEBUncalibRecHitsChi2Gt100adc_)       meEBUncalibRecHitsChi2Gt100adc_       -> Fill(uncalibRecHit->chi2());
     }
 
       // supermodule maps
       int ic = EBid.ic();              
       int ie = (ic-1)/20 + 1;          
       int ip = (ic-1)%20 + 1;
       int ism = EBid.ism();
       float xie = ie - 0.5;
       float xip = ip - 0.5;      
       if( meEBUncalibRecHitPedMap_[ism-1] ) meEBUncalibRecHitPedMap_[ism-1]->Fill(xie, xip, uncalibRecHit->pedestal());
       if( meEBUncalibRecHitAmplMap_[ism-1] ) meEBUncalibRecHitAmplMap_[ism-1]->Fill(xie, xip, uncalibRecHit->amplitude());
       
       if ( ! skipDigis ) { 
         // find the rechit corresponding digi and the max sample
         EBDigiCollection::const_iterator myDigi = EBDigi->find(EBid);
         // int sMax = -1; // UNUSED
         double eMax = 0.;
         if (myDigi != EBDigi->end()){
           for (unsigned int sample = 0 ; sample < myDigi->size(); ++sample){
             EcalMGPASample thisSample = (*myDigi)[sample];
             double analogSample = thisSample.adc();
             if ( eMax < analogSample ){
               eMax = analogSample;
               // sMax = sample; // UNUSED
             }
           }
         }
         else
           continue;
 
         // ratio uncalibratedRecHit amplitude + ped / max energy digi  
         const EcalPedestals* myped = ecalPeds.product();
         EcalPedestalsMap::const_iterator it=myped->getMap().find( EBid );
         if( it != myped->getMap().end() ){
           
           if (eMax > (*it).mean_x1 + 5 * (*it).rms_x1 && eMax != 0 ) {//only real signal RecHit
             
             if ( meEBUncalibRecHitMaxSampleRatio_ ) meEBUncalibRecHitMaxSampleRatio_->Fill( (uncalibRecHit->amplitude()+uncalibRecHit->pedestal())/eMax);
             if ( meEBUncalibRecHitMaxSampleRatioGt100adc_ && (uncalibRecHit->amplitude()>100) ) meEBUncalibRecHitMaxSampleRatioGt100adc_->Fill( (uncalibRecHit->amplitude()+uncalibRecHit->pedestal())/eMax);
             LogDebug("EcalRecHitsTaskInfo") << "barrel, eMax = " << eMax << " Amplitude = " << uncalibRecHit->amplitude()+uncalibRecHit->pedestal();  
           }
           else
             continue;
         }
         else
           continue;
       }
       
     }  // loop over the UncalibratedRecHitCollection
 }

void EcalBarrelRecHitsValidation::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Analyze.

Implements DQMEDAnalyzer.

Definition at line 52 of file EcalBarrelRecHitsValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile2D(), trackerHits::histo, mps_fire::i, and DQMStore::IBooker::setCurrentFolder().

                                                                                                              { 
   Char_t histo[200];
    
   ibooker.setCurrentFolder("EcalRecHitsV/EcalBarrelRecHitsTask");
    
   sprintf (histo, "EB Occupancy" );  
   meEBUncalibRecHitsOccupancy_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
   
   sprintf (histo, "EB Amplitude" );
   meEBUncalibRecHitsAmplitude_ = ibooker.book1D(histo, histo, 201, -20., 4000.);
   
   sprintf (histo, "EB Pedestal" );
   meEBUncalibRecHitsPedestal_ = ibooker.book1D(histo, histo, 50, 190., 210.);
   
   sprintf (histo, "EB Jitter" );
   meEBUncalibRecHitsJitter_ = ibooker.book1D(histo, histo, 100, 0., 100.);
       
   sprintf (histo, "EB Chi2" );
   meEBUncalibRecHitsChi2_ = ibooker.book1D(histo, histo, 100, 18000., 22000.);
 
   sprintf (histo, "EB RecHit Max Sample Ratio"); 
   meEBUncalibRecHitMaxSampleRatio_ = ibooker.book1D(histo, histo, 120, 0.90, 1.05);
 
   sprintf (histo, "EB Occupancy gt 100 adc counts" );
   meEBUncalibRecHitsOccupancyGt100adc_ = ibooker.book2D(histo, histo, 170, -85., 85., 360, 0., 360.);
   
   sprintf (histo, "EB Amplitude gt 100 adc counts" );
   meEBUncalibRecHitsAmplitudeGt100adc_ = ibooker.book1D(histo, histo, 200, 0., 4000.);
 
   sprintf (histo, "EB Pedestal gt 100 adc counts" );
   meEBUncalibRecHitsPedestalGt100adc_ = ibooker.book1D(histo, histo, 50, 190., 210.);
 
   sprintf (histo, "EB Jitter gt 100 adc counts" );
   meEBUncalibRecHitsJitterGt100adc_ = ibooker.book1D(histo, histo, 100, 0., 100.);
 
   sprintf (histo, "EB Chi2 gt 100 adc counts" );
   meEBUncalibRecHitsChi2Gt100adc_ = ibooker.book1D(histo, histo, 100, 18000., 22000.);
 
   sprintf (histo, "EB RecHit Max Sample Ratio gt 100 adc counts"); 
   meEBUncalibRecHitMaxSampleRatioGt100adc_ = ibooker.book1D(histo, histo, 120, 0.90, 1.05);
   
   sprintf (histo, "EB Amplitude Full Map");
   meEBUncalibRecHitsAmpFullMap_ = ibooker.bookProfile2D(histo, histo, 170, -85., 85., 360, 0., 360., 200, 0., 4000.);
 
   sprintf (histo, "EB Pedestal Full Map");
   meEBUncalibRecHitsPedFullMap_ = ibooker.bookProfile2D(histo, histo, 170, -85., 85., 360, 0., 360., 50, 194., 201.);
 
   for (int i=0; i<36 ; i++) 
     {
       sprintf(histo, "EB Amp SM%02d", i+1);
       meEBUncalibRecHitAmplMap_[i] = ibooker.bookProfile2D(histo, histo, 85, 0., 85., 20, 0., 20., 200, 0., 4000.);
       
       sprintf(histo, "EB Ped SM%02d", i+1);
       meEBUncalibRecHitPedMap_[i] = ibooker.bookProfile2D(histo, histo, 85, 0., 85., 20, 0., 20., 50, 194., 201.);
     }
 }