scarrazza / apfel Goto Github PK

View Code? Open in Web Editor NEW

12.0 12.0 9.0 33.15 MB

A PDF Evolution Library

Home Page: http://apfel.mi.infn.it/

License: GNU General Public License v3.0

C++ 12.40% C 0.44% Python 0.50% Shell 10.85% Fortran 75.04% Makefile 0.14% SWIG 0.01% CMake 0.62%

apfel's People

Contributors

Stargazers

Watchers

Forkers

claria gitter-badger mariaubiali dnapoletano vbertone amsimone adrianneschauss neogyk andriish veprbl

apfel's Issues

EVOL: VFNS with scale variations

Dear @vbertone,
while doing our usual benchmarks with the LHA tables and EKO, we noticed that there is a discrepancy between Apfel and the reported values in case of NNLO with scale variations.
We believe that this is due to the strong coupling at the the threshold scale, which is not computed correctly.
In particular here, we propose the change:

asthUp(inf)   = a_QCD(  m2th(inf) * ( 1d0 + eps ) )
asthDown(inf) = a_QCD(  m2th(inf) * ( 1d0 - eps ) )

Since a_QCD takes factorization scales as input and m2th(inf) are already in factorization scales.
With this fix, we are able to find good agreement.

DIS: TMC for F3

We believe TMC for F3 is wrong:
apfel/src/DIS/F3light.f

Lines 40 to 41 in a6bd3c1

c1 = x**2 / xi**2 / tau

c2 = 4d0 * rhop * x**3 / tau**1.5d0
APFEL Implementation corresponds to https://arxiv.org/pdf/0808.1231.pdf sec. 3.8
we believe Eq. 24 of http://iopscience.iop.org/0954-3899/35/5/053101 to be correct
see also discussion in sec. C1 (ii) in that paper: they claim a typo in the original paper by Georgi & Politzer
in esp. c2 should be 2 * rho * x^3 / tau^1.5

can not import in Python

I installed APFEL and tried to use it in Python

import apfel

and see the following error

ModuleNotFoundError: No module named 'apfel'

However, import in C++ has no issues

#include "APFEL/APFEL.h"

Could you please point out to me what I should do to make APFEL importable in Python? Thanks!

DIS: TMC for FL

We believe TMC for FL is wrong

apfel/src/DIS/FLlight.f

Lines 56 to 58 in a6bd3c1

    
                  FLlight = FLlight + w_int_gen(n,alpha,xi)  
        
           1           * ( FL(3,0,alpha)  
        
           2           + c1 * F2(3,0,alpha) + c2 * I2(3,0,alpha) )

APFEL Implementation corresponds to https://arxiv.org/pdf/0808.1231.pdf sec. 3.8
we believe Eq. 26 of http://iopscience.iop.org/0954-3899/35/5/053101 to be correct
APFEL is not producing physical results for this - see log below

log

FLlight with theory=109 using CT14llo_NF6
           x          Q2         APFEL        yadism  yadism_error    rel_err[%]
0   0.000317   90.000000 -1.096388e-09  6.453450e-09           0.0   -688.610175
1   0.001007   90.000000 -7.971752e-09  4.682564e-08           0.0   -687.394585
2   0.003195   90.000000 -5.319054e-08  3.368383e-07           0.0   -733.267388
3   0.010138   90.000000 -3.578034e-07  2.480502e-06           0.0   -793.258370
4   0.032169   90.000000 -3.500474e-06  1.894879e-05           0.0   -641.320747
5   0.102077   90.000000 -4.737209e-05  1.355193e-04           0.0   -386.074071
6   0.304545   90.000000 -5.363571e-04  5.690540e-04           0.0   -206.096106
7   0.536364   90.000000 -6.862291e-04  3.825738e-04           0.0   -155.750147
8   0.768182   90.000000 -1.603016e-04  4.102670e-05           0.0   -125.593442
9   1.000000   90.000000 -1.977082e-07  5.276186e-09           0.0   -102.668674
10  0.010000    2.000000 -4.734504e-07  7.321240e-05           0.0 -15563.583255
11  0.010000    3.088904 -1.401287e-06  5.034189e-05           0.0  -3692.547885
12  0.010000    4.770665 -1.653054e-06  3.449513e-05           0.0  -2186.751469
13  0.010000    7.368063 -1.562449e-06  2.353910e-05           0.0  -1606.551240
14  0.010000   11.379620 -1.330316e-06  1.599546e-05           0.0  -1302.380041
15  0.010000   17.575279 -1.065192e-06  1.082591e-05           0.0  -1116.334043
16  0.010000   27.144176 -8.190154e-07  7.299937e-06           0.0   -991.306470
17  0.010000   41.922880 -6.119055e-07  4.905680e-06           0.0   -901.705498
18  0.010000   64.747880 -4.474956e-07  3.286587e-06           0.0   -834.440182
19  0.010000  100.000000 -3.219011e-07  2.195792e-06           0.0   -782.132479

corresponding theory

{'ID': 22,
 'PTO': 0,
 'FNS': 'FFNS',
 'DAMP': 0,
 'IC': 0,
 'ModEv': 'EXA',
 'XIR': 1.0,
 'XIF': 1.0,
 'NfFF': 3,
 'MaxNfAs': 6,
 'MaxNfPdf': 6,
 'Q0': 1,
 'alphas': 0.11800000000000001,
 'Qref': 91.2,
 'QED': 0,
 'alphaqed': 0.007496251999999999,
 'Qedref': 1.777,
 'SxRes': 0,
 'SxOrd': 'LL',
 'HQ': 'POLE',
 'mc': 2,
 'Qmc': 2,
 'kcThr': 1.0,
 'mb': 4,
 'Qmb': 4,
 'kbThr': 1.0,
 'mt': 173.07,
 'Qmt': 173.07,
 'ktThr': 1.0,
 'CKM': '0.97428 0.22530 0.003470 0.22520 0.97345 0.041000 0.00862 0.04030 0.999152',
 'MZ': 91.1876,
 'MW': 80.398,
 'GF': 1.1663787e-05,
 'SIN2TW': 0.23126,
 'TMC': 2,
 'MP': 0.938,
 'Comments': 'LO baseline for small-x res',
 'global_nx': 0,
 'EScaleVar': 1,
 '_modify_time': '2020-05-28 16:55:02.232461'}

Evol: NNLO Polarized

While doing the benchmark of NNLO polarized evolution with Eko,
I found two things which are not fully consistent:

P_{ns,s} is not implemented here, but it's available from this source ref eq 23, 25.
This should have a minor effect on the Valence distribution.
As far as I understand, VFNS evolution is not possible since polarized matching conditions are not implemented.

DIS: IC in F3charm in NC

we believe the IC contributions for F3charm in the NC case are not present

it is discussed in the documentation but not implemented in the code:

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 92 to 150 in fdfd068

    
                          SC2mNC(igrid,ixi,3,0,beta,alpha) = ( 2d0 - eta ) * fL 
        
                          SCLmNC(igrid,ixi,3,0,beta,alpha) = 4d0 * ( 1d0 - eta ) 
        
                1              / ( 2d0 - eta ) / factL * fL 
        
           * 
        
           *     NLO 
        
           * 
        
                          if(ipt.ge.1)then 
        
           * 
        
                             bound = alpha-inter_degree(igrid) 
        
                             if(alpha.lt.inter_degree(igrid)) bound = 0 
        
           * 
        
                             c = max(xg(igrid,beta), 
        
                1                 xg(igrid,beta)/xg(igrid,alpha+1)) / eta 
        
                             d = min(eta,xg(igrid,beta)/xg(igrid,bound)) / eta 
        
           * 
        
                             if(c.ge.1d0) cycle 
        
           * 
        
                             walpha = alpha 
        
                             wbeta  = beta 
        
           * 
        
                             k  = 3 
        
                             wl = 1 
        
           * 
        
                             sf   = 1 
        
                             C2RS = fact2 * eta * dgauss(integrandsICm,c,d,eps) 
        
                             C2L  = fact2 * eta * c21ICL(c) 
        
           * 
        
                             sf   = 2 
        
                             CLRS = eta * dgauss(integrandsICm,c,d,eps) 
        
                             CLL  = eta * cL1ICL(c) 
        
           * 
        
                             if(MassScheme(1:5).eq."FONLL")then 
        
                                sf = 2 
        
           *     Gluon 
        
                                k  = 1 
        
                                gDIC = dgauss(integrandsICm,c,d,eps) 
        
           *     Put these terms in the LO gluon slot (e.g. SC2mNC(igrid,ixi,1,0,beta,alpha)) 
        
           *     which is empty. 
        
                                SC2mNC(igrid,ixi,1,0,beta,alpha) = 
        
                1                    - ( 2d0 - eta ) * lnF * gDIC 
        
                                SCLmNC(igrid,ixi,1,0,beta,alpha) = 
        
                1                    - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) / factL 
        
                2                    * lnF * gDIC 
        
           *     Non-singlet 
        
                                k  = 3 
        
                                wl = 2 
        
                                nsDIC = dgauss(integrandsICm,c,d,eps) 
        
           * 
        
                                C2RS = C2RS - ( 2d0 - eta ) * nsDIC 
        
                                C2L  = C2L - ( 2d0 - eta ) * DICc(xi,c) 
        
           * 
        
                                CLRS = CLRS - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) 
        
                1                    * nsDIC 
        
                                CLL  = CLL - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) 
        
                1                    * DICc(xi,c) 
        
                             endif 
        
           * 
        
                             SC2mNC(igrid,ixi,3,1,beta,alpha) = C2RS + C2L * fL 
        
                             SCLmNC(igrid,ixi,3,1,beta,alpha) = CLRS + CLL * fL

