-
- fPtPara = new TF1("Pt-Parametrization",fPtParaFunc,fPtMin,fPtMax,0);
- fYPara = new TF1("Y -Parametrization",fYParaFunc,fYMin,fYMax,0);
- TF1* PtPara = new TF1("Pt-Parametrization",fPtParaFunc,0,15,0);
- TF1* YPara = new TF1("Y -Parametrization",fYParaFunc,-6,6,0);
+ char name[256];
+ snprintf(name, 256, "pt-parameterisation for %s", GetName());
+
+ if (fPtPara) fPtPara->Delete();
+ fPtPara = new TF1(name, fPtParaFunc, fPtMin, fPtMax,0);
+ gROOT->GetListOfFunctions()->Remove(fPtPara);
+// Set representation precision to 10 MeV
+ Int_t npx= Int_t((fPtMax - fPtMin) / fDeltaPt);
+
+ fPtPara->SetNpx(npx);
+
+ snprintf(name, 256, "y-parameterisation for %s", GetName());
+ if (fYPara) fYPara->Delete();
+ fYPara = new TF1(name, fYParaFunc, fYMin, fYMax, 0);
+ gROOT->GetListOfFunctions()->Remove(fYPara);
+
+
+ snprintf(name, 256, "pt-for-%s", GetName());
+ TF1 ptPara(name ,fPtParaFunc, 0, 15, 0);
+ snprintf(name, 256, "y-for-%s", GetName());
+ TF1 yPara(name, fYParaFunc, -6, 6, 0);