DTParametrizedDriftAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalMuon_DTParametrizedDriftAlgo_H
#define RecoLocalMuon_DTParametrizedDriftAlgo_H

#include "RecoLocalMuon/DTRecHit/interface/DTRecHitBaseAlgo.h"
#include "FWCore/Utilities/interface/ESGetToken.h"

class MagneticField;
class IdealMagneticFieldRecord;

class DTParametrizedDriftAlgo : public DTRecHitBaseAlgo {
public:
  DTParametrizedDriftAlgo(const edm::ParameterSet& config, edm::ConsumesCollector);

  ~DTParametrizedDriftAlgo() override;

  // Operations

  void setES(const edm::EventSetup& setup) override;

  bool compute(const DTLayer* layer,
               const DTDigi& digi,
               LocalPoint& leftPoint,
               LocalPoint& rightPoint,
               LocalError& error) const override;

  bool compute(const DTLayer* layer,
               const DTRecHit1D& recHit1D,
               const float& angle,
               DTRecHit1D& newHit1D) const override;

  bool compute(const DTLayer* layer,
               const DTRecHit1D& recHit1D,
               const float& angle,
               const GlobalPoint& globPos,
               DTRecHit1D& newHit1D) const override;

private:
  // Interpolate parametrization function
  const bool interpolate;

  // Times below MinTime (ns) are considered as coming from previous BXs.
  const float minTime;

  // Times above MaxTime (ns) are considered as coming from following BXs
  const float maxTime;

  // Do the actual work.
  virtual bool compute(const DTLayer* layer,
                       const DTWireId& wireId,
                       const float digiTime,
                       const float& angle,
                       const GlobalPoint& globPos,
                       LocalPoint& leftPoint,
                       LocalPoint& rightPoint,
                       LocalError& error,
                       int step) const;

  // Interface to the method which does the actual work suited for 2nd and 3rd steps
  virtual bool compute(const DTLayer* layer,
                       const DTWireId& wireId,
                       const float digiTime,
                       const float& angle,
                       const GlobalPoint& globPos,
                       DTRecHit1D& newHit1D,
                       int step) const;

  // Switch on/off the verbosity
  const bool debug;

  // Pointer to the magnetic field (read from ES once per event)
  const MagneticField* magField;
  const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> magFieldToken_;
};
#endif