track_associator Namespace Reference

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 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)
double	trackAssociationChi2 (const reco::TrackBase &track, const TrackingParticle &trackingParticle, const MagneticField &magfield, const reco::BeamSpot &bs)

Variables
double	invalidChi2 = 10000000000.

Function Documentation

const TrackingRecHit* track_associator::getHitFromIter

(

trackingRecHit_iterator

iter

)

Definition at line 14 of file trackHitsToClusterRefs.h.

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

{ return &(**iter); }

template<typename iter >

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

Definition at line 19 of file trackHitsToClusterRefs.h.

References SiStripRecHit1D::cluster(), Phase2TrackerRecHit1D::cluster(), SiStripRecHit2D::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().

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

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.

References vertices_cff::chi2, reco::deltaPhi(), f, objects.autophobj::float, invalidChi2, LogDebug, funct::sin(), and reco::trackingParametersAtClosestApproachToBeamSpot().

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

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

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.

References ALCARECOTkAlJpsiMuMu_cff::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 63 of file trackAssociationChi2.cc.

References reco::TrackBase::covariance(), reco::TrackBase::parameters(), and trackAssociationChi2().

                                                      {
    auto cov = track.covariance();
    cov.Invert();

    return trackAssociationChi2(track.parameters(), cov, trackingParticle, magfield, bs);
  }

Variable Documentation

double track_associator::invalidChi2 = 10000000000.

Definition at line 7 of file trackAssociationChi2.cc.

Referenced by trackAssociationChi2().