CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | 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 (PileUpEventPrincipal const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
 
template<typename GEOM >
void accumulate (edm::Handle< edm::PCaloHitContainer > const &hits, int bxCrossing, const GEOM *geom, CLHEP::HepRandomEngine *hre)
 
void beginRun (const edm::EventSetup &es)
 actions at the start/end of run More...
 
std::string digiCollection ()
 
void endRun ()
 
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...
 
bool producesEEDigis ()
 
bool producesHEbackDigis ()
 
bool producesHEfrontDigis ()
 
 ~HGCDigitizer ()
 

Static Public Member Functions

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

Private Member Functions

void resetSimHitDataAccumulator ()
 

Private Attributes

std::array< double, 3 > averageOccupancies_
 
double bxTime_
 
std::vector< float > cce_
 
std::string digiCollection_
 
int digitizationType_
 
double ev_per_eh_pair_
 
const HcalGeometrygHcal_
 
const HGCalGeometrygHGCal_
 
std::string hitCollection_
 
int maxSimHitsAccTime_
 
ForwardSubdetector mySubDet_
 
uint32_t nEvents_
 
float refSpeed_
 
std::unique_ptr< hgc::HGCSimHitDataAccumulatorsimHitAccumulator_
 
std::unique_ptr< HGCEEDigitizertheHGCEEDigitizer_
 
std::unique_ptr< HGCHEbackDigitizertheHGCHEbackDigitizer_
 
std::unique_ptr< HGCHEfrontDigitizertheHGCHEfrontDigitizer_
 
float tofDelay_
 
std::unordered_set< DetIdvalidIds_
 
uint32_t verbosity_
 

Detailed Description

Definition at line 28 of file HGCDigitizer.h.

Member Typedef Documentation

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

Definition at line 36 of file HGCDigitizer.h.

Constructor & Destructor Documentation

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

Definition at line 118 of file HGCDigitizer.cc.

References bxTime_, cce_, edm::ConsumesCollector::consumes(), digiCollection_, digitizationType_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HGCEE, HGCHEB, HGCHEF, hitCollection_, maxSimHitsAccTime_, mySubDet_, AlCaHLTBitMon_QueryRunRegistry::string, edm::swap(), groupFilesInBlocks::temp, theHGCEEDigitizer_, theHGCHEbackDigitizer_, theHGCHEfrontDigitizer_, tofDelay_, validIds_, and verbosity_.

119  :
122  refSpeed_(0.1*CLHEP::c_light), //[CLHEP::c_light]=mm/ns convert to cm/ns
123  averageOccupancies_(occupancyGuesses),
124  nEvents_(1)
125 {
126  //configure from cfg
127  hitCollection_ = ps.getParameter< std::string >("hitCollection");
128  digiCollection_ = ps.getParameter< std::string >("digiCollection");
129  maxSimHitsAccTime_ = ps.getParameter< uint32_t >("maxSimHitsAccTime");
130  bxTime_ = ps.getParameter< double >("bxTime");
131  digitizationType_ = ps.getParameter< uint32_t >("digitizationType");
132  verbosity_ = ps.getUntrackedParameter< uint32_t >("verbosity",0);
133  tofDelay_ = ps.getParameter< double >("tofDelay");
134 
135  std::unordered_set<DetId>().swap(validIds_);
136 
137  iC.consumes<std::vector<PCaloHit> >(edm::InputTag("g4SimHits",hitCollection_));
138  const auto& myCfg_ = ps.getParameter<edm::ParameterSet>("digiCfg");
139 
140  if( myCfg_.existsAs<std::vector<double> >( "chargeCollectionEfficiencies" ) ) {
141  cce_.clear();
142  const auto& temp = myCfg_.getParameter<std::vector<double> >("chargeCollectionEfficiencies");
143  for( double cce : temp ) {
144  cce_.push_back(cce);
145  }
146  } else {
147  std::vector<float>().swap(cce_);
148  }
149 
150  if(hitCollection_.find("HitsEE")!=std::string::npos) {
152  theHGCEEDigitizer_=std::unique_ptr<HGCEEDigitizer>(new HGCEEDigitizer(ps) );
153  }
154  if(hitCollection_.find("HitsHEfront")!=std::string::npos)
155  {
157  theHGCHEfrontDigitizer_=std::unique_ptr<HGCHEfrontDigitizer>(new HGCHEfrontDigitizer(ps) );
158  }
159  if(hitCollection_.find("HcalHits")!=std::string::npos)
160  {
162  theHGCHEbackDigitizer_=std::unique_ptr<HGCHEbackDigitizer>(new HGCHEbackDigitizer(ps) );
163  }
164 }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
uint32_t nEvents_
Definition: HGCDigitizer.h:116
T getUntrackedParameter(std::string const &, T const &) const
std::vector< float > cce_
Definition: HGCDigitizer.h:118
ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
int digitizationType_
Definition: HGCDigitizer.h:84
std::string hitCollection_
Definition: HGCDigitizer.h:81
void swap(Association< C > &lhs, Association< C > &rhs)
Definition: Association.h:116
std::string digiCollection_
Definition: HGCDigitizer.h:81
std::unordered_map< uint32_t, HGCCellInfo > HGCSimHitDataAccumulator
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
std::unique_ptr< HGCHEbackDigitizer > theHGCHEbackDigitizer_
Definition: HGCDigitizer.h:94
int maxSimHitsAccTime_
Definition: HGCDigitizer.h:87
std::unique_ptr< HGCHEfrontDigitizer > theHGCHEfrontDigitizer_
Definition: HGCDigitizer.h:95
std::unique_ptr< HGCEEDigitizer > theHGCEEDigitizer_
Definition: HGCDigitizer.h:93
std::array< double, 3 > averageOccupancies_
Definition: HGCDigitizer.h:115
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:89
double bxTime_
Definition: HGCDigitizer.h:88
uint32_t verbosity_
Definition: HGCDigitizer.h:106
HGCDigitizer::~HGCDigitizer ( )
inline

Definition at line 33 of file HGCDigitizer.h.

33 { }

Member Function Documentation

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

handle SimHit accumulation

Definition at line 244 of file HGCDigitizer.cc.

References Exception, edm::Event::getByLabel(), gHcal_, gHGCal_, hitCollection_, hfClusterShapes_cfi::hits, and edm::HandleBase::isValid().

Referenced by accumulate(), and orderByDetIdThenTime().

244  {
245 
246  //get inputs
248  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits);
249  if( !hits.isValid() ){
250  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
251  return;
252  }
253 
254  //accumulate in-time the main event
255  if( nullptr != gHGCal_ ) {
256  accumulate(hits, 0, gHGCal_, hre);
257  } else if( nullptr != gHcal_ ) {
258  accumulate(hits, 0, gHcal_, hre);
259  } else {
260  throw cms::Exception("BadConfiguration")
261  << "HGCDigitizer is not producing EE, FH, or BH digis!";
262  }
263 }
std::string hitCollection_
Definition: HGCDigitizer.h:81
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool isValid() const
Definition: HandleBase.h:74
void accumulate(edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
handle SimHit accumulation
const HcalGeometry * gHcal_
Definition: HGCDigitizer.h:100
void HGCDigitizer::accumulate ( PileUpEventPrincipal const &  e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 266 of file HGCDigitizer.cc.

References accumulate(), PileUpEventPrincipal::bunchCrossing(), Exception, PileUpEventPrincipal::getByLabel(), gHcal_, gHGCal_, hitCollection_, hfClusterShapes_cfi::hits, and edm::HandleBase::isValid().

266  {
267 
268  //get inputs
270  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits);
271  if( !hits.isValid() ){
272  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
273  return;
274  }
275 
276  //accumulate for the simulated bunch crossing
277  if( nullptr != gHGCal_ ) {
278  accumulate(hits, e.bunchCrossing(), gHGCal_, hre);
279  } else if ( nullptr != gHcal_ ) {
280  accumulate(hits, e.bunchCrossing(), gHcal_, hre);
281  } else {
282  throw cms::Exception("BadConfiguration")
283  << "HGCDigitizer is not producing EE, FH, or BH digis!";
284  }
285 }
std::string hitCollection_
Definition: HGCDigitizer.h:81
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool isValid() const
Definition: HandleBase.h:74
void accumulate(edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre)
handle SimHit accumulation
const HcalGeometry * gHcal_
Definition: HGCDigitizer.h:100
template<typename GEOM >
void HGCDigitizer::accumulate ( edm::Handle< edm::PCaloHitContainer > const &  hits,
int  bxCrossing,
const GEOM *  geom,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 289 of file HGCDigitizer.cc.

References bxTime_, cce_, ALCARECOTkAlJpsiMuMu_cff::charge, TauDecayModes::dec, f, HGCEE, HGCHEB, HGCHEF, mps_fire::i, createfilelist::int, hgcalDigitizer_cfi::keV2fC, mySubDet_, orderByDetIdThenTime(), producesEEDigis(), refSpeed_, simHitAccumulator_, HGCalUncalibRecHit_cfi::tdcOnset, theHGCEEDigitizer_, theHGCHEbackDigitizer_, theHGCHEfrontDigitizer_, tofDelay_, validIds_, and verbosity_.

292  {
293  if( nullptr == geom ) return;
294 
295 
296 
297  //configuration to apply for the computation of time-of-flight
298  bool weightToAbyEnergy(false);
299  float tdcOnset(0.f),keV2fC(0.f);
300  switch( mySubDet_ ) {
302  weightToAbyEnergy = theHGCEEDigitizer_->toaModeByEnergy();
303  tdcOnset = theHGCEEDigitizer_->tdcOnset();
304  keV2fC = theHGCEEDigitizer_->keV2fC();
305  break;
307  weightToAbyEnergy = theHGCHEfrontDigitizer_->toaModeByEnergy();
308  tdcOnset = theHGCHEfrontDigitizer_->tdcOnset();
309  keV2fC = theHGCHEfrontDigitizer_->keV2fC();
310  break;
312  weightToAbyEnergy = theHGCHEbackDigitizer_->toaModeByEnergy();
313  tdcOnset = theHGCHEbackDigitizer_->tdcOnset();
314  keV2fC = theHGCHEbackDigitizer_->keV2fC();
315  break;
316  default:
317  break;
318  }
319 
320  //create list of tuples (pos in container, RECO DetId, time) to be sorted first
321  int nchits=(int)hits->size();
322  std::vector< HGCCaloHitTuple_t > hitRefs;
323  hitRefs.reserve(nchits);
324  for(int i=0; i<nchits; ++i) {
325  const auto& the_hit = hits->at(i);
326 
327  DetId id = simToReco(geom,the_hit.id());
328 
329  if (verbosity_>0) {
330  if (producesEEDigis())
331  edm::LogInfo("HGCDigitizer") << " i/p " << std::hex << the_hit.id() << " o/p " << id.rawId() << std::dec << std::endl;
332  else
333  edm::LogInfo("HGCDigitizer") << " i/p " << std::hex << the_hit.id() << " o/p " << id.rawId() << std::dec << std::endl;
334  }
335 
336  if( 0 != id.rawId() ) {
337  hitRefs.emplace_back( i, id.rawId(), (float)the_hit.time() );
338  }
339  }
340  std::sort(hitRefs.begin(),hitRefs.end(),this->orderByDetIdThenTime);
341 
342  //loop over sorted hits
343  nchits = hitRefs.size();
344  for(int i=0; i<nchits; ++i) {
345  const int hitidx = std::get<0>(hitRefs[i]);
346  const uint32_t id = std::get<1>(hitRefs[i]);
347 
348  //get the data for this cell, if not available then we skip it
349 
350  if( !validIds_.count(id) ) continue;
351  HGCSimHitDataAccumulator::iterator simHitIt = simHitAccumulator_->emplace(id,HGCCellInfo()).first;
352 
353  if(id==0) continue; // to be ignored at RECO level
354 
355  const float toa = std::get<2>(hitRefs[i]);
356  const PCaloHit &hit=hits->at( hitidx );
357  const float charge = hit.energy()*1e6*keV2fC*getCCE(geom,id,cce_);
358 
359  //distance to the center of the detector
360  const float dist2center( getPositionDistance(geom,id) );
361 
362  //hit time: [time()]=ns [centerDist]=cm [refSpeed_]=cm/ns + delay by 1ns
363  //accumulate in 15 buckets of 25ns (9 pre-samples, 1 in-time, 5 post-samples)
364  const float tof = toa-dist2center/refSpeed_+tofDelay_ ;
365  const int itime= std::floor( tof/bxTime_ ) + 9;
366 
367  //no need to add bx crossing - tof comes already corrected from the mixing module
368  //itime += bxCrossing;
369  //itime += 9;
370 
371  if(itime<0 || itime>14) continue;
372 
373  //check if time index is ok and store energy
374  if(itime >= (int)simHitIt->second.hit_info[0].size() ) continue;
375 
376  (simHitIt->second).hit_info[0][itime] += charge;
377  float accCharge=(simHitIt->second).hit_info[0][itime];
378 
379  //time-of-arrival (check how to be used)
380  if(weightToAbyEnergy) (simHitIt->second).hit_info[1][itime] += charge*tof;
381  else if((simHitIt->second).hit_info[1][itime]==0) {
382  if( accCharge>tdcOnset)
383  {
384  //extrapolate linear using previous simhit if it concerns to the same DetId
385  float fireTDC=tof;
386  if(i>0)
387  {
388  uint32_t prev_id = std::get<1>(hitRefs[i-1]);
389  if(prev_id==id)
390  {
391  float prev_toa = std::get<2>(hitRefs[i-1]);
392  float prev_tof(prev_toa-dist2center/refSpeed_+tofDelay_);
393  //float prev_charge = std::get<3>(hitRefs[i-1]);
394  float deltaQ2TDCOnset = tdcOnset-((simHitIt->second).hit_info[0][itime]-charge);
395  float deltaQ = charge;
396  float deltaT = (tof-prev_tof);
397  fireTDC = deltaT*(deltaQ2TDCOnset/deltaQ)+prev_tof;
398  }
399  }
400 
401  (simHitIt->second).hit_info[1][itime]=fireTDC;
402  }
403  }
404  }
405  hitRefs.clear();
406 }
std::vector< float > cce_
Definition: HGCDigitizer.h:118
ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
bool producesEEDigis()
Definition: HGCDigitizer.h:67
double f[11][100]
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
std::unique_ptr< HGCHEbackDigitizer > theHGCHEbackDigitizer_
Definition: HGCDigitizer.h:94
Definition: DetId.h:18
std::unique_ptr< HGCHEfrontDigitizer > theHGCHEfrontDigitizer_
Definition: HGCDigitizer.h:95
std::unique_ptr< HGCEEDigitizer > theHGCEEDigitizer_
Definition: HGCDigitizer.h:93
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:89
static bool orderByDetIdThenTime(const HGCCaloHitTuple_t &a, const HGCCaloHitTuple_t &b)
Definition: HGCDigitizer.h:37
double bxTime_
Definition: HGCDigitizer.h:88
uint32_t verbosity_
Definition: HGCDigitizer.h:106
void HGCDigitizer::beginRun ( const edm::EventSetup es)

actions at the start/end of run

Definition at line 409 of file HGCDigitizer.cc.

References Exception, DetId::Forward, relativeConstraints::geom, edm::EventSetup::get(), CaloGeometry::getSubdetectorGeometry(), gHcal_, gHGCal_, DetId::Hcal, HcalEndcap, HGCEE, HGCHEF, hitCollection_, producesEEDigis(), producesHEbackDigis(), producesHEfrontDigis(), validIds_, and verbosity_.

Referenced by digiCollection().

410 {
411  //get geometry
413  es.get<CaloGeometryRecord>().get(geom);
414 
415  gHGCal_ = nullptr;
416  gHcal_ = nullptr;
417 
418  if( producesEEDigis() ) gHGCal_ = dynamic_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCEE));
419  if( producesHEfrontDigis() ) gHGCal_ = dynamic_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCHEF));
420  if( producesHEbackDigis() ) gHcal_ = dynamic_cast<const HcalGeometry*>(geom->getSubdetectorGeometry(DetId::Hcal, HcalEndcap));
421 
422  int nadded(0);
423  //valid ID lists
424  if( nullptr != gHGCal_ ) {
425  getValidDetIds( gHGCal_, validIds_ );
426  } else if( nullptr != gHcal_ ) {
427  getValidDetIds( gHcal_, validIds_ );
428  } else {
429  throw cms::Exception("BadConfiguration")
430  << "HGCDigitizer is not producing EE, FH, or BH digis!";
431  }
432 
433  if (verbosity_ > 0)
434  edm::LogInfo("HGCDigitizer")
435  << "Added " << nadded << ":" << validIds_.size()
436  << " detIds without " << hitCollection_
437  << " in first event processed" << std::endl;
438 }
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:45
std::string hitCollection_
Definition: HGCDigitizer.h:81
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool producesEEDigis()
Definition: HGCDigitizer.h:67
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
bool producesHEfrontDigis()
Definition: HGCDigitizer.h:68
bool producesHEbackDigis()
Definition: HGCDigitizer.h:69
const T & get() const
Definition: EventSetup.h:55
const HcalGeometry * gHcal_
Definition: HGCDigitizer.h:100
uint32_t verbosity_
Definition: HGCDigitizer.h:106
std::string HGCDigitizer::digiCollection ( )
inline

