ROOTのTSpectrumクラスを使ってpeak searchをする

PeakSearchByTSpectrum.md

ROOTのTSpectrumクラスを使ってpeak searchをする

TSpectrum を使うと簡単にpeak searchができた。

• How to Find Peaks in Histograms ? | ROOT

の中のサンプルプログラムを動かしたら大体分かる。

TSpectrumインスタンスの生成

ピークの最大数を適当に決めてTSpecrumのインスタンスをつくる。

Int_t maxpeaks = 10;
TSpectrum *s = new TSpectrum(maxpeaks);

ちなみにこれだとresolutionは1になる。今はピークが探せさえすればいいので気にしない。あとからTSpectum::SetResolutionで変えれる。resolution=1はピーク間の距離は3 sigmaに対応するのだとか。詳しくはTSpectrum::SetResolutionを参照。

ピークを探す

TSpectrum::Searchを使う。h1というTH1Fオブジェクトがあったら、

s->Search(h1, 1, "new");

をするだけでいい。めっちゃ簡単。2つ目の引数はsigma。詳しくはマニュアルって書いてあるのであとで読む。3つ目はオプション。4つ目は省略してるけどthresholdを設定する。ピークサーチが終わるとfNPeaks, fPositionXなどの変数に情報が入ってるので、TSpectrum::GetNPeaks, TSpectrum::GetPositionXとかで取ってくる。

Float_t *xpeaks = s->GetPositionX();

あとはresolutionやthresholdをいじってoptimizeしたり、それはピークじゃないってやつをfitできるかどうかで除けばよいのかな。あとsmoothingしたりも。

コンパイル

コンパイルするときのオプションに、-lSpectrumが必要。これは

root-config --libs

しても出てこないので自分でMakefileに書いておく必要がある(?)。

• RootTalk • View topic - TSpectrum error on the version installed on my laptop

に書いてあったけどもっと分かりやすいところに書いて欲しい。