{ 
ifstream cfile("Appearance.txt");
string dummyline;
for(int iline=1; iline<18; iline++){
getline(cfile,dummyline);
}
float a1, a2;
const int NPOINT = 83;
float energy[NPOINT];
float event[NPOINT];

for(int iline = 0; iline<NPOINT; iline++){
cfile >> a1 >> a2;
energy[iline] = a1;
event[iline] = a2;
cout<< "a1" <<  a1  << "a2"<<  a2  <<endl;
}
Int_t ci; //color
TGraph *eventrate = new TGraph(NPOINT,energy,event);
eventrate ->GetXaxis()->SetRangeUser(0,1.2);
eventrate ->GetYaxis()->SetRangeUser(0,250);
eventrate ->GetXaxis()->SetTitle("Reconstructed Energy E (GeV)");
eventrate ->GetYaxis()->SetTitle("Number of events");
eventrate ->SetLineWidth(3);
ci = TColor::GetColor("#D55E00");
eventrate->SetLineColor(ci);

eventrate->Draw("hist");
}