SignalSimulation/Drift.h at master · venturia/SignalSimulation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#ifndef DRIFT_H
#define DRIFT_H

#include <string>

class Mobility;
class TF1;
class TF2;
class Sensor;

class Drift {

 public:

  Drift();
  //  Drift(const char* name, const Mobility* mob, const TF1* efield, const double& bfield, const double thickness);
  Drift(const char* name, const Mobility* mob, const Sensor* sens);
  ~Drift();

  void SetLineColor(Color_t color);
  void SetLineStyle(Style_t style);

  TF1* MobilityFunc() const;
  TF1* tanLAFunc() const;
  TF1* parallelVdriftFunc() const;
  TF1* invParallelVdriftFunc() const;
  TF1* DriftTimeFunc() const;
  TF1* FromDriftTimeFunc() const;
  TF1* ToDriftTimeFunc() const;
  TF2* GenericDriftTimeFunc() const;
  TF2* GenericLorentzDriftFunc() const;

 private:

  double getEfield(const double *x, const double *p) const;
  double getMobility(const double *x, const double *p) const;
  double getTanLA(const double *x, const double *p) const;
  double getParallelVdrift(const double *x, const double *p) const;
  double getInvParallelVdrift(const double *x, const double *p) const;
  double getDriftTime(const double *x, const double *p) const;
  double getFromDriftTime(const double *x, const double *p) const;
  double getToDriftTime(const double *x, const double *p) const;
  double getGenericDriftTime(const double *x, const double *p) const;
  double getGenericLorentzDrift(const double *x, const double *p) const;

  double genericDriftTime(const double x, const double x0) const;
  double genericLorentzDrift(const double x, const double x0) const;

  std::string _name;
  const Mobility* _mob;

  const Sensor* _sens;

  TF1* _mobilityfunc;
  TF1* _tanLAfunc;
  TF1* _parallelVdriftfunc;
  TF1* _invParallelVdriftfunc;
  TF1* _driftTimefunc;
  TF1* _fromDriftTimefunc;
  TF1* _toDriftTimefunc;
  TF2* _genericDriftTimefunc;
  TF2* _genericLorentzDriftfunc;

};


#endif // DRIFT_H