SignalSimulation/Signal.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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#ifndef SIGNAL_H
#define SIGNAL_H

#include <string>

class Drift;
class Sensor;
class Efield;

class Signal {

 public:

  Signal();
  Signal(const char* name, const Sensor* sens, const Drift* hole, const Drift* elect, const Efield* efield, TF1* pulse);
  ~Signal();

  TF1* SignalFractionFunc() const;
  TF2* SignalFractionVsTimeFunc() const;
  TF2* CurrentSignalVsTimeFunc() const;
  TF1* TotalSignalFunc() const;
  TF1* WeightedPositionFunc() const;
  TF2* WeightedPositionVsTimeFunc() const;
  TF2* WeightedLADriftVsTimeFunc() const;
  TF1* ChargeBarycenterFunc() const;
  TF1* LADriftBarycenterFunc() const;

  void SaveFunctions(const int npsig=0, const int npcurrx=0, const int npcurry=0, const int nplax=0, const int nplay=0) const;
  std::string getName() const;

 private:

  double getCurrentSignal(const double z, const double t0) const;
  double getSignalFraction(const double z, const double t0) const;
  double getConvolutedLADrift(const double z, const double t0) const;

  double signalToBeConvoluted(const double x, const double x0, const double t0, const Drift* drift) const;
  double laDriftToBeConvoluted(const double x, const double x0, const double t0, const Drift* drift) const;

  double getHoleSignalToBeConvoluted(const double* x, const double* p) const;
  double getElectronSignalToBeConvoluted(const double* x, const double* p) const;
  double getHoleLADriftToBeConvoluted(const double* x, const double* p) const;
  double getElectronLADriftToBeConvoluted(const double* x, const double* p) const;

  double getSignalFraction1D(const double *x, const double *p) const;
  double getSignalFraction2D(const double *x, const double *p) const;
  double getCurrentSignal2D(const double *x, const double *p) const;
  double getTotalSignal(const double *x, const double *p) const;
  double getWeightedPosition(const double *x, const double *p) const;
  double getWeightedPosition2D(const double *x, const double *p) const;
  double getWeightedLADrift(const double *x, const double *p) const;
  double getWeightedLADrift2D(const double *x, const double *p) const;
  double getChargeBarycenter(const double *x, const double *p) const;
  double getLADriftBarycenter(const double *x, const double *p) const;

  std::string _name;
  const Drift* _holed;
  const Drift* _electd;
  const Sensor* _sens;
  const Efield* _efield;

  TF1* _signalFractionfunc;
  TF1* _signalFractionfuncaux;
  TF2* _signalFractionVsTimefunc;
  TF2* _currentSignalVsTimefunc;
  TF1* _totalSignalfunc;
  TF1* _weightedPositionfunc;
  TF1* _weightedPositionfuncaux;
  TF2* _weightedPositionVsTimefunc;
  TF1* _weightedLADriftfuncaux;
  TF2* _weightedLADriftVsTimefunc;
  TF1* _chargeBarycenterfunc;
  TF1* _laDriftBarycenterfunc;

  TF1* _holeToBeConvolutedfunc;
  TF1* _electToBeConvolutedfunc;
  TF1* _holeLADriftToBeConvolutedfunc;
  TF1* _electLADriftToBeConvolutedfunc;

  TF1* _decoPulse;
  bool _pulseowner;

};

#endif // SIGNAL_H