CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
TrackTransformer Class Referencefinal

#include <TrackTransformer.h>

Inheritance diagram for TrackTransformer:
TrackTransformerBase

Public Member Functions

TransientTrackingRecHit::ConstRecHitContainer getTransientRecHits (const reco::TransientTrack &track) const
 
const MagneticFieldmagneticField () const
 the magnetic field More...
 
std::unique_ptr< TrajectoryFitter > const & refitter () const
 the refitter used to refit the reco::Track More...
 
void setServices (const edm::EventSetup &) override
 set the services needed by the TrackTransformer More...
 
std::unique_ptr< TrajectorySmoother > const & smoother () const
 the smoother used to smooth the trajectory which came from the refitting step More...
 
edm::ESHandle< GlobalTrackingGeometrytrackingGeometry () const
 the tracking geometry More...
 
 TrackTransformer (const edm::ParameterSet &)
 Constructor. More...
 
std::vector< Trajectorytransform (const reco::Track &) const override
 Convert a reco::Track into Trajectory. More...
 
std::vector< Trajectorytransform (const reco::TrackRef &) const
 Convert a reco::TrackRef into Trajectory. More...
 
std::vector< Trajectorytransform (const reco::TransientTrack &, TransientTrackingRecHit::ConstRecHitContainer &) const
 Convert a reco::TrackRef into Trajectory, refit with a new set of hits. More...
 
 ~TrackTransformer () override
 Destructor. More...
 
- Public Member Functions inherited from TrackTransformerBase
 TrackTransformerBase ()
 Constructor. More...
 
virtual ~TrackTransformerBase ()
 Destructor. More...
 

Static Public Member Functions

static void fillPSetDescription (edm::ParameterSetDescription &descriptions, bool DoPredictionsOnly=false, const std::string &Fitter="KFFitterForRefitInsideOut", const std::string &Smoother="KFSmootherForRefitInsideOut", const std::string &Propagator="SmartPropagatorAnyRK", const std::string &RefitDirection="alongMomentum", bool RefitRPCHits=true, const std::string &TrackerRecHitBuilder="WithTrackAngle", const std::string &MuonRecHitBuilder="MuonRecHitBuilder", const std::string &MTDRecHitBuilder="MTDRecHitBuilder")
 fillDescriptions More...
 

Private Member Functions

RefitDirection::GeometricalDirection checkRecHitsOrdering (TransientTrackingRecHit::ConstRecHitContainer const &) const
 
edm::ESHandle< Propagator > const & propagator () const
 

Private Attributes

TkClonerImpl hitCloner
 
unsigned long long theCacheId_GTG = 0
 
unsigned long long theCacheId_MG = 0
 
unsigned long long theCacheId_TC = 0
 
unsigned long long theCacheId_TRH = 0
 
const bool theDoPredictionsOnly
 
std::unique_ptr< TrajectoryFittertheFitter
 
const std::string theFitterName
 
edm::ESHandle< MagneticFieldtheMGField
 
bool theMtdAvailable
 
edm::ESHandle< TransientTrackingRecHitBuildertheMTDRecHitBuilder
 
const std::string theMTDRecHitBuilderName
 
edm::ESHandle< TransientTrackingRecHitBuildertheMuonRecHitBuilder
 
const std::string theMuonRecHitBuilderName
 
edm::ESHandle< PropagatorthePropagator
 
const std::string thePropagatorName
 
const RefitDirection theRefitDirection
 
const bool theRPCInTheFit
 
std::unique_ptr< TrajectorySmoothertheSmoother
 
const std::string theSmootherName
 
edm::ESHandle< TransientTrackingRecHitBuildertheTrackerRecHitBuilder
 
const std::string theTrackerRecHitBuilderName
 
edm::ESHandle< GlobalTrackingGeometrytheTrackingGeometry
 

Detailed Description

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Author
R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Rec hits are ordered in rank of increasing |z| assuing muons are coming from the top of the detector.

