SiTrivialInduceChargeOnStrips Class Reference

#include <SiTrivialInduceChargeOnStrips.h>

Inheritance diagram for SiTrivialInduceChargeOnStrips:

Public Member Functions
void	induce (SiChargeCollectionDrifter::collection_type collection_points, const StripGeomDetUnit &det, std::vector< double > &localAmplitudes, size_t &recordMinAffectedStrip, size_t &recordMaxAffectedStrip)
	SiTrivialInduceChargeOnStrips (const edm::ParameterSet &conf, double g)
virtual	~SiTrivialInduceChargeOnStrips ()
Public Member Functions inherited from SiInduceChargeOnStrips
virtual	~SiInduceChargeOnStrips ()

Private Member Functions
double	chargeDeposited (size_t strip, size_t Nstrips, double amplitude, double chargeSpread, double chargePosition) const

Static Private Member Functions
static unsigned int	indexOf (const std::string &)
static unsigned int	typeOf (const StripGeomDetUnit &)

Private Attributes
const double	geVperElectron
const double	Nsigma
std::vector< std::vector < double > >	signalCoupling

Static Private Attributes
static const int	Ntypes = 14
static const std::string	type [] = { "IB1", "IB2","OB1","OB2","W1a","W2a","W3a","W1b","W2b","W3b","W4","W5","W6","W7"}

Detailed Description

Definition at line 7 of file SiTrivialInduceChargeOnStrips.h.

Constructor & Destructor Documentation

SiTrivialInduceChargeOnStrips::SiTrivialInduceChargeOnStrips	(	const edm::ParameterSet &	conf,
		double	g
	)

Definition at line 42 of file SiTrivialInduceChargeOnStrips.cc.

References edm::ParameterSet::getParameter(), i, alignBH_cfg::mode, Ntypes, and signalCoupling.

  : Nsigma(3.), geVperElectron(g)  {
  std::string mode = conf.getParameter<bool>("APVpeakmode") ? "Peak" : "Dec";
  for(int i=0; i<Ntypes; i++)
    signalCoupling.push_back(conf.getParameter<std::vector<double> >("CouplingConstant"+mode+type[i]));
}

virtual SiTrivialInduceChargeOnStrips::~SiTrivialInduceChargeOnStrips ( )

inlinevirtual

Definition at line 10 of file SiTrivialInduceChargeOnStrips.h.

{}

Member Function Documentation

double SiTrivialInduceChargeOnStrips::chargeDeposited ( size_t  strip, size_t  Nstrips, double  amplitude, double  chargeSpread, double  chargePosition  ) const

inlineprivate

Definition at line 89 of file SiTrivialInduceChargeOnStrips.cc.

References geVperElectron.

Referenced by induce().

                                                                                                                  {
  double integralUpToStrip = (strip == 0)         ? 0. : ( ROOT::Math::normal_cdf(   strip, chargeSpread, chargePosition) );
  double integralUpToNext  = (strip+1 == Nstrips) ? 1. : ( ROOT::Math::normal_cdf( strip+1, chargeSpread, chargePosition) );
  double percentOfSignal = integralUpToNext - integralUpToStrip;
  
  return percentOfSignal * amplitude / geVperElectron;
}

unsigned int SiTrivialInduceChargeOnStrips::indexOf ( const std::string &  t )

inlinestaticprivate

Definition at line 26 of file SiTrivialInduceChargeOnStrips.cc.

References spr::find(), Ntypes, and type.

Referenced by typeOf().

{ return std::find( type, type + Ntypes, t) - type;}

void SiTrivialInduceChargeOnStrips::induce ( SiChargeCollectionDrifter::collection_type  collection_points, const StripGeomDetUnit &  det, std::vector< double > &  localAmplitudes, size_t &  recordMinAffectedStrip, size_t &  recordMaxAffectedStrip  )

virtual

Implements SiInduceChargeOnStrips.

Definition at line 51 of file SiTrivialInduceChargeOnStrips.cc.

