@@ -47,7 +47,7 @@ Double_t MassFitFunction::fitLD(Double_t *x, Double_t *par) {
47
47
}
48
48
Double_t MassFitFunction::fitParabola ( Double_t *x, Double_t *par ){
49
49
50
- Double_t chi = (x[0 ] - par[0 ]) / (par[1 ]+1 . );
50
+ Double_t chi = (x[0 ] - par[0 ]) / (par[1 ]+ 0.001 );
51
51
Double_t fV = chi*chi + par[2 ] ;
52
52
return fV ;
53
53
}
@@ -79,9 +79,18 @@ Double_t MassFitFunction::fitSG(Double_t *x, Double_t *par) {
79
79
80
80
Double_t MassFitFunction::fitData (Double_t *x, Double_t *par) {
81
81
82
- // return fitBW(x,par) + fitGS(x, &par[3]) ;
83
- return fitSG (x,par) + fitLD (x, &par[6 ]) + fitLD (x, &par[9 ]) ;
84
- // return fitSG(x,par) + fitLD(x, &par[6]) ;
82
+ Double_t gs = TMath::Gaus (x[0 ],par[1 ],par[2 ]);
83
+
84
+ Double_t A0 = log ( x[0 ] ) - par[4 ] ;
85
+ Double_t A1 = (-1 .*par[5 ]*A0*A0) ;
86
+ Double_t LG_Val = par[3 ]*exp ( A1 ) / x[0 ] ;
87
+
88
+ Double_t ld1_Val = TMath::Landau (x[0 ],par[6 ],par[7 ]) ;
89
+
90
+ Double_t sg_Val = gs + LG_Val + (par[8 ]*ld1_Val) ;
91
+ Double_t fitV = par[0 ]*sg_Val ;
92
+
93
+ return fitV ;
85
94
}
86
95
87
96
Double_t MassFitFunction::fitData1 (Double_t *x, Double_t *par) {
0 commit comments