n3956.net

よく使う物理定数
名称	記号	数値	単位
真空中の光速	$$c$$	2.99792458E+08	$$\mathrm{m}\ \mathrm{s}^{-1}$$
真空の透磁率	$$\mu_{0}$$	1.25664E-06	$$\mathrm{N}\ \mathrm{A}^{-2}$$
真空の誘電率	$$\varepsilon_{0}$$	8.85419E-12	$$\mathrm{F}\ \mathrm{m}^{-1}$$
万有引力定数	$$G$$	6.67428E-11	$$\mathrm{N}\ \mathrm{m}^{2}\ \mathrm{kg}^{-2}$$
Planck 定数	$$h$$	6.62607E-34	$$\mathrm{J}\ \mathrm{s}$$
		4.13567E-15	$$\mathrm{eV}\ \mathrm{s}$$
Dirac 定数	$$\hbar=\frac{h}{2\pi}$$	1.05457E-34	$$\mathrm{J}\ \mathrm{s}$$
		6.58212E-16	$$\mathrm{eV}\ \mathrm{s}$$
素電荷	$$e$$	1.60218E-19	$$\mathrm{C}$$
ボーア磁子	$$\mu_{\mathrm{B}}=\frac{e\hbar}{2m_{\mathrm{e}}}$$	9.27401E-24	$$\mathrm{J}\ \mathrm{T}^{-1}$$
核磁子	$$\mu_{\mathrm{N}}=\frac{e\hbar}{2m_{\mathrm{p}}}$$	5.05078E-27	$$\mathrm{J}\ \mathrm{T}^{-1}$$
電子の質量	$$m_{\mathrm{e}}$$	9.10938E-31	$$\mathrm{kg}$$
	$$m_{\mathrm{e}}c^{2}$$	0.5109989	$$\mathrm{MeV}$$
陽子の質量	$$m_{\mathrm{p}}$$	1.67262E-27	$$\mathrm{kg}$$
	$$m_{\mathrm{p}}c^{2}$$	938.272	$$\mathrm{MeV}$$
中性子の質量	$$m_{\mathrm{n}}$$	1.67493E-27	$$\mathrm{kg}$$
	$$m_{\mathrm{n}}c^{2}$$	939.565	$$\mathrm{MeV}$$
電子の磁気モーメント	$$\mu_{\mathrm{e}}$$	-9.28476E-24	$$\mathrm{J}\ \mathrm{T}^{-1}$$
陽子の磁気モーメント	$$\mu_{\mathrm{p}}$$	1.41061E-26	$$\mathrm{J}\ \mathrm{T}^{-1}$$
中性子の磁気モーメント	$$\mu_{\mathrm{n}}$$	-9.66236E-27	$$\mathrm{J}\ \mathrm{T}^{-1}$$
中性子の磁気回転比	$$\gamma_{\mathrm{n}}=\frac{g_{\mathrm{n}}\mu_{\mathrm{N}}}{\hbar}$$	1.83247E+8	$$\mathrm{T}^{-1}\ \mathrm{s}^{-1}$$ $$\mathrm{Hz}\ \mathrm{T}^{-1}$$
		1.83247E+4	$$\mathrm{G}^{-1}\ \mathrm{s}^{-1}$$ $$\mathrm{Hz}\ \mathrm{G}^{-1}$$
統一原子質量単位	$$[\mathrm{u}]$$	1.66054E-27	$$\mathrm{kg}$$
Avogadro 定数	$$N_{\mathrm{A}}$$	6.02214E+23	$$\mathrm{mol}^{-1}$$
Boltzmann 定数	$$k_{\mathrm{B}}$$	1.38065E-23	$$\mathrm{J}\ \mathrm{K}^{-1}$$

emcas で認識されるパスを通常の端末と同じにする

投稿日: 2013-10-6 作成者: n3956

問題：
例えば emacs で tex を文章をコンパイルするとき（例えば　yatex などの機能で C-c C-t j する），No such directory or file. bad interpreter と出る．

（emacs からではなく）通常のOSの機能として開いた端末からでは latex にパスを通していることは確認（例えば ~/.bashrc に
PATH=$PATH:/usr/local/texlive/??/bin
と書いている）．

emacs で呼ばれる bash が .bashrc の設定を読んでいないことが原因なので、emacs に通常使っている PATH を教える必要がある．
それには以下のような設定を .emacs に書き加えればよい。
参考　http://stackoverflow.com/a/6415812

