#include <EcalTBValidation.h>

Inheritance diagram for EcalTBValidation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	EcalTBValidation (const edm::ParameterSet &)

	~EcalTBValidation () override

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

Private Attributes
int	data_

edm::EDGetTokenT< EBDigiCollection >	digi_Token_

std::string	digiCollection_

std::string	digiProducer_

edm::EDGetTokenT< EcalTBEventHeader >	eventHeader_Token_

std::string	eventHeaderCollection_

std::string	eventHeaderProducer_

edm::EDGetTokenT< EBUncalibratedRecHitCollection >	hit_Token_

std::string	hitCollection_

std::string	hitProducer_

edm::EDGetTokenT< EcalTBHodoscopeRecInfo >	hodoRec_Token_

std::string	hodoRecInfoCollection_

std::string	hodoRecInfoProducer_

MonitorElement *	meETBampltdc_

MonitorElement *	meETBe1e25_

MonitorElement *	meETBe1e25vsX_

MonitorElement *	meETBe1e25vsY_

MonitorElement *	meETBe1e9_

MonitorElement *	meETBe1e9vsX_

MonitorElement *	meETBe1e9vsY_

MonitorElement *	meETBe1vsX_

MonitorElement *	meETBe1vsY_

MonitorElement *	meETBe1x1_

MonitorElement *	meETBe1x1_center_

MonitorElement *	meETBe3x3_

MonitorElement *	meETBe3x3_center_

MonitorElement *	meETBe5x5_

MonitorElement *	meETBe5x5_center_

MonitorElement *	meETBe9e25_

MonitorElement *	meETBe9e25vsX_

MonitorElement *	meETBe9e25vsY_

MonitorElement *	meETBhodoX_

MonitorElement *	meETBhodoY_

MonitorElement *	meETBShape_

MonitorElement *	meETBxib_

edm::EDGetTokenT< EcalTBTDCRecInfo >	tdcRec_Token_

std::string	tdcRecInfoCollection_

std::string	tdcRecInfoProducer_

bool	verbose_

int	xtalInBeam_

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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 19 of file EcalTBValidation.h.

Constructor & Destructor Documentation

◆ EcalTBValidation()

EcalTBValidation::EcalTBValidation ( const edm::ParameterSet & config )

explicit

Definition at line 15 of file EcalTBValidation.cc.

References data_, digi_Token_, digiCollection_, digiProducer_, eventHeader_Token_, eventHeaderCollection_, eventHeaderProducer_, hit_Token_, hitCollection_, hitProducer_, hodoRec_Token_, hodoRecInfoCollection_, hodoRecInfoProducer_, ProducerED_cfi::InputTag, meETBampltdc_, meETBe1e25_, meETBe1e25vsX_, meETBe1e25vsY_, meETBe1e9_, meETBe1e9vsX_, meETBe1e9vsY_, meETBe1vsX_, meETBe1vsY_, meETBe1x1_, meETBe1x1_center_, meETBe3x3_, meETBe3x3_center_, meETBe5x5_, meETBe5x5_center_, meETBe9e25_, meETBe9e25vsX_, meETBe9e25vsY_, meETBhodoX_, meETBhodoY_, meETBShape_, meETBxib_, AlCaHLTBitMon_QueryRunRegistry::string, tdcRec_Token_, tdcRecInfoCollection_, tdcRecInfoProducer_, verbose_, and xtalInBeam_.

                                                                 {
   data_ = config.getUntrackedParameter<int>("data", -1000);
   xtalInBeam_ = config.getUntrackedParameter<int>("xtalInBeam", -1000);
 
   digiCollection_ = config.getParameter<std::string>("digiCollection");
   digiProducer_ = config.getParameter<std::string>("digiProducer");
   hitCollection_ = config.getParameter<std::string>("hitCollection");
   hitProducer_ = config.getParameter<std::string>("hitProducer");
   hodoRecInfoCollection_ = config.getParameter<std::string>("hodoRecInfoCollection");
   hodoRecInfoProducer_ = config.getParameter<std::string>("hodoRecInfoProducer");
   tdcRecInfoCollection_ = config.getParameter<std::string>("tdcRecInfoCollection");
   tdcRecInfoProducer_ = config.getParameter<std::string>("tdcRecInfoProducer");
   eventHeaderCollection_ = config.getParameter<std::string>("eventHeaderCollection");
   eventHeaderProducer_ = config.getParameter<std::string>("eventHeaderProducer");
 
   digi_Token_ = consumes<EBDigiCollection>(edm::InputTag(digiProducer_, digiCollection_));
   hit_Token_ = consumes<EBUncalibratedRecHitCollection>(edm::InputTag(hitProducer_, hitCollection_));
   hodoRec_Token_ = consumes<EcalTBHodoscopeRecInfo>(edm::InputTag(hodoRecInfoProducer_, hodoRecInfoCollection_));
   tdcRec_Token_ = consumes<EcalTBTDCRecInfo>(edm::InputTag(tdcRecInfoProducer_, tdcRecInfoCollection_));
   eventHeader_Token_ = consumes<EcalTBEventHeader>(edm::InputTag(eventHeaderProducer_));
   // rootfile_              =
   // config.getUntrackedParameter<std::string>("rootfile","EcalTBValidation.root");
 
   // verbosity...
   verbose_ = config.getUntrackedParameter<bool>("verbose", false);
 
   meETBxib_ = nullptr;
   meETBampltdc_ = nullptr;
   meETBShape_ = nullptr;
   meETBhodoX_ = nullptr;
   meETBhodoY_ = nullptr;
   meETBe1x1_ = nullptr;
   meETBe3x3_ = nullptr;
   meETBe5x5_ = nullptr;
   meETBe1e9_ = nullptr;
   meETBe1e25_ = nullptr;
   meETBe9e25_ = nullptr;
   meETBe1x1_center_ = nullptr;
   meETBe3x3_center_ = nullptr;
   meETBe5x5_center_ = nullptr;
   meETBe1vsX_ = nullptr;
   meETBe1vsY_ = nullptr;
   meETBe1e9vsX_ = nullptr;
   meETBe1e9vsY_ = nullptr;
   meETBe1e25vsX_ = nullptr;
   meETBe1e25vsY_ = nullptr;
   meETBe9e25vsX_ = nullptr;
   meETBe9e25vsY_ = nullptr;
 }

