src/Level1/Wigner.cc
// ______________________________________________________________________
//          REDUCED WIGNER ROTATION MATRIX ELEMENTS: RANK 1/2 
// ______________________________________________________________________
REDUCED WIGNER ROTATION MATRIX ELEMENTS: RANK 1/2
REDUCED WIGNER ROTATION MATRIX ELEMENTS : RANK 1 / 2
double d1half(int n, double beta )
  {
  double r = 0;
  beta = beta*PI/360.;
    switch (n)
      {
      case -1:				
        r = -sin(beta);			
        break;				
      case 1:				
        r = cos(beta);			
        break;				
      default:
        Wigner_error(0);
        cout << "1/2," << n << "/2";
        Wigner_fatality(10);
        break;
      }
  return r;
  }
double d1half ( int n , double beta ) { double r = 0 ; beta = beta * PI / 360. ; switch ( n ) { case -1 : r = - sin ( beta ) ; break ; case 1 : r = cos ( beta ) ; break ; default : Wigner_error ( 0 ) ; cout << " " << n << " " ; Wigner_fatality ( 10 ) ; break ; } return r ; }
______________________________________________________________________
REDUCED WIGNER ROTATION MATRIX ELEMENTS: RANK 1/2
______________________________________________________________________
______________________________________________________________________ REDUCED WIGNER ROTATION MATRIX ELEMENTS : RANK 1 / 2 ______________________________________________________________________
// Input		n     : momentum index // 			beta  : Euler angle // Output		r     : rank 1/2, m = 1/2 reduced Wigner //			        rotation matrix element // Note 		      : Euler angle input in degrees //  (1/2) // d        (beta) = - sin(beta/2) //  1/2,-1/2 //  (1/2) // d       (beta) = cos(beta/2) //  1/2,1/2