sparsfit.h File Reference

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)

Function Documentation

void fit_histo ( class TH1F *  histo,
Double_t *  params 
)
void fit_histo ( class TH1F *  histo  ) 

Definition at line 11 of file theFit.h.

References fit_histo().

Referenced by fit_histo().

00011                                   {
00012    Double_t params[9]={0};
00013    Double_t parerrs[9]={0};
00014    fit_histo(histo,params,parerrs,kTRUE);
00015 }

void fit_histo_spars0 ( class TH1F *  histo,
Double_t *  params 
)
void fit_histo_spars0 ( class TH1F *  histo  ) 
void fit_histo_spars10 ( class TH1F *  histo,
Double_t *  params 
)
void fit_histo_spars10 ( class TH1F *  histo  ) 
void fit_histo_spars20 ( class TH1F *  histo,
Double_t *  params 
)
Double_t my_gamma ( Double_t *  x,
Double_t *  par 
)

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 }

Double_t myfped_1gauss ( Double_t *  x,
Double_t *  par 
)

Definition at line 288 of file theFit.h.

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 }

Double_t mypefit ( Double_t *  x,
Double_t *  par 
)

Definition at line 249 of file theFit.h.

References MuELoss::a.

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 }

Double_t mypefit_noped ( Double_t *  x,
Double_t *  par 
)
Double_t mypefit_noped_ordyn ( Double_t *  x,
Double_t *  par 
)

Generated on 2 Nov 2017 for loon by  doxygen 1.6.1