(defun set-exec-path-from-shell-PATH ()
  (let((path-from-shell (replace-regexp-in-string
                         "[ \t\n]*$"
                         ""
                         (shell-command-to-string "$SHELL --login -i -c 'echo $PATH'"))))
    (setenv "PATH" path-from-shell)
    (setq eshell-path-env path-from-shell) ; for eshell users
    (setq exec-path (split-string path-from-shell path-separator))))

(defun set-exec-path-from-shell-PATH ()

(let((path-from-shell (replace-regexp-in-string

"[ \t\n]*$"

(shell-command-to-string "$SHELL --login -i -c 'echo $PATH'"))))

(setenv "PATH" path-from-shell)

(setq eshell-path-env path-from-shell) ; for eshell users

(setq exec-path (split-string path-from-shell path-separator))))

余談ですが，シェルスクリプトの 1 行目に書いてある
#!/bin/sh
などのインタプリタ指定のための行はシェバン (shebang) というらしい．

カテゴリー: 未分類 | コメントする

ROOT のヒストグラムの内容を python のリストとして得る

投稿日: 2013-9-17 作成者: n3956

基本的に ROOT のヒストグラムならROOTの機能 (Draw) で描画して絵を得るのが素直ですが，python の方でいろいろと処理したい場合に便利？

ビンの下端の値をリストとして得る
[testh.GetBinLowEdge(i+1) for i in range(testh.GetNbinsX())]
ビンの中身をリストとして得る
[testh.GetBinContent(i+1) for i in range(testh.GetNbinsX())]
使用例

>>> testh = TH1F("testh", "testh title", 20, -1.0, 1.0)
>>> testh.FillRandom("gaus", 100000) # mean=0, sigma=1 のガウス分布を 100000 サンプル

>>> [testh.GetBinLowEdge(i+1) for i in range(testh.GetNbinsX())] 
[-1.0, -0.9, -0.8, -0.7, -0.6,
 -0.5, -0.3999999999999999, -0.29999999999999993, -0.19999999999999996,
 -0.09999999999999998, 0.0, 0.10000000000000009, 0.20000000000000018,
 0.30000000000000004, 0.40000000000000013, 0.5, 0.6000000000000001,
 0.7000000000000002, 0.8, 0.9000000000000001]
>>> [testh.GetBinContent(i+1) for i in range(testh.GetNbinsX())] 
[0.0,
 1.0,
 20.0,
 117.0,
 444.0,
 1670.0,
 4427.0,
 9067.0,
 15108.0,
 18890.0,
 19146.0,
 15026.0,
 9332.0,
 4509.0,
 1651.0,
 458.0,
 111.0,
 22.0,
 1.0,
 0.0]

>>> testh = TH1F("testh", "testh title", 20, -1.0, 1.0)

>>> testh.FillRandom("gaus", 100000) # mean=0, sigma=1 のガウス分布を 100000 サンプル

>>> [testh.GetBinLowEdge(i+1) for i in range(testh.GetNbinsX())]

[-1.0, -0.9, -0.8, -0.7, -0.6,

-0.5, -0.3999999999999999, -0.29999999999999993, -0.19999999999999996,

-0.09999999999999998, 0.0, 0.10000000000000009, 0.20000000000000018,

0.30000000000000004, 0.40000000000000013, 0.5, 0.6000000000000001,

0.7000000000000002, 0.8, 0.9000000000000001]

>>> [testh.GetBinContent(i+1) for i in range(testh.GetNbinsX())]

[0.0,

1.0,

20.0,

117.0,

444.0,

1670.0,

4427.0,

9067.0,

15108.0,

18890.0,

19146.0,

15026.0,

9332.0,

4509.0,

1651.0,

458.0,

111.0,

22.0,

1.0,

0.0]

ROOT の Draw と matplotlib でプロットしてみる

>>> xp = [testh.GetBinCenter(i+1) for i in range(testh.GetNbinsX())]  # ビンの中心値のリスト
>>> yp = [testh.GetBinContent(i+1) for i in range(testh.GetNbinsX())] # ビンの値のリスト
>>> plot(xp, yp, 'o-') # matplotlib でプロット
>>> testh.Draw()       # pyROOT で Draw

>>> xp = [testh.GetBinCenter(i+1) for i in range(testh.GetNbinsX())] # ビンの中心値のリスト

>>> yp = [testh.GetBinContent(i+1) for i in range(testh.GetNbinsX())] # ビンの値のリスト

>>> plot(xp, yp, 'o-') # matplotlib でプロット

>>> testh.Draw() # pyROOT で Draw

カテゴリー: matplotlib, pyROOT, python, ROOT | コメントする

pyROOT tips

投稿日: 2013-9-16 作成者: n3956

Legend の作成

