#include <StripCPEgeometric.h>

Inheritance diagram for StripCPEgeometric:

Classes
class	WrappedCluster

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

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

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

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

virtual LocalValues	localParameters (const SiStripCluster &, const GeomDetUnit &) const

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

	StripCPEgeometric (edm::ParameterSet &conf, const MagneticField &mag, const TrackerGeometry &geom, const SiStripLorentzAngle &LorentzAngle, const SiStripBackPlaneCorrection &BackPlaneCorrection, const SiStripConfObject &confObj, const SiStripLatency &latency)

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

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

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const TrajectoryStateOnSurface &tsos) const

virtual LocalValues	localParameters (const SiStripCluster &, const GeomDetUnit &) const

virtual LocalValues	localParameters (const SiStripCluster &cluster, const GeomDetUnit &gd, const LocalTrajectoryParameters &) const

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

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

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

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

Public Member Functions inherited from StripClusterParameterEstimator
virtual void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const

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 Member Functions
bool	ambiguousSize (const WrappedCluster &, const stats_t< float > &) const

stats_t< float >	geometric_position (const WrappedCluster &, const stats_t< float > &) const

stats_t< float >	offset_from_firstStrip (const std::vector< stats_t< float > > &, const stats_t< float > &) const

bool	useNMinusOne (const WrappedCluster &, const stats_t< float > &) const

bool	useNPlusOne (const WrappedCluster &, const stats_t< float > &) const

Private Attributes
const float	maybe_noise_threshold

const float	minimum_uncertainty_squared

const float	noise_threshold

const float	scaling_squared

const float	tan_diffusion_angle

const float	thickness_rel_err2

Additional Inherited Members
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

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

Constructor & Destructor Documentation

◆ StripCPEgeometric()

StripCPEgeometric::StripCPEgeometric	(	edm::ParameterSet &	conf,
		const MagneticField &	mag,
		const TrackerGeometry &	geom,
		const SiStripLorentzAngle &	LorentzAngle,
		const SiStripBackPlaneCorrection &	BackPlaneCorrection,
		const SiStripConfObject &	confObj,
		const SiStripLatency &	latency
	)

inline

Definition at line 14 of file StripCPEgeometric.h.

       : StripCPE(conf, mag, geom, LorentzAngle, BackPlaneCorrection, confObj, latency),
         tan_diffusion_angle(conf.getParameter<double>("TanDiffusionAngle")),
         thickness_rel_err2(pow(conf.getParameter<double>("ThicknessRelativeUncertainty"), 2)),
         noise_threshold(conf.getParameter<double>("NoiseThreshold")),
         maybe_noise_threshold(conf.getParameter<double>("MaybeNoiseThreshold")),
         scaling_squared(pow(conf.getParameter<double>("UncertaintyScaling"), 2)),
         minimum_uncertainty_squared(pow(conf.getParameter<double>("MinimumUncertainty"), 2)) {}

Member Function Documentation

◆ ambiguousSize()

bool StripCPEgeometric::ambiguousSize	(	const WrappedCluster &	wc,
		const stats_t< float > &	proj
	)		const

inlineprivate

Definition at line 96 of file StripCPEgeometric.cc.

References maybe_noise_threshold, StripCPEgeometric::WrappedCluster::N, amptDefault_cfi::proj, and StripCPEgeometric::WrappedCluster::smallerEdgeStrip().

Referenced by offset_from_firstStrip().

                                                                                                        {
   return proj() < wc.N - 1 && wc.smallerEdgeStrip()() > 0 && wc.smallerEdgeStrip().sigmaFrom(0) < maybe_noise_threshold;
 }

◆ geometric_position()

stats_t< float > StripCPEgeometric::geometric_position	(	const WrappedCluster &	wc,
		const stats_t< float > &	proj
	)		const

private

Definition at line 63 of file StripCPEgeometric.cc.

References StripCPEgeometric::WrappedCluster::eta(), StripCPEgeometric::WrappedCluster::middle(), StripCPEgeometric::WrappedCluster::N, funct::pow(), amptDefault_cfi::proj, scaling_squared, and x.

Referenced by offset_from_firstStrip().

                                                                                        {
   const stats_t<float> x = wc.middle() + 0.5 * proj * wc.eta();
   return wc.N == 1 ? stats_t<float>(x(), pow(1 - 0.82 * proj(), 2) / 12)
                    : stats_t<float>(x(), scaling_squared * x.error2());
 }

◆ localParameters() [1/7]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 40 of file StripClusterParameterEstimator.h.

                                                                               {
     return localParameters(cluster, gd);
   }