Definition at line 70 of file HGCDigitizer.h.

References beginRun(), digiCollection_, and endRun().

Referenced by finalizeEvent().

70 { return digiCollection_; }
std::string digiCollection_
Definition: HGCDigitizer.h:81
void HGCDigitizer::endRun ( )

Definition at line 441 of file HGCDigitizer.cc.

References edm::swap(), and validIds_.

Referenced by digiCollection().

442 {
443  std::unordered_set<DetId>().swap(validIds_);
444 }
void swap(Association< C > &lhs, Association< C > &rhs)
Definition: Association.h:116
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
void HGCDigitizer::finalizeEvent ( edm::Event e,
edm::EventSetup const &  c,
CLHEP::HepRandomEngine *  hre 
)

Definition at line 188 of file HGCDigitizer.cc.

References averageOccupancies_, digiCollection(), digitizationType_, gHcal_, gHGCal_, HGCEE, HGCHEB, HGCHEF, training_settings::idx, hpstanc_transforms::max, eostools::move(), mySubDet_, nEvents_, producesEEDigis(), producesHEbackDigis(), producesHEfrontDigis(), edm::Event::put(), simHitAccumulator_, theHGCEEDigitizer_, theHGCHEbackDigitizer_, theHGCHEfrontDigitizer_, and validIds_.

