#include <StripCPEfromTrackAngle.h>

Inheritance diagram for StripCPEfromTrackAngle:

Public Member Functions
float	legacyStripErrorSquared (const unsigned N, const float uProj) const

StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &, const GeomDetUnit &, const LocalTrajectoryParameters &) const

	StripCPEfromTrackAngle (edm::ParameterSet &conf, const MagneticField &mag, const TrackerGeometry &geom, const SiStripLorentzAngle &lorentz, const SiStripBackPlaneCorrection &backPlaneCorrection, const SiStripConfObject &confObj, const SiStripLatency &latency)

float	stripErrorSquared (const unsigned N, const float uProj, const SiStripDetId::SubDetector loc) const

Public Member Functions inherited from StripCPE
LocalVector	driftDirection (const StripGeomDetUnit *det) const

StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const

	StripCPE (edm::ParameterSet &conf, const MagneticField &, const TrackerGeometry &, const SiStripLorentzAngle &, const SiStripBackPlaneCorrection &, const SiStripConfObject &, const SiStripLatency &)

Public Member Functions inherited from StripClusterParameterEstimator
virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd) const

virtual VLocalValues	localParametersV (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual	~StripClusterParameterEstimator ()

Private Attributes
const float	maxChgOneMIP

float	mHC_P [4][2]

std::map < SiStripDetId::SubDetector, float >	mHC_P0

std::map < SiStripDetId::SubDetector, float >	mHC_P1

float	mLC_P [3]

const bool	useLegacyError

Additional Inherited Members
Public Types inherited from StripClusterParameterEstimator
typedef std::pair< LocalPoint, LocalError >	LocalValues

typedef std::vector< LocalValues >	VLocalValues

Protected Member Functions inherited from StripCPE
Param const &	param (const GeomDetUnit &det) const

Protected Attributes inherited from StripCPE
const SiStripBackPlaneCorrection &	BackPlaneCorrectionMap_

const TrackerGeometry &	geom_

const SiStripLorentzAngle &	LorentzAngleMap_

const MagneticField &	magfield_

const bool	peakMode_

std::vector< float >	xtalk1

std::vector< float >	xtalk2

Detailed Description

Definition at line 6 of file StripCPEfromTrackAngle.h.

Constructor & Destructor Documentation

StripCPEfromTrackAngle::StripCPEfromTrackAngle	(	edm::ParameterSet &	conf,
		const MagneticField &	mag,
		const TrackerGeometry &	geom,
		const SiStripLorentzAngle &	lorentz,
		const SiStripBackPlaneCorrection &	backPlaneCorrection,
		const SiStripConfObject &	confObj,
		const SiStripLatency &	latency
	)

inline

Definition at line 34 of file StripCPEfromTrackAngle.h.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), mHC_P, mLC_P, SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, and SiStripDetId::TOB.

   : StripCPE(conf, mag, geom, lorentz, backPlaneCorrection, confObj, latency )
   , useLegacyError(conf.existsAs<bool>("useLegacyError") ? conf.getParameter<bool>("useLegacyError") : true)
   , maxChgOneMIP(conf.existsAs<float>("maxChgOneMIP") ? conf.getParameter<double>("maxChgOneMIP") : -6000.)
   {
     mLC_P[0] = conf.existsAs<double>("mLC_P0") ? conf.getParameter<double>("mLC_P0") : -.326;
     mLC_P[1] = conf.existsAs<double>("mLC_P1") ? conf.getParameter<double>("mLC_P1") :  .618;
     mLC_P[2] = conf.existsAs<double>("mLC_P2") ? conf.getParameter<double>("mLC_P2") :  .300;
 
     mHC_P[SiStripDetId::TIB - 3][0] = conf.existsAs<double>("mTIB_P0") ? conf.getParameter<double>("mTIB_P0") : -.742  ;
     mHC_P[SiStripDetId::TIB - 3][1] = conf.existsAs<double>("mTIB_P1") ? conf.getParameter<double>("mTIB_P1") :  .202  ;
     mHC_P[SiStripDetId::TID - 3][0] = conf.existsAs<double>("mTID_P0") ? conf.getParameter<double>("mTID_P0") : -1.026 ;
     mHC_P[SiStripDetId::TID - 3][1] = conf.existsAs<double>("mTID_P1") ? conf.getParameter<double>("mTID_P1") :  .253  ;
     mHC_P[SiStripDetId::TOB - 3][0] = conf.existsAs<double>("mTOB_P0") ? conf.getParameter<double>("mTOB_P0") : -1.427 ;
     mHC_P[SiStripDetId::TOB - 3][1] = conf.existsAs<double>("mTOB_P1") ? conf.getParameter<double>("mTOB_P1") :  .433  ;
     mHC_P[SiStripDetId::TEC - 3][0] = conf.existsAs<double>("mTEC_P0") ? conf.getParameter<double>("mTEC_P0") : -1.885 ;
     mHC_P[SiStripDetId::TEC - 3][1] = conf.existsAs<double>("mTEC_P1") ? conf.getParameter<double>("mTEC_P1") :  .471  ;  
   }

