#include <Validation/HGCalValidation/plugins/HGCalDigiValidation.cc>

Inheritance diagram for HGCalDigiValidation:

Classes
struct	digiInfo

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

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

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

	HGCalDigiValidation (const edm::ParameterSet &)

	~HGCalDigiValidation ()

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

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

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

virtual 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

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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 &	itemsToGetFromEvent () 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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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)

Private Member Functions
template<class T1 , class T2 >
void	digiValidation (const T1 &detId, const T2 *geom, int, uint16_t, double)

void	fillDigiInfo (digiInfo &hinfo)

void	fillDigiInfo ()

void	fillOccupancyMap (std::map< int, int > &OccupancyMap, int layer)

Private Attributes
std::vector< MonitorElement * >	ADC_

std::vector< MonitorElement * >	charge_

std::vector< MonitorElement * >	DigiOccupancy_Minus_

std::vector< MonitorElement * >	DigiOccupancy_Plus_

std::vector< MonitorElement * >	DigiOccupancy_XY_

edm::EDGetToken	digiSource_

bool	ifHCAL_

int	layers_

MonitorElement *	MeanDigiOccupancy_Minus_

MonitorElement *	MeanDigiOccupancy_Plus_

std::string	nameDetector_

std::map< int, int >	OccupancyMap_minus_

std::map< int, int >	OccupancyMap_plus_

int	SampleIndx_

int	verbosity_

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

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

Description: Validates SimHits of High Granularity Calorimeter Implementation: [Notes on implementation]

Definition at line 43 of file HGCalDigiValidation.h.

Constructor & Destructor Documentation

HGCalDigiValidation::HGCalDigiValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 51 of file HGCalDigiValidation.cc.

References digiSource_, Exception, edm::ParameterSet::getParameter(), ifHCAL_, nameDetector_, and groupFilesInBlocks::temp.

                                                                        :
   nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
   verbosity_(iConfig.getUntrackedParameter<int>("Verbosity",0)),
   SampleIndx_(iConfig.getUntrackedParameter<int>("SampleIndx",5)) {
   ifHCAL_   = iConfig.getParameter<bool>("ifHCAL");
   auto temp = iConfig.getParameter<edm::InputTag>("DigiSource");
   if( nameDetector_ == "HGCalEESensitive" ) {
     digiSource_    = consumes<HGCEEDigiCollection>(temp);
   } else if ( nameDetector_ == "HGCalHESiliconSensitive" ||
               nameDetector_ == "HGCalHEScintillatorSensitive" ) {
     digiSource_    = consumes<HGCHEDigiCollection>(temp);
   } else if ( nameDetector_ == "HCal" ) {
     if (ifHCAL_) digiSource_ = consumes<QIE11DigiCollection>(temp);
     else         digiSource_ = consumes<HGCBHDigiCollection>(temp);
   } else {
     throw cms::Exception("BadHGCDigiSource")
       << "HGCal DetectorName given as " << nameDetector_ << " must be: "
       << "\"HGCalHESiliconSensitive\", \"HGCalHESiliconSensitive\", "
       << "\"HGCalHEScintillatorSensitive\", or \"HCal\"!"; 
   }  
 }

HGCalDigiValidation::~HGCalDigiValidation ( )

Definition at line 74 of file HGCalDigiValidation.cc.

74 { }

Member Function Documentation

void HGCalDigiValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 76 of file HGCalDigiValidation.cc.

