Page MenuHomeHEPForge
Files F7879916

No OneTemporary
Actions

View Options

Index: trunk/src/lhapdf/LHAPDFWrap.cpp
===================================================================
--- trunk/src/lhapdf/LHAPDFWrap.cpp (revision 5917)
+++ trunk/src/lhapdf/LHAPDFWrap.cpp (revision 5918)
@@ -1,77 +1,77 @@
#include "LHAPDF/LHAPDF.h"
#include "LHAPDF/Info.h"
#include "LHAPDF/Exceptions.h"
using namespace LHAPDF;
using namespace std;
extern "C" {
LHAPDF::PDF* lhapdf_init_pdf (char* setname, const int imem) {
LHAPDF::PDF* pdf = LHAPDF::mkPDF(setname, imem);
return pdf;
}
void lhapdf_pdf_delete (LHAPDF::PDF* pdf) {
delete pdf;
}
double lhapdf_pdf_getxmin (const LHAPDF::PDF* pdf) {
return pdf->info().get_entry_as<double>("XMin");
}
double lhapdf_pdf_getxmax (const LHAPDF::PDF* pdf) {
return pdf->info().get_entry_as<double>("XMax");
}
double lhapdf_pdf_getq2min (const LHAPDF::PDF* pdf) {
return LHAPDF::sqr(pdf->info().get_entry_as<double>("QMin"));
}
double lhapdf_pdf_getq2max (const LHAPDF::PDF* pdf) {
return LHAPDF::sqr(pdf->info().get_entry_as<double>("QMax"));
}
bool lhapdf_has_photon (const LHAPDF::PDF* pdf) {
return pdf->hasFlavor(22);
}
/// Get xf(x) values for common partons from PDF pdf
// Evaluate for the 13 LHAPDF5 standard partons
- void lhapdf_evolvepdfm (const LHAPDF::PDF* pdf, const double& x, const double& q, double* fxq) {
+ void lhapdf_evolvepdfm (const LHAPDF::PDF* pdf, const double x, const double q, double* fxq) {
for (size_t i = 0; i < 13; ++i) {
fxq[i] = pdf->xfxQ(i-6, x, q);
}
}
/// Get xfx values from current PDF, including an extra photon flavor
- void lhapdf_evolvepdfphotonm (const LHAPDF::PDF* pdf, const double& x, const double& q, double* fxq, double& photonfxq) {
+ void lhapdf_evolvepdfphotonm (const LHAPDF::PDF* pdf, const double x, const double q, double* fxq, double photonfxq) {
lhapdf_evolvepdfm (pdf, x, q, fxq);
photonfxq = pdf->xfxQ(22, x, q);
}
- void lhapdf_evolvepdfpm (const LHAPDF::PDF* pdf, const double& x, const double& q, const double& s, const int& scheme, double& fxq) {
+ void lhapdf_evolvepdfpm (const LHAPDF::PDF* pdf, const double x, const double q, const double s, const int scheme, double fxq) {
throw LHAPDF::NotImplementedError("Photon structure function are not yet supported");
}
- double lhapdf_getqmass (const LHAPDF::PDF* pdf, const int& nf) {
+ double lhapdf_getqmass (const LHAPDF::PDF* pdf, const int nf) {
if (nf*nf == 1) return pdf->info().get_entry_as<double>("MDown");
else if (nf*nf == 4) return pdf->info().get_entry_as<double>("MUp");
else if (nf*nf == 9) return pdf->info().get_entry_as<double>("MStrange");
else if (nf*nf == 16) return pdf->info().get_entry_as<double>("MCharm");
else if (nf*nf == 25) return pdf->info().get_entry_as<double>("MBottom");
else if (nf*nf == 36) return pdf->info().get_entry_as<double>("MTop");
else throw LHAPDF::UserError("Trying to get quark mass for invalid quark ID #" + LHAPDF::to_str(nf));
}
- int lhapdf_numpdfm (const LHAPDF::PDF* pdf, int& numpdf) {
+ int lhapdf_numpdfm (const LHAPDF::PDF* pdf, int numpdf) {
numpdf = pdf->info().get_entry_as<int>("NumMembers");
if (numpdf > 1) numpdf-=1;
return numpdf;
}
- double lhapdf_alphaspdf (const LHAPDF::PDF* pdf, const double& q) {
+ double lhapdf_alphaspdf (const LHAPDF::PDF* pdf, const double q) {
return pdf->alphasQ(q);
}
}
Index: trunk/share/tests/ref-output/sf_lhapdf6_2.ref
===================================================================
--- trunk/share/tests/ref-output/sf_lhapdf6_2.ref (revision 5917)
+++ trunk/share/tests/ref-output/sf_lhapdf6_2.ref (revision 5918)
@@ -1,193 +1,193 @@
* Test output: sf_lhapdf_2
* Purpose: initialize and fill test structure function object
* Initialize configuration data
* Initialize structure-function object
LHAPDF data:
flavor = f(2212)
prefix = <empty (non-verbose version)>
- file = cteq6ll.LHpdf
+ file = CT10
member = 0
- x(min) = 1.000000000000E-06
+ x(min) = 1.000000000000E-08
x(max) = 1.000000000000E+00
Q(min) = 1.300000000000E+00
- Q(max) = 1.000000000000E+04
+ Q(max) = 1.000000000000E+05
invert = F
mask = T T T T T T * T * T T T T T T
photon mask = T
hoppet_b = F
SF instance: [initialized]
beam = 1
incoming = 1
radiated = 2
outgoing = 3
Interaction: 1
Incoming:
Particle 1
[momentum undefined]
mask [fch] = [FFgT]
internal links: X => 2 3
Outgoing:
Particle 2
[momentum undefined]
mask [fch] = [FFgT]
internal links: 1 => X
Particle 3
[momentum undefined]
mask [fch] = [FFgT]
internal links: 1 => X
State matrix: norm = 1.000000000000E+00
[f(2212)]
[f(92) c(1 )]
[f(-6) c(-1 )] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-5) c(-1 )] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-4) c(-1 )] => ME(3) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-3) c(-1 )] => ME(4) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-2) c(-1 )] => ME(5) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-1) c(-1 )] => ME(6) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(93) c(1 -2)]
[f(21) c(2 -1)] => ME(7) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-92) c(-1 )]
[f(1) c(1 )] => ME(8) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(2) c(1 )] => ME(9) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(3) c(1 )] => ME(10) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(4) c(1 )] => ME(11) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(5) c(1 )] => ME(12) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(6) c(1 )] => ME(13) = ( 0.000000000000E+00, 0.000000000000E+00)
* Initialize incoming momentum with E=500
E = 5.000000000000E+02
P = 0.000000000000E+00 0.000000000000E+00 5.000000000000E+02
* Set kinematics for x=0.5
LHAPDF data:
flavor = f(2212)
prefix = <empty (non-verbose version)>
- file = cteq6ll.LHpdf
+ file = CT10
member = 0
- x(min) = 1.000000000000E-06
+ x(min) = 1.000000000000E-08
x(max) = 1.000000000000E+00
Q(min) = 1.300000000000E+00
- Q(max) = 1.000000000000E+04
+ Q(max) = 1.000000000000E+05
invert = F
mask = T T T T T T * T * T T T T T T
photon mask = T
hoppet_b = F
SF parameters:
x = 5.0000000000E-01
SF instance: [kinematics set]
beam = 1
incoming = 1
radiated = 2
outgoing = 3
Interaction: 1
Incoming:
Particle 1
E = 5.000000000000E+02
P = 0.000000000000E+00 0.000000000000E+00 5.000000000000E+02
mask [fch] = [FFgT]
internal links: X => 2 3
Outgoing:
Particle 2
E = 2.500000000000E+02
P = 0.000000000000E+00 0.000000000000E+00 2.500000000000E+02
mask [fch] = [FFgT]
internal links: 1 => X
Particle 3
E = 2.500000000000E+02
P = 0.000000000000E+00 0.000000000000E+00 2.500000000000E+02
mask [fch] = [FFgT]
internal links: 1 => X
State matrix: norm = 1.000000000000E+00
[f(2212)]
[f(92) c(1 )]
[f(-6) c(-1 )] => ME(1) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-5) c(-1 )] => ME(2) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-4) c(-1 )] => ME(3) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-3) c(-1 )] => ME(4) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-2) c(-1 )] => ME(5) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-1) c(-1 )] => ME(6) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(93) c(1 -2)]
[f(21) c(2 -1)] => ME(7) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(-92) c(-1 )]
[f(1) c(1 )] => ME(8) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(2) c(1 )] => ME(9) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(3) c(1 )] => ME(10) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(4) c(1 )] => ME(11) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(5) c(1 )] => ME(12) = ( 0.000000000000E+00, 0.000000000000E+00)
[f(6) c(1 )] => ME(13) = ( 0.000000000000E+00, 0.000000000000E+00)
x = 0.5000000
f = 1.0000000
* Recover x from momenta
x = 0.5000000
* Evaluate for Q = 100 GeV
LHAPDF data:
flavor = f(2212)
prefix = <empty (non-verbose version)>
- file = cteq6ll.LHpdf
+ file = CT10
member = 0
- x(min) = 1.000000000000E-06
+ x(min) = 1.000000000000E-08
x(max) = 1.000000000000E+00
Q(min) = 1.300000000000E+00
- Q(max) = 1.000000000000E+04
+ Q(max) = 1.000000000000E+05
invert = F
mask = T T T T T T * T * T T T T T T
photon mask = T
hoppet_b = F
SF parameters:
x = 5.0000000000E-01
Q = 1.0000000000E+02
SF instance: [evaluated]
beam = 1
incoming = 1
radiated = 2
outgoing = 3
Interaction: 2
Incoming:
Particle 1
E = 5.00000000E+02
P = 0.00000000E+00 0.00000000E+00 5.00000000E+02
mask [fch] = [FFgT]
internal links: X => 2 3
Outgoing:
Particle 2
E = 2.50000000E+02
P = 0.00000000E+00 0.00000000E+00 2.50000000E+02
mask [fch] = [FFgT]
internal links: 1 => X
Particle 3
E = 2.50000000E+02
P = 0.00000000E+00 0.00000000E+00 2.50000000E+02
mask [fch] = [FFgT]
internal links: 1 => X
State matrix: norm = 1.0000000000E+00
[f(2212)]
[f(92) c(1 )]
[f(-6) c(-1 )] => ME(1) = ( 0.0000000000E+00, 0.0000000000E+00)
- [f(-5) c(-1 )] => ME(2) = ( 1.1259893344E-04, 0.0000000000E+00)
- [f(-4) c(-1 )] => ME(3) = ( 2.2557712257E-04, 0.0000000000E+00)
- [f(-3) c(-1 )] => ME(4) = ( 7.1301796617E-04, 0.0000000000E+00)
- [f(-2) c(-1 )] => ME(5) = ( 1.9333939912E-03, 0.0000000000E+00)
- [f(-1) c(-1 )] => ME(6) = ( 9.5496447194E-04, 0.0000000000E+00)
+ [f(-5) c(-1 )] => ME(2) = ( 3.0097859775E-04, 0.0000000000E+00)
+ [f(-4) c(-1 )] => ME(3) = ( 7.2905446546E-04, 0.0000000000E+00)
+ [f(-3) c(-1 )] => ME(4) = ( 1.0563432274E-03, 0.0000000000E+00)
+ [f(-2) c(-1 )] => ME(5) = ( 1.8057880542E-03, 0.0000000000E+00)
+ [f(-1) c(-1 )] => ME(6) = ( 7.7501337532E-04, 0.0000000000E+00)
[f(93) c(1 -2)]
- [f(21) c(2 -1)] => ME(7) = ( 1.1698522829E-02, 0.0000000000E+00)
+ [f(21) c(2 -1)] => ME(7) = ( 2.9267411026E-02, 0.0000000000E+00)
[f(-92) c(-1 )]
- [f(1) c(1 )] => ME(8) = ( 6.0550918734E-02, 0.0000000000E+00)
- [f(2) c(1 )] => ME(9) = ( 2.3995914039E-01, 0.0000000000E+00)
- [f(3) c(1 )] => ME(10) = ( 7.1301796617E-04, 0.0000000000E+00)
- [f(4) c(1 )] => ME(11) = ( 2.2557712257E-04, 0.0000000000E+00)
- [f(5) c(1 )] => ME(12) = ( 1.1259893344E-04, 0.0000000000E+00)
+ [f(1) c(1 )] => ME(8) = ( 5.7527640865E-02, 0.0000000000E+00)
+ [f(2) c(1 )] => ME(9) = ( 2.1108632693E-01, 0.0000000000E+00)
+ [f(3) c(1 )] => ME(10) = ( 1.0563432274E-03, 0.0000000000E+00)
+ [f(4) c(1 )] => ME(11) = ( 7.2905446546E-04, 0.0000000000E+00)
+ [f(5) c(1 )] => ME(12) = ( 3.0097859775E-04, 0.0000000000E+00)
[f(6) c(1 )] => ME(13) = ( 0.0000000000E+00, 0.0000000000E+00)
* Cleanup
* Test output end: sf_lhapdf_2
Index: trunk/share/tests/ref-output/sf_lhapdf6_3.ref
===================================================================
--- trunk/share/tests/ref-output/sf_lhapdf6_3.ref (revision 5917)
+++ trunk/share/tests/ref-output/sf_lhapdf6_3.ref (revision 5918)
@@ -1,18 +1,18 @@
* Test output: sf_lhapdf_3
* Purpose: initialize and evaluate alpha_s
* Initialize configuration data
* Initialize qcd object
QCD parameters (lhapdf):
- PDF set = cteq6ll.LHpdf
+ PDF set = CT10
PDF member = 1
* Evaluate for Q = 100
- alpha = 0.12791
+ alpha = 0.11639
* Cleanup
* Test output end: sf_lhapdf_3

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 9:22 PM (1 d, 51 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3806263
Default Alt Text
(11 KB)

Event Timeline