EvtBTo4piCP.cc File Reference

#include "EvtGenBase/EvtPatches.hh"
#include <stdlib.h>
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtGenKine.hh"
#include "EvtGenBase/EvtCPUtil.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenModels/EvtBTo4piCP.hh"
#include "EvtGenBase/EvtId.hh"
#include "EvtGenBase/EvtConst.hh"
#include <string>

Go to the source code of this file.

Functions
EvtComplex	EvtAmpA2 (const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
EvtComplex	EvtAmpA1 (const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)

Function Documentation

EvtAmpA1()

EvtComplex EvtAmpA1	(	const EvtVector4R &	p4pi1,
		const EvtVector4R &	p4pi2,
		const EvtVector4R &	p4pi3,
		const EvtVector4R &	p4pi4 )

Definition at line 95 of file EvtBTo4piCP.cc.

                                                              {
 
  //added by Lange Jan4,2000
  static EvtId A1M=EvtPDL::getId("a_1-");
  static EvtId RHO0=EvtPDL::getId("rho0");
 
  EvtVector4R p4a1,p4rho,p4b;
 
  p4rho=p4pi1+p4pi2;
 
  p4a1=p4rho+p4pi3;
 
  p4b=p4a1+p4pi4;
 
  EvtVector4R p4b_a1,p4rho_a1,p4pi1_a1,p4a1_a1;
 
  p4b_a1=boostTo(p4b,p4a1);
  p4rho_a1=boostTo(p4rho,p4a1);
  p4pi1_a1=boostTo(p4pi1,p4a1);
  p4a1_a1=boostTo(p4a1,p4a1);
 
  EvtVector4R p4pi1_rho;
 
  p4pi1_rho=boostTo(p4pi1_a1,p4rho_a1);
 
  EvtVector4R vb,vrho,vpi,t;
 
  vb=p4b_a1/p4b_a1.d3mag();
  vrho=p4rho_a1/p4rho_a1.d3mag();
  vpi=p4pi1_rho/p4pi1_rho.d3mag();
 
  t.set(1.0,0.0,0.0,0.0);
 
  EvtComplex amp_a1;
 
  double bwm_a1=EvtPDL::getMeanMass(A1M);
  double gamma_a1=EvtPDL::getWidth(A1M);
  //  double bwm_a2=EvtPDL::getMeanMass(A2M);
  //  double gamma_a2=EvtPDL::getWidth(A2M);
  double bwm_rho=EvtPDL::getMeanMass(RHO0);
  double gamma_rho=EvtPDL::getWidth(RHO0);
 
  amp_a1=(sqrt(gamma_a1/EvtConst::twoPi)/
    ((p4a1).mass()-bwm_a1-EvtComplex(0.0,0.5*gamma_a1)))*
         (sqrt(gamma_rho/EvtConst::twoPi)/
    ((p4rho).mass()-bwm_rho-EvtComplex(0.0,0.5*gamma_rho)));
 
  return amp_a1*
    (vb.get(1)*vpi.get(1)+vb.get(2)*vpi.get(2)+vb.get(3)*vpi.get(3));
 
}

Referenced by EvtBTo4piCP::decay().

EvtAmpA2()

EvtComplex EvtAmpA2	(	const EvtVector4R &	p4pi1,
		const EvtVector4R &	p4pi2,
		const EvtVector4R &	p4pi3,
		const EvtVector4R &	p4pi4 )

Definition at line 36 of file EvtBTo4piCP.cc.

                                                              {
 
  //added by Lange Jan4,2000
  static EvtId A2M=EvtPDL::getId("a_2-");
  static EvtId RHO0=EvtPDL::getId("rho0");
 
  EvtVector4R p4a2,p4rho,p4b;
 
  p4rho=p4pi1+p4pi2;
 
  p4a2=p4rho+p4pi3;
 
  p4b=p4a2+p4pi4;
 
  EvtVector4R p4b_a2,p4rho_a2,p4pi1_a2,p4a2_a2;
 
  p4b_a2=boostTo(p4b,p4a2);
  p4rho_a2=boostTo(p4rho,p4a2);
  p4pi1_a2=boostTo(p4pi1,p4a2);
  p4a2_a2=boostTo(p4a2,p4a2);
 
  EvtVector4R p4pi1_rho;
 
  p4pi1_rho=boostTo(p4pi1_a2,p4rho_a2);
 
  EvtVector4R vb,vrho,vpi,t;
 
  vb=p4b_a2/p4b_a2.d3mag();
  vrho=p4rho_a2/p4rho_a2.d3mag();
  vpi=p4pi1_rho/p4pi1_rho.d3mag();
 
  t.set(1.0,0.0,0.0,0.0);
 
  //  EvtComplex amp_a1,amp_a2;
  EvtComplex amp_a2;
 
  //  double bwm_a1=EvtPDL::getMeanMass(A1M);
  //  double gamma_a1=EvtPDL::getWidth(A1M);
  double bwm_a2=EvtPDL::getMeanMass(A2M);
  double gamma_a2=EvtPDL::getWidth(A2M);
  double bwm_rho=EvtPDL::getMeanMass(RHO0);
  double gamma_rho=EvtPDL::getWidth(RHO0);
 
  amp_a2=(sqrt(gamma_a2/EvtConst::twoPi)/
    ((p4a2).mass()-bwm_a2-EvtComplex(0.0,0.5*gamma_a2)))*
         (sqrt(gamma_rho/EvtConst::twoPi)/
    ((p4rho).mass()-bwm_rho-EvtComplex(0.0,0.5*gamma_rho)));
 
  return amp_a2*
    (vb.get(1)*vrho.get(1)+vb.get(2)*vrho.get(2)+vb.get(3)*vrho.get(3))*
    (
     vpi.get(1)*(vb.get(2)*vrho.get(3)-vb.get(3)*vrho.get(2))+
     vpi.get(2)*(vb.get(3)*vrho.get(1)-vb.get(1)*vrho.get(3))+
     vpi.get(3)*(vb.get(1)*vrho.get(2)-vb.get(2)*vrho.get(1))
     );
 
}

Referenced by EvtBTo4piCP::decay().