TEvtGen/EvtGenModels/EvtSLPole.cxx

   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/EvtSLPole.hh"
  29 #include "EvtGenModels/EvtSLPoleFF.hh"
  30 #include "EvtGenBase/EvtSemiLeptonicScalarAmp.hh"
  31 #include "EvtGenBase/EvtSemiLeptonicVectorAmp.hh"
  32 #include "EvtGenBase/EvtSemiLeptonicTensorAmp.hh"
  33 #include <string>
  34
  35 EvtSLPole::~EvtSLPole() {}
  36
  37 std::string EvtSLPole::getName(){
  38
  39   return "SLPOLE";
  40
  41 }
  42
  43
  44 EvtDecayBase* EvtSLPole::clone(){
  45
  46   return new EvtSLPole;
  47
  48 }
  49
  50 void EvtSLPole::decay( EvtParticle *p ){
  51
  52   p->initializePhaseSpace(getNDaug(),getDaugs(),_resetDaughterTree);
  53   calcamp->CalcAmp(p,_amp2,SLPoleffmodel);
  54   return;
  55 }
  56
  57 void EvtSLPole::initProbMax(){
  58
  59 EvtId parnum,mesnum,lnum,nunum;
  60
  61 parnum = getParentId();
  62 mesnum = getDaug(0);
  63 lnum = getDaug(1);
  64 nunum = getDaug(2);
  65
  66 double mymaxprob = calcamp->CalcMaxProb(parnum,mesnum,
  67                            lnum,nunum,SLPoleffmodel);
  68
  69 setProbMax(mymaxprob);
  70
  71 }
  72
  73
  74 void EvtSLPole::init(){
  75
  76   checkNDaug(3);
  77
  78   //We expect the parent to be a scalar
  79   //and the daughters to be X lepton neutrino
  80
  81   checkSpinParent(EvtSpinType::SCALAR);
  82   checkSpinDaughter(1,EvtSpinType::DIRAC);
  83   checkSpinDaughter(2,EvtSpinType::NEUTRINO);
  84
  85   EvtSpinType::spintype mesontype=EvtPDL::getSpinType(getDaug(0));
  86
  87   SLPoleffmodel = new EvtSLPoleFF(getNArg(),getArgs());
  88
  89   if ( mesontype==EvtSpinType::SCALAR ) {
  90     calcamp = new EvtSemiLeptonicScalarAmp;
  91   }
  92   if ( mesontype==EvtSpinType::VECTOR ) {
  93     calcamp = new EvtSemiLeptonicVectorAmp;
  94   }
  95   if ( mesontype==EvtSpinType::TENSOR ) {
  96     calcamp = new EvtSemiLeptonicTensorAmp;
  97   }
  98
  99   _resetDaughterTree=false;
 100   if ( getArgStr(getNArg()-1) == "true") _resetDaughterTree=true;
 101
 102 }
 103