EvtdFunctionSingle Class Reference

#include <EvtdFunctionSingle.hh>

List of all members.

Public Member Functions
	EvtdFunctionSingle ()
	~EvtdFunctionSingle ()
void	init (int j, int m1, int m2)
double	d (int j, int m1, int m2, double theta)
Private Member Functions
int	fact (int n)
Private Attributes
int	_j
int	_m1
int	_m2
double *	_coef
int	_kmin
int	_kmax

---

Detailed Description

Definition at line 24 of file EvtdFunctionSingle.hh.

---

Constructor & Destructor Documentation

EvtdFunctionSingle::EvtdFunctionSingle

(

)

Definition at line 29 of file EvtdFunctionSingle.cc.

References _coef, _j, _kmax, _kmin, _m1, and _m2.

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

EvtdFunctionSingle::~EvtdFunctionSingle

(

)

Definition at line 39 of file EvtdFunctionSingle.cc.

References _coef.

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

---

Member Function Documentation

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

Definition at line 81 of file EvtdFunctionSingle.cc.

References _coef, _j, _kmax, _kmin, _m1, _m2, cos(), and sin().

Referenced by EvtdFunction::d().

                                                             {

  //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;
  
}

int EvtdFunctionSingle::fact

(

int

n

)

[private]

Definition at line 104 of file EvtdFunctionSingle.cc.

Referenced by init().

                                 {

  assert(n>=0);

  int f=1;

  int k;
  for(k=2;k<=n;k++) f*=k;

  return f;

}

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

Definition at line 44 of file EvtdFunctionSingle.cc.

References _coef, _j, _kmax, _kmin, _m1, _m2, abs, fact(), and sign().

Referenced by EvtdFunction::d().

                                                {

  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;
  }

}

---

Member Data Documentation

double* EvtdFunctionSingle::_coef [private]

Definition at line 43 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), init(), and ~EvtdFunctionSingle().

int EvtdFunctionSingle::_j [private]

Definition at line 39 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), and init().

int EvtdFunctionSingle::_kmax [private]

Definition at line 47 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), and init().

int EvtdFunctionSingle::_kmin [private]

Definition at line 46 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), and init().

int EvtdFunctionSingle::_m1 [private]

Definition at line 40 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), and init().

int EvtdFunctionSingle::_m2 [private]

Definition at line 41 of file EvtdFunctionSingle.hh.

Referenced by d(), EvtdFunctionSingle(), and init().

---