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

#include "T3piCrossPart.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	Rrho (double Q2, bool IsIsospinBreak)
complex_t	Amp (double s)
complex_t	AmpOmegaPhi (double s)
complex_t	F3pi (double s, double Q02, double Qp2, double Qm2)

---

Typedef Documentation

typedef std::complex<double> complex_t

Definition at line 8 of file T3piCrossPart.C.

---

Function Documentation

complex_t Amp

(

double

s

)

[inline]

Definition at line 55 of file T3piCrossPart.C.

References alpha, exp(), rb::Geepsi, rb::Geepsip, rb::Geepsipp, rb::Gpsi, rb::Gpsip, rb::Gpsipp, I, M_PI, rb::Mpsi, rb::Mpsi2, rb::Mpsip, rb::Mpsip2, rb::Mpsipp, and rb::Mpsipp2.

Referenced by EvtSVSCPiso::decay(), F3pi(), FK0K0Eta(), FKcKcEta(), FKcKcP0(), FRhoEta(), and T3piEtaCrossPart::SetJ().

                              {
  const complex_t I(0,1);
  const double Apsi = 4*alpha*Gpsi/(3*Geepsi);
  const double phipsi = -M_PI/2;
  const double Apsip = alpha*Gpsip/(3*Geepsip);
  const double phipsip = -M_PI/2;
  const double Apsipp = alpha*Gpsipp/(3*Geepsipp);
  const double phipsipp = -M_PI/2;
  complex_t res = 3*sqrt(s)/alpha*
    (
     Geepsi*(1 + Apsi*exp(I*phipsi))/(s - Mpsi2 + I*Mpsi*Gpsi) +
     Geepsip*(1 + Apsip*exp(I*phipsip))/(s - Mpsip2 + I*Mpsip*Gpsip) +
     Geepsipp*(1 + Apsipp*exp(I*phipsipp))/(s - Mpsipp2 + I*Mpsipp*Gpsipp));
  return res;
}

complex_t AmpOmegaPhi

(

double

s

)

[inline]

Definition at line 71 of file T3piCrossPart.C.

References rb::Gomega, rb::Gphi, rb::Momega, rb::Momega2, rb::Mphi, and rb::Mphi2.

                                      {
  //  const double eta = 3.4*M_PI/180;
  //  const double theta = eta + asin(1/sqrt(3));
  //  const double comega = sin(theta)*cos(eta);
  //  const double cphi = -cos(theta)*sin(eta);
  const double comega = 6.23657803940740774e-01;
  const double cphi = -4.63075347511148622e-02;
  complex_t Romega = Momega2/complex_t(s-Momega2,Gomega*Momega);
  complex_t Rphi = Mphi2/complex_t(s-Mphi2,Gphi*Mphi);
  return (comega*Romega + cphi*Rphi);
}

complex_t F3pi	(	double	s,
		double	Q02,
		double	Qp2,
		double	Qm2
	)

Definition at line 83 of file T3piCrossPart.C.

References Amp(), H, and Rrho().

Referenced by T3piCrossPart::SetJ().

                                                            {
  // Rrho has rho-omega mixing otherwise use Rrho
  complex_t H = Rrho(Q02,true) + Rrho(Qm2,false) + Rrho(Qp2,false);
  return (5e-3*(1/s))*H*(1+Amp(s));
  //  const double f_pi3 = 93*93*93;// [MeV^3]
  //  return complex_t(1000/(s*sqrt(s)),0);
  //  const double alphaK = 0.5;
  //  return 2e5/(s*s)*(1+Amp(s))*(1-3*alphaK-alphaK*H);
  //  return 0.325*1e-9*Amp(s)*H;
  //  return sqrt(3)/(4*M_PI*M_PI*f_pi3)*(AmpOmegaPhi(s)+Amp(s))*(1-3*alphaK-alphaK*H);
  //  return sqrt(3)/(4*M_PI*M_PI*f_pi3)*(1)/(s)*(1-3*alphaK-alphaK*H);
}

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

Definition at line 18 of file T3piCrossPart.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 10 of file T3piCrossPart.C.

                                                                {
  return y + x;
}

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

Definition at line 14 of file T3piCrossPart.C.

                                                                {
  return -(y - x);
}

complex_t Rrho	(	double	Q2,
		bool	IsIsospinBreak
	)			[inline]

Definition at line 25 of file T3piCrossPart.C.

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

Referenced by F3pi(), FRhoEta(), J5pi(), and T4piCrossPart::SetJ().

                                                     {
  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));
  return prop_rho;
  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));
  complex_t prop_rhoppp = Mrhoppp2/(Q2 - Mrhoppp2 + I*(Gamma(Q2,Mrhoppp2,Grhoppp)*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);

  const complex_t xi(0.05,0.0);
  if(IsIsospinBreak){
    // rho-omega mixing
    const complex_t delta(0.002,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 + xi*prop_rhoppp)/(1+beta+gamma+xi);
  } else {
    return (prop_rho + beta*prop_rhop  + gamma*prop_rhopp + xi*prop_rhoppp)/(1+beta+gamma+xi);
  }
}

---