Track.hh

Go to the documentation of this file.

#ifndef G_TRACK_H
#define G_TRACK_H
 
#include <array>
#include <cmath>
#include <string>
 
namespace Garfield {
 
class Sensor;
class ViewDrift;

 
class Track {
 public:
  Track() = delete;
  Track(const std::string& name);
  virtual ~Track() {}

  virtual void SetParticle(const std::string& part);

  void SetEnergy(const double e);
  void SetBetaGamma(const double bg);
  void SetBeta(const double beta);
  void SetGamma(const double gamma);
  void SetMomentum(const double p);
  void SetKineticEnergy(const double ekin);

  double GetEnergy() const { return m_energy; }
  double GetBetaGamma() const { return sqrt(m_beta2 / (1. - m_beta2)); }
  double GetBeta() const { return sqrt(m_beta2); }
  double GetGamma() const { return sqrt(1. / (1. - m_beta2)); }
  double GetMomentum() const { return m_mass * sqrt(m_beta2 / (1. - m_beta2)); }
  double GetKineticEnergy() const { return m_energy - m_mass; }

  double GetCharge() const { return m_q; }
  double GetMass() const { return m_mass; }

  void SetSensor(Sensor* s);

  virtual bool NewTrack(const double x0, const double y0, const double z0,
                        const double t0, const double dx0, const double dy0,
                        const double dz0) = 0;

  virtual double GetClusterDensity() { return 0.; }
  virtual double GetStoppingPower() { return 0.; }

  void EnablePlotting(ViewDrift* viewer);
  void DisablePlotting();

  void EnableDebugging() { m_debug = true; }
  void DisableDebugging() { m_debug = false; }
 
 protected:
  std::string m_className = "Track";
 
  double m_q = -1.;
  int m_spin = 1;
  double m_mass;
  double m_energy = 0.;
  double m_beta2 = 1.;
  bool m_isElectron = false;
  std::string m_particleName = "mu-";
 
  Sensor* m_sensor = nullptr;
 
  bool m_isChanged = true;
 
  ViewDrift* m_viewer = nullptr;
 
  bool m_debug = false;
 
  size_t m_plotId = 0;
  void PlotNewTrack(const double x0, const double y0, const double z0);
  void PlotCluster(const double x0, const double y0, const double z0);
 
  static std::array<double, 3> StepBfield(const double dt, const double qoverm,
                                          const double vmag, double bx,
                                          double by, double bz,
                                          std::array<double, 3>& dir);
};
}  // namespace Garfield
 
#endif