{ 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(); }