#include <EvtVVS.hh>

Public Member Functions
	VVS (EvtVector4R pd1, EvtVector4R pd2, EvtVector4R pd3, double ResonanceMass, double ResonanceWidth, double r1, double r2, double phase1, double phase2)
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 34 of file EvtVVS.hh.

Constructor & Destructor Documentation

◆ VVS()

VVS::VVS	(	EvtVector4R	pd1,
		EvtVector4R	pd2,
		EvtVector4R	pd3,
		double	ResonanceMass,
		double	ResonanceWidth,
		double	r1,
		double	r2,
		double	phase1,
		double	phase2 )

inline

Definition at line 36 of file EvtVVS.hh.

                                                                                   {
    ////////////////////////////////////////////////////////////////////////////////////
    //     For the decay 1-- -->R(1--) +Pseudoscalar-->3 Pseudoscalars
    //  pd1,pd2,pd3: four momentum vector for three daughters
    //   ResonanceMass,  ResonanceWidth : resonance parameter for intermediate state R(1--)
    //  r1,r2      : the relative coupling strength for two resonaces R R.P.T another R
    //  phase1,phase2: the relatev phase of two phases R.P.T another R
    ///////////////////////////////////////////////////////////////////////////////////////////
    _pd[0]  = pd1;
    _pd[1]  = pd2;
    _pd[2]  = pd3;
    _res[0] = ResonanceMass;
    _res[1] = ResonanceWidth;
    _par[0] = r1;
    _par[1] = phase1;
    _par[2] = r2;
    _par[3] = phase2;
  }

Member Function Documentation

◆ amps()

double VVS::amps ( )

Definition at line 71 of file EvtVVS.cc.

                 {
  double r1     = _par[0];
  double phase1 = _par[1];
  double r2     = _par[2];
  double phase2 = _par[3];
  double temp   = 0.0;
  int    m;
  for ( m = -1; m <= 1; m += 2 )
    temp = temp + pow( abs( r1 * exp( EvtComplex( 0, phase1 ) ) * amps1( m, 0, 1 ) +
                            r2 * exp( EvtComplex( 0, phase2 ) ) * amps1( m, 0, 2 ) +
                            amps1( m, 1, 2 ) ),
                       2. );
  return temp;
}

Referenced by EvtRhoPi::AmplitudeSquare().

◆ amps1()

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

Definition at line 47 of file EvtVVS.cc.

                                           {
  double      mrho = _res[0];
  double      wrho = _res[1];
  double      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 VVS::Fij	(	int	i,
		int	j,
		double	r )

Definition at line 37 of file EvtVVS.cc.

                                        {
  double mpsi = ( _pd[0] + _pd[1] + _pd[2] ).mass();
  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;
  std::cerr << __FILE__ << ":" << __LINE__ << ": should not reach here" << std::endl;
  exit( 1 );
}

Referenced by amps1().

◆ R00()

double VVS::R00 ( double r )

Definition at line 35 of file EvtVVS.cc.

35{ return r; }