質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.46%
C++

C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

Q&A

解決済

1回答

1896閲覧

ROOTで2次元ヒストグラムの画素値をファイルに出力したい

instantuser

総合スコア11

C++

C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

0グッド

0クリップ

投稿2020/04/21 15:41

ROOTの2次元ヒストグラムの値をCSVファイルとして出力したいです。

以下のようなコードを書いて図が出力するところまで確認しているのですが、
k = h2->GetBinContent(i) の時点で値がすべて0となっています。

C++

1 2 3void NtuplePout(){ 4 TH2D *h2 = new TH2D("name","name",10,-300,300,10,0,3); 5 Coincidences->Draw("sinogramTheta:sinogramS>>h2","","colz"); 6 7 int x = h2->GetXaxis()->GetNbins()+2; 8 int y = h2->GetYaxis()->GetNbins()+2; 9 int k = 0; 10 11 ofstream outf("test.txt"); 12 13 14 for(int i=1; i<=x*y; i++){ 15 k = h2->GetBinContent(i); 16 if(k>0){ 17 printf("**********************"); 18 } 19 20 printf("%d",k); 21 outf<<k; 22 23 if(i%x==0){ 24 printf("\n"); 25 outf<<"\n"; 26 }else{ 27 printf(","); 28 outf<<","; 29 } 30 } 31 outf.close(); 32 return; 33}

何かおかしいところはございますか?

コマンド
root outout.root NtuplePout.C

結果
イメージ説明
test.txtも全て0

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

コンストラクタの指定方法が間違っているのでは。

以下のようにしてみては。

TH2D *h2 = new TH2D("name", "name", 10, -300.0, 300.0, 10, 0.0, 3.0);

投稿2020/04/21 19:49

TaroToyotomi

総合スコア1430

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

instantuser

2020/04/22 01:32

結果は変わらなかったです。。。 確かにh2に正しく格納されていないように思います(範囲が-300~300になっていない) 最初の部分を以下のように変更すると指定した範囲で表示されるのですが //TH2D *h2 = new TH2D("name","name",10,-300.0,300.0,10,0.0,3.0); Coincidences->Draw("sinogramTheta:sinogramS>>h2(128,-300.,300.,128,0.,3.2)"); 以下のようなエラーがでます /home/gate/GateContrib/imaging/myPET/NtuplePout.C:5:11: error: use of undeclared identifier 'h2' int x = h2->GetXaxis()->GetNbins(); ^ /home/gate/GateContrib/imaging/myPET/NtuplePout.C:6:11: error: use of undeclared identifier 'h2' int y = h2->GetYaxis()->GetNbins(); ^ /home/gate/GateContrib/imaging/myPET/NtuplePout.C:13:9: error: use of undeclared identifier 'h2' k = h2->GetBinContent(i);
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.46%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問