References ecalMGPA::adc(), HcalCoderDb::adc2fC(), RecoTauCleanerPlugins::charge, HGCSample::data(), HcalDetId::depth(), digiSource_, digiValidation(), edm::EventID::event(), fillDigiInfo(), relativeConstraints::geom, motionPolicyChamber::geom0, geometryDiff::geom1, edm::EventSetup::get(), edm::Event::getByToken(), HcalEndcap, QIE11DataFrame::id(), edm::EventBase::id(), ifHCAL_, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), HGCEEDetId::layer(), HGCHEDetId::layer(), nameDetector_, OccupancyMap_minus_, OccupancyMap_plus_, HcalCalibrations::pedestal(), edm::ESHandle< class >::product(), SampleIndx_, HcalDetId::subdet(), HGCSample::toa(), and verbosity_.

                                                { 
   OccupancyMap_plus_.clear();
   OccupancyMap_minus_.clear();
   
   const HGCalGeometry* geom0(0);
   const CaloGeometry*  geom1(0);
   if (nameDetector_ == "HCal") {
     edm::ESHandle<CaloGeometry> geom;
     iSetup.get<CaloGeometryRecord>().get(geom);
     if (!geom.isValid()) edm::LogWarning("HGCalValidation") << "Cannot get valid HGCalGeometry Object for " << nameDetector_;
     geom1 = geom.product();
   } else {
     edm::ESHandle<HGCalGeometry> geom;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, geom);
     if (!geom.isValid()) edm::LogWarning("HGCalValidation") << "Cannot get valid HGCalGeometry Object for " << nameDetector_;
     geom0 = geom.product();
   }
 
   unsigned int ntot(0), nused(0);
   if (nameDetector_ == "HGCalEESensitive") {
     //HGCalEE
     edm::Handle<HGCEEDigiCollection> theHGCEEDigiContainers;
     iEvent.getByToken(digiSource_, theHGCEEDigiContainers);
     if (theHGCEEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                     << theHGCEEDigiContainers->size() 
                     << " element(s)";
       
       for (HGCEEDigiCollection::const_iterator it =theHGCEEDigiContainers->begin();
        it !=theHGCEEDigiContainers->end(); ++it) {
     ntot++; nused++;
     HGCEEDetId detId     = (it->id());
     int        layer     = detId.layer();
     HGCSample  hgcSample = it->sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom0, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "HGCEEDigiCollection handle does not exist !!!";
     }
   } else if ((nameDetector_ == "HGCalHESiliconSensitive") || 
          (nameDetector_ == "HGCalHEScintillatorSensitive")) {
     //HGCalHE
     edm::Handle<HGCHEDigiCollection> theHGCHEDigiContainers;
     iEvent.getByToken(digiSource_, theHGCHEDigiContainers);
     if (theHGCHEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                     << theHGCHEDigiContainers->size()
                     << " element(s)";
       
       for (HGCHEDigiCollection::const_iterator it =theHGCHEDigiContainers->begin();
        it !=theHGCHEDigiContainers->end(); ++it) {
     ntot++; nused++;
     HGCHEDetId detId     = (it->id());
     int        layer     = detId.layer();
     HGCSample  hgcSample = it->sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom0, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "HGCHEDigiCollection handle does not exist !!!";
     }
   } else if ((nameDetector_ == "HCal") && (!ifHCAL_)) {
     //HGCalBH
     edm::Handle<HGCBHDigiCollection> theHGCBHDigiContainers;
     iEvent.getByToken(digiSource_, theHGCBHDigiContainers);
     if (theHGCBHDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                     << theHGCBHDigiContainers->size()
                     << " element(s)";
       
       for (HGCBHDigiCollection::const_iterator it =theHGCBHDigiContainers->begin();
        it !=theHGCBHDigiContainers->end(); ++it) {
     ntot++; nused++;
     HcalDetId  detId     = (it->id());
     int        layer     = detId.depth();
     HGCSample  hgcSample = it->sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom1, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "HGCBHDigiCollection handle does not exist !!!";
     }
   } else if (nameDetector_ == "HCal") {
     //HE
     edm::Handle<QIE11DigiCollection>  theHEDigiContainers;
     iEvent.getByToken(digiSource_, theHEDigiContainers);
     if (theHEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                     << theHEDigiContainers->size()
                     << " element(s)";
       edm::ESHandle<HcalDbService> conditions;
       iSetup.get<HcalDbRecord > ().get(conditions);
 
       for (QIE11DigiCollection::const_iterator it =theHEDigiContainers->begin();
        it !=theHEDigiContainers->end(); ++it) {
     QIE11DataFrame df(*it);
     HcalDetId detId  = (df.id());
     ntot++;
     if (detId.subdet() == HcalEndcap) {
       nused++;
       HcalCalibrations calibrations = conditions->getHcalCalibrations(detId);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder(detId);
       const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
       HcalCoderDb coder(*channelCoder, *shape);
       CaloSamples tool;
       coder.adc2fC(df, tool);
       int       layer  = detId.depth();
       uint16_t  adc    = (df)[SampleIndx_].adc();
       int       capid  = (df)[SampleIndx_].capid();
       double    charge = (tool[SampleIndx_] - calibrations.pedestal(capid));
       digiValidation(detId, geom1, layer, adc, charge);
     }
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "HGCHEDigiCollection handle does not exist !!!";
     }
   } else {
     edm::LogWarning("HGCalValidation") << "invalid detector name !! " 
                        << nameDetector_;
   }
   edm::LogInfo("HGCalValidation") << "Event " << iEvent.id().event()
                   << " with " << ntot << " total and "
                   << nused << " used digis";
 }

void HGCalDigiValidation::bookHistograms	(	DQMStore::IBooker &	iB,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 296 of file HGCalDigiValidation.cc.

References ADC_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), charge_, DigiOccupancy_Minus_, DigiOccupancy_Plus_, DigiOccupancy_XY_, layers_, MeanDigiOccupancy_Minus_, MeanDigiOccupancy_Plus_, nameDetector_, and DQMStore::IBooker::setCurrentFolder().

                                            {
   
   iB.setCurrentFolder("HGCalDigiV/"+nameDetector_);
 
   std::ostringstream histoname;
   for (int ilayer = 0; ilayer < layers_; ilayer++ ) {
     histoname.str(""); histoname << "charge_"<< "layer_" << ilayer;
     charge_.push_back(iB.book1D(histoname.str().c_str(),"charge_",100,-25,25));
       
     histoname.str(""); histoname << "ADC_" << "layer_" << ilayer;
     ADC_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy",200,0,1000));
       
     histoname.str(""); histoname << "DigiOccupancy_XY_" << "layer_" << ilayer;
     DigiOccupancy_XY_.push_back(iB.book2D(histoname.str().c_str(), "DigiOccupancy", 50, -500, 500, 50, -500, 500));
       
     histoname.str(""); histoname << "DigiOccupancy_Plus_" << "layer_" << ilayer;
     DigiOccupancy_Plus_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy +z", 100, 0, 1000));
     histoname.str(""); histoname << "DigiOccupancy_Minus_" << "layer_" << ilayer;
     DigiOccupancy_Minus_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy -z", 100, 0, 1000));
   }
 
   histoname.str(""); histoname << "SUMOfDigiOccupancy_Plus";
   MeanDigiOccupancy_Plus_ = iB.book1D( histoname.str().c_str(), "SUMOfDigiOccupancy_Plus", layers_, -0.5, layers_-0.5);
   histoname.str(""); histoname << "SUMOfRecDigiOccupancy_Minus";
   MeanDigiOccupancy_Minus_ = iB.book1D( histoname.str().c_str(), "SUMOfDigiOccupancy_Minus", layers_, -0.5,layers_-0.5);
 }

template<class T1 , class T2 >

void HGCalDigiValidation::digiValidation	(	const T1 &	detId,
		const T2 *	geom,
		int	layer,
		uint16_t	adc,
		double	charge
	)

private

Definition at line 218 of file HGCalDigiValidation.cc.

References ecalMGPA::adc(), HGCalDigiValidation::digiInfo::adc, RecoTauCleanerPlugins::charge, HGCalDigiValidation::digiInfo::charge, fillDigiInfo(), fillOccupancyMap(), HGCalDigiValidation::digiInfo::layer, OccupancyMap_minus_, OccupancyMap_plus_, verbosity_, HGCalDigiValidation::digiInfo::x, HGCalDigiValidation::digiInfo::y, and HGCalDigiValidation::digiInfo::z.

Referenced by analyze().

                                                              {
   
   if (verbosity_>1) edm::LogInfo("HGCalValidation") << detId;
   DetId id1 = DetId(detId.rawId());
   GlobalPoint global1 = geom->getPosition(id1);
   
   if (verbosity_>1) 
     edm::LogInfo("HGCalValidation") << " adc = "         <<  adc
                     << " charge = "      <<  charge;
   
   digiInfo   hinfo; 
   hinfo.x       =  global1.x();
   hinfo.y       =  global1.y();
   hinfo.z       =  global1.z();
   hinfo.adc     =  adc;
   hinfo.charge  =  charge; //charges[0];
   hinfo.layer   =  layer;
   
   if (verbosity_>1) 
     edm::LogInfo("HGCalValidation") << "gx =  "  << hinfo.x
                     << " gy = "  << hinfo.y
                     << " gz = "  << hinfo.z;
   
   fillDigiInfo(hinfo);
 
   if (global1.eta() > 0)  fillOccupancyMap(OccupancyMap_plus_, layer -1);
   else                    fillOccupancyMap(OccupancyMap_minus_, layer -1);
   
 }

void HGCalDigiValidation::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 276 of file HGCalDigiValidation.cc.

References edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), HGCalDDDConstants::layers(), layers_, nameDetector_, and verbosity_.

                                                    {
   if (nameDetector_ == "HCal") {
     edm::ESHandle<HcalDDDRecConstants> pHRNDC;
     iSetup.get<HcalRecNumberingRecord>().get( pHRNDC );
     const HcalDDDRecConstants *hcons   = &(*pHRNDC);
     layers_ = hcons->getMaxDepth(1);
   } else {
     edm::ESHandle<HGCalDDDConstants>  pHGDC;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, pHGDC);
     const HGCalDDDConstants & hgcons_ = (*pHGDC);
     layers_ = hgcons_.layers(true);
   }
   
   if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << "current DQM directory:  "
                     << "HGCalDigiV/" << nameDetector_ 
                     << "  layer = "<< layers_;
 }  

void HGCalDigiValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 327 of file HGCalDigiValidation.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                      {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

void HGCalDigiValidation::fillDigiInfo ( digiInfo & hinfo )

private

Definition at line 254 of file HGCalDigiValidation.cc.

References HGCalDigiValidation::digiInfo::adc, ADC_, HGCalDigiValidation::digiInfo::charge, charge_, DigiOccupancy_XY_, HGCalDigiValidation::digiInfo::layer, HGCalDigiValidation::digiInfo::x, and HGCalDigiValidation::digiInfo::y.

                                                       {
   int ilayer = hinfo.layer -1;
   charge_.at(ilayer)->Fill(hinfo.charge);
   DigiOccupancy_XY_.at(ilayer)->Fill(hinfo.x, hinfo.y);
   ADC_.at(ilayer)->Fill(hinfo.adc);
 }

void HGCalDigiValidation::fillDigiInfo ( )

private

Definition at line 261 of file HGCalDigiValidation.cc.

References DigiOccupancy_Minus_, DigiOccupancy_Plus_, OccupancyMap_minus_, and OccupancyMap_plus_.

Referenced by analyze(), and digiValidation().

                                        {
   for (auto itr = OccupancyMap_plus_.begin(); 
        itr != OccupancyMap_plus_.end(); ++itr) {
     int layer = (*itr).first;
     int occupancy = (*itr).second;
     DigiOccupancy_Plus_.at(layer)->Fill(occupancy);
   }
   for (auto itr = OccupancyMap_minus_.begin(); 
        itr != OccupancyMap_minus_.end(); ++itr) {
     int layer = (*itr).first;
     int occupancy = (*itr).second;
     DigiOccupancy_Minus_.at(layer)->Fill(occupancy);
   }
 }

void HGCalDigiValidation::fillOccupancyMap	(	std::map< int, int > &	OccupancyMap,
		int	layer
	)

private

Definition at line 249 of file HGCalDigiValidation.cc.

Referenced by digiValidation().

                                                                                     {
   if (OccupancyMap.find(layer) != OccupancyMap.end()) OccupancyMap[layer] ++;
   else                                                OccupancyMap[layer] = 1;
 }