TEvtGen/EvtGenModels/EvtKKLambdaC.cpp

   1 //--------------------------------------------------------------------------
   2 //
   3 // Environment:
   4 //      This software is part of the EvtGen package developed jointly
   5 //      for the BaBar and CLEO collaborations.  If you use all or part
   6 //      of it, please give an appropriate acknowledgement.
   7 //
   8 // Copyright Information: See EvtGen/COPYRIGHT
   9 //      Copyright (C) 1998      Caltech, UCSB
  10 //
  11 // Module: EvtSLPole.cc
  12 //
  13 // Description: Routine to implement semileptonic decays according
  14 //              to light cone sum rules
  15 //
  16 // Modification history:
  17 //
  18 //    DJL       April 23, 1998       Module created
  19 //
  20 //------------------------------------------------------------------------
  21 //
  22 #include "EvtGenBase/EvtPatches.hh"
  23 #include <stdlib.h>
  24 #include "EvtGenBase/EvtParticle.hh"
  25 #include "EvtGenBase/EvtGenKine.hh"
  26 #include "EvtGenBase/EvtPDL.hh"
  27 #include "EvtGenBase/EvtReport.hh"
  28 #include "EvtGenModels/EvtKKLambdaC.hh"
  29 #include "EvtGenModels/EvtKKLambdaCFF.hh"
  30 #include "EvtGenBase/EvtSemiLeptonicBaryonAmp.hh"
  31 #include <string>
  32
  33 EvtKKLambdaC::~EvtKKLambdaC() {}
  34
  35 std::string EvtKKLambdaC::getName(){
  36
  37   return "KK_LAMBDAC_SL";
  38
  39 }
  40
  41
  42 EvtDecayBase* EvtKKLambdaC::clone(){
  43
  44   return new EvtKKLambdaC;
  45
  46 }
  47
  48 void EvtKKLambdaC::decay( EvtParticle *p ){
  49
  50   p->initializePhaseSpace(getNDaug(),getDaugs());
  51
  52   _calcamp->CalcAmp(p,_amp2,_ffmodel);
  53   return;
  54 }
  55
  56 void EvtKKLambdaC::initProbMax(){
  57
  58 EvtId parnum,mesnum,lnum,nunum;
  59
  60 parnum = getParentId();
  61 mesnum = getDaug(0);
  62 lnum = getDaug(1);
  63 nunum = getDaug(2);
  64
  65 //double mymaxprob = _calcamp->CalcMaxProb(parnum,mesnum,
  66 //                           lnum,nunum,_ffmodel);
  67  double mymaxprob=100.;
  68 setProbMax(mymaxprob);
  69
  70 }
  71
  72
  73 void EvtKKLambdaC::init(){
  74
  75   checkNDaug(3);
  76
  77   //We expect the parent to be a dirac
  78   //and the daughters to be dirac lepton neutrino
  79
  80   checkSpinParent(EvtSpinType::DIRAC);
  81   checkSpinDaughter(0,EvtSpinType::DIRAC);
  82   checkSpinDaughter(1,EvtSpinType::DIRAC);
  83   checkSpinDaughter(2,EvtSpinType::NEUTRINO);
  84
  85   _ffmodel = new EvtKKLambdaCFF(getNArg(),getArgs());
  86
  87   _calcamp = new EvtSemiLeptonicBaryonAmp;
  88
  89 }
  90