#include <RPixLinearChargeCollectionDrifter.h>

Public Member Functions
std::vector< RPixSignalPoint >	Drift (const std::vector< RPixEnergyDepositUnit > &energy_deposition)

	RPixLinearChargeCollectionDrifter (const edm::ParameterSet &params, uint32_t det_id, const PPSPixelTopology &ppt)

Private Member Functions
double	getSigma_ (double z)

Private Attributes
std::vector< double >	charge_cloud_sigmas_vect_

uint32_t	det_id_

double	det_thickness_

double	GeV_per_electron_

std::vector< RPixSignalPoint >	temp_

int	verbosity_

Detailed Description

Definition at line 9 of file RPixLinearChargeCollectionDrifter.h.

Constructor & Destructor Documentation

◆ RPixLinearChargeCollectionDrifter()

RPixLinearChargeCollectionDrifter::RPixLinearChargeCollectionDrifter	(	const edm::ParameterSet &	params,
		uint32_t	det_id,
		const PPSPixelTopology &	ppt
	)

Definition at line 3 of file RPixLinearChargeCollectionDrifter.cc.

References charge_cloud_sigmas_vect_, det_id_, det_thickness_, PPSPixelTopology::getThickness(), GeV_per_electron_, submitPVValidationJobs::params, and verbosity_.

                                                                                                   {
   verbosity_ = params.getParameter<int>("RPixVerbosity");
 
   GeV_per_electron_ = params.getParameter<double>("RPixGeVPerElectron");
   charge_cloud_sigmas_vect_ = params.getParameter<std::vector<double> >("RPixInterSmearing");
   det_thickness_ = ppt.getThickness();
   det_id_ = det_id;
 }

Member Function Documentation

◆ Drift()

std::vector< RPixSignalPoint > RPixLinearChargeCollectionDrifter::Drift ( const std::vector< RPixEnergyDepositUnit > & energy_deposition )

Definition at line 14 of file RPixLinearChargeCollectionDrifter.cc.

References det_id_, EcalCondDBWriter_cfi::Energy, getSigma_(), GeV_per_electron_, mps_fire::i, PixelTestBeamValidation_cfi::Position, temp_, verbosity_, x, y, and z.

                                                                {
   // convert an energy deposit in a point and in a charge of electrons n=E/3.61 (eV)
   temp_.resize(energy_deposition.size());
   for (unsigned int i = 0; i < energy_deposition.size(); i++) {
     temp_[i].setPosition(LocalPoint(energy_deposition[i].Position().x(), energy_deposition[i].Position().y()));
     temp_[i].setSigma(getSigma_(energy_deposition[i].Position().z()));
     temp_[i].setCharge(energy_deposition[i].Energy() / GeV_per_electron_);
     if (verbosity_ > 1) {
       edm::LogInfo("PPS") << "RPixLinearChargeCollectionDrifter " << det_id_ << " :" << temp_[i].Position() << " "
                           << temp_[i].Sigma() << " " << temp_[i].Charge();
     }
   }
   return temp_;
 }

◆ getSigma_()

double RPixLinearChargeCollectionDrifter::getSigma_ ( double z )

private

Definition at line 29 of file RPixLinearChargeCollectionDrifter.cc.

References reco::ceil(), charge_cloud_sigmas_vect_, det_thickness_, createfilelist::int, and z.

Referenced by Drift().

                                                             {
   if (charge_cloud_sigmas_vect_.size() == 1)
     return charge_cloud_sigmas_vect_[0];
 
   double factor = (z / det_thickness_) * (charge_cloud_sigmas_vect_.size() - 1);
   double lo_i = floor(factor);
   double hi_i = ceil(factor);
   if (lo_i == hi_i) {
     return charge_cloud_sigmas_vect_[(int)factor];
   } else {
     double lo_weight = hi_i - factor;
     double hi_weight = factor - lo_i;
 
     return charge_cloud_sigmas_vect_[(int)lo_i] * lo_weight + charge_cloud_sigmas_vect_[(int)hi_i] * hi_weight;
   }
 }