#include <EcalEndcapRecHitsValidation.h>

Inheritance diagram for EcalEndcapRecHitsValidation:

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

	~EcalEndcapRecHitsValidation () 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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
edm::ESGetToken< EcalPedestals, EcalPedestalsRcd >	ecalPeds

edm::EDGetTokenT< EEDigiCollection >	EEdigiCollection_token_

edm::EDGetTokenT< EEUncalibratedRecHitCollection >	EEuncalibrechitCollection_token_

MonitorElement *	meEEUncalibRecHitMaxSampleRatio_

MonitorElement *	meEEUncalibRecHitMaxSampleRatioGt60adc_

MonitorElement *	meEEUncalibRecHitsAmpFullMap_

MonitorElement *	meEEUncalibRecHitsAmplitude_

MonitorElement *	meEEUncalibRecHitsAmplitudeGt60adc_

MonitorElement *	meEEUncalibRecHitsChi2_

MonitorElement *	meEEUncalibRecHitsChi2Gt60adc_

MonitorElement *	meEEUncalibRecHitsJitter_

MonitorElement *	meEEUncalibRecHitsJitterGt60adc_

MonitorElement *	meEEUncalibRecHitsOccupancyMinus_

MonitorElement *	meEEUncalibRecHitsOccupancyMinusGt60adc_

MonitorElement *	meEEUncalibRecHitsOccupancyPlus_

MonitorElement *	meEEUncalibRecHitsOccupancyPlusGt60adc_

MonitorElement *	meEEUncalibRecHitsPedestal_

MonitorElement *	meEEUncalibRecHitsPedestalGt60adc_

MonitorElement *	meEEUncalibRecHitsPedFullMap_

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 EcalEndcapRecHitsValidation.h.

Constructor & Destructor Documentation

◆ EcalEndcapRecHitsValidation()

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

Constructor.

Definition at line 16 of file EcalEndcapRecHitsValidation.cc.

                                                                                : ecalPeds(esConsumes()) {
   // ----------------------
   EEdigiCollection_token_ = consumes<EEDigiCollection>(ps.getParameter<edm::InputTag>("EEdigiCollection"));
   EEuncalibrechitCollection_token_ =
       consumes<EEUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("EEuncalibrechitCollection"));
  
   // ----------------------
   // verbosity switch
   verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
  
   // ----------------------
   meEEUncalibRecHitsOccupancyPlus_ = nullptr;
   meEEUncalibRecHitsOccupancyMinus_ = nullptr;
   meEEUncalibRecHitsAmplitude_ = nullptr;
   meEEUncalibRecHitsPedestal_ = nullptr;
   meEEUncalibRecHitsJitter_ = nullptr;
   meEEUncalibRecHitsChi2_ = nullptr;
   meEEUncalibRecHitMaxSampleRatio_ = nullptr;
   meEEUncalibRecHitsOccupancyPlusGt60adc_ = nullptr;
   meEEUncalibRecHitsOccupancyMinusGt60adc_ = nullptr;
   meEEUncalibRecHitsAmplitudeGt60adc_ = nullptr;
   meEEUncalibRecHitsPedestalGt60adc_ = nullptr;
   meEEUncalibRecHitsJitterGt60adc_ = nullptr;
   meEEUncalibRecHitsChi2Gt60adc_ = nullptr;
   meEEUncalibRecHitMaxSampleRatioGt60adc_ = nullptr;
   meEEUncalibRecHitsAmpFullMap_ = nullptr;
   meEEUncalibRecHitsPedFullMap_ = nullptr;
 }

◆ ~EcalEndcapRecHitsValidation()

EcalEndcapRecHitsValidation::~EcalEndcapRecHitsValidation ( )

override

Destructor.

Definition at line 45 of file EcalEndcapRecHitsValidation.cc.

