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
19
#include "
CepGen/Core/Exception.h
"
20
#include "
CepGen/Integration/FunctionIntegrand.h
"
21
22
namespace
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
Exception.h
CG_FATAL
#define CG_FATAL(mod)
Definition
Exception.h:61
FunctionIntegrand.h
CG_DEBUG_LOOP
#define CG_DEBUG_LOOP(mod)
Definition
Message.h:224
cepgen::FunctionIntegrand::FunctionIntegrand
FunctionIntegrand(size_t, const std::function< double(const std::vector< double > &)> &)
Definition
FunctionIntegrand.cpp:23
cepgen::FunctionIntegrand::eval
double eval(const std::vector< double > &) override
Compute the integrand for a given coordinates set.
Definition
FunctionIntegrand.cpp:26
cepgen::FunctionIntegrand::size
size_t size() const override
Phase space dimension.
Definition
FunctionIntegrand.h:33
cepgen
Common namespace for this Monte Carlo generator.
Definition
CommandLineHandler.cpp:36
CepGen
Integration
FunctionIntegrand.cpp
Generated on Mon Jul 29 2024 for CepGen by
1.9.7