ROOT/クレブシュ-ゴルダン係数の計算

クレブシュ-ゴルダン係数の計算†

参考サイト†

• クレブシュ-ゴルダン係数 - KobaWiki

• ROOT/MathMore - KobaWiki

• クレブシュ・ゴルダン計算機 - Wolfram|Alpha

• ROOT; 特殊関数

コード†

• ROOTのMathMoreというライブラリを使う。最小インストールでは入ってないらしいので注意。sahoには入ってる。

• ClebschGordan.h

  #include "TROOT.h"
  #include "Math/SpecFunc.h"
  
  double ClebschGordan(float j1,float m1,float j2,float m2,float j3,float m3){
    int two_ja=round(2*j1);
    int two_jb=round(2*j2);
    int two_jc=round(2*j3);
    int two_ma=round(2*m1);
    int two_mb=round(2*m2);
    int two_mc=round(2*m3);
  
    return pow(-1,(two_ja-two_jb+two_mc)/2)*sqrt(two_jc+1)
      *ROOT::Math::wigner_3j(two_ja,two_jb,two_jc,two_ma,two_mb,-two_mc);
  }

• calc.cxx

  #include <stdio.h>
  #include "ClebschGordan.h"
  
  int main(){
    printf("%f\n",ClebschGordan(2,2,2,-2,1,0));
    return 0;
  }

• コンパイル

  g++ calc.cxx -o calc `root-config --cflags --libs` -lMathMore

• 実行結果

  $ ./calc 
  0.632456

• 念の為、クレブシュ・ゴルダン計算機 で結果を確認してみたが、合ってた。

Topへ戻る