$Date: 2009/01/15 author R. Bellan - CERN ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch modified for zed ordering by N. Kypreos - UF nicho.nosp@m.las..nosp@m.theod.nosp@m.ore..nosp@m.kypre.nosp@m.os@c.nosp@m.ern.c.nosp@m.h

This class takes a reco::Track and refits the rechits inside it. The final result is a Trajectory refitted and smoothed. To make the refitting (and the smoothing) the usual KF tools are used.

CAVEAT: till now (it will be changed in the near future) the class stores the pointers to the services, therefore EACH event the setServices(const edm::EventSetup&) method MUST be called in the code in which the TrackTransformer is used.

Author
R. Bellan - CERN ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 44 of file TrackTransformer.h.

Constructor & Destructor Documentation

TrackTransformer::TrackTransformer ( const edm::ParameterSet parameterSet)
explicit

Constructor.

Definition at line 31 of file TrackTransformer.cc.

32  : theRPCInTheFit(parameterSet.getParameter<bool>("RefitRPCHits")),
33  theDoPredictionsOnly(parameterSet.getParameter<bool>("DoPredictionsOnly")),
34  theRefitDirection(parameterSet.getParameter<string>("RefitDirection")),
35  theFitterName(parameterSet.getParameter<string>("Fitter")),
36  theSmootherName(parameterSet.getParameter<string>("Smoother")),
37  thePropagatorName(parameterSet.getParameter<string>("Propagator")),
38  theTrackerRecHitBuilderName(parameterSet.getParameter<string>("TrackerRecHitBuilder")),
39  theMuonRecHitBuilderName(parameterSet.getParameter<string>("MuonRecHitBuilder")),
40  theMTDRecHitBuilderName(parameterSet.getParameter<string>("MTDRecHitBuilder")) {}
T getParameter(std::string const &) const
const std::string theTrackerRecHitBuilderName
const RefitDirection theRefitDirection
const bool theRPCInTheFit
const std::string thePropagatorName
const std::string theFitterName
const std::string theMuonRecHitBuilderName
const std::string theSmootherName
const std::string theMTDRecHitBuilderName
const bool theDoPredictionsOnly
TrackTransformer::~TrackTransformer ( )
override

Destructor.

Definition at line 43 of file TrackTransformer.cc.

43 {}

Member Function Documentation

RefitDirection::GeometricalDirection TrackTransformer::checkRecHitsOrdering ( TransientTrackingRecHit::ConstRecHitContainer const &  recHits) const
private

Definition at line 145 of file TrackTransformer.cc.

References dqmdumpme::first, GlobalTrackingGeometry::idToDet(), RefitDirection::insideOut, dqmdumpme::last, LogDebug, PV3DBase< T, PVType, FrameType >::mag2(), RefitDirection::outsideIn, position, trackingGeometry(), and RefitDirection::undetermined.

Referenced by transform().

146  {
147  if (!recHits.empty()) {
148  GlobalPoint first = trackingGeometry()->idToDet(recHits.front()->geographicalId())->position();
149  GlobalPoint last = trackingGeometry()->idToDet(recHits.back()->geographicalId())->position();
150 
151  // maybe perp2?
152  auto rFirst = first.mag2();
153  auto rLast = last.mag2();
154  if (rFirst < rLast)
156  if (rFirst > rLast)
158  }
159  LogDebug("Reco|TrackingTools|TrackTransformer") << "Impossible to determine the rechits order" << endl;
161 }
#define LogDebug(id)
T mag2() const
Definition: PV3DBase.h:63
static int position[264][3]
Definition: ReadPGInfo.cc:289
const GeomDet * idToDet(DetId) const override
edm::ESHandle< GlobalTrackingGeometry > trackingGeometry() const
the tracking geometry
void TrackTransformer::fillPSetDescription ( edm::ParameterSetDescription descriptions,
bool  DoPredictionsOnly = false,
const std::string &  Fitter = "KFFitterForRefitInsideOut",
const std::string &  Smoother = "KFSmootherForRefitInsideOut",
const std::string &  Propagator = "SmartPropagatorAnyRK",
const std::string &  RefitDirection = "alongMomentum",
bool  RefitRPCHits = true,
const std::string &  TrackerRecHitBuilder = "WithTrackAngle",
const std::string &  MuonRecHitBuilder = "MuonRecHitBuilder",
const std::string &  MTDRecHitBuilder = "MTDRecHitBuilder" 
)
static

