TEvtGen/EvtGenModels/EvtBcPsiNPi.cpp

   1 //--------------------------------------------------------------------------
   2 //
   3 // Environment:
   4 //      This software is part of the EvtGen package. If you use all or part
   5 //      of it, please give an appropriate acknowledgement.
   6 //
   7 // Copyright Information: See EvtGen/COPYRIGHT
   8 //
   9 // Module: EvtGenModels/EvtBcBsNPi.hh
  10 //
  11 // Description: Decay model for Bc -> J/psi + npi
  12 //
  13 // Modification history:
  14 //
  15 //    A.Berezhnoy, A.Likhoded, A.Luchinsky  April 2011   Module created
  16 //
  17 //------------------------------------------------------------------------
  18
  19 #include "EvtGenBase/EvtPatches.hh"
  20 #include "EvtGenModels/EvtBcPsiNPi.hh"
  21 #include "EvtGenBase/EvtSpinType.hh"
  22
  23 EvtBcPsiNPi::EvtBcPsiNPi() {
  24
  25   _beta=-0.108; _mRho=0.775; _gammaRho=0.149;
  26   _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4;
  27
  28   FA0_N=5.9; FA0_c1= 0.049; FA0_c2= 0.0015;
  29   FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0;
  30   FAp_N=-0.074; FAp_c1= 0.049; FAp_c2= 0.0015;
  31   FV_N=0.11; FV_c1= 0.049; FV_c2= 0.0015;
  32
  33 }
  34
  35 EvtBcPsiNPi::~EvtBcPsiNPi() {
  36
  37 }
  38
  39 std::string EvtBcPsiNPi::getName() {
  40
  41   return "BC_PSI_NPI";
  42
  43 }
  44
  45 EvtDecayBase* EvtBcPsiNPi::clone() {
  46
  47   return new EvtBcPsiNPi;
  48
  49 }
  50
  51 void EvtBcPsiNPi::init() {
  52
  53   checkNArg(0);
  54
  55   // check spins
  56   checkSpinParent(EvtSpinType::SCALAR);
  57   checkSpinDaughter(0,EvtSpinType::VECTOR);
  58   // the others are scalar
  59   for (int i=1; i<=(getNDaug()-1);i++) {
  60     checkSpinDaughter(i,EvtSpinType::SCALAR);
  61   }
  62
  63 }
  64
  65 void EvtBcPsiNPi::initProbMax() {
  66
  67   setProbMax(100.);
  68   if( getNDaug() == 2 ) {
  69     setProbMax(330.);
  70   } else if( getNDaug() == 3 ) {
  71     setProbMax(11000.); // checked with 30k events
  72   } else if( getNDaug() == 4 ) {
  73     setProbMax(36000.);
  74   }
  75
  76 }