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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::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_

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_

Detailed Description

Definition at line 41 of file HGCalRecHitValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 93 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)),
   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\"!"; 
   }
 }

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 125 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(), HGCalDetId::layer(), HGCScintillatorDetId::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();
   if (verbosity_>0) 
     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 &
	)

override

Definition at line 308 of file HGCalRecHitValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DEFINE_FWK_MODULE, energy_, EtaPhi_Minus_, EtaPhi_Plus_, firstLayer_, HitOccupancy_Minus_, HitOccupancy_Plus_, createfilelist::int, 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 il=0; il<layers_; ++il) {
     int ilayer = firstLayer_ + (int)(il);
     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_",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);
 }

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

override

Definition at line 291 of file HGCalRecHitValidation.cc.

References HGCalDDDConstants::firstLayer(), firstLayer_, 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);
     firstLayer_ = hgcons_.firstLayer();
   }
 }

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

static

Definition at line 116 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 265 of file HGCalRecHitValidation.cc.

References HitOccupancy_Minus_, HitOccupancy_Plus_, OccupancyMap_minus, and OccupancyMap_plus.

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

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

void HGCalRecHitValidation::fillHitsInfo ( HitsInfo & hits )

private

Definition at line 281 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;
   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 259 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 226 of file HGCalRecHitValidation.cc.

References HGCalRecHitValidation::HitsInfo::energy, HGCalRecHitValidation::HitsInfo::eta, PV3DBase< T, PVType, FrameType >::eta(), fillHitsInfo(), fillOccupancyMap(), firstLayer_, 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(), 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  = 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);
       
 }      