EvtHelSys Class Reference

#include <EvtHelSys.hh>

Public Member Functions
virtual	~EvtHelSys ()
	EvtHelSys ()
	EvtHelSys (const EvtVector4R &p4p, const EvtVector4R &p4d)
double	getHelAng (int i)
double	Angles (EvtVector4R, int)
EvtVector4R	checkparent ()
EvtVector4R	checkdaug ()
EvtVector4R	checkst (int i)
EvtVector4R	Helrotate (EvtVector4R p1, double phi, double theta)

Friends
double	djmn (int j, int m, int n, double theta)
double	djmn (double j, double m, double n, double theta)
EvtComplex	Djmn (int j, int m, int n, double phi, double theta, double gamma)
EvtComplex	Djmn (double j, double m, double n, double phi, double theta, double gamma)

Detailed Description

Definition at line 43 of file EvtHelSys.hh.

Constructor & Destructor Documentation

~EvtHelSys()

EvtHelSys::~EvtHelSys ( )

virtual

Definition at line 43 of file EvtHelSys.cc.

                      {
}

EvtHelSys() [1/2]

EvtHelSys::EvtHelSys

(

)

Definition at line 51 of file EvtHelSys.cc.

                       {
 }

EvtHelSys() [2/2]

EvtHelSys::EvtHelSys	(	const EvtVector4R &	p4p,
		const EvtVector4R &	p4d )

Definition at line 46 of file EvtHelSys.cc.

                                                                     {
   _p4parent=p4p;   // parent momentum in its parent CM system
   _p4daug  =p4d;   // daughter momentum in its parent CM system
 }

Member Function Documentation

Angles()

double EvtHelSys::Angles	(	EvtVector4R	bbst,
		int	i )

Definition at line 108 of file EvtHelSys.cc.

                                              {
 
double rxy=sqrt(pow(bbst.get(1),2.)+pow(bbst.get(2),2.));
  if(bbst.d3mag()<=1e-10) {
      if(i==0) return 0.;
      else if (i==1) return 0.;
      else if (i==2) return 0.;
      else {cout<<"Angles(i): i<=2"<<endl;abort();}
   }
  else if(rxy<=1e-10){
      if(i==0) return bbst.d3mag();
      if(i==1) return 0.;
      if(i==2) return 0.;
      else {cout<<"Angles(i): i<=2"<<endl;abort();}
   }
   else {
   double theta=acos(bbst.get(3)/bbst.d3mag());
   double csphi=bbst.get(1)/bbst.d3mag()/sin(theta);
   if(fabs(csphi)>1.0) csphi=csphi/fabs(csphi);
   double phi=acos(csphi);
   if(bbst.get(2)<0.0) phi=2*3.1415926-phi;
      if(i==0) return bbst.d3mag();
      else if (i==1) return theta;
      else if (i==2) return phi;
      else {cout<<"Angles(i): i<=2"<<endl;abort();}
    }
    
} 

Referenced by rhopi::amps1(), rhopifull::amps1(), VVS::amps1(), and getHelAng().

checkdaug()

EvtVector4R EvtHelSys::checkdaug

(

)

Definition at line 86 of file EvtHelSys.cc.

{ return _p4daug;}

checkparent()

EvtVector4R EvtHelSys::checkparent

(

)

Definition at line 85 of file EvtHelSys.cc.

{ return _p4parent;}

checkst()

EvtVector4R EvtHelSys::checkst

(

int

i

)

Definition at line 87 of file EvtHelSys.cc.

                                   {
getHelAng(1);
if(i==0) return _bp4p;      //parent momentum used to boost the daughter to the CM sys.  
if(i==1) return _rotatep4p; //the parent momentum in Hel system by rotation
if(i==2) return _rotatep4d; //the daughter momentum in Hel. system by rotation
if(i==3) return _bst;       //_bst:daughter momentum in helicity system
}

getHelAng()

double EvtHelSys::getHelAng

(

int

i

)

