CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
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::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::HGCSimHitDataAccumulator
simHitAccumulator_
 
std::unique_ptr< HGCEEDigitizertheHGCEEDigitizer_
 
std::unique_ptr
< HGCHEbackDigitizer
theHGCHEbackDigitizer_
 
std::unique_ptr
< HGCHEfrontDigitizer
theHGCHEfrontDigitizer_
 
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 100 of file HGCDigitizer.cc.

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

101  :
104  refSpeed_(0.1*CLHEP::c_light), //[CLHEP::c_light]=mm/ns convert to cm/ns
105  averageOccupancies_(occupancyGuesses),
106  nEvents_(1)
107 {
108  //configure from cfg
109  hitCollection_ = ps.getParameter< std::string >("hitCollection");
110  digiCollection_ = ps.getParameter< std::string >("digiCollection");
111  maxSimHitsAccTime_ = ps.getParameter< uint32_t >("maxSimHitsAccTime");
112  bxTime_ = ps.getParameter< double >("bxTime");
113  digitizationType_ = ps.getParameter< uint32_t >("digitizationType");
114  verbosity_ = ps.getUntrackedParameter< uint32_t >("verbosity",0);
115  tofDelay_ = ps.getParameter< double >("tofDelay");
116 
117  std::unordered_set<DetId>().swap(validIds_);
118 
119  iC.consumes<std::vector<PCaloHit> >(edm::InputTag("g4SimHits",hitCollection_));
120 
121  if(hitCollection_.find("HitsEE")!=std::string::npos) {
123  theHGCEEDigitizer_=std::unique_ptr<HGCEEDigitizer>(new HGCEEDigitizer(ps) );
124  }
125  if(hitCollection_.find("HitsHEfront")!=std::string::npos)
126  {
128  theHGCHEfrontDigitizer_=std::unique_ptr<HGCHEfrontDigitizer>(new HGCHEfrontDigitizer(ps) );
129  }
130  if(hitCollection_.find("HcalHits")!=std::string::npos)
131  {
133  theHGCHEbackDigitizer_=std::unique_ptr<HGCHEbackDigitizer>(new HGCHEbackDigitizer(ps) );
134  }
135 }
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
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 215 of file HGCDigitizer.cc.

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

Referenced by accumulate().

215  {
216 
217  //get inputs
219  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits);
220  if( !hits.isValid() ){
221  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
222  return;
223  }
224 
225  //accumulate in-time the main event
226  if( nullptr != gHGCal_ ) {
227  accumulate(hits, 0, gHGCal_, hre);
228  } else if( nullptr != gHcal_ ) {
229  accumulate(hits, 0, gHcal_, hre);
230  } else {
231  throw cms::Exception("BadConfiguration")
232  << "HGCDigitizer is not producing EE, FH, or BH digis!";
233  }
234 }
std::string hitCollection_
Definition: HGCDigitizer.h:81
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool isValid() const
Definition: HandleBase.h:75
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 237 of file HGCDigitizer.cc.

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

237  {
238 
239  //get inputs
241  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits);
242  if( !hits.isValid() ){
243  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
244  return;
245  }
246 
247  //accumulate for the simulated bunch crossing
248  if( nullptr != gHGCal_ ) {
249  accumulate(hits, e.bunchCrossing(), gHGCal_, hre);
250  } else if ( nullptr != gHcal_ ) {
251  accumulate(hits, e.bunchCrossing(), gHcal_, hre);
252  } else {
253  throw cms::Exception("BadConfiguration")
254  << "HGCDigitizer is not producing EE, FH, or BH digis!";
255  }
256 }
std::string hitCollection_
Definition: HGCDigitizer.h:81
const HGCalGeometry * gHGCal_
Definition: HGCDigitizer.h:99
bool isValid() const
Definition: HandleBase.h:75
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 260 of file HGCDigitizer.cc.

References bxTime_, RecoTauCleanerPlugins::charge, TauDecayModes::dec, f, HGCEE, HGCHEB, HGCHEF, i, mySubDet_, orderByDetIdThenTime(), producesEEDigis(), refSpeed_, simHitAccumulator_, theHGCEEDigitizer_, theHGCHEbackDigitizer_, theHGCHEfrontDigitizer_, tofDelay_, validIds_, and verbosity_.