from ROOT import TLegend
def setLegend(leg, entries):
    if isinstance(entries, (list, tuple)):
        for ent in entries:
            leg.AddEntry(ent, ent.GetName())

from ROOT import TLegend

def setLegend(leg, entries):

if isinstance(entries, (list, tuple)):

for ent in entries:

leg.AddEntry(ent, ent.GetName())

使い方

leg = TLegend(0.65, 0.65, 0.85, 0.85)
setLegend(leg, (h1, h2, h3))

1 2	leg = TLegend(0.65, 0.65, 0.85, 0.85) setLegend(leg, (h1, h2, h3))

複数のカットを and で連結する

from ROOT import TCut
def andCut(cuts):
    if isinstance(cuts, (list, tuple)):
        cutnames = [var.GetTitle() for var in cuts]
        resCut = TCut( "&&".join(cutnames) )
    else:
        resCut = cuts
    return resCut

from ROOT import TCut

def andCut(cuts):

if isinstance(cuts, (list, tuple)):

cutnames = [var.GetTitle() for var in cuts]

resCut = TCut( "&&".join(cutnames) )

else:

resCut = cuts

return resCut

and(&&)で結んだカットを定義

cut1  = TCut("-1<x")
cut2  = TCut("x<1")
cut12 = TCut( "(%s)&&(%s)" %(cut1.GetTitle(), cut2.GetTitle()) )

cut1 = TCut("-1<x")

cut2 = TCut("x<1")

cut12 = TCut( "(%s)&&(%s)" %(cut1.GetTitle(), cut2.GetTitle()) )

ROOT の TCut オブジェクトからカットを表現する文字列を得るときは、GetTitle() メソッドを使えば良い。

カテゴリー: pyROOT, python, ROOT | コメントする

pyOpenGL を Windows8 (64bit) にインストールする

投稿日: 2013-8-17 作成者: n3956

環境は OS：Windows8 Pro 64 bit，Python: Python2.7 (64bit)

> pip install "pyopengl"
pip が使えなければ，PyPI の pyOpenGL のページ
https://pypi.python.org/pypi/PyOpenGL
から PyOpenGL-3.0.2.win-amd64.exe をダウンロードし，インストールウィザードに従う．もしくは，ソース PyOpenGL-3.0.2.zip をダウンロードして解凍した後，当該フォルダでコマンドプロンプトから
> python setup.py install
とする．これで python から import OpenGL できるようになる．
しかし，このままでは適当なサンプルスクリプト（後述）を動かそうとすると以下のエラーが出た．

 Traceback (most recent call last):
 File "C:\*\***.py", line xx, in
 setup()
 File "C:\*\***.py", line xx, in setup
 glutInit(sys.argv)
 File "C:\*\lib\site-packages\OpenGL\GLUT\special.py", line 324, in glutInit
 _base_glutInit( ctypes.byref(count), holder )
 TypeError: 'NoneType' object is not callable

Traceback (most recent call last):

File "C:\*\***.py", line xx, in

setup()

File "C:\*\***.py", line xx, in setup

glutInit(sys.argv)

File "C:\*\lib\site-packages\OpenGL\GLUT\special.py", line 324, in glutInit

_base_glutInit( ctypes.byref(count), holder )

TypeError: 'NoneType' object is not callable

実行時に引数は指定していないので，sys.argv の中身はこのスクリプト teapot.py のパスのはず．どうも DLL が読み込まれていないことが原因らしい．解凍したソースの C:\***\Downloads\PyOpenGL-3.0.2\OpenGL\DLLS にある glut32.dll と glut64.dll を C:\Windows\System32 フォルダへコピーしたらこのエラーはでなくなった．

テストに使ったサンプルスクリプトは以下（こちら ”Pythonでゲーム作りますが何か？” を参考にした）．

#!/usr/bin/env python
#coding:utf-8

from OpenGL.GL import *
from OpenGL.GLU import *
from OpenGL.GLUT import *
import sys

def main():
    glutInit(sys.argv)
    glutInitDisplayMode(GLUT_RGB | GLUT_SINGLE | GLUT_DEPTH)
    glutInitWindowSize(300, 300)     # window size
    glutInitWindowPosition(100, 100) # window position
    glutCreateWindow("teapot")       # show window
    glutDisplayFunc(display)         # draw callback function
    glutReshapeFunc(reshape)         # resize callback function
    init(300, 300)
    glutMainLoop()

def init(width, height):
    """ initialize """
    glClearColor(0.0, 0.0, 0.0, 1.0)
    glEnable(GL_DEPTH_TEST) # enable shading

    glMatrixMode(GL_PROJECTION)
    glLoadIdentity()
    ##set perspective 
    gluPerspective(45.0, float(width)/float(height), 0.1, 100.0) 