apfel/src/DIS/ComputeDISOperators.f

Lines 469 to 488 in fdfd068

    
                           if(IntrinsicCharm)then 
        
                              if(Nf_FF.lt.4)then 
        
                                 do pt=0,ipt_max_IC 
        
                                    CLnsp(4) = CLnsp(4) + as(pt) 
        
           1                             * ( c0(4) 
        
           2                             * SCLm0NC(jgrid,ixi(4), 
        
           3                             3,pt,alpha,beta) 
        
           4                             + c1(4) 
        
           5                             * SCLm0NC(jgrid,ixi(4)+1, 
        
           6                             3,pt,alpha,beta) ) 
        
                                    C2nsp(4) = C2nsp(4) + as(pt) 
        
           1                             * ( c0(4) 
        
           2                             * SC2m0NC(jgrid,ixi(4), 
        
           3                             3,pt,alpha,beta) 
        
           4                             + c1(4) 
        
           5                             * SC2m0NC(jgrid,ixi(4)+1, 
        
           6                             3,pt,alpha,beta) ) 
        
                                 enddo 
        
                              endif 
        
                           endif

DIS: F2charm is non-zero at low Q2 in FFNS4 for uonly PDF

We find F2charm beeing non-zero when using the uonly PDF (which does what it name says uonly.tar.gz)
We believe this to be wrong

F2charm with theory=99 using uonly
      x        Q2     APFEL  yadism  yadism_error  rel_err[%]
0  0.01  2.000000  0.000260     0.0           0.0      -100.0
1  0.01  2.444444  0.000304     0.0           0.0      -100.0
2  0.01  2.888889  0.000336     0.0           0.0      -100.0
3  0.01  3.333333  0.000356     0.0           0.0      -100.0
4  0.01  3.777778  0.000366     0.0           0.0      -100.0
5  0.01  4.222222  0.000000     0.0           0.0         NaN
6  0.01  4.666667  0.000000     0.0           0.0         NaN
7  0.01  5.111111  0.000000     0.0           0.0         NaN
8  0.01  5.555556  0.000000     0.0           0.0         NaN
9  0.01  6.000000  0.000000     0.0           0.0         NaN

theory 99

{'ID': 22,
 'PTO': 0,
 'FNS': 'FFNS',
 'DAMP': 0,
 'IC': 0,
 'ModEv': 'EXA',
 'XIR': 1.0,
 'XIF': 1.0,
 'NfFF': 4,
 'MaxNfAs': 6,
 'MaxNfPdf': 6,
 'Q0': 1,
 'alphas': 0.11800000000000001,
 'Qref': 91.2,
 'QED': 0,
 'alphaqed': 0.007496251999999999,
 'Qedref': 1.777,
 'SxRes': 0,
 'SxOrd': 'LL',
 'HQ': 'POLE',
 'mc': 2,
 'Qmc': 2,
 'kcThr': 1.0,
 'mb': 4,
 'Qmb': 4,
 'kbThr': 1.0,
 'mt': 173.07,
 'Qmt': 173.07,
 'ktThr': 1.0,
 'CKM': '0.97428 0.22530 0.003470 0.22520 0.97345 0.041000 0.00862 0.04030 0.999152',
 'MZ': 91.1876,
 'MW': 80.398,
 'GF': 1.1663787e-05,
 'SIN2TW': 0.23126,
 'TMC': 0,
 'MP': 0.938,
 'Comments': 'LO baseline for small-x res',
 'global_nx': 0,
 'EScaleVar': 1,
 '_modify_time': '2020-05-28 16:55:02.164765'}

might be related to #22

DIS: CC FONLL-A scale variations

We believe to have found a bug in APFEL:
when computing F2charm/F3charm in FONLL-A scheme (i.e. NLO) for xiF!=0 there is the scale variations for the massive-zero coefficient missing
the actual coefficient is the same as in the NC zero-mass (light) case - so it amounts to a missing 2 in the difference of the scale-variated and non-scale-variated terms (when using a fact-scale independent PDF)
note that our result is continuous in Q2 (at ~16GeV^2), but APFEL not

non-scale variated

           x           Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000317    90.000000  0.001365  0.001365  1.052304e-11   -0.003618
1   0.001007    90.000000  0.003636  0.003636  2.957269e-11   -0.011747
2   0.003195    90.000000  0.009640  0.009638  1.129206e-10   -0.014704
3   0.010138    90.000000  0.025498  0.025494  3.063871e-10   -0.015907
4   0.032169    90.000000  0.067195  0.067179  8.827141e-10   -0.024268
5   0.102077    90.000000  0.173154  0.173063  2.515006e-09   -0.052814
6   0.304545    90.000000  0.376150  0.376114  4.893392e-09   -0.009677
7   0.536364    90.000000  0.460563  0.457903  7.999385e-09   -0.577581
8   0.768182    90.000000  0.354356  0.363529  5.272257e-09    2.588783
9   1.000000    90.000000  0.000000  0.000000  0.000000e+00         NaN
10  0.001000     4.000000  0.007303  0.007293  1.185785e-11   -0.130942
11  0.001000     5.348935  0.006370  0.006368  2.650673e-11   -0.030171
12  0.001000     7.152778  0.005681  0.005677  4.638424e-11   -0.057296
13  0.001000     9.564937  0.005165  0.005157  1.354225e-11   -0.145549
14  0.001000    12.790557  0.004764  0.004761  1.129663e-11   -0.057178
15  0.001000    17.103966  0.004459  0.004458  9.166199e-12   -0.012656
16  0.001000    22.872003  0.004228  0.004228  1.874557e-11    0.005051
17  0.001000    30.585218  0.004046  0.004046  1.615744e-11   -0.010532
18  0.001000    40.899589  0.003899  0.003899  1.228940e-11   -0.010231
19  0.001000    54.692316  0.003778  0.003778  1.823467e-11   -0.007239
20  0.001000    73.136417  0.003677  0.003677  1.288867e-11   -0.001218
21  0.001000    97.800494  0.003592  0.003592  1.575276e-11   -0.011739
22  0.001000   130.782133  0.003519  0.003518  1.637083e-11   -0.014816
23  0.001000   174.886299  0.003454  0.003454  1.217347e-11   -0.015418
24  0.001000   233.863883  0.003397  0.003396  1.171409e-11   -0.014889
25  0.001000   312.730705  0.003345  0.003345  1.807173e-11   -0.016899
26  0.001000   418.194092  0.003299  0.003298  1.244429e-11   -0.018640
27  0.001000   559.223306  0.003256  0.003255  8.366275e-12   -0.018543
28  0.001000   747.812346  0.003216  0.003216  1.052677e-11   -0.018331
29  0.001000  1000.000000  0.003179  0.003179  1.092562e-11   -0.018418

scale variated

           x           Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000317    90.000000  0.001856  0.001615  6.113991e-12  -12.974212
1   0.001007    90.000000  0.004953  0.004306  2.015077e-11  -13.052670
2   0.003195    90.000000  0.013050  0.011373  5.598218e-11  -12.849692
3   0.010138    90.000000  0.033837  0.029727  1.880561e-10  -12.144662
4   0.032169    90.000000  0.085403  0.076386  5.067254e-10  -10.558479
5   0.102077    90.000000  0.202139  0.187519  1.417134e-09   -7.232290
6   0.304545    90.000000  0.372621  0.373441  2.981850e-09    0.220166
7   0.536364    90.000000  0.386012  0.418997  4.154874e-09    8.545153
8   0.768182    90.000000  0.239398  0.305492  2.403995e-09   27.608191
9   1.000000    90.000000  0.000000  0.000000  0.000000e+00         NaN
10  0.001000     4.000000  0.009186  0.009174  5.671767e-12   -0.130288
11  0.001000     5.348935  0.007951  0.007948  2.224786e-11   -0.039398
12  0.001000     7.152778  0.007037  0.007033  2.030852e-11   -0.055678
13  0.001000     9.564937  0.006348  0.006343  6.949554e-12   -0.089371
14  0.001000    12.790557  0.005820  0.005816  4.434129e-12   -0.070409
15  0.001000    17.103966  0.006214  0.005413  1.974105e-11  -12.893534
16  0.001000    22.872003  0.005873  0.005106  2.391705e-11  -13.057690
17  0.001000    30.585218  0.005600  0.004862  2.468689e-11  -13.166203
18  0.001000    40.899589  0.005374  0.004665  2.013577e-11  -13.195396
19  0.001000    54.692316  0.005185  0.004502  2.422744e-11  -13.169462
20  0.001000    73.136417  0.005024  0.004366  1.795298e-11  -13.105272
21  0.001000    97.800494  0.004886  0.004250  2.000938e-11  -13.025702
22  0.001000   130.782133  0.004765  0.004149  1.841735e-11  -12.922172
23  0.001000   174.886299  0.004657  0.004061  1.918520e-11  -12.802721
24  0.001000   233.863883  0.004560  0.003982  1.439796e-11  -12.673503
25  0.001000   312.730705  0.004472  0.003911  1.593528e-11  -12.541315
26  0.001000   418.194092  0.004392  0.003847  1.555100e-11  -12.405532
27  0.001000   559.223306  0.004317  0.003788  1.470485e-11  -12.267555
28  0.001000   747.812346  0.004248  0.003733  1.484172e-11  -12.129212
29  0.001000  1000.000000  0.004183  0.003682  1.069974e-11  -11.992034

