//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package. If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//
// Module: EvtGenModels/EvtBcBsStarNPi.hh
//
// Description: Decay model for Bc -> Bs* + npi
//
// Modification history:
//
//    A.Berezhnoy, A.Likhoded, A.Luchinsky  April 2011   Module created
//
//------------------------------------------------------------------------

#include "EvtGenBase/EvtPatches.hh"

#include "EvtGenModels/EvtBcBsStarNPi.hh"
#include "EvtGenBase/EvtSpinType.hh"

EvtBcBsStarNPi::EvtBcBsStarNPi() {

  _beta=-0.108; _mRho=0.775; _gammaRho=0.149;
  _mRhopr=1.364; _gammaRhopr=0.400; _mA1=1.23; _gammaA1=0.4;

  FA0_N=8.1; FA0_c1=0.30; FA0_c2=0.069;
  FAm_N=0.0; FAm_c1=0.0; FAm_c2=0.0;
  FAp_N=0.15; FAp_c1=0.30; FAp_c2=0.069;
  FV_N= 1.08; FV_c1=0.30; FV_c2=0.069; 

}

EvtBcBsStarNPi::~EvtBcBsStarNPi() {

}

std::string EvtBcBsStarNPi::getName() {

  return "BC_BSSTAR_NPI";

}

EvtDecayBase* EvtBcBsStarNPi::clone() {

  return new EvtBcBsStarNPi;

}

void EvtBcBsStarNPi::init() {

  checkNArg(0);

  // check spins
  checkSpinParent(EvtSpinType::SCALAR);
  checkSpinDaughter(0,EvtSpinType::VECTOR);
  // the others are scalar
  for (int i=1; i<=(getNDaug()-1);i++) {
    checkSpinDaughter(i,EvtSpinType::SCALAR);
  }

}

void EvtBcBsStarNPi::initProbMax() {

  if ( getNDaug() == 2 ) {
    setProbMax(100.);
  } else if( getNDaug() == 3 ) {
    setProbMax(40000.);
  } else if( getNDaug() == 4 ) {
    setProbMax(620.);  // checked, 30k events
  }

}