/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/Mcgpj/Mcgpj-00-01-04/src/code3pi/src/T2piCrossPart.C File Reference

#include "T2piCrossPart.h"
#include "TConstant.h"
#include <complex>
#include <iostream>

Go to the source code of this file.

Typedefs
typedef std::complex< double >	complex_t
Functions
complex_t	operator+ (const double &x, const complex_t &y)
complex_t	operator- (const double &x, const complex_t &y)
double	Gamma (const double Q2, const double M2, const double G)
complex_t	Rrho0 (double Q2)
complex_t	F2pi (double s)

---

Typedef Documentation

typedef std::complex<double> complex_t

Definition at line 7 of file T2piCrossPart.C.

---

Function Documentation

complex_t F2pi

(

double

s

)

Definition at line 46 of file T2piCrossPart.C.

Referenced by T2piCrossPart::SetJ().

                        {
  return complex_t(1,0);
  //  return Rrho0(s);
}

double Gamma	(	const double	Q2,
		const double	M2,
		const double	G
	)			[inline]

Definition at line 17 of file T2piCrossPart.C.

References rb::mpi2.

                                                                     {
  double beta2 = (Q2-4*mpi2)/(M2-4*mpi2);
  double beta = sqrt(beta2);
  return G*M2/Q2*beta2*beta;  
}

complex_t operator+	(	const double &	x,
		const complex_t &	y
	)			[inline]

Definition at line 9 of file T2piCrossPart.C.

                                                                {
  return y + x;
}

complex_t operator-	(	const double &	x,
		const complex_t &	y
	)			[inline]

Definition at line 13 of file T2piCrossPart.C.

                                                                {
  return -(y - x);
}

complex_t Rrho0

(

double

Q2

)

[inline]

Definition at line 23 of file T2piCrossPart.C.

References rb::Gamma(), rb::Gomega, rb::Grho, rb::Grhop, rb::Grhopp, I, rb::Momega, rb::Momega2, rb::mpi2, rb::Mrho2, rb::Mrhop2, and rb::Mrhopp2.

                                 {
  if(Q2<4*mpi2) return 0;
  const complex_t I(0,1);
  double Q = sqrt(Q2);
  complex_t prop_rho = Mrho2/(Q2 - Mrho2 + I*(Gamma(Q2,Mrho2,Grho)*Q));
  complex_t prop_rhop = Mrhop2/(Q2 - Mrhop2 + I*(Gamma(Q2,Mrhop2,Grhop)*Q));
  complex_t prop_rhopp = Mrhopp2/(Q2 - Mrhopp2 + I*(Gamma(Q2,Mrhopp2,Grhopp)*Q));
  // rho(770), rho(1450), rho(1700) in 2-pi inv mass
  // beta and gamma from Ph.D. thesis of Fedor Ignatov (ignatov@inp.nsk.su)

  //  const complex_t beta(0.41*cos(-2.02),0.41*sin(-2.02));
  const complex_t beta(-0.178042,-0.369325);

  //  const complex_t gamma(0.074*cos(-3.53),0.074*sin(-3.53));
  const complex_t gamma(-0.068488,0.0280249);

  // rho-omega mixing
  const complex_t delta(0.0015,0);
  complex_t prop_omega = Momega2/(Q2 - Momega2 + I*(Gomega*Momega));

  return (prop_rho*(1-delta*Q2/Momega2*prop_omega) + beta*prop_rhop  + gamma*prop_rhopp)/(1+beta+gamma);
}

---