#include <MRHChi2MeasurementEstimator.h>

Inheritance diagram for MRHChi2MeasurementEstimator:

Public Member Functions
MRHChi2MeasurementEstimator *	clone () const override

template<unsigned int N>
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &tsos, const TrackingRecHit &aRecHit) const

std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &tsos, const TrackingRecHit &aRecHit) const override

	MRHChi2MeasurementEstimator (double maxChi2, double nSigma=3.)

Public Member Functions inherited from Chi2MeasurementEstimatorBase
template<typename... Args>
	Chi2MeasurementEstimatorBase (double maxChi2, double nSigma, float maxDisp, Args &&... args)

	Chi2MeasurementEstimatorBase (double maxChi2, double nSigma=3., float maxDisp=std::numeric_limits< float >::max())

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.

17 : Chi2MeasurementEstimatorBase(maxChi2, nSigma) {}

Referenced by clone().

Member Function Documentation

◆ clone()

MRHChi2MeasurementEstimator* MRHChi2MeasurementEstimator::clone ( void ) const

inlineoverridevirtual

Implements MeasurementEstimator.

Definition at line 23 of file MRHChi2MeasurementEstimator.h.

23 { return new MRHChi2MeasurementEstimator(*this); }

References MRHChi2MeasurementEstimator().

◆ estimate() [1/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 36 of file MRHChi2MeasurementEstimator.cc.

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

References gen::C, makeMuonMisalignmentScenario::components, 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.

◆ estimate() [2/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 12 of file MRHChi2MeasurementEstimator.cc.

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

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

Public Member Functions