TrackDegrade.hh

Go to the documentation of this file.

#ifndef G_TRACK_DEGRADE_H
#define G_TRACK_DEGRADE_H
 
#include <array>
#include <vector>
 
#include "Garfield/Track.hh"
 
namespace Garfield {
class Medium;

 
class TrackDegrade : public Track {
 public:
  struct Electron {
    double x = 0.;
    double y = 0.;
    double z = 0.;
    double t = 0.;
    double energy = 0.;
    double dx = 0.;
    double dy = 0.;
    double dz = 0.;
  };
 
  struct Excitation {
    double x = 0.;
    double y = 0.;
    double z = 0.;
    double t = 0.;
    double energy = 0.;
  };
 
  struct Cluster {
    double x, y, z, t;
    std::vector<Electron> deltaElectrons;
    std::vector<Electron> electrons;
    std::vector<Excitation> excitations;
  };

  TrackDegrade() : TrackDegrade(nullptr) {}
  TrackDegrade(Sensor* sensor);
  virtual ~TrackDegrade() {}
 
  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(Medium* medium, const bool verbose = false);

  void SetThresholdEnergy(const double eth);
  void StoreExcitations(const bool on = true, const double thr = 4.) {
    m_storeExcitations = on;
    m_ethrExc = std::max(thr, 1.e-3);
  }

  void EnableBremsstrahlung(const bool on = true) { m_bremsStrahlung = on; }
  void EnableFullCascade(const bool on = true) { m_fullCascade = on; }
 
  void SetParticle(const std::string& particle) override;
 
 protected:
  std::vector<Cluster> m_clusters;
 
  bool m_penning = true;
  bool m_bremsStrahlung = true;
  bool m_fullCascade = true;
  bool m_storeExcitations = false;
  // Energy threshold for storing excitations.
  double m_ethrExc = 4.;
  // Energy threshold for tracking electrons.
  double m_ethr = 2.;
 
  double m_pressure = -1.;
  double m_temperature = -1.;
  std::string m_mediumName = "";
  unsigned int m_nGas = 0;
 
  double m_dedx = -1.;
  double m_clusterDensity = -1.;
 
  std::array<double, 6> m_rPenning;
  std::array<double, 6> m_dPenning;
 
  std::pair<std::vector<Electron>, std::vector<Excitation> >
  TransportDeltaElectron(const double x0, const double y0, const double z0,
                         const double t0, const double e0, const double dx,
                         const double dy, const double dz);
 
  void SetupPenning(Medium* medium, std::array<double, 6>& rP,
                    std::array<double, 6>& dP);
  bool IsInside(const double x, const double y, const double z);
};
}  // namespace Garfield
 
#endif