diff

           x           Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000317    90.000000 -0.000491 -0.000250  1.663703e-11  -49.071786
1   0.001007    90.000000 -0.001316 -0.000670  4.972346e-11  -49.083270
2   0.003195    90.000000 -0.003410 -0.001735  1.689027e-10  -49.129734
3   0.010138    90.000000 -0.008339 -0.004234  4.944432e-10  -49.229239
4   0.032169    90.000000 -0.018208 -0.009207  1.389440e-09  -49.434442
5   0.102077    90.000000 -0.028984 -0.014457  3.932141e-09  -50.122870
6   0.304545    90.000000  0.003530  0.002673  7.875242e-09  -24.274422
7   0.536364    90.000000  0.074551  0.038906  1.215426e-08  -47.813207
8   0.768182    90.000000  0.114958  0.058038  7.676251e-09  -49.513885
9   1.000000    90.000000  0.000000  0.000000  0.000000e+00    0.000000
10  0.001000     4.000000 -0.001883 -0.001881  1.752961e-11   -0.127752
11  0.001000     5.348935 -0.001581 -0.001580  4.875459e-11   -0.076566
12  0.001000     7.152778 -0.001356 -0.001355  6.669276e-11   -0.048901
13  0.001000     9.564937 -0.001183 -0.001185  2.049181e-11    0.155818
14  0.001000    12.790557 -0.001056 -0.001055  1.573076e-11   -0.130089
15  0.001000    17.103966 -0.001755 -0.000955  2.890725e-11  -45.617532
16  0.001000    22.872003 -0.001645 -0.000878  4.266262e-11  -46.627284
17  0.001000    30.585218 -0.001553 -0.000817  4.084433e-11  -47.432279
18  0.001000    40.899589 -0.001475 -0.000766  3.242517e-11  -48.050363
19  0.001000    54.692316 -0.001407 -0.000724  4.246211e-11  -48.519009
20  0.001000    73.136417 -0.001347 -0.000689  3.084166e-11  -48.878320
21  0.001000    97.800494 -0.001294 -0.000658  3.576214e-11  -49.153844
22  0.001000   130.782133 -0.001246 -0.000631  3.478819e-11  -49.362926
23  0.001000   174.886299 -0.001203 -0.000607  3.135867e-11  -49.518829
24  0.001000   233.863883 -0.001163 -0.000586  2.611204e-11  -49.635516
25  0.001000   312.730705 -0.001127 -0.000567  3.400701e-11  -49.724830
26  0.001000   418.194092 -0.001093 -0.000549  2.799528e-11  -49.791142
27  0.001000   559.223306 -0.001061 -0.000532  2.307112e-11  -49.841461
28  0.001000   747.812346 -0.001032 -0.000517  2.536849e-11  -49.878594
29  0.001000  1000.000000 -0.001004 -0.000503  2.162536e-11  -49.906513

APFEL banner

WARNING: FONLL-A is a VFN scheme ... setting VFNS PDF evolution WARNING: FONLL-A is NLO scheme ... setting NLO perturbative order WARNING: if there are external grids they cannot be locked ... unlocking subgrids

Welcome to
/// //// //// ///_/ _/
_/ / / / / / /
//// //// /// /// _/
_/ _/ _/ _/ _/ _/
/ / / / //// ////
_____v3.0.4 A PDF Evolution Library, arXiv:1310.1394
Authors: V. Bertone, S. Carrazza, J. Rojo

Report of the evolution parameters:

QCD evolution
Space-like evolution (PDFs)
Unpolarized evolution
Evolution scheme: VFNS at N1LO
Solution of the DGLAP equation: 'exactalpha' with maximum 6 active flavours
Solution of the coupling equations: 'exact' with maximum 6 active flavours
Coupling reference value:

AlphaQCD( 91.2000 GeV) = 0.118000
Pole heavy quark masses:
Mc = 2.0000 GeV
Mb = 4.0000 GeV
Mt = 173.0700 GeV
The matching thresholds coincide with the physical masses
muR / muF = 1.0000

Allowed evolution range [ 1.0000 : 10000.0000 ] GeV
Fast evolution enabled

Initialization of the evolution completed in 5.254 s

Report of the electroweak parameters:

Mass of the Z = 91.188 GeV
Mass of the W = 80.398 GeV
Mass of the proton = 0.9380 GeV
sin^2(thetaW) = 0.2313
GFermi = 1.16638E-05
| 0.9743 0.2253 0.0035 |
CKM = | 0.2252 0.9735 0.0410 |
| 0.0086 0.0403 0.9992 |
Z propagator correction = 0.00000

Report of the DIS parameters:

Computation in the FONLL-A mass scheme
Charged Current (CC) process
Scattering electron - proton
muR / Q = 1.0000
muF / Q = 0.5000
(No scale variation in the evolution)
Target Mass corrections disabled
FONLL damping factor disabled for all heavy quarks
Intrinsic charm disabled

Release 3.0.7?

Dear @scarrazza, @vbertone ,

it has been almost 2 years since the last release has been made.
Would that make sense to bump versions in configure/CMakeLists.txt and make a new release?

Best regards,

Andrii

Please remove ApfelGUI from Readme.md

As ApfelGUI has been removed, please remove any mentions from the Readme, and preferably also from the homepage.

Related: #3

DIS: charm pdf in F3charm at CC

We believe there is a bug in combination of CC coefficient functions with the respective PDFs
When computing F3charm with the conly pdf (which is defined by c(x)=x(1-x) and 0 otherwise ) APFEL is 0, while we're not (see table below)
looking at the PDG expression for the LO expressions (=parton-model) the charm pdf can contribute to F3 when scattering electrons or anti-neutrinos (what we are doing here) (see also Eq. 51 of CC_observables_def.pdf)
keep in mind that the APFEL definition of F3charm relies on the seperation of the CKM matrix

output:

------------------------------F3charm------------------------------

           x           Q2  APFEL    yadism  yadism_error  rel_err[%]
0   0.000317    90.000000    0.0  0.001465  6.438219e-12         inf
1   0.001007    90.000000    0.0  0.003751  9.234736e-12         inf
2   0.003195    90.000000    0.0  0.009510  4.041613e-11         inf
3   0.010138    90.000000    0.0  0.024013  1.321150e-10         inf
4   0.032169    90.000000    0.0  0.060758  4.250026e-10         inf
5   0.102077    90.000000    0.0  0.153750  1.063813e-09         inf
6   0.304545    90.000000    0.0  0.346403  1.647991e-09         inf
7   0.536364    90.000000    0.0  0.450932  2.779282e-09         inf
8   0.768182    90.000000    0.0  0.395607  6.832703e-10         inf
9   1.000000    90.000000    0.0  0.000000  0.000000e+00         NaN
10  0.001000     4.000000    0.0  0.008233  9.639907e-12         inf
11  0.001000     7.387505    0.0  0.005930  2.837073e-11         inf
12  0.001000    13.643806    0.0  0.004952  1.346871e-11         inf
13  0.001000    25.198421    0.0  0.004389  1.252950e-11         inf
14  0.001000    46.538363    0.0  0.004015  1.220218e-11         inf
15  0.001000    85.950595    0.0  0.003747  1.481200e-11         inf
16  0.001000   158.740105    0.0  0.003545  1.100350e-11         inf
17  0.001000   293.173318    0.0  0.003386  1.215494e-11         inf
18  0.001000   541.454816    0.0  0.003258  1.348017e-11         inf
19  0.001000  1000.000000    0.0  0.003152  1.028605e-11         inf
--------------------------------------------------------------------

