EvtEulerAngles Class Reference

#include <EvtEulerAngles.hh>

Public Member Functions
virtual	~EvtEulerAngles ()
	EvtEulerAngles ()
	EvtEulerAngles (const EvtVector3R &Yaxis, const EvtVector3R &Zaxis)
	EvtEulerAngles (const EvtVector4R &Pyaxis, const EvtVector4R &Pzaxis)
void	EulerAngles ()
double	getAlpha ()
double	getBeta ()
double	getGamma ()
virtual	~EvtEulerAngles ()
	EvtEulerAngles ()
	EvtEulerAngles (const EvtVector3R &Yaxis, const EvtVector3R &Zaxis)
	EvtEulerAngles (const EvtVector4R &Pyaxis, const EvtVector4R &Pzaxis)
void	EulerAngles ()
double	getAlpha ()
double	getBeta ()
double	getGamma ()

Detailed Description

Definition at line 32 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.hh.

Constructor & Destructor Documentation

~EvtEulerAngles() [1/2]

EvtEulerAngles::~EvtEulerAngles ( )

virtual

Definition at line 31 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

{}

EvtEulerAngles() [1/6]

EvtEulerAngles::EvtEulerAngles

(

)

Definition at line 46 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                              {
}

EvtEulerAngles() [2/6]

EvtEulerAngles::EvtEulerAngles	(	const EvtVector3R &	Yaxis,
		const EvtVector3R &	Zaxis
	)

Definition at line 32 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                                                                                  {
  _Yaxis=Yaxis;
  _Zaxis=Zaxis;
   EulerAngles();
}

EvtEulerAngles() [3/6]

EvtEulerAngles::EvtEulerAngles	(	const EvtVector4R &	Pyaxis,
		const EvtVector4R &	Pzaxis
	)

Definition at line 38 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                                                                                     {
  for (int i=1;i<4;i++){
    _Yaxis.set(i-1,Pyaxis.get(i));
    _Zaxis.set(i-1,Pzaxis.get(i));
  }
  EulerAngles();
}

~EvtEulerAngles() [2/2]

virtual EvtEulerAngles::~EvtEulerAngles ( )

virtual

EvtEulerAngles() [4/6]

EvtEulerAngles::EvtEulerAngles

(

)

EvtEulerAngles() [5/6]

EvtEulerAngles::EvtEulerAngles	(	const EvtVector3R &	Yaxis,
		const EvtVector3R &	Zaxis
	)

EvtEulerAngles() [6/6]

EvtEulerAngles::EvtEulerAngles	(	const EvtVector4R &	Pyaxis,
		const EvtVector4R &	Pzaxis
	)

Member Function Documentation

EulerAngles() [1/2]

void EvtEulerAngles::EulerAngles

(

)

Definition at line 61 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                                {
  // to calculate Euler angles with y-convention, i.e. R=R(Z, alpha).R(Y,beta).R(Z,gamma)
  double pi=3.1415926;
    _ry=_Yaxis.d3mag();
    _rz=_Zaxis.d3mag();
 
    if(_ry==0 ||_rz==0) {
      report(ERROR,"") << "Euler angle calculation specified by zero modules of the axis!"<<endl;
      report(ERROR,"EvtGen") << "Will terminate execution!"<<endl;
      ::abort();
    }
   double tolerance=1e-10;
   bool Y1is0=fabs(_Yaxis.get(0))<tolerance;
   bool Y2is0=fabs(_Yaxis.get(1))<tolerance;
   bool Y3is0=fabs(_Yaxis.get(2))<tolerance;
   bool Z1is0=fabs(_Zaxis.get(0))<tolerance;
   bool Z2is0=fabs(_Zaxis.get(1))<tolerance;
   bool Z3is0=fabs(_Zaxis.get(2))<tolerance;
 
   if(Y1is0 && Y3is0 && Z1is0 && Z2is0 ){
     _alpha=0; _beta=0; _gamma=0;  
   return;
   }
 
   if( Z1is0 && Z2is0 && !Y2is0){
      _alpha=0; _beta=0;
      _gamma=acos(_Yaxis.get(0)/_ry);
    if(_Yaxis.get(1)<0) _gamma=2*pi - _gamma;
      return;
   }
 
  // For general case to calculate Euler angles  
  // to calculate beta  
 
   if(Z1is0 && Z2is0) {
     _beta=0;
   } else{ _beta =acos(_Zaxis.get(2)/_rz);}
 
   //to calculate alpha
 
   if(_beta==0){
     _alpha=0.0;
   }else {
     double cosalpha=_Zaxis.get(0)/_rz/sin(_beta);
     if(fabs(cosalpha)>1.0) cosalpha=cosalpha/fabs(cosalpha);
     _alpha=acos(cosalpha);
     if(_Zaxis.get(1)<0.0) _alpha=2*pi - _alpha;
   }
 
   //to calculate gamma, alpha=0 and beta=0 case has been calculated, so only alpha !=0 and beta !=0 case left
 
   double singamma=_Yaxis.get(2)/_ry/sin(_beta);
   double cosgamma=(-_Yaxis.get(0)/_ry-cos(_alpha)*cos(_beta)*singamma)/sin(_alpha);
  if(fabs(singamma)>1.0) singamma=singamma/fabs(singamma);
  _gamma=asin(singamma);
  if(singamma>0 && cosgamma<0 ) _gamma=pi - _gamma;  // _gamma>0
  if(singamma<0 && cosgamma<0 ) _gamma=pi - _gamma;  //_gamma<0
  if(singamma<0 && cosgamma>0 ) _gamma=2*pi + _gamma; //_gamma<0
 
  
}

Referenced by EvtEulerAngles().

EulerAngles() [2/2]

void EvtEulerAngles::EulerAngles

(

)

getAlpha() [1/2]

double EvtEulerAngles::getAlpha

(

)

Definition at line 49 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                               {
  return _alpha;
}

Referenced by EvtJTO3P::decay().

getAlpha() [2/2]

double EvtEulerAngles::getAlpha

(

)

getBeta() [1/2]

double EvtEulerAngles::getBeta

(

)

Definition at line 53 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                              {
  return _beta;
}

Referenced by EvtJTO3P::decay().

getBeta() [2/2]

double EvtEulerAngles::getBeta

(

)

getGamma() [1/2]

double EvtEulerAngles::getGamma

(

)

Definition at line 57 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc.

                               {
  return _gamma;
}

Referenced by EvtJTO3P::decay().

getGamma() [2/2]

double EvtEulerAngles::getGamma

(

)

---

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

• bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.hh
• BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.hh
• bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc
• BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtEulerAngles.cc