Computing observables in the FONLL scheme

The actual implementation of FONLL observables requires the computation of the partonic coefficients functions using different FFNS settings with yadism and the possibility to compute coexisting flavor number PDFs, which is provided by the EKO library. In this tutorial we explain, using an example, how observables can be calculated using a fixed input PDF by evolving it to the three different flavor number schemes.

As a first step we define the define the inputs needed for a computation of a DIS observable. In the following example we will compute predictions for \(F_{2c}\).

[82]:
from eko.interpolation import lambertgrid

# Define a DIS observable runcard:
observablecard = {
    # Process type: "EM", "NC", "CC"
    "prDIS": "NC",
    # Projectile: "electron", "positron", "neutrino", "antineutrino"
    "ProjectileDIS": "electron",
    # Scattering target: "proton", "neutron", "isoscalar", "lead", "iron", "neon" or "marble"
    "TargetDIS": "proton",
    # Interpolation: if True use log interpolation
    "interpolation_is_log": True,
    # Interpolation: polynomial degree, 1 = linear, ...
    "interpolation_polynomial_degree": 4,
    # Interpolation: xgrid values
    # here we use an equally log-linear spaced grid
    "interpolation_xgrid": lambertgrid(30),
    # Observables configurations
    "observables": {
      "F2_charm": [
        {
            "y": 0.,
            "x": 1e-4,
            "Q2": 15.0**2,
        },
        {
            "y": 0.,
            "x": 1e-3,
            "Q2": 15.0**2,
        },
        # Add here the kinematics of other datapoints
      ],
      # Potentially include observables other than F2_charm,
      # each of them has to be: TYPE_heaviness, where heaviness can take:
      # "charm", "bottom", "total" or "light".
    },
    # Projectile polarization faction, float from 0 to 1.
    "PolarizationDIS": 0.0,
    # Exchanged boson propagator correction
    "PropagatorCorrection": 0.0,
    # Restrict boson coupling to a single parton ? Monte Carlo PID or None for all partons
    "NCPositivityCharge": None,
}

# Define the theory parameters:
theorycard = {
    "CKM": "0.97428 0.22530 0.003470 0.22520 0.97345 0.041000 0.00862 0.04030 0.999152", # CKM matrix elements
    "GF": 1.1663787e-05, # [GeV^-2] Fermi coupling constant
    "MP": 0.938, # [GeV] proton mass
    "MW": 80.398, # [GeV] W boson mass
    "MZ": 91.1876, # [GeV] Z boson mass
    "mc": 1.51, # [GeV] charm mass
    "mb": 4.92, # [GeV] bottom mass
    "mt": 172.5, # [GeV] top mass
    "kcThr": 1.0, # ratio of the charm matching scale over the charm mass
    "kbThr": 1.0, # ratio of the bottom matching scale over the bottom mass
    "ktThr": 1.0, # ratio of the top matching scale over the top mass
    "alphas": 0.118, # alphas value at the reference scale
    "Q0": 1.65, # [GeV] reference scale for the flavor patch determination
    "nf0": 4, # number of active flavors at the Q0 reference scale
    "PTO": 1, # perturbative order in alpha_s: 0 = LO (alpha_s^0), 1 = NLO (alpha_s^1) ...
    "alphaqed": 0.007496252, # alpha_em value
    "Qref": 91.2, # [GeV] reference scale for the alphas value
    "nfref": 5, # number of active flavors at the reference scale Qref
    "QED": 0, # QED correction to running of strong coupling: 0 = disabled, 1 = allowed
    "XIF": 1.0, # ratio of factorization scale over the hard scattering scale
    "XIR": 1.0, # ratio of renormalization scale over the hard scattering scale
    "TMC": 1, # include target mass corrections: 0 = disabled, 1 = leading twist, 2 = higher twist approximated, 3 = higher twist exact
    "n3lo_cf_variation": 0, # N3LO coefficient functions variation: -1 = lower bound, 0 = central , 1 = upper bound
    "MaxNfAs": 5, # maximum number of flavors in running of strong coupling
    "HQ": "POLE", # heavy quark mass scheme (only for DGLAP evolution)
    "MaxNfPdf": 5, # maximum number of flavors in running of PDFs (needed by pineko)
    "IC": 1, # 0 = perturbative charm only, 1 = intrinsic charm allowed
    "ModEv": "EXA",  # alphas evolution method, "EXA" = "exact", "TRN" = "truncated"
}

