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 23 of file trackHitsToClusterRefs.h.

References SiStripRecHit1D::cluster(), Phase2TrackerRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), end, TrackingRecHit::geographicalId(), getHitFromIter(), edm::Ref< C, T, F >::isNonnull(), TrackingRecHit::isValid(), 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 (rhit->isValid()) {
        int subdetid = rhit->geographicalId().subdetId();
        if (subdetid==PixelSubdetector::PixelBarrel||subdetid==PixelSubdetector::PixelEndcap) {
          const SiPixelRecHit* pRHit = dynamic_cast<const SiPixelRecHit*>(rhit);
          if (!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 64 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().