include/EvtGenModels/EvtVubAC.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 EVTVUBAC_HH
0022 #define EVTVUBAC_HH
0023
0024 #include "EvtGenBase/EvtDecayIncoherent.hh"
0025
0026 #include <vector>
0027
0028 class EvtParticle;
0029
0030 // Analytic Coupling Model (based on hep-ph/0608047 by Aglietti, Ferrera and Ricciardi)
0031
0032 class EvtVubAC : public EvtDecayIncoherent {
0033   public:
0034     std::string getName() override;
0035
0036     EvtDecayBase* clone() override;
0037
0038     void initProbMax() override;
0039
0040     void init() override;
0041
0042     void decay( EvtParticle* Bmeson ) override;
0043
0044   private:
0045     // Input parameters
0046     double mB;
0047
0048     double alphaSmZ;
0049     double alphaSmB;
0050     double c;
0051     double q;
0052     double k;
0053
0054     double CF;
0055     double CA;
0056
0057     double beta0;
0058
0059     std::vector<double> gvars;
0060
0061     double rate( double u, double w, double xb );
0062     double wreg( double w );
0063     double alphaS( double Q );
0064     double PolyLog( double v, double z );
0065     double ureg( double u );
0066     double ularge( double u );
0067     double Coeff( double u, double w, double xb );
0068     double Coeff1( double w, double xb );
0069     double Coeff0( double w, double xb );
0070     double Sigma( double x1, double x2 );
0071     double max( double ub, double lb );
0072     double d1( double u, double w, double xb );
0073     double d( double u, double w, double xb );
0074     double f( double w );
0075     double Lambda2( double x, double alphaSmZ );
0076     int Bisect( double x1, double x2, double precision, double& root,
0077                 const double alphaSmZ );
0078     double FindRoot( const double alphaSmZ );
0079 };
0080
0081 #endif