Referenced by orderByDetIdThenTime().

189 {
190 
191  const CaloSubdetectorGeometry* theGeom = ( nullptr == gHGCal_ ?
192  static_cast<const CaloSubdetectorGeometry*>(gHcal_) :
193  static_cast<const CaloSubdetectorGeometry*>(gHGCal_) );
194 
195  ++nEvents_;
197  switch(mySubDet_) {
199  idx = 0;
200  break;
202  idx = 1;
203  break;
205  idx = 2;
206  break;
207  default:
208  break;
209  }
210  // release memory for unfilled parts of hash table
211  if( validIds_.size()*averageOccupancies_[idx] > simHitAccumulator_->size() ) {
212  simHitAccumulator_->reserve(simHitAccumulator_->size());
213  }
214  //update occupancy guess
215  const double thisOcc = simHitAccumulator_->size()/((double)validIds_.size());
217 
218  if( producesEEDigis() )
219  {
220  std::unique_ptr<HGCEEDigiCollection> digiResult(new HGCEEDigiCollection() );
221  theHGCEEDigitizer_->run(digiResult,*simHitAccumulator_,theGeom,validIds_,digitizationType_, hre);
222  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " EE hits";
223  e.put(std::move(digiResult),digiCollection());
224  }
225  if( producesHEfrontDigis())
226  {
227  std::unique_ptr<HGCHEDigiCollection> digiResult(new HGCHEDigiCollection() );
229  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " HE front hits";
230  e.put(std::move(digiResult),digiCollection());
231  }
232  if( producesHEbackDigis() )
233  {
234  std::unique_ptr<HGCBHDigiCollection> digiResult(new HGCBHDigiCollection() );
236  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " HE back hits";
237  e.put(std::move(digiResult),digiCollection());
238  }
239 
241 }
uint32_t nEvents_
Definition: HGCDigitizer.h:116
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
edm::SortedCollection< HGCEEDataFrame > HGCEEDigiCollection
ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
int digitizationType_
Definition: HGCDigitizer.h:84
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool producesEEDigis()
Definition: HGCDigitizer.h:67
std::unordered_map< uint32_t, HGCCellInfo > HGCSimHitDataAccumulator
edm::SortedCollection< HGCBHDataFrame > HGCBHDigiCollection
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
bool producesHEfrontDigis()
Definition: HGCDigitizer.h:68
std::unique_ptr< HGCHEbackDigitizer > theHGCHEbackDigitizer_
Definition: HGCDigitizer.h:94
bool producesHEbackDigis()
Definition: HGCDigitizer.h:69
std::unique_ptr< HGCHEfrontDigitizer > theHGCHEfrontDigitizer_
Definition: HGCDigitizer.h:95
std::unique_ptr< HGCEEDigitizer > theHGCEEDigitizer_
Definition: HGCDigitizer.h:93
std::array< double, 3 > averageOccupancies_
Definition: HGCDigitizer.h:115
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:89
std::string digiCollection()
Definition: HGCDigitizer.h:70
edm::SortedCollection< HGCHEDataFrame > HGCHEDigiCollection
const HcalGeometry * gHcal_
Definition: HGCDigitizer.h:100
def move(src, dest)
Definition: eostools.py:510
void HGCDigitizer::initializeEvent ( edm::Event const &  e,
edm::EventSetup const &  c 
)