# Define the PDF evolution grid and other DGLAP evolution parameters
operatorcard = {
    # Evolution starting scale
    'mu0': theorycard["Q0"],
    # Evolution final scales: add here all the points needed in the evolved PDF
    # with the corresponding active number of flavor (Q_i, nf_i)
    # (we will update this entry later)
    'mugrid': [(10.0, 5),(15.0, 5),(20.0, 5),],
    # specify an interpolation grid
    'xgrid': lambertgrid(30),
    # DGLAP evolution settings
    'configs': {
        # evolution method: 'iterate-exact', 'iterate-expanded', 'truncated' ... see https://eko.readthedocs.io/en/latest/theory/DGLAP.html
        'evolution_method': 'iterate-exact',
        # other evolution parameters. See https://eko.readthedocs.io/en/latest/theory/DGLAP.html
        'ev_op_max_order': [10, 0],
        'ev_op_iterations': 10,
        # Interpolation: polynomial degree, 1 = linear, ...
        'interpolation_polynomial_degree': 4,
        # Interpolation: if True use log interpolation
        'interpolation_is_log': True,
        # Scale variations: "None", "expanded", "exponentiated"
        'scvar_method': None,
        # Backward evolution matching inversion: "None", "expanded", "exact"
        'inversion_method': None,
        # Number of parallel cores to use: 0 = all, 1 = one, -1 = all but one ...
        'n_integration_cores': 1,
        # Polarized evolution ?
        'polarized': False,
        # Time like evolution ?
        'time_like': False
    },
    # Debug options
    'debug': {
        'skip_singlet': False,
        'skip_non_singlet': False
    }
}

Then we compute the necessary coexisting flavor number PDFs starting from a given PDF defined for some value of \(Q\) and corresponding \(n_f\) using EKO (for more detailed information on the use of EKO we refer to the documentation ). Here, we show how to evolve NNPDF40_nnlo_as_01180 to produce three PDFs in FFNS3, FFNS4, and FFNS5 respectively.

[83]:
import lhapdf
import pathlib
from eko.io import runcards
from ekobox.evol_pdf import evolve_pdfs


# Define the PDF for which we want to make the prediction
# Here we take the replica 0 of NNPDF4.0 - remember to have it available via LHAPDF.
# You can install it if necessary via `lhapdf get NNPDF40_nnlo_as_01180`.
input_pdf = lhapdf.mkPDF("NNPDF40_nlo_as_01180", 0)

# Build eko-like runcards objects
theory_card = runcards.Legacy(theorycard, {}).new_theory
operators_card = runcards.OperatorCard.from_dict(operatorcard)

# Iterate the flavor numbers
for nf in [3, 4, 5]:
    # Update the nf value in the mugrid
    operators_card.mugrid = [(mu, nf) for (mu,_) in operators_card.mugrid]

    eko_name = f"eko_nf{nf}.tar"
    pathlib.Path(eko_name).unlink(missing_ok=True)
    # Evolve the LHAPDF grid and store the output
    evolve_pdfs(
        initial_PDF_list=[input_pdf],
        theory_card=theory_card,
        operators_card=operators_card,
        store_path=eko_name,  # where to store the eko
        install=True, # Install in the evolved PDF in the LHAPDF directory
        name=f"NNPDF40_evolved_nf{nf}",  # name of the evolved pdf
    )
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_nlo_as_01180/NNPDF40_nlo_as_01180_0000.dat
NNPDF40_nlo_as_01180 PDF set, member #0, version 1; LHAPDF ID = 331700
[15:30:11] INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=2.7224999999999997 matchings.py:72
                    @ 4                                                                                            
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     {'polarized': False}                                                                 grid.py:93
           INFO     {'time_like': False}                                                                 grid.py:94
           INFO     Prepare threshold operator                                                          grid.py:139
           INFO     Evolution: computing operators 2.722500e+00 -> 2.280100e+00, nf=4               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.722500e+00 -> 2.280100e+00                         __init__.py:924
           INFO     Evolution: a_s distance: 2.605852e-02 -> 2.728162e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:30:20] INFO     Evolution: Total time 8.296076 s                                                __init__.py:982
           INFO     Matching: order: (1, 0), backward method: None                   operator_matrix_element.py:327
