{
  TString outputfilename;
  outputfilename = "Spatial.txt";
  ofstream outputfile;
  outputfile.open(outputfilename.Data(),ios::app);

//=========Macro generated from canvas: c1/c1
//=========  (Wed Apr 15 11:54:19 2009) by ROOT version5.20/00
/*   TCanvas *c1 = new TCanvas("c1", "c1",9,44,700,500);
   c1->Range(-0.9858439,-1.117922,0.9858439,1.117922);
   TView *view = TView::CreateView(1);
   view->SetRange(-160,-180,-40,100,20,80);
   c1->SetBorderSize(2);
   c1->SetFrameFillColor(0);
*/
 
   TH1 *Spatial = new TH3F("Spatial","Distribution spatiale de puissance",261,-160,100,201,-180,20,121,-40,80);
   Spatial->SetBinContent(47236,8.651172);

   Spatial->SetEntries(1.172108e+07);
   

   Double_t binc(0.);
   Double_t binc2(0.);
   Int_t binno(0.), binno2(0.);
   for(int x=1; x < 262; x++)
     {
       for(int y=1; y < 202; y++)
	 {
	   for(int z=1; z < 192; z++) 
	     {
	       binno = Spatial->GetBin(x,y,z);
	       //cout << binno << endl;
	       //binc  = (Double_t)Spatial -> GetBinContent(x,y,z);
	       binc2 = (Double_t)Spatial -> GetBinContent(binno);

	       if(binno == 47236) cout << "The Binc : " << binc2 << endl;
	       //cout << binc << endl;
	       if(binc =! 0. && binno2==0)
		 {
		   cout << (x) << "\t" << (y) << "\t" << (z) << "\t" << binc << endl;
		   cout << (x-160) << "\t" << (y-180) << "\t" << (z-40) << "\t" << binc << endl;
		   //outputfile << (x-160) << "\t" << (y-180) << "\t" << (z-40) << "\t" << binc << endl;
		   binno2++;
		 }
	       //binno++;
	     }
	 }
     }
   outputfile.close();
}