HGCalSimHitValidation Class Reference

Inheritance diagram for HGCalSimHitValidation:

Classes
struct	energysum
struct	hitsinfo

Public Member Functions
	HGCalSimHitValidation (const edm::ParameterSet &)
	~HGCalSimHitValidation () 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 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
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	analyzeHits (std::vector< PCaloHit > &hits)
bool	defineGeometry (edm::ESTransientHandle< cms::DDCompactView > &ddViewH)
bool	defineGeometry (edm::ESTransientHandle< DDCompactView > &ddViewH)
void	fillHitsInfo (std::pair< hitsinfo, energysum > hit_, unsigned int itimeslice, double esum)
void	fillOccupancyMap (std::map< int, int > &OccupancyMap, int layer)

Private Attributes
std::string	caloHitSource_
std::vector< MonitorElement * >	energy_ [maxTime_]
std::vector< MonitorElement * >	EtaPhi_Minus_
std::vector< MonitorElement * >	EtaPhi_Plus_
int	firstLayer_
bool	fromDDD_
const HcalDDDRecConstants *	hcons_
bool	heRebuild_
const HGCalDDDConstants *	hgcons_
std::vector< MonitorElement * >	HitOccupancy_Minus_
std::vector< MonitorElement * >	HitOccupancy_Plus_
unsigned int	layers_
MonitorElement *	MeanHitOccupancy_Minus_
MonitorElement *	MeanHitOccupancy_Plus_
std::string	nameDetector_
unsigned int	nTimes_
bool	symmDet_
bool	testNumber_
std::vector< double >	times_
edm::EDGetTokenT< edm::HepMCProduct >	tok_hepMC_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_
std::map< uint32_t, HepGeom::Transform3D >	transMap_
int	verbosity_

Static Private Attributes
static const unsigned int	maxTime_ = 6

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

Definition at line 55 of file HGCalSimHitValidation.cc.

Constructor & Destructor Documentation

HGCalSimHitValidation()

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

explicit

Definition at line 113 of file HGCalSimHitValidation.cc.

    : nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
      caloHitSource_(iConfig.getParameter<std::string>("CaloHitSource")),
      times_(iConfig.getParameter<std::vector<double> >("TimeSlices")),
      verbosity_(iConfig.getUntrackedParameter<int>("Verbosity", 0)),
      testNumber_(iConfig.getUntrackedParameter<bool>("TestNumber", true)),
      symmDet_(true),
      fromDDD_(iConfig.getUntrackedParameter<bool>("fromDDD", true)),
      firstLayer_(1) {
  heRebuild_ = (nameDetector_ == "HCal") ? true : false;
  tok_hepMC_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
  tok_hits_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", caloHitSource_));
  nTimes_ = (times_.size() > maxTime_) ? maxTime_ : times_.size();
}

References caloHitSource_, heRebuild_, maxTime_, nameDetector_, nTimes_, times_, tok_hepMC_, and tok_hits_.

~HGCalSimHitValidation()

HGCalSimHitValidation::~HGCalSimHitValidation ( )

inlineoverride

Definition at line 76 of file HGCalSimHitValidation.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 140 of file HGCalSimHitValidation.cc.

                                                                                       {
  //Generator input
  if (verbosity_ > 0) {
    edm::Handle<edm::HepMCProduct> evtMC;
    iEvent.getByToken(tok_hepMC_, evtMC);
    if (!evtMC.isValid()) {
      edm::LogVerbatim("HGCalValidation") << "no HepMCProduct found";
    } else {
      const HepMC::GenEvent* myGenEvent = evtMC->GetEvent();
      unsigned int k(0);
      for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
           ++p, ++k) {
        edm::LogVerbatim("HGCalValidation") << "Particle[" << k << "] with pt " << (*p)->momentum().perp() << " eta "
                                            << (*p)->momentum().eta() << " phi " << (*p)->momentum().phi();
      }
    }
  }
 
  //Now the hits
  edm::Handle<edm::PCaloHitContainer> theCaloHitContainers;
  iEvent.getByToken(tok_hits_, theCaloHitContainers);
  if (theCaloHitContainers.isValid()) {
    if (verbosity_ > 0)
      edm::LogVerbatim("HGCalValidation") << " PcalohitItr = " << theCaloHitContainers->size();
    std::vector<PCaloHit> caloHits;
    caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
    if (heRebuild_ && testNumber_) {
      for (unsigned int i = 0; i < caloHits.size(); ++i) {
        unsigned int id_ = caloHits[i].id();
        HcalDetId hid = HcalHitRelabeller::relabel(id_, hcons_);
        if (hid.subdet() != int(HcalEndcap))
          hid = HcalDetId(HcalEmpty, hid.ieta(), hid.iphi(), hid.depth());
        caloHits[i].setID(hid.rawId());
        if (verbosity_ > 0)
          edm::LogVerbatim("HGCalValidation") << "Hit[" << i << "] " << hid;
      }
    }
    analyzeHits(caloHits);
  } else if (verbosity_ > 0) {
    edm::LogVerbatim("HGCalValidation") << "PCaloHitContainer does not exist!";
  }
}

