trajectoryStateTransform Namespace Reference

Functions
FreeTrajectoryState	initialFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
FreeTrajectoryState	initialFreeStateL1TTrack (const TTTrack< Ref_Phase2TrackerDigi_ > &tk, const MagneticField *field, bool withErr=false)
FreeTrajectoryState	innerFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
TrajectoryStateOnSurface	innerStateOnSurface (const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
FreeTrajectoryState	outerFreeState (const reco::Track &tk, const MagneticField *field, bool withErr=true)
TrajectoryStateOnSurface	outerStateOnSurface (const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
PTrajectoryStateOnDet	persistentState (const TrajectoryStateOnSurface &ts, unsigned int detid)
TrajectoryStateOnSurface	transientState (const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)

Function Documentation

initialFreeState()

FreeTrajectoryState trajectoryStateTransform::initialFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Construct a FreeTrajectoryState from the reco::Track innermost or outermost state, does not require access to tracking geometry

Definition at line 58 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), reco::TrackBase::covariance(), submitPVResolutionJobs::err, reco::TrackBase::momentum(), and reco::TrackBase::vertex().

Referenced by TrackDetectorAssociator::associate(), reco::CandidatePtrTransientTrack::CandidatePtrTransientTrack(), TrackAssociatorByPositionImpl::getState(), pat::PATIsolatedTrackProducer::getTrackDetMatchInfo(), reco::GsfTransientTrack::GsfTransientTrack(), TSGForRoadSearch::IPfts(), TSGFromPropagation::outerTkState(), FastTSGFromPropagation::outerTkState(), TSGForOIFromL2::produce(), OutsideInMuonSeeder::produce(), TSGForOIDNN::produce(), PropagateToMuon::startingState(), MatcherUsingTracksAlgorithm::startingState(), and reco::TrackTransientTrack::TrackTransientTrack().

                                                                                                      {
    Basic3DVector<float> pos(tk.vertex());
    GlobalPoint gpos(pos);
    Basic3DVector<float> mom(tk.momentum());
    GlobalVector gmom(mom);
    GlobalTrajectoryParameters par(gpos, gmom, tk.charge(), field);
    if (!withErr)
      return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err(tk.covariance());
    return FreeTrajectoryState(par, err);
  }

initialFreeStateL1TTrack()

FreeTrajectoryState trajectoryStateTransform::initialFreeStateL1TTrack	(	const TTTrack< Ref_Phase2TrackerDigi_ > &	tk,
		const MagneticField *	field,
		bool	withErr = false
	)

Definition at line 70 of file TrajectoryStateTransform.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, MillePedeFileConverter_cfg::e, TTTrack< T >::momentum(), TTTrack< T >::POCA(), and TTTrack< T >::rInv().

Referenced by SeedGeneratorFromL1TTracksEDProducer::produce().

                                                             {
    Basic3DVector<float> pos(tk.POCA());
    GlobalPoint gpos(pos);
    GlobalVector gmom = tk.momentum();
    int charge = tk.rInv() > 0.f ? 1 : -1;

    GlobalTrajectoryParameters par(gpos, gmom, charge, field);
    if (!withErr)
      return FreeTrajectoryState(par);
    AlgebraicSymMatrix55 mat = AlgebraicMatrixID();
    mat *= 1e-8;
    return FreeTrajectoryState(par, mat);
  }

innerFreeState()

FreeTrajectoryState trajectoryStateTransform::innerFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 86 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), submitPVResolutionJobs::err, reco::Track::extra(), reco::Track::innerMomentum(), and reco::Track::innerPosition().

Referenced by TrackDetectorAssociator::associate(), innerStateOnSurface(), EgammaHLTPixelMatchElectronAlgo::process(), PropagateToMuon::startingState(), and MatcherUsingTracksAlgorithm::startingState().

                                                                                                    {
    Basic3DVector<float> pos(tk.innerPosition());
    GlobalPoint gpos(pos);
    Basic3DVector<float> mom(tk.innerMomentum());
    GlobalVector gmom(mom);
    GlobalTrajectoryParameters par(gpos, gmom, tk.charge(), field);
    if (!withErr)
      return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err(tk.extra()->innerStateCovariance());
    return FreeTrajectoryState(par, err);
  }

