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 39 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 81 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 44 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