HGCalRecHitValidation Class Reference

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::HGCalRecHitValidation ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 84 of file HGCalRecHitValidation.cc.

    : 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\"!";
  }
}

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

~HGCalRecHitValidation()

HGCalRecHitValidation::~HGCalRecHitValidation ( )

inlineoverride

Definition at line 50 of file HGCalRecHitValidation.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 110 of file HGCalRecHitValidation.cc.

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

References fillHitsInfo(), motionPolicyChamber::geom0, geometry_token_, edm::EventSetup::getHandle(), iEvent, 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< T >::product(), recHitSource_, recHitValidation(), edm::SortedCollection< T, SORT >::size(), and verbosity_.

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

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

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

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

fillDescriptions()

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

static

Definition at line 102 of file HGCalRecHitValidation.cc.

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

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

fillHitsInfo() [1/2]

void HGCalRecHitValidation::fillHitsInfo ( )

private

Definition at line 190 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_, phase1PixelTopology::layer, OccupancyMap_minus, and OccupancyMap_plus.

Referenced by analyze(), and recHitValidation().

fillHitsInfo() [2/2]

void HGCalRecHitValidation::fillHitsInfo ( HitsInfo &  hits )

private

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

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

References phase1PixelTopology::layer.

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 153 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, phase1PixelTopology::layer, 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 77 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

EtaPhi_Minus_

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

private

Definition at line 76 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

EtaPhi_Plus_

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

private

Definition at line 75 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

firstLayer_

int HGCalRecHitValidation::firstLayer_

private

Definition at line 71 of file HGCalRecHitValidation.cc.

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

geometry_beginRun_token_

edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalRecHitValidation::geometry_beginRun_token_

private

Definition at line 67 of file HGCalRecHitValidation.cc.

Referenced by dqmBeginRun().

geometry_token_

edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalRecHitValidation::geometry_token_

private

Definition at line 66 of file HGCalRecHitValidation.cc.

Referenced by analyze().

HitOccupancy_Minus_

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

private

Definition at line 79 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

HitOccupancy_Plus_

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

private

Definition at line 78 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

layers_

unsigned int HGCalRecHitValidation::layers_

private

Definition at line 70 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms(), and dqmBeginRun().

MeanHitOccupancy_Minus_

MonitorElement* HGCalRecHitValidation::MeanHitOccupancy_Minus_

private

Definition at line 81 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms().

MeanHitOccupancy_Plus_

MonitorElement* HGCalRecHitValidation::MeanHitOccupancy_Plus_

private

Definition at line 80 of file HGCalRecHitValidation.cc.

Referenced by bookHistograms().

nameDetector_

std::string HGCalRecHitValidation::nameDetector_

private

Definition at line 65 of file HGCalRecHitValidation.cc.

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

OccupancyMap_minus

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

private

Definition at line 73 of file HGCalRecHitValidation.cc.

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

OccupancyMap_plus

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

private

Definition at line 72 of file HGCalRecHitValidation.cc.

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

recHitSource_

edm::EDGetToken HGCalRecHitValidation::recHitSource_

private

Definition at line 68 of file HGCalRecHitValidation.cc.

Referenced by analyze(), and HGCalRecHitValidation().

verbosity_

int HGCalRecHitValidation::verbosity_

private

Definition at line 69 of file HGCalRecHitValidation.cc.

Referenced by analyze(), and recHitValidation().