Line data    Source code

       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/EvtBcBsStarNPi.hh
      10             : //
      11             : // Description: Decay model for Bc -> Bs* + 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             : 
      21             : #include "EvtGenModels/EvtBcBsStarNPi.hh"
      22             : #include "EvtGenBase/EvtSpinType.hh"
      23             : 
      24           0 : EvtBcBsStarNPi::EvtBcBsStarNPi() {
      25             : 
      26           0 :   _beta=-0.108; _mRho=0.775; _gammaRho=0.149;
      27           0 :   _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4;
      28             : 
      29           0 :   FA0_N=8.1; FA0_c1=0.30; FA0_c2=0.069;
      30           0 :   FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0;
      31           0 :   FAp_N=0.15; FAp_c1=0.30; FAp_c2=0.069;
      32           0 :   FV_N= 1.08; FV_c1=0.30; FV_c2=0.069; 
      33             : 
      34           0 : }
      35             : 
      36           0 : EvtBcBsStarNPi::~EvtBcBsStarNPi() {
      37             : 
      38           0 : }
      39             : 
      40             : std::string EvtBcBsStarNPi::getName() {
      41             : 
      42           0 :   return "BC_BSSTAR_NPI";
      43             : 
      44             : }
      45             : 
      46             : EvtDecayBase* EvtBcBsStarNPi::clone() {
      47             : 
      48           0 :   return new EvtBcBsStarNPi;
      49             : 
      50           0 : }
      51             : 
      52             : void EvtBcBsStarNPi::init() {
      53             : 
      54           0 :   checkNArg(0);
      55             : 
      56             :   // check spins
      57           0 :   checkSpinParent(EvtSpinType::SCALAR);
      58           0 :   checkSpinDaughter(0,EvtSpinType::VECTOR);
      59             :   // the others are scalar
      60           0 :   for (int i=1; i<=(getNDaug()-1);i++) {
      61           0 :     checkSpinDaughter(i,EvtSpinType::SCALAR);
      62             :   }
      63             : 
      64           0 : }
      65             : 
      66             : void EvtBcBsStarNPi::initProbMax() {
      67             : 
      68           0 :   if ( getNDaug() == 2 ) {
      69           0 :     setProbMax(100.);
      70           0 :   } else if( getNDaug() == 3 ) {
      71           0 :     setProbMax(40000.);
      72           0 :   } else if( getNDaug() == 4 ) {
      73           0 :     setProbMax(620.);  // checked, 30k events
      74           0 :   }
      75             : 
      76           0 : }