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

#include <iostream>
#include "T5piCrossPart.h"
#include "TConstant.h"
#include "TUtil.h"
#include <complex>

Go to the source code of this file.

Typedefs
typedef std::complex< double >	complex_t
Functions
complex_t	Rrho (double s, bool IsIsospinBreak)
complex_t	bw_om (double s)
complex_t	bw_b1 (double s)
TLorentzVectorC	J5pi (const TLorentzVector &q1, const TLorentzVector &q2, const TLorentzVector &q3, const TLorentzVector &q4, const TLorentzVector &q5)

---

Typedef Documentation

typedef std::complex<double> complex_t

Definition at line 8 of file T5piCrossPart.C.

---

Function Documentation

complex_t bw_b1

(

double

s

)

[inline]

Definition at line 24 of file T5piCrossPart.C.

References rb::Gb1, rb::Mb1, and rb::Mb12.

Referenced by J5pi().

                                {
  complex_t Db1 = complex_t(Mb12-s,-Mb1*Gb1);
  complex_t bw = Mb12/Db1;
  return bw;
}

complex_t bw_om

(

double

s

)

[inline]

Definition at line 18 of file T5piCrossPart.C.

References rb::Gomega, rb::Momega, and rb::Momega2.

                                {
  complex_t Domega = complex_t(Momega2-s,-Momega*Gomega);
  complex_t bw = Momega2/Domega;
  return bw;
}

TLorentzVectorC J5pi	(	const TLorentzVector &	q1,
		const TLorentzVector &	q2,
		const TLorentzVector &	q3,
		const TLorentzVector &	q4,
		const TLorentzVector &	q5
	)

Definition at line 80 of file T5piCrossPart.C.

References bw_b1(), bw_om(), rb::PseudoScalars3(), and Rrho().

Referenced by T5piCrossPart::SetJ().

                                              {
  TLorentzVector qom = q3 + q4 + q5;
  TLorentzVector qp = q3 + q5;
  TLorentzVector qm = q4 + q5;
  TLorentzVector q0 = q3 + q4;
  TLorentzVector qb1_0 = qom + q2;
  TLorentzVector qb1_1 = qom + q1;
  complex_t pom = bw_om(qom.M2());
  pom *= Rrho(qp.M2(),false) + Rrho(qm.M2(),false) + Rrho(q0.M2(),true);

  TLorentzVector tom = PseudoScalars3(q3,q4,q5);

  complex_t pb1_0 = bw_b1(qb1_0.M2());
  complex_t pb1_1 = bw_b1(qb1_1.M2());

  complex_t prop_0 = pom*pb1_0;
  complex_t prop_1 = pom*pb1_1;
  double    q1q2   =    q1*q2;
  double    q1tom  =    q1*tom;
  double qb1_0qom  = qb1_0*qom;
  double qb1_0tom  = qb1_0*tom;
  double qb1_0q2   = qb1_0*q2;
  double    q1qom  =    q1*qom;
  TLorentzVector res0
    = qom*((   q1q2 )*(qb1_0tom) - (qb1_0q2 )*(   q1tom))
    +  q2*((qb1_0qom)*(   q1tom) - (   q1qom)*(qb1_0tom))
    + tom*((   q1qom)*(qb1_0q2 ) - (qb1_0qom)*(   q1q2 ));

  double    q2tom  =    q2*tom;
  double qb1_1qom  = qb1_1*qom;
  double qb1_1tom  = qb1_1*tom;
  double qb1_1q1   = qb1_1*q1;
  double    q2qom  =    q2*qom;
  TLorentzVector res1
    = qom*((   q1q2 )*(qb1_1tom) - (qb1_1q1 )*(   q2tom))
    +  q1*((qb1_1qom)*(   q2tom) - (   q2qom)*(qb1_1tom))
    + tom*((   q2qom)*(qb1_1q1 ) - (qb1_1qom)*(   q1q2 ));
  
  return TLorentzVectorC(res0.X()*prop_0 + res1.X()*prop_1,
                         res0.Y()*prop_0 + res1.Y()*prop_1,
                         res0.Z()*prop_0 + res1.Z()*prop_1,
                         res0.T()*prop_0 + res1.T()*prop_1);
}

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

Definition at line 12 of file T5piCrossPart.C.

References rb::Grho, rb::Mrho, and rb::Mrho2.

                                                    {
  complex_t Drho = complex_t(Mrho2-s, -Mrho*Grho);
  complex_t prop_rho = Mrho2/Drho;
  return prop_rho;
}

---