StripCPEgeometric Class Reference

#include <StripCPEgeometric.h>

Inheritance diagram for StripCPEgeometric:

Classes
class	WrappedCluster

Public Member Functions
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
virtual StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, AlgoParam const &ap) const
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const override
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
	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 void	localParameters (AClusters const &clusters, ALocalValues &retValues, const GeomDetUnit &gd, const LocalTrajectoryParameters &ltp) const
virtual LocalValues	localParameters (const SiStripCluster &, const GeomDetUnit &) const
virtual StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, AlgoParam const &ap) const
StripClusterParameterEstimator::LocalValues	localParameters (const SiStripCluster &cl, const GeomDetUnit &) const override
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
	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.

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

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

Referenced by offset_from_firstStrip().

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.

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

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

Referenced by offset_from_firstStrip().

localParameters() [1/7]

virtual void StripClusterParameterEstimator::localParameters

inline

Definition at line 32 of file StripClusterParameterEstimator.h.

                                                                           {}

localParameters() [2/7]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 37 of file StripClusterParameterEstimator.h.

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

localParameters() [3/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.

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

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(), digitizers_cfi::strip, tan_diffusion_angle, thickness_rel_err2, HLT_FULL_cff::track, reco::InverseCrosstalkMatrix::unfold(), LocalTrajectoryParameters::vector(), and StripCPE::xtalk1.

localParameters() [4/7]

virtual StripClusterParameterEstimator::LocalValues StripCPE::localParameters

inline

Definition at line 52 of file StripCPE.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 40 of file StripClusterParameterEstimator.h.

                                                                              {
    return localParameters(cluster, gd);
  }

localParameters() [7/7]

virtual LocalValues StripClusterParameterEstimator::localParameters

inline

Definition at line 45 of file StripClusterParameterEstimator.h.

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

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.

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

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, useNMinusOne(), and useNPlusOne().

Referenced by localParameters().

useNMinusOne()

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

inlineprivate

Definition at line 74 of file StripCPEgeometric.cc.

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

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

useNPlusOne()

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

inlineprivate

Definition at line 70 of file StripCPEgeometric.cc.

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

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

Referenced by offset_from_firstStrip().

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