HGCalHitValidation Class Reference

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

Inheritance diagram for HGCalHitValidation:

Public Member Functions
	HGCalHitValidation (const edm::ParameterSet &)
	~HGCalHitValidation () 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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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 &)

Protected Types
typedef std::tuple< float, float, float, float >	HGCHitTuple

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
template<class T1 >
void	analyzeHGCalRecHit (T1 const &theHits, std::map< unsigned int, HGCHitTuple > const &hitRefs)
void	analyzeHGCalSimHit (edm::Handle< std::vector< PCaloHit >> const &simHits, int idet, MonitorElement *hist, std::map< unsigned int, HGCHitTuple > &)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
edm::EDGetTokenT< HGChebRecHitCollection >	bhRecHitTokeng_
edm::EDGetTokenT< HBHERecHitCollection >	bhRecHitTokenh_
edm::InputTag	bhSimHitSource
edm::EDGetTokenT< std::vector< PCaloHit > >	bhSimHitToken_
edm::EDGetTokenT< HGCeeRecHitCollection >	eeRecHitToken_
edm::InputTag	eeSimHitSource
edm::EDGetTokenT< std::vector< PCaloHit > >	eeSimHitToken_
edm::EDGetTokenT< HGChefRecHitCollection >	fhRecHitToken_
edm::InputTag	fhSimHitSource
edm::EDGetTokenT< std::vector< PCaloHit > >	fhSimHitToken_
std::vector< std::string >	geometrySource_
const HcalDDDRecConstants *	hcConr_
const HcalDDDSimConstants *	hcCons_
const CaloSubdetectorGeometry *	hcGeometry_
MonitorElement *	hebdEtaVsEta
MonitorElement *	hebdPhiVsPhi
MonitorElement *	hebdzVsZ
MonitorElement *	hebEnRec
MonitorElement *	hebEnSim
MonitorElement *	hebEnSimRec
MonitorElement *	hebRecVsSimX
MonitorElement *	hebRecVsSimY
MonitorElement *	hebRecVsSimZ
MonitorElement *	heedxVsX
MonitorElement *	heedyVsY
MonitorElement *	heedzVsZ
MonitorElement *	heeEnRec
MonitorElement *	heeEnSim
MonitorElement *	heeEnSimRec
MonitorElement *	heeRecVsSimX
MonitorElement *	heeRecVsSimY
MonitorElement *	heeRecVsSimZ
MonitorElement *	hefdxVsX
MonitorElement *	hefdyVsY
MonitorElement *	hefdzVsZ
MonitorElement *	hefEnRec
MonitorElement *	hefEnSim
MonitorElement *	hefEnSimRec
MonitorElement *	hefRecVsSimX
MonitorElement *	hefRecVsSimY
MonitorElement *	hefRecVsSimZ
std::vector< const HGCalDDDConstants * >	hgcCons_
std::vector< const HGCalGeometry * >	hgcGeometry_
std::vector< int >	ietaExcludeBH_
bool	ifHCAL_
bool	ifHCALsim_

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 61 of file HGCalHitValidation.cc.

Member Typedef Documentation

HGCHitTuple

typedef std::tuple<float, float, float, float> HGCalHitValidation::HGCHitTuple

protected

Definition at line 68 of file HGCalHitValidation.cc.

Constructor & Destructor Documentation

HGCalHitValidation()

HGCalHitValidation::HGCalHitValidation ( const edm::ParameterSet &  cfg )

explicit

Definition at line 116 of file HGCalHitValidation.cc.

                                                                 {
  geometrySource_ = cfg.getUntrackedParameter<std::vector<std::string>>("geometrySource");
  eeSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("eeSimHitSource"));
  fhSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("fhSimHitSource"));
  bhSimHitToken_ = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("bhSimHitSource"));
  eeRecHitToken_ = consumes<HGCeeRecHitCollection>(cfg.getParameter<edm::InputTag>("eeRecHitSource"));
  fhRecHitToken_ = consumes<HGChefRecHitCollection>(cfg.getParameter<edm::InputTag>("fhRecHitSource"));
  ietaExcludeBH_ = cfg.getParameter<std::vector<int>>("ietaExcludeBH");
  ifHCAL_ = cfg.getParameter<bool>("ifHCAL");
  ifHCALsim_ = cfg.getParameter<bool>("ifHCALsim");
  if (ifHCAL_)
    bhRecHitTokenh_ = consumes<HBHERecHitCollection>(cfg.getParameter<edm::InputTag>("bhRecHitSource"));
  else
    bhRecHitTokeng_ = consumes<HGChebRecHitCollection>(cfg.getParameter<edm::InputTag>("bhRecHitSource"));
 