References analyzeHits(), HcalDetId::depth(), edm::HepMCProduct::GetEvent(), HcalEmpty, HcalEndcap, hcons_, heRebuild_, mps_fire::i, edm::HandleBase::id(), HcalDetId::ieta(), iEvent, createfilelist::int, HcalDetId::iphi(), edm::HandleBase::isValid(), dqmdumpme::k, AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), HcalHitRelabeller::relabel(), HcalDetId::subdet(), testNumber_, tok_hepMC_, tok_hits_, and verbosity_.

analyzeHits()

void HGCalSimHitValidation::analyzeHits ( std::vector< PCaloHit > &  hits )

private

Definition at line 183 of file HGCalSimHitValidation.cc.

                                                                 {
  std::map<int, int> OccupancyMap_plus, OccupancyMap_minus;
  OccupancyMap_plus.clear();
  OccupancyMap_minus.clear();
 
  std::map<uint32_t, std::pair<hitsinfo, energysum> > map_hits;
  map_hits.clear();
 
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " << hits.size() << " PcaloHit elements";
  unsigned int nused(0);
  for (unsigned int i = 0; i < hits.size(); i++) {
    double energy = hits[i].energy();
    double time = hits[i].time();
    uint32_t id_ = hits[i].id();
    int cell, sector, subsector(0), layer, zside;
    int subdet(0), cell2(0), type(0);
    if (heRebuild_) {
      HcalDetId detId = HcalDetId(id_);
      subdet = detId.subdet();
      if (subdet != static_cast<int>(HcalEndcap))
        continue;
      cell = detId.ietaAbs();
      sector = detId.iphi();
      subsector = 1;
      layer = detId.depth();
      zside = detId.zside();
    } else if (hgcons_->waferHexagon8()) {
      HGCSiliconDetId detId = HGCSiliconDetId(id_);
      subdet = ForwardEmpty;
      cell = detId.cellU();
      cell2 = detId.cellV();
      sector = detId.waferU();
      subsector = detId.waferV();
      type = detId.type();
      layer = detId.layer();
      zside = detId.zside();
    } else if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
      HGCalTestNumbering::unpackSquareIndex(id_, zside, layer, sector, subsector, cell);
    } else if (hgcons_->tileTrapezoid()) {
      HGCScintillatorDetId detId = HGCScintillatorDetId(id_);
      subdet = ForwardEmpty;
      sector = detId.ietaAbs();
      cell = detId.iphi();
      subsector = 1;
      type = detId.type();
      layer = detId.layer();
      zside = detId.zside();
    } else {
      HGCalTestNumbering::unpackHexagonIndex(id_, subdet, zside, layer, sector, type, cell);
    }
    nused++;
    if (verbosity_ > 1)
      edm::LogVerbatim("HGCalValidation")
          << "Detector " << nameDetector_ << " zside = " << zside << " sector|wafer = " << sector << ":" << subsector
          << " type = " << type << " layer = " << layer << " cell = " << cell << ":" << cell2 << " energy = " << energy
          << " energyem = " << hits[i].energyEM() << " energyhad = " << hits[i].energyHad() << " time = " << time;
 
    HepGeom::Point3D<float> gcoord;
    if (heRebuild_) {
      std::pair<double, double> etaphi = hcons_->getEtaPhi(subdet, zside * cell, sector);
      double rz = hcons_->getRZ(subdet, zside * cell, layer);
      if (verbosity_ > 2)
        edm::LogVerbatim("HGCalValidation") << "i/p " << subdet << ":" << zside << ":" << cell << ":" << sector << ":"
                                            << layer << " o/p " << etaphi.first << ":" << etaphi.second << ":" << rz;
      gcoord = HepGeom::Point3D<float>(rz * cos(etaphi.second) / cosh(etaphi.first),
                                       rz * sin(etaphi.second) / cosh(etaphi.first),
                                       rz * tanh(etaphi.first));
    } else if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
      std::pair<float, float> xy = hgcons_->locateCell(cell, layer, subsector, false);
      const HepGeom::Point3D<float> lcoord(xy.first, xy.second, 0);
      int subs = (symmDet_ ? 0 : subsector);
      id_ = HGCalTestNumbering::packSquareIndex(zside, layer, sector, subs, 0);
      gcoord = (transMap_[id_] * lcoord);
    } else {
      std::pair<float, float> xy;
      if (hgcons_->waferHexagon8()) {
        xy = hgcons_->locateCell(layer, sector, subsector, cell, cell2, false, true);
      } else if (hgcons_->tileTrapezoid()) {
        xy = hgcons_->locateCellTrap(layer, sector, cell, false);
      } else {
        xy = hgcons_->locateCell(cell, layer, sector, false);
      }
      double zp = hgcons_->waferZ(layer, false);
      if (zside < 0)
        zp = -zp;
      float xp = (zp < 0) ? -xy.first : xy.first;
      gcoord = HepGeom::Point3D<float>(xp, xy.second, zp);
    }
    double tof = (gcoord.mag() * CLHEP::mm) / CLHEP::c_light;
    if (verbosity_ > 1)
      edm::LogVerbatim("HGCalValidation")
          << std::hex << id_ << std::dec << " global coordinate " << gcoord << " time " << time << ":" << tof;
    time -= tof;
 
    energysum esum;
    hitsinfo hinfo;
    if (map_hits.count(id_) != 0) {
      hinfo = map_hits[id_].first;
      esum = map_hits[id_].second;
    } else {
      hinfo.x = gcoord.x();
      hinfo.y = gcoord.y();
      hinfo.z = gcoord.z();
      hinfo.sector = sector;
      hinfo.sector2 = subsector;
      hinfo.cell = cell;
      hinfo.cell2 = cell2;
      hinfo.type = type;
      hinfo.layer = layer - firstLayer_;
      hinfo.phi = gcoord.getPhi();
      hinfo.eta = gcoord.getEta();
    }
    esum.etotal += energy;
    for (unsigned int k = 0; k < nTimes_; ++k) {
      if (time > 0 && time < times_[k])
        esum.eTime[k] += energy;
    }
 
    if (verbosity_ > 1)
      edm::LogVerbatim("HGCalValidation") << " -----------------------   gx = " << hinfo.x << " gy = " << hinfo.y
                                          << " gz = " << hinfo.z << " phi = " << hinfo.phi << " eta = " << hinfo.eta;
    map_hits[id_] = std::pair<hitsinfo, energysum>(hinfo, esum);
  }
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " << map_hits.size()
                                        << " detector elements being hit";
 
  std::map<uint32_t, std::pair<hitsinfo, energysum> >::iterator itr;
  for (itr = map_hits.begin(); itr != map_hits.end(); ++itr) {
    hitsinfo hinfo = (*itr).second.first;
    energysum esum = (*itr).second.second;
    int layer = hinfo.layer;
    double eta = hinfo.eta;
 
    for (unsigned int itimeslice = 0; itimeslice < nTimes_; itimeslice++) {
      fillHitsInfo((*itr).second, itimeslice, esum.eTime[itimeslice]);
    }
 
    if (eta > 0.0)
      fillOccupancyMap(OccupancyMap_plus, layer);
    else
      fillOccupancyMap(OccupancyMap_minus, layer);
  }
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << "With map:used:total " << hits.size() << "|" << nused << "|"
                                        << map_hits.size() << " hits";
 
  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 hgcalTopologyTester_cfi::cell2, HGCSiliconDetId::cellU(), HGCSiliconDetId::cellV(), funct::cos(), TauDecayModes::dec, HcalDetId::depth(), HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, fastmath::etaphi(), HGCalSimHitValidation::energysum::eTime, HGCalSimHitValidation::energysum::etotal, fillHitsInfo(), fillOccupancyMap(), firstLayer_, ForwardEmpty, HGCalDDDConstants::geomMode(), HcalDDDRecConstants::getEtaPhi(), HcalDDDRecConstants::getRZ(), HcalEndcap, hcons_, heRebuild_, hgcons_, HitOccupancy_Minus_, HitOccupancy_Plus_, hfClusterShapes_cfi::hits, mps_fire::i, HGCScintillatorDetId::ietaAbs(), HcalDetId::ietaAbs(), HGCScintillatorDetId::iphi(), HcalDetId::iphi(), dqmdumpme::k, HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellTrap(), nameDetector_, nTimes_, HGCalTestNumbering::packSquareIndex(), funct::sin(), HGCalGeometryMode::Square, HcalDetId::subdet(), symmDet_, HGCalDDDConstants::tileTrapezoid(), ntuplemaker::time, times_, transMap_, HGCScintillatorDetId::type(), HGCSiliconDetId::type(), HGCalTestNumbering::unpackHexagonIndex(), HGCalTestNumbering::unpackSquareIndex(), verbosity_, HGCalDDDConstants::waferHexagon8(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferV(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xy, HGCScintillatorDetId::zside(), HGCSiliconDetId::zside(), HcalDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 486 of file HGCalSimHitValidation.cc.

                                                                                                   {
  iB.setCurrentFolder("HGCAL/HGCalSimHitsV/" + nameDetector_);
 
  std::ostringstream histoname;
  for (unsigned int il = 0; il < layers_; ++il) {
    int ilayer = firstLayer_ + static_cast<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(), "HitOccupancy_Plus", 501, -0.5, 500.5));
    histoname.str("");
    histoname << "HitOccupancy_Minus_layer_" << istr1;
    HitOccupancy_Minus_.push_back(iB.book1D(histoname.str().c_str(), "HitOccupancy_Minus", 501, -0.5, 500.5));
 
    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));
 
    for (unsigned int itimeslice = 0; itimeslice < nTimes_; itimeslice++) {
      histoname.str("");
      histoname << "energy_time_" << itimeslice << "_layer_" << istr1;
      energy_[itimeslice].push_back(iB.book1D(histoname.str().c_str(), "energy_", 100, 0, 0.1));
    }
  }
 
  MeanHitOccupancy_Plus_ = iB.book1D("MeanHitOccupancy_Plus", "MeanHitOccupancy_Plus", layers_, 0.5, layers_ + 0.5);
  MeanHitOccupancy_Minus_ = iB.book1D("MeanHitOccupancy_Minus", "MeanHitOccupancy_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_, layers_, MeanHitOccupancy_Minus_, MeanHitOccupancy_Plus_, nameDetector_, nTimes_, pi, and dqm::implementation::NavigatorBase::setCurrentFolder().

defineGeometry() [1/2]

bool HGCalSimHitValidation::defineGeometry ( edm::ESTransientHandle< cms::DDCompactView > &  ddViewH )

private

Definition at line 416 of file HGCalSimHitValidation.cc.

                                                                                          {
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << "Initialize HGCalDDDConstants (DD4hep) for " << nameDetector_ << " : "
                                        << hgcons_;
 
  if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
    const cms::DDCompactView& cview = *ddViewH;
    const cms::DDFilter filter("Volume", nameDetector_);
    cms::DDFilteredView fv(cview, filter);
 
    while (fv.firstChild()) {
      const auto& name = fv.name();
      int isd = (name.find(nameDetector_) == std::string::npos) ? -1 : 1;
      if (isd > 0) {
        const auto& copy = fv.copyNos();
        int nsiz = static_cast<int>(copy.size());
        int lay = (nsiz > 0) ? copy[0] : -1;
        int sec = (nsiz > 1) ? copy[1] : -1;
        int zp = (nsiz > 3) ? copy[3] : -1;
        if (zp != 1)
          zp = -1;
        const auto& pars = fv.parameters();
        int subs = (pars[6] > 0 ? 1 : 0);
        symmDet_ = (pars[6] == 0 ? true : false);
        uint32_t id = HGCalTestNumbering::packSquareIndex(zp, lay, sec, subs, 0);
        DD3Vector x, y, z;
        fv.rotation().GetComponents(x, y, z);
        const CLHEP::HepRep3x3 rotation(x.X(), y.X(), z.X(), x.Y(), y.Y(), z.Y(), x.Z(), y.Z(), z.Z());
        const CLHEP::HepRotation hr(rotation);
        const CLHEP::Hep3Vector h3v(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
        const HepGeom::Transform3D ht3d(hr, h3v);
        transMap_.insert(std::make_pair(id, ht3d));
        if (verbosity_ > 2)
          edm::LogVerbatim("HGCalValidation") << HGCalDetId(id) << " Transform using " << h3v << " and " << hr;
      }
    }
    if (verbosity_ > 0)
      edm::LogVerbatim("HGCalValidation") << "Finds " << transMap_.size() << " elements and SymmDet_ = " << symmDet_;
  }
  return true;
}

References filterCSVwithJSON::copy, cms::DDFilteredView::copyNos(), ALCARECOTkAlBeamHalo_cff::filter, cms::DDFilteredView::firstChild(), HGCalDDDConstants::geomMode(), hgcons_, Skims_PA_cff::name, cms::DDFilteredView::name(), nameDetector_, HGCalTestNumbering::packSquareIndex(), cms::DDFilteredView::parameters(), idealTransformation::rotation, cms::DDFilteredView::rotation(), fileinputsource_cfi::sec, HGCalGeometryMode::Square, symmDet_, cms::DDFilteredView::translation(), transMap_, funct::true, verbosity_, x, y, and z.

defineGeometry() [2/2]

bool HGCalSimHitValidation::defineGeometry ( edm::ESTransientHandle< DDCompactView > &  ddViewH )

private

Definition at line 368 of file HGCalSimHitValidation.cc.

                                                                                       {
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << "Initialize HGCalDDDConstants (DDD) for " << nameDetector_ << " : "
                                        << hgcons_;
 
  if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
    const DDCompactView& cview = *ddViewH;
    std::string attribute = "Volume";
    std::string value = nameDetector_;
 
    DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0)};
    DDFilteredView fv(cview, filter);
    bool dodet = fv.firstChild();
 
    while (dodet) {
      const DDSolid& sol = fv.logicalPart().solid();
      const std::string& name = sol.name().fullname();
      int isd = (name.find(nameDetector_) == std::string::npos) ? -1 : 1;
      if (isd > 0) {
        std::vector<int> copy = fv.copyNumbers();
        int nsiz = static_cast<int>(copy.size());
        int lay = (nsiz > 0) ? copy[nsiz - 1] : -1;
        int sec = (nsiz > 1) ? copy[nsiz - 2] : -1;
        int zp = (nsiz > 3) ? copy[nsiz - 4] : -1;
        if (zp != 1)
          zp = -1;
        const DDTrap& trp = static_cast<DDTrap>(sol);
        int subs = (trp.alpha1() > 0 ? 1 : 0);
        symmDet_ = (trp.alpha1() == 0 ? true : false);
        uint32_t id = HGCalTestNumbering::packSquareIndex(zp, lay, sec, subs, 0);
        DD3Vector x, y, z;
        fv.rotation().GetComponents(x, y, z);
        const CLHEP::HepRep3x3 rotation(x.X(), y.X(), z.X(), x.Y(), y.Y(), z.Y(), x.Z(), y.Z(), z.Z());
        const CLHEP::HepRotation hr(rotation);
        const CLHEP::Hep3Vector h3v(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
        const HepGeom::Transform3D ht3d(hr, h3v);
        transMap_.insert(std::make_pair(id, ht3d));
        if (verbosity_ > 2)
          edm::LogVerbatim("HGCalValidation") << HGCalDetId(id) << " Transform using " << h3v << " and " << hr;
      }
      dodet = fv.next();
    }
    if (verbosity_ > 0)
      edm::LogVerbatim("HGCalValidation") << "Finds " << transMap_.size() << " elements and SymmDet_ = " << symmDet_;
  }
  return true;
}