innerStateOnSurface()

TrajectoryStateOnSurface trajectoryStateTransform::innerStateOnSurface	(	const reco::Track &	tk,
		const TrackingGeometry &	geom,
		const MagneticField *	field,
		bool	withErr = true
	)

Construct a TrajectoryStateOnSurface from the reco::Track innermost or outermost state, requires access to tracking geometry

Definition at line 110 of file TrajectoryStateTransform.cc.

References reco::Track::extra(), relativeConstraints::geom, and innerFreeState().

Referenced by ConversionProducer::checkPhi(), GlobalMuonTrackMatcher::convertToTSOSMuHit(), GlobalMuonTrackMatcher::convertToTSOSTkHit(), OutsideInMuonSeeder::doDebug(), MuonShowerInformationFiller::getCompatibleDets(), NuclearTrackCorrector::getInitialState(), GlobalCosmicMuonTrajectoryBuilder::getTransientRecHits(), GlobalTrajectoryBuilderBase::getTransientRecHits(), ConversionSeedFilter::getTSOS(), reco::CandidatePtrTransientTrack::innermostMeasurementState(), reco::TrackTransientTrack::innermostMeasurementState(), TSGFromPropagation::innerState(), FastTSGFromPropagation::innerState(), CalibrationTrackSelectorFromDetIdList::makeCandidate(), reco::modules::CosmicTrackSplitter::makeCandidate(), reco::modules::TrackerTrackHitFilter::makeCandidate(), GlobalCosmicMuonTrajectoryBuilder::match(), TrackMerger::merge(), TSGForOIFromL2::produce(), OutsideInMuonSeeder::produce(), TSGForOIDNN::produce(), TkConvValidator::recalculateMomentumAtFittedVertex(), MatcherUsingTracksAlgorithm::targetState(), and GlobalCosmicMuonTrajectoryBuilder::trajectories().

                                                             {
    const Surface& surface = geom.idToDet(DetId(tk.extra()->innerDetId()))->surface();
    return TrajectoryStateOnSurface(innerFreeState(tk, field, withErr), surface);
  }

outerFreeState()

