#include <FTLDigitizer.h>

Inheritance diagram for ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >:

Public Member Functions
void	accumulate (edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *hre) override
	handle SimHit accumulation More...

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

void	accumulate (edm::Handle< edm::PSimHitContainer > const &hits, int bxCrossing, CLHEP::HepRandomEngine *hre) override

void	beginRun (const edm::EventSetup &es) override
	actions at the start/end of run More...

void	endRun () override

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

	FTLDigitizer (const edm::ParameterSet &config, edm::ProducesCollector producesCollector, edm::ConsumesCollector &iC)

void	initializeEvent (edm::Event const &e, edm::EventSetup const &c) override
	actions at the start/end of event More...

	~FTLDigitizer () override

Public Member Functions inherited from FTLDigitizerBase
	FTLDigitizerBase (const edm::ParameterSet &config, edm::ProducesCollector producesCollector, edm::ConsumesCollector &iC)

const std::string &	name () const

virtual	~FTLDigitizerBase ()

Private Member Functions
void	resetSimHitDataAccumulator ()

Private Attributes
const double	bxTime_

edm::ESHandle< FastTimeDDDConstants >	dddFTL_

SensorPhysics	deviceSim_

ElectronicsSim	electronicsSim_

edm::ESWatcher< IdealGeometryRecord >	idealGeomWatcher_

const int	maxSimHitsAccTime_

FTLSimHitDataAccumulator	simHitAccumulator_

const float	tofDelay_

std::unordered_set< DetId >	validIds_

Additional Inherited Members
Protected Attributes inherited from FTLDigitizerBase
const std::string	digiCollection_

const edm::InputTag	inputSimHits_

const ForwardSubdetector	mySubDet_

const float	refSpeed_

const uint32_t	verbosity_

Detailed Description

template<class SensorPhysics, class ElectronicsSim>
class ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >

Definition at line 50 of file FTLDigitizer.h.

Constructor & Destructor Documentation