263  {
264  if( nullptr == geom ) return;
265 
266 
267 
268  //configuration to apply for the computation of time-of-flight
269  bool weightToAbyEnergy(false);
270  float tdcOnset(0.f),keV2fC(0.f);
271  switch( mySubDet_ ) {
273  weightToAbyEnergy = theHGCEEDigitizer_->toaModeByEnergy();
274  tdcOnset = theHGCEEDigitizer_->tdcOnset();
275  keV2fC = theHGCEEDigitizer_->keV2fC();
276  break;
278  weightToAbyEnergy = theHGCHEfrontDigitizer_->toaModeByEnergy();
279  tdcOnset = theHGCHEfrontDigitizer_->tdcOnset();
280  keV2fC = theHGCHEfrontDigitizer_->keV2fC();
281  break;
283  weightToAbyEnergy = theHGCHEbackDigitizer_->toaModeByEnergy();
284  tdcOnset = theHGCHEbackDigitizer_->tdcOnset();
285  keV2fC = theHGCHEbackDigitizer_->keV2fC();
286  break;
287  default:
288  break;
289  }
290 
291  //create list of tuples (pos in container, RECO DetId, time) to be sorted first
292  int nchits=(int)hits->size();
293  std::vector< HGCCaloHitTuple_t > hitRefs;
294  hitRefs.reserve(nchits);
295  for(int i=0; i<nchits; ++i) {
296  const auto& the_hit = hits->at(i);
297 
298  DetId id = simToReco(geom,the_hit.id());
299 
300  if (verbosity_>0) {
301  if (producesEEDigis())
302  edm::LogInfo("HGCDigitizer") << " i/p " << std::hex << the_hit.id() << std::dec << " o/p " << id.rawId() << std::endl;
303  else
304  edm::LogInfo("HGCDigitizer") << " i/p " << std::hex << the_hit.id() << std::dec << " o/p " << id.rawId() << std::endl;
305  }
306 
307  if( 0 != id.rawId() ) {
308  hitRefs.emplace_back( i, id.rawId(), (float)the_hit.time() );
309  }
310  }
311  std::sort(hitRefs.begin(),hitRefs.end(),this->orderByDetIdThenTime);
312 
313  //loop over sorted hits
314  nchits = hitRefs.size();
315  for(int i=0; i<nchits; ++i) {
316  const int hitidx = std::get<0>(hitRefs[i]);
317  const uint32_t id = std::get<1>(hitRefs[i]);
318 
319  //get the data for this cell, if not available then we skip it
320 
321  if( !validIds_.count(id) ) continue;
322  HGCSimHitDataAccumulator::iterator simHitIt = simHitAccumulator_->emplace(id,HGCCellInfo()).first;
323 
324  if(id==0) continue; // to be ignored at RECO level
325 
326  const float toa = std::get<2>(hitRefs[i]);
327  const PCaloHit &hit=hits->at( hitidx );
328  const float charge = hit.energy()*1e6*keV2fC;
329 
330  //distance to the center of the detector
331  const float dist2center( getPositionDistance(geom,id) );
332 
333  //hit time: [time()]=ns [centerDist]=cm [refSpeed_]=cm/ns + delay by 1ns
334  //accumulate in 15 buckets of 25ns (9 pre-samples, 1 in-time, 5 post-samples)
335  const float tof = toa-dist2center/refSpeed_+tofDelay_ ;
336  const int itime= std::floor( tof/bxTime_ ) + 9;
337 
338  //no need to add bx crossing - tof comes already corrected from the mixing module
339  //itime += bxCrossing;
340  //itime += 9;
341 
342  if(itime<0 || itime>14) continue;
343 
344  //check if time index is ok and store energy
345  if(itime >= (int)simHitIt->second.hit_info[0].size() ) continue;
346 
347  (simHitIt->second).hit_info[0][itime] += charge;
348  float accCharge=(simHitIt->second).hit_info[0][itime];
349 
350  //time-of-arrival (check how to be used)
351  if(weightToAbyEnergy) (simHitIt->second).hit_info[1][itime] += charge*tof;
352  else if((simHitIt->second).hit_info[1][itime]==0) {
353  if( accCharge>tdcOnset)
354  {
355  //extrapolate linear using previous simhit if it concerns to the same DetId
356  float fireTDC=tof;
357  if(i>0)
358  {
359  uint32_t prev_id = std::get<1>(hitRefs[i-1]);
360  if(prev_id==id)
361  {
362  float prev_toa = std::get<2>(hitRefs[i-1]);
363  float prev_tof(prev_toa-dist2center/refSpeed_+tofDelay_);
364  //float prev_charge = std::get<3>(hitRefs[i-1]);
365  float deltaQ2TDCOnset = tdcOnset-((simHitIt->second).hit_info[0][itime]-charge);
366  float deltaQ = charge;
367  float deltaT = (tof-prev_tof);
368  fireTDC = deltaT*(deltaQ2TDCOnset/deltaQ)+prev_tof;
369  }
370  }
371 
372  (simHitIt->second).hit_info[1][itime]=fireTDC;
373  }
374  }
375  }
376  hitRefs.clear();
377 }
int i
Definition: DBlmapReader.cc:9
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 380 of file HGCDigitizer.cc.

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

