EvtBTo3hCP.hh

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/>.     *
***********************************************************************/
 
#ifndef EVTBTO3HCP_HH
#define EVTBTO3HCP_HH
 
#include "EvtGenBase/EvtDecayAmp.hh"
#include "EvtGenBase/EvtVector4R.hh"
 
class EvtParticle;
 
class EvtBTo3hCP {
  public:
    void EvtKpipi( double alpha, double beta, int iset, EvtVector4R& p_K_plus,
                   EvtVector4R& p_pi_minus, EvtVector4R& p_gamma_1,
                   EvtVector4R& p_gamma_2, double& Real_B0, double& Imag_B0,
                   double& Real_B0bar, double& Imag_B0bar );
 
    void Evt3pi( double alpha, int iset, EvtVector4R& p_K_plus,
                 EvtVector4R& p_pi_minus, EvtVector4R& p_gamma_1,
                 EvtVector4R& p_gamma_2, double& Real_B0, double& Imag_B0,
                 double& Real_B0bar, double& Imag_B0bar );
 
    void Evt3piMPP( double alpha, int iset, EvtVector4R& p_p1,
                    EvtVector4R& p_p2, EvtVector4R& p_p3, double& Real_B0,
                    double& Imag_B0, double& Real_B0bar, double& Imag_B0bar );
 
    void Evt3piP00( double alpha, int iset, EvtVector4R& p_p1,
                    EvtVector4R& p_p1_gamma1, EvtVector4R& p_p1_gamma2,
                    EvtVector4R& p_p2_gamma1, EvtVector4R& p_p2_gamma2,
                    double& Real_B0, double& Imag_B0, double& Real_B0bar,
                    double& Imag_B0bar );
 
  private:
    void setConstants( double balpha, double bbeta );
    int computeKpipi( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3,
                      double& real_B0, double& imag_B0, double& real_B0bar,
                      double& imag_B0bar, int set );
    int compute3pi( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3,
                    double& real_B0, double& imag_B0, double& real_B0bar,
                    double& imag_B0bar, int set );
    int compute3piMPP( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3,
                       double& real_B0, double& imag_B0, double& real_B0bar,
                       double& imag_B0bar, int set );
    int compute3piP00( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3,
                       double& real_B0, double& imag_B0, double& real_B0bar,
                       double& imag_B0bar, int set );
 
    // Modes are : 0 = Kpipi, 1 = 3pi, 2 = MPP, 3 = P00
    void firstStep( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3, int mode );
    void generateSqMasses_Kpipi( double& m12, double& m13, double& m23,
                                 double MB2, double m1sq, double m2sq,
                                 double m3sq );
    void generateSqMasses_3pi( double& m12, double& m13, double& m23, double MB2,
                               double m1sq, double m2sq, double m3sq );
    void generateSqMasses_3piMPP( double& m12, double& m13, double& m23,
                                  double MB2, double m1sq, double m2sq,
                                  double m3sq );
    void generateSqMasses_3piP00( double& m12, double& m13, double& m23,
                                  double MB2, double m1sq, double m2sq,
                                  double m3sq );
 
    void rotation( EvtVector4R& p, int newRot );
    void gammaGamma( EvtVector4R& p, EvtVector4R& pgamma1, EvtVector4R& pgamma2 );
    EvtComplex BreitWigner( EvtVector4R& p1, EvtVector4R& p2, EvtVector4R& p3,
                            int& ierr, double Mass = 0, double Width = 0 );
    EvtComplex EvtRBW( double s, double Am2, double Gam, double Am2Min );
    EvtComplex EvtCRhoF_W( double s );
    EvtComplex EvtcBW_KS( double s, double Am2, double Gam );
    EvtComplex EvtcBW_GS( double s, double Am2, double Gam );
    double d( double AmRho2 );
    double k( double s );
    double Evtfs( double s, double AmRho2, double GamRho );
    double h( double s );
    double dh_ds( double s );
 
    EvtComplex m_Mat_S1, m_Mat_S2, m_Mat_S3, m_Mat_S4, m_Mat_S5, m_Nat_S1,
        m_Nat_S2, m_Nat_S3, m_Nat_S4, m_Nat_S5, m_MatKstarp, m_MatKstar0,
        m_MatKrho, m_NatKstarp, m_NatKstar0, m_NatKrho;
    double m_alphaCP = 1.365;
    double m_betaCP = 0.362;
    double m_MA2 = 27.927981186;
    double m_MB2 = 27.929242450;
    double m_MC2 = 28.153482608;
    double m_pi = 3.141592653;
    double m_Mass_rho = 0.770;
    double m_Gam_rho = 0.150;
    double m_M_B = 5.2794;
    double m_M_pip = 0.13957;
    double m_M_pim = 0.13957;
    double m_M_pi0 = 0.134976;
    double m_M_Kp = 0.49368;
    double m_Mass_Kstarp = 0.8916;
    double m_Mass_Kstar0 = 0.8961;
    double m_Gam_Kstarp = 0.0498;
    double m_Gam_Kstar0 = 0.0505;
 
    double m_rotMatrix[3][3];
    double m_factor_max = 1;
};
 
#endif