Inheritance diagram for HGCalRecHitValidation:

Classes
struct	energysum

struct	HitsInfo

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

	HGCalRecHitValidation (const edm::ParameterSet &)

	~HGCalRecHitValidation () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) 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 edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	fillHitsInfo ()

void	fillHitsInfo (HitsInfo &hits)

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

template<class T1 , class T2 >
void	recHitValidation (DetId &detId, int layer, const T1 *geom, T2 it)

Private Attributes
std::vector< MonitorElement * >	energy_

std::vector< MonitorElement * >	EtaPhi_Minus_

std::vector< MonitorElement * >	EtaPhi_Plus_

std::vector< MonitorElement * >	HitOccupancy_Minus_

std::vector< MonitorElement * >	HitOccupancy_Plus_

bool	ifHCAL_

unsigned int	layers_

MonitorElement *	MeanHitOccupancy_Minus_

MonitorElement *	MeanHitOccupancy_Plus_

std::string	nameDetector_

std::map< int, int >	OccupancyMap_minus

std::map< int, int >	OccupancyMap_plus

edm::EDGetToken	recHitSource_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 41 of file HGCalRecHitValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 92 of file HGCalRecHitValidation.cc.

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

                                                                           :
   nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
   ifHCAL_(iConfig.getParameter<bool>("ifHCAL")),
   verbosity_(iConfig.getUntrackedParameter<int>("Verbosity",0)) {
 
   auto temp = iConfig.getParameter<edm::InputTag>("RecHitSource");
   if (nameDetector_ == "HGCalEESensitive" || 
       nameDetector_ == "HGCalHESiliconSensitive" ||
       nameDetector_ == "HGCalHEScintillatorSensitive") {
     recHitSource_    = consumes<HGCRecHitCollection>(temp);
   } else if (nameDetector_ == "HCal") {
     if (ifHCAL_) recHitSource_ = consumes<HBHERecHitCollection>(temp);
     else         recHitSource_ = consumes<HGChebRecHitCollection>(temp);
   } else {
     throw cms::Exception("BadHGCRecHitSource")
       << "HGCal DetectorName given as " << nameDetector_ << " must be: "
       << "\"HGCalHESiliconSensitive\", \"HGCalHESiliconSensitive\", "
       << "\"HGCalHEScintillatorSensitive\", or \"HCal\"!"; 
   }
 }

HGCalRecHitValidation::~HGCalRecHitValidation ( )

inlineoverride

Definition at line 60 of file HGCalRecHitValidation.cc.

References analyze(), bookHistograms(), dqmBeginRun(), fillDescriptions(), fillHitsInfo(), fillOccupancyMap(), relativeConstraints::geom, hfClusterShapes_cfi::hits, and recHitValidation().

60 {}

Member Function Documentation

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

override

Definition at line 123 of file HGCalRecHitValidation.cc.

References HcalDetId::depth(), edm::EventID::event(), fillHitsInfo(), relativeConstraints::geom, motionPolicyChamber::geom0, HGCalTopology::geomMode(), edm::EventSetup::get(), edm::Event::getByToken(), HcalEndcap, HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, edm::EventBase::id(), ifHCAL_, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), HGCScintillatorDetId::layer(), HGCalDetId::layer(), HGCSiliconDetId::layer(), ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, nameDetector_, OccupancyMap_minus, OccupancyMap_plus, convertSQLiteXML::ok, edm::Handle< T >::product(), edm::ESHandle< T >::product(), recHitSource_, recHitValidation(), edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), HGCalGeometry::topology(), HGCalGeometryMode::Trapezoid, and verbosity_.