fillDescriptions

Definition at line 45 of file TrackTransformer.cc.

References edm::ParameterSetDescription::add(), trackingRecoMaterialAnalyzer_cfi::DoPredictionsOnly, HLTSiStripMonitoring_cff::Fitter, trackingRecoMaterialAnalyzer_cfi::MTDRecHitBuilder, trackingRecoMaterialAnalyzer_cfi::MuonRecHitBuilder, LaserSeedGenerator_cfi::Propagator, trackingRecoMaterialAnalyzer_cfi::RefitDirection, MuonAlignmentFromReference_cff::RefitRPCHits, HLTSiStripMonitoring_cff::Smoother, AlCaHLTBitMon_QueryRunRegistry::string, and trackingRecoMaterialAnalyzer_cfi::TrackerRecHitBuilder.

Referenced by L3MuonProducer::fillDescriptions(), and TrackExtenderWithMTDT< TrackCollection >::fillDescriptions().

54  {
55  desc.add<bool>("DoPredictionsOnly", DoPredictionsOnly);
56  desc.add<std::string>("Fitter", Fitter);
57  desc.add<std::string>("Smoother", Smoother);
58  desc.add<std::string>("Propagator", Propagator);
59  desc.add<std::string>("RefitDirection", RefitDirection);
60  desc.add<bool>("RefitRPCHits", RefitRPCHits);
61  desc.add<std::string>("TrackerRecHitBuilder", TrackerRecHitBuilder);
62  desc.add<std::string>("MuonRecHitBuilder", MuonRecHitBuilder);
63  desc.add<std::string>("MTDRecHitBuilder", MTDRecHitBuilder);
64 }
TransientTrackingRecHit::ConstRecHitContainer TrackTransformer::getTransientRecHits ( const reco::TransientTrack track) const

Definition at line 117 of file TrackTransformer.cc.

References TransientTrackingRecHitBuilder::build(), watchdog::const, Exception, FastTime, DetId::Forward, LogTrace, DetId::Muon, edm::ESHandle< T >::product(), reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), mps_fire::result, theMtdAvailable, theMTDRecHitBuilder, theMuonRecHitBuilder, theRPCInTheFit, theTrackerRecHitBuilder, and DetId::Tracker.

Referenced by transform().

