「絶望」の編集履歴(バックアップ)一覧はこちら
絶望 - (2007/12/03 (月) 19:34:46) の1つ前との変更点
追加された行は緑色になります。
削除された行は赤色になります。
今現在(12/2、4:30)アップされてるzetubou.Cの情報です。
バイナリから直で読む形式になってます。
例えば、gfitE(範囲min,範囲max,const,mean,sigma)
と打つことでフィットができます。変数を足せばダブルガウシアン+直線とかも可。
gfitE(範囲min,範囲max,const1,mean1,sigma1,const2,mean2,sigma2,slope,intercept)
ここでフィットに定義したガウシアンはconst×normalizeされたガウシアンという形になってます。rootのpredefinedされたガウシアンとは違うので注意。
ただソースがだいぶややこしくなったので、万が一使いたい場合は一旦ソースをチェックされてから使うことをお薦めします。
今現在(12/3、19:30)アップされてるzetubou.Cの情報です。
バイナリから直で読む形式になってます。
"dat2e("****.dat" ,"title",binの数(区切りの数),範囲min,範囲max)"
でEカウンター側のヒスト。
"dat2de("****.dat" ,"title",binの数(区切りの数),範囲min,範囲max)"
でdEカウンター側のヒスト。
カットをかけたい場合。
dat2e("filename", "title" ,int nbinx,int xmin,int xmax,int dEmin,int dEmax)でカットをかけれます。
ただの散布図+プロジェクションx,y
dat2scat("filename", "title")
散布図+プロジェクション+範囲限定散布図
dat2scat("filename","title",Int_t nbinx,Int_t xmin, Int_t xmax,Int_t nbiny,Int_t ymin, Int_t ymax)
てなカンジで。
例えば、gfitE(範囲min,範囲max,const,mean,sigma)
と打つことでシングルガウシアンでフィットができます。const,mean,sigmaはパラメーターの初期値。経験上meanの値をちゃんと定めればあとは割と適当で大丈夫。backgroundの上にピークがあるぞ、という場合は、
gfitE(範囲min,範囲max,const,mean,sigma、slope,intercept)
と2つ変数を足すことでガウシアン+線形でfit可。こいつでbackgroundを消したことにします。
変数を足せばダブルガウシアン+直線とかも可。
gfitE(範囲min,範囲max,const1,mean1,sigma1,const2,mean2,sigma2,slope,intercept)
ここでフィットに定義したガウシアンはconst×normalizeされたガウシアンという形になってます。rootのpredefinedされたガウシアンとは違うので注意。
χ^2/NDFがコンソール上に現れるのでこいつが1に近かったらfit成功。2,3…とかなら範囲や初期値を変えてみてください。どうにもよくならない場合もあるけど。
ただソースがだいぶややこしくなったので、万が一使いたい場合は一旦ソースをチェックされてから使うことをお薦めします。