References DDTrap::alpha1(), filterCSVwithJSON::copy, DDFilteredView::copyNumbers(), ALCARECOTkAlBeamHalo_cff::filter, DDFilteredView::firstChild(), DDName::fullname(), HGCalDDDConstants::geomMode(), hgcons_, DDFilteredView::logicalPart(), Skims_PA_cff::name, DDBase< N, C >::name(), nameDetector_, DDFilteredView::next(), HGCalTestNumbering::packSquareIndex(), idealTransformation::rotation, DDFilteredView::rotation(), fileinputsource_cfi::sec, DDLogicalPart::solid(), HGCalGeometryMode::Square, AlCaHLTBitMon_QueryRunRegistry::string, symmDet_, DDFilteredView::translation(), transMap_, funct::true, verbosity_, x, y, and z.

Referenced by dqmBeginRun().

dqmBeginRun()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 459 of file HGCalSimHitValidation.cc.

                                                                                  {
  if (heRebuild_) {
    edm::ESHandle<HcalDDDRecConstants> pHRNDC;
    iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
    hcons_ = &(*pHRNDC);
    layers_ = hcons_->getMaxDepth(1);
  } else {
    edm::ESHandle<HGCalDDDConstants> pHGDC;
    iSetup.get<IdealGeometryRecord>().get(nameDetector_, pHGDC);
    hgcons_ = &(*pHGDC);
    layers_ = hgcons_->layers(false);
    firstLayer_ = hgcons_->firstLayer();
    if (fromDDD_) {
      edm::ESTransientHandle<DDCompactView> pDD;
      iSetup.get<IdealGeometryRecord>().get(pDD);
      defineGeometry(pDD);
    } else {
      edm::ESTransientHandle<cms::DDCompactView> pDD;
      iSetup.get<IdealGeometryRecord>().get(pDD);
      defineGeometry(pDD);
    }
  }
  if (verbosity_ > 0)
    edm::LogVerbatim("HGCalValidation") << nameDetector_ << " defined with " << layers_ << " Layers with first at "
                                        << firstLayer_;
}