Referenced by ~HGCalRecHitValidation().

                                                  {
   OccupancyMap_plus.clear();
   OccupancyMap_minus.clear();
 
   bool ok(true);
   unsigned int ntot(0), nused(0);
   if (nameDetector_ == "HCal") {
     edm::ESHandle<CaloGeometry> geom;
     iSetup.get<CaloGeometryRecord>().get(geom);
     if (!geom.isValid()) {
       edm::LogVerbatim("HGCalValidation") << "Cannot get valid HGCalGeometry "
                       << "Object for " << nameDetector_;
     } else {
       const CaloGeometry* geom0 = geom.product();
 
       if (ifHCAL_) {
     edm::Handle<HBHERecHitCollection> hbhecoll;
     iEvent.getByToken(recHitSource_, hbhecoll);
     if (hbhecoll.isValid()) {
       if (verbosity_>0) 
         edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << hbhecoll->size() 
                         << " element(s)";
       for (const auto & it : *(hbhecoll.product())) {
         DetId detId = it.id();
         ntot++;
         if (detId.subdetId() == HcalEndcap) {
           nused++;
           int   layer = HcalDetId(detId).depth();
           recHitValidation(detId, layer, geom0, &it);
         }
       }
     } else {
       ok = false;
       edm::LogVerbatim("HGCalValidation") << "HBHERecHitCollection Handle does not "
                           << "exist !!!";
     }
       } else {
     edm::Handle<HGChebRecHitCollection> hbhecoll;
     iEvent.getByToken(recHitSource_, hbhecoll);
     if (hbhecoll.isValid()) {
       if (verbosity_>0) 
         edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << hbhecoll->size() 
                         << " element(s)";
       for (const auto & it : *(hbhecoll.product())) {
         DetId detId = it.id();
         ntot++; nused++;
         int   layer = HcalDetId(detId).depth();
         recHitValidation(detId, layer, geom0, &it);
       }
     } else {
       ok = false;
       edm::LogVerbatim("HGCalValidation") << "HGChebRecHitCollection Handle does not "
                           << "exist !!!";
     }
       }
     }
   } else {
     edm::ESHandle<HGCalGeometry> geom;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, geom);
     if (!geom.isValid()) {
       edm::LogVerbatim("HGCalValidation") << "Cannot get valid HGCalGeometry Object for "
                        << nameDetector_;
     } else {
       const HGCalGeometry* geom0 = geom.product();
       HGCalGeometryMode::GeometryMode mode = geom0->topology().geomMode();
       int geomType = (((mode == HGCalGeometryMode::Hexagon8) ||
                (mode == HGCalGeometryMode::Hexagon8Full)) ? 1 :
               ((mode == HGCalGeometryMode::Trapezoid) ? 2 : 0));
       
       edm::Handle<HGCRecHitCollection> theRecHitContainers;
       iEvent.getByToken(recHitSource_, theRecHitContainers);
       if (theRecHitContainers.isValid()) {
     if (verbosity_>0) 
       edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                           << theRecHitContainers->size()
                           << " element(s)";
     for (const auto &  it : *(theRecHitContainers.product())) {
       ntot++; nused++;
       DetId detId = it.id();
       int layer   = ((geomType == 0) ? HGCalDetId(detId).layer() :
              ((geomType == 1) ? HGCSiliconDetId(detId).layer() :
               HGCScintillatorDetId(detId).layer()));
       recHitValidation(detId, layer, geom0, &it);
     }
       } else {
     ok = false;
     edm::LogVerbatim("HGCalValidation") << "HGCRecHitCollection Handle does not "
                         << "exist !!!";
       }
     }
   }
   if (ok) fillHitsInfo();
   edm::LogVerbatim("HGCalValidation") << "Event " << iEvent.id().event()
                       << " with " << ntot << " total and "
                       << nused << " used recHits";
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 303 of file HGCalRecHitValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DEFINE_FWK_MODULE, energy_, EtaPhi_Minus_, EtaPhi_Plus_, HitOccupancy_Minus_, HitOccupancy_Plus_, layers_, MeanHitOccupancy_Minus_, MeanHitOccupancy_Plus_, nameDetector_, pi, and DQMStore::IBooker::setCurrentFolder().

Referenced by ~HGCalRecHitValidation().

                                              {
 
   iB.setCurrentFolder("HGCAL/HGCalRecHitsV/"+nameDetector_);
   std::ostringstream histoname;
   for (unsigned int ilayer = 0; ilayer < layers_; ilayer++ ) {
     histoname.str(""); histoname << "HitOccupancy_Plus_layer_" << ilayer;
     HitOccupancy_Plus_.push_back(iB.book1D( histoname.str().c_str(), "RecHitOccupancy_Plus", 100, 0, 10000));
     histoname.str(""); histoname << "HitOccupancy_Minus_layer_" << ilayer;
     HitOccupancy_Minus_.push_back(iB.book1D( histoname.str().c_str(), "RecHitOccupancy_Minus", 100, 0, 10000));
 
     histoname.str(""); histoname << "EtaPhi_Plus_" << "layer_" << ilayer;
     EtaPhi_Plus_.push_back(iB.book2D(histoname.str().c_str(), "Occupancy", 31, 1.45, 3.0, 72, -CLHEP::pi, CLHEP::pi));
     histoname.str(""); histoname << "EtaPhi_Minus_" << "layer_" << ilayer;
     EtaPhi_Minus_.push_back(iB.book2D(histoname.str().c_str(), "Occupancy", 31, -3.0, -1.45, 72, -CLHEP::pi, CLHEP::pi));
       
     histoname.str(""); histoname << "energy_layer_" << ilayer; 
     energy_.push_back(iB.book1D(histoname.str().c_str(),"energy_",100,0,0.002));
   }//loop over layers ends here 
 
   histoname.str(""); histoname << "SUMOfRecHitOccupancy_Plus";
   MeanHitOccupancy_Plus_= iB.book1D( histoname.str().c_str(), "SUMOfRecHitOccupancy_Plus", layers_, -0.5, layers_-0.5);
   histoname.str(""); histoname << "SUMOfRecHitOccupancy_Minus";
   MeanHitOccupancy_Minus_ = iB.book1D( histoname.str().c_str(), "SUMOfRecHitOccupancy_Minus", layers_, -0.5,layers_-0.5);
 }

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

