double ffit(double *x, double *par) { double y = 0.; double x0 = x[0]/0.8954; y = par[0] + par[1]*x0 + par[2]*x0*x0 + par[3]*x0*x0*x0 + par[4]*x0*x0*x0*x0 + par[5]*x0*x0*x0*x0*x0 + par[6]*x0*x0*x0*x0*x0*x0; if (x0 > 116.) { y = par[7] + par[8]*x0 + par[9]*x0*x0; } return y; } void RjRb_dbratio_short() { gStyle->SetTitleAlign (13); gStyle->SetTitleX (0.02); gStyle->SetTitleSize (0.07,"t"); gStyle->SetTitleSize (0.07,"xyz"); gStyle->SetTitleOffset (1.5,"y"); gStyle->SetTitleOffset (0.9,"x"); gStyle->SetPadRightMargin (0.1); gStyle->SetPadLeftMargin (0.25); gStyle->SetPadBottomMargin(0.16); gStyle->SetPadTopMargin(0.05); double xmin = 10.; // allowed function range double xmax = 150.; // allowed function range TF1 *func = new TF1("func",ffit,xmin,xmax,10); func->SetParameter(0, 3.26609e-01); func->SetParameter(1, 6.40266e-02); func->SetParameter(2, -2.38273e-03); func->SetParameter(3, 4.62234e-05); func->SetParameter(4, -5.02656e-07); func->SetParameter(5, 2.87526e-09); func->SetParameter(6, -6.70040e-12); //--- func->SetParameter(7, 9.57575e-01+5.02e-4); func->SetParameter(8, 2.70604e-04); func->SetParameter(9, -7.08136e-06 ); func->SetName("RjRb_DMC"); func->Draw(); func->SetTitle(""); func->GetYaxis()->SetRangeUser(0.4,1.4); func->GetXaxis()->SetLimits(0.0,140.0); func->GetXaxis()->SetRangeUser(0.0,140.0); func->GetXaxis()->SetNdivisions(408); lgrid100 = new TLine(func->GetXaxis()->GetXmin(),1.00,func->GetXaxis()->GetXmax(),1.00); lgrid100->SetLineStyle(2); lgrid100->SetLineColor(1); lgrid100->SetLineWidth(1); lgrid100->Draw(); }