References defineGeometry(), HGCalDDDConstants::firstLayer(), firstLayer_, fromDDD_, edm::EventSetup::get(), get, HcalDDDRecConstants::getMaxDepth(), hcons_, heRebuild_, hgcons_, HGCalDDDConstants::layers(), layers_, nameDetector_, and verbosity_.

fillDescriptions()

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

static

Definition at line 128 of file HGCalSimHitValidation.cc.

                                                                                       {
  edm::ParameterSetDescription desc;
  std::vector<double> times = {25.0, 1000.0};
  desc.add<std::string>("DetectorName", "HGCalEESensitive");
  desc.add<std::string>("CaloHitSource", "HGCHitsEE");
  desc.add<std::vector<double> >("TimeSlices", times);
  desc.addUntracked<int>("Verbosity", 0);
  desc.addUntracked<bool>("TestNumber", true);
  desc.addUntracked<bool>("fromDDD", true);
  descriptions.add("hgcalSimHitValidationEE", desc);
}

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

fillHitsInfo()

void HGCalSimHitValidation::fillHitsInfo ( std::pair< hitsinfo, energysum >  hit_, unsigned int  itimeslice, double  esum  )

private

Definition at line 351 of file HGCalSimHitValidation.cc.

                                                                                                                {
  unsigned int ilayer = hits.first.layer;
  if (ilayer < layers_) {
    energy_[itimeslice].at(ilayer)->Fill(esum);
    if (itimeslice == 0) {
      EtaPhi_Plus_.at(ilayer)->Fill(hits.first.eta, hits.first.phi);
      EtaPhi_Minus_.at(ilayer)->Fill(hits.first.eta, hits.first.phi);
    }
  } else {
    if (verbosity_ > 0)
      edm::LogVerbatim("HGCalValidation")
          << "Problematic Hit for " << nameDetector_ << " at sector " << hits.first.sector << ":" << hits.first.sector2
          << " layer " << hits.first.layer << " cell " << hits.first.cell << ":" << hits.first.cell2 << " energy "
          << hits.second.etotal;
  }
}