APFEL setup:

   
 WARNING: FFNS is a FFN scheme
          ... setting NF = 3 FFNS PDF evolution
 WARNING: if there are external grids they cannot be locked
          ... unlocking subgrids
 
 Welcome to 
      _/_/_/    _/_/_/_/   _/_/_/_/   _/_/_/_/   _/
    _/    _/   _/    _/   _/         _/         _/
   _/_/_/_/   _/_/_/_/   _/_/_/     _/_/_/     _/
  _/    _/   _/         _/         _/         _/
 _/    _/   _/         _/         _/_/_/_/   _/_/_/_/
 _____v3.0.4 A PDF Evolution Library, arXiv:1310.1394
      Authors: V. Bertone, S. Carrazza, J. Rojo
 
 Report of the evolution parameters:
   
 QCD evolution
 Space-like evolution (PDFs)
 Unpolarized evolution
 Evolution scheme: FFNS with 3 active flavours at N1LO
 Solution of the DGLAP equation: 'exactalpha'
 Solution of the coupling equations: 'exact'
 Coupling reference value:
 - AlphaQCD( 91.2000 GeV) =  0.118000
 Pole heavy quark masses:
 - Mc =   2.0000 GeV
 - Mb =   4.0000 GeV
 - Mt = 173.0700 GeV
 The matching thresholds coincide with the physical masses
 muR / muF =  1.0000
  
 Allowed evolution range [   1.0000 :  10000.0000 ] GeV
 Fast evolution enabled
 
 Initialization of the evolution completed in   1.309 s
  
 Report of the electroweak parameters:
   
 Mass of the Z = 91.188 GeV
 Mass of the W = 80.398 GeV
 Mass of the proton = 0.9380 GeV
 sin^2(thetaW) = 0.2313
 GFermi = 1.16638E-05
       | 0.9743 0.2253 0.0035 |
 CKM = | 0.2252 0.9735 0.0410 |
       | 0.0086 0.0403 0.9992 |
 Z propagator correction = 0.00000
   
 Report of the DIS parameters:
   
 Computation in the FFNS3 mass scheme
 Charged Current (CC) process
 Scattering electron - proton   
 muR / Q =  0.5000
 muF / Q =  1.0000
 (No scale variation in the evolution)
 Target Mass corrections disabled
 Intrinsic charm disabled

Certificate for apfel.hepforge.org has expired

http://apfel.hepforge.org redirects to https://apfel.hepforge.org, but the latter has an expired certificate.

DIS: IC FONLL-A q Matching for NC FL

we believe the NLO IC contributions to FL for FONLL (massive-zero coefficient) is not correctly normalized

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 143 to 146 in fdfd068

    
                           CLRS = CLRS - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) 
        
           1                    * nsDIC 
        
                           CLL  = CLL - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) 
        
           1                    * DICc(xi,c)

these lines should contain an additional division by factL as is present for the gluon:

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 132 to 134 in fdfd068

    
                           SCLmNC(igrid,ixi,1,0,beta,alpha) = 
        
           1                    - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) / factL 
        
           2                    * lnF * gDIC

this is just the LO as prefactor

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 93 to 94 in fdfd068

SCLmNC(igrid,ixi,3,0,beta,alpha) = 4d0 * ( 1d0 - eta )

1 / ( 2d0 - eta ) / factL * fL

Benchmarking to yadism with the error fixed gives:

Log added, hash=3594cca65cabea4b2cec85345986c0fb978b9fedaf5b2dd69ff53368680a9301
 ────────── 
  F2_charm  
 ────────── 
           x         Q2    yadism  yadism_error     APFEL  percent_error
0   0.000100  20.000000  0.000132 -1.629179e-13  0.000132       0.183309
1   0.000278  20.000000  0.000316  1.282237e-12  0.000315       0.112155
2   0.000774  20.000000  0.000749 -5.303289e-12  0.000749       0.020648
3   0.002154  20.000000  0.001766 -5.497683e-12  0.001766       0.008649
4   0.005995  20.000000  0.004159 -2.203334e-11  0.004160      -0.028787
5   0.016681  20.000000  0.009780 -4.266334e-11  0.009806      -0.257337
6   0.046416  20.000000  0.022883 -1.142255e-10  0.023073      -0.824161
7   0.129155  20.000000  0.051691 -6.739887e-11  0.052715      -1.942159
8   0.359381  20.000000  0.095410 -1.148124e-09  0.105458      -9.527575
9   1.000000  20.000000  0.000000  0.000000e+00  0.000000            NaN
10  0.000100   2.280100  0.000190 -5.409928e-13  0.000190      -0.070432
11  0.000278   2.280100  0.000454 -1.627653e-12  0.000454      -0.045346
12  0.000774   2.280100  0.001075 -8.679910e-13  0.001076      -0.130771
13  0.002154   2.280100  0.002523 -8.150368e-12  0.002525      -0.107948
14  0.005995   2.280100  0.005879 -4.668564e-11  0.005880      -0.017036
15  0.016681   2.280100  0.013619 -3.032851e-11  0.013582       0.272896
16  0.046416   2.280100  0.031256 -9.510437e-11  0.030975       0.906939
17  0.129155   2.280100  0.069233 -2.617881e-10  0.067295       2.880376
18  0.359381   2.280100  0.115947 -1.927640e-09  0.116128      -0.155885
19  1.000000   2.280100  0.000000  0.000000e+00  0.000000            NaN
20  0.100000   4.000000  0.059420  1.094726e-10  0.059094       0.550973
21  0.100000   4.783253  0.056225 -3.421848e-10  0.056151       0.131524
22  0.100000   5.719877  0.053445 -3.893411e-10  0.053563      -0.220346
23  0.100000   6.839904  0.051039 -4.499993e-10  0.051319      -0.544447
24  0.100000   8.179247  0.048972 -2.546579e-10  0.049383      -0.831340
25  0.100000   9.780851  0.047206 -2.433610e-10  0.047706      -1.048973
26  0.100000  11.696071  0.045706 -2.289117e-10  0.046287      -1.255947
27  0.100000  13.986316  0.044438 -2.061357e-10  0.045081      -1.425556
28  0.100000  16.725021  0.043373 -1.307333e-10  0.044064      -1.569260
29  0.100000  20.000000  0.042481 -4.404693e-10  0.043211      -1.688722
30  0.001000   4.000000  0.001330 -1.039740e-11  0.001330       0.022644
31  0.001000   4.783253  0.001251  2.227644e-12  0.001250       0.083949
32  0.001000   5.719877  0.001184 -1.421386e-12  0.001183       0.087022
33  0.001000   6.839904  0.001128 -1.325243e-12  0.001126       0.134046
34  0.001000   8.179247  0.001080 -2.238019e-12  0.001078       0.133931
35  0.001000   9.780851  0.001039 -3.464261e-12  0.001036       0.270599
36  0.001000  11.696071  0.001005 -2.260645e-12  0.001005      -0.024660
37  0.001000  13.986316  0.000975 -3.595089e-12  0.000977      -0.168079
38  0.001000  16.725021  0.000950 -3.651491e-12  0.000949       0.080485
39  0.001000  20.000000  0.000928 -2.905575e-12  0.000927       0.078009
 ────────── 
  FL_charm  
 ────────── 
           x         Q2    yadism  yadism_error     APFEL  percent_error
0   0.000100  20.000000  0.000005 -5.098901e-14  0.000005      -0.230613
1   0.000278  20.000000  0.000017  2.898225e-13  0.000017      -0.450185
2   0.000774  20.000000  0.000055 -1.642976e-12  0.000056      -0.444158
3   0.002154  20.000000  0.000174 -1.713579e-12  0.000174      -0.395645
4   0.005995  20.000000  0.000534 -6.902627e-12  0.000537      -0.479797
5   0.016681  20.000000  0.001609 -1.336231e-11  0.001620      -0.701185
6   0.046416  20.000000  0.004668 -3.577503e-11  0.004740      -1.517500
7   0.129155  20.000000  0.012471 -2.107481e-11  0.012829      -2.787246
8   0.359381  20.000000  0.025555 -3.595873e-10  0.028757     -11.135603
9   1.000000  20.000000  0.000000  0.000000e+00  0.000000            NaN
10  0.000100   2.280100  0.000075 -4.328476e-13  0.000075      -0.302434
11  0.000278   2.280100  0.000196 -1.302162e-12  0.000197      -0.229494
12  0.000774   2.280100  0.000510 -6.883614e-13  0.000511      -0.278590
13  0.002154   2.280100  0.001318 -6.521384e-12  0.001321      -0.231433
14  0.005995   2.280100  0.003384 -3.734854e-11  0.003385      -0.035319
15  0.016681   2.280100  0.008610 -2.426286e-11  0.008579       0.364472
16  0.046416   2.280100  0.021500 -7.608356e-11  0.021281       1.028513
17  0.129155   2.280100  0.050964 -2.094306e-10  0.049420       3.124914
18  0.359381   2.280100  0.089296 -1.542112e-09  0.089411      -0.128247
19  1.000000   2.280100  0.000000  0.000000e+00  0.000000            NaN
20  0.100000   4.000000  0.035471  6.754984e-11  0.035258       0.602447
21  0.100000   4.783253  0.031047 -2.244614e-10  0.031010       0.119739
22  0.100000   5.719877  0.027081 -2.392777e-10  0.027160      -0.289815
23  0.100000   6.839904  0.023550 -2.571488e-10  0.023704      -0.653177
24  0.100000   8.179247  0.020424 -1.342561e-10  0.020645      -1.068999
25  0.100000   9.780851  0.017678 -1.174287e-10  0.017936      -1.437156
26  0.100000  11.696071  0.015280 -1.002940e-10  0.015556      -1.770667
27  0.100000  13.986316  0.013202 -8.136334e-11  0.013479      -2.053995
28  0.100000  16.725021  0.011412 -4.613363e-11  0.011680      -2.291399
29  0.100000  20.000000  0.009881 -1.379583e-10  0.010132      -2.475503
30  0.001000   4.000000  0.000453 -7.227605e-12  0.000453      -0.018019
31  0.001000   4.783253  0.000373  8.304087e-13  0.000373       0.035504
32  0.001000   5.719877  0.000305 -8.726530e-13  0.000305       0.196615
33  0.001000   6.839904  0.000249 -7.211695e-13  0.000249       0.301937
34  0.001000   8.179247  0.000203 -1.179184e-12  0.000203       0.169805
35  0.001000   9.780851  0.000165 -1.671592e-12  0.000165       0.028966
36  0.001000  11.696071  0.000134 -9.904684e-13  0.000134      -0.084462
37  0.001000  13.986316  0.000109 -1.419030e-12  0.000110      -0.148850
38  0.001000  16.725021  0.000090 -1.288593e-12  0.000090      -0.254437
39  0.001000  20.000000  0.000074 -9.100151e-13  0.000074      -0.271499

