DTParametrizedDriftAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalMuon_DTParametrizedDriftAlgo_H
#define RecoLocalMuon_DTParametrizedDriftAlgo_H

#include "RecoLocalMuon/DTRecHit/interface/DTRecHitBaseAlgo.h"

class MagneticField;

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

  ~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;


};
#endif