FreeTrajectoryState trajectoryStateTransform::outerFreeState	(	const reco::Track &	tk,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 98 of file TrajectoryStateTransform.cc.

References reco::TrackBase::charge(), submitPVResolutionJobs::err, reco::Track::extra(), reco::Track::outerMomentum(), and reco::Track::outerPosition().

Referenced by TICLDumper::analyze(), TrackDetectorAssociator::associate(), ticl::LinkingAlgoByDirectionGeometric::linkTracksters(), ticl::SeedingRegionByTracks::makeRegions(), outerStateOnSurface(), PropagateToMuon::startingState(), and MatcherUsingTracksAlgorithm::startingState().

                                                                                                    {
    Basic3DVector<float> pos(tk.outerPosition());
    GlobalPoint gpos(pos);
    Basic3DVector<float> mom(tk.outerMomentum());
    GlobalVector gmom(mom);
    GlobalTrajectoryParameters par(gpos, gmom, tk.charge(), field);
    if (!withErr)
      return FreeTrajectoryState(par);
    CurvilinearTrajectoryError err(tk.extra()->outerStateCovariance());
    return FreeTrajectoryState(par, err);
  }

outerStateOnSurface()

TrajectoryStateOnSurface trajectoryStateTransform::outerStateOnSurface	(	const reco::Track &	tk,
		const TrackingGeometry &	geom,
		const MagneticField *	field,
		bool	withErr = true
	)

Definition at line 118 of file TrajectoryStateTransform.cc.

References reco::Track::extra(), relativeConstraints::geom, and outerFreeState().

Referenced by GlobalMuonTrackMatcher::convertToTSOSMuHit(), GlobalMuonTrackMatcher::convertToTSOSTk(), GlobalMuonTrackMatcher::convertToTSOSTkHit(), MuonShowerInformationFiller::getCompatibleDets(), NuclearTrackCorrector::getInitialState(), ConversionProducer::getTrackImpactPosition(), CalibrationTrackSelectorFromDetIdList::makeCandidate(), reco::modules::CosmicTrackSplitter::makeCandidate(), reco::modules::TrackerTrackHitFilter::makeCandidate(), GlobalCosmicMuonTrajectoryBuilder::match(), TrackMerger::merge(), reco::CandidatePtrTransientTrack::outermostMeasurementState(), reco::TrackTransientTrack::outermostMeasurementState(), HGCalTrackCollectionProducer::produce(), MatcherUsingTracksAlgorithm::targetState(), and GlobalCosmicMuonTrajectoryBuilder::trajectories().

                                                             {
    const Surface& surface = geom.idToDet(DetId(tk.extra()->outerDetId()))->surface();
    return TrajectoryStateOnSurface(outerFreeState(tk, field, withErr), surface);
  }

persistentState()

PTrajectoryStateOnDet trajectoryStateTransform::persistentState	(	const TrajectoryStateOnSurface &	ts,
		unsigned int	detid
	)

should check if corresponds to m

should check if corresponds to m

Definition at line 14 of file TrajectoryStateTransform.cc.

References TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::hasError(), mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), visualization-live-secondInstance_cfg::m, LocalTrajectoryError::matrix(), PV3DBase< T, PVType, FrameType >::perp(), DiDispStaMuonMonitor_cfi::pt, and TrajectoryStateOnSurface::surfaceSide().

Referenced by SeedFromConsecutiveHitsCreator::buildSeed(), SeedFromGenericPairOrTriplet::buildSeed(), SeedForPhotonConversion1Leg::buildSeed(), SeedForPhotonConversionFromQuadruplets::buildSeed(), SeedFromNuclearInteraction::construct(), MkFitOutputConverter::convertCandidates(), RPCSeedPattern::createFakeSeed(), MuonSeedFromRecHits::createSeed(), MuonSeedCreator::createSeed(), OutInConversionSeedFinder::createSeed(), SeedingOTEDProducer::createSeed(), RPCSeedPattern::createSeed(), InOutConversionSeedFinder::createSeed(), TSGFromPropagation::createSeed(), FastTSGFromPropagation::createSeed(), OutsideInMuonSeeder::doLayer(), SeedGeneratorFromL1TTracksEDProducer::findSeedsOnLayer(), SeedFromProtoTrack::init(), CalibrationTrackSelectorFromDetIdList::makeCandidate(), reco::modules::CosmicTrackSplitter::makeCandidate(), reco::modules::TrackerTrackHitFilter::makeCandidate(), CosmicSeedCreator::makeSeed(), SETSeedFinder::makeSeed(), TSGForOIDNN::makeSeedsFromHitDoublets(), TSGForOIFromL2::makeSeedsFromHits(), TSGForOIDNN::makeSeedsFromHits(), TSGForOIFromL2::makeSeedsWithoutHits(), TSGForOIDNN::makeSeedsWithoutHits(), TrackMerger::merge(), GlobalTrackerMuonAlignment::muonFitter(), MuonReSeeder::produce(), SETMuonSeedProducer::produce(), SeedProducerFromSoAT< TrackerTraits >::produce(), TrackCandidateProducer::produce(), TrajectorySeedFromMuonProducer::produce(), L2MuonSeedGeneratorFromL1T::produce(), ConvBremSeedProducer::produce(), L2MuonSeedGeneratorFromL1TkMu::produce(), L2MuonSeedGenerator::produce(), TrackerToMuonPropagator::produce(), CSCOverlapsTrackPreparation::produce(), cms::CosmicTrackFinder::produce(), cms::CkfTrackCandidateMakerBase::produceBase(), TSGForRoadSearch::pushTrajectorySeed(), TrackInfoProducerAlgorithm::run(), SeedGeneratorForCRack::seeds(), SeedGeneratorForCosmics::seeds(), SimpleCosmicBONSeeder::seeds(), GlobalTrackerMuonAlignment::trackFitter(), InOutConversionTrackFinder::tracks(), OutInConversionTrackFinder::tracks(), StandAloneMuonTrajectoryBuilder::trajectories(), and CosmicMuonSeedGenerator::tsosToSeed().

                                                                                                {
    int surfaceSide = static_cast<int>(ts.surfaceSide());
    auto pt = ts.globalMomentum().perp();

    if (ts.hasError()) {
      AlgebraicSymMatrix55 const& m = ts.localError().matrix();

      int dim = 5;  
      float localErrors[15];

      int k = 0;
      for (int i = 0; i < dim; i++) {
        for (int j = 0; j <= i; j++) {
          localErrors[k++] = m(i, j);
        }
      }
      return PTrajectoryStateOnDet(ts.localParameters(), pt, localErrors, detid, surfaceSide);
    }
    return PTrajectoryStateOnDet(ts.localParameters(), pt, detid, surfaceSide);
  }

