#include "Math/SpecFunc.h"
#include <iostream>

Double_t Clebsch_Gordan(Int_t two_ja, Int_t two_jb, Int_t two_jc,
			Int_t two_ma, Int_t two_mb, Int_t two_mc ) {
  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);
}

void calc_Clebsch_Gordan() {
  Int_t two_j1 = 3;
  Int_t two_j2 = 5;
  Int_t two_j3 = 4;
  Int_t two_m1 = 1;
  Int_t two_m2 = -1;
  Int_t two_m3 = 0;
  std::cout << "C("
	    << two_j1/2. << ", " << two_j2/2. << ", " << two_j3/2. << "; "
	    << two_m1/2. << ", " << two_m2/2. << ", " << two_m3/2. << ") = "
	    << Clebsch_Gordan(two_j1,two_j2,two_j3,two_m1,two_m2,two_m3) << std::endl;
}