template<class SensorPhysics , class ElectronicsSim >

ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::FTLDigitizer	(	const edm::ParameterSet &	config,
		edm::ProducesCollector	producesCollector,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 52 of file FTLDigitizer.h.

         : FTLDigitizerBase(config, producesCollector, iC),
           deviceSim_(config.getParameterSet("DeviceSimulation")),
           electronicsSim_(config.getParameterSet("ElectronicsSimulation")),
           maxSimHitsAccTime_(config.getParameter<uint32_t>("maxSimHitsAccTime")),
           bxTime_(config.getParameter<double>("bxTime")),
           tofDelay_(config.getParameter<double>("tofDelay")) {}

template<class SensorPhysics , class ElectronicsSim >

ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::~FTLDigitizer ( )

inlineoverride

Definition at line 60 of file FTLDigitizer.h.

References HltBtagPostValidation_cff::c, MillePedeFileConverter_cfg::e, and hfClusterShapes_cfi::hits.

60 {}

Member Function Documentation

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::accumulate	(	edm::Event const &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	hre
	)

overridevirtual

handle SimHit accumulation

Implements FTLDigitizerBase.

Definition at line 105 of file FTLDigitizer.h.

References edm::Event::getByLabel(), and FastTrackerRecHitCombiner_cfi::simHits.

                                                                                           {
     edm::Handle<edm::PSimHitContainer> simHits;
     e.getByLabel(inputSimHits_, simHits);
     accumulate(simHits, 0, hre);
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::accumulate	(	PileUpEventPrincipal const &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	hre
	)

overridevirtual

Implements FTLDigitizerBase.

Definition at line 114 of file FTLDigitizer.h.

References PileUpEventPrincipal::bunchCrossing(), PileUpEventPrincipal::getByLabel(), and FastTrackerRecHitCombiner_cfi::simHits.

                                                                                           {
     edm::Handle<edm::PSimHitContainer> simHits;
     e.getByLabel(inputSimHits_, simHits);
     accumulate(simHits, e.bunchCrossing(), hre);
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::accumulate	(	edm::Handle< edm::PSimHitContainer > const &	hits,
		int	bxCrossing,
		CLHEP::HepRandomEngine *	hre
	)

overridevirtual

Implements FTLDigitizerBase.

Definition at line 123 of file FTLDigitizer.h.

References ALCARECOTkAlJpsiMuMu_cff::charge, TauDecayModes::dec, f, mps_fire::i, createfilelist::int, ftl_digitizer::FTLHelpers::orderByDetIdThenTime(), and HGCalUncalibRecHit_cfi::tdcOnset.

                                                                                           {
     using namespace FTLHelpers;
     //configuration to apply for the computation of time-of-flight
     bool weightToAbyEnergy(false);
     float tdcOnset(0.f);
 
     //create list of tuples (pos in container, RECO DetId, time) to be sorted first
     int nchits = (int)hits->size();
     std::vector<FTLCaloHitTuple_t> hitRefs;
     hitRefs.reserve(nchits);
     for (int i = 0; i < nchits; ++i) {
       const auto& the_hit = hits->at(i);
 
       DetId id = (validIds_.count(the_hit.detUnitId()) ? the_hit.detUnitId() : 0);
 
       if (verbosity_ > 0) {
         edm::LogInfo("HGCDigitizer") << " i/p " << std::hex << the_hit.detUnitId() << std::dec << " o/p " << id.rawId()
                                      << std::endl;
       }
 
       if (0 != id.rawId()) {
         hitRefs.emplace_back(i, id.rawId(), the_hit.tof());
       }
     }
     std::sort(hitRefs.begin(), hitRefs.end(), FTLHelpers::orderByDetIdThenTime);
 
     //loop over sorted hits
     nchits = hitRefs.size();
     for (int i = 0; i < nchits; ++i) {
       const int hitidx = std::get<0>(hitRefs[i]);
       const uint32_t id = std::get<1>(hitRefs[i]);
 
       //get the data for this cell, if not available then we skip it
 
       if (!validIds_.count(id))
         continue;
       auto simHitIt = simHitAccumulator_.emplace(id, FTLCellInfo()).first;
 
       if (id == 0)
         continue;  // to be ignored at RECO level
 
       const float toa = std::get<2>(hitRefs[i]);
       const PSimHit& hit = hits->at(hitidx);
       const float charge = deviceSim_.getChargeForHit(hit);
 
       //distance to the center of the detector
       const float dist2center(0.1f * hit.entryPoint().mag());
 
       //hit time: [time()]=ns  [centerDist]=cm [refSpeed_]=cm/ns + delay by 1ns
       //accumulate in 15 buckets of 25ns (9 pre-samples, 1 in-time, 5 post-samples)
       const float tof = toa - dist2center / refSpeed_ + tofDelay_;
       const int itime = std::floor(tof / bxTime_) + 9;
 
       if (itime < 0 || itime > 14)
         continue;
 
       //check if time index is ok and store energy
       if (itime >= (int)simHitIt->second.hit_info[0].size())
         continue;
 
       (simHitIt->second).hit_info[0][itime] += charge;
       float accCharge = (simHitIt->second).hit_info[0][itime];
 
       //time-of-arrival (check how to be used)
       if (weightToAbyEnergy)
         (simHitIt->second).hit_info[1][itime] += charge * tof;
       else if ((simHitIt->second).hit_info[1][itime] == 0) {
         if (accCharge > tdcOnset) {
           //extrapolate linear using previous simhit if it concerns to the same DetId
           float fireTDC = tof;
           if (i > 0) {
             uint32_t prev_id = std::get<1>(hitRefs[i - 1]);
             if (prev_id == id) {
               float prev_toa = std::get<2>(hitRefs[i - 1]);
               float prev_tof(prev_toa - dist2center / refSpeed_ + tofDelay_);
               float deltaQ2TDCOnset = tdcOnset - ((simHitIt->second).hit_info[0][itime] - charge);
               float deltaQ = charge;
               float deltaT = (tof - prev_tof);
               fireTDC = deltaT * (deltaQ2TDCOnset / deltaQ) + prev_tof;
             }
           }
 
           (simHitIt->second).hit_info[1][itime] = fireTDC;
         }
       }
     }
     hitRefs.clear();
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::beginRun ( const edm::EventSetup & es )

overridevirtual

actions at the start/end of run

Get DDD constants

Implements FTLDigitizerBase.

Definition at line 235 of file FTLDigitizer.h.

References edm::EventSetup::get(), LEDCalibrationChannels::iphi, edm::swap(), and ecaldqm::zside().

                                                                                     {
     if (idealGeomWatcher_.check(es)) {
       es.get<IdealGeometryRecord>().get(dddFTL_);
       {  // force scope for the temporary nameless unordered_set
         std::unordered_set<DetId>().swap(validIds_);
       }
       // recalculate valid detids
       for (int zside = -1; zside <= 1; zside += 2) {
         for (unsigned type = 1; type <= 2; ++type) {
           for (unsigned izeta = 0; izeta < 1 << 10; ++izeta) {
             for (unsigned iphi = 0; iphi < 1 << 10; ++iphi) {
               if (dddFTL_->isValidXY(type, izeta, iphi)) {
                 validIds_.emplace(FastTimeDetId(type, izeta, iphi, zside));
               }
             }
           }
         }
       }
       validIds_.reserve(validIds_.size());
     }
     deviceSim_.getEventSetup(es);
     electronicsSim_.getEventSetup(es);
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::endRun ( )

inlineoverridevirtual

Implements FTLDigitizerBase.

Definition at line 81 of file FTLDigitizer.h.

81 {}

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::finalizeEvent	(	edm::Event &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	hre
	)

overridevirtual

Implements FTLDigitizerBase.

Definition at line 221 of file FTLDigitizer.h.

References apvshotsanalyzer_cfi::digiCollection, eostools::move(), and edm::Event::put().

                                                                                              {
     auto digiCollection = std::make_unique<FTLDigiCollection>();
 
     electronicsSim_.run(simHitAccumulator_, *digiCollection);
 
     e.put(std::move(digiCollection), digiCollection_);
 
     //release memory for next event
     resetSimHitDataAccumulator();
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::initializeEvent	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

overridevirtual

actions at the start/end of event

Implements FTLDigitizerBase.

Definition at line 215 of file FTLDigitizer.h.

                                                                                                              {
     deviceSim_.getEvent(e);
     electronicsSim_.getEvent(e);
   }

template<class SensorPhysics , class ElectronicsSim >

void ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::resetSimHitDataAccumulator ( )

inlineprivate

Definition at line 84 of file FTLDigitizer.h.

84 { FTLSimHitDataAccumulator().swap(simHitAccumulator_); }

ftl_digitizer::FTLSimHitDataAccumulator

std::unordered_map< uint32_t, FTLCellInfo > FTLSimHitDataAccumulator

Definition: FTLDigitizerTypes.h:21

ftl_digitizer::FTLDigitizer::simHitAccumulator_

FTLSimHitDataAccumulator simHitAccumulator_

Definition: FTLDigitizer.h:93

Member Data Documentation

template<class SensorPhysics , class ElectronicsSim >

const double ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::bxTime_

private

Definition at line 92 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

edm::ESHandle<FastTimeDDDConstants> ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::dddFTL_

private

Definition at line 101 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

SensorPhysics ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::deviceSim_

private

Definition at line 87 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

ElectronicsSim ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::electronicsSim_

private

Definition at line 88 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

edm::ESWatcher<IdealGeometryRecord> ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::idealGeomWatcher_

private

Definition at line 100 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

const int ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::maxSimHitsAccTime_

private

Definition at line 91 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

FTLSimHitDataAccumulator ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::simHitAccumulator_

private

Definition at line 93 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

const float ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::tofDelay_

private

Definition at line 96 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

std::unordered_set<DetId> ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >::validIds_

private

Definition at line 99 of file FTLDigitizer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<class SensorPhysics, class ElectronicsSim> class ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<class SensorPhysics, class ElectronicsSim>
class ftl_digitizer::FTLDigitizer< SensorPhysics, ElectronicsSim >