TrackBichsel.hh

Go to the documentation of this file.

#ifndef G_TRACK_BICHSEL_H
#define G_TRACK_BICHSEL_H
 
#include <array>
#include <vector>
 
#include "Garfield/FundamentalConstants.hh"
#include "Garfield/Track.hh"
 
namespace Garfield {

 
class TrackBichsel : public Track {
 public:
  struct Cluster {
    double x, y, z, t;
    double energy;
  };

  TrackBichsel() : TrackBichsel(nullptr) {}
  TrackBichsel(Sensor* sensor);
  virtual ~TrackBichsel() {}
 
  bool NewTrack(const double x0, const double y0, const double z0,
                const double t0, const double dx0, const double dy0,
                const double dz0) override;
  const std::vector<Cluster>& GetClusters() const { return m_clusters; }
 
  double GetClusterDensity() override;
  double GetStoppingPower() override;
 
  bool Initialise();
  bool ComputeCrossSection();
 
 private:
  constexpr static size_t NEnergyBins = 1250;
  std::array<double, NEnergyBins + 1> m_E;

  std::array<double, NEnergyBins> m_dfdE;
  std::array<double, NEnergyBins> m_eps1;
  std::array<double, NEnergyBins> m_eps2;
  std::array<double, NEnergyBins> m_int;
  std::array<double, NEnergyBins> m_k1;
 
  constexpr static size_t NCdfBins = 10000;
  std::array<double, NCdfBins> m_tab;

  double m_density;
  double m_conv = 0.0092456;
 
  bool m_initialised = false;

  double m_imfp = 0.;
  double m_dEdx = 0.;

  double m_speed = SpeedOfLight;
 
  std::vector<Cluster> m_clusters;
};
}  // namespace Garfield
 
#endif