{
c1 = new TCanvas("c1","A Simple Graph with error bars",0,0,700,500);
//c1->SetGrid();
c1->GetFrame()->SetFillColor(21);
c1->GetFrame()->SetBorderSize(12);
gPad->SetLogx();
double q1[100], q2[100], q3[100], y1[100], y2[100], y3[100];
int k, i=0; int j=0; int m=0;
float q, a, ainv;
//read Standard model
ifstream fin("nosusy.dat");
while (fin >> k >> q >> a >> ainv ) {
if(k==1){ i++; q1[i] = q ;y1[i]=ainv;}
if(k==2){ j++; q2[j] = q ;y2[j]=ainv;}
if(k==3){ m++; q3[m] = q ;y3[m]=ainv;}
}
f1 = new TGraphErrors(i,q1,y1,0,0);
f2 = new TGraphErrors(j,q2,y2,0,0);
f3 = new TGraphErrors(m,q3,y3,0,0);
f3->SetMaximum(70); 
f3->SetMinimum(0);
f3->GetXaxis()->SetTitle("Interaction Energy (GeV)");
f3->GetYaxis()->SetTitle("1/#alpha _{i}"); 
f3->GetYaxis()->SetTitleOffset(1.0);
f3->Draw("AL");
f1->Draw("SAME L"); 
f2->Draw("SAME L"); 
//
TText *t2=new TText(1.E13,46,"no SUSY"); 
t2->SetTextSize(0.04); t2->Draw(); 
TLatex *u1=new TLatex(1.E5,57,"1/#alpha_{1}");
TLatex *u2=new TLatex(1.E5,35,"1/#alpha_{2}");
TLatex *u3=new TLatex(1.E5,19,"1/#alpha_{3}");
u1->SetTextAngle(-15); u1->SetTextSize(0.05); u1->Draw();
u2->SetTextAngle(+15); u2->SetTextSize(0.05); u2->Draw();
u3->SetTextAngle(30); u3->SetTextSize(0.05); u3->Draw();
c1->Print("guts_nosusy.png");
 }