actions at the start/end of event

Definition at line 167 of file HGCDigitizer.cc.

References averageOccupancies_, HGCEE, HGCHEB, HGCHEF, training_settings::idx, hpstanc_transforms::max, mySubDet_, simHitAccumulator_, and validIds_.

Referenced by orderByDetIdThenTime().

168 {
169  // reserve memory for a full detector
171  switch(mySubDet_) {
173  idx = 0;
174  break;
176  idx = 1;
177  break;
179  idx = 2;
180  break;
181  default:
182  break;
183  }
184  simHitAccumulator_->reserve( averageOccupancies_[idx]*validIds_.size() );
185 }
ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
std::unordered_set< DetId > validIds_
Definition: HGCDigitizer.h:98
std::array< double, 3 > averageOccupancies_
Definition: HGCDigitizer.h:115
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:89
static bool HGCDigitizer::orderByDetIdThenTime ( const HGCCaloHitTuple_t a,
const HGCCaloHitTuple_t b 
)
inlinestatic

Definition at line 37 of file HGCDigitizer.h.

References accumulate(), EnergyCorrector::c, MillePedeFileConverter_cfg::e, finalizeEvent(), relativeConstraints::geom, hfClusterShapes_cfi::hits, and initializeEvent().

Referenced by accumulate().