#ifdef EDM_ML_DEBUG
  edm::LogInfo("HGCalValid") << "Exclude the following " << ietaExcludeBH_.size() << " ieta values from BH plots (BH "
                             << ifHCAL_ << ") ";
  for (unsigned int k = 0; k < ietaExcludeBH_.size(); ++k)
    edm::LogInfo("HGCalValid") << " [" << k << "] " << ietaExcludeBH_[k];
#endif
}

References bhRecHitTokeng_, bhRecHitTokenh_, bhSimHitToken_, looper::cfg, eeRecHitToken_, eeSimHitToken_, fhRecHitToken_, fhSimHitToken_, geometrySource_, ietaExcludeBH_, ifHCAL_, ifHCALsim_, and dqmdumpme::k.

~HGCalHitValidation()

HGCalHitValidation::~HGCalHitValidation ( )

override

Definition at line 139 of file HGCalHitValidation.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 234 of file HGCalHitValidation.cc.

                                                                                    {
  std::map<unsigned int, HGCHitTuple> eeHitRefs, fhHitRefs, bhHitRefs;
 
  //Accesing ee simhits
  edm::Handle<std::vector<PCaloHit>> eeSimHits;
  iEvent.getByToken(eeSimHitToken_, eeSimHits);
 
  if (eeSimHits.isValid()) {
    analyzeHGCalSimHit(eeSimHits, 0, heeEnSim, eeHitRefs);
#ifdef EDM_ML_DEBUG
    for (std::map<unsigned int, HGCHitTuple>::iterator itr = eeHitRefs.begin(); itr != eeHitRefs.end(); ++itr) {
      int idx = std::distance(eeHitRefs.begin(), itr);
      edm::LogInfo("HGCalValid") << "EEHit[" << idx << "] " << std::hex << itr->first << std::dec << "; Energy "
                                 << std::get<0>(itr->second) << "; Position (" << std::get<1>(itr->second) << ", "
                                 << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ")";
    }
#endif
  } else {
    edm::LogVerbatim("HGCalValid") << "No EE SimHit Found ";
  }
 
  //Accesing fh simhits
  edm::Handle<std::vector<PCaloHit>> fhSimHits;
  iEvent.getByToken(fhSimHitToken_, fhSimHits);
  if (fhSimHits.isValid()) {
    analyzeHGCalSimHit(fhSimHits, 1, hefEnSim, fhHitRefs);
#ifdef EDM_ML_DEBUG
    for (std::map<unsigned int, HGCHitTuple>::iterator itr = fhHitRefs.begin(); itr != fhHitRefs.end(); ++itr) {
      int idx = std::distance(fhHitRefs.begin(), itr);
      edm::LogInfo("HGCalValid") << "FHHit[" << idx << "] " << std::hex << itr->first << std::dec << "; Energy "
                                 << std::get<0>(itr->second) << "; Position (" << std::get<1>(itr->second) << ", "
                                 << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ")";
    }
#endif
  } else {
    edm::LogVerbatim("HGCalValid") << "No FH SimHit Found ";
  }
 
  //Accessing bh simhits
  edm::Handle<std::vector<PCaloHit>> bhSimHits;
  iEvent.getByToken(bhSimHitToken_, bhSimHits);
  if (bhSimHits.isValid()) {
    if (ifHCALsim_) {
      for (std::vector<PCaloHit>::const_iterator simHit = bhSimHits->begin(); simHit != bhSimHits->end(); ++simHit) {
        int subdet, z, depth, eta, phi, lay;
        HcalTestNumbering::unpackHcalIndex(simHit->id(), subdet, z, depth, eta, phi, lay);
 
        if (subdet == static_cast<int>(HcalEndcap)) {
          HcalCellType::HcalCell cell = hcCons_->cell(subdet, z, lay, eta, phi);
          double zp = cell.rz / 10;
 
          HcalDDDRecConstants::HcalID idx = hcConr_->getHCID(subdet, eta, phi, lay, depth);
          int sign = (z == 0) ? (-1) : (1);
          zp *= sign;
          HcalDetId id = HcalDetId(HcalEndcap, sign * idx.eta, idx.phi, idx.depth);
 
          float energy = simHit->energy();
          float energySum(energy);
          if (bhHitRefs.count(id.rawId()) != 0)
            energySum += std::get<0>(bhHitRefs[id.rawId()]);
          hebEnSim->Fill(energy);
          if (std::find(ietaExcludeBH_.begin(), ietaExcludeBH_.end(), idx.eta) == ietaExcludeBH_.end()) {
            bhHitRefs[id.rawId()] = std::make_tuple(energySum, cell.eta, cell.phi, zp);
#ifdef EDM_ML_DEBUG
            edm::LogInfo("HGCalValid") << "Accept " << id << std::endl;
          } else {
            edm::LogInfo("HGCalValid") << "Reject " << id << std::endl;
#endif
          }
        }
      }
    } else {
      analyzeHGCalSimHit(bhSimHits, 2, hebEnSim, bhHitRefs);
    }
#ifdef EDM_ML_DEBUG
    for (std::map<unsigned int, HGCHitTuple>::iterator itr = bhHitRefs.begin(); itr != bhHitRefs.end(); ++itr) {
      int idx = std::distance(bhHitRefs.begin(), itr);
      edm::LogInfo("HGCalValid") << "BHHit[" << idx << "] " << std::hex << itr->first << std::dec << "; Energy "
                                 << std::get<0>(itr->second) << "; Position (" << std::get<1>(itr->second) << ", "
                                 << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ")";
    }
#endif
  } else {
    edm::LogVerbatim("HGCalValid") << "No BH SimHit Found ";
  }
 
  //accessing EE Rechit information
  edm::Handle<HGCeeRecHitCollection> eeRecHit;
  iEvent.getByToken(eeRecHitToken_, eeRecHit);
  if (eeRecHit.isValid()) {
    const HGCeeRecHitCollection* theHits = (eeRecHit.product());
    for (auto it = theHits->begin(); it != theHits->end(); ++it) {
      double energy = it->energy();
      heeEnRec->Fill(energy);
      std::map<unsigned int, HGCHitTuple>::const_iterator itr = eeHitRefs.find(it->id().rawId());
      if (itr != eeHitRefs.end()) {
        GlobalPoint xyz = hgcGeometry_[0]->getPosition(it->id());
        heeRecVsSimX->Fill(std::get<1>(itr->second), xyz.x());
        heeRecVsSimY->Fill(std::get<2>(itr->second), xyz.y());
        heeRecVsSimZ->Fill(std::get<3>(itr->second), xyz.z());
        heedxVsX->Fill(std::get<1>(itr->second), (xyz.x() - std::get<1>(itr->second)));
        heedyVsY->Fill(std::get<2>(itr->second), (xyz.y() - std::get<2>(itr->second)));
        heedzVsZ->Fill(std::get<3>(itr->second), (xyz.z() - std::get<3>(itr->second)));
        heeEnSimRec->Fill(std::get<0>(itr->second), energy);
#ifdef EDM_ML_DEBUG
        edm::LogInfo("HGCalValid") << "EEHit: " << std::hex << it->id().rawId() << std::dec << " Sim ("
                                   << std::get<0>(itr->second) << ", " << std::get<1>(itr->second) << ", "
                                   << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ") Rec ("
                                   << energy << ", " << xyz.x() << ", " << xyz.y() << ", " << xyz.z() << ")";
#endif
      }
    }
  } else {
    edm::LogVerbatim("HGCalValid") << "No EE RecHit Found ";
  }
 
  //accessing FH Rechit information
  edm::Handle<HGChefRecHitCollection> fhRecHit;
  iEvent.getByToken(fhRecHitToken_, fhRecHit);
  if (fhRecHit.isValid()) {
    const HGChefRecHitCollection* theHits = (fhRecHit.product());
    for (auto it = theHits->begin(); it != theHits->end(); ++it) {
      double energy = it->energy();
      hefEnRec->Fill(energy);
      std::map<unsigned int, HGCHitTuple>::const_iterator itr = fhHitRefs.find(it->id().rawId());
      if (itr != fhHitRefs.end()) {
        GlobalPoint xyz = hgcGeometry_[1]->getPosition(it->id());
 
        hefRecVsSimX->Fill(std::get<1>(itr->second), xyz.x());
        hefRecVsSimY->Fill(std::get<2>(itr->second), xyz.y());
        hefRecVsSimZ->Fill(std::get<3>(itr->second), xyz.z());
        hefdxVsX->Fill(std::get<1>(itr->second), (xyz.x() - std::get<1>(itr->second)));
        hefdyVsY->Fill(std::get<2>(itr->second), (xyz.y() - std::get<2>(itr->second)));
        hefdzVsZ->Fill(std::get<3>(itr->second), (xyz.z() - std::get<3>(itr->second)));
        hefEnSimRec->Fill(std::get<0>(itr->second), energy);
#ifdef EDM_ML_DEBUG
        edm::LogInfo("HGCalValid") << "FHHit: " << std::hex << it->id().rawId() << std::dec << " Sim ("
                                   << std::get<0>(itr->second) << ", " << std::get<1>(itr->second) << ", "
                                   << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ") Rec ("
                                   << energy << "," << xyz.x() << ", " << xyz.y() << ", " << xyz.z() << ")";
#endif
      }
    }
  } else {
    edm::LogVerbatim("HGCalValid") << "No FH RecHit Found ";
  }
 
  //accessing BH Rechit information
  if (ifHCAL_) {
    edm::Handle<HBHERecHitCollection> bhRecHit;
    iEvent.getByToken(bhRecHitTokenh_, bhRecHit);
    if (bhRecHit.isValid()) {
      const HBHERecHitCollection* theHits = (bhRecHit.product());
      analyzeHGCalRecHit(theHits, bhHitRefs);
    } else {
      edm::LogVerbatim("HGCalValid") << "No BH RecHit Found ";
    }
  } else {
    edm::Handle<HGChebRecHitCollection> bhRecHit;
    iEvent.getByToken(bhRecHitTokeng_, bhRecHit);
    if (bhRecHit.isValid()) {
      const HGChebRecHitCollection* theHits = (bhRecHit.product());
      analyzeHGCalRecHit(theHits, bhHitRefs);
    } else {
      edm::LogVerbatim("HGCalValid") << "No BH RecHit Found ";
    }
  }
}