Definition at line 54 of file EvtHelSys.cc.

                                    {
  EvtVector4R b_p4p,rp4p, rp4d,boostdaug; 
  EvtVector3R GetHelAng;
  while (_p4parent.d3mag()!=0) {
 
    //  b_p4p=-1 * _p4parent;   //boost from Lab to HEL sys. required to reverse mom.Vec.
    //  b_p4p.set(0,_p4parent.get(0));  
    //  _bp4p=b_p4p;
 
// first to rotate the mother and daugher momentum to the helicity system
  double theta=Angles(_p4parent,1);
  double phi=Angles(_p4parent,2);
 
 
  rp4p=Helrotate(_p4parent,phi,theta);
  rp4d=Helrotate(_p4daug,  phi,theta);
 
// then boos to the CM system
//  EvtVector4R r_p4p=-1*rp4p;  //boost from Lab to HEL sys. required to reverse mom.Vec.
//  r_p4p.set(0,rp4p.get(0));
//  boostdaug=boostTo(rp4d,r_p4p);
 
  _rotatep4p=rp4p;
  _rotatep4d=rp4d;
  _bst=rp4d;
  //  _bp4p=r_p4p;
  return Angles(_bst,i);  //_bst:daughter momentum in helicity system, i=0==>|_bst|;i=1,2==>(theta,phi)
  }
 return Angles(_p4daug,i);
}

Referenced by rhopi::amps1(), rhopifull::amps1(), VVS::amps1(), EvtPsi3Sdecay::AngSam(), EvtConExc::baryon_sampling(), checkst(), EvtAngSam::decay(), EvtAngSamX::decay(), EvtAV2GV::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtDeBD::decay(), EvtHelPPJ::decay(), EvtHypWK::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJPE::decay(), EvtP2GC0::decay(), EvtP2GC1::decay(), EvtP2GC2::decay(), EvtSPL::decay(), EvtT2GV::decay(), EvtEvalHelAmp::evalAmp(), and EvtConExc::meson_sampling().

Helrotate()

EvtVector4R EvtHelSys::Helrotate	(	EvtVector4R	p1,
		double	phi,
		double	theta )

Definition at line 96 of file EvtHelSys.cc.

                                                                        {
EvtVector4R Rp;
double cp=cos(phi);
double sp=sin(phi);
double ct=cos(theta);
double st=sin(theta);
double t=p1.get(0),x=p1.get(1), y=p1.get(2),   z=p1.get(3);
double xp=x*cp*ct+y*sp*ct-z*st,yp= -x*sp+y*cp,zp=x*cp*st+y*sp*st+z*ct;
      Rp.set(t,xp,yp,zp);
 return Rp;
}

Referenced by getHelAng().

Friends And Related Symbol Documentation

Djmn [1/2]

EvtComplex Djmn ( double j, double m, double n, double phi, double theta, double gamma )

friend

Definition at line 165 of file EvtHelSys.cc.

                                                                                    {
int j2=(int)(j*2*1.1),m2=(int)(m*2*1.1),n2=(int)(n*2*1.1);
EvtComplex gp(cos(-phi*m  ), -sin(phi*m));
EvtComplex gm(cos(-gamma*n), -sin(gamma*n)); 
double  tp3=EvtdFunction::d(j2,m2,n2,theta);
EvtComplex temp=gp * tp3 * gm;
return temp;
}

Djmn [2/2]

EvtComplex Djmn ( int j, int m, int n, double phi, double theta, double gamma )

friend

Definition at line 151 of file EvtHelSys.cc.

                                                                           {
int j2=j*2,m2=m*2,n2=n*2;
EvtComplex gp(cos(-phi*m  ), -sin(phi*m));
EvtComplex gm(cos(-gamma*n), -sin(gamma*n));
double  tp3=EvtdFunction::d(j2,m2,n2,theta);
 
//EvtComplex temp=wignerD(j2,m2,n2,phi,theta,gamma); //wignerD is corrected by pingrg, 2007,04,28, it gives the same result as this definition
 
EvtComplex temp=gp * tp3 * gm;
 
return temp;
}

djmn [1/2]

double djmn ( double j, double m, double n, double theta )

friend

Definition at line 144 of file EvtHelSys.cc.

                                                       {  
int j2=(int)(j*2*1.1),m2=(int)(m*2*1.1),n2=(int)(n*2*1.1);
double temp=EvtdFunction::d(j2,m2,n2,theta);
return temp;
} 

djmn [2/2]

double djmn ( int j, int m, int n, double theta )

friend

Definition at line 138 of file EvtHelSys.cc.

                                              {  
int j2=j*2,m2=m*2,n2=n*2;
double temp=EvtdFunction::d(j2,m2,n2,theta);
return temp;
}

---

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

• EvtHelSys.hh
• EvtHelSys.cc