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/EvtBcBsNPi.hh
      10             : //
      11             : // Description: Decay model for Bc -> Bs + npi
      12             : //
      13             : // Modification history:
      14             : //
      15             : //    A.Berezhnoy, A.Likhoded, A.Luchinsky  July 2011   Module created
      16             : //
      17             : //------------------------------------------------------------------------
      18             : 
      19             : #include "EvtGenBase/EvtPatches.hh"
      20             : #include "EvtGenBase/EvtSpinType.hh"
      21             : #include "EvtGenModels/EvtBcBsNPi.hh"
      22             : 
      23           0 : EvtBcBsNPi::EvtBcBsNPi() {
      24             : 
      25           0 :   _beta=-0.108; _mRho=0.775; _gammaRho=0.149;
      26           0 :   _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4;
      27             :   //            Fp_N=1.3; Fp_c1=0.30; Fp_c2=0.069;
      28           0 :   Fp_N=3*1.3; Fp_c1=0.30; Fp_c2=0.069;
      29           0 :   Fm_N=0.0; Fm_c1=0.0; Fm_c2=0.0;
      30             : 
      31           0 : }
      32             : 
      33           0 : EvtBcBsNPi::~EvtBcBsNPi() {
      34           0 : }
      35             : 
      36             : std::string EvtBcBsNPi::getName() {
      37             : 
      38           0 :   return "BC_BS_NPI";
      39             : 
      40             : }
      41             : 
      42             : EvtDecayBase* EvtBcBsNPi::clone() {
      43             : 
      44           0 :   return new EvtBcBsNPi;
      45             : 
      46           0 : }
      47             : 
      48             : void EvtBcBsNPi::init() {
      49             : 
      50           0 :   checkNArg(0);
      51             : 
      52             :   // check spins
      53           0 :   checkSpinParent(EvtSpinType::SCALAR);
      54           0 :   checkSpinDaughter(0,EvtSpinType::SCALAR);
      55             :   // the others are scalar
      56           0 :   for (int i=1; i<=(getNDaug()-1);i++) {
      57           0 :     checkSpinDaughter(i,EvtSpinType::SCALAR);
      58             :   }
      59             : 
      60           0 : }
      61             : 
      62             : void EvtBcBsNPi::initProbMax() {
      63             : 
      64           0 :   if ( getNDaug() == 2 ) {
      65           0 :     setProbMax(250.);
      66           0 :   } else if ( getNDaug() == 3 ) {
      67           0 :     setProbMax(25000.);// checked at 30k events
      68           0 :   } else if( getNDaug() == 4 ) {
      69           0 :     setProbMax(45000.); // checked at 30k events
      70           0 :   }
      71             : 
      72           0 : }