def display():
    """ display """
    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)
    glMatrixMode(GL_MODELVIEW)
    glLoadIdentity()
    ##set camera
    gluLookAt(0.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0)    
    ##draw a teapot
    glColor3f(1.0, 0.0, 0.0)
    glutWireTeapot(1.0)   # wireframe
#    glutSolidTeapot(1.0)  # solid
    glFlush()  # enforce OpenGL command

def reshape(width, height):
    """callback function resize window"""
    glViewport(0, 0, width, height)
    glMatrixMode(GL_PROJECTION)
    glLoadIdentity()
    gluPerspective(45.0, float(width)/float(height), 0.1, 100.0)

if __name__ == "__main__":
    main()

#!/usr/bin/env python

#coding:utf-8

from OpenGL.GL import *

from OpenGL.GLU import *

from OpenGL.GLUT import *

import sys

def main():

glutInit(sys.argv)

glutInitDisplayMode(GLUT_RGB | GLUT_SINGLE | GLUT_DEPTH)

glutInitWindowSize(300, 300) # window size

glutInitWindowPosition(100, 100) # window position

glutCreateWindow("teapot") # show window

glutDisplayFunc(display) # draw callback function

glutReshapeFunc(reshape) # resize callback function

init(300, 300)

glutMainLoop()

def init(width, height):

""" initialize """

glClearColor(0.0, 0.0, 0.0, 1.0)

glEnable(GL_DEPTH_TEST) # enable shading

glMatrixMode(GL_PROJECTION)

glLoadIdentity()

##set perspective

gluPerspective(45.0, float(width)/float(height), 0.1, 100.0)

def display():

""" display """

glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)

glMatrixMode(GL_MODELVIEW)

glLoadIdentity()

##set camera

gluLookAt(0.0, 1.0, 5.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0)

##draw a teapot

glColor3f(1.0, 0.0, 0.0)

glutWireTeapot(1.0) # wireframe

# glutSolidTeapot(1.0) # solid

glFlush() # enforce OpenGL command

def reshape(width, height):

"""callback function resize window"""

glViewport(0, 0, width, height)

glMatrixMode(GL_PROJECTION)

glLoadIdentity()

gluPerspective(45.0, float(width)/float(height), 0.1, 100.0)

if __name__ == "__main__":

main()

> python teapot.py

http://sourceforge.net/projects/pyopengl/files/
にある PyOpenGL-Demo をインストールすると色々なデモが見られます．例えば
> python .\PyOpenGL-Demo\GLUT\gears.py

25828 frames in 5.0 seconds = 5165.600 FPS
26693 frames in 5.0 seconds = 5338.600 FPS
26626 frames in 5.0 seconds = 5325.200 FPS
26851 frames in 5.0 seconds = 5370.200 FPS
26447 frames in 5.0 seconds = 5289.400 FPS
...

25828 frames in 5.0 seconds = 5165.600 FPS

26693 frames in 5.0 seconds = 5338.600 FPS

26626 frames in 5.0 seconds = 5325.200 FPS

26851 frames in 5.0 seconds = 5370.200 FPS

26447 frames in 5.0 seconds = 5289.400 FPS

...

カテゴリー: pyOpenGL, python | コメントする

Ratings
- リストの操作 (0.00)
- 三角関数に関する公式 (0.00)
- WSLでGUIがサポートされるようになった（Windows10でも） (0.00)
StatPress TopPosts
- Page: Home
- /blog/
- matplotlib のプロットの単純な例
- matplotlib で等高線グラフ
- pyOpenGL を Windows8 (64bit) にインストールする
Calendar
2025年7月

月火水木金土日

« 5月

1 2 3 4 5 6

7 8 9 10 11 12 13

14 15 16 17 18 19 20

21 22 23 24 25 26 27

28 29 30 31

よく使う物理定数

emcas で認識されるパスを通常の端末と同じにする

ROOT のヒストグラムの内容を python のリストとして得る

pyROOT tips

Legend の作成

複数のカットを and で連結する

and(&&)で結んだカットを定義

pyOpenGL を Windows8 (64bit) にインストールする

最近の投稿

最近のコメント

Archive

カテゴリー

メタ情報

Ratings

StatPress TopPosts

Calendar

月	火	水	木	金	土	日
« 5月
	1	2	3	4	5	6
7	8	9	10	11	12	13
14	15	16	17	18	19	20
21	22	23	24	25	26	27
28	29	30	31