#include <BTLDeviceSim.h>

Public Member Functions
	BTLDeviceSim (const edm::ParameterSet &pset)

void	getEvent (const edm::Event &evt)

void	getEventSetup (const edm::EventSetup &evt)

void	getHitsResponse (const std::vector< std::tuple< int, uint32_t, float > > &hitRefs, const edm::Handle< edm::PSimHitContainer > &hits, mtd_digitizer::MTDSimHitDataAccumulator simHitAccumulator, CLHEP::HepRandomEngine hre)

Private Attributes
const float	bxTime_

const float	LightCollEff_

const float	LightCollTime_

const float	LightYield_

const float	PDE_

const float	refSpeed_

const float	smearLightCollTime_

Detailed Description

Definition at line 18 of file BTLDeviceSim.h.

Constructor & Destructor Documentation

BTLDeviceSim::BTLDeviceSim ( const edm::ParameterSet & pset )

Definition at line 11 of file BTLDeviceSim.cc.

                                                       : 
   refSpeed_( 0.1*CLHEP::c_light ),
   bxTime_(pset.getParameter<double>("bxTime") ),
   LightYield_(pset.getParameter<double>("LightYield")),
   LightCollEff_(pset.getParameter<double>("LightCollectionEff")),
   LightCollTime_(pset.getParameter<double>("LightCollectionTime")),
   smearLightCollTime_(pset.getParameter<double>("smearLightCollectionTime")),
   PDE_(pset.getParameter<double>("PhotonDetectionEff")) { }

Member Function Documentation

void BTLDeviceSim::getEvent ( const edm::Event & evt )

inline

Definition at line 24 of file BTLDeviceSim.h.

24 { }

void BTLDeviceSim::getEventSetup ( const edm::EventSetup & evt )

inline

Definition at line 26 of file BTLDeviceSim.h.

References hfClusterShapes_cfi::hits.

26 { }

void BTLDeviceSim::getHitsResponse	(	const std::vector< std::tuple< int, uint32_t, float > > &	hitRefs,
		const edm::Handle< edm::PSimHitContainer > &	hits,
		mtd_digitizer::MTDSimHitDataAccumulator *	simHitAccumulator,
		CLHEP::HepRandomEngine *	hre
	)

Definition at line 21 of file BTLDeviceSim.cc.

References bxTime_, PSimHit::energyLoss(), PSimHit::entryPoint(), DetId::Forward, LightCollEff_, LightCollTime_, LightYield_, PV3DBase< T, PVType, FrameType >::mag(), PDE_, refSpeed_, and smearLightCollTime_.

                                              {
 
   //loop over sorted simHits
   const int nchits = hitRefs.size();
   for(int ihit=0; ihit<nchits; ++ihit) {
 
     const int hitidx   = std::get<0>(hitRefs[ihit]);
     const uint32_t id  = std::get<1>(hitRefs[ihit]);
     const MTDDetId detId(id);
     const PSimHit &hit = hits->at( hitidx );     
     
     // --- Safety check on the detector ID
     if ( detId.det()!=DetId::Forward || detId.mtdSubDetector()!=1 ) continue;
 
     // --- Store the detector element ID as a key of the MTDSimHitDataAccumulator map
     auto simHitIt = simHitAccumulator->emplace(id,mtd_digitizer::MTDCellInfo()).first;
 
     // --- Get the simHit energy and convert it from MeV to photo-electrons
     float Npe = 1000.*hit.energyLoss()*LightYield_*LightCollEff_*PDE_;
 
     // --- Get the distance to the center of the detector
     const float dist2center = 0.1f*hit.entryPoint().mag();
 
     // --- Get the simHit time of arrival and correct for the nominal time of flight
     float toa = std::get<2>(hitRefs[ihit]) - dist2center/refSpeed_ + LightCollTime_;
 
     if ( smearLightCollTime_ > 0. )
       toa += CLHEP::RandGaussQ::shoot(hre, 0., smearLightCollTime_);
 
     // --- Accumulate in 15 buckets of 25 ns (9 pre-samples, 1 in-time, 5 post-samples)
     const int itime = std::floor( toa/bxTime_ ) + 9;
 
     if(itime<0 || itime>14) continue;     
 
     // --- Check if the time index is ok and accumulate the energy
     if(itime >= (int)simHitIt->second.hit_info[0].size() ) continue;
     (simHitIt->second).hit_info[0][itime] += Npe;
 
 
     // --- Store the time of the first SimHit
     if( (simHitIt->second).hit_info[1][itime] == 0 )
       (simHitIt->second).hit_info[1][itime] = toa;
 
   } // ihit loop
 
 }