References energy_, EtaPhi_Minus_, EtaPhi_Plus_, hfClusterShapes_cfi::hits, layers_, nameDetector_, and verbosity_.

Referenced by analyzeHits().

fillOccupancyMap()

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

private

Definition at line 343 of file HGCalSimHitValidation.cc.

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

Referenced by analyzeHits().

Member Data Documentation

caloHitSource_

std::string HGCalSimHitValidation::caloHitSource_

private

Definition at line 93 of file HGCalSimHitValidation.cc.

Referenced by HGCalSimHitValidation().

energy_

std::vector<MonitorElement*> HGCalSimHitValidation::energy_[maxTime_]

private

Definition at line 109 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

EtaPhi_Minus_

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

private

Definition at line 106 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

EtaPhi_Plus_

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

private

Definition at line 106 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

firstLayer_

int HGCalSimHitValidation::firstLayer_

private

Definition at line 102 of file HGCalSimHitValidation.cc.

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

fromDDD_

bool HGCalSimHitValidation::fromDDD_

private

Definition at line 98 of file HGCalSimHitValidation.cc.

Referenced by dqmBeginRun().

hcons_

const HcalDDDRecConstants* HGCalSimHitValidation::hcons_

private

Definition at line 95 of file HGCalSimHitValidation.cc.

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