◆ localParameters() [2/7]

virtual StripClusterParameterEstimator::LocalValues StripCPE::localParameters

inline

Definition at line 52 of file StripCPE.h.

                                                                                                  {
     return std::make_pair(LocalPoint(), LocalError());
   }

◆ localParameters() [3/7]

virtual void StripClusterParameterEstimator::localParameters

inline

Definition at line 32 of file StripClusterParameterEstimator.h.

35 {}

◆ localParameters() [4/7]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 37 of file StripClusterParameterEstimator.h.

                                                                                        {
     return std::make_pair(LocalPoint(), LocalError());
   }

◆ localParameters() [5/7]

StripClusterParameterEstimator::LocalValues StripCPE::localParameters

override

Definition at line 62 of file StripCPE.cc.

                                                                                                     {
   StripCPE::Param const& p = param(det);
   const float barycenter = cluster.barycenter();
   const float fullProjection =
       p.coveredStrips(p.drift + LocalVector(0, 0, -p.thickness), p.topology->localPosition(barycenter));
   const float strip = barycenter - 0.5f * (1.f - p.backplanecorrection) * fullProjection;
 
   return std::make_pair(p.topology->localPosition(strip), p.topology->localError(strip, 1.f / 12.f));
 }

◆ localParameters() [6/7]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 45 of file StripClusterParameterEstimator.h.

                                                                                   {
     return localParameters(cluster, gd, tsos.localParameters());
   }

◆ localParameters() [7/7]

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

overridevirtual

Reimplemented from StripClusterParameterEstimator.

Definition at line 6 of file StripCPEgeometric.cc.

References SiStripCluster::firstStrip(), stats_t< T >::from_relative_uncertainty2(), SiStripPI::max, minimum_uncertainty_squared, LocalTrajectoryParameters::momentum(), offset_from_firstStrip(), AlCaHLTBitMon_ParallelJobs::p, StripCPE::param(), LocalTrajectoryParameters::position(), nano_mu_digi_cff::strip, tan_diffusion_angle, thickness_rel_err2, HLT_2023v12_cff::track, reco::InverseCrosstalkMatrix::unfold(), LocalTrajectoryParameters::vector(), and StripCPE::xtalk1.

                                                                                                        {
   StripCPE::Param const& p = param(det);
 
   const LocalPoint& pos = ltp.position();
   LocalVector track = ltp.momentum();
   track *= (track.z() < 0)   ? fabs(p.thickness / track.z())
            : (track.z() > 0) ? -fabs(p.thickness / track.z())
                              : p.maxLength / track.mag();
 
   const float fullProjection = p.coveredStrips(track + p.drift, pos);
   stats_t<float> projection;
   {
     const float absProj = fabs(fullProjection);
     const float minProj = 2 * p.thickness * tan_diffusion_angle / p.topology->localPitch(pos);
     const float projection_rel_err2 = thickness_rel_err2 + p.pitch_rel_err2;
     projection = stats_t<float>::from_relative_uncertainty2(std::max(absProj, minProj), projection_rel_err2);
   }
 
   const std::vector<stats_t<float> > Q =
       reco::InverseCrosstalkMatrix::unfold(cluster, xtalk1[static_cast<int>(p.moduleGeom)]);
   const stats_t<float> strip = cluster.firstStrip() + offset_from_firstStrip(Q, projection);
 
   const float corrected =
       strip() - 0.5 * (1 - p.backplanecorrection) * fullProjection + 0.5 * p.coveredStrips(track, ltp.position());
   const float error2 = std::max(strip.error2(), minimum_uncertainty_squared);
 
   return std::make_pair(p.topology->localPosition(corrected, ltp.vector()),
                         p.topology->localError(corrected, error2, ltp.vector()));
 }

◆ offset_from_firstStrip()

stats_t< float > StripCPEgeometric::offset_from_firstStrip	(	const std::vector< stats_t< float > > &	Q,
		const stats_t< float > &	proj
	)		const

private

Definition at line 37 of file StripCPEgeometric.cc.

References StripCPEgeometric::WrappedCluster::addSuppressedEdgeStrip(), ambiguousSize(), StripCPEgeometric::WrappedCluster::centroid(), StripCPEgeometric::WrappedCluster::deformed(), StripCPEgeometric::WrappedCluster::dropSmallerEdgeStrip(), stats_t< T >::error2(), nano_mu_digi_cff::float, geometric_position(), SiStripPI::max, StripCPEgeometric::WrappedCluster::maxProjection(), StripCPEgeometric::WrappedCluster::middle(), StripCPEgeometric::WrappedCluster::N, funct::pow(), amptDefault_cfi::proj, useNMinusOne(), and useNPlusOne().

