#include <EcalTBValidation.h>

Inheritance diagram for EcalTBValidation:

Public Member Functions
virtual 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 ()

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

Detailed Description

Definition at line 21 of file EcalTBValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file EcalTBValidation.cc.

References data_, digi_Token_, digiCollection_, digiProducer_, eventHeader_Token_, eventHeaderCollection_, eventHeaderProducer_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hit_Token_, hitCollection_, hitProducer_, hodoRec_Token_, hodoRecInfoCollection_, hodoRecInfoProducer_, 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_ = 0;
   meETBampltdc_ = 0;
   meETBShape_ = 0;
   meETBhodoX_ = 0;
   meETBhodoY_ = 0;
   meETBe1x1_ = 0;
   meETBe3x3_ = 0;
   meETBe5x5_ = 0;
   meETBe1e9_ = 0;
   meETBe1e25_ = 0;
   meETBe9e25_ = 0;
   meETBe1x1_center_ = 0;
   meETBe3x3_center_ = 0;
   meETBe5x5_center_ = 0;
   meETBe1vsX_ = 0;
   meETBe1vsY_ = 0;
   meETBe1e9vsX_ = 0;
   meETBe1e9vsY_ = 0;
   meETBe1e25vsX_ = 0;
   meETBe1e25vsY_ = 0;
   meETBe9e25vsX_ = 0;
   meETBe9e25vsY_ = 0;
 
 }

EcalTBValidation::~EcalTBValidation ( )

Definition at line 69 of file EcalTBValidation.cc.

69 {}

Member Function Documentation

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

overridevirtual

Definition at line 123 of file EcalTBValidation.cc.

References EcalMGPASample::adc(), CustomPhysics_cfi::amplitude, MessageLogger_cfi::cerr, data_, digi_Token_, digiCollection_, photonPostprocessing_cfi::eMax, edm::DataFrameContainer::end(), EBDetId::ETAPHIMODE, eventHeader_Token_, eventHeaderProducer_, MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::SortedCollection< T, SORT >::find(), edm::DataFrameContainer::find(), hit_Token_, hitCollection_, hodoRec_Token_, hodoRecInfoCollection_, EBDetId::ieta(), cuy::ii, 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(), EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), edm::Handle< T >::product(), simplePhotonAnalyzer_cfi::sample, EcalDataFrame::sample(), EcalDataFrame::size(), EBDetId::SMCRYSTALMODE, EcalTBEventHeader::tableIsMoving(), tdcRec_Token_, tdcRecInfoCollection_, EBDetId::validDetId(), and xtalInBeam_.

                                                                                   {
 
   using namespace edm;
   using namespace cms;
 
   // digis
   const EBDigiCollection* theDigis=0;
   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=0;  
   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=0;  
   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=0;
   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=0;
   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);
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 71 of file EcalTBValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::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_, DQMStore::IBooker::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 );
 
 }