f2.SetParameter(4,0.9) f2.SetParameter(5,0.1) maxval = f.GetMaximumX() print maxval gr.Fit('f2','V') reject = False gr.Draw("APSAME") c.Update() fl = TF1("fl",fline,xmin,0.8,3) fl.SetParameters(f.GetParameters()) gr.GetListOfFunctions().Add(fl) gROOT.GetListOfFunctions().Remove(fl) fr = TF1("fr",fline,1.3,xmax,3) fr.SetParameters(f.GetParameters()) gr.GetListOfFunctions().Add(fr) gROOT.GetListOfFunctions().Remove(fr) gr.Draw("APSAME") c.Update() #raw_input() fnew = TF1('fnew', func ,xmin, xmax,6) fnew.Add(fl,-1)