GamR  0.0.0
GammaROOT
Loading...
Searching...
No Matches
Peak.hh
Go to the documentation of this file.
1#ifndef GAMR_TOOLKIT_PEAK_HH
2#define GAMR_TOOLKIT_PEAK_HH
3
4#include <algorithm>
5#include <vector>
6#include <map>
7#include <sstream>
8#include <string>
9
10/* ROOT */
11#include <Rtypes.h>
12#include <TBrowser.h>
13#include <TCanvas.h>
14#include <TH1.h>
15#include <TH2.h>
16#include <TNamed.h>
17
18#include <TF1.h>
19#include <TMath.h>
20
21#include "Gate.hh"
22
23namespace GamR {
28
31 namespace TK {
32
33 enum PeakType { Gaussian, StepGaussian, OneTailGaussian, OneTailStepGaussian, TwoTailGaussian, TwoTailStepGaussian};
34
35 //Peak: something with area, amplitude, a centroid, and a width.
36 class Peak {
37 public:
38 virtual double GetArea() { return 0; }
39 virtual double GetAreaError() { return 0; }
40 virtual double GetAmp() { return 0; }
41 virtual double GetAmpError() { return 0; }
42 virtual double GetCent() { return 0; }
43 virtual double GetCentError() { return 0; }
44 virtual double GetWidth() { return 0; }
45 virtual double GetWidthError() { return 0; }
46
47 virtual void SetAmp(double amp) {}
48 virtual void SetAmpError(double error) {}
49 virtual void SetCent(double cent) {}
50 virtual void SetCentError(double error) {}
51 virtual void SetWidth(double width) {}
52 virtual void SetWidthError(double error) {}
53 };
54
55 //Basic peak : no functional form
56 class BPeak : public Peak, public TNamed, public TAttLine {
57 private:
58 double Area;
59 double AreaError;
60 double Counts;
61 double CountsError;
62 double Amp;
63 double AmpError;
64 double Cent;
65 double CentError;
66 double Width;
67 double WidthError;
68
69 GamR::TK::Gate mPeak;
70 std::vector<GamR::TK::Gate > mBackgrounds;
71
72 TF1 *linear;
73 TF1 *backlow;
74 TF1 *backhigh;
75
76 double Lowest;
77 double Highest;
78
79 void Set(TH1 *hist, GamR::TK::Gate peak, TF1 *background); //actually calculates area
80
81 void Paint(Option_t *option = "same");
82 Int_t DistancetoPrimitive(Int_t px, Int_t py);
83 void ExecuteEvent(Int_t event, Int_t px, Int_t py);
84
85 //helper class
86 class LFit {
87 TH1D *back;
88 public:
89 LFit(TH1D *background) { back = background; }
90 double operator() (double *x, double *p) { return back->GetBinContent(back->FindBin(x[0])) + p[0] + x[0]*p[1]; }
91 };
92
93 public :
94 BPeak();
95 ~BPeak();
96 double GetArea() { return Area; }
97 double GetAreaError() { return AreaError; }
98 double GetCounts() { return Counts; }
99 double GetCountsError() { return CountsError; }
100 double GetAmp() { return Amp; }
101 double GetAmpError() { return AmpError; }
102 double GetCent() { return Cent; }
103 double GetCentError() { return CentError; }
104 double GetWidth() { return Width; }
105 double GetWidthError() { return WidthError; }
106 void Set(TH1 *hist, double x1, double y1, double x2, double y2);
107 void Set(TH1 *hist, GamR::TK::Gate peak, std::vector<GamR::TK::Gate > background, Option_t *foption="", Option_t *option="sl");
108 void Set(TH1 *hist, GamR::TK::Gate peak, GamR::TK::Gate background, Option_t *foption="", Option_t *option="sl");
109
110 void SetCont(TH1 *hist, GamR::TK::Gate peak, std::vector<GamR::TK::Gate > background, Option_t *foption="", Option_t *option=""); //for continuous
111 void SetPoint(TH1 *hist, GamR::TK::Gate peak, std::vector<GamR::TK::Gate > background, Option_t *foption="", Option_t *option=""); //for point
112 void SetStep(TH1 *hist, GamR::TK::Gate peak, std::vector<GamR::TK::Gate > background, Option_t *foption="", Option_t *option=""); //for step
113
114
115 void Set(TH1 *hist, Option_t *foption="", Option_t *option="sl");
116 void Set(TH2 *hist, Option_t *foption="", Option_t *option="sl");
117 void Print();
118
119 void Delete(Option_t *option) { this->~BPeak(); } // *MENU*
120
121 ClassDef(BPeak, 1);
122 };
123
124 //"Function Peak" - a peak with a pre-defined shape
125 class FPeak : public Peak {
126 protected :
127 //int iTails1;
128 //int iTails2;
129 //int iStep;
130
131 TF1 *fPeakFunc;
132 double fLow;
133 double fHigh;
134 PeakType fPeakType;
135 //double fArea;
136 //double fAreaError;
137
138 public :
139 FPeak() {}
140 FPeak(double Low, double High, PeakType Type) : fLow(Low), fHigh(High), fPeakType(Type) {}
141 virtual ~FPeak() {};
142 //Peak(int Tails1, int Tails2, int Step) { iTails1 = Tails1; iTails2 = Tails2; iStep = Step; }
143 TF1 *GetFunc() { return fPeakFunc; }
144 PeakType GetPeakType() { return fPeakType; }
145
146 virtual TF1 *GetStepFunc() { return NULL; }
147 virtual TF1 *GetGaussFunc() { return NULL; }
148 virtual void SetParameters(std::vector<double> params) {}
149 virtual void SetParErrors(std::vector<double> errors) {}
150
151 //void SetFunc(TF1 *func) { fPeakFunc = func; }
152 //void SetArea(double area) { fArea = area; }
153 //void SetAreaError(double area_error) { fAreaError = area_error; }
154 };
155
156 class GaussianPeak : public FPeak {
157 public :
158 GaussianPeak () {};
159 GaussianPeak (double Low, double High);
160 ~GaussianPeak() {};
161 double GetArea() { return sqrt(2 * TMath::Pi()) * fPeakFunc->GetParameter(0) *
162 fabs(fPeakFunc->GetParameter(2)); }
163 double GetAreaError() { return sqrt(2 * TMath::Pi() *
164 (pow(fPeakFunc->GetParError(0) * fPeakFunc->GetParameter(2), 2) +
165 pow(fPeakFunc->GetParError(2) * fPeakFunc->GetParameter(0), 2))); }
166 double GetAmp() { return fPeakFunc->GetParameter(0); }
167 double GetAmpError() { return fPeakFunc->GetParError(0); }
168 double GetCent() { return fPeakFunc->GetParameter(1); }
169 double GetCentError() { return fPeakFunc->GetParError(1); }
170 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
171 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
172
173 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
174 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
175 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
176 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
177 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
178 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
179 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
180 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
181
182 TF1 *GetStepFunc() { return NULL; }
183 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
184 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
185 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
186 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
187 return gaussfunc;
188 }
189 };
190
191 class StepGaussianPeak : public FPeak {
192 public:
193 StepGaussianPeak () {};
194 StepGaussianPeak (double Low, double High);
195 ~StepGaussianPeak() {};
196 double GetArea() { return sqrt(2 * TMath::Pi()) * fPeakFunc->GetParameter(0) *
197 fabs(fPeakFunc->GetParameter(2)); }
198 double GetAreaError() { return sqrt(2 * TMath::Pi() *
199 (pow(fPeakFunc->GetParError(0) * fPeakFunc->GetParameter(2), 2) +
200 pow(fPeakFunc->GetParError(2) * fPeakFunc->GetParameter(0), 2))); }
201 double GetAmp() { return fPeakFunc->GetParameter(0); }
202 double GetAmpError() { return fPeakFunc->GetParError(0); }
203 double GetCent() { return fPeakFunc->GetParameter(1); }
204 double GetCentError() { return fPeakFunc->GetParError(1); }
205 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
206 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
207
208 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
209 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
210 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
211 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
212 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
213 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
214 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
215 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
216
217 TF1 *GetStepFunc() { TF1 *stepfunc = new TF1("StepFunc", "[0]*([3]/100)*std::erfc((x-[1])/(sqrt(2)*[2]))", fLow, fHigh);
218 stepfunc->SetParameter(0, fPeakFunc->GetParameter(0));
219 stepfunc->SetParameter(1, fPeakFunc->GetParameter(1));
220 stepfunc->SetParameter(2, fPeakFunc->GetParameter(2));
221 stepfunc->SetParameter(3, fPeakFunc->GetParameter(3));
222 return stepfunc;
223 }
224
225 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
226 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
227 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
228 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
229 return gaussfunc;
230 }
231
232 };
233
234 class OneTailGaussianPeak : public FPeak {
235 public:
236 OneTailGaussianPeak ();
237 OneTailGaussianPeak (double Low, double High);
238 ~OneTailGaussianPeak();
239 double GetArea() { return fPeakFunc->Integral(fLow, fHigh); }
240 double GetAreaError();
241 double GetAmp() { return fPeakFunc->GetParameter(0); }
242 double GetAmpError() { return fPeakFunc->GetParError(0); }
243 double GetCent() { return fPeakFunc->GetParameter(1); }
244 double GetCentError() { return fPeakFunc->GetParError(1); }
245 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
246 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
247
248 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
249 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
250 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
251 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
252 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
253 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
254 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
255 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
256
257 TF1 *GetStepFunc() { return NULL; }
258 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*(1-[3]/100)*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
259 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
260 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
261 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
262 gaussfunc->SetParameter(3, fPeakFunc->GetParameter(3));
263 return gaussfunc;
264 }
265 };
266
267 class OneTailStepGaussianPeak : public FPeak {
268 public:
269 OneTailStepGaussianPeak () {};
270 OneTailStepGaussianPeak (double Low, double High);
271 ~OneTailStepGaussianPeak() {};
272 double GetArea();
273 double GetAreaError();
274 double GetAmp() { return fPeakFunc->GetParameter(0); }
275 double GetAmpError() { return fPeakFunc->GetParError(0); }
276 double GetCent() { return fPeakFunc->GetParameter(1); }
277 double GetCentError() { return fPeakFunc->GetParError(1); }
278 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
279 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
280
281 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
282 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
283 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
284 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
285 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
286 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
287 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
288 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
289
290 TF1 *GetStepFunc() { TF1 *stepfunc = new TF1("StepFunc", "[0]*([3]/100)*std::erfc((x-[1])/(sqrt(2)*[2]))", fLow, fHigh);
291 stepfunc->SetParameter(0, fPeakFunc->GetParameter(0));
292 stepfunc->SetParameter(1, fPeakFunc->GetParameter(1));
293 stepfunc->SetParameter(2, fPeakFunc->GetParameter(2));
294 stepfunc->SetParameter(3, fPeakFunc->GetParameter(5));
295 return stepfunc;
296 }
297 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*(1-[3]/100)*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
298 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
299 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
300 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
301 gaussfunc->SetParameter(3, fPeakFunc->GetParameter(3));
302 return gaussfunc;
303 }
304 };
305
306 class TwoTailGaussianPeak : public FPeak {
307 public:
308 TwoTailGaussianPeak () {};
309 TwoTailGaussianPeak (double Low, double High);
310 ~TwoTailGaussianPeak() {};
311 double GetArea() { return fPeakFunc->Integral(fLow, fHigh); }
312 double GetAreaError();
313 double GetAmp() { return fPeakFunc->GetParameter(0); }
314 double GetAmpError() { return fPeakFunc->GetParError(0); }
315 double GetCent() { return fPeakFunc->GetParameter(1); }
316 double GetCentError() { return fPeakFunc->GetParError(1); }
317 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
318 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
319
320 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
321 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
322 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
323 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
324 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
325 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
326 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
327 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
328
329 TF1 *GetStepFunc() { return NULL; }
330 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*(1 - [3]/100 - [4]/100)*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
331 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
332 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
333 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
334 gaussfunc->SetParameter(3, fPeakFunc->GetParameter(3));
335 gaussfunc->SetParameter(4, fPeakFunc->GetParameter(5));
336 return gaussfunc;
337 }
338 };
339
340 class TwoTailStepGaussianPeak : public FPeak {
341 public:
342 TwoTailStepGaussianPeak () {};
343 TwoTailStepGaussianPeak (double Low, double High);
344 ~TwoTailStepGaussianPeak() {};
345 double GetArea();
346 double GetAreaError();
347 double GetAmp() { return fPeakFunc->GetParameter(0); }
348 double GetAmpError() { return fPeakFunc->GetParError(0); }
349 double GetCent() { return fPeakFunc->GetParameter(1); }
350 double GetCentError() { return fPeakFunc->GetParError(1); }
351 double GetWidth() { return 2.3548*fPeakFunc->GetParameter(2); }
352 double GetWidthError() { return 2.3548*fPeakFunc->GetParError(2); }
353
354 void SetParameters(std::vector<double> params) { for (int i=0; i<(int)params.size(); ++i) { fPeakFunc->SetParameter(i, params[i]); } }
355 void SetParErrors(std::vector<double> errors) { for (int i=0; i<(int)errors.size(); ++i) { fPeakFunc->SetParError(i, errors[i]); } }
356 void SetAmp(double amp) { fPeakFunc->SetParameter(0, amp); }
357 void SetAmpError(double error) { fPeakFunc->SetParError(0, error); }
358 void SetCent(double cent) { fPeakFunc->SetParameter(1, cent); }
359 void SetCentError(double error) { fPeakFunc->SetParError(1, error); }
360 void SetWidth(double width) { fPeakFunc->SetParameter(2, width/2.3548); }
361 void SetWidthError(double error) { fPeakFunc->SetParError(2, error/2.3548); }
362
363 TF1 *GetStepFunc() { TF1 *stepfunc = new TF1("StepFunc", "[0]*([3]/100)*std::erfc((x-[1])/(sqrt(2)*[2]))", fLow, fHigh);
364 stepfunc->SetParameter(0, fPeakFunc->GetParameter(0));
365 stepfunc->SetParameter(1, fPeakFunc->GetParameter(1));
366 stepfunc->SetParameter(2, fPeakFunc->GetParameter(2));
367 stepfunc->SetParameter(3, fPeakFunc->GetParameter(7));
368 return stepfunc;
369 }
370 TF1 *GetGaussFunc() { TF1 *gaussfunc = new TF1("GaussFunc", "[0]*(1 - [3]/100 - [4]/100)*exp(-pow(x-[1],2)/(2*pow([2],2)))", fLow, fHigh);
371 gaussfunc->SetParameter(0, fPeakFunc->GetParameter(0));
372 gaussfunc->SetParameter(1, fPeakFunc->GetParameter(1));
373 gaussfunc->SetParameter(2, fPeakFunc->GetParameter(2));
374 gaussfunc->SetParameter(3, fPeakFunc->GetParameter(3));
375 gaussfunc->SetParameter(4, fPeakFunc->GetParameter(5));
376 return gaussfunc;
377 }
378 };
379
380 BPeak* FitBPeak(TCanvas *canvas = gPad->GetCanvas(), Option_t *foption="", Option_t *option="sl");
381 BPeak* ClickBPeak(TCanvas *canvas = gPad->GetCanvas());
382 }
383}
384
385#endif
386
py
void py(TVirtualPad *canvas)
Definition Abbreviations.cc:195
px
void px(TVirtualPad *canvas)
Definition Abbreviations.cc:194
GamR::TK::OneTailStepGaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.cc:766
GamR::TK::BPeak::GetAmp
double GetAmp()
Definition Peak.hh:100
GamR::TK::GaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:183
GamR::TK::TwoTailStepGaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:347
GamR::TK::Peak::GetWidth
virtual double GetWidth()
Definition Peak.hh:44
GamR::TK::OneTailStepGaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:279
GamR::TK::TwoTailGaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:320
GamR::TK::TwoTailStepGaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:351
GamR::TK::StepGaussianPeak::~StepGaussianPeak
~StepGaussianPeak()
Definition Peak.hh:195
GamR::TK::GaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:180
GamR::TK::OneTailGaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:252
GamR::TK::OneTailGaussianPeak::GetCent
double GetCent()
Definition Peak.hh:243
GamR::TK::Peak::GetAmpError
virtual double GetAmpError()
Definition Peak.hh:41
GamR::TK::StepGaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:205
GamR::TK::FPeak::~FPeak
virtual ~FPeak()
Definition Peak.hh:141
GamR::TK::FPeak::fPeakFunc
TF1 * fPeakFunc
Definition Peak.hh:131
GamR::TK::OneTailStepGaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:282
GamR::TK::GaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:169
GamR::TK::GaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:178
GamR::TK::OneTailGaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:248
GamR::TK::OneTailStepGaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:290
GamR::TK::StepGaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:212
GamR::TK::TwoTailGaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:321
GamR::TK::FPeak::GetStepFunc
virtual TF1 * GetStepFunc()
Definition Peak.hh:146
GamR::TK::GaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.hh:163
GamR::TK::Peak::GetCentError
virtual double GetCentError()
Definition Peak.hh:43
GamR::TK::OneTailGaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:242
GamR::TK::Peak::GetWidthError
virtual double GetWidthError()
Definition Peak.hh:45
GamR::TK::TwoTailGaussianPeak::GetCent
double GetCent()
Definition Peak.hh:315
GamR::TK::StepGaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:210
GamR::TK::BPeak::GetCent
double GetCent()
Definition Peak.hh:102
GamR::TK::OneTailStepGaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:275
GamR::TK::PeakType
PeakType
Definition Peak.hh:33
GamR::TK::BPeak::GetCentError
double GetCentError()
Definition Peak.hh:103
GamR::TK::Peak::SetAmp
virtual void SetAmp(double amp)
Definition Peak.hh:47
GamR::TK::OneTailStepGaussianPeak::OneTailStepGaussianPeak
OneTailStepGaussianPeak()
Definition Peak.hh:269
GamR::TK::TwoTailGaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:325
GamR::TK::OneTailStepGaussianPeak::GetCent
double GetCent()
Definition Peak.hh:276
GamR::TK::TwoTailGaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:317
GamR::TK::GaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:166
GamR::TK::TwoTailStepGaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:361
GamR::TK::OneTailGaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:257
GamR::TK::OneTailStepGaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:281
GamR::TK::FPeak::fHigh
double fHigh
Definition Peak.hh:133
GamR::TK::TwoTailGaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:322
GamR::TK::Peak::GetAreaError
virtual double GetAreaError()
Definition Peak.hh:39
GamR::TK::TwoTailStepGaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:355
GamR::TK::Peak::GetArea
virtual double GetArea()
Definition Peak.hh:38
GamR::TK::OneTailGaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:244
GamR::TK::BPeak::GetAreaError
double GetAreaError()
Definition Peak.hh:97
GamR::TK::TwoTailGaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:314
GamR::TK::Peak::GetCent
virtual double GetCent()
Definition Peak.hh:42
GamR::TK::TwoTailStepGaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:357
GamR::TK::GaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:179
GamR::TK::OneTailStepGaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:277
GamR::TK::OneTailGaussianPeak::GetArea
double GetArea()
Definition Peak.hh:239
GamR::TK::OneTailStepGaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:285
GamR::TK::BPeak::SetCont
void SetCont(TH1 *hist, GamR::TK::Gate peak, std::vector< GamR::TK::Gate > background, Option_t *foption="", Option_t *option="")
Definition Peak.cc:146
GamR::TK::OneTailGaussianPeak::OneTailGaussianPeak
OneTailGaussianPeak()
Definition Peak.cc:906
GamR::TK::TwoTailGaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:327
GamR::TK::TwoTailStepGaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:370
GamR::TK::GaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:177
GamR::TK::TwoTailGaussianPeak::GetArea
double GetArea()
Definition Peak.hh:311
GamR::TK::StepGaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:214
GamR::TK::TwoTailStepGaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:350
GamR::TK::GaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:167
GamR::TK::GaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:182
GamR::TK::TwoTailGaussianPeak::~TwoTailGaussianPeak
~TwoTailGaussianPeak()
Definition Peak.hh:310
GamR::TK::TwoTailStepGaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.cc:846
GamR::TK::OneTailStepGaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:284
GamR::TK::StepGaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:215
GamR::TK::StepGaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.hh:198
GamR::TK::OneTailGaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:250
GamR::TK::BPeak::SetStep
void SetStep(TH1 *hist, GamR::TK::Gate peak, std::vector< GamR::TK::Gate > background, Option_t *foption="", Option_t *option="")
Definition Peak.cc:372
GamR::TK::GaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:171
GamR::TK::OneTailStepGaussianPeak::GetArea
double GetArea()
Definition Peak.cc:753
GamR::TK::GaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:176
GamR::TK::TwoTailStepGaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:356
GamR::TK::GaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:175
GamR::TK::StepGaussianPeak::StepGaussianPeak
StepGaussianPeak()
Definition Peak.hh:193
GamR::TK::BPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:105
GamR::TK::OneTailGaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:249
GamR::TK::OneTailGaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:241
GamR::TK::Peak::SetWidth
virtual void SetWidth(double width)
Definition Peak.hh:51
GamR::TK::TwoTailStepGaussianPeak::TwoTailStepGaussianPeak
TwoTailStepGaussianPeak()
Definition Peak.hh:342
GamR::TK::StepGaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:217
GamR::TK::GaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:170
GamR::TK::OneTailGaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:245
GamR::TK::BPeak::GetCounts
double GetCounts()
Definition Peak.hh:98
GamR::TK::BPeak::ClassDef
ClassDef(BPeak, 1)
GamR::TK::StepGaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:204
GamR::TK::StepGaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:213
GamR::TK::TwoTailGaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:326
GamR::TK::BPeak::LFit::LFit
LFit(TH1D *background)
Definition Peak.hh:89
GamR::TK::FPeak::FPeak
FPeak()
Definition Peak.hh:139
GamR::TK::BPeak::Delete
void Delete(Option_t *option)
Definition Peak.hh:119
GamR::TK::OneTailGaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:253
GamR::TK::TwoTailStepGaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:363
GamR::TK::Peak::SetAmpError
virtual void SetAmpError(double error)
Definition Peak.hh:48
GamR::TK::OneTailGaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:251
GamR::TK::TwoTailGaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:330
GamR::TK::Peak::SetCentError
virtual void SetCentError(double error)
Definition Peak.hh:50
GamR::TK::TwoTailStepGaussianPeak::GetArea
double GetArea()
Definition Peak.cc:833
GamR::TK::OneTailGaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:246
GamR::TK::OneTailGaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.cc:723
GamR::TK::GaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:173
GamR::TK::BPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:101
GamR::TK::TwoTailStepGaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:360
GamR::TK::BPeak::GetWidth
double GetWidth()
Definition Peak.hh:104
GamR::TK::OneTailStepGaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:286
GamR::TK::BPeak::Print
void Print()
Definition Peak.cc:672
GamR::TK::OneTailGaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:258
GamR::TK::BPeak::BPeak
BPeak()
Definition Peak.cc:14
GamR::TK::GaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:174
GamR::TK::OneTailGaussianPeak::~OneTailGaussianPeak
~OneTailGaussianPeak()
Definition Peak.cc:907
GamR::TK::OneTailGaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:254
GamR::TK::OneTailStepGaussianPeak::SetAmp
void SetAmp(double amp)
Definition Peak.hh:283
GamR::TK::BPeak::LFit::operator()
double operator()(double *x, double *p)
Definition Peak.hh:90
GamR::TK::TwoTailGaussianPeak::GetStepFunc
TF1 * GetStepFunc()
Definition Peak.hh:329
GamR::TK::StepGaussianPeak::GetCent
double GetCent()
Definition Peak.hh:203
GamR::TK::Peak::SetWidthError
virtual void SetWidthError(double error)
Definition Peak.hh:52
GamR::TK::StepGaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:202
GamR::TK::FPeak::fPeakType
PeakType fPeakType
Definition Peak.hh:134
GamR::TK::OneTailStepGaussianPeak::~OneTailStepGaussianPeak
~OneTailStepGaussianPeak()
Definition Peak.hh:271
GamR::TK::TwoTailGaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:323
GamR::TK::StepGaussianPeak::SetParErrors
void SetParErrors(std::vector< double > errors)
Definition Peak.hh:209
GamR::TK::TwoTailGaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:313
GamR::TK::ClickBPeak
BPeak * ClickBPeak(TCanvas *canvas)
Definition Peak.cc:870
GamR::TK::OneTailStepGaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:297
GamR::TK::Peak::GetAmp
virtual double GetAmp()
Definition Peak.hh:40
GamR::TK::BPeak::GetCountsError
double GetCountsError()
Definition Peak.hh:99
GamR::TK::Peak::SetCent
virtual void SetCent(double cent)
Definition Peak.hh:49
GamR::TK::GaussianPeak::~GaussianPeak
~GaussianPeak()
Definition Peak.hh:160
GamR::TK::GaussianPeak::GaussianPeak
GaussianPeak()
Definition Peak.hh:158
GamR::TK::FPeak::GetGaussFunc
virtual TF1 * GetGaussFunc()
Definition Peak.hh:147
GamR::TK::TwoTailGaussianPeak::GetCentError
double GetCentError()
Definition Peak.hh:316
GamR::TK::GaussianPeak::GetArea
double GetArea()
Definition Peak.hh:161
GamR::TK::GaussianPeak::GetCent
double GetCent()
Definition Peak.hh:168
GamR::TK::TwoTailGaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:318
GamR::TK::TwoTailStepGaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:354
GamR::TK::StepGaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:201
GamR::TK::StepGaussianPeak::SetAmpError
void SetAmpError(double error)
Definition Peak.hh:211
GamR::TK::FPeak::GetPeakType
PeakType GetPeakType()
Definition Peak.hh:144
GamR::TK::OneTailStepGaussianPeak::GetWidth
double GetWidth()
Definition Peak.hh:278
GamR::TK::BPeak::SetPoint
void SetPoint(TH1 *hist, GamR::TK::Gate peak, std::vector< GamR::TK::Gate > background, Option_t *foption="", Option_t *option="")
Definition Peak.cc:227
GamR::TK::OneTailGaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:255
GamR::TK::TwoTailStepGaussianPeak::~TwoTailStepGaussianPeak
~TwoTailStepGaussianPeak()
Definition Peak.hh:344
GamR::TK::BPeak::GetArea
double GetArea()
Definition Peak.hh:96
GamR::TK::FPeak::SetParErrors
virtual void SetParErrors(std::vector< double > errors)
Definition Peak.hh:149
GamR::TK::TwoTailStepGaussianPeak::GetAmpError
double GetAmpError()
Definition Peak.hh:348
GamR::TK::FPeak::FPeak
FPeak(double Low, double High, PeakType Type)
Definition Peak.hh:140
GamR::TK::StepGaussianPeak::GetGaussFunc
TF1 * GetGaussFunc()
Definition Peak.hh:225
GamR::TK::StepGaussianPeak::GetArea
double GetArea()
Definition Peak.hh:196
GamR::TK::FPeak::SetParameters
virtual void SetParameters(std::vector< double > params)
Definition Peak.hh:148
GamR::TK::TwoTailStepGaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:358
GamR::TK::TwoTailGaussianPeak::TwoTailGaussianPeak
TwoTailGaussianPeak()
Definition Peak.hh:308
GamR::TK::StepGaussianPeak::SetParameters
void SetParameters(std::vector< double > params)
Definition Peak.hh:208
GamR::TK::TwoTailStepGaussianPeak::SetCentError
void SetCentError(double error)
Definition Peak.hh:359
GamR::TK::FPeak::GetFunc
TF1 * GetFunc()
Definition Peak.hh:143
GamR::TK::BPeak::~BPeak
~BPeak()
Definition Peak.cc:15
GamR::TK::StepGaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:206
GamR::TK::OneTailStepGaussianPeak::GetAmp
double GetAmp()
Definition Peak.hh:274
GamR::TK::TwoTailGaussianPeak::GetAreaError
double GetAreaError()
Definition Peak.cc:797
GamR::TK::FPeak::fLow
double fLow
Definition Peak.hh:132
GamR::TK::FitBPeak
BPeak * FitBPeak(TCanvas *canvas, Option_t *foption, Option_t *option)
Definition Peak.cc:859
GamR::TK::OneTailStepGaussianPeak::SetWidthError
void SetWidthError(double error)
Definition Peak.hh:288
GamR::TK::OneTailStepGaussianPeak::SetWidth
void SetWidth(double width)
Definition Peak.hh:287
GamR::TK::TwoTailStepGaussianPeak::GetCent
double GetCent()
Definition Peak.hh:349
GamR::TK::TwoTailStepGaussianPeak::GetWidthError
double GetWidthError()
Definition Peak.hh:352
GamR::TK::TwoTailGaussianPeak::SetCent
void SetCent(double cent)
Definition Peak.hh:324
GamR::TK::OneTailStepGaussian
@ OneTailStepGaussian
Definition Peak.hh:33
GamR::TK::StepGaussian
@ StepGaussian
Definition Peak.hh:33
GamR::TK::Gaussian
@ Gaussian
Definition Peak.hh:33
GamR::TK::TwoTailGaussian
@ TwoTailGaussian
Definition Peak.hh:33
GamR::TK::OneTailGaussian
@ OneTailGaussian
Definition Peak.hh:33
GamR::TK::TwoTailStepGaussian
@ TwoTailStepGaussian
Definition Peak.hh:33
GamR
Definition Gain.cc:19