4

rootpy と matplotlib を使用して ROOT 2D ヒストグラムをプロットしようとしています。

これに使用するコードは次のとおりです。

from rootpy.io import File
from rootpy.plotting import Hist
import rootpy.plotting.root2matplotlib as rplt
import matplotlib.pyplot as plt
inputFile = File('mydata.root', 'read')
h_response = inputFile.myfolder.response

plt.figure(figsize=(16, 10), dpi=100)
rplt.hist(h_response, label='response matrix')
h_response.Draw()
plt.xlabel('reconstructed $E_{\mathrm{T}}^{miss}$')
plt.ylabel('Generated $E_{\mathrm{T}}^{miss}$')
plt.title('Response Matrix')
plt.savefig('ResponseMatrix.png')

ただし、これによりエラーメッセージが表示されます。

Traceback (most recent call last):
  File "/storage/Dropbox/Workspace/Analysis/DailyPythonScripts/src/unfolding.py", line 66, in <module>
    rplt.hist(h_response, label='response matrix')
  File "/usr/local/lib/python2.7/dist-packages/rootpy-0.7.0_a0-py2.7-linux-x86_64.egg/rootpy/plotting/root2matplotlib.py", line 140, in hist
    snap_zero=snap_zero)
  File "/usr/local/lib/python2.7/dist-packages/rootpy-0.7.0_a0-py2.7-linux-x86_64.egg/rootpy/plotting/root2matplotlib.py", line 82, in _set_bounds
    ywidth = ymax - ymin
TypeError: unsupported operand type(s) for -: 'list' and 'list'

明らかに、私は間違った rootpy2matplotlib モジュールを使用しているので、調べてみました: モジュールが提供するもの: hist、bar、および errorbar 関数 - >= 2D に固有のものではありません。

何か不足していますか?簡単な回避策はありますか?

PS: この質問に「rootpy」タグを付けたいのですが、できません。この質問は非常に具体的であるため、申し訳ありません。

4

2 に答える 2

4

rootpy の root2matplotlib インターフェイスは、2D ROOT ヒストグラムをプロットするための hist2d、imshow、contour 関数を提供するようになりました。次の例を参照してください。

https://github.com/rootpy/rootpy/blob/master/examples/plotting/plot_matplotlib_hist2d.py

from matplotlib import pyplot as plt
from rootpy.plotting import root2matplotlib as rplt
from rootpy.plotting import Hist2D
import numpy as np

a = Hist2D(100, -3, 3, 100, 0, 6)
a.fill_array(np.random.multivariate_normal(
    mean=(0, 3),
    cov=np.arange(4).reshape(2, 2),
    size=(1E6,)))

fig, (ax1, ax2, ax3) = plt.subplots(nrows=1, ncols=3, figsize=(15, 5))

ax1.set_title('hist2d')
rplt.hist2d(a, axes=ax1)

ax2.set_title('imshow')
im = rplt.imshow(a, axes=ax2)

ax3.set_title('contour')
rplt.contour(a, axes=ax3)

fig.subplots_adjust(right=0.8)
cbar_ax = fig.add_axes([0.85, 0.15, 0.05, 0.7])
fig.colorbar(im, cax=cbar_ax)

plt.show()

ここに画像の説明を入力

于 2013-09-18T08:00:22.000 に答える
3

TH2Dヒストグラムを含むROOTファイルでスクリプトを試しました。すべてがうまくいった。

/ opt / rootpy#cat version.txtは私に:0.7.0を与えます

私がチェックした場合

/usr/local/lib/python2.7/dist-packages/rootpy-dev-py2.7.egg/rootpy/ploting/root2matplotlib.py

表示されるエラーメッセージと比較すると、異なるバージョンのrootpyを使用しているように見えます。

rootpyの最新バージョンを試してください。

于 2013-03-04T14:19:42.650 に答える