[15:30:31] INFO     Matching: Total time 11.472626 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.280100e+00 -> 1.000000e+02, nf=3               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.280100e+00 -> 1.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.728162e-02 -> 1.330271e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:30:39] INFO     Evolution: Total time 8.060058 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.280100e+00 -> 2.250000e+02, nf=3               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.280100e+00 -> 2.250000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.728162e-02 -> 1.202975e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:30:47] INFO     Evolution: Total time 7.896101 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.280100e+00 -> 4.000000e+02, nf=3               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.280100e+00 -> 4.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.728162e-02 -> 1.126961e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:30:55] INFO     Evolution: Total time 7.954777 s                                                __init__.py:982
[15:30:56] INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
Overwriting old PDF installation
install_pdf NNPDF40_evolved_nf3
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=2.7224999999999997 matchings.py:72
                    @ 4                                                                                            
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     {'polarized': False}                                                                 grid.py:93
           INFO     {'time_like': False}                                                                 grid.py:94
           INFO     Evolution: computing operators 2.722500e+00 -> 1.000000e+02, nf=4               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.722500e+00 -> 1.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.605852e-02 -> 1.391613e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:04] INFO     Evolution: Total time 7.748763 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.722500e+00 -> 2.250000e+02, nf=4               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.722500e+00 -> 2.250000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.605852e-02 -> 1.263129e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:12] INFO     Evolution: Total time 7.965763 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.722500e+00 -> 4.000000e+02, nf=4               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.722500e+00 -> 4.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 2.605852e-02 -> 1.185893e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:20] INFO     Evolution: Total time 7.973417 s                                                __init__.py:982
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
Overwriting old PDF installation
install_pdf NNPDF40_evolved_nf4
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=2.7224999999999997 matchings.py:72
                    @ 4                                                                                            
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     {'polarized': False}                                                                 grid.py:93
           INFO     {'time_like': False}                                                                 grid.py:94
           INFO     Prepare threshold operator                                                          grid.py:139
           INFO     Evolution: computing operators 2.722500e+00 -> 2.420640e+01, nf=4               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.722500e+00 -> 2.420640e+01                         __init__.py:924
           INFO     Evolution: a_s distance: 2.605852e-02 -> 1.697192e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:28] INFO     Evolution: Total time 7.918660 s                                                __init__.py:982
           INFO     Matching: order: (1, 0), backward method: None                   operator_matrix_element.py:327
[15:31:41] INFO     Matching: Total time 13.032816 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.420640e+01 -> 1.000000e+02, nf=5               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.420640e+01 -> 1.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 1.697192e-02 -> 1.415523e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:49] INFO     Evolution: Total time 8.292476 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.420640e+01 -> 2.250000e+02, nf=5               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.420640e+01 -> 2.250000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 1.697192e-02 -> 1.293886e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:31:58] INFO     Evolution: Total time 8.344759 s                                                __init__.py:982
           INFO     Evolution: computing operators 2.420640e+01 -> 4.000000e+02, nf=5               __init__.py:917
           INFO     Evolution: µ_R^2 distance: 2.420640e+01 -> 4.000000e+02                         __init__.py:924
           INFO     Evolution: a_s distance: 1.697192e-02 -> 1.219852e-02                           __init__.py:936
           INFO     Evolution: order: (2, 0), solution strategy: iterate-exact, use fhmruvv: False  __init__.py:946
[15:32:05] INFO     Evolution: Total time 7.532989 s                                                __init__.py:982
[15:32:06] INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
Overwriting old PDF installation
install_pdf NNPDF40_evolved_nf5

Next, we compute the coefficient functions using yadism and convolve them with the three PDFs generated in the script above. The example script works for FONLL-A and FONLL-C without damping. FONLL-B requires more care as one has to compute the same coefficient functions but at different perturbative orders. Also note that the script below just combines various FFNSs, but does not produce a full VFNS. For example, one might want to introduce here a dumping function to turn off e.g. \(n_f=4\) contributions below the charm matching scale.

[84]:
import yadism