Member Function Documentation

float StripCPEfromTrackAngle::legacyStripErrorSquared	(	const unsigned	N,
		const float	uProj
	)		const

Definition at line 16 of file StripCPEfromTrackAngle.cc.

References constexpr, f, and myMath::fast_expf().

Referenced by localParameters().

                                                                                                {
   if( (float(N)-uProj) > 3.5f )
     return float(N*N)/12.f;
   else {
     static constexpr float P1=-0.339;
     static constexpr float P2=0.90;
     static constexpr float P3=0.279;
     const float uerr = P1*uProj*vdt::fast_expf(-uProj*P2)+P3;
     return uerr*uerr;
   }
 }

StripClusterParameterEstimator::LocalValues StripCPEfromTrackAngle::localParameters	(	const SiStripCluster &	cluster,
		const GeomDetUnit &	det,
		const LocalTrajectoryParameters &	ltp
	)		const

virtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 29 of file StripCPEfromTrackAngle.cc.

References funct::abs(), SiStripCluster::amplitudes(), StripCPE::Param::backplanecorrection, SiStripCluster::barycenter(), siStripClusterTools::chargePerCM(), StripCPE::Param::coveredStrips(), StripCPE::Param::drift, f, GeomDet::geographicalId(), SiStripCluster::isMerged(), legacyStripErrorSquared(), StripTopology::localError(), StripTopology::localPosition(), maxChgOneMIP, StripCPE::Param::maxLength, LocalTrajectoryParameters::momentum(), N, AlCaHLTBitMon_ParallelJobs::p, StripCPE::param(), LocalTrajectoryParameters::position(), stripErrorSquared(), SiStripDetId::subDetector(), StripCPE::Param::thickness, StripCPE::Param::topology, useLegacyError, and LocalTrajectoryParameters::vector().

                                                                                                                     {
   
   StripCPE::Param const & p = param(det);
   SiStripDetId ssdid = SiStripDetId( det.geographicalId() );  
  
   LocalVector track = ltp.momentum();
   track *= 
     (track.z()<0) ?  std::abs(p.thickness/track.z()) : 
     (track.z()>0) ? -std::abs(p.thickness/track.z()) :  
                          p.maxLength/track.mag() ;
 
   const unsigned N = cluster.amplitudes().size();
   const float fullProjection = p.coveredStrips( track+p.drift, ltp.position());
   float uerr2;
   if (useLegacyError) {
     uerr2 = legacyStripErrorSquared(N,std::abs(fullProjection));
   } else if (maxChgOneMIP < 0.0) {
     uerr2 = cluster.isMerged() ? legacyStripErrorSquared(N,std::abs(fullProjection)) : stripErrorSquared( N, std::abs(fullProjection),ssdid.subDetector() );
   } else {
     float dQdx = siStripClusterTools::chargePerCM(ssdid, cluster, ltp);
     uerr2 = dQdx > maxChgOneMIP ? legacyStripErrorSquared(N,std::abs(fullProjection)) : stripErrorSquared( N, std::abs(fullProjection),ssdid.subDetector() );
   }
   const float strip = cluster.barycenter() -  0.5f*(1.f-p.backplanecorrection) * fullProjection
     + 0.5f*p.coveredStrips(track, ltp.position());
 
   return std::make_pair( p.topology->localPosition(strip, ltp.vector()),
              p.topology->localError(strip, uerr2, ltp.vector()) );
 }

float StripCPEfromTrackAngle::stripErrorSquared	(	const unsigned	N,
		const float	uProj,
		const SiStripDetId::SubDetector	loc
	)		const

Definition at line 7 of file StripCPEfromTrackAngle.cc.

References Exception, myMath::fast_expf(), mHC_P, mLC_P, SiStripDetId::UNKNOWN, and x.

Referenced by localParameters().

                                                                                                                              {
   if( loc == SiStripDetId::UNKNOWN)
     throw cms::Exception("StripCPEfromTrackAngle::stripErrorSquared", "Incompatible sub-detector.");
 
   auto fun = [&] (float x)  -> float { return mLC_P[0]*x*vdt::fast_expf(-x*mLC_P[1])+mLC_P[2];};
   auto uerr = (N <= 4) ?  fun(uProj) :  mHC_P[loc-3][0]+float(N)*mHC_P[loc-3][1];
   return uerr*uerr;
 }