StripCPEfromTrackAngle Class Reference

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

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

AClusters

using StripCPEfromTrackAngle::AClusters = StripClusterParameterEstimator::AClusters

Definition at line 31 of file StripCPEfromTrackAngle.h.

AlgoParam

using StripCPEfromTrackAngle::AlgoParam = StripCPE::AlgoParam

Definition at line 30 of file StripCPEfromTrackAngle.h.

ALocalValues

using StripCPEfromTrackAngle::ALocalValues = StripClusterParameterEstimator::ALocalValues

Definition at line 32 of file StripCPEfromTrackAngle.h.

Member Enumeration Documentation

Algo

enum StripCPEfromTrackAngle::Algo

strongprivate

Enumerator
legacy
mergeCK
chargeCK

Definition at line 25 of file StripCPEfromTrackAngle.h.

{ legacy, mergeCK, chargeCK };

Constructor & Destructor Documentation

StripCPEfromTrackAngle()

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

legacyStripErrorSquared()

float StripCPEfromTrackAngle::legacyStripErrorSquared	(	const unsigned	N,
		const float	uProj
	)		const

Definition at line 40 of file StripCPEfromTrackAngle.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), f, myMath::fast_expf(), ALCARECOEcalPhiSym_cff::float, N, 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;
  }
}

localParameters() [1/9]

virtual LocalValues StripClusterParameterEstimator::localParameters

inlineprivate

Definition at line 45 of file StripClusterParameterEstimator.h.

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

localParameters() [2/9]

virtual void StripClusterParameterEstimator::localParameters

inlineprivate

Definition at line 32 of file StripClusterParameterEstimator.h.

                                                                           {}

localParameters() [3/9]

virtual LocalValues StripClusterParameterEstimator::localParameters

inlineprivate

Definition at line 37 of file StripClusterParameterEstimator.h.

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

localParameters() [4/9]

virtual LocalValues StripClusterParameterEstimator::localParameters

inlineprivate

Definition at line 40 of file StripClusterParameterEstimator.h.

                                                                              {
    return localParameters(cluster, gd);
  }

localParameters() [5/9]

virtual StripClusterParameterEstimator::LocalValues StripCPE::localParameters

inlineprivate

Definition at line 52 of file StripCPE.h.

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

localParameters() [6/9]

StripClusterParameterEstimator::LocalValues StripCPE::localParameters

overrideprivate

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() [7/9]

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(), bsc_activity_cfg::clusters, alignCSCRings::corr, EcalPhiSymFlatTableProducers_cfi::fill, StripCPE::getAlgoParam(), mps_fire::i, legacy, legacyStripErrorSquared(), isotrackTrainRegressor::loc, m_algo, maxChgOneMIP, mergeCK, N, AlCaHLTBitMon_ParallelJobs::p, nano_mu_digi_cff::strip, stripErrorSquared(), mitigatedMETSequence_cff::U, 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;
  }
}

localParameters() [8/9]

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, isotrackTrainRegressor::loc, StripCPE::AlgoParam::ltp, m_algo, maxChgOneMIP, mergeCK, N, StripCPE::AlgoParam::p, AlCaHLTBitMon_ParallelJobs::p, SiStripCluster::size(), nano_mu_digi_cff::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()));
}

localParameters() [9/9]

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

stripErrorSquared()

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(), ALCARECOEcalPhiSym_cff::float, isotrackTrainRegressor::loc, mHC_P, mLC_P, N, 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

m_algo

Algo StripCPEfromTrackAngle::m_algo

private

Definition at line 27 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

maxChgOneMIP

const float StripCPEfromTrackAngle::maxChgOneMIP

private

Definition at line 23 of file StripCPEfromTrackAngle.h.

Referenced by localParameters().

mHC_P

float StripCPEfromTrackAngle::mHC_P[4][2]

private

Definition at line 12 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

mHC_P0

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

private

Definition at line 15 of file StripCPEfromTrackAngle.h.

mHC_P1

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

private

Definition at line 16 of file StripCPEfromTrackAngle.h.

mLC_P

float StripCPEfromTrackAngle::mLC_P[3]

private

Definition at line 11 of file StripCPEfromTrackAngle.h.

Referenced by StripCPEfromTrackAngle(), and stripErrorSquared().

useLegacyError

const bool StripCPEfromTrackAngle::useLegacyError

private

Definition at line 19 of file StripCPEfromTrackAngle.h.