include/EvtGenModels/EvtKStopizmumu.hh

0001
0002 /***********************************************************************
0003 * Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
0004 *                                                                      *
0005 * This file is part of EvtGen.                                         *
0006 *                                                                      *
0007 * EvtGen is free software: you can redistribute it and/or modify       *
0008 * it under the terms of the GNU General Public License as published by *
0009 * the Free Software Foundation, either version 3 of the License, or    *
0010 * (at your option) any later version.                                  *
0011 *                                                                      *
0012 * EvtGen is distributed in the hope that it will be useful,            *
0013 * but WITHOUT ANY WARRANTY; without even the implied warranty of       *
0014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
0015 * GNU General Public License for more details.                         *
0016 *                                                                      *
0017 * You should have received a copy of the GNU General Public License    *
0018 * along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
0019 ***********************************************************************/
0020
0021 #ifndef EVTKSTOPIZMUMU_HH
0022 #define EVTKSTOPIZMUMU_HH
0023
0024 #include "EvtGenBase/EvtDecayAmp.hh"
0025
0026 #include <string>
0027
0028 class EvtParticle;
0029
0030 // Description: Routine to implement KS -> pi0 mu mu; see JHEP08(1998)004.
0031
0032 class EvtKStopizmumu : public EvtDecayAmp {
0033   public:
0034     std::string getName() override { return "KS_PI0MUMU"; }
0035
0036     EvtDecayBase* clone() override { return new EvtKStopizmumu; }
0037
0038     void init() override;
0039
0040     void initProbMax() override { setProbMax( 1.0e-10 ); }
0041
0042     void decay( EvtParticle* p ) override;
0043
0044     double F_z( const double& z, const double& rvsq );
0045     EvtComplex G_z( const double& z );
0046     double Wpol_z( const double& z, const double& as, const double& bs );
0047     EvtComplex chi_z( const double& z, const double& rpisq );
0048     EvtComplex Wpipi_z( const double& z, const double& alpha_s,
0049                         const double& beta_s, const double& rvsq,
0050                         const double& rpisq, const double& z0 );
0051 };
0052
0053 #endif    //EVTKTOPIZMUMU_HH