00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 #include "EvtGenBase/EvtPatches.hh"
00023 #include <stdlib.h>
00024 #include "EvtGenBase/EvtParticle.hh"
00025 #include "EvtGenBase/EvtGenKine.hh"
00026 #include "EvtGenBase/EvtPDL.hh"
00027 #include "EvtGenBase/EvtReport.hh"
00028 #include "EvtGenModels/EvtSLPole.hh"
00029 #include "EvtGenModels/EvtSLPoleFF.hh"
00030 #include "EvtGenBase/EvtSemiLeptonicScalarAmp.hh"
00031 #include "EvtGenBase/EvtSemiLeptonicVectorAmp.hh"
00032 #include "EvtGenBase/EvtSemiLeptonicTensorAmp.hh"
00033 #include <string>
00034
00035 EvtSLPole::~EvtSLPole() {}
00036
00037 void EvtSLPole::getName(std::string& model_name){
00038
00039 model_name="SLPOLE";
00040
00041 }
00042
00043
00044 EvtDecayBase* EvtSLPole::clone(){
00045
00046 return new EvtSLPole;
00047
00048 }
00049
00050 void EvtSLPole::decay( EvtParticle *p ){
00051
00052 p->initializePhaseSpace(getNDaug(),getDaugs());
00053
00054 calcamp->CalcAmp(p,_amp2,SLPoleffmodel);
00055 return;
00056 }
00057
00058 void EvtSLPole::initProbMax(){
00059
00060 EvtId parnum,mesnum,lnum,nunum;
00061
00062 parnum = getParentId();
00063 mesnum = getDaug(0);
00064 lnum = getDaug(1);
00065 nunum = getDaug(2);
00066
00067 double mymaxprob = calcamp->CalcMaxProb(parnum,mesnum,
00068 lnum,nunum,SLPoleffmodel);
00069
00070 setProbMax(mymaxprob);
00071
00072 }
00073
00074
00075 void EvtSLPole::init(){
00076
00077 checkNDaug(3);
00078
00079
00080
00081
00082 checkSpinParent(EvtSpinType::SCALAR);
00083 checkSpinDaughter(1,EvtSpinType::DIRAC);
00084 checkSpinDaughter(2,EvtSpinType::NEUTRINO);
00085
00086 EvtSpinType::spintype mesontype=EvtPDL::getSpinType(getDaug(0));
00087
00088 SLPoleffmodel = new EvtSLPoleFF(getNArg(),getArgs());
00089
00090 if ( mesontype==EvtSpinType::SCALAR ) {
00091 calcamp = new EvtSemiLeptonicScalarAmp;
00092 }
00093 if ( mesontype==EvtSpinType::VECTOR ) {
00094 calcamp = new EvtSemiLeptonicVectorAmp;
00095 }
00096 if ( mesontype==EvtSpinType::TENSOR ) {
00097 calcamp = new EvtSemiLeptonicTensorAmp;
00098 }
00099
00100 }
00101