#include <TPaveText.h>
#include <TList.h>
#include <TFitResult.h>
-#include <../hist/hist/src/TF1Helper.h>
+//#include <../hist/hist/src/TF1Helper.h>
#include <AliLog.h>
// by the user in its macro
fHistDataPM->Fit(fFuncSigBack, fFitOpt.Data(), "", fFitMin, fFitMax);
TFitResultPtr pmFitPtr = fHistDataPM->Fit(fFuncSigBack, fFitOpt.Data(), "", fFitMin, fFitMax);
- TFitResult *pmFitResult = pmFitPtr.Get();
+ //TFitResult *pmFitResult = pmFitPtr.Get(); // used only with TF1Helper
fFuncSignal->SetParameters(fFuncSigBack->GetParameters());
fFuncBackground->SetParameters(fFuncSigBack->GetParameters()+fFuncSignal->GetNpar());
Double_t bknd = fFuncBackground->Eval(m);
Double_t ebknd = 0;
for(Int_t iPar=fFuncSignal->GetNpar(); iPar<fFuncSigBack->GetNpar(); iPar++) {
+/* TF1Helper problem on alien compilation
for(Int_t jPar=iPar; jPar<fFuncSigBack->GetNpar(); jPar++) {
TF1 gradientIpar("gradientIpar",
ROOT::TF1Helper::TGradientParFunction(iPar-fFuncSignal->GetNpar(),fFuncBackground),0,0,0);
gradientIpar.Eval(m)*gradientJpar.Eval(m)*
(iPar==jPar ? 1.0 : 2.0);
}
+*/
}
Double_t signal = pm-bknd;
Double_t error = TMath::Sqrt(epm*epm+ebknd);
fValues(0) = fFuncSignal->Integral(fIntMin, fIntMax)/fHistDataPM->GetBinWidth(1);
fErrors(0) = 0;
for(Int_t iPar=0; iPar<fFuncSignal->GetNpar(); iPar++) {
+/* TF1Helper problem on alien compilation
for(Int_t jPar=iPar; jPar<fFuncSignal->GetNpar(); jPar++) {
TF1 gradientIpar("gradientIpar",
ROOT::TF1Helper::TGradientParFunction(iPar,fFuncSignal),0,0,0);
gradientIpar.Integral(fIntMin,fIntMax)*gradientJpar.Integral(fIntMin,fIntMax)*
(iPar==jPar ? 1.0 : 2.0);
}
+*/
}
// background
fValues(1) = fFuncBackground->Integral(fIntMin, fIntMax)/fHistDataPM->GetBinWidth(1);
fErrors(1) = 0;
for(Int_t iPar=fFuncSignal->GetNpar(); iPar<fFuncSigBack->GetNpar(); iPar++) {
+/* TF1Helper problem on alien compilation
for(Int_t jPar=iPar; jPar<fFuncSigBack->GetNpar(); jPar++) {
TF1 gradientIpar("gradientIpar",
ROOT::TF1Helper::TGradientParFunction(iPar-fFuncSignal->GetNpar(),fFuncBackground),0,0,0);
gradientIpar.Integral(fIntMin, fIntMax)*gradientJpar.Integral(fIntMin, fIntMax)*
(iPar==jPar ? 1.0 : 2.0);
}
+*/
}
}
else {
Double_t epm = fHistDataPM->GetBinError(iBin);
Double_t epp = 0;
for(Int_t iPar=0; iPar<funcClonePP->GetNpar(); iPar++) {
+/* TF1Helper problem on alien compilation
for(Int_t jPar=iPar; jPar<funcClonePP->GetNpar(); jPar++) {
TF1 gradientIpar("gradientIpar",
ROOT::TF1Helper::TGradientParFunction(iPar,funcClonePP),0,0,0);
gradientIpar.Eval(m)*gradientJpar.Eval(m)*
(iPar==jPar ? 1.0 : 2.0);
}
+*/
}
Double_t emm = 0;
for(Int_t iPar=0; iPar<funcCloneMM->GetNpar(); iPar++) {
+/* TF1Helper problem on alien compilation
for(Int_t jPar=iPar; jPar<funcCloneMM->GetNpar(); jPar++) {
TF1 gradientIpar("gradientIpar",
ROOT::TF1Helper::TGradientParFunction(iPar,funcCloneMM),0,0,0);
gradientIpar.Eval(m)*gradientJpar.Eval(m)*
(iPar==jPar ? 1.0 : 2.0);
}
+*/
}
Double_t signal = pm-2.0*TMath::Sqrt(pp*mm);