References analyzeHGCalRecHit(), analyzeHGCalSimHit(), bhRecHitTokeng_, bhRecHitTokenh_, hgcalTriggerNtuples_cfi::bhSimHits, bhSimHitToken_, HcalDDDSimConstants::cell(), TauDecayModes::dec, LEDCalibrationChannels::depth, HLT_FULL_cff::distance, eeRecHitToken_, hgcalTriggerNtuples_cfi::eeSimHits, eeSimHitToken_, HCALHighEnergyHPDFilter_cfi::energy, CastorDataFrameFilter_impl::energySum(), HcalCellType::HcalCell::eta, PVValHelper::eta, fhRecHitToken_, hgcalTriggerNtuples_cfi::fhSimHits, fhSimHitToken_, dqm::impl::MonitorElement::Fill(), spr::find(), HcalDDDRecConstants::getHCID(), HcalEndcap, hcConr_, hcCons_, hebEnSim, heedxVsX, heedyVsY, heedzVsZ, heeEnRec, heeEnSim, heeEnSimRec, heeRecVsSimX, heeRecVsSimY, heeRecVsSimZ, hefdxVsX, hefdyVsY, hefdzVsZ, hefEnRec, hefEnSim, hefEnSimRec, hefRecVsSimX, hefRecVsSimY, hefRecVsSimZ, hgcGeometry_, heavyIonCSV_trainingSettings::idx, ietaExcludeBH_, iEvent, ifHCAL_, ifHCALsim_, edm::HandleBase::isValid(), phi, HcalCellType::HcalCell::phi, edm::Handle< T >::product(), HcalCellType::HcalCell::rz, Validation_hcalonly_cfi::sign, rpcPointValidation_cfi::simHit, bscTrigger_cfi::theHits, HcalTestNumbering::unpackHcalIndex(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

analyzeHGCalRecHit()

template<class T1 >

void HGCalHitValidation::analyzeHGCalRecHit ( T1 const &  theHits, std::map< unsigned int, HGCHitTuple > const &  hitRefs  )

protected

Definition at line 440 of file HGCalHitValidation.cc.

                                                                                                               {
  for (auto it = theHits->begin(); it != theHits->end(); ++it) {
    DetId id = it->id();
    if (id.det() == DetId::Hcal and id.subdetId() == (int)(HcalEndcap)) {
      double energy = it->energy();
      hebEnRec->Fill(energy);
      GlobalPoint xyz = hcGeometry_->getGeometry(id)->getPosition();
 
      std::map<unsigned int, HGCHitTuple>::const_iterator itr = hitRefs.find(id.rawId());
      if (itr != hitRefs.end()) {
        float ang3 = xyz.phi().value();  // returns the phi in radians
        double fac = sinh(std::get<1>(itr->second));
        double pT = std::get<3>(itr->second) / fac;
        double xp = pT * cos(std::get<2>(itr->second));
        double yp = pT * sin(std::get<2>(itr->second));
        hebRecVsSimX->Fill(xp, xyz.x());
        hebRecVsSimY->Fill(yp, xyz.y());
        hebRecVsSimZ->Fill(std::get<3>(itr->second), xyz.z());
        hebdEtaVsEta->Fill(std::get<1>(itr->second), (xyz.eta() - std::get<1>(itr->second)));
        hebdPhiVsPhi->Fill(std::get<2>(itr->second), (ang3 - std::get<2>(itr->second)));
        hebdzVsZ->Fill(std::get<3>(itr->second), (xyz.z() - std::get<3>(itr->second)));
        hebEnSimRec->Fill(std::get<0>(itr->second), energy);
 
#ifdef EDM_ML_DEBUG
        edm::LogInfo("HGCalValid") << "BHHit: " << std::hex << id.rawId() << std::dec << " Sim ("
                                   << std::get<0>(itr->second) << ", " << std::get<1>(itr->second) << ", "
                                   << std::get<2>(itr->second) << ", " << std::get<3>(itr->second) << ") Rec ("
                                   << energy << ", " << xyz.x() << ", " << xyz.y() << ", " << xyz.z() << ")\n";
#endif
      }
    }
  }
}

References funct::cos(), TauDecayModes::dec, HCALHighEnergyHPDFilter_cfi::energy, PV3DBase< T, PVType, FrameType >::eta(), dqm::impl::MonitorElement::Fill(), CaloSubdetectorGeometry::getGeometry(), DetId::Hcal, HcalEndcap, hcGeometry_, hebdEtaVsEta, hebdPhiVsPhi, hebdzVsZ, hebEnRec, hebEnSimRec, hebRecVsSimX, hebRecVsSimY, hebRecVsSimZ, PV3DBase< T, PVType, FrameType >::phi(), PVValHelper::pT, funct::sin(), bscTrigger_cfi::theHits, Geom::Phi< T1, Range >::value(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

analyzeHGCalSimHit()

void HGCalHitValidation::analyzeHGCalSimHit ( edm::Handle< std::vector< PCaloHit >> const &  simHits, int  idet, MonitorElement *  hist, std::map< unsigned int, HGCHitTuple > &  hitRefs  )

protected

Definition at line 403 of file HGCalHitValidation.cc.

                                                                                        {
  const HGCalTopology& hTopo = hgcGeometry_[idet]->topology();
  for (std::vector<PCaloHit>::const_iterator simHit = simHits->begin(); simHit != simHits->end(); ++simHit) {
    int subdet, zside, layer, wafer, celltype, cell;
    HGCalTestNumbering::unpackHexagonIndex(simHit->id(), subdet, zside, layer, wafer, celltype, cell);
    std::pair<float, float> xy = hgcCons_[idet]->locateCell(cell, layer, wafer, false);
    float zp = hgcCons_[idet]->waferZ(layer, false);
    if (zside < 0)
      zp = -zp;
    float xp = (zp < 0) ? -xy.first / 10 : xy.first / 10;
    float yp = xy.second / 10.0;
 
    //skip this hit if after ganging it is not valid
    std::pair<int, int> recoLayerCell = hgcCons_[idet]->simToReco(cell, layer, wafer, hTopo.detectorType());
    cell = recoLayerCell.first;
    layer = recoLayerCell.second;
 
    //skip this hit if after ganging it is not valid
    if (layer < 0 || cell < 0) {
    } else {
      //assign the RECO DetId
      HGCalDetId id = HGCalDetId((ForwardSubdetector)(subdet), zside, layer, celltype, wafer, cell);
      float energy = simHit->energy();
 
      float energySum(energy);
      if (hitRefs.count(id.rawId()) != 0)
        energySum += std::get<0>(hitRefs[id.rawId()]);
      hitRefs[id.rawId()] = std::make_tuple(energySum, xp, yp, zp);
      hist->Fill(energy);
    }
  }
}

References HGCalTopology::detectorType(), HCALHighEnergyHPDFilter_cfi::energy, CastorDataFrameFilter_impl::energySum(), hgcCons_, hgcGeometry_, compare::hist, rpcPointValidation_cfi::simHit, FastTrackerRecHitCombiner_cfi::simHits, HGCalTestNumbering::unpackHexagonIndex(), geometryCSVtoXML::xy, and ecaldqm::zside().

Referenced by analyze().

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 149 of file HGCalHitValidation.cc.

                                                                                                {
  iB.setCurrentFolder("HGCAL/HGCalSimHitsV/HitValidation");
 
  //initiating histograms
  heedzVsZ = iB.book2D("heedzVsZ", "", 720, -360, 360, 100, -0.1, 0.1);
  heedyVsY = iB.book2D("heedyVsY", "", 400, -200, 200, 100, -0.02, 0.02);
  heedxVsX = iB.book2D("heedxVsX", "", 400, -200, 200, 100, -0.02, 0.02);
  heeRecVsSimZ = iB.book2D("heeRecVsSimZ", "", 720, -360, 360, 720, -360, 360);
  heeRecVsSimY = iB.book2D("heeRecVsSimY", "", 400, -200, 200, 400, -200, 200);
  heeRecVsSimX = iB.book2D("heeRecVsSimX", "", 400, -200, 200, 400, -200, 200);
 
  hefdzVsZ = iB.book2D("hefdzVsZ", "", 820, -410, 410, 100, -0.1, 0.1);
  hefdyVsY = iB.book2D("hefdyVsY", "", 400, -200, 200, 100, -0.02, 0.02);
  hefdxVsX = iB.book2D("hefdxVsX", "", 400, -200, 200, 100, -0.02, 0.02);
  hefRecVsSimZ = iB.book2D("hefRecVsSimZ", "", 820, -410, 410, 820, -410, 410);
  hefRecVsSimY = iB.book2D("hefRecVsSimY", "", 400, -200, 200, 400, -200, 200);
  hefRecVsSimX = iB.book2D("hefRecVsSimX", "", 400, -200, 200, 400, -200, 200);
 
  hebdzVsZ = iB.book2D("hebdzVsZ", "", 1080, -540, 540, 100, -1.0, 1.0);
  hebdPhiVsPhi = iB.book2D("hebdPhiVsPhi", "", M_PI * 100, -0.5, M_PI + 0.5, 200, -0.2, 0.2);
  hebdEtaVsEta = iB.book2D("hebdEtaVsEta", "", 1000, -5, 5, 200, -0.1, 0.1);
  hebRecVsSimZ = iB.book2D("hebRecVsSimZ", "", 1080, -540, 540, 1080, -540, 540);
  hebRecVsSimY = iB.book2D("hebRecVsSimY", "", 400, -200, 200, 400, -200, 200);
  hebRecVsSimX = iB.book2D("hebRecVsSimX", "", 400, -200, 200, 400, -200, 200);
 
  heeEnRec = iB.book1D("heeEnRec", "", 1000, 0, 10);
  heeEnSim = iB.book1D("heeEnSim", "", 1000, 0, 0.01);
  heeEnSimRec = iB.book2D("heeEnSimRec", "", 200, 0, 0.002, 200, 0, 0.2);
 
  hefEnRec = iB.book1D("hefEnRec", "", 1000, 0, 10);
  hefEnSim = iB.book1D("hefEnSim", "", 1000, 0, 0.01);
  hefEnSimRec = iB.book2D("hefEnSimRec", "", 200, 0, 0.001, 200, 0, 0.5);
 
  hebEnRec = iB.book1D("hebEnRec", "", 1000, 0, 15);
  hebEnSim = iB.book1D("hebEnSim", "", 1000, 0, 0.01);
  hebEnSimRec = iB.book2D("hebEnSimRec", "", 200, 0, 0.02, 200, 0, 4);
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), hebdEtaVsEta, hebdPhiVsPhi, hebdzVsZ, hebEnRec, hebEnSim, hebEnSimRec, hebRecVsSimX, hebRecVsSimY, hebRecVsSimZ, heedxVsX, heedyVsY, heedzVsZ, heeEnRec, heeEnSim, heeEnSimRec, heeRecVsSimX, heeRecVsSimY, heeRecVsSimZ, hefdxVsX, hefdyVsY, hefdzVsZ, hefEnRec, hefEnSim, hefEnSimRec, hefRecVsSimX, hefRecVsSimY, hefRecVsSimZ, M_PI, and dqm::implementation::NavigatorBase::setCurrentFolder().

dqmBeginRun()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 187 of file HGCalHitValidation.cc.

                                                                                    {
  //initiating hgc Geometry
  for (size_t i = 0; i < geometrySource_.size(); i++) {
    if (geometrySource_[i].find("Hcal") != std::string::npos) {
      edm::ESHandle<HcalDDDSimConstants> pHSNDC;
      iSetup.get<HcalSimNumberingRecord>().get(pHSNDC);
      if (pHSNDC.isValid()) {
        hcCons_ = pHSNDC.product();
        hgcCons_.push_back(nullptr);
      } else {
        edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDSimConstants: " << geometrySource_[i] << std::endl;
      }
      edm::ESHandle<HcalDDDRecConstants> pHRNDC;
      iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
      if (pHRNDC.isValid()) {
        hcConr_ = pHRNDC.product();
      } else {
        edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDRecConstants: " << geometrySource_[i] << std::endl;
      }
      edm::ESHandle<CaloGeometry> caloG;
      iSetup.get<CaloGeometryRecord>().get(caloG);
      if (caloG.isValid()) {
        const CaloGeometry* geo = caloG.product();
        hcGeometry_ = geo->getSubdetectorGeometry(DetId::Hcal, HcalBarrel);
        hgcGeometry_.push_back(nullptr);
      } else {
        edm::LogWarning("HGCalValid") << "Cannot initiate HcalGeometry for " << geometrySource_[i] << std::endl;
      }
    } else {
      edm::ESHandle<HGCalDDDConstants> hgcCons;
      iSetup.get<IdealGeometryRecord>().get(geometrySource_[i], hgcCons);
      if (hgcCons.isValid()) {
        hgcCons_.push_back(hgcCons.product());
      } else {
        edm::LogWarning("HGCalValid") << "Cannot initiate HGCalDDDConstants for " << geometrySource_[i] << std::endl;
      }
      edm::ESHandle<HGCalGeometry> hgcGeom;
      iSetup.get<IdealGeometryRecord>().get(geometrySource_[i], hgcGeom);
      if (hgcGeom.isValid()) {
        hgcGeometry_.push_back(hgcGeom.product());
      } else {
        edm::LogWarning("HGCalValid") << "Cannot initiate HGCalGeometry for " << geometrySource_[i] << std::endl;
      }
    }
  }
}

References spr::find(), geometrySource_, edm::EventSetup::get(), get, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, hcConr_, hcCons_, hcGeometry_, hgcCons_, hgcGeometry_, mps_fire::i, edm::ESHandleBase::isValid(), and edm::ESHandle< T >::product().

fillDescriptions()

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

static

Definition at line 141 of file HGCalHitValidation.cc.

                                                                                    {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

Member Data Documentation

bhRecHitTokeng_

edm::EDGetTokenT<HGChebRecHitCollection> HGCalHitValidation::bhRecHitTokeng_

private

Definition at line 98 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

bhRecHitTokenh_

edm::EDGetTokenT<HBHERecHitCollection> HGCalHitValidation::bhRecHitTokenh_

private

Definition at line 99 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

bhSimHitSource

edm::InputTag HGCalHitValidation::bhSimHitSource

private

Definition at line 92 of file HGCalHitValidation.cc.

bhSimHitToken_

edm::EDGetTokenT<std::vector<PCaloHit> > HGCalHitValidation::bhSimHitToken_

private

Definition at line 95 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

eeRecHitToken_

edm::EDGetTokenT<HGCeeRecHitCollection> HGCalHitValidation::eeRecHitToken_

private

Definition at line 96 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

eeSimHitSource

edm::InputTag HGCalHitValidation::eeSimHitSource

private

Definition at line 92 of file HGCalHitValidation.cc.

eeSimHitToken_

edm::EDGetTokenT<std::vector<PCaloHit> > HGCalHitValidation::eeSimHitToken_

private

Definition at line 93 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

fhRecHitToken_

edm::EDGetTokenT<HGChefRecHitCollection> HGCalHitValidation::fhRecHitToken_

private

Definition at line 97 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

fhSimHitSource

edm::InputTag HGCalHitValidation::fhSimHitSource

private

Definition at line 92 of file HGCalHitValidation.cc.

fhSimHitToken_

edm::EDGetTokenT<std::vector<PCaloHit> > HGCalHitValidation::fhSimHitToken_

private

Definition at line 94 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

geometrySource_

std::vector<std::string> HGCalHitValidation::geometrySource_

private

Definition at line 87 of file HGCalHitValidation.cc.

Referenced by dqmBeginRun(), and HGCalHitValidation().

hcConr_

const HcalDDDRecConstants* HGCalHitValidation::hcConr_

private

Definition at line 85 of file HGCalHitValidation.cc.

Referenced by analyze(), and dqmBeginRun().

hcCons_

const HcalDDDSimConstants* HGCalHitValidation::hcCons_

private

Definition at line 84 of file HGCalHitValidation.cc.

Referenced by analyze(), and dqmBeginRun().

hcGeometry_

const CaloSubdetectorGeometry* HGCalHitValidation::hcGeometry_

private

Definition at line 86 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and dqmBeginRun().

hebdEtaVsEta

MonitorElement * HGCalHitValidation::hebdEtaVsEta

private

Definition at line 104 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebdPhiVsPhi

MonitorElement * HGCalHitValidation::hebdPhiVsPhi

private

Definition at line 104 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebdzVsZ

MonitorElement* HGCalHitValidation::hebdzVsZ

private

Definition at line 104 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebEnRec

MonitorElement* HGCalHitValidation::hebEnRec

private

Definition at line 111 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebEnSim

MonitorElement * HGCalHitValidation::hebEnSim

private

Definition at line 111 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hebEnSimRec

MonitorElement * HGCalHitValidation::hebEnSimRec

private

Definition at line 109 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebRecVsSimX

MonitorElement * HGCalHitValidation::hebRecVsSimX

private

Definition at line 108 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebRecVsSimY

MonitorElement * HGCalHitValidation::hebRecVsSimY

private

Definition at line 108 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

hebRecVsSimZ

MonitorElement* HGCalHitValidation::hebRecVsSimZ

private

Definition at line 108 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalRecHit(), and bookHistograms().

heedxVsX

MonitorElement * HGCalHitValidation::heedxVsX

private

Definition at line 102 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heedyVsY

MonitorElement * HGCalHitValidation::heedyVsY

private

Definition at line 102 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heedzVsZ

MonitorElement* HGCalHitValidation::heedzVsZ

private

Definition at line 102 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeEnRec

MonitorElement* HGCalHitValidation::heeEnRec

private

Definition at line 113 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeEnSim

MonitorElement * HGCalHitValidation::heeEnSim

private

Definition at line 113 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeEnSimRec

MonitorElement* HGCalHitValidation::heeEnSimRec

private

Definition at line 109 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeRecVsSimX

MonitorElement * HGCalHitValidation::heeRecVsSimX

private

Definition at line 106 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeRecVsSimY

MonitorElement * HGCalHitValidation::heeRecVsSimY

private

Definition at line 106 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

heeRecVsSimZ

MonitorElement* HGCalHitValidation::heeRecVsSimZ

private

Definition at line 106 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefdxVsX

MonitorElement * HGCalHitValidation::hefdxVsX

private

Definition at line 103 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefdyVsY

MonitorElement * HGCalHitValidation::hefdyVsY

private

Definition at line 103 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefdzVsZ

MonitorElement* HGCalHitValidation::hefdzVsZ

private

Definition at line 103 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefEnRec

MonitorElement* HGCalHitValidation::hefEnRec

private

Definition at line 112 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefEnSim

MonitorElement * HGCalHitValidation::hefEnSim

private

Definition at line 112 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefEnSimRec

MonitorElement * HGCalHitValidation::hefEnSimRec

private

Definition at line 109 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefRecVsSimX

MonitorElement * HGCalHitValidation::hefRecVsSimX

private

Definition at line 107 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefRecVsSimY

MonitorElement * HGCalHitValidation::hefRecVsSimY

private

Definition at line 107 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hefRecVsSimZ

MonitorElement* HGCalHitValidation::hefRecVsSimZ

private

Definition at line 107 of file HGCalHitValidation.cc.

Referenced by analyze(), and bookHistograms().

hgcCons_

std::vector<const HGCalDDDConstants*> HGCalHitValidation::hgcCons_

private

Definition at line 82 of file HGCalHitValidation.cc.

Referenced by analyzeHGCalSimHit(), and dqmBeginRun().

hgcGeometry_

std::vector<const HGCalGeometry*> HGCalHitValidation::hgcGeometry_

private

Definition at line 83 of file HGCalHitValidation.cc.

Referenced by analyze(), analyzeHGCalSimHit(), and dqmBeginRun().

ietaExcludeBH_

std::vector<int> HGCalHitValidation::ietaExcludeBH_

private

Definition at line 88 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

ifHCAL_

bool HGCalHitValidation::ifHCAL_

private

Definition at line 89 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().

ifHCALsim_

bool HGCalHitValidation::ifHCALsim_

private

Definition at line 90 of file HGCalHitValidation.cc.

Referenced by analyze(), and HGCalHitValidation().