45 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 103 of file EcalEndcapRecHitsValidation.cc.

                                                                              {
   const EEUncalibratedRecHitCollection *EEUncalibRecHit = nullptr;
   Handle<EEUncalibratedRecHitCollection> EcalUncalibRecHitEE;
   e.getByToken(EEuncalibrechitCollection_token_, EcalUncalibRecHitEE);
   if (EcalUncalibRecHitEE.isValid()) {
     EEUncalibRecHit = EcalUncalibRecHitEE.product();
   } else {
     return;
   }
  
   bool skipDigis = false;
   const EEDigiCollection *EEDigi = nullptr;
   Handle<EEDigiCollection> EcalDigiEE;
   e.getByToken(EEdigiCollection_token_, EcalDigiEE);
   if (EcalDigiEE.isValid()) {
     EEDigi = EcalDigiEE.product();
   } else {
     skipDigis = true;
   }
  
   // ----------------------
   // loop over UncalibRecHits
   for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit = EEUncalibRecHit->begin();
        uncalibRecHit != EEUncalibRecHit->end();
        ++uncalibRecHit) {
     EEDetId EEid = EEDetId(uncalibRecHit->id());
  
     int mySide = EEid.zside();
  
     // general checks
     if (mySide > 0) {
       if (meEEUncalibRecHitsOccupancyPlus_)
         meEEUncalibRecHitsOccupancyPlus_->Fill(EEid.ix(), EEid.iy());
     }
     if (mySide < 0) {
       if (meEEUncalibRecHitsOccupancyMinus_)
         meEEUncalibRecHitsOccupancyMinus_->Fill(EEid.ix(), EEid.iy());
     }
     if (meEEUncalibRecHitsAmplitude_)
       meEEUncalibRecHitsAmplitude_->Fill(uncalibRecHit->amplitude());
     if (meEEUncalibRecHitsPedestal_)
       meEEUncalibRecHitsPedestal_->Fill(uncalibRecHit->pedestal());
     if (meEEUncalibRecHitsJitter_)
       meEEUncalibRecHitsJitter_->Fill(uncalibRecHit->jitter());
     if (meEEUncalibRecHitsChi2_)
       meEEUncalibRecHitsChi2_->Fill(uncalibRecHit->chi2());
     if (meEEUncalibRecHitsAmpFullMap_)
       meEEUncalibRecHitsAmpFullMap_->Fill(EEid.ix(), EEid.iy(), uncalibRecHit->amplitude());
     if (meEEUncalibRecHitsPedFullMap_)
       meEEUncalibRecHitsPedFullMap_->Fill(EEid.ix(), EEid.iy(), uncalibRecHit->pedestal());
  
     // general checks, with threshold at 60 ADC counts
     if (uncalibRecHit->amplitude() > 60) {
       if (mySide > 0) {
         if (meEEUncalibRecHitsOccupancyPlusGt60adc_)
           meEEUncalibRecHitsOccupancyPlusGt60adc_->Fill(EEid.ix(), EEid.iy());
       }
       if (mySide < 0) {
         if (meEEUncalibRecHitsOccupancyMinusGt60adc_)
           meEEUncalibRecHitsOccupancyMinusGt60adc_->Fill(EEid.ix(), EEid.iy());
       }
       if (meEEUncalibRecHitsAmplitudeGt60adc_)
         meEEUncalibRecHitsAmplitudeGt60adc_->Fill(uncalibRecHit->amplitude());
       if (meEEUncalibRecHitsPedestalGt60adc_)
         meEEUncalibRecHitsPedestalGt60adc_->Fill(uncalibRecHit->pedestal());
       if (meEEUncalibRecHitsJitterGt60adc_)
         meEEUncalibRecHitsJitterGt60adc_->Fill(uncalibRecHit->jitter());
       if (meEEUncalibRecHitsChi2Gt60adc_)
         meEEUncalibRecHitsChi2Gt60adc_->Fill(uncalibRecHit->chi2());
     }
  
     if (!skipDigis) {
       // Find the rechit corresponding digi
       EEDigiCollection::const_iterator myDigi = EEDigi->find(EEid);
       // int sMax    = -1; // UNUSED
       double eMax = 0.;
       if (myDigi != EEDigi->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 = &c.getData(ecalPeds);
       EcalPedestalsMap::const_iterator it = myped->getMap().find(EEid);
       if (it != myped->getMap().end()) {
         if (eMax > (*it).mean_x1 + 5 * (*it).rms_x1 && eMax != 0) {  // only real signal RecHit
  
           if (meEEUncalibRecHitMaxSampleRatio_) {
             meEEUncalibRecHitMaxSampleRatio_->Fill((uncalibRecHit->amplitude() + uncalibRecHit->pedestal()) / eMax);
           }
  
           if (meEEUncalibRecHitMaxSampleRatioGt60adc_ && (uncalibRecHit->amplitude() > 60)) {
             meEEUncalibRecHitMaxSampleRatioGt60adc_->Fill((uncalibRecHit->amplitude() + uncalibRecHit->pedestal()) /
                                                           eMax);
           }
  
           LogDebug("EcalRecHitsTaskInfo")
               << "endcap, eMax = " << eMax << " Amplitude = " << uncalibRecHit->amplitude() + uncalibRecHit->pedestal();
         } else
           continue;
       } else
         continue;
     }
   }  // loop over the UncalibratedRecHitCollection
 }

