#include <StripCPEfromTrackAngle.h>

Inheritance diagram for StripCPEfromTrackAngle:

Public Types
using	AClusters = StripClusterParameterEstimator::AClusters

using	AlgoParam = StripCPE::AlgoParam

using	ALocalValues = StripClusterParameterEstimator::ALocalValues

Public Types inherited from StripClusterParameterEstimator
using	AClusters = DynArray< SiStripCluster const * >

using	ALocalValues = DynArray< LocalValues >

using	LocalValues = std::pair< LocalPoint, LocalError >

typedef std::vector< LocalValues >	VLocalValues

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

void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const override

StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, AlgoParam const &ap) const override

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

	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

AlgoParam	getAlgoParam (const GeomDetUnit &det, const LocalTrajectoryParameters &ltp) const

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

	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 Types
enum	Algo { Algo::legacy, Algo::mergeCK, Algo::chargeCK }

Private Attributes
Algo	m_algo

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
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.

Member Typedef Documentation

using StripCPEfromTrackAngle::AClusters = StripClusterParameterEstimator::AClusters

Definition at line 33 of file StripCPEfromTrackAngle.h.

using StripCPEfromTrackAngle::AlgoParam = StripCPE::AlgoParam

Definition at line 32 of file StripCPEfromTrackAngle.h.

using StripCPEfromTrackAngle::ALocalValues = StripClusterParameterEstimator::ALocalValues

Definition at line 34 of file StripCPEfromTrackAngle.h.

Member Enumeration Documentation

enum StripCPEfromTrackAngle::Algo

strongprivate

Enumerator
legacy
mergeCK
chargeCK

Definition at line 27 of file StripCPEfromTrackAngle.h.

27 { legacy, mergeCK, chargeCK };

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
	)

Definition at line 7 of file StripCPEfromTrackAngle.cc.

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.),
     m_algo(useLegacyError ? Algo::legacy : ( maxChgOneMIP<0 ? Algo::mergeCK : Algo::chargeCK))
   {
     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 40 of file StripCPEfromTrackAngle.cc.

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

Referenced by localParameters().

                                                                                                {
   if unlikely( (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;
   }
 }

void StripCPEfromTrackAngle::localParameters	(	AClusters const &	clusters,
		ALocalValues &	retValues,
		const GeomDetUnit &	gd,
		const LocalTrajectoryParameters &	ltp
	)		const

overridevirtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 54 of file StripCPEfromTrackAngle.cc.

References chargeCK, siStripClusterTools::chargePerCM(), corr, lumiContext::fill, StripCPE::getAlgoParam(), i, legacy, legacyStripErrorSquared(), m_algo, maxChgOneMIP, mergeCK, N, AlCaHLTBitMon_ParallelJobs::p, DynArray< T >::size(), stripErrorSquared(), and LocalTrajectoryParameters::vector().

Referenced by localParameters().

                                                                                                                                                                       {
 
 
   auto const & par = getAlgoParam(det,ltp);
   auto const & p = par.p;
   auto loc = par.loc;
   auto corr = par.corr;
   auto afp = par.afullProjection;
 
   auto fill = [&](unsigned int i, float uerr2) {
     const float strip = clusters[i]->barycenter() + corr;
     retValues[i].first  = p.topology->localPosition(strip, ltp.vector());
     retValues[i].second = p.topology->localError(strip, uerr2, ltp.vector());
   };
 
   
   switch (m_algo) {
   case Algo::chargeCK :
     for (auto i=0U; i< clusters.size(); ++i) {
       auto dQdx = siStripClusterTools::chargePerCM(*clusters[i], ltp, p.invThickness);
       auto N = clusters[i]->amplitudes().size();
       auto uerr2 = dQdx > maxChgOneMIP ? legacyStripErrorSquared(N,afp) : stripErrorSquared( N, afp,loc );
       fill(i, uerr2);
     }
     break;
   case Algo::legacy :
     for (auto i=0U; i< clusters.size(); ++i)  {
       auto N = clusters[i]->amplitudes().size();
       auto uerr2 = legacyStripErrorSquared(N,afp);
       fill(i, uerr2);
     }
     break;
   case Algo::mergeCK :
     for (auto i=0U; i< clusters.size(); ++i) {
       auto N = clusters[i]->amplitudes().size();
       auto uerr2 = clusters[i]->isMerged() ? legacyStripErrorSquared(N,afp) : stripErrorSquared( N,afp,loc );
       fill(i, uerr2);
     }
     break;
   }
 
 
 }

StripClusterParameterEstimator::LocalValues StripCPEfromTrackAngle::localParameters	(	const SiStripCluster &	cl,
		AlgoParam const &	ap
	)		const

overridevirtual

Reimplemented from StripCPE.

Definition at line 99 of file StripCPEfromTrackAngle.cc.

References StripCPE::AlgoParam::afullProjection, SiStripCluster::amplitudes(), SiStripCluster::barycenter(), chargeCK, siStripClusterTools::chargePerCM(), corr, StripCPE::AlgoParam::corr, SiStripCluster::isMerged(), legacy, legacyStripErrorSquared(), StripCPE::AlgoParam::loc, StripCPE::AlgoParam::ltp, m_algo, maxChgOneMIP, mergeCK, N, StripCPE::AlgoParam::p, AlCaHLTBitMon_ParallelJobs::p, and stripErrorSquared().

                                                                                                    {
   auto const & p = par.p;
   auto const & ltp = par.ltp;
   auto loc = par.loc;
   auto corr = par.corr;
   auto afp = par.afullProjection;
 
   float uerr2=0;
   
   auto N = cluster.amplitudes().size();
 
   switch (m_algo) {
   case Algo::chargeCK :
     {
       auto dQdx = siStripClusterTools::chargePerCM(cluster, ltp, p.invThickness);
       uerr2 = dQdx > maxChgOneMIP ? legacyStripErrorSquared(N,afp) : stripErrorSquared( N, afp,loc );
     }
     break;
   case Algo::legacy :
     uerr2 = legacyStripErrorSquared(N,afp);
     break;
   case Algo::mergeCK :
     uerr2 = cluster.isMerged() ? legacyStripErrorSquared(N,afp) : stripErrorSquared( N, afp,loc );
     break;
   }
   
   const float strip = cluster.barycenter() + corr;
   
   return std::make_pair( p.topology->localPosition(strip, ltp.vector()),
              p.topology->localError(strip, uerr2, ltp.vector()) );
 }

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

overridevirtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 134 of file StripCPEfromTrackAngle.cc.

References StripCPE::getAlgoParam(), and localParameters().

                                                                                                                                           {
   
   auto const & par = getAlgoParam(det,ltp);
   return localParameters(cluster,par);
 }

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

Definition at line 33 of file StripCPEfromTrackAngle.cc.

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

Referenced by localParameters().

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

Member Data Documentation

Algo StripCPEfromTrackAngle::m_algo

private

Definition at line 29 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

const float StripCPEfromTrackAngle::maxChgOneMIP

private

Definition at line 25 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

float StripCPEfromTrackAngle::mHC_P[4][2]

private

Definition at line 14 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

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

private

Definition at line 17 of file StripCPEfromTrackAngle.h.

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

private

Definition at line 18 of file StripCPEfromTrackAngle.h.

float StripCPEfromTrackAngle::mLC_P[3]

private

Definition at line 13 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

const bool StripCPEfromTrackAngle::useLegacyError

private

Definition at line 21 of file StripCPEfromTrackAngle.h.

Public Types

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation