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 -> 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           0 : EvtBcPsiNPi::EvtBcPsiNPi() {
      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             : 
      28           0 :   FA0_N=5.9; FA0_c1= 0.049; FA0_c2= 0.0015;
      29           0 :   FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0;
      30           0 :   FAp_N=-0.074; FAp_c1= 0.049; FAp_c2= 0.0015;
      31           0 :   FV_N=0.11; FV_c1= 0.049; FV_c2= 0.0015; 
      32             : 
      33           0 : }
      34             : 
      35           0 : EvtBcPsiNPi::~EvtBcPsiNPi() {
      36             : 
      37           0 : }
      38             : 
      39             : std::string EvtBcPsiNPi::getName() {
      40             : 
      41           0 :   return "BC_PSI_NPI";
      42             : 
      43             : }
      44             : 
      45             : EvtDecayBase* EvtBcPsiNPi::clone() {
      46             : 
      47           0 :   return new EvtBcPsiNPi;
      48             : 
      49           0 : }
      50             : 
      51             : void EvtBcPsiNPi::init() {
      52             : 
      53           0 :   checkNArg(0);
      54             : 
      55             :   // check spins
      56           0 :   checkSpinParent(EvtSpinType::SCALAR);
      57           0 :   checkSpinDaughter(0,EvtSpinType::VECTOR);
      58             :   // the others are scalar
      59           0 :   for (int i=1; i<=(getNDaug()-1);i++) {
      60           0 :     checkSpinDaughter(i,EvtSpinType::SCALAR);
      61             :   }
      62             : 
      63           0 : }
      64             : 
      65             : void EvtBcPsiNPi::initProbMax() {
      66             : 
      67           0 :   setProbMax(100.);
      68           0 :   if( getNDaug() == 2 ) {
      69           0 :     setProbMax(330.);
      70           0 :   } else if( getNDaug() == 3 ) {
      71           0 :     setProbMax(11000.); // checked with 30k events
      72           0 :   } else if( getNDaug() == 4 ) {
      73           0 :     setProbMax(36000.);
      74           0 :   }
      75             : 
      76           0 : }