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();    
}