---

MuonBestMeasurementFinder Class Reference

Algorithmic class to get best measurement from a list of TM the chi2 cut for the MeasurementEstimator is huge since should not be used. More...

#include <RecoMuon/TrackingTools/interface/MuonBestMeasurementFinder.h>

List of all members.

Public Member Functions
const MeasurementEstimator *	estimator () const
TrajectoryMeasurement *	findBestMeasurement (std::vector< TrajectoryMeasurement > &measC, const Propagator *propagator)
	return the Tm with the best chi2: no cut applied.
std::pair< double, int >	lookAtSubRecHits (TrajectoryMeasurement *measurement, const Propagator *propagator)
	MuonBestMeasurementFinder ()
	Constructor.
virtual	~MuonBestMeasurementFinder ()
	Destructor.
Private Types
typedef std::vector < TrajectoryMeasurement * >	TMContainer
typedef TMContainer::iterator	TMIterator
Private Attributes
const MeasurementEstimator *	theEstimator

---

Detailed Description

Algorithmic class to get best measurement from a list of TM the chi2 cut for the MeasurementEstimator is huge since should not be used.

The aim of this class is to return the "best" measurement according to the chi2, but without any cut. The decision whether to use or not the measurement is taken in the caller class. The evaluation is made (in hard-code way) with the granularity = 1. Where the granularity is the one defined in the MuonTrajectoyUpdatorClass.

Date

2007/12/19 15:46:03

Revision

1.6

Author:

R. Bellan - INFN Torino <riccardo.bellan@cern.ch>

S. Lacaprara - INFN Legnaro <stefano.lacaprara@pd.infn.it>

The aim of this class is to return the "best" measurement according to the chi2, but without any cut. The decision whether to use or not the measurement is taken in the caller class. The evaluation is made (in hard-code way) with the granularity = 1. Where the granularity is the one defined in the MuonTrajectoyUpdatorClass.

Date

2007/12/19 15:46:40

Revision

1.14

Author:

R. Bellan - INFN Torino <riccardo.bellan@cern.ch>

S. Lacaprara - INFN Legnaro <stefano.lacaprara@pd.infn.it>

Definition at line 26 of file MuonBestMeasurementFinder.h.

---

Member Typedef Documentation

typedef std::vector<TrajectoryMeasurement*> MuonBestMeasurementFinder::TMContainer [private]

Definition at line 27 of file MuonBestMeasurementFinder.h.

typedef TMContainer::iterator MuonBestMeasurementFinder::TMIterator [private]

Definition at line 28 of file MuonBestMeasurementFinder.h.

---

Constructor & Destructor Documentation

MuonBestMeasurementFinder::MuonBestMeasurementFinder

(

)

Constructor.

Definition at line 31 of file MuonBestMeasurementFinder.cc.

References Chi2MeasurementEstimatorESProducer_cfi::Chi2MeasurementEstimator, and theEstimator.

                                                    {
  
  theEstimator = new Chi2MeasurementEstimator(100000.);
}

MuonBestMeasurementFinder::~MuonBestMeasurementFinder

(

)

[virtual]

Destructor.

Definition at line 36 of file MuonBestMeasurementFinder.cc.

References theEstimator.

                                                     {
  delete theEstimator;
}

---

Member Function Documentation

const MeasurementEstimator* MuonBestMeasurementFinder::estimator

(

void

)

const [inline]

Definition at line 47 of file MuonBestMeasurementFinder.h.

References theEstimator.

Referenced by lookAtSubRecHits().

{ return theEstimator;}

TrajectoryMeasurement * MuonBestMeasurementFinder::findBestMeasurement	(	std::vector< TrajectoryMeasurement > &	measC,
		const Propagator *	propagator
	)

return the Tm with the best chi2: no cut applied.

Definition at line 41 of file MuonBestMeasurementFinder.cc.

References lat::endl(), LogTrace, and lookAtSubRecHits().