# Define the settings for the different contributions to the FONLL observable
settings_list = [
    ["FONLL-FFNS", 3, "full", "NNPDF40_evolved_nf3"],
    ["FONLL-FFN0", 3, "full", "NNPDF40_evolved_nf3"],
    ["FONLL-FFNS", 4, "massless", "NNPDF40_evolved_nf4"],
    ["FONLL-FFN0", 4, "full", "NNPDF40_evolved_nf4"],
    ["FONLL-FFNS", 5, "full", "NNPDF40_evolved_nf5"],
]

single_scheme_result = []
# Iterate parts
for fns, nf, fonllparts, pdfname in settings_list:
    # Replace the relevant settings in the theory card
    theorycard["FNS"] = fns
    theorycard["NfFF"] = nf
    theorycard["FONLLParts"] = fonllparts
    pdf = lhapdf.mkPDF(pdfname)

    # Compute the observable in each of the schemes defined in settings_list
    out = yadism.run_yadism(theorycard, observablecard)

    # Store the results in single_scheme_result
    single_scheme_result.append(out.apply_pdf(pdf))

# Extract the different terms form single_scheme_result and combine them to
# obtain the coexisting flavor number scheme FONLL result
result = {}
for observable_name, kinematic_points in single_scheme_result[0].items():
    result[observable_name] = []
    for i in range(len(kinematic_points)):
        result[observable_name].append(
            single_scheme_result[0][observable_name][i]["result"]  # FONLL-FFNS3
            - single_scheme_result[1][observable_name][i]["result"]  # FONLL-FFN03
            + single_scheme_result[2][observable_name][i]["result"]  # FONLL-FFNS4
            - single_scheme_result[3][observable_name][i]["result"]  # FONLL-FFN04
            + single_scheme_result[4][observable_name][i]["result"]  # FONLL-FFNS5
        )
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_evolved_nf3/NNPDF40_evolved_nf3_0000.dat
NNPDF40_evolved_nf3 PDF set, member #0, version 1
[15:32:10] INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01 5.07600900e-02  coupling_constants.py:53
                    1.20409000e-05] [5.07150400e-02 9.47604903e-01 1.68100000e-03]                                 
                    [7.43044000e-05 1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                            
                    0.23121, 'MW2': 6463.838403999999}                                                             
           INFO     {'process': 'NC', 'projectilePID': 11, 'polarization': 0.0,            coupling_constants.py:54
                    'propagatorCorrection': 0.0, 'nc_pos_charge': None}                                            
           INFO     RenScaleVar: True, FactScaleVar: True                                    scale_variations.py:62
           INFO     Atlas [0.00e+00 - inf - inf - inf - inf], ref=2.7224999999999997 @ 4            matchings.py:72
           INFO     PTO: 1, PTO@evolution: 1, process: NC                                             runner.py:142
           INFO     FNS: FONLL-FFNS, NfFF: 3                                                          runner.py:149
           INFO     projectile: electron, target: {Z: 1, A: 1}                                        runner.py:150
                                      ┌────────────────────────────────────┐
                                      │                                    │
                                      │ __     __       _ _                │
                                      │ \ \   / /      | (_)               │
                                      │  \ \_/ /_ _  __| |_ ___ _ __ ___   │
                                      │   \   / _` |/ _` | / __| '_ ` _ \  │
                                      │    | | (_| | (_| | \__ \ | | | | | │
                                      │    |_|\__,_|\__,_|_|___/_| |_| |_| │
                                      │                                    │
                                      └────────────────────────────────────┘
                                                       Plan
F2_charm at 2 pts
                                                    Calculation
yadism took off! please stay tuned ...
[15:32:10] INFO     computing P_qq_0 - took: 0.161964 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”>[15:32:10]</span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qq_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.161964</span> s ]8;id=580267;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=808659;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{[15:32:10]}textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qq_0 - took: textcolor{ansi-cyan-intense}{textbf{0.161964}} s ]8;id=580267;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=808659;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

[15:32:10] INFO  computing P_qq_0 - took: 0.161964 s ]8;id=580267;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=808659;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–

           INFO     computing P_qg_0 - took: 0.134157 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”> </span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qg_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.134157</span> s ]8;id=209988;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=643495;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{ }textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qg_0 - took: textcolor{ansi-cyan-intense}{textbf{0.134157}} s ]8;id=209988;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=643495;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

  INFO  computing P_qg_0 - took: 0.134157 s ]8;id=209988;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=643495;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–