381 {
382  //get geometry
384  es.get<CaloGeometryRecord>().get(geom);
385 
386  gHGCal_ = nullptr;
387  gHcal_ = nullptr;
388 
389  if( producesEEDigis() ) gHGCal_ = dynamic_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCEE));
390  if( producesHEfrontDigis() ) gHGCal_ = dynamic_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCHEF));
391  if( producesHEbackDigis() ) gHcal_ = dynamic_cast<const HcalGeometry*>(geom->getSubdetectorGeometry(DetId::Hcal, HcalEndcap));
392 
393  int nadded(0);
394  //valid ID lists
395  if( nullptr != gHGCal_ ) {
396  getValidDetIds( gHGCal_, validIds_ );
397  } else if( nullptr != gHcal_ ) {
398  getValidDetIds( gHcal_, validIds_ );
399  } else {
400  throw cms::Exception("BadConfiguration")
401  << "HGCDigitizer is not producing EE, FH, or BH digis!";
402  }
403 
404  if (verbosity_ > 0)
405  edm::LogInfo("HGCDigitizer")
406  << "Added " << nadded << ":" << validIds_.size()
407  << " detIds without " << hitCollection_
408  << " in first event processed" << std::endl;
409 }
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:56
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 digiCollection_.

Referenced by finalizeEvent().

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

Definition at line 412 of file HGCDigitizer.cc.

References edm::swap(), and validIds_.

413 {
414  std::unordered_set<DetId>().swap(validIds_);
415 }
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 159 of file HGCDigitizer.cc.

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

160 {
161 
162  const CaloSubdetectorGeometry* theGeom = ( nullptr == gHGCal_ ?
163  static_cast<const CaloSubdetectorGeometry*>(gHcal_) :
164  static_cast<const CaloSubdetectorGeometry*>(gHGCal_) );
165 
166  ++nEvents_;
167  unsigned idx = std::numeric_limits<unsigned>::max();
168  switch(mySubDet_) {
170  idx = 0;
171  break;
173  idx = 1;
174  break;
176  idx = 2;
177  break;
178  default:
179  break;
180  }
181  // release memory for unfilled parts of hash table
182  if( validIds_.size()*averageOccupancies_[idx] > simHitAccumulator_->size() ) {
183  simHitAccumulator_->reserve(simHitAccumulator_->size());
184  }
185  //update occupancy guess
186  const double thisOcc = simHitAccumulator_->size()/((double)validIds_.size());
187  averageOccupancies_[idx] = (averageOccupancies_[idx]*(nEvents_-1) + thisOcc)/nEvents_;
188 
189  if( producesEEDigis() )
190  {
191  std::unique_ptr<HGCEEDigiCollection> digiResult(new HGCEEDigiCollection() );
192  theHGCEEDigitizer_->run(digiResult,*simHitAccumulator_,theGeom,validIds_,digitizationType_, hre);
193  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " EE hits";
194  e.put(std::move(digiResult),digiCollection());
195  }
196  if( producesHEfrontDigis())
197  {
198  std::unique_ptr<HGCHEDigiCollection> digiResult(new HGCHEDigiCollection() );
200  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " HE front hits";
201  e.put(std::move(digiResult),digiCollection());
202  }
203  if( producesHEbackDigis() )
204  {
205  std::unique_ptr<HGCBHDigiCollection> digiResult(new HGCBHDigiCollection() );
207  edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() << " HE back hits";
208  e.put(std::move(digiResult),digiCollection());
209  }
210 
212 }
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
def move
Definition: eostools.py:510
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
void HGCDigitizer::initializeEvent ( edm::Event const &  e,
edm::EventSetup const &  c 
)

actions at the start/end of event

Definition at line 138 of file HGCDigitizer.cc.

References averageOccupancies_, HGCEE, HGCHEB, HGCHEF, bookConverter::max, mySubDet_, simHitAccumulator_, and validIds_.

139 {
140  // reserve memory for a full detector
141  unsigned idx = std::numeric_limits<unsigned>::max();
142  switch(mySubDet_) {
144  idx = 0;
145  break;
147  idx = 1;
148  break;
150  idx = 2;
151  break;
152  default:
153  break;
154  }
155  simHitAccumulator_->reserve( averageOccupancies_[idx]*validIds_.size() );
156 }
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.

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 418 of file HGCDigitizer.cc.

References simHitAccumulator_.

419 {
420  for( HGCSimHitDataAccumulator::iterator it = simHitAccumulator_->begin(); it!=simHitAccumulator_->end(); it++)
421  {
422  it->second.hit_info[0].fill(0.);
423  it->second.hit_info[1].fill(0.);
424  }
425 }
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::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().