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

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

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

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

Referenced by accumulate(), and orderByDetIdThenTime().

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

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

234  {
235 
236  //get inputs
238  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits);
239  if( !hits.isValid() ){
240  edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
241  return;
242  }
243 
244  //accumulate for the simulated bunch crossing
245  if( nullptr != gHGCal_ ) {
246  accumulate(hits, e.bunchCrossing(), gHGCal_, hre);
247  } else if ( nullptr != gHcal_ ) {
248  accumulate(hits, e.bunchCrossing(), gHcal_, hre);
249  } else {
250  throw cms::Exception("BadConfiguration")
251  << "HGCDigitizer is not producing EE, FH, or BH digis!";
252  }
253 }
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 257 of file HGCDigitizer.cc.

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

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

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

Referenced by finalizeEvent().

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

Definition at line 409 of file HGCDigitizer.cc.

References edm::swap(), and validIds_.

Referenced by digiCollection().

410 {
411  std::unordered_set<DetId>().swap(validIds_);
412 }
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 156 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().

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

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

Referenced by orderByDetIdThenTime().

136 {
137  // reserve memory for a full detector
139  switch(mySubDet_) {
141  idx = 0;
142  break;
144  idx = 1;
145  break;
147  idx = 2;
148  break;
149  default:
150  break;
151  }
152  simHitAccumulator_->reserve( averageOccupancies_[idx]*validIds_.size() );
153 }
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 415 of file HGCDigitizer.cc.

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

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