Public Member Functions
	rhopifull (EvtVector4R pd1, EvtVector4R pd2, EvtVector4R pd3)
double	Fij (int i, int j, double r)
double	R00 (double r)
EvtComplex	amps1 (int m, int i, int j)
double	amps ()

Detailed Description

Definition at line 117 of file UserDIY.cc.

Constructor & Destructor Documentation

◆ rhopifull()

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

inline

Definition at line 119 of file UserDIY.cc.

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

Member Function Documentation

◆ amps()

double rhopifull::amps ( )

Definition at line 172 of file UserDIY.cc.

                       {
  double temp = 0.0;
  int    m;
  for ( m = -1; m <= 1; m += 2 )
    temp = temp + pow( abs( amps1( m, 0, 1 ) + amps1( m, 0, 2 ) + amps1( m, 1, 2 ) ), 2. );
  return temp;
}

Referenced by EvtDIY::AmplitudeSquare().

◆ amps1()

EvtComplex rhopifull::amps1	(	int	m,
		int	i,
		int	j )

Definition at line 150 of file UserDIY.cc.

                                                 {
  double      mrho = 0.771, wrho = 0.1492, s;
  EvtComplex  img( 0.0, 1.0 ), dpro;
  EvtVector4R prho;
  prho = _pd[i] + _pd[j];
  s    = prho.mass2();
  dpro = s - mrho * mrho + img * sqrt( s ) * wrho;
  EvtHelSys angles( prho, _pd[i] ), labAngles;
  double    theta, phi, ct1, st1, phi1, st, ct;
  double    rpp = angles.getHelAng( 0 );
  theta         = angles.getHelAng( 1 );
  phi           = angles.getHelAng( 2 );
  ct1           = labAngles.Angles( prho, 1 );
  phi1          = labAngles.Angles( prho, 2 );
  int        lamb;
  EvtComplex temp( 0.0, 0.0 );
  for ( lamb = -1; lamb <= 1; lamb++ )
    temp = temp + Fij( lamb, 0, prho.d3mag() ) * Djmn( 1, m, lamb, phi1, ct1, 0.0 ) / dpro *
                      R00( rpp ) * Djmn( 1, lamb, 0, phi, theta, 0.0 );
  return temp;
}

Referenced by amps().

◆ Fij()

double rhopifull::Fij	(	int	i,
		int	j,
		double	r )

Definition at line 138 of file UserDIY.cc.

                                              {
  double mpi = 0.1395, mpsi = 3.096916;
  double temp = mpsi * r;
  if ( i == 0 && j == 0 ) return 0;
  if ( i == 1 && j == 0 ) return temp;
  if ( i == -1 && j == 0 ) return -temp;
 
  cerr << __FILE__ << ":" << __LINE__ << ": "
       << "Should not reach here!" << endl;
  abort();
}

Referenced by amps1().

◆ R00()

double rhopifull::R00 ( double r )

Definition at line 133 of file UserDIY.cc.

                                {
  double mpi = 0.1395;
  return r;
}