]>
Commit | Line | Data |
---|---|---|
1 | /************************************************************************** | |
2 | * Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: Per Thomas Hille for the ALICE HLT Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | #include "AliHLTPHOSRawAnalyzer.h" | |
17 | ||
18 | ||
19 | AliHLTPHOSRawAnalyzer:: AliHLTPHOSRawAnalyzer(): AliHLTPHOSBase(), | |
20 | fFloatDataPtr(0), | |
21 | fIntDataPtr(0), | |
22 | fSampleFrequency(10), | |
23 | fDTofGuess(-1), | |
24 | fDAmplGuess(-1), | |
25 | fTau(2), | |
26 | fDTof(99999), | |
27 | fDAmpl(99999), | |
28 | fStartIndex(0) | |
29 | { | |
30 | fIntDataPtr = new UInt_t[1008]; | |
31 | } | |
32 | ||
33 | AliHLTPHOSRawAnalyzer::~AliHLTPHOSRawAnalyzer() | |
34 | { | |
35 | delete[] fIntDataPtr; | |
36 | } | |
37 | ||
38 | ||
39 | ||
40 | /** | |
41 | * Main constructor | |
42 | * @param dtaPtr Data array for wich a subarray will be taken to perform the fit | |
43 | * @param fs the sampling frequency in entities of MHz. Needed in order to calculate physical time | |
44 | **/ | |
45 | AliHLTPHOSRawAnalyzer::AliHLTPHOSRawAnalyzer(double * /*dtaPtr*/, double fs): AliHLTPHOSBase(), | |
46 | fFloatDataPtr(0), | |
47 | fIntDataPtr(0), | |
48 | fSampleFrequency(10), | |
49 | fDTofGuess(-1), | |
50 | fDAmplGuess(-1), | |
51 | fTau(2), | |
52 | fDTof(99999), | |
53 | fDAmpl(99999), | |
54 | fStartIndex(0) | |
55 | { | |
56 | fSampleFrequency = fs; | |
57 | } //end | |
58 | ||
59 | ||
60 | /** | |
61 | * Attemps to level the basline to zero. | |
62 | * The baseline will be calculated from the pretrigger samples and subtracted from the | |
63 | * data array. | |
64 | * If pretrigger samples are not present then the basline correction will be incorrect. | |
65 | * @param dataPtr array for wich to correct the basline | |
66 | * @param N the number of pretrigger samples used to calculate the baseline. | |
67 | **/ | |
68 | void | |
69 | AliHLTPHOSRawAnalyzer::BaselineCorrection(double * /*dataPtr*/, int /*N*/) | |
70 | { | |
71 | ||
72 | } //end BaselineCorrection | |
73 | ||
74 | ||
75 | /** | |
76 | * Shifts the baseline with the amount given by baselineValue | |
77 | * If pretrigger samples are not present then the basline correction will be incorrect. | |
78 | * @param dataPtr array for wich to correct the basline | |
79 | * @param baselineValue the basline value to subtract.. | |
80 | **/ | |
81 | void | |
82 | AliHLTPHOSRawAnalyzer::BaselineCorrection(double * /*dataPtr*/, double baselineValue) | |
83 | { | |
84 | ||
85 | } //end BaslineCorrection | |
86 | ||
87 | ||
88 | /** | |
89 | * Gives the timing in entities of sample indexes | |
90 | * Physical time is found by multiplying with the sampling intervall (Ts). | |
91 | **/ | |
92 | float | |
93 | AliHLTPHOSRawAnalyzer::GetTiming() const | |
94 | { | |
95 | return fDTof; | |
96 | } //end GetTiming | |
97 | ||
98 | ||
99 | /** | |
100 | * Gives the time in entities of ADC channels (quantization levels). | |
101 | * Absolute enrgy is found by multiplying with offline calibration constants. | |
102 | **/ | |
103 | float | |
104 | AliHLTPHOSRawAnalyzer::GetEnergy() const | |
105 | { | |
106 | return fDAmpl; | |
107 | } //end GetEnergy | |
108 | ||
109 | ||
110 | /** | |
111 | * Set data array. Overrides data data array set in the constructor. | |
112 | **/ | |
113 | void | |
114 | AliHLTPHOSRawAnalyzer::SetData(const UInt_t *data) | |
115 | { | |
116 | fIntDataPtr = data; | |
117 | } | |
118 | ||
119 | /** | |
120 | * Set data array. Overrides data data array set in the constructor. | |
121 | **/ | |
122 | void | |
123 | AliHLTPHOSRawAnalyzer::SetData(const double *data) | |
124 | { | |
125 | fFloatDataPtr = data; | |
126 | } | |
127 | ||
128 | ||
129 | ||
130 | void | |
131 | AliHLTPHOSRawAnalyzer::SetSampleFreq(double freq) | |
132 | { | |
133 | fSampleFrequency = freq; | |
134 | } | |
135 | ||
136 | int | |
137 | AliHLTPHOSRawAnalyzer::FindStartIndex(double /*treshold*/) | |
138 | { | |
139 | // cout << "Find Start index not yet implemented" << endl; | |
140 | return 0; | |
141 | } //end FindStartIndex | |
142 | ||
143 | ||
144 | /** | |
145 | * This function applies only to the Chi and Least mean square fit. An initial guess is made | |
146 | * based on the average of the first 5 samples and the first value exeeding this value. | |
147 | **/ | |
148 | void | |
149 | AliHLTPHOSRawAnalyzer::MakeInitialGuess() | |
150 | { | |
151 | // cout << "Make initial guess not yet implemeted" << endl; | |
152 | } | |
153 | ||
154 | ||
155 | /** | |
156 | * This function applies only to the Chi and Least mean square fit. An initial guess is made | |
157 | * based on the average of the first 5 samples and the first value exeeding threshold + this value. | |
158 | * @param treshold The index of the first value above treshold is ntaken to be the first value. | |
159 | **/ | |
160 | void | |
161 | AliHLTPHOSRawAnalyzer::MakeInitialGuess(int /*treshold*/) | |
162 | { | |
163 | // cout << "Make initial guess not yet implemeted" << endl; | |
164 | } | |
165 | ||
166 | ||
167 | void | |
168 | AliHLTPHOSRawAnalyzer::SetStartIndex(int index) | |
169 | { | |
170 | fStartIndex = index; | |
171 | } | |
172 | ||
173 | ||
174 | ||
175 | void | |
176 | AliHLTPHOSRawAnalyzer::SetTVector(Double_t * /*tVector*/, Int_t /*size*/) | |
177 | { | |
178 | // cout <<"ERROR: AliHLTPHOSRawAnalyzer::SetTVector: You cannot set the peakfindervector here, must be set in derived class peakfinder"<<endl; | |
179 | } | |
180 | ||
181 | ||
182 | ||
183 | void | |
184 | AliHLTPHOSRawAnalyzer::SetAVector(Double_t * /*aVector*/, Int_t /*size*/) | |
185 | { | |
186 | // cout <<"ERROR: AliHLTPHOSRawAnalyzer::SetAVector: You cannot set the peakfindervector here, must be set in derived class peakfinder"<<endl; | |
187 | } | |
188 | ||
189 |