transientState()

TrajectoryStateOnSurface trajectoryStateTransform::transientState	(	const PTrajectoryStateOnDet &	ts,
		const Surface *	surface,
		const MagneticField *	field
	)

Definition at line 35 of file TrajectoryStateTransform.cc.

References PTrajectoryStateOnDet::error(), PTrajectoryStateOnDet::hasError(), mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, visualization-live-secondInstance_cfg::m, PTrajectoryStateOnDet::parameters(), and PTrajectoryStateOnDet::surfaceSide().

Referenced by MuonSeedsAnalyzer::analyze(), TrackBuildingAnalyzer::analyze(), TestHits::analyze(), TestSmoothHits::analyze(), SeedMultiplicityAnalyzer::analyze(), TrackerSeedCleaner::clean(), MkFitSeedConverter::convertSeeds(), TrackingNtuple::fillCandidates(), TrackingNtuple::fillSeeds(), MuonCkfTrajectoryBuilder::findCompatibleMeasurements(), BaseCkfTrajectoryBuilder::findStateAndLayers(), MuonSeedTrack::getSeedTSOS(), MuonTrackAnalyzer::getSeedTSOS(), ConversionSeedFilter::getTSOS(), SeedMatcher::matchRecHitCombinations(), CtfSpecialSeedGenerator::postCheck(), TrackProducerWithSCAssociation::produce(), ConversionSeedFilterCharge::produce(), FakeTrackProducer< T >::produce(), TrackCandidateProducer::produce(), TrackFromSeedProducer::produce(), BaseCkfTrajectoryBuilder::seedMeasurements(), MuonSeedCleaner::SeedMomentum(), MuonSeedCleaner::SeedPosition(), SeedTransformer::seedTransientState(), SeedToTrackProducer::seedTransientState(), CosmicTrajectoryBuilder::startingTSOS(), CRackTrajectoryBuilder::startingTSOS(), FastTSGFromPropagation::trackerSeeds(), and CosmicMuonTrajectoryBuilder::trajectories().

                                                                      {
    AlgebraicSymMatrix55 m;
    bool errInv = true;
    if (ts.hasError()) {
      errInv = false;
      int dim = 5;
      int k = 0;
      for (int i = 0; i < dim; i++) {
        for (int j = 0; j <= i; j++) {
          m(i, j) = ts.error(k++);  // NOTE: here we do a cast float => double.
        }
      }
    }

    return TrajectoryStateOnSurface(ts.parameters(),
                                    errInv ? LocalTrajectoryError(InvalidError()) : LocalTrajectoryError(m),
                                    *surface,
                                    field,
                                    static_cast<SurfaceSide>(ts.surfaceSide()));
  }