#include <StripCPEgeometric.h>

Inheritance diagram for StripCPEgeometric:

Classes
class	WrappedCluster

Public Member Functions
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

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

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

inlineprivate

Definition at line 95 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;
 }

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

private

Definition at line 62 of file StripCPEgeometric.cc.

References stats_t< T >::error2(), 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());
 }

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 StripCPE::Param::backplanecorrection, StripCPE::Param::coveredStrips(), StripCPE::Param::drift, stats_t< T >::error2(), SiStripCluster::firstStrip(), stats_t< T >::from_relative_uncertainty2(), StripTopology::localError(), StripTopology::localPitch(), StripTopology::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), SiStripPI::max, StripCPE::Param::maxLength, minimum_uncertainty_squared, StripCPE::Param::moduleGeom, LocalTrajectoryParameters::momentum(), offset_from_firstStrip(), AlCaHLTBitMon_ParallelJobs::p, StripCPE::param(), StripCPE::Param::pitch_rel_err2, LocalTrajectoryParameters::position(), class-composition::Q, digitizers_cfi::strip, tan_diffusion_angle, StripCPE::Param::thickness, thickness_rel_err2, StripCPE::Param::topology, HLT_2018_cff::track, reco::InverseCrosstalkMatrix::unfold(), LocalTrajectoryParameters::vector(), StripCPE::xtalk1, and PV3DBase< T, PVType, FrameType >::z().

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

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

private

Definition at line 36 of file StripCPEgeometric.cc.

References StripCPEgeometric::WrappedCluster::addSuppressedEdgeStrip(), ambiguousSize(), StripCPEgeometric::WrappedCluster::centroid(), StripCPEgeometric::WrappedCluster::deformed(), StripCPEgeometric::WrappedCluster::dropSmallerEdgeStrip(), stats_t< T >::error2(), dqmMemoryStats::float, geometric_position(), SiStripPI::max, StripCPEgeometric::WrappedCluster::maxProjection(), StripCPEgeometric::WrappedCluster::middle(), StripCPEgeometric::WrappedCluster::N, funct::pow(), amptDefault_cfi::proj, class-composition::Q, 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);
 }

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

inlineprivate

Definition at line 73 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;
 }

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

inlineprivate

Definition at line 69 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;
 }