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 &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

	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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

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

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_

int	firstLayer_

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

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

Constructor & Destructor Documentation

◆ HGCalRecHitValidation()

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

explicit

Definition at line 91 of file HGCalRecHitValidation.cc.

     : nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
       ifHCAL_(iConfig.getParameter<bool>("ifHCAL")),
       verbosity_(iConfig.getUntrackedParameter<int>("Verbosity", 0)),
       firstLayer_(1) {
   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\"!";
   }
 }

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

◆ ~HGCalRecHitValidation()

HGCalRecHitValidation::~HGCalRecHitValidation ( )

inlineoverride

Definition at line 58 of file HGCalRecHitValidation.cc.

58 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 121 of file HGCalRecHitValidation.cc.

                                                                                        {
   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();
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << "Event " << iEvent.id().event() << " with " << ntot << " total and " << nused
                                         << " used recHits";
 }

References HcalDetId::depth(), fillHitsInfo(), relativeConstraints::geom, motionPolicyChamber::geom0, edm::EventSetup::get(), get, HcalEndcap, HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, iEvent, ifHCAL_, edm::HandleBase::isValid(), HGCalDetId::layer(), HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, nameDetector_, OccupancyMap_minus, OccupancyMap_plus, convertSQLiteXML::ok, edm::Handle< T >::product(), recHitSource_, recHitValidation(), edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), HGCalGeometryMode::Trapezoid, and verbosity_.

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 292 of file HGCalRecHitValidation.cc.

                                                                                                    {
   iB.setCurrentFolder("HGCAL/HGCalRecHitsV/" + nameDetector_);
   std::ostringstream histoname;
   for (unsigned int il = 0; il < layers_; ++il) {
     int ilayer = firstLayer_ + (int)(il);
     auto istr1 = std::to_string(ilayer);
     while (istr1.size() < 2) {
       istr1.insert(0, "0");
     }
     histoname.str("");
     histoname << "HitOccupancy_Plus_layer_" << istr1;
     HitOccupancy_Plus_.push_back(iB.book1D(histoname.str().c_str(), "RecHitOccupancy_Plus", 100, 0, 10000));
     histoname.str("");
     histoname << "HitOccupancy_Minus_layer_" << istr1;
     HitOccupancy_Minus_.push_back(iB.book1D(histoname.str().c_str(), "RecHitOccupancy_Minus", 100, 0, 10000));
  
     histoname.str("");
     histoname << "EtaPhi_Plus_"
               << "layer_" << istr1;
     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_" << istr1;
     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_" << istr1;
     energy_.push_back(iB.book1D(histoname.str().c_str(), "energy_", 500, 0, 1));
   }  //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);
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), energy_, EtaPhi_Minus_, EtaPhi_Plus_, firstLayer_, HitOccupancy_Minus_, HitOccupancy_Plus_, createfilelist::int, layers_, MeanHitOccupancy_Minus_, MeanHitOccupancy_Plus_, nameDetector_, pi, and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 277 of file HGCalRecHitValidation.cc.

                                                                                   {
   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);
     firstLayer_ = hgcons_.firstLayer();
   }
 }

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

◆ fillDescriptions()

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

static

Definition at line 112 of file HGCalRecHitValidation.cc.

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

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

◆ fillHitsInfo() [1/2]

void HGCalRecHitValidation::fillHitsInfo ( )

private

Definition at line 255 of file HGCalRecHitValidation.cc.

                                          {
   for (auto const& itr : OccupancyMap_plus) {
     int layer = itr.first;
     int occupancy = itr.second;
     HitOccupancy_Plus_.at(layer)->Fill(occupancy);
   }
  
   for (auto const& itr : OccupancyMap_minus) {
     int layer = itr.first;
     int occupancy = itr.second;
     HitOccupancy_Minus_.at(layer)->Fill(occupancy);
   }
 }

References HitOccupancy_Minus_, HitOccupancy_Plus_, OccupancyMap_minus, and OccupancyMap_plus.

Referenced by analyze(), and recHitValidation().

◆ fillHitsInfo() [2/2]

void HGCalRecHitValidation::fillHitsInfo ( HitsInfo & hits )

private

Definition at line 269 of file HGCalRecHitValidation.cc.

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

References energy_, EtaPhi_Minus_, EtaPhi_Plus_, and hfClusterShapes_cfi::hits.

