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

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

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

virtual void	beginStream (edm::StreamID id) final

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

	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
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::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 44 of file HGCalRecHitValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 41 of file HGCalRecHitValidation.cc.

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

                                                       :
   nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
   verbosity_(iConfig.getUntrackedParameter<int>("Verbosity",0)) {
   auto temp = iConfig.getParameter<edm::InputTag>("RecHitSource");
   ifHCAL_   = iConfig.getParameter<bool>("ifHCAL");
   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 ( )

Definition at line 62 of file HGCalRecHitValidation.cc.

62 { }

Member Function Documentation

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

override

Definition at line 64 of file HGCalRecHitValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), edm::EventID::event(), fillHitsInfo(), relativeConstraints::geom, motionPolicyChamber::geom0, edm::EventSetup::get(), edm::Event::getByToken(), HcalEndcap, HGCEE, edm::EventBase::id(), ifHCAL_, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), HGCEEDetId::layer(), HGCHEDetId::layer(), nameDetector_, OccupancyMap_minus, OccupancyMap_plus, convertSQLiteXML::ok, edm::ESHandle< T >::product(), recHitSource_, recHitValidation(), edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), and verbosity_.

                                                  {
   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::LogWarning("HGCalValidation") << "Cannot get valid HGCalGeometry Object for " << nameDetector_;
     const CaloGeometry* geom0 = geom.product();
 
     if (ifHCAL_) {
       edm::Handle<HBHERecHitCollection> hbhecoll;
       iEvent.getByToken(recHitSource_, hbhecoll);
       if (hbhecoll.isValid()) {
     if (verbosity_>0) 
       edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                       << hbhecoll->size() << " element(s)";
     for (HBHERecHitCollection::const_iterator it=hbhecoll->begin(); 
          it != hbhecoll->end(); ++it) {
       DetId detId = it->id();
       ntot++;
       if (detId.subdetId() == HcalEndcap) {
         nused++;
         int   layer = HcalDetId(detId).depth();
         recHitValidation(detId, layer, geom0, it);
       }
     }
       } else {
     ok = false;
     edm::LogWarning("HGCalValidation") << "HBHERecHitCollection Handle does not exist !!!";
       }
     } else {
       edm::Handle<HGChebRecHitCollection> hbhecoll;
       iEvent.getByToken(recHitSource_, hbhecoll);
       if (hbhecoll.isValid()) {
     if (verbosity_>0) 
       edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                       << hbhecoll->size() << " element(s)";
     for (HGChebRecHitCollection::const_iterator it=hbhecoll->begin(); 
          it != hbhecoll->end(); ++it) {
       DetId detId = it->id();
       ntot++; nused++;
       int   layer = HcalDetId(detId).depth();
       recHitValidation(detId, layer, geom0, it);
     }
       } else {
     ok = false;
     edm::LogWarning("HGCalValidation") << "HGChebRecHitCollection Handle does not exist !!!";
       }
     }
   } else {
     edm::ESHandle<HGCalGeometry> geom;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, geom);
     if (!geom.isValid()) edm::LogWarning("HGCalValidation") << "Cannot get valid HGCalGeometry Object for " << nameDetector_;
     const HGCalGeometry* geom0 = geom.product();
 
     edm::Handle<HGCRecHitCollection> theRecHitContainers;
     iEvent.getByToken(recHitSource_, theRecHitContainers);
     if (theRecHitContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogInfo("HGCalValidation") << nameDetector_ << " with " 
                     << theRecHitContainers->size()
                     << " element(s)";
       for (HGCRecHitCollection::const_iterator it=theRecHitContainers->begin();
        it !=theRecHitContainers->end(); ++it) {
     ntot++; nused++;
     DetId detId = it->id();
     int layer   = (detId.subdetId() == HGCEE) ? (HGCEEDetId(detId).layer()) : (HGCHEDetId(detId).layer());
     recHitValidation(detId, layer, geom0, it);
       }
     } else {
       ok = false;
       edm::LogWarning("HGCalValidation") << "HGCRecHitCollection Handle does not exist !!!";
     }
   }
   if (ok) fillHitsInfo();
   edm::LogWarning("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 228 of file HGCalRecHitValidation.cc.

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

                                              {
 
   iB.setCurrentFolder("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 212 of file HGCalRecHitValidation.cc.

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

                                                    {
 
   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 257 of file HGCalRecHitValidation.cc.

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

                                                                                        {
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

void HGCalRecHitValidation::fillHitsInfo ( )

private

Definition at line 186 of file HGCalRecHitValidation.cc.

References HitOccupancy_Minus_, HitOccupancy_Plus_, OccupancyMap_minus, and OccupancyMap_plus.

Referenced by analyze(), and recHitValidation().

                                          { 
 
   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 202 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 181 of file HGCalRecHitValidation.cc.

Referenced by recHitValidation().

                                                                                      {
   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 149 of file HGCalRecHitValidation.cc.

References HGCalRecHitValidation::HitsInfo::energy, HGCalRecHitValidation::HitsInfo::eta, PV3DBase< T, PVType, FrameType >::eta(), fillHitsInfo(), fillOccupancyMap(), HGCalRecHitValidation::HitsInfo::layer, 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().

                                                 {
 
   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  = layer;
   hinfo.phi    = global.phi();
   hinfo.eta    = global.eta();
       
   if (verbosity_>1) 
     edm::LogInfo("HGCalValidation") << " --------------------------   gx = "
                     << globalx << " gy = "  << globaly   
                     << " gz = " << globalz << " phi = " 
                     << hinfo.phi << " eta = " << hinfo.eta;
       
   fillHitsInfo(hinfo);
       
   if (hinfo.eta > 0)  fillOccupancyMap(OccupancyMap_plus, layer -1);
   else                fillOccupancyMap(OccupancyMap_minus, layer -1);
       
 }      