00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021 #include "EvtGenBase/EvtPatches.hh"
00022 #include <stdlib.h>
00023 #include "EvtGenBase/EvtParticle.hh"
00024 #include "EvtGenBase/EvtGenKine.hh"
00025 #include "EvtGenBase/EvtPDL.hh"
00026 #include "EvtGenBase/EvtVector4R.hh"
00027 #include "EvtGenBase/EvtVector4C.hh"
00028 #include "EvtGenBase/EvtTensor4C.hh"
00029 #include "EvtGenBase/EvtTensor3C.hh"
00030 #include "EvtGenBase/EvtReport.hh"
00031 #include "EvtGenBase/EvtComplex.hh"
00032 #include "EvtGenModels/EvtSVPCP.hh"
00033 #include "EvtGenBase/EvtId.hh"
00034 #include "EvtGenBase/EvtCPUtil.hh"
00035 #include <string>
00036 #include "EvtGenBase/EvtConst.hh"
00037
00038 EvtSVPCP::~EvtSVPCP() {}
00039
00040 void EvtSVPCP::getName(std::string& model_name){
00041
00042 model_name="SVP_CP";
00043
00044 }
00045
00046
00047 EvtDecayBase* EvtSVPCP::clone(){
00048
00049 return new EvtSVPCP;
00050
00051 }
00052
00053 void EvtSVPCP::initProbMax(){
00054
00055 setProbMax(2*(getArg(3)*getArg(3)+getArg(5)*getArg(5)));
00056
00057 }
00058
00059
00060 void EvtSVPCP::init(){
00061
00062
00063 checkNArg(7);
00064 checkNDaug(2);
00065
00066 checkSpinParent(EvtSpinType::SCALAR);
00067
00068 checkSpinDaughter(0,EvtSpinType::VECTOR);
00069 checkSpinDaughter(1,EvtSpinType::PHOTON);
00070
00071 }
00072
00073 void EvtSVPCP::decay( EvtParticle *p ){
00074
00075 static EvtId B0=EvtPDL::getId("B0");
00076 static EvtId B0B=EvtPDL::getId("anti-B0");
00077
00078 double t;
00079 EvtId other_b;
00080
00081 EvtCPUtil::OtherB(p,t,other_b);
00082
00083 EvtComplex G1P,G1M, G1_T_even, G1_T_odd;
00084
00085 double norm = getArg(3)*getArg(3)+getArg(5)*getArg(5);
00086
00087 G1P=EvtComplex(getArg(3)*cos(getArg(4))/norm,getArg(3)*sin(getArg(4))/norm);
00088 G1M=EvtComplex(getArg(5)*cos(getArg(6))/norm,getArg(5)*sin(getArg(6))/norm);
00089
00090 G1_T_even = (G1P+G1M)/sqrt(2.0);
00091 G1_T_odd = (G1P-G1M)/sqrt(2.0);
00092
00093 EvtComplex lambda_km =EvtComplex(cos(-2*getArg(0)),sin(-2*getArg(0)));
00094
00095 double cdmt=cos(getArg(1)*t/(2*EvtConst::c));
00096 double sdmt=sin(getArg(1)*t/(2*EvtConst::c));
00097
00098 EvtComplex cG1_T_even,cG1_T_odd;
00099
00100 if (other_b==B0B){
00101 cG1_T_even = G1_T_even*(cdmt+lambda_km*EvtComplex(0.0,getArg(2)*sdmt));
00102 cG1_T_odd = G1_T_odd*(cdmt-lambda_km*EvtComplex(0.0,getArg(2)*sdmt));
00103 }
00104 if (other_b==B0){
00105 cG1_T_even = G1_T_even*(cdmt+(1.0/lambda_km)*EvtComplex(0.0,getArg(2)*sdmt));
00106 cG1_T_odd =-G1_T_odd*(cdmt-(1.0/lambda_km)*EvtComplex(0.0,getArg(2)*sdmt));
00107 }
00108
00109 EvtComplex hp, hm, h0;
00110
00111
00112
00113
00114 hm=(cG1_T_even-cG1_T_odd)/sqrt(2.0);
00115 hp=(cG1_T_even+cG1_T_odd)/sqrt(2.0);
00116 h0=EvtComplex(0.0,0.0);
00117
00118 EvtParticle *v1,*ph;
00119
00120 p->initializePhaseSpace(getNDaug(),getDaugs());
00121 v1 = p->getDaug(0);
00122 ph = p->getDaug(1);
00123 EvtVector4R momv1 = v1->getP4();
00124 EvtVector4R momph = ph->getP4();
00125
00126 EvtTensor4C d,g;
00127
00128 g.setdiag(1.0,-1.0,-1.0,-1.0);
00129
00130 EvtVector4R v,vp;
00131
00132 v=momv1/momv1.d3mag();
00133 vp=(momv1+momph)/(momv1+momph).mass();
00134
00135 d=((1.0/sqrt(3.0))*(h0-(hp+hm))*(-1.0/sqrt(3.0)))*g+
00136 ((1.0/sqrt(2.0))*(hp-hm)*EvtComplex(0.0,1.0)*(sqrt(1.0/2.0)))*dual(directProd(v,vp))+
00137 (sqrt(2.0/3.0)*(h0+0.5*(hp+hm))*sqrt(3.0/2.0))*(directProd(v,v)+(1.0/3.0)*g);
00138
00139 EvtVector4C ep0,ep1,ep2;
00140
00141 ep0=d.cont1(v1->eps(0).conj());
00142 ep1=d.cont1(v1->eps(1).conj());
00143 ep2=d.cont1(v1->eps(2).conj());
00144
00145 EvtVector4C ep20,ep21,ep22;
00146
00147 ep20=ph->epsParentPhoton(0).conj();
00148 ep21=ph->epsParentPhoton(1).conj();
00149
00150 vertex(0,0,ep0*ep20);
00151 vertex(0,1,ep0*ep21);
00152
00153 vertex(1,0,ep1*ep20);
00154 vertex(1,1,ep1*ep21);
00155
00156 vertex(2,0,ep2*ep20);
00157 vertex(2,1,ep2*ep21);
00158
00159
00160 return ;
00161
00162 }
00163