118  {
120  auto tkbuilder = static_cast<TkTransientTrackingRecHitBuilder const*>(theTrackerRecHitBuilder.product());
121 
122  for (auto hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
123  if ((*hit)->isValid()) {
124  if ((*hit)->geographicalId().det() == DetId::Tracker) {
125  result.emplace_back((**hit).cloneForFit(*tkbuilder->geometry()->idToDet((**hit).geographicalId())));
126  } else if ((*hit)->geographicalId().det() == DetId::Muon) {
127  if ((*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit) {
128  LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
129  continue;
130  }
131  result.push_back(theMuonRecHitBuilder->build(&**hit));
132  } else if ((*hit)->geographicalId().det() == DetId::Forward && (*hit)->geographicalId().subdetId() == FastTime) {
133  if (theMtdAvailable)
134  result.push_back(theMTDRecHitBuilder->build(&**hit));
135  else
136  throw cms::Exception("TrackTransformer") << "MTD hit encountered but MTD not available!";
137  }
138  }
139  }
140 
141  return result;
142 }
const bool theRPCInTheFit
edm::ESHandle< TransientTrackingRecHitBuilder > theMTDRecHitBuilder
edm::ESHandle< TransientTrackingRecHitBuilder > theMuonRecHitBuilder
edm::ESHandle< TransientTrackingRecHitBuilder > theTrackerRecHitBuilder
#define LogTrace(id)
virtual RecHitPointer build(const TrackingRecHit *p) const =0
build a tracking rechit from an existing rechit
trackingRecHit_iterator recHitsEnd() const
last iterator to RecHits
std::vector< ConstRecHitPointer > ConstRecHitContainer
T const * product() const
Definition: ESHandle.h:86
trackingRecHit_iterator recHitsBegin() const
first iterator to RecHits
const MagneticField* TrackTransformer::magneticField ( ) const
inline

the magnetic field

Definition at line 76 of file TrackTransformer.h.

Referenced by transform().

76 { return &*theMGField; }
edm::ESHandle< MagneticField > theMGField
edm::ESHandle<Propagator> const& TrackTransformer::propagator ( ) const
inlineprivate

Definition at line 115 of file TrackTransformer.h.

Referenced by transform().

115 { return thePropagator; }
edm::ESHandle< Propagator > thePropagator
std::unique_ptr<TrajectoryFitter> const& TrackTransformer::refitter ( ) const
inline

the refitter used to refit the reco::Track

Definition at line 85 of file TrackTransformer.h.

85 { return theFitter; }
std::unique_ptr< TrajectoryFitter > theFitter
void TrackTransformer::setServices ( const edm::EventSetup setup)
overridevirtual

set the services needed by the TrackTransformer

Implements TrackTransformerBase.

Definition at line 66 of file TrackTransformer.cc.

References TrajectorySmoother::clone(), TrajectoryFitter::clone(), watchdog::const, edm::EventSetup::get(), hitCloner, edm::ESHandleBase::isValid(), LogTrace, metname, edm::ESHandle< T >::product(), AlCaHLTBitMon_QueryRunRegistry::string, theCacheId_GTG, theCacheId_MG, theCacheId_TC, theCacheId_TRH, theFitter, theFitterName, theMGField, theMtdAvailable, theMTDRecHitBuilder, theMTDRecHitBuilderName, theMuonRecHitBuilder, theMuonRecHitBuilderName, thePropagator, thePropagatorName, theSmoother, theSmootherName, theTrackerRecHitBuilder, theTrackerRecHitBuilderName, and theTrackingGeometry.

Referenced by TrackingRecoMaterialAnalyser::analyze(), MuonKinkFinder::init(), MuonReSeeder::produce(), TrackerToMuonPropagator::produce(), CSCOverlapsAlignmentAlgorithm::run(), and GlobalTrajectoryBuilderBase::setEvent().

66  {
67  const std::string metname = "Reco|TrackingTools|TrackTransformer";
68 
71  setup.get<TrajectoryFitter::Record>().get(theFitterName, aFitter);
72  setup.get<TrajectoryFitter::Record>().get(theSmootherName, aSmoother);
73  theFitter = aFitter->clone();
74  theSmoother.reset(aSmoother->clone());
75 
76  unsigned long long newCacheId_TC = setup.get<TrackingComponentsRecord>().cacheIdentifier();
77 
78  if (newCacheId_TC != theCacheId_TC) {
79  LogTrace(metname) << "Tracking Component changed!";
80  theCacheId_TC = newCacheId_TC;
82  }
83 
84  // Global Tracking Geometry
85  unsigned long long newCacheId_GTG = setup.get<GlobalTrackingGeometryRecord>().cacheIdentifier();
86  if (newCacheId_GTG != theCacheId_GTG) {
87  LogTrace(metname) << "GlobalTrackingGeometry changed!";
88  theCacheId_GTG = newCacheId_GTG;
90  }
91 
92  // Magfield Field
93  unsigned long long newCacheId_MG = setup.get<IdealMagneticFieldRecord>().cacheIdentifier();
94  if (newCacheId_MG != theCacheId_MG) {
95  LogTrace(metname) << "Magnetic Field changed!";
96  theCacheId_MG = newCacheId_MG;
98  }
99 
100  // Transient Rechit Builders
101  unsigned long long newCacheId_TRH = setup.get<TransientRecHitRecord>().cacheIdentifier();
102  if (newCacheId_TRH != theCacheId_TRH) {
103  theCacheId_TRH = newCacheId_TRH;
104  LogTrace(metname) << "TransientRecHitRecord changed!";
110  }
111  theFitter->setHitCloner(&hitCloner);
112  theSmoother->setHitCloner(&hitCloner);
113 }
const std::string theTrackerRecHitBuilderName
const std::string metname
unsigned long long theCacheId_TRH
virtual TrajectorySmoother * clone() const =0
edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry
edm::ESHandle< TransientTrackingRecHitBuilder > theMTDRecHitBuilder
virtual std::unique_ptr< TrajectoryFitter > clone() const =0
edm::ESHandle< Propagator > thePropagator
unsigned long long theCacheId_MG
edm::ESHandle< TransientTrackingRecHitBuilder > theMuonRecHitBuilder
std::unique_ptr< TrajectoryFitter > theFitter
unsigned long long theCacheId_TC
const std::string thePropagatorName
const std::string theFitterName
edm::ESHandle< MagneticField > theMGField
unsigned long long theCacheId_GTG
edm::ESHandle< TransientTrackingRecHitBuilder > theTrackerRecHitBuilder
#define LogTrace(id)
const std::string theMuonRecHitBuilderName
TkClonerImpl hitCloner
const std::string theSmootherName
const std::string theMTDRecHitBuilderName
T get() const
Definition: EventSetup.h:73
bool isValid() const
Definition: ESHandle.h:44
T const * product() const
Definition: ESHandle.h:86
std::unique_ptr< TrajectorySmoother > theSmoother
std::unique_ptr<TrajectorySmoother> const& TrackTransformer::smoother ( ) const
inline

the smoother used to smooth the trajectory which came from the refitting step

Definition at line 88 of file TrackTransformer.h.

References HLT_2018_cff::track.

88 { return theSmoother; }
std::unique_ptr< TrajectorySmoother > theSmoother
edm::ESHandle<GlobalTrackingGeometry> TrackTransformer::trackingGeometry ( ) const
inline

the tracking geometry

Definition at line 79 of file TrackTransformer.h.

Referenced by checkRecHitsOrdering(), and transform().

79 { return theTrackingGeometry; }
edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry
vector< Trajectory > TrackTransformer::transform ( const reco::Track newTrack) const
overridevirtual

Convert a reco::Track into Trajectory.

Convert Tracks into Trajectories.

Implements TrackTransformerBase.

Definition at line 164 of file TrackTransformer.cc.

References getTransientRecHits(), magneticField(), metname, AlCaHLTBitMon_QueryRunRegistry::string, HLT_2018_cff::track, and trackingGeometry().

Referenced by TrackingRecoMaterialAnalyser::analyze(), GlobalTrajectoryBuilderBase::build(), CSCPairResidualsConstraint::dphidzFromTrack(), MuonKinkFinder::fillTrkKink(), MuonReSeeder::produce(), TrackerToMuonPropagator::produce(), and transform().

164  {
165  const std::string metname = "Reco|TrackingTools|TrackTransformer";
166 
168 
169  auto recHitsForReFit = getTransientRecHits(track);
170  return transform(track, recHitsForReFit);
171 }
const MagneticField * magneticField() const
the magnetic field
const std::string metname
TransientTrackingRecHit::ConstRecHitContainer getTransientRecHits(const reco::TransientTrack &track) const
std::vector< Trajectory > transform(const reco::Track &) const override
Convert a reco::Track into Trajectory.
edm::ESHandle< GlobalTrackingGeometry > trackingGeometry() const
the tracking geometry
vector< Trajectory > TrackTransformer::transform ( const reco::TrackRef track) const

Convert a reco::TrackRef into Trajectory.

Definition at line 115 of file TrackTransformer.cc.

References transform().

115 { return transform(*track); }
std::vector< Trajectory > transform(const reco::Track &) const override
Convert a reco::Track into Trajectory.
vector< Trajectory > TrackTransformer::transform ( const reco::TransientTrack track,
TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit 
) const

Convert a reco::TrackRef into Trajectory, refit with a new set of hits.

Convert Tracks into Trajectories with a given set of hits.

Definition at line 174 of file TrackTransformer.cc.

References alongMomentum, anyDirection, checkRecHitsOrdering(), RefitDirection::geometricalDirection(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), reco::Track::innerDetId(), reco::TransientTrack::innermostMeasurementState(), RefitDirection::insideOut, TrajectoryStateOnSurface::isValid(), LogTrace, metname, oppositeToMomentum, reco::Track::outerDetId(), reco::TransientTrack::outermostMeasurementState(), RefitDirection::outsideIn, AlCaHLTBitMon_ParallelJobs::p, position, Propagator::propagate(), RefitDirection::propagationDirection(), propagator(), Trajectory::push(), groupFilesInBlocks::reverse, SurveyInfoScenario_cff::seed, AlCaHLTBitMon_QueryRunRegistry::string, TrajectoryStateOnSurface::surface(), theDoPredictionsOnly, theFitter, theRefitDirection, theSmoother, reco::TransientTrack::track(), FastTrackerRecHitMaskProducer_cfi::trajectories, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

175  {
176  const std::string metname = "Reco|TrackingTools|TrackTransformer";
177 
178  if (recHitsForReFit.size() < 2)
179  return vector<Trajectory>();
180 
181  // 8 cases are foreseen:
182  // [RH = rec hit order, P = momentum dir, FD = fit direction. IO/OI = inside-out/outside-in, AM/OM = along momentum/opposite to momentum]
183  // (1) RH IO | P IO | FD AM ---> Start from IN
184  // (2) RH IO | P IO | FD OM ---> Reverse RH and start from OUT
185  // (3) RH IO | P OI | FD AM ---> Reverse RH and start from IN
186  // (4) RH IO | P OI | FD OM ---> Start from OUT
187  // (5) RH OI | P IO | FD AM ---> Reverse RH and start from IN
188  // (6) RH OI | P IO | FD OM ---> Start from OUT
189  // (7) RH OI | P OI | FD AM ---> Start from IN
190  // (8) RH OI | P OI | FD OM ---> Reverse RH and start from OUT
191  //
192  // *** Rules: ***
193  // -A- If RH-FD agree (IO-AM,OI-OM) do not reverse the RH
194  // -B- If FD along momentum start from innermost state, otherwise use outermost
195 
196  // Other special cases can be handled:
197  // (1 bis) RH IO | P IO | GFD IO => FD AM ---> Start from IN
198  // (2 bis) RH IO | P IO | GFD OI => FD OM ---> Reverse RH and start from OUT
199  // (3 bis) RH IO | P OI | GFD OI => FD AM ---> Reverse RH and start from OUT
200  // (4 bis) RH IO | P OI | GFD IO => FD OM ---> Start from IN
201  // (5 bis) RH OI | P IO | GFD IO => FD AM ---> Reverse RH and start from IN
202  // (6 bis) RH OI | P IO | GFD OI => FD OM ---> Start from OUT
203  // (7 bis) RH OI | P OI | GFD OI => FD AM ---> Start from OUT
204  // (8 bis) RH OI | P OI | GFD IO => FD OM ---> Reverse RH and start from IN
205  //
206  // *** Additional rule: ***
207  // -A0- If P and GFD agree, then FD is AM otherwise is OM
208  // -A00- rechit must be ordered as GFD in order to handle the case of cosmics
209  // -B0- The starting state is decided by GFD
210 
211  // Determine the RH order
213  checkRecHitsOrdering(recHitsForReFit); // FIXME change nome of the *type* --> RecHit order!
214  LogTrace(metname) << "RH order (0-insideOut, 1-outsideIn): " << recHitsOrder;
215 
217 
218  // Apply rule -A0-
219  if (propagationDirection == anyDirection) {
222  RefitDirection::GeometricalDirection p = (momentum.x() * position.x() > 0 || momentum.y() * position.y() > 0)
225 
227  LogTrace(metname) << "P (0-insideOut, 1-outsideIn): " << p;
228  LogTrace(metname) << "FD (0-OM, 1-AM, 2-ANY): " << propagationDirection;
229  }
230  // -A0-
231 
232  // Apply rule -A-
234  if ((recHitsOrder == RefitDirection::insideOut && propagationDirection == oppositeToMomentum) ||
235  (recHitsOrder == RefitDirection::outsideIn && propagationDirection == alongMomentum))
236  reverse(recHitsForReFit.begin(), recHitsForReFit.end());
237  }
238  // -A-
239  // Apply rule -A00-
240  else {
241  // reorder the rechit as defined in theRefitDirection.geometricalDirection();
242  if (theRefitDirection.geometricalDirection() != recHitsOrder)
243  reverse(recHitsForReFit.begin(), recHitsForReFit.end());
244  }
245  // -A00-
246 
247  // Apply rule -B-
249  unsigned int innerId = track.track().innerDetId();
251  if (propagationDirection == oppositeToMomentum) {
252  innerId = track.track().outerDetId();
253  firstTSOS = track.outermostMeasurementState();
254  }
255  } else { // if(theRefitDirection.propagationDirection() == anyDirection)
256  // Apply rule -B0-
258  innerId = track.track().outerDetId();
259  firstTSOS = track.outermostMeasurementState();
260  }
261  // -B0-
262  }
263  // -B-
264 
265  if (!firstTSOS.isValid()) {
266  LogTrace(metname) << "Error wrong initial state!" << endl;
267  return vector<Trajectory>();
268  }
269 
271 
272  if (recHitsForReFit.front()->geographicalId() != DetId(innerId)) {
273  LogTrace(metname) << "Propagation occured" << endl;
274  firstTSOS = propagator()->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
275  if (!firstTSOS.isValid()) {
276  LogTrace(metname) << "Propagation error!" << endl;
277  return vector<Trajectory>();
278  }
279  }
280 
281  if (theDoPredictionsOnly) {
282  Trajectory aTraj(seed, propagationDirection);
283  TrajectoryStateOnSurface predTSOS = firstTSOS;
284  for (auto const& hit : recHitsForReFit) {
285  predTSOS = propagator()->propagate(predTSOS, hit->det()->surface());
286  if (predTSOS.isValid())
287  aTraj.push(TrajectoryMeasurement(predTSOS, hit));
288  }
289  return vector<Trajectory>(1, aTraj);
290  }
291 
292  auto const& trajectories = theFitter->fit(seed, recHitsForReFit, firstTSOS);
293 
294  if (trajectories.empty()) {
295  LogTrace(metname) << "No Track refitted!" << endl;
296  return trajectories;
297  }
298 
299  auto const& trajectoryBW = trajectories.front();
300 
301  auto const& trajectoriesSM = theSmoother->trajectories(trajectoryBW);
302 
303  if (trajectoriesSM.empty()) {
304  LogTrace(metname) << "No Track smoothed!" << endl;
305  }
306 
307  return trajectoriesSM;
308 }
PropagationDirection propagationDirection() const
const std::string metname
const RefitDirection theRefitDirection
GeometricalDirection geometricalDirection() const
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
T y() const
Definition: PV3DBase.h:60
GlobalPoint globalPosition() const
PropagationDirection
TrajectoryStateOnSurface innermostMeasurementState() const
RefitDirection::GeometricalDirection checkRecHitsOrdering(TransientTrackingRecHit::ConstRecHitContainer const &) const
const SurfaceType & surface() const
std::unique_ptr< TrajectoryFitter > theFitter
unsigned int outerDetId() const
DetId of the detector on which surface the outermost state is located.
Definition: Track.h:79
TrajectoryStateOnSurface outermostMeasurementState() const
edm::ESHandle< Propagator > const & propagator() const
#define LogTrace(id)
Definition: DetId.h:17
const Track & track() const
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
const bool theDoPredictionsOnly
GlobalVector globalMomentum() const
static int position[264][3]
Definition: ReadPGInfo.cc:289
T x() const
Definition: PV3DBase.h:59
unsigned int innerDetId() const
DetId of the detector on which surface the innermost state is located.
Definition: Track.h:82
std::unique_ptr< TrajectorySmoother > theSmoother

