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)
