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 59 of file UserDIY.cc.

Constructor & Destructor Documentation

◆ rhopi()

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

inline

Definition at line 61 of file UserDIY.cc.

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

Member Function Documentation

◆ amps()

double rhopi::amps ( )

Definition at line 106 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 75 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 80 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;
}