Referenced by StandAloneMuonFilter::findBestMeasurements(), and CosmicMuonTrajectoryBuilder::findBestMeasurements().

                                                                            {
  
  const std::string metname = "Muon|RecoMuon|MuonBestMeasurementFinder";

  TMContainer validMeasurements;

  TrajectoryMeasurement* bestMeasurement=0;

  // consider only valid TM
  int NumValidMeas=0;
  for ( vector<TrajectoryMeasurement>::iterator measurement = measC.begin(); 
        measurement!= measC.end(); ++measurement ) {
    if ((*measurement).recHit()->isValid()) {
      ++NumValidMeas;
      bestMeasurement = &(*measurement);
      validMeasurements.push_back( &(*measurement) );
    }
  }

  // If we have just one (or zero) valid meas, return it at once 
  // (or return null measurement)
  if(NumValidMeas<=1) {
    LogTrace(metname) << "MuonBestMeasurement: just " << NumValidMeas
                      << " valid measurement ";
    return bestMeasurement;
  }

  TMIterator measurement;
  double minChi2PerNDoF=1.E6;

  // if there are more than one valid measurement, then sort them.
  for ( measurement = validMeasurements.begin(); measurement!= validMeasurements.end(); measurement++ ) {

    TransientTrackingRecHit::ConstRecHitPointer muonRecHit = (*measurement)->recHit();
    
    // FIXME!! FIXME !! FIXME !!
    pair<double,int> chi2Info = lookAtSubRecHits(*measurement, propagator);

    double chi2PerNDoF = chi2Info.first/chi2Info.second;
    LogTrace(metname) << " The measurement has a chi2/npts " << chi2PerNDoF << " with dof = " << chi2Info.second
                      << " \n Till now the best chi2 is " << minChi2PerNDoF;
    
    if ( chi2PerNDoF && chi2PerNDoF<minChi2PerNDoF ) {
      minChi2PerNDoF = chi2PerNDoF;     
      bestMeasurement = *measurement;
    }
    
  }
  LogTrace(metname)<<"The final best chi2 is "<<minChi2PerNDoF<<endl;
  return bestMeasurement;
}

pair< double, int > MuonBestMeasurementFinder::lookAtSubRecHits	(	TrajectoryMeasurement *	measurement,
		const Propagator *	propagator
	)

Definition at line 95 of file MuonBestMeasurementFinder.cc.

References lat::endl(), MeasurementEstimator::estimate(), estimator(), TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::isValid(), LogTrace, TrajectoryMeasurement::predictedState(), Propagator::propagate(), TrajectoryMeasurement::recHit(), and HLT_VtxMuL3::result.

Referenced by findBestMeasurement().

                                                                                          {
  
  const std::string metname = "Muon|RecoMuon|MuonBestMeasurementFinder";

  unsigned int npts=0;
  double thisChi2 = 0.;

  TransientTrackingRecHit::ConstRecHitPointer muonRecHit = measurement->recHit();

  
  // ask for the 2D-segments/2D-rechit
  TransientTrackingRecHit::ConstRecHitContainer rhits_list = muonRecHit->transientHits();
  
  LogTrace(metname)<<"Number of rechits in the measurement rechit: "<<rhits_list.size()<<endl;
  
  // loop over them
  for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator rhit = rhits_list.begin(); 
       rhit!= rhits_list.end(); ++rhit)
    if ((*rhit)->isValid() ) {
      LogTrace(metname)<<"Rechit dimension: "<<(*rhit)->dimension()<<endl;
      npts+=(*rhit)->dimension();
      
      TrajectoryStateOnSurface predState;
      
      if (!( (*rhit)->geographicalId() == muonRecHit->geographicalId() ) ){
        predState = propagator->propagate(*measurement->predictedState().freeState(),
                                          (*rhit)->det()->surface()); 
      }
      else predState = measurement->predictedState();  
      
      if ( predState.isValid() ) { 
        std::pair<bool,double> sing_chi2 = estimator()->estimate( predState, *((*rhit).get()));
        thisChi2 += sing_chi2.second ;
        LogTrace(metname) << " single incremental chi2: " << sing_chi2.second;
      }
    }
  
  pair<double,int> result = pair<double,int>(thisChi2,npts);
  return result;
}

---

Member Data Documentation

const MeasurementEstimator* MuonBestMeasurementFinder::theEstimator [private]

Definition at line 53 of file MuonBestMeasurementFinder.h.

Referenced by estimator(), MuonBestMeasurementFinder(), and ~MuonBestMeasurementFinder().

---

The documentation for this class was generated from the following files:

• RecoMuon/TrackingTools/interface/MuonBestMeasurementFinder.h
• RecoMuon/TrackingTools/src/MuonBestMeasurementFinder.cc

---