Referenced by localParameters().

                                                                                            {
   WrappedCluster wc(Q);
   if (useNPlusOne(wc, proj))
     wc.addSuppressedEdgeStrip();
   else
     while (useNMinusOne(wc, proj))
       wc.dropSmallerEdgeStrip();
 
   if (proj() < wc.N - 2)
     return stats_t<float>(wc.middle(), pow(wc.N - proj(), 2) / 12.);
   if (wc.deformed())
     return stats_t<float>(wc.centroid()(), 1 / 12.);
   if (proj > wc.maxProjection())
     return stats_t<float>(wc.centroid()(), 1 / 12.);
 
   if (ambiguousSize(wc, proj)) {
     const stats_t<float> probably = geometric_position(wc, proj);
     wc.dropSmallerEdgeStrip();
     const stats_t<float> maybe = geometric_position(wc, proj);
     return stats_t<float>(probably(),
                           std::max(probably.error2(), float(maybe.error2() + pow(probably() - maybe(), 2) / 12)));
   }
   return geometric_position(wc, proj);
 }

◆ useNMinusOne()

bool StripCPEgeometric::useNMinusOne	(	const WrappedCluster &	wc,
		const stats_t< float > &	proj
	)		const

inlineprivate

Definition at line 74 of file StripCPEgeometric.cc.

References StripCPEgeometric::WrappedCluster::dropSmallerEdgeStrip(), StripCPEgeometric::WrappedCluster::eta(), StripCPEgeometric::WrappedCluster::maxProjection(), StripCPEgeometric::WrappedCluster::N, noise_threshold, amptDefault_cfi::proj, StripCPEgeometric::WrappedCluster::sign(), and StripCPEgeometric::WrappedCluster::smallerEdgeStrip().

Referenced by offset_from_firstStrip().

                                                                                                       {
   if (proj() > wc.N - 1)
     return false;
   if (wc.smallerEdgeStrip() < 0)
     return true;
   if (proj() < wc.N - 3)
     return false;
   if (proj() < wc.N - 2)
     return true;
   if (wc.eta().sigmaFrom(0) < 3)
     return false;
 
   WrappedCluster wcTest(wc);
   wcTest.dropSmallerEdgeStrip();
   if (proj >= wcTest.maxProjection())
     return false;
   if (wc.sign() * wc.eta()() > 1. / (wc.N - 1))
     return true;
 
   return wc.smallerEdgeStrip().sigmaFrom(0) < noise_threshold;
 }

◆ useNPlusOne()

bool StripCPEgeometric::useNPlusOne	(	const WrappedCluster &	wc,
		const stats_t< float > &	proj
	)		const

inlineprivate

Definition at line 70 of file StripCPEgeometric.cc.

References StripCPEgeometric::WrappedCluster::maxProjection(), StripCPEgeometric::WrappedCluster::N, and amptDefault_cfi::proj.

Referenced by offset_from_firstStrip().

                                                                                                      {
   return wc.maxProjection() < proj && proj() < wc.N + 1;
 }

Member Data Documentation

◆ maybe_noise_threshold

const float StripCPEgeometric::maybe_noise_threshold

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by ambiguousSize().

◆ minimum_uncertainty_squared

const float StripCPEgeometric::minimum_uncertainty_squared

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by localParameters().

◆ noise_threshold

const float StripCPEgeometric::noise_threshold

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by useNMinusOne().

◆ scaling_squared

const float StripCPEgeometric::scaling_squared

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by geometric_position().

◆ tan_diffusion_angle

const float StripCPEgeometric::tan_diffusion_angle

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by localParameters().

◆ thickness_rel_err2

const float StripCPEgeometric::thickness_rel_err2

private

Definition at line 30 of file StripCPEgeometric.h.

Referenced by localParameters().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ StripCPEgeometric()

Member Function Documentation

◆ ambiguousSize()

◆ geometric_position()

◆ localParameters() [1/7]

◆ localParameters() [2/7]

◆ localParameters() [3/7]

◆ localParameters() [4/7]

◆ localParameters() [5/7]

◆ localParameters() [6/7]

◆ localParameters() [7/7]

◆ offset_from_firstStrip()

◆ useNMinusOne()

◆ useNPlusOne()

Member Data Documentation

◆ maybe_noise_threshold

◆ minimum_uncertainty_squared

◆ noise_threshold

◆ scaling_squared

◆ tan_diffusion_angle

◆ thickness_rel_err2