took 2.99 s
[15:32:13] INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_evolved_nf3/NNPDF40_evolved_nf3_0000.dat
NNPDF40_evolved_nf3 PDF set, member #0, version 1
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01 5.07600900e-02  coupling_constants.py:53
                    1.20409000e-05] [5.07150400e-02 9.47604903e-01 1.68100000e-03]                                 
                    [7.43044000e-05 1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                            
                    0.23121, 'MW2': 6463.838403999999}                                                             
           INFO     {'process': 'NC', 'projectilePID': 11, 'polarization': 0.0,            coupling_constants.py:54
                    'propagatorCorrection': 0.0, 'nc_pos_charge': None}                                            
           INFO     RenScaleVar: True, FactScaleVar: True                                    scale_variations.py:62
           INFO     Atlas [0.00e+00 - inf - inf - inf - inf], ref=2.7224999999999997 @ 4            matchings.py:72
           INFO     PTO: 1, PTO@evolution: 1, process: NC                                             runner.py:142
           INFO     FNS: FONLL-FFN0, NfFF: 3                                                          runner.py:149
           INFO     projectile: electron, target: {Z: 1, A: 1}                                        runner.py:150
                                      ┌────────────────────────────────────┐
                                      │                                    │
                                      │ __     __       _ _                │
                                      │ \ \   / /      | (_)               │
                                      │  \ \_/ /_ _  __| |_ ___ _ __ ___   │
                                      │   \   / _` |/ _` | / __| '_ ` _ \  │
                                      │    | | (_| | (_| | \__ \ | | | | | │
                                      │    |_|\__,_|\__,_|_|___/_| |_| |_| │
                                      │                                    │
                                      └────────────────────────────────────┘
                                                       Plan
F2_charm at 2 pts
                                                    Calculation
yadism took off! please stay tuned ...
[15:32:13] INFO     computing P_qq_0 - took: 0.165385 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”>[15:32:13]</span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qq_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.165385</span> s ]8;id=470632;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=159092;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{[15:32:13]}textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qq_0 - took: textcolor{ansi-cyan-intense}{textbf{0.165385}} s ]8;id=470632;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=159092;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

[15:32:13] INFO  computing P_qq_0 - took: 0.165385 s ]8;id=470632;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=159092;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–

           INFO     computing P_qg_0 - took: 0.134486 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”> </span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qg_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.134486</span> s ]8;id=259115;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=905741;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{ }textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qg_0 - took: textcolor{ansi-cyan-intense}{textbf{0.134486}} s ]8;id=259115;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=905741;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

  INFO  computing P_qg_0 - took: 0.134486 s ]8;id=259115;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=905741;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–



took 1.57 s
[15:32:14] INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_evolved_nf4/NNPDF40_evolved_nf4_0000.dat
NNPDF40_evolved_nf4 PDF set, member #0, version 1
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01 5.07600900e-02  coupling_constants.py:53
                    1.20409000e-05] [5.07150400e-02 9.47604903e-01 1.68100000e-03]                                 
                    [7.43044000e-05 1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                            
                    0.23121, 'MW2': 6463.838403999999}                                                             
           INFO     {'process': 'NC', 'projectilePID': 11, 'polarization': 0.0,            coupling_constants.py:54
                    'propagatorCorrection': 0.0, 'nc_pos_charge': None}                                            
           INFO     RenScaleVar: True, FactScaleVar: True                                    scale_variations.py:62
           INFO     Atlas [0.00e+00 - 0.00e+00 - inf - inf - inf], ref=2.7224999999999997 @ 4       matchings.py:72
           INFO     PTO: 1, PTO@evolution: 1, process: NC                                             runner.py:142
           INFO     FNS: FONLL-FFNS, NfFF: 4                                                          runner.py:149
           INFO     projectile: electron, target: {Z: 1, A: 1}                                        runner.py:150
                                      ┌────────────────────────────────────┐
                                      │                                    │
                                      │ __     __       _ _                │
                                      │ \ \   / /      | (_)               │
                                      │  \ \_/ /_ _  __| |_ ___ _ __ ___   │
                                      │   \   / _` |/ _` | / __| '_ ` _ \  │
                                      │    | | (_| | (_| | \__ \ | | | | | │
                                      │    |_|\__,_|\__,_|_|___/_| |_| |_| │
                                      │                                    │
                                      └────────────────────────────────────┘
                                                       Plan
F2_charm at 2 pts
                                                    Calculation
yadism took off! please stay tuned ...
[15:32:15] INFO     computing P_qq_0 - took: 0.162115 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”>[15:32:15]</span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qq_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.162115</span> s ]8;id=624537;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=491586;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{[15:32:15]}textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qq_0 - took: textcolor{ansi-cyan-intense}{textbf{0.162115}} s ]8;id=624537;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=491586;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