38  {
39  unsigned int detId_a(std::get<1>(a)), detId_b(std::get<1>(b));
40 
41  if(detId_a<detId_b) return true;
42  if(detId_a>detId_b) return false;
43 
44  double time_a(std::get<2>(a)), time_b(std::get<2>(b));
45  if(time_a<time_b) return true;
46 
47  return false;
48  }
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121
bool HGCDigitizer::producesEEDigis ( )
inline

Definition at line 67 of file HGCDigitizer.h.

References HGCEE, and mySubDet_.

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

ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
bool HGCDigitizer::producesHEbackDigis ( )
inline

Definition at line 69 of file HGCDigitizer.h.

References HGCHEB, and mySubDet_.

Referenced by beginRun(), and finalizeEvent().

ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
bool HGCDigitizer::producesHEfrontDigis ( )
inline

Definition at line 68 of file HGCDigitizer.h.

References HGCHEF, and mySubDet_.

Referenced by beginRun(), and finalizeEvent().

ForwardSubdetector mySubDet_
Definition: HGCDigitizer.h:103
void HGCDigitizer::resetSimHitDataAccumulator ( )
private

Definition at line 447 of file HGCDigitizer.cc.

References relativeConstraints::geom, hfClusterShapes_cfi::hits, and simHitAccumulator_.

