#include <RPLinearChargeCollectionDrifter.h>

Public Member Functions
simromanpot::charge_induced_on_surface	Drift (const simromanpot::energy_path_distribution &energy_deposition)

	RPLinearChargeCollectionDrifter (const edm::ParameterSet &params, RPDetId det_id)

Private Member Functions
double	getSigma (double z)

Private Attributes
std::vector< double >	charge_cloud_sigmas_vect_

RPDetId	det_id_

double	det_thickness_

double	GeV_per_electron_

int	verbosity_

Detailed Description

Definition at line 9 of file RPLinearChargeCollectionDrifter.h.

Constructor & Destructor Documentation

RPLinearChargeCollectionDrifter::RPLinearChargeCollectionDrifter	(	const edm::ParameterSet &	params,
		RPDetId	det_id
	)

Definition at line 7 of file RPLinearChargeCollectionDrifter.cc.

References charge_cloud_sigmas_vect_, det_id_, det_thickness_, RPTopology::DetThickness(), edm::ParameterSet::getParameter(), GeV_per_electron_, and verbosity_.

                                                                                                               {
   verbosity_ = params.getParameter<int>("RPVerbosity");
   GeV_per_electron_ = params.getParameter<double>("RPGeVPerElectron");
   charge_cloud_sigmas_vect_ = params.getParameter<std::vector<double> >("RPInterStripSmearing");
   det_thickness_ = RPTopology().DetThickness();
   det_id_ = det_id;
 }

Member Function Documentation

simromanpot::charge_induced_on_surface RPLinearChargeCollectionDrifter::Drift ( const simromanpot::energy_path_distribution & energy_deposition )

Definition at line 15 of file RPLinearChargeCollectionDrifter.cc.

References det_id_, getSigma(), GeV_per_electron_, mps_fire::i, PixelTestBeamValidation_cfi::Position, verbosity_, x, y, and z.

                                                                   {
   simromanpot::charge_induced_on_surface temp_;
   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()));  //befor charge_cloud_sigma_ used, now a vector of sigmas;
     temp_[i].setCharge(energy_deposition[i].Energy() / GeV_per_electron_);
     if (verbosity_) {
       edm::LogInfo("RPLinearChargeCollectionDrifter")
           << det_id_ << " :" << temp_[i].Position() << " " << temp_[i].Sigma() << " " << temp_[i].Charge() << "\n";
     }
   }
   return temp_;
 }

double RPLinearChargeCollectionDrifter::getSigma ( double z )

private

Definition at line 31 of file RPLinearChargeCollectionDrifter.cc.

References reco::ceil(), charge_cloud_sigmas_vect_, and det_thickness_.

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;
   }
 }