CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
HGCDigitizer Class Reference

#include <HGCDigitizer.h>

Public Types

typedef std::tuple< int, uint32_t, float > HGCCaloHitTuple_t
 

Public Member Functions

void accumulate (edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 handle SimHit accumulation More...
 
void accumulate (edm::Handle< edm::PCaloHitContainer > const &hits, int bxCrossing, const HGCalGeometry *geom, CLHEP::HepRandomEngine *hre)
 
void accumulate (PileUpEventPrincipal const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 
void accumulate_forPreMix (const PHGCSimAccumulator &simAccumulator, const bool minbiasFlag)
 
void accumulate_forPreMix (edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 
void accumulate_forPreMix (edm::Handle< edm::PCaloHitContainer > const &hits, int bxCrossing, const HGCalGeometry *geom, CLHEP::HepRandomEngine *hre)
 
void accumulate_forPreMix (PileUpEventPrincipal const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 
std::string digiCollection ()
 
void finalizeEvent (edm::Event &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 
 HGCDigitizer (const edm::ParameterSet &ps, edm::ConsumesCollector &iC)
 
void initializeEvent (edm::Event const &e, edm::EventSetup const &c)
 actions at the start/end of event More...
 
 ~HGCDigitizer ()
 

Static Public Member Functions

static bool orderByDetIdThenTime (const HGCCaloHitTuple_t &a, const HGCCaloHitTuple_t &b)
 

Private Member Functions

void checkPosition (const HGCalDigiCollection *digis) const
 
uint32_t getType () const
 
void resetPUSimHitDataAccumulator ()
 
void resetSimHitDataAccumulator ()
 

Private Attributes

std::array< double, 4 > averageOccupancies_
 
double bxTime_
 
std::vector< float > cce_
 
std::string digiCollection_
 
int digitizationType_
 
double ev_per_eh_pair_
 
const edm::ESGetToken< CaloGeometry, CaloGeometryRecordgeomToken_
 
edm::ESWatcher< CaloGeometryRecordgeomWatcher_
 
const HGCalGeometrygHGCal_ = nullptr
 
std::string hitCollection_
 
std::unordered_map< uint32_t, bool > hitOrder_monitor
 
std::unordered_map< uint32_t, std::vector< std::pair< float, float > > > hitRefs_bx0
 
int maxSimHitsAccTime_
 
uint32_t nEvents_
 
std::unordered_map< uint32_t, std::vector< std::tuple< float, float, float > > > PhitRefs_bx0
 
bool premixStage1_
 
double premixStage1MaxCharge_
 
double premixStage1MinCharge_
 
std::unique_ptr< hgc::HGCPUSimHitDataAccumulatorpusimHitAccumulator_
 
float refSpeed_
 
std::unique_ptr< hgc::HGCSimHitDataAccumulatorsimHitAccumulator_
 
std::unique_ptr< HGCDigitizerBasetheDigitizer_
 
float tofDelay_
 
std::unordered_set< DetIdvalidIds_
 
uint32_t verbosity_
 

Static Private Attributes

static const unsigned int maxBx_ = 14
 
static const unsigned int thisBx_ = 9
 

Detailed Description

Definition at line 29 of file HGCDigitizer.h.

Member Typedef Documentation

◆ HGCCaloHitTuple_t

typedef std::tuple<int, uint32_t, float> HGCDigitizer::HGCCaloHitTuple_t

Definition at line 35 of file HGCDigitizer.h.

Constructor & Destructor Documentation

◆ HGCDigitizer()

HGCDigitizer::HGCDigitizer ( const edm::ParameterSet ps,
edm::ConsumesCollector iC 
)

Definition at line 231 of file HGCDigitizer.cc.

234  geomToken_(iC.esConsumes()),
235  refSpeed_(0.1 * CLHEP::c_light), //[CLHEP::c_light]=mm/ns convert to cm/ns
236  averageOccupancies_(occupancyGuesses),
237  nEvents_(1) {
238  //configure from cfg
239 
240  hitCollection_ = ps.getParameter<std::string>("hitCollection");
241  digiCollection_ = ps.getParameter<std::string>("digiCollection");
242  maxSimHitsAccTime_ = ps.getParameter<uint32_t>("maxSimHitsAccTime");
243  bxTime_ = ps.getParameter<double>("bxTime");
244  digitizationType_ = ps.getParameter<uint32_t>("digitizationType");
245  verbosity_ = ps.getUntrackedParameter<uint32_t>("verbosity", 0);
246  tofDelay_ = ps.getParameter<double>("tofDelay");
247  premixStage1_ = ps.getParameter<bool>("premixStage1");
248  premixStage1MinCharge_ = ps.getParameter<double>("premixStage1MinCharge");
249  premixStage1MaxCharge_ = ps.getParameter<double>("premixStage1MaxCharge");
250  iC.consumes<std::vector<PCaloHit>>(edm::InputTag("g4SimHits", hitCollection_));
251  const auto& myCfg_ = ps.getParameter<edm::ParameterSet>("digiCfg");
252 
253  if (myCfg_.existsAs<edm::ParameterSet>("chargeCollectionEfficiencies")) {
254  cce_.clear();
255  const auto& temp = myCfg_.getParameter<edm::ParameterSet>("chargeCollectionEfficiencies")
256  .getParameter<std::vector<double>>("values");
257  for (double cce : temp) {
258  cce_.emplace_back(cce);
259  }
260  } else {
261  std::vector<float>().swap(cce_);
262  }
263 
264  auto const& pluginName = ps.getParameter<std::string>("digitizer");
265  theDigitizer_ = HGCDigitizerPluginFactory::get()->create(pluginName, ps);
266 }

References bxTime_, cce_, edm::ConsumesCollector::consumes(), digiCollection_, digitizationType_, get, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hitCollection_, HLT_FULL_cff::InputTag, maxSimHitsAccTime_, premixStage1_, premixStage1MaxCharge_, premixStage1MinCharge_, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, theDigitizer_, tofDelay_, and verbosity_.

◆ ~HGCDigitizer()

HGCDigitizer::~HGCDigitizer ( )
inline

Definition at line 32 of file HGCDigitizer.h.

32 {}

Member Function Documentation

◆ accumulate() [1/3]

void HGCDigitizer::accumulate ( edm::Event const &  e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

handle SimHit accumulation

Definition at line 365 of file HGCDigitizer.cc.

365  {
366  //get inputs
368  e.getByLabel(edm::InputTag("g4SimHits", hitCollection_), hits);
369  if (!hits.isValid()) {
370  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
371  return;
372  }
373 
374  //accumulate in-time the main event
375  if (nullptr != gHGCal_) {
376  accumulate(hits, 0, gHGCal_, hre);
377  } else {
378  throw cms::Exception("BadConfiguration") << "HGCDigitizer is not producing EE, FH, or BH digis!";
379  }
380 }

References MillePedeFileConverter_cfg::e, Exception, gHGCal_, hitCollection_, and hfClusterShapes_cfi::hits.

Referenced by HGCDigiProducer::accumulate(), and accumulate().

◆ accumulate() [2/3]

void HGCDigitizer::accumulate ( edm::Handle< edm::PCaloHitContainer > const &  hits,
int  bxCrossing,
const HGCalGeometry geom,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 545 of file HGCDigitizer.cc.

548  {
549  if (nullptr == geom)
550  return;
551 
552  //configuration to apply for the computation of time-of-flight
553  auto weightToAbyEnergy = theDigitizer_->toaModeByEnergy();
554  auto tdcForToAOnset = theDigitizer_->tdcForToAOnset();
555  auto keV2fC = theDigitizer_->keV2fC();
556 
557  //create list of tuples (pos in container, RECO DetId, time) to be sorted first
558  int nchits = (int)hits->size();
559 
560  std::vector<HGCCaloHitTuple_t> hitRefs;
561  hitRefs.reserve(nchits);
562  for (int i = 0; i < nchits; ++i) {
563  const auto& the_hit = hits->at(i);
564  DetId id = simToReco(geom, the_hit.id());
565 
566  if (verbosity_ > 0) {
567  edm::LogVerbatim("HGCDigitizer") << "HGCDigitizer::i/p " << std::hex << the_hit.id() << " o/p " << id.rawId()
568  << std::dec;
569  }
570 
571  if (0 != id.rawId()) {
572  hitRefs.emplace_back(i, id.rawId(), (float)the_hit.time());
573  }
574  }
575 
576  std::sort(hitRefs.begin(), hitRefs.end(), this->orderByDetIdThenTime);
577  //loop over sorted hits
578  nchits = hitRefs.size();
579  for (int i = 0; i < nchits; ++i) {
580  const int hitidx = std::get<0>(hitRefs[i]);
581  const uint32_t id = std::get<1>(hitRefs[i]);
582 
583  //get the data for this cell, if not available then we skip it
584 
585  if (!validIds_.count(id))
586  continue;
587  HGCSimHitDataAccumulator::iterator simHitIt = simHitAccumulator_->emplace(id, HGCCellInfo()).first;
588 
589  if (id == 0)
590  continue; // to be ignored at RECO level
591 
592  const float toa = std::get<2>(hitRefs[i]);
593  const PCaloHit& hit = hits->at(hitidx);
594  const float charge = hit.energy() * 1e6 * keV2fC;
595 
596  //distance to the center of the detector
597  const float dist2center(getPositionDistance(geom, id));
598 
599  //hit time: [time()]=ns [centerDist]=cm [refSpeed_]=cm/ns + delay by 1ns
600  //accumulate in 15 buckets of 25ns (9 pre-samples, 1 in-time, 5 post-samples)
601  const float tof = toa - dist2center / refSpeed_ + tofDelay_;
602  const int itime = std::floor(tof / bxTime_) + 9;
603 
604  //no need to add bx crossing - tof comes already corrected from the mixing module
605  //itime += bxCrossing;
606  //itime += 9;
607 
608  if (itime < 0 || itime > (int)maxBx_)
609  continue;
610 
611  //check if time index is ok and store energy
612  if (itime >= (int)simHitIt->second.hit_info[0].size())
613  continue;
614 
615  (simHitIt->second).hit_info[0][itime] += charge;
616 
617  //for time-of-arrival: save the time-sorted list of timestamps with cumulative charge just above threshold
618  //working version with pileup only for in-time hits
619  int waferThickness = getCellThickness(geom, id);
620  bool orderChanged = false;
621  if (itime == (int)thisBx_) {
622  //if start empty => just add charge and time
623  if (hitRefs_bx0[id].empty()) {
624  hitRefs_bx0[id].emplace_back(charge, tof);
625 
626  } else if (tof <= hitRefs_bx0[id].back().second) {
627  //find position to insert new entry preserving time sorting
628  std::vector<std::pair<float, float>>::iterator findPos =
629  std::upper_bound(hitRefs_bx0[id].begin(),
630  hitRefs_bx0[id].end(),
631  std::pair<float, float>(0.f, tof),
632  [](const auto& i, const auto& j) { return i.second <= j.second; });
633 
634  std::vector<std::pair<float, float>>::iterator insertedPos = findPos;
635  if (findPos->second == tof) {
636  //just merge timestamps with exact timing
637  findPos->first += charge;
638  } else {
639  //insert new element cumulating the charge
640  insertedPos = hitRefs_bx0[id].insert(findPos,
641  (findPos == hitRefs_bx0[id].begin())
642  ? std::pair<float, float>(charge, tof)
643  : std::pair<float, float>((findPos - 1)->first + charge, tof));
644  }
645 
646  //cumulate the charge of new entry for all elements that follow in the sorted list
647  //and resize list accounting for cases when the inserted element itself crosses the threshold
648  for (std::vector<std::pair<float, float>>::iterator step = insertedPos; step != hitRefs_bx0[id].end(); ++step) {
649  if (step != insertedPos)
650  step->first += charge;
651  // resize the list stopping with the first timeStamp with cumulative charge above threshold
652  if (step->first > tdcForToAOnset[waferThickness - 1] && step->second != hitRefs_bx0[id].back().second) {
653  hitRefs_bx0[id].resize(std::upper_bound(hitRefs_bx0[id].begin(),
654  hitRefs_bx0[id].end(),
655  std::pair<float, float>(0.f, step->second),
656  [](const auto& i, const auto& j) { return i.second < j.second; }) -
657  hitRefs_bx0[id].begin());
658  for (auto stepEnd = step + 1; stepEnd != hitRefs_bx0[id].end(); ++stepEnd)
659  stepEnd->first += charge;
660  break;
661  }
662  }
663 
664  orderChanged = true;
665  } else {
666  //add new entry at the end of the list
667  if (hitRefs_bx0[id].back().first <= tdcForToAOnset[waferThickness - 1]) {
668  hitRefs_bx0[id].emplace_back(hitRefs_bx0[id].back().first + charge, tof);
669  }
670  }
671  }
672  float accChargeForToA = hitRefs_bx0[id].empty() ? 0.f : hitRefs_bx0[id].back().first;
673  //now compute the firing ToA through the interpolation of the consecutive time-stamps at threshold
674  if (weightToAbyEnergy)
675  (simHitIt->second).hit_info[1][itime] += charge * tof;
676  else if (accChargeForToA > tdcForToAOnset[waferThickness - 1] &&
677  ((simHitIt->second).hit_info[1][itime] == 0 || orderChanged == true)) {
678  float fireTDC = hitRefs_bx0[id].back().second;
679  if (hitRefs_bx0[id].size() > 1) {
680  float chargeBeforeThr = (hitRefs_bx0[id].end() - 2)->first;
681  float tofchargeBeforeThr = (hitRefs_bx0[id].end() - 2)->second;
682 
683  float deltaQ = accChargeForToA - chargeBeforeThr;
684  float deltaTOF = fireTDC - tofchargeBeforeThr;
685  fireTDC = (tdcForToAOnset[waferThickness - 1] - chargeBeforeThr) * deltaTOF / deltaQ + tofchargeBeforeThr;
686  }
687  (simHitIt->second).hit_info[1][itime] = fireTDC;
688  }
689  }
690  hitRefs.clear();
691 }

References bxTime_, ALCARECOTkAlJpsiMuMu_cff::charge, TauDecayModes::dec, relativeConstraints::empty, mps_fire::end, f, first, relativeConstraints::geom, hitRefs_bx0, hfClusterShapes_cfi::hits, mps_fire::i, triggerObjects_cff::id, createfilelist::int, dqmiolumiharvest::j, hgcalTriggerNtuples_cfi::keV2fC, maxBx_, orderByDetIdThenTime(), refSpeed_, edm::second(), simHitAccumulator_, findQualityFiles::size, jetUpdater_cfi::sort, theDigitizer_, thisBx_, tofDelay_, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, validIds_, trackerHitRTTI::vector, and verbosity_.

◆ accumulate() [3/3]

void HGCDigitizer::accumulate ( PileUpEventPrincipal const &  e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 402 of file HGCDigitizer.cc.

404  {
405  //get inputs
407  e.getByLabel(edm::InputTag("g4SimHits", hitCollection_), hits);
408 
409  if (!hits.isValid()) {
410  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
411  return;
412  }
413 
414  //accumulate for the simulated bunch crossing
415  if (nullptr != gHGCal_) {
416  accumulate(hits, e.bunchCrossing(), gHGCal_, hre);
417  } else {
418  throw cms::Exception("BadConfiguration") << "HGCDigitizer is not producing EE, FH, or BH digis!";
419  }
420 }

References accumulate(), MillePedeFileConverter_cfg::e, Exception, gHGCal_, hitCollection_, and hfClusterShapes_cfi::hits.

◆ accumulate_forPreMix() [1/4]

void HGCDigitizer::accumulate_forPreMix ( const PHGCSimAccumulator simAccumulator,
const bool  minbiasFlag 
)

Definition at line 692 of file HGCDigitizer.cc.

692  {
693  //configuration to apply for the computation of time-of-flight
694  auto weightToAbyEnergy = theDigitizer_->toaModeByEnergy();
695  auto tdcForToAOnset = theDigitizer_->tdcForToAOnset();
696 
697  if (nullptr != gHGCal_) {
698  loadSimHitAccumulator_forPreMix(*simHitAccumulator_,
700  gHGCal_,
701  hitRefs_bx0,
702  simAccumulator,
705  !weightToAbyEnergy,
706  tdcForToAOnset,
707  minbiasFlag,
709  thisBx_);
710  }
711 }

References gHGCal_, hitOrder_monitor, hitRefs_bx0, premixStage1MaxCharge_, premixStage1MinCharge_, pusimHitAccumulator_, simHitAccumulator_, theDigitizer_, and thisBx_.

◆ accumulate_forPreMix() [2/4]

void HGCDigitizer::accumulate_forPreMix ( edm::Event const &  e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 343 of file HGCDigitizer.cc.

345  {
346  //get inputs
347 
349  e.getByLabel(edm::InputTag("g4SimHits", hitCollection_), hits);
350  if (!hits.isValid()) {
351  edm::LogError("HGCDigitizer") << " @ accumulate_minbias : can't find " << hitCollection_
352  << " collection of g4SimHits";
353  return;
354  }
355 
356  //accumulate in-time the main event
357  if (nullptr != gHGCal_) {
359  } else {
360  throw cms::Exception("BadConfiguration") << "HGCDigitizer is not producing EE, FH, or BH digis!";
361  }
362 }

References MillePedeFileConverter_cfg::e, Exception, gHGCal_, hitCollection_, and hfClusterShapes_cfi::hits.

Referenced by HGCDigiProducer::accumulate(), accumulate_forPreMix(), PreMixingHGCalWorker::addPileups(), and PreMixingHGCalWorker::addSignals().

◆ accumulate_forPreMix() [3/4]

void HGCDigitizer::accumulate_forPreMix ( edm::Handle< edm::PCaloHitContainer > const &  hits,
int  bxCrossing,
const HGCalGeometry geom,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 423 of file HGCDigitizer.cc.

426  {
427  if (nullptr == geom)
428  return;
429 
430  auto keV2fC = theDigitizer_->keV2fC();
431  auto tdcForToAOnset = theDigitizer_->tdcForToAOnset();
432 
433  int nchits = (int)hits->size();
434  int count_thisbx = 0;
435  std::vector<HGCCaloHitTuple_t> hitRefs;
436  hitRefs.reserve(nchits);
437  for (int i = 0; i < nchits; ++i) {
438  const auto& the_hit = hits->at(i);
439  DetId id = simToReco(geom, the_hit.id());
440  // to be written the verbosity block
441  if (id.rawId() != 0) {
442  hitRefs.emplace_back(i, id.rawId(), (float)the_hit.time());
443  }
444  }
445  std::sort(hitRefs.begin(), hitRefs.end(), this->orderByDetIdThenTime);
446 
447  nchits = hitRefs.size();
448  for (int i = 0; i < nchits; ++i) {
449  const int hitidx = std::get<0>(hitRefs[i]);
450  const uint32_t id = std::get<1>(hitRefs[i]);
451  if (!validIds_.count(id))
452  continue;
453 
454  if (id == 0)
455  continue;
456 
457  const float toa = std::get<2>(hitRefs[i]);
458  const PCaloHit& hit = hits->at(hitidx);
459  const float charge = hit.energy() * 1e6 * keV2fC; // * getCCE(geom, id, cce_);
460 
461  const float dist2center(getPositionDistance(geom, id));
462  const float tof = toa - dist2center / refSpeed_ + tofDelay_;
463  const int itime = std::floor(tof / bxTime_) + 9;
464 
465  if (itime < 0 || itime > (int)maxBx_)
466  continue;
467 
468  if (itime >= (int)(maxBx_ + 1))
469  continue;
470 
471  int waferThickness = getCellThickness(geom, id);
472  if (itime == (int)thisBx_) {
473  ++count_thisbx;
474  if (PhitRefs_bx0[id].empty()) {
475  PhitRefs_bx0[id].emplace_back(charge, charge, tof);
476  } else if (tof > std::get<2>(PhitRefs_bx0[id].back())) {
477  PhitRefs_bx0[id].emplace_back(charge, charge + std::get<1>(PhitRefs_bx0[id].back()), tof);
478  } else if (tof == std::get<2>(PhitRefs_bx0[id].back())) {
479  std::get<0>(PhitRefs_bx0[id].back()) += charge;
480  std::get<1>(PhitRefs_bx0[id].back()) += charge;
481  } else {
482  //find position to insert new entry preserving time sorting
483  auto findPos = std::upper_bound(PhitRefs_bx0[id].begin(),
484  PhitRefs_bx0[id].end(),
485  hit_timeStamp(charge, 0.f, tof),
486  [](const auto& i, const auto& j) { return std::get<2>(i) <= std::get<2>(j); });
487 
488  auto insertedPos = findPos;
489 
490  if (tof == std::get<2>(*(findPos - 1))) {
491  std::get<0>(*(findPos - 1)) += charge;
492  std::get<1>(*(findPos - 1)) += charge;
493 
494  } else {
495  insertedPos = PhitRefs_bx0[id].insert(findPos,
496  (findPos == PhitRefs_bx0[id].begin())
497  ? hit_timeStamp(charge, charge, tof)
498  : hit_timeStamp(charge, charge + std::get<1>(*(findPos - 1)), tof));
499  }
500  //cumulate the charge of new entry for all elements that follow in the sorted list
501  //and resize list accounting for cases when the inserted element itself crosses the threshold
502 
503  for (auto step = insertedPos; step != PhitRefs_bx0[id].end(); ++step) {
504  if (step != insertedPos)
505  std::get<1>(*(step)) += charge;
506 
507  // resize the list stopping with the first timeStamp with cumulative charge above threshold
508  if (std::get<1>(*step) > tdcForToAOnset[waferThickness - 1] &&
509  std::get<2>(*step) != std::get<2>(PhitRefs_bx0[id].back())) {
510  PhitRefs_bx0[id].resize(
511  std::upper_bound(PhitRefs_bx0[id].begin(),
512  PhitRefs_bx0[id].end(),
513  hit_timeStamp(charge, 0.f, std::get<2>(*step)),
514  [](const auto& i, const auto& j) { return std::get<2>(i) < std::get<2>(j); }) -
515  PhitRefs_bx0[id].begin());
516  for (auto stepEnd = step + 1; stepEnd != PhitRefs_bx0[id].end(); ++stepEnd)
517  std::get<1>(*stepEnd) += charge;
518  break;
519  }
520  }
521  }
522  }
523  }
524 
525  for (const auto& hitCollection : PhitRefs_bx0) {
526  const uint32_t detectorId = hitCollection.first;
527  auto simHitIt = pusimHitAccumulator_->emplace(detectorId, HGCCellHitInfo()).first;
528 
529  for (const auto& hit_timestamp : PhitRefs_bx0[detectorId]) {
530  (simHitIt->second).PUhit_info[1][thisBx_].push_back(std::get<2>(hit_timestamp));
531  (simHitIt->second).PUhit_info[0][thisBx_].push_back(std::get<0>(hit_timestamp));
532  }
533  }
534 
535  if (nchits == 0) {
536  HGCPUSimHitDataAccumulator::iterator simHitIt = pusimHitAccumulator_->emplace(0, HGCCellHitInfo()).first;
537  (simHitIt->second).PUhit_info[1][9].push_back(0.0);
538  (simHitIt->second).PUhit_info[0][9].push_back(0.0);
539  }
540  hitRefs.clear();
541  PhitRefs_bx0.clear();
542 }

References bxTime_, ALCARECOTkAlJpsiMuMu_cff::charge, relativeConstraints::empty, mps_fire::end, f, relativeConstraints::geom, egammaRechitFilter_cfi::hitCollection, hfClusterShapes_cfi::hits, mps_fire::i, triggerObjects_cff::id, createfilelist::int, dqmiolumiharvest::j, hgcalTriggerNtuples_cfi::keV2fC, maxBx_, orderByDetIdThenTime(), PhitRefs_bx0, pusimHitAccumulator_, refSpeed_, jetUpdater_cfi::sort, theDigitizer_, thisBx_, tofDelay_, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, and validIds_.

◆ accumulate_forPreMix() [4/4]

void HGCDigitizer::accumulate_forPreMix ( PileUpEventPrincipal const &  e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 383 of file HGCDigitizer.cc.

385  {
387  e.getByLabel(edm::InputTag("g4SimHits", hitCollection_), hits);
388 
389  if (!hits.isValid()) {
390  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
391  return;
392  }
393 
394  if (nullptr != gHGCal_) {
395  accumulate_forPreMix(hits, e.bunchCrossing(), gHGCal_, hre);
396  } else {
397  throw cms::Exception("BadConfiguration") << "HGCDigitizer is not producing EE, FH, or BH digis!";
398  }
399 }

References accumulate_forPreMix(), MillePedeFileConverter_cfg::e, Exception, gHGCal_, hitCollection_, and hfClusterShapes_cfi::hits.

◆ checkPosition()

void HGCDigitizer::checkPosition ( const HGCalDigiCollection digis) const
private

Definition at line 742 of file HGCDigitizer.cc.

742  {
743  const double tol(0.5);
744  if (nullptr != gHGCal_) {
745  for (const auto& digi : *(digis)) {
746  const DetId& id = digi.id();
747  const GlobalPoint& global = gHGCal_->getPosition(id);
748  double r = global.perp();
749  double z = std::abs(global.z());
750  std::pair<double, double> zrange = gHGCal_->topology().dddConstants().rangeZ(true);
751  std::pair<double, double> rrange = gHGCal_->topology().dddConstants().rangeR(z, true);
752  bool ok = ((r >= rrange.first) && (r <= rrange.second) && (z >= zrange.first) && (z <= zrange.second));
753  std::string ck = (((r < rrange.first - tol) || (r > rrange.second + tol) || (z < zrange.first - tol) ||
754  (z > zrange.second + tol))
755  ? "***** ERROR *****"
756  : "");
757  bool val = gHGCal_->topology().valid(id);
758  if ((!ok) || (!val)) {
759  if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
760  edm::LogVerbatim("HGCDigitizer") << "Check " << HGCSiliconDetId(id) << " " << global << " R " << r << ":"
761  << rrange.first << ":" << rrange.second << " Z " << z << ":" << zrange.first
762  << ":" << zrange.second << " Flag " << ok << ":" << val << " " << ck;
763  } else if (id.det() == DetId::HGCalHSc) {
764  edm::LogVerbatim("HGCDigitizer") << "Check " << HGCScintillatorDetId(id) << " " << global << " R " << r << ":"
765  << rrange.first << ":" << rrange.second << " Z " << z << ":" << zrange.first
766  << ":" << zrange.second << " Flag " << ok << ":" << val << " " << ck;
767  } else if ((id.det() == DetId::Forward) && (id.subdetId() == static_cast<int>(HFNose))) {
768  edm::LogVerbatim("HGCDigitizer") << "Check " << HFNoseDetId(id) << " " << global << " R " << r << ":"
769  << rrange.first << ":" << rrange.second << " Z " << z << ":" << zrange.first
770  << ":" << zrange.second << " Flag " << ok << ":" << val << " " << ck;
771  } else {
772  edm::LogVerbatim("HGCDigitizer")
773  << "Check " << std::hex << id.rawId() << std::dec << " " << id.det() << ":" << id.subdetId() << " "
774  << global << " R " << r << ":" << rrange.first << ":" << rrange.second << " Z " << z << ":"
775  << zrange.first << ":" << zrange.second << " Flag " << ok << ":" << val << " " << ck;
776  }
777  }
778  }
779  }
780 }

References funct::abs(), HGCalTopology::dddConstants(), TauDecayModes::dec, DetId::Forward, HGCalGeometry::getPosition(), gHGCal_, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, convertSQLiteXML::ok, PV3DBase< T, PVType, FrameType >::perp(), alignCSCRings::r, HGCalDDDConstants::rangeR(), HGCalDDDConstants::rangeZ(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalGeometry::topology(), heppy_batch::val, HGCalTopology::valid(), z, PV3DBase< T, PVType, FrameType >::z(), and TkClusParameters_cff::zrange.

Referenced by finalizeEvent().

◆ digiCollection()

std::string HGCDigitizer::digiCollection ( )
inline

◆ finalizeEvent()

void HGCDigitizer::finalizeEvent ( edm::Event e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 299 of file HGCDigitizer.cc.

299  {
300  hitRefs_bx0.clear();
301  PhitRefs_bx0.clear();
302  hitOrder_monitor.clear();
303 
304  const CaloSubdetectorGeometry* theGeom = static_cast<const CaloSubdetectorGeometry*>(gHGCal_);
305 
306  ++nEvents_;
307 
308  unsigned idx = getType();
309  // release memory for unfilled parts of hash table
310  if (validIds_.size() * averageOccupancies_[idx] > simHitAccumulator_->size()) {
311  simHitAccumulator_->reserve(simHitAccumulator_->size());
312  pusimHitAccumulator_->reserve(simHitAccumulator_->size());
313  }
314  //update occupancy guess
315  const double thisOcc = simHitAccumulator_->size() / ((double)validIds_.size());
317 
318  if (premixStage1_) {
319  auto simRecord = std::make_unique<PHGCSimAccumulator>();
320 
321  if (!pusimHitAccumulator_->empty()) {
322  saveSimHitAccumulator_forPreMix(
324  }
325 
326  e.put(std::move(simRecord), digiCollection());
327 
328  } else {
329  auto digiResult = std::make_unique<HGCalDigiCollection>();
330  theDigitizer_->run(digiResult, *simHitAccumulator_, theGeom, validIds_, digitizationType_, hre);
331  edm::LogVerbatim("HGCDigitizer") << "HGCDigitizer:: finalize event - produced " << digiResult->size()
332  << " hits in det/subdet " << theDigitizer_->det() << "/"
333  << theDigitizer_->subdet();
334 #ifdef EDM_ML_DEBUG
335  checkPosition(&(*digiResult));
336 #endif
337  e.put(std::move(digiResult), digiCollection());
338  }
339 
342 }

References averageOccupancies_, checkPosition(), digiCollection(), digitizationType_, MillePedeFileConverter_cfg::e, getType(), gHGCal_, hitOrder_monitor, hitRefs_bx0, heavyIonCSV_trainingSettings::idx, eostools::move(), nEvents_, PhitRefs_bx0, premixStage1_, premixStage1MaxCharge_, premixStage1MinCharge_, pusimHitAccumulator_, simHitAccumulator_, theDigitizer_, and validIds_.

Referenced by HGCDigiProducer::finalizeEvent(), and PreMixingHGCalWorker::put().

◆ getType()

uint32_t HGCDigitizer::getType ( void  ) const
private

Definition at line 721 of file HGCDigitizer.cc.

721  {
723  switch (theDigitizer_->det()) {
724  case DetId::HGCalEE:
725  idx = 0;
726  break;
727  case DetId::HGCalHSi:
728  idx = 1;
729  break;
730  case DetId::HGCalHSc:
731  idx = 2;
732  break;
733  case DetId::Forward:
734  idx = 3;
735  break;
736  default:
737  break;
738  }
739  return idx;
740 }

References DetId::Forward, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, heavyIonCSV_trainingSettings::idx, SiStripPI::max, and theDigitizer_.

Referenced by finalizeEvent(), and initializeEvent().

◆ initializeEvent()

void HGCDigitizer::initializeEvent ( edm::Event const &  e,
edm::EventSetup const &  c 
)

actions at the start/end of event

Definition at line 269 of file HGCDigitizer.cc.

269  {
270  if (geomWatcher_.check(es)) {
271  std::unordered_set<DetId>().swap(validIds_);
272 
273  //get geometry
274  CaloGeometry const& geom = es.getData(geomToken_);
275 
276  gHGCal_ =
277  dynamic_cast<const HGCalGeometry*>(geom.getSubdetectorGeometry(theDigitizer_->det(), theDigitizer_->subdet()));
278 
279  int nadded(0);
280  //valid ID lists
281  if (nullptr != gHGCal_) {
282  getValidDetIds(gHGCal_, validIds_);
283  } else {
284  throw cms::Exception("BadConfiguration") << "HGCDigitizer is not producing EE, FH, or BH digis!";
285  }
286 
287  if (verbosity_ > 0)
288  edm::LogInfo("HGCDigitizer") << "Added " << nadded << ":" << validIds_.size() << " detIds without "
289  << hitCollection_ << " in first event processed" << std::endl;
290  }
291 
292  // reserve memory for a full detector
293  unsigned idx = getType();
296 }

References averageOccupancies_, edm::ESWatcher< T >::check(), Exception, relativeConstraints::geom, geomToken_, geomWatcher_, edm::EventSetup::getData(), getType(), gHGCal_, hitCollection_, heavyIonCSV_trainingSettings::idx, pusimHitAccumulator_, simHitAccumulator_, theDigitizer_, validIds_, and verbosity_.

Referenced by PreMixingHGCalWorker::initializeEvent(), and HGCDigiProducer::initializeEvent().

◆ orderByDetIdThenTime()

static bool HGCDigitizer::orderByDetIdThenTime ( const HGCCaloHitTuple_t a,
const HGCCaloHitTuple_t b 
)
inlinestatic

Definition at line 36 of file HGCDigitizer.h.

36  {
37  unsigned int detId_a(std::get<1>(a)), detId_b(std::get<1>(b));
38 
39  if (detId_a < detId_b)
40  return true;
41  if (detId_a > detId_b)
42  return false;
43 
44  double time_a(std::get<2>(a)), time_b(std::get<2>(b));
45  if (time_a < time_b)
46  return true;
47 
48  return false;
49  }

References a, and b.

Referenced by accumulate(), and accumulate_forPreMix().

◆ resetPUSimHitDataAccumulator()

void HGCDigitizer::resetPUSimHitDataAccumulator ( )
private

◆ resetSimHitDataAccumulator()

void HGCDigitizer::resetSimHitDataAccumulator ( )
private

Definition at line 714 of file HGCDigitizer.cc.

714  {
715  for (HGCSimHitDataAccumulator::iterator it = simHitAccumulator_->begin(); it != simHitAccumulator_->end(); it++) {
716  it->second.hit_info[0].fill(0.);
717  it->second.hit_info[1].fill(0.);
718  }
719 }

References simHitAccumulator_.

Member Data Documentation

◆ averageOccupancies_

std::array<double, 4> HGCDigitizer::averageOccupancies_
private

Definition at line 122 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and initializeEvent().

◆ bxTime_

double HGCDigitizer::bxTime_
private

Definition at line 95 of file HGCDigitizer.h.

Referenced by accumulate(), accumulate_forPreMix(), and HGCDigitizer().

◆ cce_

std::vector<float> HGCDigitizer::cce_
private

Definition at line 128 of file HGCDigitizer.h.

Referenced by HGCDigitizer().

◆ digiCollection_

std::string HGCDigitizer::digiCollection_
private

Definition at line 80 of file HGCDigitizer.h.

Referenced by digiCollection(), and HGCDigitizer().

◆ digitizationType_

int HGCDigitizer::digitizationType_
private

Definition at line 83 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

◆ ev_per_eh_pair_

double HGCDigitizer::ev_per_eh_pair_
private

Definition at line 95 of file HGCDigitizer.h.

◆ geomToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HGCDigitizer::geomToken_
private

Definition at line 107 of file HGCDigitizer.h.

Referenced by initializeEvent().

◆ geomWatcher_

edm::ESWatcher<CaloGeometryRecord> HGCDigitizer::geomWatcher_
private

Definition at line 108 of file HGCDigitizer.h.

Referenced by initializeEvent().

◆ gHGCal_

const HGCalGeometry* HGCDigitizer::gHGCal_ = nullptr
private

◆ hitCollection_

std::string HGCDigitizer::hitCollection_
private

Definition at line 80 of file HGCDigitizer.h.

Referenced by accumulate(), accumulate_forPreMix(), HGCDigitizer(), and initializeEvent().

◆ hitOrder_monitor

std::unordered_map<uint32_t, bool> HGCDigitizer::hitOrder_monitor
private

Definition at line 131 of file HGCDigitizer.h.

Referenced by accumulate_forPreMix(), and finalizeEvent().

◆ hitRefs_bx0

std::unordered_map<uint32_t, std::vector<std::pair<float, float> > > HGCDigitizer::hitRefs_bx0
private

Definition at line 129 of file HGCDigitizer.h.

Referenced by accumulate(), accumulate_forPreMix(), and finalizeEvent().

◆ maxBx_

const unsigned int HGCDigitizer::maxBx_ = 14
staticprivate

Definition at line 126 of file HGCDigitizer.h.

Referenced by accumulate(), and accumulate_forPreMix().

◆ maxSimHitsAccTime_

int HGCDigitizer::maxSimHitsAccTime_
private

Definition at line 94 of file HGCDigitizer.h.

Referenced by HGCDigitizer().

◆ nEvents_

uint32_t HGCDigitizer::nEvents_
private

Definition at line 123 of file HGCDigitizer.h.

Referenced by finalizeEvent().

◆ PhitRefs_bx0

std::unordered_map<uint32_t, std::vector<std::tuple<float, float, float> > > HGCDigitizer::PhitRefs_bx0
private

Definition at line 130 of file HGCDigitizer.h.

Referenced by accumulate_forPreMix(), and finalizeEvent().

◆ premixStage1_

bool HGCDigitizer::premixStage1_
private

Definition at line 86 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

◆ premixStage1MaxCharge_

double HGCDigitizer::premixStage1MaxCharge_
private

Definition at line 91 of file HGCDigitizer.h.

Referenced by accumulate_forPreMix(), finalizeEvent(), and HGCDigitizer().

◆ premixStage1MinCharge_

double HGCDigitizer::premixStage1MinCharge_
private

Definition at line 89 of file HGCDigitizer.h.

Referenced by accumulate_forPreMix(), finalizeEvent(), and HGCDigitizer().

◆ pusimHitAccumulator_

std::unique_ptr<hgc::HGCPUSimHitDataAccumulator> HGCDigitizer::pusimHitAccumulator_
private

Definition at line 97 of file HGCDigitizer.h.

Referenced by accumulate_forPreMix(), finalizeEvent(), and initializeEvent().

◆ refSpeed_

float HGCDigitizer::refSpeed_
private

Definition at line 116 of file HGCDigitizer.h.

Referenced by accumulate(), and accumulate_forPreMix().

◆ simHitAccumulator_

std::unique_ptr<hgc::HGCSimHitDataAccumulator> HGCDigitizer::simHitAccumulator_
private

◆ theDigitizer_

std::unique_ptr<HGCDigitizerBase> HGCDigitizer::theDigitizer_
private

◆ thisBx_

const unsigned int HGCDigitizer::thisBx_ = 9
staticprivate

Definition at line 127 of file HGCDigitizer.h.

Referenced by accumulate(), and accumulate_forPreMix().

◆ tofDelay_

float HGCDigitizer::tofDelay_
private

Definition at line 119 of file HGCDigitizer.h.

Referenced by accumulate(), accumulate_forPreMix(), and HGCDigitizer().

◆ validIds_

std::unordered_set<DetId> HGCDigitizer::validIds_
private

Definition at line 109 of file HGCDigitizer.h.

Referenced by accumulate(), accumulate_forPreMix(), finalizeEvent(), and initializeEvent().

◆ verbosity_

uint32_t HGCDigitizer::verbosity_
private

Definition at line 113 of file HGCDigitizer.h.

Referenced by accumulate(), HGCDigitizer(), and initializeEvent().

edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57
HGCDigitizer::thisBx_
static const unsigned int thisBx_
Definition: HGCDigitizer.h:127
pfDeepBoostedJetPreprocessParams_cfi.upper_bound
upper_bound
Definition: pfDeepBoostedJetPreprocessParams_cfi.py:16
mps_fire.i
i
Definition: mps_fire.py:428
hgc_digi::HGCCellHitInfo
Definition: HGCDigitizerTypes.h:25
HGCalTopology::dddConstants
const HGCalDDDConstants & dddConstants() const
Definition: HGCalTopology.h:98
TkClusParameters_cff.zrange
zrange
Definition: TkClusParameters_cff.py:7
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
HGCalDDDConstants::rangeZ
std::pair< double, double > rangeZ(bool reco) const
Definition: HGCalDDDConstants.cc:1011
HGCDigitizer::premixStage1MinCharge_
double premixStage1MinCharge_
Definition: HGCDigitizer.h:89
HGCDigitizer::simHitAccumulator_
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:96
HGCDigitizer::digiCollection
std::string digiCollection()
Definition: HGCDigitizer.h:76
HGCDigitizer::gHGCal_
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:110
step
step
Definition: StallMonitor.cc:94
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
HGCDigitizer::nEvents_
uint32_t nEvents_
Definition: HGCDigitizer.h:123
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
HGCDigitizer::maxBx_
static const unsigned int maxBx_
Definition: HGCDigitizer.h:126
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HGCDigitizer::refSpeed_
float refSpeed_
Definition: HGCDigitizer.h:116
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
HGCDigitizer::theDigitizer_
std::unique_ptr< HGCDigitizerBase > theDigitizer_
Definition: HGCDigitizer.h:104
HGCDigitizer::premixStage1_
bool premixStage1_
Definition: HGCDigitizer.h:86
HFNoseDetId
Definition: HFNoseDetId.h:22
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
edm::ConsumesCollector::esConsumes
auto esConsumes()
Definition: ConsumesCollector.h:97
edm::Handle
Definition: AssociativeIterator.h:50
groupFilesInBlocks.temp
list temp
Definition: groupFilesInBlocks.py:142
HGCSiliconDetId
Definition: HGCSiliconDetId.h:22
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
DetId
Definition: DetId.h:17
DetId::HGCalHSi
Definition: DetId.h:33
DetId::HGCalEE
Definition: DetId.h:32
HGCDigitizer::accumulate_forPreMix
void accumulate_forPreMix(edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
Definition: HGCDigitizer.cc:343
CaloGeometry
Definition: CaloGeometry.h:21
HFNose
Definition: ForwardSubdetector.h:11
mps_fire.end
end
Definition: mps_fire.py:242
DDAxes::z
relativeConstraints.geom
geom
Definition: relativeConstraints.py:72
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:55
HGCDigitizer::pusimHitAccumulator_
std::unique_ptr< hgc::HGCPUSimHitDataAccumulator > pusimHitAccumulator_
Definition: HGCDigitizer.h:97
egammaRechitFilter_cfi.hitCollection
hitCollection
Definition: egammaRechitFilter_cfi.py:10
HGCDigitizer::orderByDetIdThenTime
static bool orderByDetIdThenTime(const HGCCaloHitTuple_t &a, const HGCCaloHitTuple_t &b)
Definition: HGCDigitizer.h:36
Point3DBase< float, GlobalTag >
b
double b
Definition: hdecay.h:118
HGCalGeometry::topology
const HGCalTopology & topology() const
Definition: HGCalGeometry.h:111
HGCalDDDConstants::rangeR
std::pair< double, double > rangeR(double z, bool reco) const
Definition: HGCalDDDConstants.cc:966
first
auto first
Definition: CAHitNtupletGeneratorKernelsImpl.h:125
hgc_digi::HGCCellInfo
Definition: HGCDigitizerTypes.h:31
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
HGCDigitizer::hitCollection_
std::string hitCollection_
Definition: HGCDigitizer.h:80
edm::ParameterSet
Definition: ParameterSet.h:47
a
double a
Definition: hdecay.h:119
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:29
HGCDigitizer::accumulate
void accumulate(edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
handle SimHit accumulation
Definition: HGCDigitizer.cc:365
HGCDigitizer::maxSimHitsAccTime_
int maxSimHitsAccTime_
Definition: HGCDigitizer.h:94
HGCDigitizer::geomToken_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > geomToken_
Definition: HGCDigitizer.h:107
HGCDigitizer::hitRefs_bx0
std::unordered_map< uint32_t, std::vector< std::pair< float, float > > > hitRefs_bx0
Definition: HGCDigitizer.h:129
createfilelist.int
int
Definition: createfilelist.py:10
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
HGCalGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Definition: HGCalGeometry.cc:198
HGCDigitizer::digitizationType_
int digitizationType_
Definition: HGCDigitizer.h:83
PCaloHit
Definition: PCaloHit.h:8
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
HGCDigitizer::PhitRefs_bx0
std::unordered_map< uint32_t, std::vector< std::tuple< float, float, float > > > PhitRefs_bx0
Definition: HGCDigitizer.h:130
get
#define get
hgc_digi::HGCSimHitDataAccumulator
std::unordered_map< uint32_t, HGCCellInfo > HGCSimHitDataAccumulator
Definition: HGCDigitizerTypes.h:38
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HGCDigitizer::validIds_
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:109
HGCDigitizer::bxTime_
double bxTime_
Definition: HGCDigitizer.h:95
alignCSCRings.r
r
Definition: alignCSCRings.py:93
hit_timeStamp
std::tuple< float, float, float > hit_timeStamp
Definition: HGCDigitizer.cc:28
HGCDigitizer::tofDelay_
float tofDelay_
Definition: HGCDigitizer.h:119
HGCDigitizer::premixStage1MaxCharge_
double premixStage1MaxCharge_
Definition: HGCDigitizer.h:91
heppy_batch.val
val
Definition: heppy_batch.py:351
eostools.move
def move(src, dest)
Definition: eostools.py:511
HGCDigitizer::verbosity_
uint32_t verbosity_
Definition: HGCDigitizer.h:113
HGCDigitizer::digiCollection_
std::string digiCollection_
Definition: HGCDigitizer.h:80
hgc_digi::HGCPUSimHitDataAccumulator
std::unordered_map< uint32_t, HGCCellHitInfo > HGCPUSimHitDataAccumulator
Definition: HGCDigitizerTypes.h:39
HGCScintillatorDetId
Definition: HGCScintillatorDetId.h:23
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
DetId::HGCalHSc
Definition: DetId.h:34
relativeConstraints.empty
bool empty
Definition: relativeConstraints.py:46
Exception
Definition: hltDiff.cc:245
HGCDigitizer::geomWatcher_
edm::ESWatcher< CaloGeometryRecord > geomWatcher_
Definition: HGCDigitizer.h:108
HGCDigitizer::checkPosition
void checkPosition(const HGCalDigiCollection *digis) const
Definition: HGCDigitizer.cc:742
CaloSubdetectorGeometry
Definition: CaloSubdetectorGeometry.h:22
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HGCDigitizer::cce_
std::vector< float > cce_
Definition: HGCDigitizer.h:128
HGCalTopology::valid
bool valid(const DetId &id) const override
Is this a valid cell id.
Definition: HGCalTopology.cc:464
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
DetId::Forward
Definition: DetId.h:30
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
HGCDigitizer::getType
uint32_t getType() const
Definition: HGCDigitizer.cc:721
TauDecayModes.dec
dec
Definition: TauDecayModes.py:142
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
edm::InputTag
Definition: InputTag.h:15
hit
Definition: SiStripHitEffFromCalibTree.cc:88
HGCDigitizer::hitOrder_monitor
std::unordered_map< uint32_t, bool > hitOrder_monitor
Definition: HGCDigitizer.h:131
hgcalTriggerNtuples_cfi.keV2fC
keV2fC
Definition: hgcalTriggerNtuples_cfi.py:11
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
HGCDigitizer::averageOccupancies_
std::array< double, 4 > averageOccupancies_
Definition: HGCDigitizer.h:122