Length mismatch in test when built with debug flags

In order to debug NNPDF/apfelcomb#9, I have changed the build.sh file from
the conda recipe
https://github.com/Zaharid/conda-recipes/tree/master/apfel
to look like:

#!/bin/bash
export CXXFLAGS="-fvisibility-inlines-hidden -std=c++17 -fmessage-length=0 -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -pipe -Og -g -Wall -Wextra -fcheck=all -fbacktrace -fimplicit-none -fvar-tracking-assignments"

export CFLAGS="-march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -pipe -Og -g -Wall -Wextra -fcheck=all -fbacktrace -fimplicit-none -fvar-tracking-assignments"

export FFLAGS="-fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -pipe -fopenmp -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -O2 -pipe -Og -g -Wall -Wextra -fcheck=all -fbacktrace -fvar-tracking-assignments"

./configure --prefix=$PREFIX 
make -j${CPU_COUNT}
make install

so with a bunch of debug flags. This then fails the test as follows:

+ apfel-config --checkapfel
 
 Checking APFEL v3.0.0  ...
At line 8 of file DIS/SetProjectileDIS.f
Fortran runtime error: Actual string length is shorter than the declared one for dummy argument 'lept' (8/12)

Error termination. Backtrace:
#0  0x7f3d76481853 in ???
#1  0x7f3d76482193 in ???
#2  0x7f3d764824c4 in ???
#3  0x7f3d7684a37b in setprojectiledis_
	at DIS/SetProjectileDIS.f:8
#4  0x7f3d768352ae in initparametersdis_
	at DIS/initParametersDIS.f:47
#5  0x7f3d76833785 in initializeapfel_dis_
	at DIS/InitializeAPFEL_DIS.f:25
#6  0x7f3d767b8be9 in checkapfel_
	at Evolution/CheckAPFEL.f:51
#7  0x7f3e2d92f803 in MAIN__
	at /home/zah/miniconda3/conda-bld/apfel_1511985431631/work/examples/CheckAPFEL.f:11
#8  0x7f3e2d92f803 in main
	at /home/zah/miniconda3/conda-bld/apfel_1511985431631/work/examples/CheckAPFEL.f:13

+ exit 0

Also I am using 3.0.0 rather than master.

APFELgui does not compile in version 2.7.0

I get failures when building the APFELgui component of APFEL, namely:

