TrackTransformer.h

Go to the documentation of this file.

#ifndef TrackingTools_TrackRefitter_TrackTransformer_H
#define TrackingTools_TrackRefitter_TrackTransformer_H

#include "TrackingTools/TrackRefitter/interface/TrackTransformerBase.h"

#include "TrackingTools/TrackRefitter/interface/RefitDirection.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/ESGetToken.h"

#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "MagneticField/Engine/interface/MagneticField.h"
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "RecoTracker/TransientTrackingRecHit/interface/TkTransientTrackingRecHitBuilder.h"

#include "DataFormats/TrackReco/interface/TrackFwd.h"

namespace reco {
  class TransientTrack;
}

class TrajectoryFitter;
class TrajectorySmoother;
class Propagator;
class TransientTrackingRecHitBuilder;
class Trajectory;
class GlobalTrackingGeometryRecord;
class IdealMagneticFieldRecord;
class TrajectoryFitterRecord;
class TrackingComponentsRecord;
class TransientRecHitRecord;

class TrackTransformer final : public TrackTransformerBase {
public:
  explicit TrackTransformer(const edm::ParameterSet&, edm::ConsumesCollector&);
  explicit TrackTransformer(const edm::ParameterSet& parameterSet, edm::ConsumesCollector&& iC)
      : TrackTransformer(parameterSet, iC) {}

  ~TrackTransformer() override;

  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");

  // Operations

  std::vector<Trajectory> transform(const reco::Track&) const override;

  std::vector<Trajectory> transform(const reco::TrackRef&) const;

  std::vector<Trajectory> transform(const reco::TransientTrack&, TransientTrackingRecHit::ConstRecHitContainer&) const;

  const MagneticField* magneticField() const { return &*theMGField; }

  edm::ESHandle<GlobalTrackingGeometry> trackingGeometry() const { return theTrackingGeometry; }

  void setServices(const edm::EventSetup&) override;

  std::unique_ptr<TrajectoryFitter> const& refitter() const { return theFitter; }

  std::unique_ptr<TrajectorySmoother> const& smoother() const { return theSmoother; }

  TransientTrackingRecHit::ConstRecHitContainer getTransientRecHits(const reco::TransientTrack& track) const;

private:
  RefitDirection::GeometricalDirection checkRecHitsOrdering(TransientTrackingRecHit::ConstRecHitContainer const&) const;

  unsigned long long theCacheId_TRH = 0;

  const bool theRPCInTheFit;

  const bool theDoPredictionsOnly;
  const RefitDirection theRefitDirection;

  edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> theTrackingGeometryToken;
  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> theMGFieldToken;
  edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
  edm::ESHandle<MagneticField> theMGField;

  const std::string theFitterName;
  edm::ESGetToken<TrajectoryFitter, TrajectoryFitterRecord> theFitterToken;
  std::unique_ptr<TrajectoryFitter> theFitter;

  const std::string theSmootherName;
  edm::ESGetToken<TrajectorySmoother, TrajectoryFitterRecord> theSmootherToken;
  std::unique_ptr<TrajectorySmoother> theSmoother;

  const std::string thePropagatorName;
  edm::ESGetToken<Propagator, TrackingComponentsRecord> thePropagatorToken;
  edm::ESHandle<Propagator> const& propagator() const { return thePropagator; }
  edm::ESHandle<Propagator> thePropagator;

  const std::string theTrackerRecHitBuilderName;
  edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> theTrackerRecHitBuilderToken;
  const TransientTrackingRecHitBuilder* theTrackerRecHitBuilder;
  TkClonerImpl hitCloner;

  const std::string theMuonRecHitBuilderName;
  edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> theMuonRecHitBuilderToken;
  edm::ESHandle<TransientTrackingRecHitBuilder> theMuonRecHitBuilder;

  const std::string theMTDRecHitBuilderName;
  edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> theMTDRecHitBuilderToken;
  bool theMtdAvailable;
  edm::ESHandle<TransientTrackingRecHitBuilder> theMTDRecHitBuilder;
};
#endif