Public Member Functions
	rhopi (EvtVector4R pd1, EvtVector4R pd2, EvtVector4R pd3)

double	F00 (double s)

double	F10 (double s)

double	amps1 (double s, int i, int j)

double	amps ()

Detailed Description

Definition at line 60 of file UserDIY.cc.

Constructor & Destructor Documentation

◆ rhopi()

rhopi::rhopi	(	EvtVector4R	pd1,
		EvtVector4R	pd2,
		EvtVector4R	pd3
	)

inline

Definition at line 62 of file UserDIY.cc.

                                                        {
 _pd[0]=pd1;
 _pd[1]=pd2;
 _pd[2]=pd3;
 }

Member Function Documentation

◆ amps()

double rhopi::amps ( )

Definition at line 104 of file UserDIY.cc.

                   {
 double temp,s12,s13,s23;
 s12=(_pd[0]+_pd[1]).mass2();
 s13=(_pd[0]+_pd[2]).mass2();
 s23=(_pd[1]+_pd[2]).mass2();
 temp=amps1(s12,0,1)+amps1(s13,0,2)+amps1(s23,1,2);
 return temp;
 } 

◆ amps1()

double rhopi::amps1	(	double	s,
		int	i,
		int	j
	)

Definition at line 87 of file UserDIY.cc.

                                          {
   double mrho=0.771,wrho=0.1492,dpro;
   EvtComplex img(0.0,1.0);
   dpro=pow(abs(s-mrho*mrho+img*sqrt(s)*wrho),2.);
   EvtVector4R prho;
   prho=_pd[i]+_pd[j];
   EvtHelSys angles(prho,_pd[i]),labAngles;
   double theta,phi,ct1,st1,phi1,st,ct,temp;
   theta=angles.getHelAng(1);
   phi  =angles.getHelAng(2);
   ct1 =labAngles.Angles(prho,1);
   phi1=labAngles.Angles(prho,2);
   st=sin(theta);ct=cos(theta);
   temp=pow(F00(s),2.)*pow(F10(s),2.)*pow(st,2.)/dpro; // *(1+pow(ct1,2.)+pow(st1,2.)*cos(2*(phi1+phi))); 
    return temp;
  }

Referenced by amps().

◆ F00()

double rhopi::F00 ( double s )

Definition at line 76 of file UserDIY.cc.

                           {
  double mpi=0.1395;
  return sqrt(s-4*mpi*mpi)/sqrt(s);
  }

Referenced by amps1().

◆ F10()

double rhopi::F10 ( double s )

Definition at line 81 of file UserDIY.cc.

                            {
    double mpi=0.1395,mpsi=3.096916;
    double tep=sqrt((mpsi*mpsi-pow(mpi+sqrt(s),2.))*(mpsi*mpsi-pow(mpi-sqrt(s),2.)));
   return tep;
  }