◆ fillOccupancyMap()

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

private

Definition at line 248 of file HGCalRecHitValidation.cc.

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

Referenced by recHitValidation().

◆ recHitValidation()

template<class T1 , class T2 >

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

private

Definition at line 218 of file HGCalRecHitValidation.cc.

                                                                                            {
   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 = layer - firstLayer_;
   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);
   else
     fillOccupancyMap(OccupancyMap_minus, hinfo.layer);
 }

References HCALHighEnergyHPDFilter_cfi::energy, PV3DBase< T, PVType, FrameType >::eta(), fillHitsInfo(), fillOccupancyMap(), firstLayer_, relativeConstraints::geom, OccupancyMap_minus, OccupancyMap_plus, PV3DBase< T, PVType, FrameType >::phi(), verbosity_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

Member Data Documentation

◆ energy_

std::vector<MonitorElement*> HGCalRecHitValidation::energy_

private

Definition at line 84 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ EtaPhi_Minus_

std::vector<MonitorElement*> HGCalRecHitValidation::EtaPhi_Minus_

private

Definition at line 83 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ EtaPhi_Plus_

std::vector<MonitorElement*> HGCalRecHitValidation::EtaPhi_Plus_

private

Definition at line 82 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ firstLayer_

int HGCalRecHitValidation::firstLayer_

private

Definition at line 78 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), dqmBeginRun(), and recHitValidation().

◆ HitOccupancy_Minus_

std::vector<MonitorElement*> HGCalRecHitValidation::HitOccupancy_Minus_

private

Definition at line 86 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ HitOccupancy_Plus_

std::vector<MonitorElement*> HGCalRecHitValidation::HitOccupancy_Plus_

private

Definition at line 85 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ ifHCAL_

bool HGCalRecHitValidation::ifHCAL_

private

Definition at line 75 of file HGCalRecHitValidation.cc.

Referenced by analyze(), and HGCalRecHitValidation().

◆ layers_

unsigned int HGCalRecHitValidation::layers_

private

Definition at line 77 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and dqmBeginRun().

◆ MeanHitOccupancy_Minus_

MonitorElement* HGCalRecHitValidation::MeanHitOccupancy_Minus_

private

Definition at line 88 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms().

◆ MeanHitOccupancy_Plus_

MonitorElement* HGCalRecHitValidation::MeanHitOccupancy_Plus_

private

Definition at line 87 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms().

◆ nameDetector_

std::string HGCalRecHitValidation::nameDetector_

private

Definition at line 73 of file HGCalRecHitValidation.cc.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and HGCalRecHitValidation().

◆ OccupancyMap_minus

std::map<int, int> HGCalRecHitValidation::OccupancyMap_minus

private

Definition at line 80 of file HGCalRecHitValidation.cc.

Referenced by analyze(), fillHitsInfo(), and recHitValidation().

◆ OccupancyMap_plus

std::map<int, int> HGCalRecHitValidation::OccupancyMap_plus

private

Definition at line 79 of file HGCalRecHitValidation.cc.

Referenced by analyze(), fillHitsInfo(), and recHitValidation().

◆ recHitSource_

edm::EDGetToken HGCalRecHitValidation::recHitSource_

private

Definition at line 74 of file HGCalRecHitValidation.cc.

Referenced by analyze(), and HGCalRecHitValidation().

◆ verbosity_

int HGCalRecHitValidation::verbosity_

private

Definition at line 76 of file HGCalRecHitValidation.cc.

Referenced by analyze(), and recHitValidation().

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalRecHitValidation()

◆ ~HGCalRecHitValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ fillDescriptions()

◆ fillHitsInfo() [1/2]

◆ fillHitsInfo() [2/2]

◆ fillOccupancyMap()

◆ recHitValidation()

Member Data Documentation

◆ energy_

◆ EtaPhi_Minus_

◆ EtaPhi_Plus_

◆ firstLayer_

◆ HitOccupancy_Minus_

◆ HitOccupancy_Plus_

◆ ifHCAL_

◆ layers_

◆ MeanHitOccupancy_Minus_

◆ MeanHitOccupancy_Plus_

◆ nameDetector_

◆ OccupancyMap_minus

◆ OccupancyMap_plus

◆ recHitSource_

◆ verbosity_