References abs, chargeDeposited(), StripTopology::localPitch(), max(), min, Nsigma, StripTopology::nstrips(), signalCoupling, StripGeomDetUnit::specificTopology(), strip(), StripTopology::strip(), and typeOf().

                                       {

  std::vector<double>& coupling = signalCoupling.at(typeOf(det));
  const StripTopology& topology = dynamic_cast<const StripTopology&>(det.specificTopology());
  size_t Nstrips =  topology.nstrips();

  for (SiChargeCollectionDrifter::collection_type::const_iterator 
     signalpoint = collection_points.begin();  signalpoint != collection_points.end();  signalpoint++ ) {
    
    //In strip coordinates:
    double chargePosition = topology.strip(signalpoint->position());
    double chargeSpread = signalpoint->sigma() / topology.localPitch(signalpoint->position());
    
    size_t fromStrip  = size_t(std::max( 0,          int(std::floor( chargePosition - Nsigma*chargeSpread))));
    size_t untilStrip = size_t(std::min( Nstrips, size_t(std::ceil( chargePosition + Nsigma*chargeSpread) )));
    for (size_t strip = fromStrip;  strip < untilStrip; strip++) {

      double chargeDepositedOnStrip = chargeDeposited( strip, Nstrips, signalpoint->amplitude(), chargeSpread, chargePosition);

      size_t affectedFromStrip  = size_t(std::max( 0, int(strip - coupling.size() + 1)));
      size_t affectedUntilStrip = size_t(std::min( Nstrips, strip + coupling.size())   );  
      for (size_t affectedStrip = affectedFromStrip;  affectedStrip < affectedUntilStrip;  affectedStrip++) {
    localAmplitudes.at( affectedStrip ) += chargeDepositedOnStrip * coupling.at(abs( affectedStrip - strip )) ;
      }

      if( affectedFromStrip  < recordMinAffectedStrip ) recordMinAffectedStrip = affectedFromStrip;
      if( affectedUntilStrip > recordMaxAffectedStrip ) recordMaxAffectedStrip = affectedUntilStrip;
    }
  }
  return;
}

unsigned int SiTrivialInduceChargeOnStrips::typeOf ( const StripGeomDetUnit &  det )

inlinestaticprivate

Definition at line 30 of file SiTrivialInduceChargeOnStrips.cc.

References edm::hlt::Exception, GeomDet::geographicalId(), errorMatrix2Lands_multiChannel::id, indexOf(), TIBDetId::layer(), TIDDetId::ring(), TECDetId::ring(), StripGeomDetUnit::specificType(), GeomDetType::subDetector(), GeomDetEnumerators::TEC, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, and GeomDetEnumerators::TOB.

Referenced by induce().

                                    {
  DetId id = det.geographicalId();
  switch (det.specificType().subDetector()) {
  case GeomDetEnumerators::TIB: {return (TIBDetId(id).layer() < 3) ? indexOf("IB1") : indexOf("IB2");}
  case GeomDetEnumerators::TOB: {return (TOBDetId(id).layer() > 4) ? indexOf("OB1") : indexOf("OB2");}
  case GeomDetEnumerators::TID: {return indexOf("W1a") -1 + TIDDetId(id).ring();} //fragile: relies on ordering of 'type'
  case GeomDetEnumerators::TEC: {return indexOf("W1b") -1 + TECDetId(id).ring();} //fragile: relies on ordering of 'type'
  default: throw cms::Exception("Invalid subdetector") << id();
  }
}

Member Data Documentation

const double SiTrivialInduceChargeOnStrips::geVperElectron

private

Definition at line 30 of file SiTrivialInduceChargeOnStrips.h.

Referenced by chargeDeposited().

const double SiTrivialInduceChargeOnStrips::Nsigma

private

Definition at line 29 of file SiTrivialInduceChargeOnStrips.h.

Referenced by induce().

const int SiTrivialInduceChargeOnStrips::Ntypes = 14

staticprivate

Definition at line 26 of file SiTrivialInduceChargeOnStrips.h.

Referenced by indexOf(), and SiTrivialInduceChargeOnStrips().

std::vector<std::vector<double> > SiTrivialInduceChargeOnStrips::signalCoupling

private

Definition at line 27 of file SiTrivialInduceChargeOnStrips.h.

Referenced by induce(), and SiTrivialInduceChargeOnStrips().

const std::string SiTrivialInduceChargeOnStrips::type = { "IB1", "IB2","OB1","OB2","W1a","W2a","W3a","W1b","W2b","W3b","W4","W5","W6","W7"}

staticprivate

Definition at line 25 of file SiTrivialInduceChargeOnStrips.h.

Referenced by indexOf(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties().