EvtdFunctionSingle Class Reference

#include <EvtdFunctionSingle.hh>

Public Member Functions
	EvtdFunctionSingle ()
	~EvtdFunctionSingle ()
void	init (int j, int m1, int m2)
double	d (int j, int m1, int m2, double theta)

Detailed Description

Definition at line 24 of file EvtdFunctionSingle.hh.

Constructor & Destructor Documentation

EvtdFunctionSingle()

EvtdFunctionSingle::EvtdFunctionSingle

(

)

Definition at line 29 of file EvtdFunctionSingle.cc.

                                       {
  _j    = 0;
  _m1   = 0;
  _m2   = 0;
  _coef = 0;
  _kmin = 0;
  _kmax = 0;
}

~EvtdFunctionSingle()

EvtdFunctionSingle::~EvtdFunctionSingle

(

)

Definition at line 38 of file EvtdFunctionSingle.cc.

                                        {
  if ( _coef != 0 ) delete[] _coef;
}

Member Function Documentation

d()

double EvtdFunctionSingle::d	(	int	j,
		int	m1,
		int	m2,
		double	theta )

Definition at line 80 of file EvtdFunctionSingle.cc.

                                                                  {
 
  // report(INFO,"EvtGen") << "j,m1,m2:"<<j<<","<<m1<<","<<m2<<" theta:"<<theta<<endl;
 
  assert( j == _j );
  assert( m1 == _m1 );
  assert( m2 == _m2 );
 
  double c2 = cos( 0.5 * theta );
  double s2 = sin( 0.5 * theta );
 
  double d = 0.0;
 
  int k;
  for ( k = _kmin; k <= _kmax; k += 2 )
  {
    d += _coef[( k - _kmin ) / 2] * pow( c2, ( 2 * _j - 2 * k + m2 - m1 ) / 2 ) *
         pow( s2, ( 2 * k - m2 + m1 ) / 2 );
  }
 
  return d;
}

Referenced by EvtdFunction::d(), and d().

init()

void EvtdFunctionSingle::init	(	int	j,
		int	m1,
		int	m2 )

Definition at line 42 of file EvtdFunctionSingle.cc.

                                                     {
 
  assert( abs( m2 ) >= abs( m1 ) );
  assert( m2 >= 0 );
 
  _j  = j;
  _m1 = m1;
  _m2 = m2;
 
  _kmin = _m2 - _m1;
  _kmax = _j - _m1;
 
  assert( _kmin <= _kmax );
 
  _coef = new double[( _kmax - _kmin ) / 2 + 1];
 
  int k;
 
  for ( k = _kmin; k <= _kmax; k += 2 )
  {
    int sign = 1;
    if ( ( k - _m2 + _m1 ) % 4 != 0 ) sign = -sign;
    // fkw change to satisfy KCC
    double fkwTmp = fact( ( _j + _m2 ) / 2 ) * fact( ( _j - _m2 ) / 2 ) *
                    fact( ( _j + _m1 ) / 2 ) * fact( ( _j - _m1 ) / 2 );
    _coef[( k - _kmin ) / 2] =
        sign * sqrt( fkwTmp ) /
        ( fact( ( _j + _m2 - k ) / 2 ) * fact( k / 2 ) * fact( ( _j - _m1 - k ) / 2 ) *
          fact( ( k - _m2 + _m1 ) / 2 ) );
 
    // fkw _coef[(k-_kmin)/2]=sign*sqrt(fact((_j+_m2)/2)*fact((_j-_m2)/2)*
    // fkw          fact((_j+_m1)/2)*fact((_j-_m1)/2))/
    // fkw  (fact((_j+_m2-k)/2)*fact(k/2)*fact((_j-_m1-k)/2)*fact((k-_m2+_m1)/2));
 
    // report(INFO,"EvtGen") << "k, coef:"<<k/2<<" "<<_coef[(k-_kmin)/2]<<endl;
  }
}

Referenced by EvtdFunction::d().

---

The documentation for this class was generated from the following files:

• EvtdFunctionSingle.hh
• EvtdFunctionSingle.cc