APFEL 4.8.0
A PDF evolution library in C++
Loading...
Searching...
No Matches
tmdbuilder.h
Go to the documentation of this file.
1//
2// APFEL++ 2017
3//
4// Author: Valerio Bertone: valerio.bertone@cern.ch
5//
6
7#pragma once
8
9#include "apfel/grid.h"
10#include "apfel/operator.h"
11#include "apfel/set.h"
12#include "apfel/dglapbuilder.h"
13#include "apfel/tabulateobject.h"
14#include "apfel/constants.h"
15
16namespace apfel
17{
21 enum JetAlgorithm: int {CONE = 0, KT = 1};
22
30 struct TmdObjects
31 {
32 double Threshold;
33 std::map<int, double> Beta;
34 std::map<int, double> GammaFq;
35 std::map<int, double> GammaFg;
36 std::map<int, double> GammaK;
37 std::map<int, std::vector<double>> KCS;
38 std::map<int, std::vector<Set<Operator>>> MatchingFunctionsPDFs;
39 std::map<int, std::vector<Set<Operator>>> MatchingFunctionsFFs;
40 std::map<std::string, std::map<int, double>> HardFactors;
41 };
42
49
58 std::map<int, TmdObjects> InitializeTmdObjects(Grid const& g,
59 std::vector<double> const& Thresholds,
60 double const& IntEps = 1e-5);
61
74 std::map<int, TmdObjects> InitializeTmdObjectsDYResScheme(Grid const& g,
75 std::vector<double> const& Thresholds,
76 double const& IntEps = 1e-5);
77
89 std::map<int, TmdObjects> InitializeTmdObjectsBM(Grid const& g,
90 std::vector<double> const& Thresholds,
91 double const& IntEps = 1e-5);
92
105 std::map<int, TmdObjects> InitializeTmdObjectsSivers(Grid const& g,
106 std::vector<double> const& Thresholds,
107 double const& IntEps = 1e-5);
108
118 std::map<int, TmdObjects> InitializeTmdObjectsg1(Grid const& g,
119 std::vector<double> const& Thresholds,
120 double const& IntEps = 1e-5);
122
132
145 std::function<Set<Distribution>(double const&, double const&, double const&)> BuildTmdPDFs(std::map<int, TmdObjects> const& TmdObj,
146 std::function<Set<Distribution>(double const&)> const& CollPDFs,
147 std::function<double(double const&)> const& Alphas,
148 int const& PerturbativeOrder,
149 double const& Ci = 1,
150 double const& IntEps = 1e-7);
151
165 std::function<Set<Distribution>(double const&, double const&, double const&)> BuildTmdFFs(std::map<int, TmdObjects> const& TmdObj,
166 std::function<Set<Distribution>(double const&)> const& CollFFs,
167 std::function<double(double const&)> const& Alphas,
168 int const& PerturbativeOrder,
169 double const& Ci = 1,
170 double const& IntEps = 1e-7);
171
187 std::function<double(double const&, double const&, double const&)> BuildTmdJet(std::map<int, TmdObjects> const& TmdObj,
188 JetAlgorithm const& JetAlgo,
189 double const& JetR,
190 std::function<double(double const&)> const& Alphas,
191 int const& PerturbativeOrder,
192 double const& CJ = 1,
193 double const& Ci = 1,
194 double const& IntEps = 1e-7);
195
206 std::function<Set<Distribution>(double const&)> MatchTmdPDFs(std::map<int, TmdObjects> const& TmdObj,
207 std::function<Set<Distribution>(double const&)> const& CollPDFs,
208 std::function<double(double const&)> const& Alphas,
209 int const& PerturbativeOrder,
210 double const& Ci = 1);
211
222 std::function<Set<Distribution>(double const&)> MatchTmdFFs(std::map<int, TmdObjects> const& TmdObj,
223 std::function<Set<Distribution>(double const&)> const& CollFFs,
224 std::function<double(double const&)> const& Alphas,
225 int const& PerturbativeOrder,
226 double const& Ci = 1);
227
242 std::function<double(double const&, double const&)> MatchTmdJet(std::map<int, TmdObjects> const& TmdObj,
243 JetAlgorithm const& JetAlgo,
244 double const& tR,
245 std::function<double(double const&)> const& Alphas,
246 int const& PerturbativeOrder,
247 double const& CJ = 1,
248 double const& Ci = 1,
249 double const& IntEps = 1e-7);
250
261 std::function<Set<Operator>(double const&)> MatchingFunctionsPDFs(std::map<int, TmdObjects> const& TmdObj,
262 std::function<double(double const&)> const& Alphas,
263 int const& PerturbativeOrder,
264 double const& Ci = 1);
265
276 std::function<Set<Operator>(double const&)> MatchingFunctionsFFs(std::map<int, TmdObjects> const& TmdObj,
277 std::function<double(double const&)> const& Alphas,
278 int const& PerturbativeOrder,
279 double const& Ci = 1);
280
294 std::function<std::vector<double>(double const&, double const&, double const&)> EvolutionFactors(std::map<int, TmdObjects> const& TmdObj,
295 std::function<double(double const&)> const& Alphas,
296 int const& PerturbativeOrder,
297 double const& Ci = 1,
298 double const& IntEps = 1e-7);
299
316 std::function<std::vector<double>(double const&, double const&, double const&)> EvolutionFactorsK(std::map<int, TmdObjects> const& TmdObj,
317 std::function<double(double const&)> const& Alphas,
318 int const& PerturbativeOrder,
319 double const& Ci = 1,
320 double const& IntEps = 1e-7);
321
334 std::function<double(double const&, double const&, double const&)> QuarkEvolutionFactor(std::map<int, TmdObjects> const& TmdObj,
335 std::function<double(double const&)> const& Alphas,
336 int const& PerturbativeOrder,
337 double const& Ci = 1,
338 double const& IntEps = 1e-7);
339
354 std::function<double(double const&, double const&, double const&)> QuarkEvolutionFactorxi(std::map<int, TmdObjects> const& TmdObj,
355 std::function<double(double const&)> const& Alphas,
356 int const& PerturbativeOrder,
357 double const& xi = 1,
358 double const& Ci = 1,
359 double const& IntEps = 1e-7);
360
373 std::function<double(double const&, double const&, double const&)> GluonEvolutionFactor(std::map<int, TmdObjects> const& TmdObj,
374 std::function<double(double const&)> const& Alphas,
375 int const& PerturbativeOrder,
376 double const& Ci = 1,
377 double const& IntEps = 1e-7);
378
389 std::function<double(double const&)> QuarkAnalyticEvolutionFactor(TmdObjects const& TmdObj,
390 double const& mu0,
391 double const& Alphas0,
392 double const& kappa,
393 double const& kappa0,
394 int const& PerturbativeOrder);
395
406 std::function<double(double const&)> GluonAnalyticEvolutionFactor(TmdObjects const& TmdObj,
407 double const& mu0,
408 double const& Alphas0,
409 double const& kappa,
410 double const& kappa0,
411 int const& PerturbativeOrder);
412
425 std::function<double(double const&, double const&)> CollinsSoperKernel(std::map<int, TmdObjects> const& TmdObj,
426 std::function<double(double const&)> const& Alphas,
427 int const& PerturbativeOrder,
428 double const& Ci = 1,
429 double const& IntEps = 1e-7);
430
440 std::function<double(double const&)> HardFactor(std::string const& Process,
441 std::map<int, TmdObjects> const& TmdObj,
442 std::function<double(double const&)> const& Alphas,
443 int const& PerturbativeOrder,
444 double const& Cf = 1);
446}
apfel::Grid
The Grid class defines ab object that is essentially a collection of "SubGrid" objects plus other glo...
Definition grid.h:22
apfel::Set
The Set template class allocates a collection of objects of type T along the ConvolutionMap and provi...
Definition set.h:22
apfel
Namespace for all APFEL++ functions and classes.
Definition alphaqcd.h:14
apfel::QuarkAnalyticEvolutionFactor
std::function< double(double const &) QuarkAnalyticEvolutionFactor)(TmdObjects const &TmdObj, double const &mu0, double const &Alphas0, double const &kappa, double const &kappa0, int const &PerturbativeOrder)
Analytic evolution factor for the quark TMD.
apfel::InitializeTmdObjectsDYResScheme
std::map< int, TmdObjects > InitializeTmdObjectsDYResScheme(Grid const &g, std::vector< double > const &Thresholds, double const &IntEps=1e-5)
The InitializeTmdObjectsDYResScheme function precomputes the perturbative coefficients required for t...
apfel::MatchTmdPDFs
std::function< Set< Distribution >(double const &) MatchTmdPDFs)(std::map< int, TmdObjects > const &TmdObj, std::function< Set< Distribution >(double const &)> const &CollPDFs, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1)
Function that returns the matched TMD PDFs in b-space.
apfel::InitializeTmdObjectsSivers
std::map< int, TmdObjects > InitializeTmdObjectsSivers(Grid const &g, std::vector< double > const &Thresholds, double const &IntEps=1e-5)
The InitializeTmdObjectsSivers function precomputes the perturbative coefficients required for the ev...
apfel::BuildTmdPDFs
std::function< Set< Distribution >(double const &, double const &, double const &) BuildTmdPDFs)(std::map< int, TmdObjects > const &TmdObj, std::function< Set< Distribution >(double const &)> const &CollPDFs, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the matched and evolved TMD PDFs in b-space as functions of the final scale and...
apfel::QuarkEvolutionFactor
std::function< double(double const &, double const &, double const &) QuarkEvolutionFactor)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the evolution factor for quarks.
apfel::MatchTmdJet
std::function< double(double const &, double const &) MatchTmdJet)(std::map< int, TmdObjects > const &TmdObj, JetAlgorithm const &JetAlgo, double const &tR, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &CJ=1, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the jet TMD in b-space at the initial scale.
apfel::BuildTmdFFs
std::function< Set< Distribution >(double const &, double const &, double const &) BuildTmdFFs)(std::map< int, TmdObjects > const &TmdObj, std::function< Set< Distribution >(double const &)> const &CollFFs, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the matched and evolved TMD FFs in b-space as functions of the final scale and ...
apfel::EvolutionFactorsK
std::function< std::vector< double >(double const &, double const &, double const &) EvolutionFactorsK)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the evolution factors for gluon and quarks. As compared to "EvolutionFactors",...
apfel::GluonAnalyticEvolutionFactor
std::function< double(double const &) GluonAnalyticEvolutionFactor)(TmdObjects const &TmdObj, double const &mu0, double const &Alphas0, double const &kappa, double const &kappa0, int const &PerturbativeOrder)
Analytic evolution factor for the gluon TMD.
apfel::HardFactor
std::function< double(double const &) HardFactor)(std::string const &Process, std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Cf=1)
Function that returns the hard factor.
apfel::MatchingFunctionsFFs
std::function< Set< Operator >(double const &) MatchingFunctionsFFs)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1)
Function that returns the mathing functions for the TMD FFs.
apfel::InitializeTmdObjects
std::map< int, TmdObjects > InitializeTmdObjects(Grid const &g, std::vector< double > const &Thresholds, double const &IntEps=1e-5)
The InitializeTmdObjects function precomputes the perturbative coefficients required for the evolutio...
apfel::MatchingFunctionsPDFs
std::function< Set< Operator >(double const &) MatchingFunctionsPDFs)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1)
Function that returns the mathing functions for the TMD PDFs.
apfel::InitializeTmdObjectsg1
std::map< int, TmdObjects > InitializeTmdObjectsg1(Grid const &g, std::vector< double > const &Thresholds, double const &IntEps=1e-5)
The InitializeTmdObjects function precomputes the perturbative coefficients required for the evolutio...
apfel::GluonEvolutionFactor
std::function< double(double const &, double const &, double const &) GluonEvolutionFactor)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the evolution factor for the gluon.
apfel::CollinsSoperKernel
std::function< double(double const &, double const &) CollinsSoperKernel)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the perturbative part of the Collins-Soper kernel for quarks.
apfel::JetAlgorithm
JetAlgorithm
Enumerator for the jet algoritms for the jet TMDs.
Definition tmdbuilder.h:21
apfel::KT
@ KT
Definition tmdbuilder.h:21
apfel::CONE
@ CONE
Definition tmdbuilder.h:21
apfel::MatchTmdFFs
std::function< Set< Distribution >(double const &) MatchTmdFFs)(std::map< int, TmdObjects > const &TmdObj, std::function< Set< Distribution >(double const &)> const &CollFFs, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1)
Function that returns the matched TMD FFs in b-space.
apfel::InitializeTmdObjectsBM
std::map< int, TmdObjects > InitializeTmdObjectsBM(Grid const &g, std::vector< double > const &Thresholds, double const &IntEps=1e-5)
The InitializeTmdObjectsBM function precomputes the perturbative coefficients required for the evolut...
apfel::QuarkEvolutionFactorxi
std::function< double(double const &, double const &, double const &) QuarkEvolutionFactorxi)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &xi=1, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the evolution factor for quarks with explicit dependence on the resummation-sca...
apfel::EvolutionFactors
std::function< std::vector< double >(double const &, double const &, double const &) EvolutionFactors)(std::map< int, TmdObjects > const &TmdObj, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the evolution factors for gluon and quarks.
apfel::BuildTmdJet
std::function< double(double const &, double const &, double const &) BuildTmdJet)(std::map< int, TmdObjects > const &TmdObj, JetAlgorithm const &JetAlgo, double const &JetR, std::function< double(double const &)> const &Alphas, int const &PerturbativeOrder, double const &CJ=1, double const &Ci=1, double const &IntEps=1e-7)
Function that returns the TMD of a jet in b-space as functions of the final scale and rapidity.
apfel::TmdObjects
Structure that contains all precomputed quantities needed to perform the TMD evolution,...
Definition tmdbuilder.h:31
apfel::TmdObjects::Beta
std::map< int, double > Beta
Definition tmdbuilder.h:33
apfel::TmdObjects::GammaFg
std::map< int, double > GammaFg
Definition tmdbuilder.h:35
apfel::TmdObjects::MatchingFunctionsFFs
std::map< int, std::vector< Set< Operator > > > MatchingFunctionsFFs
Definition tmdbuilder.h:39
apfel::TmdObjects::GammaFq
std::map< int, double > GammaFq
Definition tmdbuilder.h:34
apfel::TmdObjects::GammaK
std::map< int, double > GammaK
Definition tmdbuilder.h:36
apfel::TmdObjects::KCS
std::map< int, std::vector< double > > KCS
Definition tmdbuilder.h:37
apfel::TmdObjects::MatchingFunctionsPDFs
std::map< int, std::vector< Set< Operator > > > MatchingFunctionsPDFs
Definition tmdbuilder.h:38
apfel::TmdObjects::HardFactors
std::map< std::string, std::map< int, double > > HardFactors
Definition tmdbuilder.h:40
apfel::TmdObjects::Threshold
double Threshold
Definition tmdbuilder.h:32
Generated on Wed Apr 3 2024 18:46:06 for APFEL by doxygen 1.10.0