cepgen is hosted by Hepforge, IPPP Durham
CepGen 1.2.5
Central exclusive processes event generator
Loading...
Searching...
No Matches
FunctionIntegrand.cpp
Go to the documentation of this file.
1/*
2 * CepGen: a central exclusive processes event generator
3 * Copyright (C) 2013-2022 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
21
22namespace cepgen {
23 FunctionIntegrand::FunctionIntegrand(size_t ndim, const std::function<double(const std::vector<double>&)>& func)
24 : function_(func), ndim_(ndim) {}
25
26 double FunctionIntegrand::eval(const std::vector<double>& x) {
27 if (x.size() != size())
28 throw CG_FATAL("FunctionIntegrand:eval")
29 << "Invalid coordinates multiplicity: expected(" << size() << ") != received(" << x.size() << ")!";
30
31 //--- calculate weight for the phase space point to probe
32 double weight = function_(x);
33
34 //--- a bit of useful debugging
35 CG_DEBUG_LOOP("FunctionIntegrand:eval")
36 << "f value for dim-" << x.size() << " point " << x << ": " << weight << ".";
37
38 return weight;
39 }
40} // namespace cepgen
#define CG_FATAL(mod)
Definition Exception.h:61
#define CG_DEBUG_LOOP(mod)
Definition Message.h:224
FunctionIntegrand(size_t, const std::function< double(const std::vector< double > &)> &)
double eval(const std::vector< double > &) override
Compute the integrand for a given coordinates set.
size_t size() const override
Phase space dimension.
Common namespace for this Monte Carlo generator.