Change of naming convention from Analyzer to RawAnalyser for classes that analyses...
[u/mrichter/AliRoot.git] / HLT / PHOS / AliHLTPHOSRawAnalyzer.cxx
1 #include "AliHLTPHOSRawAnalyzer.h"
2 #include <iostream>
3
4 using std::cout;
5 using std::endl;
6
7
8
9 AliHLTPHOSRawAnalyzer:: AliHLTPHOSRawAnalyzer():fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
10 {
11
12
13 }
14
15 AliHLTPHOSRawAnalyzer::~AliHLTPHOSRawAnalyzer()
16 {
17
18 }
19
20 AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(AliHLTPHOSRawAnalyzer const&):fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
21 {
22
23 }
24
25 /**
26 * Main constructor
27 * @param dataPtr Data array for wich a subarray will be taken to perform the fit
28 * @param fs the sampling frequency in entities of MHz. Needed in order to calculate physical time
29 **/
30 AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(double *dtaPtr, double fs):fFloatDataPtr(0), fSampleFrequency(10), fTau(2), fDTof(99999), fDAmpl(99999), n(99999)
31 {
32   fFloatDataPtr = dtaPtr;  
33   fSampleFrequency = fs;
34 } //end  
35
36
37 /**
38 * Attemps to level the basline to zero.
39 * The baseline will be calculated from the pretrigger samples and subtracted from the 
40 * data array. 
41 * If pretrigger samples are not present then the basline correction will be incorrect. 
42 * @param dataPtr array for wich to correct the basline 
43 * @param N the number of pretrigger samples used to calculate the baseline.
44 **/
45 void 
46 AliHLTPHOSRawAnalyzer::BaselineCorrection(double *dataPtr, int N)
47 {
48   fFloatDataPtr = dataPtr;  
49   n=N;
50   cout << "Baseline correction not yet implemeted" << endl;
51 } //end BaselineCorrection
52
53
54 /**
55 * Shifts the basline with the amount given by baselineValue
56 * If pretrigger samples are not present then the basline correction will be incorrect. 
57 * @param dataPtr array for wich to correct the basline 
58 * @param BaslineValue the basline value to subtract..
59 **/
60 void 
61 AliHLTPHOSRawAnalyzer::BaselineCorrection(double *dataPtr, double baselineValue)
62 {
63   fFloatDataPtr = dataPtr;   
64   printf("\nbaselineValue = %f\n", baselineValue);
65   cout << "Baseline correction not yet implemeted" << endl;
66 } //end BaslineCorrection
67
68
69 /**
70  * Gives the timing in entities of sample indexes
71  * Physical time is found by multiplying  with the sampling intervall (Ts).
72  **/
73 float
74 AliHLTPHOSRawAnalyzer::GetTiming()
75 {
76   return fDTof;
77 } //end GetTiming
78
79
80 /**
81  * Gives the time in entities of ADC channels (quantization levels).  
82  * Absolute enrgy is found by multiplying with offline calibration constants.
83  **/
84 float
85 AliHLTPHOSRawAnalyzer::GetEnergy()
86 {
87   return fDAmpl;
88 } //end GetEnergy
89
90
91 /**
92  * Set data array. Overrides data data array set in the constructor.
93  **/
94 void 
95 AliHLTPHOSRawAnalyzer::SetData(double *data)
96 {
97   //  double *dta;
98   //  dta = data;
99   fFloatDataPtr = data;
100   //  cout << "Set data not yet implemented" << endl;
101 }
102
103 void 
104 AliHLTPHOSRawAnalyzer::SetSampleFreq(double freq)
105 {
106
107 }
108
109 int 
110 AliHLTPHOSRawAnalyzer::FindStartIndex(double treshold)
111 {
112   printf("\ntreshold = %f \n", treshold);
113   cout << "Find Start index not yet implemented" << endl;
114   return 0;
115 } //end FindStartIndex
116
117
118 /**
119  * This function applies only to the Chi and Least mean square fit. An initial guess is made
120  * based on the average of the first 5 samples and the first value exeeding this value.
121  **/
122 void 
123 AliHLTPHOSRawAnalyzer::MakeInitialGuess()
124 {
125   cout << "Make initial guess not yet implemeted" << endl;
126 }
127
128
129 /**
130  * This function applies only to the Chi and Least mean square fit. An initial guess is made
131  * based on the average of the first 5 samples and the first value exeeding threshold + this value.
132  * @param treshold The index of the first value above treshold is ntaken to be the first value.
133  **/
134 void 
135 AliHLTPHOSRawAnalyzer::MakeInitialGuess(int treshold)
136 {
137   printf("\ntreshold = %d\n", treshold);
138   cout << "Make initial guess not yet implemeted" << endl;  
139 }