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 ()

Detailed Description

Definition at line 32 of file EvtEulerAngles.hh.

Constructor & Destructor Documentation

~EvtEulerAngles()

EvtEulerAngles::~EvtEulerAngles ( )

virtual

Definition at line 31 of file EvtEulerAngles.cc.

{}

EvtEulerAngles() [1/3]

EvtEulerAngles::EvtEulerAngles

(

)

Definition at line 47 of file EvtEulerAngles.cc.

{}

EvtEulerAngles() [2/3]

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

Definition at line 32 of file EvtEulerAngles.cc.

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

EvtEulerAngles() [3/3]

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

Definition at line 38 of file EvtEulerAngles.cc.

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

Member Function Documentation

EulerAngles()

void EvtEulerAngles::EulerAngles

(

)

Definition at line 55 of file 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(), and EvtEulerAngles().

getAlpha()

double EvtEulerAngles::getAlpha

(

)

Definition at line 49 of file EvtEulerAngles.cc.

{ return _alpha; }

Referenced by EvtJTO3P::decay().

getBeta()

double EvtEulerAngles::getBeta

(

)

Definition at line 51 of file EvtEulerAngles.cc.

{ return _beta; }

Referenced by EvtJTO3P::decay().

getGamma()

double EvtEulerAngles::getGamma

(

)

Definition at line 53 of file EvtEulerAngles.cc.

{ return _gamma; }

Referenced by EvtJTO3P::decay().

---

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

• EvtEulerAngles.hh
• EvtEulerAngles.cc