#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) |
EvtComplex EvtAmpA1 | ( | const EvtVector4R & | p4pi1, | |
const EvtVector4R & | p4pi2, | |||
const EvtVector4R & | p4pi3, | |||
const EvtVector4R & | p4pi4 | |||
) |
Definition at line 95 of file EvtBTo4piCP.cc.
References boostTo(), EvtVector4R::d3mag(), EvtVector4R::get(), EvtPDL::getId(), EvtPDL::getMeanMass(), EvtPDL::getWidth(), mass, t(), and EvtConst::twoPi.
Referenced by EvtBTo4piCP::decay().
00096 { 00097 00098 //added by Lange Jan4,2000 00099 static EvtId A1M=EvtPDL::getId("a_1-"); 00100 static EvtId RHO0=EvtPDL::getId("rho0"); 00101 00102 EvtVector4R p4a1,p4rho,p4b; 00103 00104 p4rho=p4pi1+p4pi2; 00105 00106 p4a1=p4rho+p4pi3; 00107 00108 p4b=p4a1+p4pi4; 00109 00110 EvtVector4R p4b_a1,p4rho_a1,p4pi1_a1,p4a1_a1; 00111 00112 p4b_a1=boostTo(p4b,p4a1); 00113 p4rho_a1=boostTo(p4rho,p4a1); 00114 p4pi1_a1=boostTo(p4pi1,p4a1); 00115 p4a1_a1=boostTo(p4a1,p4a1); 00116 00117 EvtVector4R p4pi1_rho; 00118 00119 p4pi1_rho=boostTo(p4pi1_a1,p4rho_a1); 00120 00121 EvtVector4R vb,vrho,vpi,t; 00122 00123 vb=p4b_a1/p4b_a1.d3mag(); 00124 vrho=p4rho_a1/p4rho_a1.d3mag(); 00125 vpi=p4pi1_rho/p4pi1_rho.d3mag(); 00126 00127 t.set(1.0,0.0,0.0,0.0); 00128 00129 EvtComplex amp_a1; 00130 00131 double bwm_a1=EvtPDL::getMeanMass(A1M); 00132 double gamma_a1=EvtPDL::getWidth(A1M); 00133 // double bwm_a2=EvtPDL::getMeanMass(A2M); 00134 // double gamma_a2=EvtPDL::getWidth(A2M); 00135 double bwm_rho=EvtPDL::getMeanMass(RHO0); 00136 double gamma_rho=EvtPDL::getWidth(RHO0); 00137 00138 amp_a1=(sqrt(gamma_a1/EvtConst::twoPi)/ 00139 ((p4a1).mass()-bwm_a1-EvtComplex(0.0,0.5*gamma_a1)))* 00140 (sqrt(gamma_rho/EvtConst::twoPi)/ 00141 ((p4rho).mass()-bwm_rho-EvtComplex(0.0,0.5*gamma_rho))); 00142 00143 return amp_a1* 00144 (vb.get(1)*vpi.get(1)+vb.get(2)*vpi.get(2)+vb.get(3)*vpi.get(3)); 00145 00146 }
EvtComplex EvtAmpA2 | ( | const EvtVector4R & | p4pi1, | |
const EvtVector4R & | p4pi2, | |||
const EvtVector4R & | p4pi3, | |||
const EvtVector4R & | p4pi4 | |||
) |
Definition at line 36 of file EvtBTo4piCP.cc.
References boostTo(), EvtVector4R::d3mag(), EvtVector4R::get(), EvtPDL::getId(), EvtPDL::getMeanMass(), EvtPDL::getWidth(), mass, t(), and EvtConst::twoPi.
Referenced by EvtBTo4piCP::decay().
00037 { 00038 00039 //added by Lange Jan4,2000 00040 static EvtId A2M=EvtPDL::getId("a_2-"); 00041 static EvtId RHO0=EvtPDL::getId("rho0"); 00042 00043 EvtVector4R p4a2,p4rho,p4b; 00044 00045 p4rho=p4pi1+p4pi2; 00046 00047 p4a2=p4rho+p4pi3; 00048 00049 p4b=p4a2+p4pi4; 00050 00051 EvtVector4R p4b_a2,p4rho_a2,p4pi1_a2,p4a2_a2; 00052 00053 p4b_a2=boostTo(p4b,p4a2); 00054 p4rho_a2=boostTo(p4rho,p4a2); 00055 p4pi1_a2=boostTo(p4pi1,p4a2); 00056 p4a2_a2=boostTo(p4a2,p4a2); 00057 00058 EvtVector4R p4pi1_rho; 00059 00060 p4pi1_rho=boostTo(p4pi1_a2,p4rho_a2); 00061 00062 EvtVector4R vb,vrho,vpi,t; 00063 00064 vb=p4b_a2/p4b_a2.d3mag(); 00065 vrho=p4rho_a2/p4rho_a2.d3mag(); 00066 vpi=p4pi1_rho/p4pi1_rho.d3mag(); 00067 00068 t.set(1.0,0.0,0.0,0.0); 00069 00070 // EvtComplex amp_a1,amp_a2; 00071 EvtComplex amp_a2; 00072 00073 // double bwm_a1=EvtPDL::getMeanMass(A1M); 00074 // double gamma_a1=EvtPDL::getWidth(A1M); 00075 double bwm_a2=EvtPDL::getMeanMass(A2M); 00076 double gamma_a2=EvtPDL::getWidth(A2M); 00077 double bwm_rho=EvtPDL::getMeanMass(RHO0); 00078 double gamma_rho=EvtPDL::getWidth(RHO0); 00079 00080 amp_a2=(sqrt(gamma_a2/EvtConst::twoPi)/ 00081 ((p4a2).mass()-bwm_a2-EvtComplex(0.0,0.5*gamma_a2)))* 00082 (sqrt(gamma_rho/EvtConst::twoPi)/ 00083 ((p4rho).mass()-bwm_rho-EvtComplex(0.0,0.5*gamma_rho))); 00084 00085 return amp_a2* 00086 (vb.get(1)*vrho.get(1)+vb.get(2)*vrho.get(2)+vb.get(3)*vrho.get(3))* 00087 ( 00088 vpi.get(1)*(vb.get(2)*vrho.get(3)-vb.get(3)*vrho.get(2))+ 00089 vpi.get(2)*(vb.get(3)*vrho.get(1)-vb.get(1)*vrho.get(3))+ 00090 vpi.get(3)*(vb.get(1)*vrho.get(2)-vb.get(2)*vrho.get(1)) 00091 ); 00092 00093 }