References EcalMGPASample::adc(), c, MillePedeFileConverter_cfg::e, ecalPeds, CollectionTags_cfi::EEDigi, EEdigiCollection_token_, CollectionTags_cfi::EEUncalibRecHit, EEuncalibrechitCollection_token_, cosmicPhotonAnalyzer_cfi::eMax, EcalCondObjectContainer< T >::end(), dqm::impl::MonitorElement::Fill(), EcalCondObjectContainer< T >::find(), EcalCondObjectContainer< T >::getMap(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), LogDebug, meEEUncalibRecHitMaxSampleRatio_, meEEUncalibRecHitMaxSampleRatioGt60adc_, meEEUncalibRecHitsAmpFullMap_, meEEUncalibRecHitsAmplitude_, meEEUncalibRecHitsAmplitudeGt60adc_, meEEUncalibRecHitsChi2_, meEEUncalibRecHitsChi2Gt60adc_, meEEUncalibRecHitsJitter_, meEEUncalibRecHitsJitterGt60adc_, meEEUncalibRecHitsOccupancyMinus_, meEEUncalibRecHitsOccupancyMinusGt60adc_, meEEUncalibRecHitsOccupancyPlus_, meEEUncalibRecHitsOccupancyPlusGt60adc_, meEEUncalibRecHitsPedestal_, meEEUncalibRecHitsPedestalGt60adc_, meEEUncalibRecHitsPedFullMap_, edm::Handle< T >::product(), simplePhotonAnalyzer_cfi::sample, and EEDetId::zside().

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 47 of file EcalEndcapRecHitsValidation.cc.

                                                                         {
   Char_t histo[200];
  
   ibooker.setCurrentFolder("EcalRecHitsV/EcalEndcapRecHitsTask");
  
   sprintf(histo, "EE+ Occupancy");
   meEEUncalibRecHitsOccupancyPlus_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
  
   sprintf(histo, "EE- Occupancy");
   meEEUncalibRecHitsOccupancyMinus_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
  
   sprintf(histo, "EE Amplitude");
   meEEUncalibRecHitsAmplitude_ = ibooker.book1D(histo, histo, 201, -20., 4000.);
  
   sprintf(histo, "EE Pedestal");
   meEEUncalibRecHitsPedestal_ = ibooker.book1D(histo, histo, 50, 190., 210.);
  
   sprintf(histo, "EE Jitter");
   meEEUncalibRecHitsJitter_ = ibooker.book1D(histo, histo, 100, 0., 100.);
  
   sprintf(histo, "EE Chi2");
   meEEUncalibRecHitsChi2_ = ibooker.book1D(histo, histo, 100, 18000., 22000.);
  
   sprintf(histo, "EE RecHit Max Sample Ratio");
   meEEUncalibRecHitMaxSampleRatio_ = ibooker.book1D(histo, histo, 120, 0.90, 1.05);
  
   sprintf(histo, "EE+ Occupancy gt 60 adc counts");
   meEEUncalibRecHitsOccupancyPlusGt60adc_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
  
   sprintf(histo, "EE- Occupancy gt 60 adc counts");
   meEEUncalibRecHitsOccupancyMinusGt60adc_ = ibooker.book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
  
   sprintf(histo, "EE Amplitude gt 60 adc counts");
   meEEUncalibRecHitsAmplitudeGt60adc_ = ibooker.book1D(histo, histo, 200, 0., 4000.);
  
   sprintf(histo, "EE Pedestal gt 60 adc counts");
   meEEUncalibRecHitsPedestalGt60adc_ = ibooker.book1D(histo, histo, 50, 190., 210.);
  
   sprintf(histo, "EE Jitter gt 60 adc counts");
   meEEUncalibRecHitsJitterGt60adc_ = ibooker.book1D(histo, histo, 100, 0., 100.);
  
   sprintf(histo, "EE Chi2 gt 60 adc counts");
   meEEUncalibRecHitsChi2Gt60adc_ = ibooker.book1D(histo, histo, 100, 18000., 22000.);
  
   sprintf(histo, "EE RecHit Max Sample Ratio gt 60 adc counts");
   meEEUncalibRecHitMaxSampleRatioGt60adc_ = ibooker.book1D(histo, histo, 120, 0.90, 1.05);
  
   sprintf(histo, "EE Amplitude Full Map");
   meEEUncalibRecHitsAmpFullMap_ = ibooker.bookProfile2D(histo, histo, 100, 0., 100., 100, 0., 100., 200, 0., 4000.);
  
   sprintf(histo, "EE Pedestal Full Map");
   meEEUncalibRecHitsPedFullMap_ = ibooker.bookProfile2D(histo, histo, 100, 0., 100., 100, 0., 100., 50, 194., 201.);
 }