MuonBestMeasurementFinder Class Reference

#include <MuonBestMeasurementFinder.h>

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. More...
std::pair< double, int >	lookAtSubRecHits (TrajectoryMeasurement *measurement, const Propagator *propagator)
	MuonBestMeasurementFinder ()
	Constructor. More...
virtual	~MuonBestMeasurementFinder ()
	Destructor. More...

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.

Author

R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

S. Lacaprara - INFN Legnaro stefa.nosp@m.no.l.nosp@m.acapr.nosp@m.ara@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t

Definition at line 24 of file MuonBestMeasurementFinder.h.

Member Typedef Documentation

typedef std::vector<TrajectoryMeasurement*> MuonBestMeasurementFinder::TMContainer

private

Definition at line 25 of file MuonBestMeasurementFinder.h.

typedef TMContainer::iterator MuonBestMeasurementFinder::TMIterator

private

Definition at line 26 of file MuonBestMeasurementFinder.h.

Constructor & Destructor Documentation

MuonBestMeasurementFinder::MuonBestMeasurementFinder

(

)

Constructor.

Definition at line 29 of file MuonBestMeasurementFinder.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator.

                                                    {
  
  theEstimator = new Chi2MeasurementEstimator(100000.);
}

MuonBestMeasurementFinder::~MuonBestMeasurementFinder ( )

virtual

Destructor.

Definition at line 34 of file MuonBestMeasurementFinder.cc.

                                                     {
  delete theEstimator;
}

Member Function Documentation

const MeasurementEstimator* MuonBestMeasurementFinder::estimator ( void  ) const

inline

Definition at line 45 of file MuonBestMeasurementFinder.h.

References theEstimator.

{ 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 39 of file MuonBestMeasurementFinder.cc.

References LogTrace, metname, TrajectoryMeasurement::recHit(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by CosmicMuonTrajectoryBuilder::findBestMeasurements(), and StandAloneMuonFilter::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 93 of file MuonBestMeasurementFinder.cc.

References gsfElectronCkfTrackCandidateMaker_cff::estimator, metname, TrajectoryMeasurement::predictedState(), TrajectoryMeasurement::recHit(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

  TransientTrackingRecHit::ConstRecHitPointer muonRecHit = measurement->recHit();
  TrajectoryStateOnSurface predState = measurement->predictedState();                          // temporarily introduced by DT
  std::pair<bool, double> sing_chi2 = estimator()->estimate( predState, *(muonRecHit.get()));  // temporarily introduced by DT
  npts = 1;                                                                                    // temporarily introduced by DT
  std::pair<double, int> result = pair<double, int>(sing_chi2.second, npts);                   // temporarily introduced by DT
  
//   // ask for the 2D-segments/2D-rechit                                                      // temporarily excluded by DT
//   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 51 of file MuonBestMeasurementFinder.h.

Referenced by estimator().