[   96s] src/pdfdialog.cxx: In member function 'void PDFDialog::DIS(double, double, double, double, const string&amp;, const string&amp;, int, const string&amp;, const string&amp;, double*, double*, double*, double*, double*, double*, double*, double*)':
[   96s] src/pdfdialog.cxx:969:9: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]          APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),irep,
[   96s]          ^
[   96s] src/pdfdialog.cxx:980:9: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]          APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),0,
[   96s]          ^
[   96s] src/pdfdialog.cxx:994:13: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]              APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),r+1,
[   96s]              ^
[   96s] src/pdfdialog.cxx:1033:9: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]          APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),0,
[   96s]          ^
[   96s] src/pdfdialog.cxx:1047:13: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]              APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),r+1,
[   96s]              ^
[   96s] src/pdfdialog.cxx:1091:13: error: 'DIS_xsec' is not a member of 'APFEL'
[   96s]              APFEL::DIS_xsec(x,qi,qf,y,0,proc,scheme,pto,pdfset.toStdString(),r+1,
[   96s]              ^
[   97s] Makefile:321: recipe for target 'src/pdfdialog.o' failed
[   97s] make: *** [src/pdfdialog.o] Error 1
[   97s] make: *** Waiting for unfinished jobs....

Would be very grateful if you could take a look at this and let me know if I am doing something wrong or suggest how to fix it. Thanks for your great software

DIS: FONLL b-PDF contributions

We believe to have found a bug in APFEL:

The FONLL prescription uses 4-flavor PDFs in the range between the charm mass and the bottom mass (using the assumtions threshold = mass)
we do not consider the case of intrinsic charm nor intrinsic bottom here
hence the bottom PDF should not contribute to F2light, since there is no 4-flavor bottom (as said, no intrinsic bottom), but we find it is contributing (see below)
this is present for all FONLL flavors
HOWEVER, this bug is not a problem in the fit, because there the bottom PDF only gets a non-zero value beyond the bottom mass and hence will not trigger this bug,
but, it is a bug for anyone else who uses APFEL to compute structure functions and who uses PDFs that where not obtained using the exact same value for the bottom PDF
the bug is most likely triggered by the fact, that APFEL uses evolution basis, which is always the full basis (and so, e.g., Singlet = u+d+s+c+b+t)

To give an explicit example:

below I attach a computation of F2light/FLlight in FONLL-A
mc = 1.95, mb=1e6 so it should be FONLL all the way
the PDF is bonly which is a dummy PDF containing only the b quark as x(1-x)

yadism + APFEL output

Computing for theory=b7b56b0, ocard=21882ed and pdf=bonly ...
LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat
bonly PDF set, member #0, version 1
[15:52:39] INFO     Inspecting runcards...                                               inspector.py:96
           INFO     Inspection completed: success ✓                                     inspector.py:100
           INFO     Interpolation: number of points = 49, polynomial degree = 4,    interpolation.py:440
                    logarithmic = True                                                                  
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01      coupling_constants.py:44
                    5.07600900e-02 1.20409000e-05] [5.07150400e-02                                      
                    9.47604903e-01 1.68100000e-03] [7.43044000e-05                                      
                    1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                                 
                    0.23126, 'MW2': 6463.838403999999}                                                  
           INFO     {'process': 'EM', 'projectilePID': 11, 'polarization': 0,   coupling_constants.py:45
                    'propagatorCorrection': 0}                                                          
           INFO     RenScaleVar: False, FactScaleVar: False                       scale_variations.py:56
           INFO     Thresholds: walls = [0, 3.8024999999999998, 1000000000000.0, 1e+16, thresholds.py:79
                    inf]                                                                                
           INFO     PTO: 1, process: EM                                                    runner.py:146
           INFO     FNS: FONLL-A, NfFF: 4                                                  runner.py:147
           INFO     Intrinsic: []                                                          runner.py:148
           INFO     projectile: electron, target: {Z: 1, A: 1}                             runner.py:149
           INFO     Thresholds: walls = [0, 0.0, 1000000000000.0, 1e+16, inf]           thresholds.py:79
           INFO     Strong Coupling:                                              strong_coupling.py:151
                    a_s(µ_R^2=8317.440000)=0.009390=0.118000/(4π)                                       
                                 ┌────────────────────────────────────┐                                 
                                 │                                    │                                 
                                 │ __     __       _ _                │                                 
                                 │ \ \   / /      | (_)               │                                 
                                 │  \ \_/ /_ _  __| |_ ___ _ __ ___   │                                 
                                 │   \   / _` |/ _` | / __| '_ ` _ \  │                                 
                                 │    | | (_| | (_| | \__ \ | | | | | │                                 
                                 │    |_|\__,_|\__,_|_|___/_| |_| |_| │                                 
                                 │                                    │                                 
                                 └────────────────────────────────────┘                                 
                                              Plan                                                  

• F2_light at 40 pts

• FL_light at 40 pts
                                          Calculation                                               

yadism took off! please stay tuned ...

/home/felix/Physik/N3PDF/DIS/yadism/src/yadism/esf/conv.py:152: IntegrationWarning: The algorithm does

not converge.  Roundoff error is detected

in the extrapolation table.  It is assumed that the requested tolerance

cannot be achieved, and that the returned result (if full_output = 1) is

the best which can be obtained.

res, err = scipy.integrate.quad(

took 4.96 s

Compute external result
WARNING: FONLL-A is a VFN scheme

... setting VFNS PDF evolution

WARNING: FONLL-A is NLO scheme

... setting NLO perturbative order
Welcome to

///    ////   ////   ///_/   _/

_/    /   /    /   /         /         /

////   ////   ///     ///     _/

_/    _/   _/         _/         _/         _/

/    /   /         /         ////   ////

_____v3.0.5 A PDF Evolution Library, arXiv:1310.1394

Authors: V. Bertone, S. Carrazza, J. Rojo
Report of the evolution parameters:
QCD evolution

Space-like evolution (PDFs)

Unpolarized evolution

Evolution scheme: VFNS at N1LO

Solution of the DGLAP equation: 'exactalpha' with maximum 6 active flavours

Solution of the coupling equations: 'exact' with maximum 6 active flavours

Coupling reference value:

AlphaQCD( 91.2000 GeV) =  0.118000

Pole heavy quark masses:
Mc =   1.9500 GeV
Mb = ******** GeV
Mt = ******** GeV

The matching thresholds coincide with the physical masses

muR / muF =  1.0000

Allowed evolution range [   1.0000 :  10000.0000 ] GeV

The internal subgrids will be locked

Fast evolution enabled
Initialization of the evolution completed in   0.536 s
Report of the electroweak parameters:
Mass of the Z = 91.188 GeV

Mass of the W = 80.398 GeV

Mass of the proton = 0.9380 GeV

sin^2(thetaW) = 0.2313

GFermi = 1.16638E-05

| 0.9743 0.2253 0.0035 |

CKM = | 0.2252 0.9735 0.0410 |

| 0.0086 0.0403 0.9992 |

Z propagator correction = 0.00000
Report of the DIS parameters:
Computation in the FONLL-A mass scheme

Electromagnetic (EM) process

Scattering electron - proton

muR / Q =  1.0000

muF / Q =  1.0000

(No scale variation in the evolution)

Target Mass corrections disabled

FONLL damping factor disabled for all heavy quarks

Intrinsic charm disabled
Initialization of the DIS module completed in   5.943 s
LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1

LHAPDF 6.3.0 loading /home/felix/local/share/LHAPDF/bonly/bonly_0000.dat

bonly PDF set, member #0, version 1
Log added, hash=06dd72d0fd24320c82219a81afe221f795fc915f3a14ca711c8da4e2349f05a4
──────────

F2_light

──────────

x          Q2  yadism  yadism_error     APFEL  percent_error

0   0.000100   20.000000     0.0           0.0  0.000053         -100.0

1   0.000275   20.000000     0.0           0.0  0.000123         -100.0

2   0.000756   20.000000     0.0           0.0  0.000285         -100.0

3   0.002080   20.000000     0.0           0.0  0.000655         -100.0

4   0.005721   20.000000     0.0           0.0  0.001494         -100.0

5   0.015733   20.000000     0.0           0.0  0.003383         -100.0

6   0.043267   20.000000     0.0           0.0  0.007593         -100.0

7   0.118993   20.000000     0.0           0.0  0.016630         -100.0

8   0.327252   20.000000     0.0           0.0  0.032793         -100.0

9   0.900000   20.000000     0.0           0.0  0.021900         -100.0

10  0.000100    4.000000     0.0           0.0  0.000068         -100.0

11  0.000275    4.000000     0.0           0.0  0.000155         -100.0

12  0.000756    4.000000     0.0           0.0  0.000351         -100.0

13  0.002080    4.000000     0.0           0.0  0.000784         -100.0

14  0.005721    4.000000     0.0           0.0  0.001720         -100.0

15  0.015733    4.000000     0.0           0.0  0.003715         -100.0

16  0.043267    4.000000     0.0           0.0  0.007880         -100.0

17  0.118993    4.000000     0.0           0.0  0.016282         -100.0

18  0.327252    4.000000     0.0           0.0  0.031179         -100.0

19  0.900000    4.000000     0.0           0.0  0.024752         -100.0

20  0.100000    4.000000     0.0           0.0  0.014415         -100.0

21  0.100000    5.719877     0.0           0.0  0.014465         -100.0

22  0.100000    8.179247     0.0           0.0  0.014507         -100.0

23  0.100000   11.696071     0.0           0.0  0.014543         -100.0

24  0.100000   16.725021     0.0           0.0  0.014573         -100.0

25  0.100000   23.916263     0.0           0.0  0.014599         -100.0

26  0.100000   34.199519     0.0           0.0  0.014623         -100.0

27  0.100000   48.904257     0.0           0.0  0.014643         -100.0

28  0.100000   69.931579     0.0           0.0  0.014661         -100.0

29  0.100000  100.000000     0.0           0.0  0.014678         -100.0

30  0.001000    4.000000     0.0           0.0  0.000439         -100.0

31  0.001000    5.719877     0.0           0.0  0.000416         -100.0

32  0.001000    8.179247     0.0           0.0  0.000396         -100.0

33  0.001000   11.696071     0.0           0.0  0.000380         -100.0

34  0.001000   16.725021     0.0           0.0  0.000366         -100.0

35  0.001000   23.916263     0.0           0.0  0.000353         -100.0

36  0.001000   34.199519     0.0           0.0  0.000343         -100.0

37  0.001000   48.904257     0.0           0.0  0.000333         -100.0

38  0.001000   69.931579     0.0           0.0  0.000324         -100.0

39  0.001000  100.000000     0.0           0.0  0.000317         -100.0
──────────

FL_light

──────────

x          Q2  yadism  yadism_error     APFEL  percent_error

0   0.000100   20.000000     0.0           0.0  0.000002         -100.0

1   0.000275   20.000000     0.0           0.0  0.000005         -100.0

2   0.000756   20.000000     0.0           0.0  0.000013         -100.0

3   0.002080   20.000000     0.0           0.0  0.000035         -100.0

4   0.005721   20.000000     0.0           0.0  0.000094         -100.0

5   0.015733   20.000000     0.0           0.0  0.000246         -100.0

6   0.043267   20.000000     0.0           0.0  0.000605         -100.0

7   0.118993   20.000000     0.0           0.0  0.001272         -100.0

8   0.327252   20.000000     0.0           0.0  0.001712         -100.0

9   0.900000   20.000000     0.0           0.0  0.000078         -100.0

10  0.000100    4.000000     0.0           0.0  0.000002         -100.0

11  0.000275    4.000000     0.0           0.0  0.000007         -100.0

12  0.000756    4.000000     0.0           0.0  0.000018         -100.0

13  0.002080    4.000000     0.0           0.0  0.000049         -100.0

14  0.005721    4.000000     0.0           0.0  0.000133         -100.0

15  0.015733    4.000000     0.0           0.0  0.000348         -100.0

16  0.043267    4.000000     0.0           0.0  0.000855         -100.0

17  0.118993    4.000000     0.0           0.0  0.001799         -100.0

18  0.327252    4.000000     0.0           0.0  0.002421         -100.0

19  0.900000    4.000000     0.0           0.0  0.000111         -100.0

20  0.100000    4.000000     0.0           0.0  0.001613         -100.0

21  0.100000    5.719877     0.0           0.0  0.001475         -100.0

22  0.100000    8.179247     0.0           0.0  0.001360         -100.0

23  0.100000   11.696071     0.0           0.0  0.001262         -100.0

24  0.100000   16.725021     0.0           0.0  0.001179         -100.0

25  0.100000   23.916263     0.0           0.0  0.001106         -100.0

26  0.100000   34.199519     0.0           0.0  0.001042         -100.0

27  0.100000   48.904257     0.0           0.0  0.000985         -100.0

28  0.100000   69.931579     0.0           0.0  0.000934         -100.0

29  0.100000  100.000000     0.0           0.0  0.000889         -100.0

30  0.001000    4.000000     0.0           0.0  0.000024         -100.0

31  0.001000    5.719877     0.0           0.0  0.000022         -100.0

32  0.001000    8.179247     0.0           0.0  0.000020         -100.0

33  0.001000   11.696071     0.0           0.0  0.000019         -100.0

34  0.001000   16.725021     0.0           0.0  0.000017         -100.0

35  0.001000   23.916263     0.0           0.0  0.000016         -100.0

36  0.001000   34.199519     0.0           0.0  0.000015         -100.0

37  0.001000   48.904257     0.0           0.0  0.000015         -100.0

38  0.001000   69.931579     0.0           0.0  0.000014         -100.0

39  0.001000  100.000000     0.0           0.0  0.000013         -100.0

DIS: structure functions in low Q2 regime seems wrong

We believe that APFEL is doing something wrong in the low Q2 regime, i.e., below some quark threshold
in FFNS3 we get a perfect agreement (or as perfect as we can get, when disabling the interpolation)
we believe that in FFNS4 there are no relevant thresholds as masses are either 0 or infinite; however APFEL is sensitive to the quark masses: here mc=2,mb=4
the worring fact is that APFEL is not a continuous function near threshold (but Yadism is) - we thus believe APFEL is wrong
also the ZM-VFNS configuration has the same kind of defect

Below we attach some logs that we created; we also attached the theory input for the first log - in all subsequent logs we only change the indicated FNS

theory 97

{'ID': 22,
 'PTO': 0,
 'FNS': 'FFNS',
 'DAMP': 0,
 'IC': 0,
 'ModEv': 'EXA',
 'XIR': 1.0,
 'XIF': 1.0,
 'NfFF': 3,
 'MaxNfAs': 6,
 'MaxNfPdf': 6,
 'Q0': 1,
 'alphas': 0.11800000000000001,
 'Qref': 91.2,
 'QED': 0,
 'alphaqed': 0.007496251999999999,
 'Qedref': 1.777,
 'SxRes': 0,
 'SxOrd': 'LL',
 'HQ': 'POLE',
 'mc': 2,
 'Qmc': 2,
 'kcThr': 1.0,
 'mb': 4,
 'Qmb': 4,
 'kbThr': 1.0,
 'mt': 173.07,
 'Qmt': 173.07,
 'ktThr': 1.0,
 'CKM': '0.97428 0.22530 0.003470 0.22520 0.97345 0.041000 0.00862 0.04030 0.999152',
 'MZ': 91.1876,
 'MW': 80.398,
 'GF': 1.1663787e-05,
 'SIN2TW': 0.23126,
 'TMC': 0,
 'MP': 0.938,
 'Comments': 'LO baseline for small-x res',
 'global_nx': 0,
 'EScaleVar': 1,
 '_modify_time': '2020-05-28 16:55:02.150509'}

FFNS3

F2total with theory=97 using CT14llo_NF6
           x          Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000100   90.000000  2.243980  2.243980           0.0   -0.000004
1   0.000317   90.000000  1.639572  1.639572           0.0   -0.000004
2   0.001007   90.000000  1.182349  1.182349           0.0   -0.000004
3   0.003195   90.000000  0.846563  0.846563           0.0   -0.000004
4   0.010138   90.000000  0.623498  0.623498           0.0   -0.000004
5   0.032169   90.000000  0.483802  0.483802           0.0   -0.000004
6   0.102077   90.000000  0.370392  0.370392           0.0   -0.000004
7   0.304545   90.000000  0.225637  0.225637           0.0   -0.000003
8   0.536364   90.000000  0.072848  0.072848           0.0   -0.000003
9   0.768182   90.000000  0.007164  0.007164           0.0   -0.000003
10  1.000000   90.000000  0.000000  0.000000           0.0         NaN
11  0.845455   90.000000  0.001729  0.001729           0.0   -0.000003
12  0.922727   90.000000  0.000148  0.000148           0.0   -0.000003
13  1.000000   90.000000  0.000000  0.000000           0.0         NaN
14  0.010000    2.000000  0.420156  0.420175           0.0    0.004658
15  0.010000    3.088904  0.446229  0.446251           0.0    0.004783
16  0.010000    4.770665  0.472257  0.472281           0.0    0.005041
17  0.010000    7.368063  0.497732  0.497759           0.0    0.005401
18  0.010000   11.379620  0.522376  0.522406           0.0    0.005827
19  0.010000   17.575279  0.546045  0.546080           0.0    0.006287
20  0.010000   27.144176  0.568668  0.568707           0.0    0.006788
21  0.010000   41.922880  0.590220  0.590263           0.0    0.007298
22  0.010000   64.747880  0.610707  0.610755           0.0    0.007827
23  0.010000  100.000000  0.630162  0.630215           0.0    0.008352

FFNS4

F2total with theory=99 using CT14llo_NF6
           x          Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000100   90.000000  3.231555  3.231554           0.0   -0.000004
1   0.000317   90.000000  2.317068  2.317068           0.0   -0.000004
2   0.001007   90.000000  1.627903  1.627903           0.0   -0.000004
3   0.003195   90.000000  1.119541  1.119541           0.0   -0.000004
4   0.010138   90.000000  0.771808  0.771808           0.0   -0.000004
5   0.032169   90.000000  0.550866  0.550866           0.0   -0.000004
6   0.102077   90.000000  0.394615  0.394615           0.0   -0.000004
7   0.304545   90.000000  0.230012  0.230012           0.0   -0.000004
8   0.536364   90.000000  0.073214  0.073214           0.0   -0.000004
9   0.768182   90.000000  0.007169  0.007169           0.0   -0.000004
10  1.000000   90.000000  0.000000  0.000000           0.0         NaN
11  0.845455   90.000000  0.001729  0.001729           0.0   -0.000004
12  0.922727   90.000000  0.000148  0.000148           0.0   -0.000004
13  1.000000   90.000000  0.000000  0.000000           0.0         NaN
14  0.010000    2.000000  0.493705  0.426851           0.0  -13.541340
15  0.010000    3.088904  0.553139  0.471357           0.0  -14.785185
16  0.010000    4.770665  0.515436  0.515463           0.0    0.005241
17  0.010000    7.368063  0.558580  0.558612           0.0    0.005756
18  0.010000   11.379620  0.600308  0.600346           0.0    0.006350
19  0.010000   17.575279  0.640385  0.640430           0.0    0.006980
20  0.010000   27.144176  0.678697  0.678749           0.0    0.007633
21  0.010000   41.922880  0.715204  0.715264           0.0    0.008301
22  0.010000   64.747880  0.749925  0.749992           0.0    0.008971
23  0.010000  100.000000  0.782909  0.782985           0.0    0.009625

FFNS5

F2total with theory=101 using CT14llo_NF6
           x          Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000100   90.000000  3.349632  3.349632           0.0   -0.000004
1   0.000317   90.000000  2.392680  2.392680           0.0   -0.000004
2   0.001007   90.000000  1.673884  1.673884           0.0   -0.000004
3   0.003195   90.000000  1.145424  1.145424           0.0   -0.000004
4   0.010138   90.000000  0.784716  0.784716           0.0   -0.000004
5   0.032169   90.000000  0.556219  0.556219           0.0   -0.000004
6   0.102077   90.000000  0.396307  0.396307           0.0   -0.000004
7   0.304545   90.000000  0.230248  0.230248           0.0   -0.000004
8   0.536364   90.000000  0.073230  0.073230           0.0   -0.000004
9   0.768182   90.000000  0.007169  0.007169           0.0   -0.000004
10  1.000000   90.000000  0.000000  0.000000           0.0         NaN
11  0.845455   90.000000  0.001729  0.001729           0.0   -0.000004
12  0.922727   90.000000  0.000148  0.000148           0.0   -0.000004
13  1.000000   90.000000  0.000000  0.000000           0.0         NaN
14  0.010000    2.000000  0.519120  0.426851           0.0  -17.774098
15  0.010000    3.088904  0.580841  0.471357           0.0  -18.849231
16  0.010000    4.770665  0.526678  0.515463           0.0   -2.129342
17  0.010000    7.368063  0.570659  0.558612           0.0   -2.110929
18  0.010000   11.379620  0.613196  0.600346           0.0   -2.095608
19  0.010000   17.575279  0.640385  0.640430           0.0    0.006980
20  0.010000   27.144176  0.680323  0.680375           0.0    0.007657
21  0.010000   41.922880  0.721307  0.721368           0.0    0.008385
22  0.010000   64.747880  0.760046  0.760115           0.0    0.009106
23  0.010000  100.000000  0.796843  0.796921           0.0    0.009806

ZM-VFNS

F2total with theory=102 using CT14llo_NF6
           x          Q2     APFEL    yadism  yadism_error  rel_err[%]
0   0.000100   90.000000  3.349632  3.349632           0.0   -0.000004
1   0.000317   90.000000  2.392680  2.392680           0.0   -0.000004
2   0.001007   90.000000  1.673884  1.673884           0.0   -0.000004
3   0.003195   90.000000  1.145424  1.145424           0.0   -0.000004
4   0.010138   90.000000  0.784716  0.784716           0.0   -0.000004
5   0.032169   90.000000  0.556219  0.556219           0.0   -0.000004
6   0.102077   90.000000  0.396307  0.396307           0.0   -0.000004
7   0.304545   90.000000  0.230248  0.230248           0.0   -0.000004
8   0.536364   90.000000  0.073230  0.073230           0.0   -0.000004
9   0.768182   90.000000  0.007169  0.007169           0.0   -0.000004
10  1.000000   90.000000  0.000000  0.000000           0.0         NaN
11  0.845455   90.000000  0.001729  0.001729           0.0   -0.000004
12  0.922727   90.000000  0.000148  0.000148           0.0   -0.000004
13  1.000000   90.000000  0.000000  0.000000           0.0         NaN
14  0.010000    2.000000  0.423493  0.420175           0.0   -0.783537
15  0.010000    3.088904  0.458781  0.446251           0.0   -2.731341
16  0.010000    4.770665  0.515436  0.515463           0.0    0.005241
17  0.010000    7.368063  0.558580  0.558612           0.0    0.005756
18  0.010000   11.379620  0.600308  0.600346           0.0    0.006350
19  0.010000   17.575279  0.640385  0.640430           0.0    0.006980
20  0.010000   27.144176  0.680323  0.680375           0.0    0.007657
21  0.010000   41.922880  0.721307  0.721368           0.0    0.008385
22  0.010000   64.747880  0.760046  0.760115           0.0    0.009106
23  0.010000  100.000000  0.796843  0.796921           0.0    0.009806

DIS: IC FONLL-A Scale variations

We believe that the order of the calls to IncludeScaleVariation and IncludeIntrinsicCharm has to be switched here:

apfel/src/DIS/initIntegralsDIS.f

Lines 52 to 60 in fdfd068

    
           * 
        
           *     Include intrinsic charm contributions 
        
           * 
        
                 if(IntrinsicCharm) call IncludeIntrinsicCharm 
        
           * 
        
           *     Scale variations 
        
           *     (only if the dynamical scale variation is disabled) 
        
           * 
        
                 if(.not.DynScVar) call IncludeScaleVariation

otherwise the assignment here

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 293 to 295 in fdfd068

    
           *     Neutral Current 
        
                                SC2m0NC(igrid,ixi,3,1,beta,alpha) = 
        
                1                    SC2zm(igrid,Nf_FF,3,1,beta,alpha)

comes too early and misses the scale variations of the light coefficient functions that are however needed to cancel the scale variation from the intrinsic counter part - compare also Appendix B of https://inspirehep.net/literature/1397009

DIS: IC FONLL-A q Matching factorization scale log

We believe the definitions of the factorization scale log in the FONLL matching coefficient functions contain a typo:

apfel/src/DIS/MassiveZeroCoefficientFunctions.f

Lines 1129 to 1131 in fdfd068

    
                 double precision DICa 
        
           * 
        
                 lnk = dlog(xi/kfacQ)                ! ln(muF2/mh2)

apfel/src/DIS/MassiveZeroCoefficientFunctions.f

Lines 1157 to 1159 in fdfd068

    
                 double precision DICb 
        
           * 
        
                 lnk = dlog(xi/kfacQ)                ! ln(muF2/mh2)

apfel/src/DIS/MassiveZeroCoefficientFunctions.f

Lines 1185 to 1187 in fdfd068

    
                 double precision DICc 
        
           * 
        
                 lnk = dlog(xi/kfacQ)                ! ln(muF2/mh2)

the correct thing to do is a product of xi and kfacQ as is done for the gluon case

apfel/src/DIS/IncludeIntrinsicCharm.f

Lines 59 to 61 in fdfd068

    
           xi     = xigrid(ixi*xistep) 
        
           lambda = 1d0 / xi 
        
           lnF    = dlog(xi*kfacQ)

see also the report of that variable:

apfel/src/DIS/ReportParametersDIS.f

Line 81 in fdfd068

write(6,"(a,f7.4)") " muF / Q = ",dsqrt(kfacQ)

EVOL: Top msbar mass running

We believe the routine computing the top mass evolution in the msbar scheme is not working properly in case more than one flavour patch is present.
In particular in file MSbarmass.f:

           ELSEIF(Q2.GE.M2TH(4).AND.NFMAXALPHA.GE.4)THEN
               EVF = DECOUP("DW",5,DLOG(K2TH(5)))
     1              / EVMASS(4,ASQ,ASIM)
     2              * DECOUP("DW",6,DLOG(K2TH(6)))
     3              / EVMASS(5,ASTH(5),ASTHM(6))
            ELSE
               EVF = DECOUP("DW",4,DLOG(K2TH(4)))
     1              / EVMASS(3,ASQ,ASIM)
     2              * DECOUP("DW",5,DLOG(K2TH(5)))
     3              / EVMASS(4,ASTH(4),ASTHM(5))
     4              * DECOUP("DW",6,DLOG(K2TH(6)))
     5              / EVMASS(5,ASTH(5),ASTHM(6))

we propose the following fix:

             ELSEIF(Q2.GE.M2TH(4).AND.NFMAXALPHA.GE.4)THEN
               EVF = DECOUP("DW",5,DLOG(K2TH(5)))
     1              / EVMASS(4,ASQ,ASTH(5))
     2              * DECOUP("DW",6,DLOG(K2TH(6)))
     3              / EVMASS(5,ASTH(5),ASIM)
            ELSE
               EVF = DECOUP("DW",4,DLOG(K2TH(4)))
     1              / EVMASS(3,ASQ,ASTH(4))
     2              * DECOUP("DW",5,DLOG(K2TH(5)))
     3              / EVMASS(4,ASTH(4),ASTHM(5))
     4              * DECOUP("DW",6,DLOG(K2TH(6)))
     5              / EVMASS(5,ASTH(5),ASIM)

We found good agreement with the Eko implementation (WIP) and the fixed apfel routine.

add topics

I suggest adding topics such as high-energy-physics, particle-physics in the About section at https://github.com/scarrazza/apfel.

	FLlight = FLlight + w_int_gen(n,alpha,xi)
	1 * ( FL(3,0,alpha)
	2 + c1 * F2(3,0,alpha) + c2 * I2(3,0,alpha) )

	SC2mNC(igrid,ixi,3,0,beta,alpha) = ( 2d0 - eta ) * fL
	SCLmNC(igrid,ixi,3,0,beta,alpha) = 4d0 * ( 1d0 - eta )
	1 / ( 2d0 - eta ) / factL * fL
	*
	* NLO
	*
	if(ipt.ge.1)then
	*
	bound = alpha-inter_degree(igrid)
	if(alpha.lt.inter_degree(igrid)) bound = 0
	*
	c = max(xg(igrid,beta),
	1 xg(igrid,beta)/xg(igrid,alpha+1)) / eta
	d = min(eta,xg(igrid,beta)/xg(igrid,bound)) / eta
	*
	if(c.ge.1d0) cycle
	*
	walpha = alpha
	wbeta = beta
	*
	k = 3
	wl = 1
	*
	sf = 1
	C2RS = fact2 * eta * dgauss(integrandsICm,c,d,eps)
	C2L = fact2 * eta * c21ICL(c)
	*
	sf = 2
	CLRS = eta * dgauss(integrandsICm,c,d,eps)
	CLL = eta * cL1ICL(c)
	*
	if(MassScheme(1:5).eq."FONLL")then
	sf = 2
	* Gluon
	k = 1
	gDIC = dgauss(integrandsICm,c,d,eps)
	* Put these terms in the LO gluon slot (e.g. SC2mNC(igrid,ixi,1,0,beta,alpha))
	* which is empty.
	SC2mNC(igrid,ixi,1,0,beta,alpha) =
	1 - ( 2d0 - eta ) * lnF * gDIC
	SCLmNC(igrid,ixi,1,0,beta,alpha) =
	1 - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta ) / factL
	2 * lnF * gDIC
	* Non-singlet
	k = 3
	wl = 2
	nsDIC = dgauss(integrandsICm,c,d,eps)
	*
	C2RS = C2RS - ( 2d0 - eta ) * nsDIC
	C2L = C2L - ( 2d0 - eta ) * DICc(xi,c)
	*
	CLRS = CLRS - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta )
	1 * nsDIC
	CLL = CLL - 4d0 * ( 1d0 - eta ) / ( 2d0 - eta )
	1 * DICc(xi,c)
	endif
	*
	SC2mNC(igrid,ixi,3,1,beta,alpha) = C2RS + C2L * fL
	SCLmNC(igrid,ixi,3,1,beta,alpha) = CLRS + CLL * fL

	if(IntrinsicCharm)then
	if(Nf_FF.lt.4)then
	do pt=0,ipt_max_IC
	CLnsp(4) = CLnsp(4) + as(pt)
	1 * ( c0(4)
	2 * SCLm0NC(jgrid,ixi(4),
	3 3,pt,alpha,beta)
	4 + c1(4)
	5 * SCLm0NC(jgrid,ixi(4)+1,
	6 3,pt,alpha,beta) )
	C2nsp(4) = C2nsp(4) + as(pt)
	1 * ( c0(4)
	2 * SC2m0NC(jgrid,ixi(4),
	3 3,pt,alpha,beta)
	4 + c1(4)
	5 * SC2m0NC(jgrid,ixi(4)+1,
	6 3,pt,alpha,beta) )
	enddo
	endif
	endif

	*
	* Include intrinsic charm contributions
	*
	if(IntrinsicCharm) call IncludeIntrinsicCharm
	*
	* Scale variations
	* (only if the dynamical scale variation is disabled)
	*
	if(.not.DynScVar) call IncludeScaleVariation

	* Neutral Current
	SC2m0NC(igrid,ixi,3,1,beta,alpha) =
	1 SC2zm(igrid,Nf_FF,3,1,beta,alpha)

	xi = xigrid(ixi*xistep)
	lambda = 1d0 / xi
	lnF = dlog(xi*kfacQ)

	c1 = x2 / xi2 / tau
	c2 = 4d0 * rhop * x3 / tau1.5d0

scarrazza / apfel Goto Github PK

apfel's People

Contributors

Stargazers

Watchers

Forkers

apfel's Issues

Recommend Projects

Recommend Topics

Recommend Org