MRHChi2MeasurementEstimator Class Reference

#include <MRHChi2MeasurementEstimator.h>

Inheritance diagram for MRHChi2MeasurementEstimator:

Public Member Functions
virtual MRHChi2MeasurementEstimator *	clone () const
std::pair< bool, double >	estimate (const TrajectoryStateOnSurface &tsos, const TrackingRecHit &aRecHit) const
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
virtual bool	estimate (const TrajectoryStateOnSurface &ts, const Plane &plane) const final
virtual 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 ( 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 17 of file MRHChi2MeasurementEstimator.h.

References estimate().

Referenced by clone().

                                                                           : 
    Chi2MeasurementEstimatorBase( maxChi2, nSigma) {}

Member Function Documentation

virtual MRHChi2MeasurementEstimator* MRHChi2MeasurementEstimator::clone ( void  ) const

inlinevirtual

Implements MeasurementEstimator.

Definition at line 27 of file MRHChi2MeasurementEstimator.h.

References MRHChi2MeasurementEstimator().

                                                     {
    return new MRHChi2MeasurementEstimator(*this);
  }

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 12 of file MRHChi2MeasurementEstimator.cc.

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

Referenced by MRHChi2MeasurementEstimator().

                                                                        {

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

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.

References patCaloMETCorrections_cff::C, makeMuonMisalignmentScenario::components, SiTrackerMultiRecHit::getAnnealingFactor(), invertPosDefMatrix(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LogDebug, LocalTrajectoryError::matrix(), packedPFCandidateRefMixer_cfi::pf, alignCSCRings::r, SiTrackerMultiRecHit::recHits(), Chi2MeasurementEstimatorBase::returnIt(), 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);

}