Propagator.h

Go to the documentation of this file.

#ifndef PROPAGATOR_H
#define PROPAGATOR_H

#include "lardata/RecoObjects/Interactor.h"
#include "lardata/RecoObjects/KETrack.h"
#include "lardata/RecoObjects/KalmanLinearAlgebra.h"
namespace detinfo {
  class DetectorPropertiesData;
}

#include <memory>
#include <optional>

namespace trkf {

  class Propagator {
  public:
    enum PropDirection { FORWARD, BACKWARD, UNKNOWN };

    Propagator(detinfo::DetectorPropertiesData const& detProp,
               double tcut,
               bool doDedx,
               const std::shared_ptr<const Interactor>& interactor);

    virtual ~Propagator();

    // Accessors.

    double getTcut() const { return fTcut; }
    bool getDoDedx() const { return fDoDedx; }
    const std::shared_ptr<const Interactor>& getInteractor() const { return fInteractor; }

    // Virtual methods.

    virtual Propagator* clone() const = 0;

    virtual std::optional<double> short_vec_prop(KTrack& trk,
                                                 const std::shared_ptr<const Surface>& psurf,
                                                 PropDirection dir,
                                                 bool doDedx,
                                                 TrackMatrix* prop_matrix = 0,
                                                 TrackError* noise_matrix = 0) const = 0;

    virtual std::optional<double> origin_vec_prop(KTrack& trk,
                                                  const std::shared_ptr<const Surface>& porient,
                                                  TrackMatrix* prop_matrix = 0) const = 0;

    std::optional<double> vec_prop(KTrack& trk,
                                   const std::shared_ptr<const Surface>& psurf,
                                   PropDirection dir,
                                   bool doDedx,
                                   TrackMatrix* prop_matrix = 0,
                                   TrackError* noise_matrix = 0) const;

    std::optional<double> lin_prop(KTrack& trk,
                                   const std::shared_ptr<const Surface>& psurf,
                                   PropDirection dir,
                                   bool doDedx,
                                   KTrack* ref = 0,
                                   TrackMatrix* prop_matrix = 0,
                                   TrackError* noise_matrix = 0) const;
    std::optional<double> err_prop(KETrack& tre,
                                   const std::shared_ptr<const Surface>& psurf,
                                   PropDirection dir,
                                   bool doDedx,
                                   KTrack* ref = 0,
                                   TrackMatrix* prop_matrix = 0) const;

    std::optional<double> noise_prop(KETrack& tre,
                                     const std::shared_ptr<const Surface>& psurf,
                                     PropDirection dir,
                                     bool doDedx,
                                     KTrack* ref = 0) const;

    std::optional<double> dedx_prop(double pinv, double mass, double s, double* deriv = 0) const;

  private:
    detinfo::DetectorPropertiesData const& fDetProp;
    double fTcut;                                  
    bool fDoDedx;                                  
    std::shared_ptr<const Interactor> fInteractor; 
  };
}

#endif