#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 &event, edm::EventSetup const &setup) final

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

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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 DQMEDAnalyzer
uint64_t	meId () const

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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 37 of file EcalBarrelRecHitsValidation.h.

Constructor & Destructor Documentation

◆ EcalBarrelRecHitsValidation()

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

Constructor.

Definition at line 16 of file EcalBarrelRecHitsValidation.cc.

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

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

◆ ~EcalBarrelRecHitsValidation()

EcalBarrelRecHitsValidation::~EcalBarrelRecHitsValidation ( )

override

Destructor.

Definition at line 47 of file EcalBarrelRecHitsValidation.cc.

47 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 107 of file EcalBarrelRecHitsValidation.cc.

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

References EcalMGPASample::adc(), HltBtagPostValidation_cff::c, MillePedeFileConverter_cfg::e, CollectionTags_cfi::EBDigi, CollectionTags_cfi::EBUncalibRecHit, cosmicPhotonAnalyzer_cfi::eMax, EcalCondObjectContainer< T >::end(), EcalCondObjectContainer< T >::find(), get, 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.

◆ bookHistograms()

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

overrideprotectedvirtual

Analyze.

Implements DQMEDAnalyzer.

Definition at line 49 of file EcalBarrelRecHitsValidation.cc.

                                                                         {
   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.);
   }
 }