以前のリビジョンの文書です
Warning: "continue" targeting switch is equivalent to "break". Did you mean to use "continue 2"? in /mnt/hep_web/hep_web/member/n-kota/dokuwiki/inc/parser/handler.php on line 1458
Warning: Declaration of syntax_plugin_note::handle($match, $state, $pos, &$handler) should be compatible with DokuWiki_Syntax_Plugin::handle($match, $state, $pos, Doku_Handler $handler) in /mnt/hep_web/hep_web/member/n-kota/dokuwiki/lib/plugins/note/syntax.php on line 79
Warning: Declaration of syntax_plugin_note::render($mode, &$renderer, $indata) should be compatible with DokuWiki_Syntax_Plugin::render($format, Doku_Renderer $renderer, $data) in /mnt/hep_web/hep_web/member/n-kota/dokuwiki/lib/plugins/note/syntax.php on line 101
Warning: preg_match(): Compilation failed: invalid range in character class at offset 3120 in /mnt/hep_web/hep_web/member/n-kota/dokuwiki/inc/parser/lexer.php on line 118
====== ヒストグラムを規格化(Scale)して、さらにエラーバーをつけたい ====== ---- ===== ヒストグラムを規格化する ===== エントリー数の違うヒストグラムの形を比べる時などに活躍。\\ TH1::Scale を使って <code c> h->Scale(ScaleFactor); </code> のようにする。\\ 形の比較(h1,h2)の時などは <code c> h2->Scale(h1->GetEntries() / h2->GetEntries()); </code> のようにする。 ===== ヒストグラムにエラーバーをつける ===== ヒストグラムにエラーバーをつけるには、TH1::Draw の"E"オプションを使えば良い。\\ 各Binのエントリー数の平方根がエラーとしてつくようだ。\\ h->GetBinContent(BinNO), h->GetBinError(BinNO)で各Binの情報を取り出して、比べてみると良い。\\ "E"にもいくつかあって、 * "E": 誤差棒( error bars)を描画する。 * "E0": データがないビンを含め、誤差棒を描画する。 * "E1": 誤差棒の端に(誤差棒に対して)垂直な線をつけて描画する。 * "E2": 誤差棒に長方形をつけて描画する。 * "E3": 縦の誤差棒の端点を通るフィル領域を描画する。 * "E4": 誤差棒の端点を通る滑らかなフィル領域を描画する。 とある。"E1"とかを使っておけばOK。小文字でも良い。\\ ただ、"E"オプションを使うと元のヒストグラムの線が表示されないのが難点（一度やってみるとわかる）。\\ なので面倒だが <code c> h->Draw(); h->Draw("same e1"); </code> のようにする。 (ヒストグラムをコピーして片方をエラー用にして、分けた状態で重ね描いても良い)\\ ===== ヒストグラムを規格化(Scale)して、さらにエラーバーをつけたい ===== 本題だが、上の二つをそのまま組み合わせただけでは、 規格化されたヒストグラムのエラーがスケール後に再計算されてしまう。\\ 例えば10倍にスケールした時、エラーバーも10倍になっていてほしいが、10倍の統計があると思って小さく(√10倍)描かれてしまう。\\ これを避けるために <code c> h->Sumw2(); h->Scale(ScaleFactor); </code> のようにスケールする前にh%%->%%Sumw2()をはさまなくてはならない。 <note tip> TH1::Sumw2 について "sum of squares of weights"のこと。\\ ヒストグラムにWeight付きでFillするとき、Fillの前にh->Sumw2()しておくことで、各Binの √(Σw_i^2)が保存される。(これは、各Binのエラーになる。Weightが全て１なら√N) また、Sumw2()をcallすることでエラーが保存されるため、 ヒストグラムを書いた後にそのヒストグラムに対して演算（四則演算、scaleなども含む）をする場合などは演算前にSumw2()をcallしておかねばならない。 (そうしないと後からエラーを計算させた時に、演算後のヒストグラムの各Binに対して単純に各Binのエントリーの平方根を返すことになってしまう。) </note> ===== サンプルコード ===== <code c_mac> // マクロ test.C のサンプルコード void test() { gStyle->SetOptStat(0); TH1D *h1 = new TH1D("h1","h1",50,-5,5); TH1D *h2 = new TH1D("h2","h2",50,-5,5); TCanvas *c = new TCanvas("c","c",800,600); h1->FillRandom("gaus",10000); h2->FillRandom("gaus",1000); TH1D *h1e = (TH1D*)h1->Clone("h1e"); TH1D *h2e = (TH1D*)h2->Clone("h2e"); h2e->Sumw2(); h2->Scale(10); h2e->Scale(10); h1->SetLineColor(2); h2->SetLineColor(4); h1e->SetLineColor(2); h2e->SetLineColor(4); h1->Draw(); h1e->Draw("same e1"); h2->Draw("same"); h2e->Draw("same e1"); } {{:ja:root:scaleerror.pdf|}} </code>
