void plotgraph(){
  gStyle->SetPaintTextFormat("5.1f");
  std::ifstream is("data.dat");
  TGraph2D* gr = new TGraph2D();
  //  gr->SetMaxIter(500000);
  int point = 0;
  double cls;
  while(is >> cls){
    double cls; int np; double mass1; double mass2; 
    double g1; double g2;
    is >> mass1;
    is >> mass2;
    is >> g1;
    is >> g2;
    is >> np;    
    if(np <= 1000){
      gr->SetPoint(point,g1,np,cls);
      point++;
    }
  }


  TCanvas* c2 = new TCanvas("c2","contour",600,500);
  c2->cd();
  TH2D* gg = gr->GetHistogram();
  gg->SetContour(1);
  gg->SetContourLevel(0, 90.);
  gg->Draw("CONT3");

}