Go to the source code of this file.
Functions | |
| void | fit_histo_spars0 (class TH1F *histo) |
| void | fit_histo_spars10 (class TH1F *histo) |
| void | fit_histo (class TH1F *histo) |
| void | fit_histo (class TH1F *histo, Double_t *params) |
| void | fit_histo_spars0 (class TH1F *histo, Double_t *params) |
| void | fit_histo_spars10 (class TH1F *histo, Double_t *params) |
| void | fit_histo_spars20 (class TH1F *histo, Double_t *params) |
| Double_t | mypefit_noped (Double_t *x, Double_t *par) |
| Double_t | mypefit_noped_ordyn (Double_t *x, Double_t *par) |
| Double_t | mypefit (Double_t *x, Double_t *par) |
| Double_t | myfped_1gauss (Double_t *x, Double_t *par) |
| Double_t | my_gamma (Double_t *x, Double_t *par) |
|
||||||||||||
|
|
|
|
Definition at line 11 of file theFit.h. References fit_histo(). 00011 {
00012 Double_t params[9]={0};
00013 Double_t parerrs[9]={0};
00014 fit_histo(histo,params,parerrs,kTRUE);
00015 }
|
|
||||||||||||
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
||||||||||||
|
|
|
||||||||||||
|
Definition at line 301 of file theFit.h. 00302 {
00303 Double_t top= pow(x[0],par[1]-1)*exp(-x[0]/par[2]);
00304 Double_t bottom= pow(par[2],par[1])*TMath::Gamma(par[1]);
00305
00306
00307
00308 return par[0]*top /bottom;
00309 }
|
|
||||||||||||
|
Definition at line 288 of file theFit.h. References s(). Referenced by fit_histo(). 00289 {
00290 static double r2pi = 2.506628275;
00291 static double r2 = 1.414213562;
00292
00293 Double_t s = r2*par[1];
00294 Double_t y = 0.5*par[2]*(TMath::Erf(Double_t(x[0]-par[0]+.5)/s)
00295 -TMath::Erf(Double_t(x[0]-par[0]-.5)/s));
00296 return y;
00297
00298 }
|
|
||||||||||||
|
Definition at line 249 of file theFit.h. References s(). Referenced by fit_histo(). 00250 {
00251 static double r2pi = 2.506628275;
00252 static double r2 = 1.414213562;
00253
00254 double a,s,y;
00255
00256 /* Pedestal p.d.f.: double Gaussian. */
00257 a = par[7];
00258 a = a*exp(-par[1]);
00259 s = r2*par[5];
00260 y = 0.5*a*(TMath::Erf(Double_t(x[0]-par[0]+.5)/s)
00261 -TMath::Erf(Double_t(x[0]-par[0]-.5)/s));
00262
00263
00264 /*Photoelectron p.d.f.: combination of Poisson and Gaussian.*/
00265 a=par[7]*exp(-par[1])*(1-par[4]);
00266 a = a/r2pi;
00267 s=par[5];
00268 for (int j=1; j<12; j++)
00269 {
00270 a = a*par[1]/double(j);
00271 s = sqrt(pow(par[3]*par[2],2)+pow(s,2));
00272 y += a/s*exp(-.5*pow(((double(x[0])-double(j)*par[2]-par[0])/s),2));
00273 }
00274
00275
00276 /* Dynode PE: assumed to be a gaussian with the mean and the width
00277 1/10 of that of single PE. */
00278 a = par[7]*par[4]*(1-exp(-par[1]));
00279 s = par[5];
00280 s = sqrt( pow(par[3]*par[2]/par[6],2) + pow(s,2) );
00281
00282 y += a/(s*r2pi)*exp(-.5*pow(((double(x[0])-par[2]/par[6]-par[0])/s),2));
00283 return y;
00284
00285 }
|
|
||||||||||||
|
|
|
||||||||||||
|
|
1.3.9.1