{

   int N=225;

   double lamdaalpha[N];
   double lamdabeta[N];
   double khi2[N],k[N], I[N],J[N];

   ifstream f("val.txt");
 
   TH2D *h = new TH2D("h","h",15,0.05,1.55,15,0.05,1.55);


   for(int i=0; i<N; i++){
      f>>k[i]>>I[i]>>J[i]>>lamdaalpha[i]>>lamdabeta[i]>>khi2[i];
      h->Fill(lamdaalpha[i],lamdabeta[i],khi2[i]);//remplissage de hitogramme
   }

   TCanvas *c1=new TCanvas("c1","c1",800,600);
   c1->SetGridx();
   c1->SetGridy();

   Double_t contours[3];
   ///contours[0] = 0.5;
   ///contours[1] = 1.;
   ///contours[2] = 2.;
   contours[0] = 60.;
   contours[1] = 100.;
   contours[2] = 120.;

   h->GetXaxis()->SetTitle("#lambda_{#alpha}");
   h->GetYaxis()->SetTitle("#lambda_{#beta}");
   h->SetContour(3, contours);

   h->Draw("CONT4"); 
}