EvtDToKSpipi0pi0.h

// Model: EvtDToKSpipi0pi0
// This file is an amplitude model for D+ -> K_S0 pi- pi0 pi0.
// The model is from the BESIII Collaboration in JHEP09 (2023) 077. DOI:  https://doi.org/10.1007/JHEP09(2023)077
//
// Permission to include these files in basf2 was generously granted by the BESIII Collaboration.
//
// Please cite the original reference for any public/published results where this model was used.
 
#pragma once
 
#include <EvtGenBase/EvtDecayProb.hh>
 
namespace Belle2 {
 
  class EvtDToKSpipi0pi0: public EvtDecayProb  {
 
  public:
    EvtDToKSpipi0pi0() {}
    virtual ~EvtDToKSpipi0pi0();
    std::string getName();
    EvtDecayBase* clone();
    void init();
    void initProbMax();
    void decay(EvtParticle* p);
 
  private:
    double calPDF(const double Ks[], const double Pip[], const double Pi01[], const double Pi02[], double& Result);

    void Com_Multi(const double a1[2], const double a2[2], double res[2]);
    void Com_Divide(const double a1[2], const double a2[2], double res[2]);
    double SCADot(const double a1[4], const double a2[4]);

    double Barrier(const double mass2, const int l, const double sa, const double sb, const double sc, const double r2);
    void calt1(const double daug1[4], const double daug2[4], double t1[4]);
    void calt2(const double daug1[4], const double daug2[4], double t2[4][4]);
    double wid(const double mass2, const double mass, const double sa, const double sb, const double sc, const double r2, const int l);
    double widl1(const double mass2, const double mass, const double sa, const double sb, const double sc, const double r2);

    void propagatorRBW(const double mass2, const double mass, const double width, const double sa, const double sb, const double sc,
                       const double r2, const int l, double prop[2]);
    void propagatorRBWl1(const double mass2, const double mass, const double width, const double sa, const double sb, const double sc,
                         const double r2, double prop[2]);
    void propagatorRBW_a1(const double mass2, const double mass, const double width, const double sa, const double sb, const double sc,
                          const double r2, const int l, double prop[2]);
    void propagatorGS(const double mass2, const double mass, const double width, const double sa, const double sb, const double sc,
                      const double r2, double prop[2]);
    void propagatorsigma500(const double sa, const double sb, const double sc, double prop[2]);

    void rhoab(const double sa, const double sb, const double sc, double res[2]);
    void rho4Pi(const double sa, double res[2]);
    void KPiSLASS(const double sa, const double sb, const double sc, double prop[2]);

    int mother_c;

    static constexpr int STDim = 4;
    double G[STDim][STDim], E[STDim][STDim][STDim][STDim];

    static constexpr int numOfParameters = 24;
    double rho[numOfParameters];
    double phi[numOfParameters];

    static constexpr int topologies = 20;
    double modetype[topologies];

    double mrho, Grho, mKstr0, GKstr0, msigma, Gsigma;
    double ma1, Ga1, mK1270, GK1270, mK1400, GK1400;

    double mD;
    double math_pi;
  };

} // Belle 2 Namespace