[15:32:15] INFO  computing P_qq_0 - took: 0.162115 s ]8;id=624537;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=491586;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–

           INFO     computing P_qg_0 - took: 0.135772 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”> </span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qg_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.135772</span> s ]8;id=161721;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=248099;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{ }textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qg_0 - took: textcolor{ansi-cyan-intense}{textbf{0.135772}} s ]8;id=161721;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=248099;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

  INFO  computing P_qg_0 - took: 0.135772 s ]8;id=161721;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=248099;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–



took 0.72 s
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_evolved_nf4/NNPDF40_evolved_nf4_0000.dat
NNPDF40_evolved_nf4 PDF set, member #0, version 1
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01 5.07600900e-02  coupling_constants.py:53
                    1.20409000e-05] [5.07150400e-02 9.47604903e-01 1.68100000e-03]                                 
                    [7.43044000e-05 1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                            
                    0.23121, 'MW2': 6463.838403999999}                                                             
           INFO     {'process': 'NC', 'projectilePID': 11, 'polarization': 0.0,            coupling_constants.py:54
                    'propagatorCorrection': 0.0, 'nc_pos_charge': None}                                            
           INFO     RenScaleVar: True, FactScaleVar: True                                    scale_variations.py:62
           INFO     Atlas [0.00e+00 - 0.00e+00 - inf - inf - inf], ref=2.7224999999999997 @ 4       matchings.py:72
           INFO     PTO: 1, PTO@evolution: 1, process: NC                                             runner.py:142
           INFO     FNS: FONLL-FFN0, NfFF: 4                                                          runner.py:149
           INFO     projectile: electron, target: {Z: 1, A: 1}                                        runner.py:150
                                      ┌────────────────────────────────────┐
                                      │                                    │
                                      │ __     __       _ _                │
                                      │ \ \   / /      | (_)               │
                                      │  \ \_/ /_ _  __| |_ ___ _ __ ___   │
                                      │   \   / _` |/ _` | / __| '_ ` _ \  │
                                      │    | | (_| | (_| | \__ \ | | | | | │
                                      │    |_|\__,_|\__,_|_|___/_| |_| |_| │
                                      │                                    │
                                      └────────────────────────────────────┘
                                                       Plan
F2_charm at 2 pts
                                                    Calculation
yadism took off! please stay tuned ...
[15:32:16] INFO     computing P_qq_0 - took: 0.162971 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”>[15:32:16]</span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qq_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.162971</span> s ]8;id=291458;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=651516;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{[15:32:16]}textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qq_0 - took: textcolor{ansi-cyan-intense}{textbf{0.162971}} s ]8;id=291458;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=651516;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

[15:32:16] INFO  computing P_qq_0 - took: 0.162971 s ]8;id=291458;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=651516;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–

           INFO     computing P_qg_0 - took: 0.140963 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”> </span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qg_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.140963</span> s ]8;id=720497;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=82799;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{ }textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qg_0 - took: textcolor{ansi-cyan-intense}{textbf{0.140963}} s ]8;id=720497;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=82799;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

  INFO  computing P_qg_0 - took: 0.140963 s ]8;id=720497;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=82799;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–



took 0.74 s
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Interpolation: number of points = 30, polynomial degree = 4, logarithmic = interpolation.py:538
                    True                                                                                           
