CosmicMuonSmoother.h

Go to the documentation of this file.

#ifndef CosmicMuonProducer_CosmicMuonSmoother_H
#define CosmicMuonProducer_CosmicMuonSmoother_H

#include "TrackingTools/PatternTools/interface/TrajectorySmoother.h"

#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
#include "RecoMuon/TransientTrackingRecHit/interface/MuonTransientTrackingRecHit.h"
#include "TrackingTools/TransientTrackingRecHit/interface/TransientTrackingRecHit.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/PatternTools/interface/Trajectory.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "RecoMuon/CosmicMuonProducer/interface/CosmicMuonUtilities.h"

class Propagator;
class KFUpdator;
class MuonServiceProxy;
class Chi2MeasurementEstimator;

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}  // namespace edm

class Trajectory;
class TrajectoryMeasurement;

typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonRecHitContainer;
typedef TransientTrackingRecHit::ConstRecHitPointer ConstRecHitPointer;
typedef TransientTrackingRecHit::ConstRecHitContainer ConstRecHitContainer;
typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer ConstMuonRecHitContainer;

class CosmicMuonSmoother : public TrajectorySmoother {
public:
  CosmicMuonSmoother(const edm::ParameterSet&, const MuonServiceProxy* service);
  ~CosmicMuonSmoother() override;

  Trajectory trajectory(const Trajectory&) const override;

  TrajectoryContainer trajectories(const Trajectory& traj) const override {
    return TrajectorySmoother::trajectories(traj);
  }

  CosmicMuonSmoother* clone() const override { return new CosmicMuonSmoother(*this); }

  virtual TrajectoryContainer trajectories(const TrajectorySeed& seed,
                                           const ConstRecHitContainer& hits,
                                           const TrajectoryStateOnSurface& firstPredTsos) const;

  const Propagator* propagatorAlong() const { return &*theService->propagator(thePropagatorAlongName); }

  const Propagator* propagatorOpposite() const { return &*theService->propagator(thePropagatorOppositeName); }

  const KFUpdator* updator() const { return theUpdator; }

  const CosmicMuonUtilities* utilities() const { return theUtilities; }

  const Chi2MeasurementEstimator* estimator() const { return theEstimator; }

  std::vector<Trajectory> fit(const Trajectory&) const;
  std::vector<Trajectory> fit(const TrajectorySeed& seed,
                              const ConstRecHitContainer& hits,
                              const TrajectoryStateOnSurface& firstPredTsos) const;

  void setHitCloner(TkCloner const* hc) override {}

private:
  std::vector<Trajectory> smooth(const std::vector<Trajectory>&) const;
  std::vector<Trajectory> smooth(const Trajectory&) const;

  TrajectoryStateOnSurface initialState(const Trajectory&) const;

  void sortHitsAlongMom(ConstRecHitContainer& hits, const TrajectoryStateOnSurface&) const;

  const KFUpdator* theUpdator;
  const Chi2MeasurementEstimator* theEstimator;
  const CosmicMuonUtilities* theUtilities;

  const MuonServiceProxy* theService;

  std::string thePropagatorAlongName;
  std::string thePropagatorOppositeName;
  double theErrorRescaling;
  std::string category_;
};
#endif