#include <SiHitDigitizer.h>

Public Member Functions
void	processHit (const PSimHit , const StripGeomDetUnit &, GlobalVector, float, std::vector< float > &, size_t &, size_t &, const TrackerTopology tTopo, CLHEP::HepRandomEngine *)

void	setChargeCollectionDrifter (SiChargeCollectionDrifter *cd)

void	setChargeDivider (SiChargeDivider *cd)

void	setInduceChargeOnStrips (SiInduceChargeOnStrips *cd)

void	setParticleDataTable (const ParticleDataTable *pdt)

	SiHitDigitizer (const edm::ParameterSet &conf)

	~SiHitDigitizer ()

Private Types
typedef GloballyPositioned< double >	Frame

Private Member Functions
LocalVector	DriftDirection (const StripGeomDetUnit *_detp, GlobalVector _bfield, float langle)

Private Attributes
const double	chargeMobility

const double	depletionVoltage

std::unique_ptr< SiChargeCollectionDrifter >	theSiChargeCollectionDrifter

std::unique_ptr< SiChargeDivider >	theSiChargeDivider

std::unique_ptr< const SiInduceChargeOnStrips >	theSiInduceChargeOnStrips

Detailed Description

Digitizes the response for a single SimHit.

Definition at line 31 of file SiHitDigitizer.h.

Member Typedef Documentation

◆ Frame

typedef GloballyPositioned<double> SiHitDigitizer::Frame

private

Definition at line 62 of file SiHitDigitizer.h.

Constructor & Destructor Documentation

◆ SiHitDigitizer()

SiHitDigitizer::SiHitDigitizer ( const edm::ParameterSet & conf )

Definition at line 16 of file SiHitDigitizer.cc.

     : depletionVoltage(conf.getParameter<double>("DepletionVoltage")),
       chargeMobility(conf.getParameter<double>("ChargeMobility")),
       theSiChargeDivider(new SiLinearChargeDivider(conf)),
       theSiChargeCollectionDrifter(new SiLinearChargeCollectionDrifter(
           CBOLTZ_over_e_SI * chargeMobility * conf.getParameter<double>("Temperature") *
               (conf.getParameter<bool>("noDiffusion") ? noDiffusionMultiplier : 1.0),
           conf.getParameter<double>("ChargeDistributionRMS"),
           depletionVoltage,
           conf.getParameter<double>("AppliedVoltage"))),
       theSiInduceChargeOnStrips(new SiTrivialInduceChargeOnStrips(conf, conf.getParameter<double>("GevPerElectron"))) {}

◆ ~SiHitDigitizer()

SiHitDigitizer::~SiHitDigitizer ( )

Definition at line 28 of file SiHitDigitizer.cc.

28 {}

Member Function Documentation

◆ DriftDirection()

LocalVector SiHitDigitizer::DriftDirection	(	const StripGeomDetUnit *	_detp,
		GlobalVector	_bfield,
		float	langle
	)

inlineprivate

Definition at line 64 of file SiHitDigitizer.h.

References mkfit::Config::Bfield, GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), GeomDet::surface(), and toLocal().

Referenced by processHit().

                                                                                                 {
     LocalVector Bfield = Frame(_detp->surface().position(), _detp->surface().rotation()).toLocal(_bfield);
     return LocalVector(-langle * Bfield.y(), langle * Bfield.x(), 1.);
   }

◆ processHit()

void SiHitDigitizer::processHit	(	const PSimHit *	hit,
		const StripGeomDetUnit &	det,
		GlobalVector	bfield,
		float	langle,
		std::vector< float > &	locAmpl,
		size_t &	firstChannelWithSignal,
		size_t &	lastChannelWithSignal,
		const TrackerTopology *	tTopo,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 30 of file SiHitDigitizer.cc.

References Surface::bounds(), chargeMobility, depletionVoltage, DriftDirection(), GeomDet::specificSurface(), theSiChargeCollectionDrifter, theSiChargeDivider, theSiInduceChargeOnStrips, and Bounds::thickness().

Referenced by DigiSimLinkAlgorithm::run().

                                                               {
   // Compute the drift direction for this det
   double moduleThickness = det.specificSurface().bounds().thickness();  // active detector thicness
   double timeNormalisation = (moduleThickness * moduleThickness) / (2. * depletionVoltage * chargeMobility);
   LocalVector driftDir = DriftDirection(&det, bfield, langle);
 
   // Fully process one SimHit
   theSiInduceChargeOnStrips->induce(
       theSiChargeCollectionDrifter->drift(theSiChargeDivider->divide(hit, driftDir, moduleThickness, det, engine),
                                           driftDir,
                                           moduleThickness,
                                           timeNormalisation),
       det,
       locAmpl,
       firstChannelWithSignal,
       lastChannelWithSignal,
       tTopo);
 }