//#include "AtlasStyle.C"


void CSvsE(const char *inputFile){

//gSystem->Load("AtlasStyle.C");
//  SetAtlasStyle();

  TCanvas *c1 = new TCanvas("c1","A Simple Graph Example",800,500);
  c1->SetGrid();
  //c1->SetLogx();
  //c1->SetLogy();
  c1->SetTickx();
  c1->SetTicky();

  TMultiGraph *multigraph = new TMultiGraph();
  multigraph->SetTitle("Cross-section vs CM Energy");
  multigraph->GetXaxis()->SetTitle("#sqrt{s}(GeV)");
//  multigraph->GetXaxis()->SetLimits(240,3210);
  multigraph->GetYaxis()->SetTitle("#sigma(pb)");
//  multigraph->GetYaxis()->SetLimits(-2,33);

  TLegend *legend = new TLegend(0.72,0.15,0.88,0.35,NULL,"brNDC");
legend->SetTextSize(0.03);

  TGraph *g1 = new TGraph(inputFile, "%lg %lg");
  g1->SetLineColor(217);
  g1->SetLineWidth(3);
  //g1->SetLineWidth(1504);
  g1->SetLineStyle(1);
  //g1->SetFillStyle(3004);
  TGraph *g2 = new TGraph(inputFile, "%lg %*lg %lg");
  g2->SetLineColor(217);
  g2->SetLineWidth(3);
  //g2->SetLineWidth(1500);
  g2->SetLineStyle(2);
  //g2->SetFillStyle(3005);
  TGraph *g3 = new TGraph(inputFile, "%lg %*lg %*lg %lg");
  g3->SetLineColor(217);
  g3->SetLineWidth(3);
  g3->SetLineStyle(3);
  
  multigraph->Add(g1,"L");
  multigraph->Add(g2,"L");
  multigraph->Add(g3,"L");

  multigraph->Draw("AP");

 TPaveText *pt = new TPaveText(400.,25.,1300.,27.5);

   pt->AddText("M=100 GeV");
   pt->SetTextColor(1);
   pt->SetTextSize(0.04);
   pt->SetShadowColor(0);
//   pt->SetLabel("Born equation");
   pt->Draw();

legend->AddEntry(g1,"Sample-1","l");
legend->AddEntry(g2,"Sample-2","l");
legend->AddEntry(g3,"Sample-3","l");
legend->Draw();

   gPad->Modified(); gPad->Update();
   multigraph->GetXaxis()->SetLimits(200., 3250.);
   multigraph->SetMinimum(-1.);
   multigraph->SetMaximum(35.);
   gPad->Modified(); gPad->Update();

c1->Print("CSvsE.png");

}