TGraphを大量に描く

ROOT 研究

TGraphをfor文使って大量に描きたいときは、二次元配列を使うとうまいことできる x[3][5] = {{…..}, {…..}, {…..}}; y[3][5] = {{…..}, {…..}, {…..}}; TGraph *g[3]; for(int i=0;i<3;i++){ g[i] = new TGraph(5, x[i], y[i]); g[i]->Draw(); }

2014-06-21

2次元ヒストグラムのピークのところでの軸の値をそれぞれもってきたいとき

研究 ROOT

二次元ヒストグラムでなんかの分布を見て、一番データが集中してるところでのx, yの値がほしいときがよくあるヒストグラムの名前をhist1とすると、 int maxbin; int maximum_xbin; int maximum_ybin; int dammyz; double maximum_x; double maximum_y; maxb…

2014-02-18

ROOT Treeからヒストグラム

研究 ROOT

t1->Draw("a[0]:a[1]"); から、 t1->Draw("a[0]:a[1]","","COLZ"); ってやると２次元ヒストがヒートマップみたいになる

2014-02-06

ROOTでコンパイル

研究 ROOT

ROOTでコンパイルするためのMakefile TARGETS=Hello ROOTFLAGS = $(shell root-config --cflags) ROOTLIBS = $(shell root-config --libs) CXXFLAGS = -Wall -O2 $(ROOTFLAGS) CXXLIBS = $(ROOTLIBS) all: $(TARGETS) Hello: Hello.o g++ -o $@ Hello.o $(CX…

2014-01-23

gslでできること一覧

gsl 物理学研究

gslでできること一覧(マニュアルより抜粋) 複素数多項式の求根法特殊関数置換組み合わせ整列線形代数 CBLAS ライブラリ高速フーリエ変換乱数数値積分乱数分布ヒストグラム統計モンテカルロ積分微分方程式シミュレーティド・アニーリングベク…

2014-01-23

gslの使い方

物理学研究 gsl

gslとは、GNU Scientific Library。いろいろめんどくさい計算すぐできる。インストールしただけではただのライブラリ。例：ベッセル関数 $J_0(x)$ の、$x=5$ での値を計算して出力する #include <stdio.h> #include <gsl/gsl_sf_bessel.h> int main (void){ double x = 5.0; double y =</gsl/gsl_sf_bessel.h></stdio.h>…

2014-01-22

Geant4でのデータ収集の仕方(2)

研究 Geant4 物理学

一番簡単なやり方、ユーザフック exampleB1などで使われる。B1では線量計算をStep、Event、Run actionを自作して求める。 1. 構造物を作り、stepping actionをするvolumeを指定 2. Stepping actionクラス内で、今のStepがスコアすべき物体中にあるかどうか…