MRHChi2MeasurementEstimator Class Reference

#include <MRHChi2MeasurementEstimator.h>

Inheritance diagram for MRHChi2MeasurementEstimator:

Public Member Functions
MRHChi2MeasurementEstimator *	clone () const override
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &tsos, const TrackingRecHit &aRecHit) const override
template<unsigned int N>
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &tsos, const TrackingRecHit &aRecHit) const
	MRHChi2MeasurementEstimator (double maxChi2, double nSigma=3.)
Public Member Functions inherited from Chi2MeasurementEstimatorBase
	Chi2MeasurementEstimatorBase (double maxChi2, double nSigma=3., float maxDisp=std::numeric_limits< float >::max())
template<typename... Args>
	Chi2MeasurementEstimatorBase (double maxChi2, double nSigma, float maxDisp, Args &&... args)
double	chiSquaredCut () const
bool	estimate (const TrajectoryStateOnSurface &ts, const Plane &plane) const final
Local2DVector	maximalLocalDisplacement (const TrajectoryStateOnSurface &ts, const Plane &plane) const final
double	nSigmaCut () const
Public Member Functions inherited from MeasurementEstimator
float	maxSagitta () const
	MeasurementEstimator ()
	MeasurementEstimator (float maxSag, float minToll, float mpt)
float	minPt2ForHitRecoveryInGluedDet () const
float	minTolerance2 () const
virtual bool	preFilter (const TrajectoryStateOnSurface &, OpaquePayload const &) const
virtual	~MeasurementEstimator ()

Additional Inherited Members
Public Types inherited from MeasurementEstimator
using	HitReturnType = std::pair< bool, double >
using	Local2DVector = Vector2DBase< float, LocalTag >
using	SurfaceReturnType = bool
Protected Member Functions inherited from Chi2MeasurementEstimatorBase
std::pair< bool, double >	returnIt (double est) const

Detailed Description

Definition at line 8 of file MRHChi2MeasurementEstimator.h.

Constructor & Destructor Documentation

MRHChi2MeasurementEstimator()

MRHChi2MeasurementEstimator::MRHChi2MeasurementEstimator ( double  maxChi2, double  nSigma = 3.  )

inlineexplicit

Construct with cuts on chi2 and nSigma. The cut on Chi2 is used to define the acceptance of RecHits. The errors of the trajectory state are multiplied by nSigma to define acceptance of Plane and maximalLocalDisplacement.

Definition at line 16 of file MRHChi2MeasurementEstimator.h.

Referenced by clone().

      : Chi2MeasurementEstimatorBase(maxChi2, nSigma) {}

Member Function Documentation

clone()

MRHChi2MeasurementEstimator* MRHChi2MeasurementEstimator::clone ( void  ) const

inlineoverridevirtual

Implements MeasurementEstimator.

Definition at line 23 of file MRHChi2MeasurementEstimator.h.

References MRHChi2MeasurementEstimator().

{ return new MRHChi2MeasurementEstimator(*this); }

estimate() [1/2]

std::pair< bool, double > MRHChi2MeasurementEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const TrackingRecHit &  hit  ) const

overridevirtual

Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible. The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.

Implements Chi2MeasurementEstimatorBase.

Definition at line 11 of file MRHChi2MeasurementEstimator.cc.

References TrackingRecHit::dimension(), Exception, LogDebug, and Skims_PA_cff::name.

                                                                                                   {
  switch (aRecHit.dimension()) {
    case 1:
      return estimate<1>(tsos, aRecHit);
    case 2:
      return estimate<2>(tsos, aRecHit);
    //avoid the not-(1D or 2D)  hit due to the final sum
    //supposing all the hits inside of a MRH have the same dimension
    case 3:
    case 4:
    case 5: {
      LogDebug("MRHChi2MeasurementEstimator") << "WARNING:The hit is not 1D either 2D:"
                                              << " does not count in the MRH Chi2 estimation.";
      double est = 0.0;
      return HitReturnType(false, est);
    }
  }
  throw cms::Exception("Rec hit of invalid dimension (not 1,2,3,4,5)")
      << "The value was " << aRecHit.dimension() << ", type is " << typeid(aRecHit).name() << "\n";
}

estimate() [2/2]

template<unsigned int N>

std::pair< bool, double > MRHChi2MeasurementEstimator::estimate ( const TrajectoryStateOnSurface &  ts, const TrackingRecHit &  hit  ) const

virtual

Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible. The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.

Implements Chi2MeasurementEstimatorBase.

Definition at line 35 of file MRHChi2MeasurementEstimator.cc.

References correctionTermsCaloMet_cff::C, makeMuonMisalignmentScenario::components, runTheMatrix::const, SiTrackerMultiRecHit::getAnnealingFactor(), invertPosDefMatrix(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LogDebug, LocalTrajectoryError::matrix(), packedPFCandidateRefMixer_cfi::pf, alignCSCRings::r, SiTrackerMultiRecHit::recHits(), Chi2MeasurementEstimatorBase::returnIt(), cms::cuda::V, LocalTrajectoryParameters::vector(), SiTrackerMultiRecHit::weight(), and x.

                                                                                                   {
  LogDebug("MRHChi2MeasurementEstimator") << "Calling MRHChi2MeasurementEstimator";
  SiTrackerMultiRecHit const& mHit = dynamic_cast<SiTrackerMultiRecHit const&>(aRecHit);
  double est = 0;

  double annealing = mHit.getAnnealingFactor();
  LogDebug("MRHChi2MeasurementEstimator") << "Current annealing factor is " << annealing;

  std::vector<const TrackingRecHit*> components = mHit.recHits();
  LogDebug("MRHChi2MeasurementEstimator") << "this hit has " << components.size() << " components";

  int iComp = 0;
  for (std::vector<const TrackingRecHit*>::const_iterator iter = components.begin(); iter != components.end();
       iter++, iComp++) {
    // define variables that will be used to setup the KfComponentsHolder
    ProjectMatrix<double, 5, N> pf;
    typename AlgebraicROOTObject<N>::Vector r, rMeas;
    typename AlgebraicROOTObject<N, N>::SymMatrix V, VMeas;
    AlgebraicVector5 x = tsos.localParameters().vector();
    const AlgebraicSymMatrix55& C = (tsos.localError().matrix());

    // setup the holder with the correct dimensions and get the values
    KfComponentsHolder holder;
    holder.template setup<N>(&r, &V, &pf, &rMeas, &VMeas, x, C);
    (**iter).getKfComponents(holder);

    r -= rMeas;
    V = V * annealing + VMeas;
    bool ierr = invertPosDefMatrix(V);
    if (!ierr) {
      edm::LogError("SiTrackerMultiRecHitUpdator")
          << "SiTrackerMultiRecHitUpdator::ComputeParameters2dim: W not valid!" << std::endl;
    }

    LogDebug("MRHChi2MeasurementEstimator") << "Hit with weight " << mHit.weight(iComp);
    est += ROOT::Math::Similarity(r, V) * mHit.weight(iComp);
  }

  return returnIt(est);
}