TEvtGen/EvtGenModels/EvtMelikhov.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: EvtMelikhov.cc
  12 //
  13 // Description: Routine to implement semileptonic B->D*lnu decays according
  14 //              to the model HQET
  15 //
  16 // Modification history:
  17 //
  18 //    DJL     April 20, 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/EvtMelikhov.hh"
  29 #include "EvtGenModels/EvtMelikhovFF.hh"
  30 #include "EvtGenBase/EvtSemiLeptonicVectorAmp.hh"
  31 #include <string>
  32
  33 EvtMelikhov::~EvtMelikhov() {}
  34
  35 std::string EvtMelikhov::getName(){
  36
  37   return "MELIKHOV";
  38 }
  39
  40
  41 EvtDecayBase* EvtMelikhov::clone(){
  42
  43   return new EvtMelikhov;
  44
  45 }
  46
  47 void EvtMelikhov::decay( EvtParticle *p ){
  48
  49   p->initializePhaseSpace(getNDaug(),getDaugs());
  50   calcamp->CalcAmp(p,_amp2,Melikhovffmodel);
  51 }
  52
  53
  54 void EvtMelikhov::init(){
  55
  56   checkNArg(1);
  57   checkNDaug(3);
  58
  59   //We expect the parent to be a scalar
  60   //and the daughters to be X lepton neutrino
  61
  62   checkSpinParent(EvtSpinType::SCALAR);
  63
  64   checkSpinDaughter(0,EvtSpinType::VECTOR);
  65   checkSpinDaughter(1,EvtSpinType::DIRAC);
  66   checkSpinDaughter(2,EvtSpinType::NEUTRINO);
  67
  68
  69   Melikhovffmodel = new EvtMelikhovFF(getArg(0));
  70   calcamp = new EvtSemiLeptonicVectorAmp;
  71
  72 }
  73