track_associator Namespace Reference

Functions
double	getChi2 (const reco::TrackBase::ParameterVector &rParameters, const reco::TrackBase::CovarianceMatrix &recoTrackCovMatrix, const Basic3DVector< double > &momAtVtx, const Basic3DVector< double > &vert, int charge, const MagneticField &, const reco::BeamSpot &)
	basic method where chi2 is computed More...

Function Documentation

double track_associator::getChi2	(	const reco::TrackBase::ParameterVector &	rParameters,
		const reco::TrackBase::CovarianceMatrix &	recoTrackCovMatrix,
		const Basic3DVector< double > &	momAtVtx,
		const Basic3DVector< double > &	vert,
		int	charge,
		const MagneticField &	magfield,
		const reco::BeamSpot &	bs
	)

basic method where chi2 is computed

Definition at line 7 of file getChi2.cc.

References HLT_FULL_cff::chi2, reco::deltaPhi(), f, LogDebug, funct::sin(), and reco::trackingParametersAtClosestApproachToBeamSpot().

Referenced by TrackGenAssociatorByChi2Impl::associateGenToReco(), TrackGenAssociatorByChi2Impl::associateRecoToGen(), TrackAssociatorByChi2Impl::associateRecoToSim(), TrackAssociatorByChi2Impl::associateSimToReco(), TrackGenAssociatorByChi2Impl::getChi2(), and TrackAssociatorByChi2Impl::getChi2().

                                         {
    
    double chi2 = 10000000000.;
    
    std::pair<bool,reco::TrackBase::ParameterVector> params = reco::trackingParametersAtClosestApproachToBeamSpot(vert, momAtVtx, charge, magfield, bs);
    if (params.first){
      reco::TrackBase::ParameterVector sParameters=params.second;
      
      reco::TrackBase::ParameterVector diffParameters = rParameters - sParameters;
      diffParameters[2] = reco::deltaPhi(diffParameters[2],0.f);
      chi2 = ROOT::Math::Dot(diffParameters * recoTrackCovMatrix, diffParameters);
      chi2 /= 5;
      
      LogDebug("TrackAssociator") << "====NEW RECO TRACK WITH PT=" << sin(rParameters[1])*float(charge)/rParameters[0] << "====\n" 
                                  << "qoverp sim: " << sParameters[0] << "\n" 
                                  << "lambda sim: " << sParameters[1] << "\n" 
                                  << "phi    sim: " << sParameters[2] << "\n" 
                                  << "dxy    sim: " << sParameters[3] << "\n" 
                                  << "dsz    sim: " << sParameters[4] << "\n" 
                                  << ": " /*<< */ << "\n" 
                                  << "qoverp rec: " << rParameters[0] << "\n" 
                                  << "lambda rec: " << rParameters[1] << "\n" 
                                  << "phi    rec: " << rParameters[2] << "\n" 
                                  << "dxy    rec: " << rParameters[3] << "\n" 
                                  << "dsz    rec: " << rParameters[4] << "\n" 
                                  << ": " /*<< */ << "\n" 
                                  << "chi2: " << chi2 << "\n";
      
    }
    return chi2;
  }