GsfChi2MeasurementEstimator.cc

Go to the documentation of this file.

#include "TrackingTools/GsfTracking/interface/GsfChi2MeasurementEstimator.h"
#include "TrackingTools/GsfTools/interface/GetComponents.h"

#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/GsfTracking/interface/PosteriorWeightsCalculator.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

std::pair<bool,double> 
GsfChi2MeasurementEstimator::estimate (const TrajectoryStateOnSurface& tsos,
                       const TrackingRecHit& hit) const {

  GetComponents comps(tsos);
  auto const & tsvec = comps();
  if (tsvec.empty()) {
    edm::LogError("GsfChi2MeasurementEstimator") 
      << "Trying to calculate chi2 of hit with respect to empty mixture!";
    return std::make_pair(false,0.);
  }

  auto const & weights = PosteriorWeightsCalculator(tsvec).weights(hit);
  if ( weights.empty() )  return std::make_pair(false,0);

  //   Chi2MeasurementEstimator est(chiSquaredCut());
  double chi2 = 0.;
  int i = 0;
  for (auto const & ts : tsvec)
    chi2 += weights[i++] * theEstimator.estimate(ts,hit).second;
  // Done - normalisation of weights is ensured 
  // by PosteriorWeightsCalculator
  return returnIt(chi2);
}