heRebuild_

bool HGCalSimHitValidation::heRebuild_

private

Definition at line 98 of file HGCalSimHitValidation.cc.

Referenced by analyze(), analyzeHits(), dqmBeginRun(), and HGCalSimHitValidation().

hgcons_

const HGCalDDDConstants* HGCalSimHitValidation::hgcons_

private

Definition at line 94 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), defineGeometry(), and dqmBeginRun().

HitOccupancy_Minus_

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

private

Definition at line 105 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and bookHistograms().

HitOccupancy_Plus_

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

private

Definition at line 105 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and bookHistograms().

layers_

unsigned int HGCalSimHitValidation::layers_

private

Definition at line 101 of file HGCalSimHitValidation.cc.

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

maxTime_

const unsigned int HGCalSimHitValidation::maxTime_ = 6

staticprivate

Definition at line 108 of file HGCalSimHitValidation.cc.

Referenced by HGCalSimHitValidation().

MeanHitOccupancy_Minus_

MonitorElement * HGCalSimHitValidation::MeanHitOccupancy_Minus_

private

Definition at line 107 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms().

MeanHitOccupancy_Plus_

MonitorElement* HGCalSimHitValidation::MeanHitOccupancy_Plus_

private

Definition at line 107 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms().

nameDetector_

