lhapdf 0.2.4

(Unofficial) Rust bindings for the LHAPDF C++ library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132

#ifndef WRAPPERS_HPP
#define WRAPPERS_HPP

#ifdef FAKE_WRAPPERS
#include "fake-lhapdf.hpp"
#else
#include <LHAPDF/LHAPDF.h>
#endif

#include <lhapdf/src/lib.rs.h>
#include <rust/cxx.h>

#include <cstdint>
#include <memory>
#include <string>
#include <vector>

#ifdef FAKE_WRAPPERS

inline std::unique_ptr<LHAPDF::PDF> pdf_with_setname_and_member(std::string const&, std::int32_t) {
    return std::unique_ptr<LHAPDF::PDF>();
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_set_and_member(LHAPDF::PDFSet const&, std::int32_t) {
    return std::unique_ptr<LHAPDF::PDF>();
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_setname_and_nmem(std::string const&) {
    return std::unique_ptr<LHAPDF::PDF>();
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_lhaid(std::int32_t) {
    return std::unique_ptr<LHAPDF::PDF>();
}

inline std::unique_ptr<LHAPDF::PDFSet> pdfset_new(std::string const&) {
    return std::unique_ptr<LHAPDF::PDFSet>();
}

inline std::unique_ptr<LHAPDF::PDFSet> pdfset_from_pdf(LHAPDF::PDF const&) {
    return std::unique_ptr<LHAPDF::PDFSet>();
}

inline void lookup_pdf_setname(std::int32_t, std::string&) {}

inline std::int32_t lookup_pdf_memberid(std::int32_t) {
    return 0;
}

inline void get_pdfset_error_type(LHAPDF::PDFSet const&, std::string&) {}

inline PdfUncertainty pdf_uncertainty(
    LHAPDF::PDFSet const&,
    rust::Slice<double const>,
    double,
    bool
) {
    PdfUncertainty result;
    return result;
}

#else

inline std::unique_ptr<LHAPDF::PDF> pdf_with_setname_and_member(
    std::string const& setname,
    std::int32_t member
) {
    return std::unique_ptr<LHAPDF::PDF>(LHAPDF::mkPDF(setname, member));
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_set_and_member(
    LHAPDF::PDFSet const& set,
    std::int32_t member
) {
    return pdf_with_setname_and_member(set.name(), member);
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_setname_and_nmem(std::string const& setname_nmem) {
    return std::unique_ptr<LHAPDF::PDF>(LHAPDF::mkPDF(setname_nmem));
}

inline std::unique_ptr<LHAPDF::PDF> pdf_with_lhaid(std::int32_t lhaid) {
    return std::unique_ptr<LHAPDF::PDF>(LHAPDF::mkPDF(lhaid));
}

inline std::unique_ptr<LHAPDF::PDFSet> pdfset_new(std::string const& setname) {
    return std::unique_ptr<LHAPDF::PDFSet>(new LHAPDF::PDFSet(setname));
}

inline std::unique_ptr<LHAPDF::PDFSet> pdfset_from_pdf(LHAPDF::PDF const& pdf) {
    return std::unique_ptr<LHAPDF::PDFSet>(new LHAPDF::PDFSet(pdf.set()));
}

inline void lookup_pdf_setname(std::int32_t lhaid, std::string& setname) {
    setname = LHAPDF::lookupPDF(lhaid).first;
}

inline std::int32_t lookup_pdf_memberid(std::int32_t lhaid) {
    return LHAPDF::lookupPDF(lhaid).second;
}

inline void get_pdfset_error_type(LHAPDF::PDFSet const& set, std::string& error_type) {
    error_type = set.errorType();
}

inline PdfUncertainty pdf_uncertainty(
    LHAPDF::PDFSet const& pdfset,
    rust::Slice<double const> values,
    double cl,
    bool alternative
) {
    std::vector<double> const vector(values.begin(), values.end());
    auto const uncertainty = pdfset.uncertainty(vector, cl, alternative);

    // convert the C++ `PDFUncertainty` to Rust's `PdfUncertainty`
    PdfUncertainty result;
    result.central = uncertainty.central;
    result.errplus = uncertainty.errplus;
    result.errminus = uncertainty.errminus;
    result.errsymm = uncertainty.errsymm;
    result.scale = uncertainty.scale;
    result.errplus_pdf = uncertainty.errplus_pdf;
    result.errminus_pdf = uncertainty.errminus_pdf;
    result.errsymm_pdf = uncertainty.errsymm_pdf;
    result.err_par = uncertainty.err_par;

    return result;
}

#endif

#endif