448 {
449  for( HGCSimHitDataAccumulator::iterator it = simHitAccumulator_->begin(); it!=simHitAccumulator_->end(); it++)
450  {
451  it->second.hit_info[0].fill(0.);
452  it->second.hit_info[1].fill(0.);
453  }
454 }
std::unique_ptr< hgc::HGCSimHitDataAccumulator > simHitAccumulator_
Definition: HGCDigitizer.h:89

Member Data Documentation

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

Definition at line 115 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and initializeEvent().

double HGCDigitizer::bxTime_
private

Definition at line 88 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

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

Definition at line 118 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

std::string HGCDigitizer::digiCollection_
private

Definition at line 81 of file HGCDigitizer.h.

Referenced by digiCollection(), and HGCDigitizer().

int HGCDigitizer::digitizationType_
private

Definition at line 84 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

double HGCDigitizer::ev_per_eh_pair_
private

Definition at line 88 of file HGCDigitizer.h.

const HcalGeometry* HGCDigitizer::gHcal_
private

Definition at line 100 of file HGCDigitizer.h.

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

const HGCalGeometry* HGCDigitizer::gHGCal_
private

Definition at line 99 of file HGCDigitizer.h.

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

std::string HGCDigitizer::hitCollection_
private

Definition at line 81 of file HGCDigitizer.h.

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

int HGCDigitizer::maxSimHitsAccTime_
private

Definition at line 87 of file HGCDigitizer.h.

Referenced by HGCDigitizer().

ForwardSubdetector HGCDigitizer::mySubDet_
private
uint32_t HGCDigitizer::nEvents_
private

Definition at line 116 of file HGCDigitizer.h.

Referenced by finalizeEvent().

float HGCDigitizer::refSpeed_
private

Definition at line 109 of file HGCDigitizer.h.

Referenced by accumulate().

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

Definition at line 93 of file HGCDigitizer.h.

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

std::unique_ptr<HGCHEbackDigitizer> HGCDigitizer::theHGCHEbackDigitizer_
private

Definition at line 94 of file HGCDigitizer.h.

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

std::unique_ptr<HGCHEfrontDigitizer> HGCDigitizer::theHGCHEfrontDigitizer_
private

Definition at line 95 of file HGCDigitizer.h.

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

float HGCDigitizer::tofDelay_
private

Definition at line 112 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

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

Definition at line 106 of file HGCDigitizer.h.

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