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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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_

edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geometry_beginRun_token_

edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geometry_token_

std::vector< MonitorElement * >	HitOccupancy_Minus_

std::vector< MonitorElement * >	HitOccupancy_Plus_

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 >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 33 of file HGCalRecHitValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 84 of file HGCalRecHitValidation.cc.

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

     : nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
       geometry_token_(esConsumes(edm::ESInputTag("", nameDetector_))),
       geometry_beginRun_token_(esConsumes<HGCalGeometry, IdealGeometryRecord, edm::Transition::BeginRun>(
           edm::ESInputTag("", nameDetector_))),
       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 {
     throw cms::Exception("BadHGCRecHitSource") << "HGCal DetectorName given as " << nameDetector_ << " must be: "
                                                << "\"HGCalHESiliconSensitive\", \"HGCalHESiliconSensitive\", or "
                                                << "\"HGCalHEScintillatorSensitive\"!";
   }
 }

HGCalRecHitValidation::~HGCalRecHitValidation ( )

inlineoverride

Definition at line 50 of file HGCalRecHitValidation.cc.

50 {}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 110 of file HGCalRecHitValidation.cc.

References edm::EventID::event(), fillHitsInfo(), motionPolicyChamber::geom0, geometry_token_, edm::Event::getByToken(), edm::EventSetup::getHandle(), edm::EventBase::id(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), HGCalDetId::layer(), HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), phase1PixelTopology::layer, nameDetector_, gpuPixelDoublets::ntot, OccupancyMap_minus, OccupancyMap_plus, convertSQLiteXML::ok, edm::Handle< T >::product(), edm::ESHandle< class >::product(), recHitSource_, recHitValidation(), HGCalTopology::tileTrapezoid(), HGCalGeometry::topology(), verbosity_, and HGCalTopology::waferHexagon8().

                                                                                        {
   OccupancyMap_plus.clear();
   OccupancyMap_minus.clear();
 
   bool ok(true);
   unsigned int ntot(0), nused(0);
   edm::ESHandle<HGCalGeometry> geomHandle = iSetup.getHandle(geometry_token_);
   if (!geomHandle.isValid()) {
     edm::LogVerbatim("HGCalValidation") << "Cannot get valid HGCalGeometry "
                                         << "Object for " << nameDetector_;
   } else {
     const HGCalGeometry* geom0 = geomHandle.product();
     int geomType = ((geom0->topology().waferHexagon8()) ? 1 : ((geom0->topology().tileTrapezoid()) ? 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 &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 225 of file HGCalRecHitValidation.cc.

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

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

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 212 of file HGCalRecHitValidation.cc.

References HGCalTopology::dddConstants(), HGCalDDDConstants::firstLayer(), firstLayer_, relativeConstraints::geom, geometry_beginRun_token_, edm::EventSetup::getHandle(), edm::ESHandleBase::isValid(), HGCalDDDConstants::layers(), layers_, nameDetector_, edm::ESHandle< class >::product(), and HGCalGeometry::topology().

                                                                                   {
   edm::ESHandle<HGCalGeometry> geomHandle = iSetup.getHandle(geometry_beginRun_token_);
   if (!geomHandle.isValid()) {
     edm::LogVerbatim("HGCalValidation") << "Cannot get valid HGCalGeometry "
                                         << "Object for " << nameDetector_;
   } else {
     const HGCalGeometry* geom = geomHandle.product();
     const HGCalDDDConstants& hgcons_ = geom->topology().dddConstants();
     layers_ = hgcons_.layers(true);
     firstLayer_ = hgcons_.firstLayer();
   }
 }

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

static

Definition at line 102 of file HGCalRecHitValidation.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

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

void HGCalRecHitValidation::fillHitsInfo ( )

private

Definition at line 190 of file HGCalRecHitValidation.cc.

References HitOccupancy_Minus_, HitOccupancy_Plus_, phase1PixelTopology::layer, OccupancyMap_minus, and OccupancyMap_plus.

Referenced by analyze(), and recHitValidation().

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

References phase1PixelTopology::layer.

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

References relval_parameters_module::energy, 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().

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