MuonChi2MeasurementEstimator.cc

Go to the documentation of this file.

#include "RecoMuon/TrackingTools/interface/MuonChi2MeasurementEstimator.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "RecoMuon/TransientTrackingRecHit/interface/MuonTransientTrackingRecHit.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Common/interface/getRef.h"
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
 

MuonChi2MeasurementEstimator::MuonChi2MeasurementEstimator(double maxChi2, double nSigma)
  :Chi2MeasurementEstimatorBase(maxChi2,nSigma), 
   theDTChi2Estimator(maxChi2, nSigma), 
   theCSCChi2Estimator(maxChi2, nSigma),
   theRPCChi2Estimator(maxChi2, nSigma){}


MuonChi2MeasurementEstimator::MuonChi2MeasurementEstimator(double dtMaxChi2, double cscMaxChi2, double rpcMaxChi2, double nSigma = 3.)
  :Chi2MeasurementEstimatorBase(dtMaxChi2,nSigma),
   theDTChi2Estimator(dtMaxChi2, nSigma), 
   theCSCChi2Estimator(cscMaxChi2, nSigma),
   theRPCChi2Estimator(rpcMaxChi2, nSigma){}


std::pair<bool,double> 
MuonChi2MeasurementEstimator::estimate(const TrajectoryStateOnSurface& tsos,
                       const TrackingRecHit& recHit) const {
  
  DetId id = recHit.geographicalId();
  
  // chi2 choise based on recHit provenance
  if(id.det() == DetId::Muon){
    if(id.subdetId() == MuonSubdetId::DT)
      return theDTChi2Estimator.estimate(tsos,recHit);
    else if(id.subdetId() == MuonSubdetId::CSC)
      return theCSCChi2Estimator.estimate(tsos,recHit);
    else if(id.subdetId() == MuonSubdetId::RPC) 
      return theRPCChi2Estimator.estimate(tsos,recHit);
    else{
      edm::LogWarning("Muon|RecoMuon|MuonChi2MeasurementEstimator")
    <<"RecHit with MuonId but not with a SubDetId neither from DT, CSC or rpc. [Use the parameters used for DTs]";
      return theDTChi2Estimator.estimate(tsos,recHit);
    }
  }
  else{
    edm::LogWarning("Muon|RecoMuon|MuonChi2MeasurementEstimator")
      <<"Rechit with a non-muon det id. [Use the parameters used for DTs]";
    return theDTChi2Estimator.estimate(tsos,recHit);
  }
}