◆ ~EcalTBValidation()

EcalTBValidation::~EcalTBValidation ( )

override

Definition at line 65 of file EcalTBValidation.cc.

65 {}

Member Function Documentation

◆ analyze()

void EcalTBValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 117 of file EcalTBValidation.cc.

References EcalMGPASample::adc(), CustomPhysics_cfi::amplitude, DMR_cfg::cerr, data_, digi_Token_, digiCollection_, cosmicPhotonAnalyzer_cfi::eMax, edm::DataFrameContainer::end(), EBDetId::ETAPHIMODE, eventHeader_Token_, eventHeaderProducer_, dqm::impl::MonitorElement::Fill(), edm::DataFrameContainer::find(), hit_Token_, hitCollection_, hodoRec_Token_, hodoRecInfoCollection_, hcalRecHitTable_cff::ieta, EBDetId::ieta(), cuy::ii, hcalRecHitTable_cff::iphi, EBDetId::iphi(), EBDetId::ism(), edm::HandleBase::isValid(), meETBampltdc_, meETBe1e25_, meETBe1e25vsX_, meETBe1e25vsY_, meETBe1e9_, meETBe1e9vsX_, meETBe1e9vsY_, meETBe1vsX_, meETBe1vsY_, meETBe1x1_, meETBe1x1_center_, meETBe3x3_, meETBe3x3_center_, meETBe5x5_, meETBe5x5_center_, meETBe9e25_, meETBe9e25vsX_, meETBe9e25vsY_, meETBhodoX_, meETBhodoY_, meETBShape_, meETBxib_, EcalTBTDCRecInfo::offset(), phits, EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), edm::Handle< T >::product(), EcalDataFrame::sample(), ecalGpuTask_cfi::sample, EcalDataFrame::size(), EBDetId::SMCRYSTALMODE, EcalTBEventHeader::tableIsMoving(), tdcRec_Token_, tdcRecInfoCollection_, bscTrigger_cfi::theHits, EBDetId::validDetId(), and xtalInBeam_.

                                                                                 {
   using namespace edm;
   using namespace cms;
 
   // digis
   const EBDigiCollection *theDigis = nullptr;
   Handle<EBDigiCollection> pdigis;
   event.getByToken(digi_Token_, pdigis);
   if (pdigis.isValid()) {
     theDigis = pdigis.product();
   } else {
     std::cerr << "Error! can't get the product " << digiCollection_.c_str() << std::endl;
     return;
   }
 
   // rechits
   const EBUncalibratedRecHitCollection *theHits = nullptr;
   Handle<EBUncalibratedRecHitCollection> phits;
   event.getByToken(hit_Token_, phits);
   if (phits.isValid()) {
     theHits = phits.product();
   } else {
     std::cerr << "Error! can't get the product " << hitCollection_.c_str() << std::endl;
     return;
   }
 
   // hodoscopes
   const EcalTBHodoscopeRecInfo *theHodo = nullptr;
   Handle<EcalTBHodoscopeRecInfo> pHodo;
   event.getByToken(hodoRec_Token_, pHodo);
   if (pHodo.isValid()) {
     theHodo = pHodo.product();
   } else {
     std::cerr << "Error! can't get the product " << hodoRecInfoCollection_.c_str() << std::endl;
     return;
   }
 
   // tdc
   const EcalTBTDCRecInfo *theTDC = nullptr;
   Handle<EcalTBTDCRecInfo> pTDC;
   event.getByToken(tdcRec_Token_, pTDC);
   if (pTDC.isValid()) {
     theTDC = pTDC.product();
   } else {
     std::cerr << "Error! can't get the product " << tdcRecInfoCollection_.c_str() << std::endl;
     return;
   }
 
   // event header
   const EcalTBEventHeader *evtHeader = nullptr;
   Handle<EcalTBEventHeader> pEventHeader;
   event.getByToken(eventHeader_Token_, pEventHeader);
   if (pEventHeader.isValid()) {
     evtHeader = pEventHeader.product();
   } else {
     std::cerr << "Error! can't get the product " << eventHeaderProducer_.c_str() << std::endl;
     return;
   }
 
   // -----------------------------------------------------------------------
   // xtal-in-beam
   EBDetId xtalInBeamId(1, xtalInBeam_, EBDetId::SMCRYSTALMODE);
   if (xtalInBeamId == EBDetId(0)) {
     return;
   }
   int xibEta = xtalInBeamId.ieta();
   int xibPhi = xtalInBeamId.iphi();
 
   // skipping events with moving table (in data)
   if (data_ && (evtHeader->tableIsMoving()))
     return;
 
   // amplitudes
   EBDetId Xtals5x5[25];
   for (unsigned int icry = 0; icry < 25; icry++) {
     unsigned int row = icry / 5;
     unsigned int column = icry % 5;
     int ieta = xtalInBeamId.ieta() + column - 2;
     int iphi = xtalInBeamId.iphi() + row - 2;
     if (EBDetId::validDetId(ieta, iphi)) {
       EBDetId tempId(ieta, iphi, EBDetId::ETAPHIMODE);
       if (tempId.ism() == 1)
         Xtals5x5[icry] = tempId;
       else
         Xtals5x5[icry] = EBDetId(0);
     } else {
       Xtals5x5[icry] = EBDetId(0);
     }
   }
 
   // matrices
   double ampl1x1 = 0.;
   double ampl3x3 = 0.;
   double ampl5x5 = 0.;
   for (unsigned int icry = 0; icry < 25; icry++) {
     if (!Xtals5x5[icry].null()) {
       double theAmpl = (theHits->find(Xtals5x5[icry]))->amplitude();
       ampl5x5 += theAmpl;
       if (icry == 12) {
         ampl1x1 = theAmpl;
       }
       if (icry == 6 || icry == 7 || icry == 8 || icry == 11 || icry == 12 || icry == 13 || icry == 16 || icry == 17 ||
           icry == 18) {
         ampl3x3 += theAmpl;
       }
     }
   }
 
   // pulse shape
   double sampleSave[10];
   for (int ii = 0; ii < 10; ++ii) {
     sampleSave[ii] = 0.0;
   }
   EBDigiCollection::const_iterator thisDigi = theDigis->find(xtalInBeamId);
   // int sMax = -1; // UNUSED
   double eMax = 0.;
   if (thisDigi != theDigis->end()) {
     EBDataFrame myDigi = (*thisDigi);
     for (int sample = 0; sample < myDigi.size(); ++sample) {
       double analogSample = myDigi.sample(sample).adc();
       sampleSave[sample] = analogSample;
       if (eMax < analogSample) {
         eMax = analogSample;
         // sMax = sample; // UNUSED
       }
     }
   }
 
   // beam profile
   double xBeam = theHodo->posX();
   double yBeam = theHodo->posY();
 
   // filling histos
   meETBxib_->Fill(xibEta, xibPhi);
   meETBhodoX_->Fill(xBeam);
   meETBhodoY_->Fill(yBeam);
   meETBampltdc_->Fill(theTDC->offset(), ampl1x1);
   meETBe1x1_->Fill(ampl1x1);
   meETBe3x3_->Fill(ampl3x3);
   meETBe5x5_->Fill(ampl5x5);
   meETBe1e9_->Fill(ampl1x1 / ampl3x3);
   meETBe1e25_->Fill(ampl1x1 / ampl5x5);
   meETBe9e25_->Fill(ampl3x3 / ampl5x5);
   meETBe1vsX_->Fill(xBeam, ampl1x1);
   meETBe1vsY_->Fill(yBeam, ampl1x1);
   meETBe1e9vsX_->Fill(xBeam, ampl1x1 / ampl3x3);
   meETBe1e9vsY_->Fill(yBeam, ampl1x1 / ampl3x3);
   meETBe1e25vsX_->Fill(xBeam, ampl1x1 / ampl5x5);
   meETBe1e25vsY_->Fill(yBeam, ampl1x1 / ampl5x5);
   meETBe9e25vsX_->Fill(xBeam, ampl3x3 / ampl5x5);
   meETBe9e25vsY_->Fill(yBeam, ampl3x3 / ampl5x5);
 
   for (int ii = 0; ii < 10; ++ii) {
     meETBShape_->Fill(double(ii) + theTDC->offset(), sampleSave[ii]);
   }
 
   if ((fabs(xBeam) < 2.5) && (fabs(yBeam) < 2.5)) {
     meETBe1x1_center_->Fill(ampl1x1);
     meETBe3x3_center_->Fill(ampl3x3);
     meETBe5x5_center_->Fill(ampl5x5);
   }
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 67 of file EcalTBValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), meETBampltdc_, meETBe1e25_, meETBe1e25vsX_, meETBe1e25vsY_, meETBe1e9_, meETBe1e9vsX_, meETBe1e9vsY_, meETBe1vsX_, meETBe1vsY_, meETBe1x1_, meETBe1x1_center_, meETBe3x3_, meETBe3x3_center_, meETBe5x5_, meETBe5x5_center_, meETBe9e25_, meETBe9e25vsX_, meETBe9e25vsY_, meETBhodoX_, meETBhodoY_, meETBShape_, meETBxib_, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                      {
   std::string hname;
   ibooker.setCurrentFolder("EcalRecHitsV/EcalTBValidationTask");
 
   hname = "xtal in beam position";
   meETBxib_ = ibooker.book2D(hname, hname, 85, 0., 85., 20, 0., 20.);
   hname = "Max Amplitude vs TDC offset";
   meETBampltdc_ = ibooker.book2D(hname, hname, 100, 0., 1., 1000, 0., 4000.);
   hname = "Beam Profile X";
   meETBhodoX_ = ibooker.book1D(hname, hname, 100, -20., 20.);
   hname = "Beam Profile Y";
   meETBhodoY_ = ibooker.book1D(hname, hname, 100, -20., 20.);
   hname = "E1x1 energy";
   meETBe1x1_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E3x3 energy";
   meETBe3x3_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E5x5 energy";
   meETBe5x5_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E1x1 energy center";
   meETBe1x1_center_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E3x3 energy center";
   meETBe3x3_center_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E5x5 energy center";
   meETBe5x5_center_ = ibooker.book1D(hname, hname, 1000, 0., 4000.);
   hname = "E1 over E9 ratio";
   meETBe1e9_ = ibooker.book1D(hname, hname, 600, 0., 1.2);
   hname = "E1 over E25 ratio";
   meETBe1e25_ = ibooker.book1D(hname, hname, 600, 0., 1.2);
   hname = "E9 over E25 ratio";
   meETBe9e25_ = ibooker.book1D(hname, hname, 600, 0., 1.2);
   hname = "E1 vs X";
   meETBe1vsX_ = ibooker.book2D(hname, hname, 80, -20, 20, 1000, 0., 4000.);
   hname = "E1 vs Y";
   meETBe1vsY_ = ibooker.book2D(hname, hname, 80, -20, 20, 1000, 0., 4000.);
   hname = "E1 over E9 vs X";
   meETBe1e9vsX_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "E1 over E9 vs Y";
   meETBe1e9vsY_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "E1 over E25 vs X";
   meETBe1e25vsX_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "E1 over E25 vs Y";
   meETBe1e25vsY_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "E9 over E25 vs X";
   meETBe9e25vsX_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "E9 over E25 vs Y";
   meETBe9e25vsY_ = ibooker.book2D(hname, hname, 80, -20, 20, 600, 0., 1.2);
   hname = "Xtal in Beam Shape";
   meETBShape_ = ibooker.book2D(hname, hname, 250, 0, 10, 350, 0, 3500);
 }