EvtBcPsiNPi.cpp

Go to the documentation of this file.

 
/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen is free software: you can redistribute it and/or modify       *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or    *
* (at your option) any later version.                                  *
*                                                                      *
* EvtGen is distributed in the hope that it will be useful,            *
* but WITHOUT ANY WARRANTY; without even the implied warranty of       *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
* GNU General Public License for more details.                         *
*                                                                      *
* You should have received a copy of the GNU General Public License    *
* along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/
 
#include "EvtGenModels/EvtBcPsiNPi.hh"
 
#include "EvtGenBase/EvtSpinType.hh"
 
EvtBcPsiNPi::EvtBcPsiNPi()
{
    m_beta = -0.108;
    m_mRho = 0.775;
    m_gammaRho = 0.149;
    m_mRhopr = 1.364;
    m_gammaRhopr = 0.400;
    m_mA1 = 1.23;
    m_gammaA1 = 0.4;
 
    m_FA0_N = 5.9;
    m_FA0_c1 = 0.049;
    m_FA0_c2 = 0.0015;
    m_FAm_N = 0.0;
    m_FAm_c1 = 0.0;
    m_FAm_c2 = 0.0;
    m_FAp_N = -0.074;
    m_FAp_c1 = 0.049;
    m_FAp_c2 = 0.0015;
    m_FV_N = 0.11;
    m_FV_c1 = 0.049;
    m_FV_c2 = 0.0015;
}
 
std::string EvtBcPsiNPi::getName() const
{
    return "BC_PSI_NPI";
}
 
EvtBcPsiNPi* EvtBcPsiNPi::clone() const
{
    return new EvtBcPsiNPi;
}
 
void EvtBcPsiNPi::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 EvtBcPsiNPi::initProbMax()
{
    setProbMax( 100. );
    if ( getNDaug() == 2 ) {
        setProbMax( 330. );
    } else if ( getNDaug() == 3 ) {
        setProbMax( 11000. );    // checked with 30k events
    } else if ( getNDaug() == 4 ) {
        setProbMax( 36000. );
    }
}