LHAPDF 6.5.0 loading /Users/giacomomagni/.conda/envs/yad_dev/share/LHAPDF/NNPDF40_evolved_nf5/NNPDF40_evolved_nf5_0000.dat
NNPDF40_evolved_nf5 PDF set, member #0, version 1
           INFO     {'MZ2': 8315.178393760001, 'CKM': CKM([[9.49221518e-01 5.07600900e-02  coupling_constants.py:53
                    1.20409000e-05] [5.07150400e-02 9.47604903e-01 1.68100000e-03]                                 
                    [7.43044000e-05 1.62409000e-03 9.98304719e-01]]), 'sin2theta_weak':                            
                    0.23121, 'MW2': 6463.838403999999}                                                             
           INFO     {'process': 'NC', 'projectilePID': 11, 'polarization': 0.0,            coupling_constants.py:54
                    'propagatorCorrection': 0.0, 'nc_pos_charge': None}                                            
           INFO     RenScaleVar: True, FactScaleVar: True                                    scale_variations.py:62
           INFO     Atlas [0.00e+00 - 0.00e+00 - 0.00e+00 - inf - inf], ref=2.7224999999999997 @ 4  matchings.py:72
           INFO     PTO: 1, PTO@evolution: 1, process: NC                                             runner.py:142
           INFO     FNS: FONLL-FFNS, NfFF: 5                                                          runner.py:149
           INFO     projectile: electron, target: {Z: 1, A: 1}                                        runner.py:150
                                      ┌────────────────────────────────────┐
                                      │                                    │
                                      │ __     __       _ _                │
                                      │ \ \   / /      | (_)               │
                                      │  \ \_/ /_ _  __| |_ ___ _ __ ___   │
                                      │   \   / _` |/ _` | / __| '_ ` _ \  │
                                      │    | | (_| | (_| | \__ \ | | | | | │
                                      │    |_|\__,_|\__,_|_|___/_| |_| |_| │
                                      │                                    │
                                      └────────────────────────────────────┘
                                                       Plan
F2_charm at 2 pts
                                                    Calculation
yadism took off! please stay tuned ...
[15:32:16] INFO     computing P_qq_0 - took: 0.162901 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”>[15:32:16]</span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qq_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.162901</span> s ]8;id=183465;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=829679;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{[15:32:16]}textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qq_0 - took: textcolor{ansi-cyan-intense}{textbf{0.162901}} s ]8;id=183465;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=829679;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

[15:32:16] INFO  computing P_qq_0 - took: 0.162901 s ]8;id=183465;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=829679;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–

           INFO     computing P_qg_0 - took: 0.135994 s                                      scale_variations.py:85
Starting... ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━   0% -:--:--
<span class=”ansi-cyan-fg”> </span><span class=”ansi-cyan-fg”> </span><span class=”ansi-blue-fg”>INFO </span> computing P_qg_0 - took: <span class=”ansi-cyan-intense-fg ansi-bold”>0.135994</span> s ]8;id=170545;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=485481;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… <span style=”color: rgb(58,58,58)”>━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━</span> <span class=”ansi-magenta-fg”> 0%</span> <span class=”ansi-cyan-fg”>-:–:–</span> </pre>

textcolor{ansi-cyan}{ }textcolor{ansi-cyan}{ }textcolor{ansi-blue}{INFO } computing P_qg_0 - took: textcolor{ansi-cyan-intense}{textbf{0.135994}} s ]8;id=170545;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pytextbackslash{}scale_variations.py]8;;textbackslash{}:]8;id=485481;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#85textbackslash{}85]8;;textbackslash{}

Starting{ldots} deftcRGB{textcolor[RGB]}expandaftertcRGBexpandafter{detokenize{58,58,58}}{━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━} textcolor{ansi-magenta}{ 0%} textcolor{ansi-cyan}{-:–:–} end{sphinxVerbatim}

  INFO  computing P_qg_0 - took: 0.135994 s ]8;id=170545;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.pyscale_variations.py]8;;:]8;id=485481;file:///Volumes/Git_Workspace/physicstools/NN3PDF/yadism/src/yadism/esf/scale_variations.py#8585]8;;

Starting… ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━  0% -:–:–



took 0.72 s
[15:32:17] INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
           INFO     Strong Coupling: a_s(µ_R^2=8317.440000)^(nf=5)=0.009390=0.118000/()          couplings.py:464
           INFO     Atlas [0.00e+00 - 2.28e+00 - 2.42e+01 - 2.98e+04 - inf], ref=8317.44 @ 5        matchings.py:72
[85]:
result
[85]:
{'F2_charm': [1.1565796072413521, 0.5407855765436533]}