Functions
const TrackingRecHit *	getHitFromIter (trackingRecHit_iterator iter)

template<typename iter >
std::vector< OmniClusterRef >	hitsToClusterRefs (iter begin, iter end)

double	trackAssociationChi2 (const reco::TrackBase &track, const TrackingParticle &trackingParticle, const MagneticField &magfield, const reco::BeamSpot &bs)

double	trackAssociationChi2 (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 More...

double	trackAssociationChi2 (const reco::TrackBase::ParameterVector &rParameters, const reco::TrackBase::CovarianceMatrix &recoTrackCovMatrix, const TrackingParticle &trackingParticle, const MagneticField &magfield, const reco::BeamSpot &bs)

Variables
constexpr double	invalidChi2 = 10000000000.

Function Documentation

◆ getHitFromIter()

const TrackingRecHit * track_associator::getHitFromIter ( trackingRecHit_iterator iter )

Definition at line 14 of file trackHitsToClusterRefs.h.

14 { return &(**iter); }

Referenced by QuickTrackAssociatorByHitsImpl::getAllSimTrackIdentifiers(), QuickTrackAssociatorByHitsImpl::getDoubleCount(), hitsToClusterRefs(), and QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters().

◆ hitsToClusterRefs()

template<typename iter >

std::vector<OmniClusterRef> track_associator::hitsToClusterRefs	(	iter	begin,
		iter	end
	)

Definition at line 19 of file trackHitsToClusterRefs.h.

                                                                     {
     std::vector<OmniClusterRef> returnValue;
     for (iter iRecHit = begin; iRecHit != end; ++iRecHit) {
       const TrackingRecHit *rhit = getHitFromIter(iRecHit);
       if (trackerHitRTTI::isFromDet(*rhit)) {
         int subdetid = rhit->geographicalId().subdetId();
         if (subdetid == PixelSubdetector::PixelBarrel || subdetid == PixelSubdetector::PixelEndcap) {
           const SiPixelRecHit *pRHit = dynamic_cast<const SiPixelRecHit *>(rhit);
           if (pRHit && !pRHit->cluster().isNonnull())
             edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                              << " file: " << __FILE__ << " line: " << __LINE__;
           returnValue.push_back(pRHit->omniClusterRef());
         } else if (subdetid == SiStripDetId::TIB || subdetid == SiStripDetId::TOB || subdetid == SiStripDetId::TID ||
                    subdetid == SiStripDetId::TEC) {
           const std::type_info &tid = typeid(*rhit);
           if (tid == typeid(SiStripMatchedRecHit2D)) {
             const SiStripMatchedRecHit2D *sMatchedRHit = dynamic_cast<const SiStripMatchedRecHit2D *>(rhit);
             if (!sMatchedRHit->monoHit().cluster().isNonnull() || !sMatchedRHit->stereoHit().cluster().isNonnull())
               edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                                << " file: " << __FILE__ << " line: " << __LINE__;
             returnValue.push_back(sMatchedRHit->monoClusterRef());
             returnValue.push_back(sMatchedRHit->stereoClusterRef());
           } else if (tid == typeid(SiStripRecHit2D)) {
             const SiStripRecHit2D *sRHit = dynamic_cast<const SiStripRecHit2D *>(rhit);
             if (!sRHit->cluster().isNonnull())
               edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                                << " file: " << __FILE__ << " line: " << __LINE__;
             returnValue.push_back(sRHit->omniClusterRef());
           } else if (tid == typeid(SiStripRecHit1D)) {
             const SiStripRecHit1D *sRHit = dynamic_cast<const SiStripRecHit1D *>(rhit);
             if (!sRHit->cluster().isNonnull())
               edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                                << " file: " << __FILE__ << " line: " << __LINE__;
             returnValue.push_back(sRHit->omniClusterRef());
           } else if (tid == typeid(Phase2TrackerRecHit1D)) {
             const Phase2TrackerRecHit1D *ph2Hit = dynamic_cast<const Phase2TrackerRecHit1D *>(rhit);
             if (!ph2Hit->cluster().isNonnull())
               edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                                << " file: " << __FILE__ << " line: " << __LINE__;
             returnValue.push_back(ph2Hit->omniClusterRef());
           } else {
             auto const &thit = static_cast<BaseTrackerRecHit const &>(*rhit);
             if (thit.isProjected()) {
             } else {
               edm::LogError("TrackAssociator") << ">>> getMatchedClusters: TrackingRecHit not associated to "
                                                   "any SiStripCluster! subdetid = "
                                                << subdetid;
             }
           }
         } else {
           edm::LogError("TrackAssociator") << ">>> getMatchedClusters: TrackingRecHit not associated to any "
                                               "cluster! subdetid = "
                                            << subdetid;
         }
       }
     }
     return returnValue;
   }

References begin, SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), Phase2TrackerRecHit1D::cluster(), SiPixelRecHit::cluster(), end, TrackingRecHit::geographicalId(), getHitFromIter(), trackerHitRTTI::isFromDet(), edm::Ref< C, T, F >::isNonnull(), SiStripMatchedRecHit2D::monoClusterRef(), SiStripMatchedRecHit2D::monoHit(), TrackerSingleRecHit::omniClusterRef(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, SiStripMatchedRecHit2D::stereoClusterRef(), SiStripMatchedRecHit2D::stereoHit(), DetId::subdetId(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, and SiStripDetId::TOB.

Referenced by QuickTrackAssociatorByHitsImpl::associateTrack(), TrackingNtuple::fillSeeds(), TrackingNtuple::fillTracks(), QuickTrackAssociatorByHitsImpl::getDoubleCount(), and QuickTrackAssociatorByHitsImpl::weightedNumberOfTrackClusters().

◆ trackAssociationChi2() [1/3]

double track_associator::trackAssociationChi2	(	const reco::TrackBase &	track,
		const TrackingParticle &	trackingParticle,
		const MagneticField &	magfield,
		const reco::BeamSpot &	bs
	)

Definition at line 63 of file trackAssociationChi2.cc.

                                                       {
     auto cov = track.covariance();
     cov.Invert();
  
     return trackAssociationChi2(track.parameters(), cov, trackingParticle, magfield, bs);
   }

References cms::cuda::bs, volumeBasedMagneticField_160812_cfi::magfield, HLT_2018_cff::track, and trackAssociationChi2().

◆ trackAssociationChi2() [2/3]

double track_associator::trackAssociationChi2	(	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 9 of file trackAssociationChi2.cc.

                                                       {
     double chi2 = invalidChi2;
  
     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;
   }

References cms::cuda::bs, ALCARECOTkAlJpsiMuMu_cff::charge, hltPixelTracks_cff::chi2, reco::deltaPhi(), f, dqmMemoryStats::float, invalidChi2, LogDebug, volumeBasedMagneticField_160812_cfi::magfield, CalibrationSummaryClient_cfi::params, funct::sin(), and reco::trackingParametersAtClosestApproachToBeamSpot().

Referenced by TrackingNtuple::fillTracks(), TrackAssociatorByChi2Impl::getChi2(), TrackGenAssociatorByChi2Impl::getChi2(), and trackAssociationChi2().

◆ trackAssociationChi2() [3/3]

double track_associator::trackAssociationChi2	(	const reco::TrackBase::ParameterVector &	rParameters,
		const reco::TrackBase::CovarianceMatrix &	recoTrackCovMatrix,
		const TrackingParticle &	trackingParticle,
		const MagneticField &	magfield,
		const reco::BeamSpot &	bs
	)

Definition at line 47 of file trackAssociationChi2.cc.

                                                       {
     const int charge = trackingParticle.charge();
     if (charge == 0)
       return invalidChi2;
  
     const auto tpMom = trackingParticle.momentum();
     Basic3DVector<double> momAtVtx(tpMom.x(), tpMom.y(), tpMom.z());
     Basic3DVector<double> vert(trackingParticle.vertex());
  
     return trackAssociationChi2(rParameters, recoTrackCovMatrix, momAtVtx, vert, charge, magfield, bs);
   }

References cms::cuda::bs, ALCARECOTkAlJpsiMuMu_cff::charge, TrackingParticle::charge(), invalidChi2, volumeBasedMagneticField_160812_cfi::magfield, TrackingParticle::momentum(), trackAssociationChi2(), and TrackingParticle::vertex().

Variable Documentation

◆ invalidChi2

constexpr double track_associator::invalidChi2 = 10000000000.

constexpr

Definition at line 7 of file trackAssociationChi2.cc.

Referenced by trackAssociationChi2().

Functions

Variables