00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #include "EvtGenBase/EvtPatches.hh"
00024 #include <stdlib.h>
00025 #include <iostream>
00026 #include <string>
00027 #include "EvtGenBase/EvtParticle.hh"
00028 #include "EvtGenBase/EvtPDL.hh"
00029 #include "EvtGenBase/EvtGenKine.hh"
00030 #include "EvtGenModels/EvtVll.hh"
00031 #include "EvtGenBase/EvtDiracSpinor.hh"
00032 #include "EvtGenBase/EvtReport.hh"
00033 #include "EvtGenBase/EvtVector4C.hh"
00034
00035 EvtVll::~EvtVll() {}
00036
00037 void EvtVll::getName(std::string& model_name){
00038
00039 model_name="VLL";
00040
00041 }
00042
00043
00044 EvtDecayBase* EvtVll::clone(){
00045
00046 return new EvtVll;
00047
00048 }
00049
00050 void EvtVll::init(){
00051
00052
00053 checkNArg(0);
00054 checkNDaug(2);
00055
00056 checkSpinParent(EvtSpinType::VECTOR);
00057
00058 checkSpinDaughter(0,EvtSpinType::DIRAC);
00059 checkSpinDaughter(1,EvtSpinType::DIRAC);
00060
00061 }
00062
00063 void EvtVll::initProbMax(){
00064
00065 setProbMax(1.0);
00066
00067 }
00068
00069 void EvtVll::decay(EvtParticle *p){
00070
00071 p->initializePhaseSpace(getNDaug(),getDaugs());
00072
00073 EvtParticle *l1, *l2;
00074 l1 = p->getDaug(0);
00075 l2 = p->getDaug(1);
00076
00077 EvtVector4C l11, l12, l21, l22;
00078 l11=EvtLeptonVCurrent(l1->spParent(0),l2->spParent(0));
00079 l12=EvtLeptonVCurrent(l1->spParent(0),l2->spParent(1));
00080 l21=EvtLeptonVCurrent(l1->spParent(1),l2->spParent(0));
00081 l22=EvtLeptonVCurrent(l1->spParent(1),l2->spParent(1));
00082
00083 EvtVector4C eps0=p->eps(0);
00084 EvtVector4C eps1=p->eps(1);
00085 EvtVector4C eps2=p->eps(2);
00086
00087 double M2=p->mass();
00088 M2*=M2;
00089 double m2=l1->mass();
00090 m2*=m2;
00091
00092 double norm=1.0/sqrt(2*M2+4*m2-4*m2*m2/M2);
00093
00094 vertex(0,0,0,norm*(eps0*l11));
00095 vertex(0,0,1,norm*(eps0*l12));
00096 vertex(0,1,0,norm*(eps0*l21));
00097 vertex(0,1,1,norm*(eps0*l22));
00098
00099 vertex(1,0,0,norm*(eps1*l11));
00100 vertex(1,0,1,norm*(eps1*l12));
00101 vertex(1,1,0,norm*(eps1*l21));
00102 vertex(1,1,1,norm*(eps1*l22));
00103
00104 vertex(2,0,0,norm*(eps2*l11));
00105 vertex(2,0,1,norm*(eps2*l12));
00106 vertex(2,1,0,norm*(eps2*l21));
00107 vertex(2,1,1,norm*(eps2*l22));
00108
00109 return;
00110
00111 }
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124