#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 override

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 31 of file StripCPEfromTrackAngle.h.

using StripCPEfromTrackAngle::AlgoParam = StripCPE::AlgoParam

Definition at line 30 of file StripCPEfromTrackAngle.h.

using StripCPEfromTrackAngle::ALocalValues = StripClusterParameterEstimator::ALocalValues

Definition at line 32 of file StripCPEfromTrackAngle.h.

Member Enumeration Documentation

enum StripCPEfromTrackAngle::Algo

strongprivate

Enumerator
legacy
mergeCK
chargeCK

Definition at line 25 of file StripCPEfromTrackAngle.h.

25 { 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 validate-o2o-wbm::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 52 of file StripCPEfromTrackAngle.cc.

References chargeCK, siStripClusterTools::chargePerCM(), alignCSCRings::corr, PVValHelper::fill(), StripCPE::getAlgoParam(), mps_fire::i, legacy, legacyStripErrorSquared(), m_algo, maxChgOneMIP, mergeCK, N, AlCaHLTBitMon_ParallelJobs::p, DynArray< T >::size(), digitizers_cfi::strip, 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 94 of file StripCPEfromTrackAngle.cc.

References StripCPE::AlgoParam::afullProjection, SiStripCluster::amplitudes(), SiStripCluster::barycenter(), chargeCK, siStripClusterTools::chargePerCM(), StripCPE::AlgoParam::corr, alignCSCRings::corr, SiStripCluster::isMerged(), legacy, legacyStripErrorSquared(), StripCPE::AlgoParam::loc, StripCPE::AlgoParam::ltp, m_algo, maxChgOneMIP, mergeCK, N, StripCPE::AlgoParam::p, AlCaHLTBitMon_ParallelJobs::p, SiStripCluster::size(), digitizers_cfi::strip, 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 125 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 32 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 27 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

const float StripCPEfromTrackAngle::maxChgOneMIP

private

Definition at line 23 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

float StripCPEfromTrackAngle::mHC_P[4][2]

private

Definition at line 12 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

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

private

Definition at line 15 of file StripCPEfromTrackAngle.h.

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

private

Definition at line 16 of file StripCPEfromTrackAngle.h.

float StripCPEfromTrackAngle::mLC_P[3]

private

Definition at line 11 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

const bool StripCPEfromTrackAngle::useLegacyError

private

Definition at line 19 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