std::string HGCalSimHitValidation::nameDetector_

private

Definition at line 93 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), bookHistograms(), defineGeometry(), dqmBeginRun(), fillHitsInfo(), and HGCalSimHitValidation().

nTimes_

unsigned int HGCalSimHitValidation::nTimes_

private

Definition at line 110 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), bookHistograms(), and HGCalSimHitValidation().

symmDet_

bool HGCalSimHitValidation::symmDet_

private

Definition at line 98 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and defineGeometry().

testNumber_

bool HGCalSimHitValidation::testNumber_

private

Definition at line 98 of file HGCalSimHitValidation.cc.

Referenced by analyze().

times_

std::vector<double> HGCalSimHitValidation::times_

private

Definition at line 96 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and HGCalSimHitValidation().

tok_hepMC_

edm::EDGetTokenT<edm::HepMCProduct> HGCalSimHitValidation::tok_hepMC_

private

Definition at line 100 of file HGCalSimHitValidation.cc.

Referenced by analyze(), and HGCalSimHitValidation().

tok_hits_

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalSimHitValidation::tok_hits_

private

Definition at line 99 of file HGCalSimHitValidation.cc.

Referenced by analyze(), and HGCalSimHitValidation().

transMap_

std::map<uint32_t, HepGeom::Transform3D> HGCalSimHitValidation::transMap_

private

Definition at line 103 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and defineGeometry().

verbosity_

int HGCalSimHitValidation::verbosity_

private

Definition at line 97 of file HGCalSimHitValidation.cc.

Referenced by analyze(), analyzeHits(), defineGeometry(), dqmBeginRun(), and fillHitsInfo().