Possibility to retrieve the sigma and peak uncertainties

diff --git a/PWGHF/vertexingHF/AliHFMassFitter.cxx b/PWGHF/vertexingHF/AliHFMassFitter.cxx
index dd7b4f5f1767c5d26b0857c237dcdc81ebaf39cb..74c5bf68d253884c88d73bc4a0b05c2731899e5b 100644 (file)
--- a/PWGHF/vertexingHF/AliHFMassFitter.cxx
+++ b/PWGHF/vertexingHF/AliHFMassFitter.cxx
@@ -66,7 +66,9 @@ AliHFMassFitter::AliHFMassFitter() :
   ffactor(1),
   fntuParam(0),
   fMass(1.865),
+  fMassErr(0.01),
   fSigmaSgn(0.012),
+  fSigmaSgnErr(0.001),
   fRawYield(0.),
   fRawYieldErr(0.),
   fSideBands(0),
@@ -102,7 +104,9 @@ AliHFMassFitter::AliHFMassFitter (const TH1F *histoToFit, Double_t minvalue, Dou
  ffactor(1),
  fntuParam(0),
  fMass(1.865),
+ fMassErr(0.01),
  fSigmaSgn(0.012),
+ fSigmaSgnErr(0.001),
  fRawYield(0.),
  fRawYieldErr(0.),
  fSideBands(0),
@@ -156,7 +160,9 @@ AliHFMassFitter::AliHFMassFitter(const AliHFMassFitter &mfit):
   ffactor(mfit.ffactor),
   fntuParam(mfit.fntuParam),
   fMass(mfit.fMass),
+  fMassErr(mfit.fMassErr),
   fSigmaSgn(mfit.fSigmaSgn),
+  fSigmaSgnErr(mfit.fSigmaSgnErr),
   fRawYield(mfit.fRawYield),
   fRawYieldErr(mfit.fRawYieldErr),
   fSideBands(mfit.fSideBands),
@@ -214,7 +220,9 @@ AliHFMassFitter& AliHFMassFitter::operator=(const AliHFMassFitter &mfit){
   ffactor= mfit.ffactor;
   fntuParam= mfit.fntuParam;
   fMass= mfit.fMass;
+  fMassErr= mfit.fMassErr;
   fSigmaSgn= mfit.fSigmaSgn;
+  fSigmaSgnErr= mfit.fSigmaSgnErr;
   fRawYield= mfit.fRawYield;
   fRawYieldErr= mfit.fRawYieldErr;
   fSideBands= mfit.fSideBands;
@@ -1259,7 +1267,9 @@ Bool_t AliHFMassFitter::MassFitter(Bool_t draw){
   cout<<"fit done"<<endl;
   //reset value of fMass and fSigmaSgn to those found from fit
   fMass=funcmass->GetParameter(fNFinalPars-2);
+  fMassErr=funcmass->GetParError(fNFinalPars-2);
   fSigmaSgn=funcmass->GetParameter(fNFinalPars-1);
+  fSigmaSgnErr=funcmass->GetParError(fNFinalPars-1);
   fRawYield=funcmass->GetParameter(fNFinalPars-3)/fhistoInvMass->GetBinWidth(1);
   fRawYieldErr=funcmass->GetParError(fNFinalPars-3)/fhistoInvMass->GetBinWidth(1);
 

diff --git a/PWGHF/vertexingHF/AliHFMassFitter.h b/PWGHF/vertexingHF/AliHFMassFitter.h
index 6d665a89409a3f38db28d8377908e9abd873d7a8..00100e0bea75d9354e5ff0c6794f8bdb5d174cdf 100644 (file)
--- a/PWGHF/vertexingHF/AliHFMassFitter.h
+++ b/PWGHF/vertexingHF/AliHFMassFitter.h
@@ -61,7 +61,9 @@ class AliHFMassFitter : public TNamed {
   void     GetTypeOfFit(Bool_t &background, Int_t &typeb) const {background = fWithBkg; typeb = ftypeOfFit4Bkg;}
   Int_t    GetReflectionSigmaFactor() const {return ffactor;} 
   Double_t GetMean() const {return fMass;}
+  Double_t GetMeanUncertainty() const {return fMassErr;}
   Double_t GetSigma()const {return fSigmaSgn;}
+  Double_t GetSigmaUncertainty()const { return fSigmaSgnErr;}
   Double_t GetRawYield()const {return fRawYield;}
   Double_t GetRawYieldError()const {return fRawYieldErr;}
   Double_t GetChiSquare() const;
@@ -135,7 +137,9 @@ class AliHFMassFitter : public TNamed {
   Int_t     ffactor;           // number to multiply to the sigma of the signal to obtain the reflected gaussian
   TNtuple*  fntuParam;         // contains fit parameters
   Double_t  fMass;             // signal gaussian mean value
+  Double_t  fMassErr;          // err signal gaussian mean value
   Double_t  fSigmaSgn;         // signal gaussian sigma
+  Double_t  fSigmaSgnErr;      // err signal gaussian sigma
   Double_t  fRawYield;         // signal gaussian integral
   Double_t  fRawYieldErr;      // err on signal gaussian integral
   Bool_t    fSideBands;        // kTRUE = only side bands considered
@@ -145,7 +149,7 @@ class AliHFMassFitter : public TNamed {
   Int_t     fcounter;          // internal counter
   TList*    fContourGraph;     // TList of TGraph containing contour plots
 
-  ClassDef(AliHFMassFitter,5); // class for invariant mass fit
+  ClassDef(AliHFMassFitter,6); // class for invariant mass fit
 };
 
 #endif