Member Data Documentation

TkClonerImpl TrackTransformer::hitCloner
private

Definition at line 120 of file TrackTransformer.h.

Referenced by setServices().

unsigned long long TrackTransformer::theCacheId_GTG = 0
private

Definition at line 96 of file TrackTransformer.h.

Referenced by setServices().

unsigned long long TrackTransformer::theCacheId_MG = 0
private

Definition at line 97 of file TrackTransformer.h.

Referenced by setServices().

unsigned long long TrackTransformer::theCacheId_TC = 0
private

Definition at line 95 of file TrackTransformer.h.

Referenced by setServices().

unsigned long long TrackTransformer::theCacheId_TRH = 0
private

Definition at line 98 of file TrackTransformer.h.

Referenced by setServices().

const bool TrackTransformer::theDoPredictionsOnly
private

Definition at line 102 of file TrackTransformer.h.

Referenced by transform().

std::unique_ptr<TrajectoryFitter> TrackTransformer::theFitter
private

Definition at line 109 of file TrackTransformer.h.

Referenced by setServices(), and transform().

const std::string TrackTransformer::theFitterName
private

Definition at line 108 of file TrackTransformer.h.

Referenced by setServices().

edm::ESHandle<MagneticField> TrackTransformer::theMGField
private

Definition at line 106 of file TrackTransformer.h.

