Functions
const TrackingRecHit *	getHitFromIter (trackingRecHit_iterator iter)

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

double	trackAssociationChi2 (const reco::TrackBase::ParameterVector &rParameters, const reco::TrackBase::CovarianceMatrix &recoTrackCovMatrix, const reco::TrackBase::ParameterVector &sParameters)
	basic method where chi2 is computed More...

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)

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

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

Variables
constexpr double	invalidChi2 = 10000000000.

Function Documentation

const TrackingRecHit * track_associator::getHitFromIter ( trackingRecHit_iterator iter )

inline

Definition at line 16 of file trackHitsToClusterRefs.h.

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

16 { return &(**iter); }

template<typename iter >

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

Definition at line 21 of file trackHitsToClusterRefs.h.

References SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), Phase2TrackerRecHit1D::cluster(), SiPixelRecHit::cluster(), VectorHit::cluster(), dataset::end, VectorHit::firstClusterRef(), 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().

                                                                     {
     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 if (tid == typeid(VectorHit)) {
             const VectorHit *vectorHit = dynamic_cast<const VectorHit *>(rhit);
             if (!vectorHit->cluster().isNonnull())
               edm::LogError("TrackAssociator") << ">>> RecHit does not have an associated cluster!"
                                                << " file: " << __FILE__ << " line: " << __LINE__;
             returnValue.push_back(vectorHit->firstClusterRef());
 
           } 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;
   }

double track_associator::trackAssociationChi2	(	const reco::TrackBase::ParameterVector &	rParameters,
		const reco::TrackBase::CovarianceMatrix &	recoTrackCovMatrix,
		const reco::TrackBase::ParameterVector &	sParameters
	)

basic method where chi2 is computed

Definition at line 9 of file trackAssociationChi2.cc.

References HLT_FULL_cff::chi2, reco::deltaPhi(), validate-o2o-wbm::f, invalidChi2, LogDebug, and funct::sin().

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

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

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
	)

Definition at line 36 of file trackAssociationChi2.cc.

References HLT_FULL_cff::chi2, invalidChi2, submitPVValidationJobs::params, trackAssociationChi2(), and reco::trackingParametersAtClosestApproachToBeamSpot().

                                                       {
     double chi2 = invalidChi2;
 
     std::pair<bool, reco::TrackBase::ParameterVector> params =
         reco::trackingParametersAtClosestApproachToBeamSpot(vert, momAtVtx, charge, magfield, bs);
     if (params.first) {
       return trackAssociationChi2(rParameters, recoTrackCovMatrix, params.second);
     }
     return chi2;
   }

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 53 of file trackAssociationChi2.cc.

References RecoTauCleanerPlugins::charge, TrackingParticle::charge(), invalidChi2, TrackingParticle::momentum(), trackAssociationChi2(), and TrackingParticle::vertex().

                                                       {
     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);
   }

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

Definition at line 69 of file trackAssociationChi2.cc.

References cms::cuda::bs, reco::TrackBase::covariance(), volumeBasedMagneticField_160812_cfi::magfield, reco::TrackBase::parameters(), and trackAssociationChi2().

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

Variable Documentation

constexpr double track_associator::invalidChi2 = 10000000000.

Definition at line 7 of file trackAssociationChi2.cc.

Referenced by trackAssociationChi2().

Functions

Variables

Function Documentation

Variable Documentation