|
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
|
#ifndef TMatacq_H
#define TMatacq_H
#include "TObject.h"
class TMatacq : public TObject {
public:
static constexpr int NMAXSAMP = 100;
static constexpr int NSPARAB = 16;
private:
int fNsamples;
int fNum_samp_bef_max;
int fNum_samp_aft_max;
int firstsample, lastsample, samplemax, presample, endsample;
int bing[101];
double nsigcut;
double level1, level2, level3;
double bong[NMAXSAMP];
double t[NSPARAB], val[NSPARAB];
double fv1[NSPARAB], fv2[NSPARAB], fv3[NSPARAB];
double bl, sigbl, val_max;
double ampl, timeatmax;
double pkval, sigpkval;
double trise;
double width20, width50, width80;
double meantrise, sigtrise;
int nevmtq0, nevmtq1, nevlasers;
int status[1200];
double comp_trise[1200], comp_peak[1200];
double slidingmean;
int nslide;
double interpolate(double);
public:
// Default Constructor, mainly for Root
TMatacq(int, int, int, int, int, int, int, int, int, int, int);
// Destructor: Does nothing
~TMatacq() override;
// Initialize
void init();
int rawPulseAnalysis(Int_t, Double_t*); // GHM
int findPeak();
int doFit();
int compute_trise();
void enterdata(Int_t);
int countBadPulses(Int_t);
void printmatacqData(Int_t, Int_t, Int_t);
void printitermatacqData(Int_t, Int_t, Int_t);
int getPt1() { return firstsample; }
int getPt2() { return lastsample; }
int getPtm() { return samplemax; }
double getBaseLine() { return bl; }
double getsigBaseLine() { return sigbl; }
double getTimpeak() { return pkval; }
double getsigTimpeak() { return sigpkval; }
double getAmpl() { return ampl; }
double getTimax() { return timeatmax; }
double getTrise() { return trise; }
double getFwhm() { return width50; }
double getWidth20() { return width20; }
double getWidth80() { return width80; }
double getSlide() { return slidingmean; }
ClassDefOverride(TMatacq, 0)
};
#endif
|