cepgen is hosted by Hepforge, IPPP Durham
CepGen 1.2.5
Central exclusive processes event generator
Loading...
Searching...
No Matches
AlphaSWebber.cpp
Go to the documentation of this file.
1/*
2 * CepGen: a central exclusive processes event generator
3 * Copyright (C) 2024 Laurent Forthomme
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19#include <cmath>
20
24
25namespace cepgen {
28 class AlphaSWebber : public Coupling {
29 public:
30 explicit AlphaSWebber(const ParametersList& params)
31 : Coupling(params),
32 nc_(steer<int>("Nc")),
33 nf_(steer<int>("nf")),
34 lambda_(steer<double>("Lambda")),
35 beta0_((11. * nc_ - 2. * nf_) / 3.),
36 prefac_(4. * M_PI / beta0_) {
37 CG_INFO("AlphaSWebber:init") << "Webber et al. alpha(S) evolution algorithm initialised with parameters:\n\t"
38 << "Nc: " << nc_ << ", nf: " << nf_ << " -> beta0: " << beta0_
39 << ", Lambda: " << lambda_ << ".";
40 }
41
43 auto desc = Coupling::description();
44 desc.setDescription("Webber alpha(S) evolution algorithm");
45 desc.add<int>("Nc", 3).setDescription("number of colours considered");
46 desc.add<int>("nf", 3).setDescription("number of fermion flavours considered");
47 desc.add<double>("Lambda", 0.25).setDescription("evolution scale (in GeV)");
48 return desc;
49 }
50
51 double operator()(double q) const override {
52 const auto mun = q * q / lambda_ / lambda_;
53 return prefac_ * (1. / std::log(mun) + 125. * (1. + 4. * mun) / (1. - mun) / std::pow(4. + mun, 4));
54 }
55
56 private:
57 const int nc_, nf_;
58 const double lambda_, beta0_;
59 const double prefac_;
60 };
61} // namespace cepgen
62
63REGISTER_ALPHAS_MODULE("webber", AlphaSWebber);
#define REGISTER_ALPHAS_MODULE(name, obj)
Add a strong coupling evolution algorithm.
#define CG_INFO(mod)
Definition Message.h:216
Simple parameterisation of the QCD running coupling at low scales .
double operator()(double q) const override
Compute for a given .
static ParametersDescription description()
AlphaSWebber(const ParametersList &params)
A generic evaluation algorithm.
Definition Coupling.h:26
A description object for parameters collection.
static ParametersDescription description()
Description of all object parameters.
Definition Steerable.cpp:42
T steer(const std::string &key) const
Retrieve a parameters as previously steered.
Definition Steerable.h:39
Common namespace for this Monte Carlo generator.