override

Definition at line 287 of file HGCalRecHitValidation.cc.

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

Referenced by ~HGCalRecHitValidation().

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

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

static

Definition at line 114 of file HGCalRecHitValidation.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~HGCalRecHitValidation().

                                                                                        {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("DetectorName","HGCalEESensitive");
   desc.add<edm::InputTag>("RecHitSource",edm::InputTag("HGCalRecHit","HGCEERecHits"));
   desc.add<bool>("ifHCAL",false);
   desc.addUntracked<int>("Verbosity",0);
   descriptions.add("hgcalRecHitValidationEE",desc);
 }

void HGCalRecHitValidation::fillHitsInfo ( )

private

Definition at line 261 of file HGCalRecHitValidation.cc.

References HitOccupancy_Minus_, HitOccupancy_Plus_, OccupancyMap_minus, and OccupancyMap_plus.

Referenced by analyze(), recHitValidation(), and ~HGCalRecHitValidation().

                                          { 
 
   for (auto itr = OccupancyMap_plus.begin() ; itr != OccupancyMap_plus.end(); ++itr) {
     int layer      = (*itr).first;
     int occupancy  = (*itr).second;
     HitOccupancy_Plus_.at(layer)->Fill(occupancy);
   }
 
   for (auto itr = OccupancyMap_minus.begin() ; itr != OccupancyMap_minus.end(); ++itr) {
     int layer      = (*itr).first;
     int occupancy  = (*itr).second;
     HitOccupancy_Minus_.at(layer)->Fill(occupancy);
   }
   
 }

void HGCalRecHitValidation::fillHitsInfo ( HitsInfo & hits )

private

Definition at line 277 of file HGCalRecHitValidation.cc.

References HGCalRecHitValidation::HitsInfo::energy, energy_, HGCalRecHitValidation::HitsInfo::eta, EtaPhi_Minus_, EtaPhi_Plus_, HGCalRecHitValidation::HitsInfo::layer, and HGCalRecHitValidation::HitsInfo::phi.

                                                        {
 
   unsigned int ilayer = hits.layer -1;
   energy_.at(ilayer)->Fill(hits.energy);
 
   EtaPhi_Plus_.at(ilayer)->Fill(hits.eta , hits.phi);
   EtaPhi_Minus_.at(ilayer)->Fill(hits.eta, hits.phi);
 
 }

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

private

Definition at line 256 of file HGCalRecHitValidation.cc.

Referenced by recHitValidation(), and ~HGCalRecHitValidation().

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

template<class T1 , class T2 >

void HGCalRecHitValidation::recHitValidation	(	DetId &	detId,
		int	layer,
		const T1 *	geom,
		T2	it
	)

private

Definition at line 224 of file HGCalRecHitValidation.cc.

References HGCalRecHitValidation::HitsInfo::energy, HGCalRecHitValidation::HitsInfo::eta, PV3DBase< T, PVType, FrameType >::eta(), fillHitsInfo(), fillOccupancyMap(), HGCalRecHitValidation::HitsInfo::layer, layers_, min(), OccupancyMap_minus, OccupancyMap_plus, HGCalRecHitValidation::HitsInfo::phi, PV3DBase< T, PVType, FrameType >::phi(), verbosity_, HGCalRecHitValidation::HitsInfo::x, PV3DBase< T, PVType, FrameType >::x(), HGCalRecHitValidation::HitsInfo::y, PV3DBase< T, PVType, FrameType >::y(), HGCalRecHitValidation::HitsInfo::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze(), and ~HGCalRecHitValidation().

                                                 {
 
   const GlobalPoint& global = geom->getPosition(detId);
   double      energy = it->energy();
 
   float globalx = global.x();
   float globaly = global.y();
   float globalz = global.z();
       
   HitsInfo   hinfo;
   hinfo.energy = energy;
   hinfo.x      = globalx;
   hinfo.y      = globaly;
   hinfo.z      = globalz;
   hinfo.layer  = std::min((unsigned)layer,layers_);
   hinfo.phi    = global.phi();
   hinfo.eta    = global.eta();
       
   if (verbosity_>1) 
     edm::LogVerbatim("HGCalValidation") << "--------------------------   gx = "
                     << globalx << " gy = "  << globaly   
                     << " gz = " << globalz << " phi = " 
                     << hinfo.phi << " eta = " << hinfo.eta << " lay = " << hinfo.layer;
       
   fillHitsInfo(hinfo);
       
   if (hinfo.eta > 0)  fillOccupancyMap(OccupancyMap_plus, hinfo.layer -1);
   else                fillOccupancyMap(OccupancyMap_minus, hinfo.layer -1);
       
 }      