void fluxbon()
  {
  TCanvas *c1 = new TCanvas("c1","c1",1000,10,1000,400);


  Int_t n = 5;
  Double_t x1[5] = {1, 10,100,1000,10000};
  
  Double_t y1[5] = {7.95719*pow(10,-61)*pow(10,30),2.66728*pow(10,-63)*pow(10,30),1.91568*pow(10,-66)*pow(10,30),5.62685*pow(10,-76)*pow(10,30),6.1552*pow(10,-74)*pow(10,30)};

  




   
   //  graph 1

    TGraph *gr1 = new TGraph(n,x1,y1);
    gr1->SetName("gr1");
    gr1->SetTitle("Flux vs K.E");

   // gr1->SetTitle("Dz_pt criteria");

    gr1->SetMarkerStyle(21);
    gr1->SetDrawOption("AP");
    gr1->SetLineColor(2);
    gr1->SetLineWidth(2);
    gr1->SetFillStyle(0);
    gr1->GetXaxis()->SetTitle("KE(MeV/nucl)");
    gr1->GetYaxis()->SetTitle("Flux");
   // gr1->GetYaxis()->SetRangeUser(pow(10,-80)*7,pow(10,-65));
   //  gr1->GetXaxis()->SetRangeUser(0,100000); 
  
   //gr1->Draw("AB");
    gr1->Draw("AC*");
// TCanvas::Update() draws the frame, after which one can change it
   c1->Update();
   //c1->GetFrame()->SetFillColor(21);
   //c1->GetFrame()->SetBorderSize(12);
   c1->Modified();
//c1->BuildLegend();
}