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