{
#include <fstream.h>
c1 = new TCanvas("c1","c1",0,0,700,500);
double mass[100], cross[100], dcross[100];
int n=0;
ifstream fin("h_gg_cross.dat");
while (fin>>mass[n]>>cross[n]>>dcross[n]) n++;
fin.close();
hgg = new TGraph(n,mass,cross);
hgg->SetMaximum(100.);
hgg->SetMinimum(1.);
gPad->SetLogy();
hgg->GetXaxis()->SetTitle("Higgs mass (GeV)");
hgg->GetYaxis()->SetTitle("Cross Section [pb]");
hgg->GetXaxis()->SetRangeUser(100,500);
c1->SetGrid();
hgg->SetLineColor(2);
hgg->Draw("AL");
leg=new TLegend(0.65,0.8,0.92,0.9,"pp #rightarrow H, 14 TeV");
leg->AddEntry(hgg,"gluon fusion","l");
leg->SetFillColor(10);
leg->SetBorderSize(0);
leg->Draw();
c1->Print("higlu.png");
}