#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::ConsumesCollector &iC, edm::ProducerBase &parent)

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::ConsumesCollector &iC, edm::ProducerBase &parent)

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 48 of file FTLDigitizer.h.

Constructor & Destructor Documentation

template<class SensorPhysics , class ElectronicsSim >

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

inline

Definition at line 52 of file FTLDigitizer.h.

                                     :
     FTLDigitizerBase(config,iC,parent),
     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 62 of file FTLDigitizer.h.

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

62 { }

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 108 of file FTLDigitizer.h.

References edm::Event::getByLabel(), and trackerHits::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 117 of file FTLDigitizer.h.

References PileUpEventPrincipal::bunchCrossing(), PileUpEventPrincipal::getByLabel(), and trackerHits::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 126 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;
       
       //no need to add bx crossing - tof comes already corrected from the mixing module
       //itime += bxCrossing;
       //itime += 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 prev_charge = std::get<3>(hitRefs[i-1]);
               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 242 of file FTLDigitizer.h.

References edm::EventSetup::get(), 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 228 of file FTLDigitizer.h.

References hgcalDigitizer_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 222 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 85 of file FTLDigitizer.h.

                                       {
       FTLSimHitDataAccumulator().swap(simHitAccumulator_);
     }

Member Data Documentation

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 95 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 104 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 90 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 91 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 103 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 94 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 96 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 99 of file FTLDigitizer.h.

template<class SensorPhysics , class ElectronicsSim >

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

private

Definition at line 102 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 >