Referenced by setServices().

bool TrackTransformer::theMtdAvailable
private

Definition at line 126 of file TrackTransformer.h.

Referenced by getTransientRecHits(), and setServices().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformer::theMTDRecHitBuilder
private

Definition at line 127 of file TrackTransformer.h.

Referenced by getTransientRecHits(), and setServices().

const std::string TrackTransformer::theMTDRecHitBuilderName
private

Definition at line 125 of file TrackTransformer.h.

Referenced by setServices().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformer::theMuonRecHitBuilder
private

Definition at line 123 of file TrackTransformer.h.

Referenced by getTransientRecHits(), and setServices().

const std::string TrackTransformer::theMuonRecHitBuilderName
private

Definition at line 122 of file TrackTransformer.h.

Referenced by setServices().

edm::ESHandle<Propagator> TrackTransformer::thePropagator
private

Definition at line 116 of file TrackTransformer.h.

Referenced by setServices().

const std::string TrackTransformer::thePropagatorName
private

Definition at line 114 of file TrackTransformer.h.

Referenced by setServices().

const RefitDirection TrackTransformer::theRefitDirection
private

Definition at line 103 of file TrackTransformer.h.

Referenced by transform().

const bool TrackTransformer::theRPCInTheFit
private

Definition at line 100 of file TrackTransformer.h.

Referenced by getTransientRecHits().

std::unique_ptr<TrajectorySmoother> TrackTransformer::theSmoother
private

Definition at line 112 of file TrackTransformer.h.

Referenced by setServices(), and transform().

const std::string TrackTransformer::theSmootherName
private

Definition at line 111 of file TrackTransformer.h.

Referenced by setServices().

edm::ESHandle<TransientTrackingRecHitBuilder> TrackTransformer::theTrackerRecHitBuilder
private

Definition at line 119 of file TrackTransformer.h.

Referenced by getTransientRecHits(), and setServices().

const std::string TrackTransformer::theTrackerRecHitBuilderName
private

Definition at line 118 of file TrackTransformer.h.

Referenced by setServices().

edm::ESHandle<GlobalTrackingGeometry> TrackTransformer::theTrackingGeometry
private

Definition at line 105 of file TrackTransformer.h.

Referenced by setServices().