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 <fstream>
00024 #include <stdio.h>
00025 #include <string>
00026 #include "EvtGenBase/EvtGenKine.hh"
00027 #include "EvtGenBase/EvtParticle.hh"
00028 #include "EvtGenBase/EvtPDL.hh"
00029 #include "EvtGenBase/EvtReport.hh"
00030 #include "EvtGenModels/EvtPi0Dalitz.hh"
00031 #include "EvtGenBase/EvtVector4C.hh"
00032 #include "EvtGenBase/EvtDiracSpinor.hh"
00033 #include "EvtGenBase/EvtVector4C.hh"
00034 #include "EvtGenBase/EvtTensor4C.hh"
00035 using std::fstream;
00036
00037 EvtPi0Dalitz::~EvtPi0Dalitz() {}
00038
00039 void EvtPi0Dalitz::getName(std::string& model_name){
00040
00041 model_name="PI0_DALITZ";
00042
00043 }
00044
00045 EvtDecayBase* EvtPi0Dalitz::clone(){
00046
00047 return new EvtPi0Dalitz;
00048
00049 }
00050
00051
00052 void EvtPi0Dalitz::initProbMax(){
00053
00054 setProbMax(3.5);
00055
00056 }
00057
00058
00059 void EvtPi0Dalitz::init(){
00060
00061
00062 checkNArg(0);
00063 checkNDaug(3);
00064
00065
00066 checkSpinParent(EvtSpinType::SCALAR);
00067
00068 checkSpinDaughter(0,EvtSpinType::DIRAC);
00069 checkSpinDaughter(1,EvtSpinType::DIRAC);
00070 checkSpinDaughter(2,EvtSpinType::PHOTON);
00071
00072 }
00073
00074
00075 void EvtPi0Dalitz::decay( EvtParticle *p){
00076
00077 EvtParticle *ep, *em, *gamma;
00078 setWeight(p->initializePhaseSpace(getNDaug(),getDaugs(),0.00000002,0,1));
00079
00080 ep=p->getDaug(0);
00081 em=p->getDaug(1);
00082 gamma=p->getDaug(2);
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092 double m2=(ep->getP4()+em->getP4()).mass2();
00093 EvtVector4R q=ep->getP4()+em->getP4();
00094
00095
00096 EvtTensor4C w,v;
00097
00098 v=2.0*(gamma->getP4()*q)*directProd(q,gamma->getP4())
00099 - (gamma->getP4()*q)*(gamma->getP4()*q)*EvtTensor4C::g()
00100 -m2*directProd(gamma->getP4(),gamma->getP4());
00101
00102 w=4.0*( directProd(ep->getP4(),em->getP4()) + directProd(em->getP4(),ep->getP4())
00103 -EvtTensor4C::g()*(ep->getP4()*em->getP4()-ep->getP4().mass2()));
00104
00105 double prob=(real(cont(v,w)))/(m2*m2);
00106 prob *=(1.0/( (0.768*0.768-m2)*(0.768*0.768-m2)
00107 +0.768*0.768*0.151*0.151));
00108
00109
00110 setProb(prob);
00111
00112 return;
00113 }
00114
00115