Intrinsic Charm

benchmarks/runners/sandbox.py

@@ -1,5 +1,6 @@
 # -*- coding: utf-8 -*-
 # pylint: skip-file
+# fmt: off
 # Compare the results with APFEL's
 import copy
 
@@ -20,29 +21,29 @@ def generate_observables():
     # defaults["interpolation_xgrid"] = np.geomspace(0.1, 1, 40).tolist()
     light_kin = []
     light_kin.extend(
-        [dict(x=x, Q2=10.0) for x in defaults["interpolation_xgrid"][3::3]]
+        [dict(x=x, Q2=90.0) for x in defaults["interpolation_xgrid"][3::3]]
         # np.linspace(1e-3, 1, 50)
     )
     # light_kin.extend([dict(x=x, Q2=90) for x in np.linspace(.8, .99, 10).tolist()])
-    light_kin.extend([dict(x=0.001, Q2=Q2) for Q2 in np.geomspace(4, 15, 10).tolist()])
+    light_kin.extend([dict(x=0.001, Q2=Q2) for Q2 in np.geomspace(4, 1e3, 10).tolist()])
     # light_kin.extend([dict(x=0.0051, Q2=Q2) for Q2 in np.geomspace(10, 1e5, 60).tolist()])
     # light_kin = [dict(x=0.001,Q2=1e4)]
     # light_kin.extend([dict(x=0.01, Q2=Q2) for Q2 in np.geomspace(500, 800, 10).tolist()])
     # light_kin.extend([dict(x=0.1, Q2=Q2) for Q2 in np.geomspace(4, 1e3, 10).tolist()])
     obs_list = [
         "F2light",
-        "F2charm",
+        #"F2charm",
         # "F2bottom",
         # "F2top",
-        "F2total",
-        "FLlight",
-        "FLcharm",
+        #"F2total",
+        # "FLlight",
+        #"FLcharm",
         # "FLbottom",
-        "FLtotal",
-        "F3light",
-        "F3charm",
+        # "FLtotal",
+        # "F3light",
+        #"F3charm",
         # "F3bottom",
-        "F3total",
+        # "F3total",
     ]
     cards = []
     card = copy.deepcopy(defaults)
@@ -80,6 +81,7 @@ def run_LO(self):
                 # "NfFF": self.db.theory_query.NfFF == 4,
                 # "TMC": self.db.theory_query.TMC == 0,
                 # "DAMP": self.db.theory_query.DAMP == 0,
+                # "IC": (self.db.theory_query.IC == 1),
             },
         )
 
@@ -89,17 +91,18 @@ def run_NLO(self):
             ["CT14llo_NF6"],
             {
                 # "XIR": self.db.theory_query.XIR == 2,
-                "NfFF": self.db.theory_query.NfFF == 4,
+                # "NfFF": self.db.theory_query.NfFF == 4,
                 # "XIF": self.db.theory_query.XIF == .7,
                 # "TMC": self.db.theory_query.TMC == 1,
-                "FNS": self.db.theory_query.FNS == "FONLL-A",
+                # "FNS": self.db.theory_query.FNS == "FONLL-A",
                 # "DAMP": self.db.theory_query.DAMP == 1,
+                # "IC": (self.db.theory_query.IC == 1),
             },
         )
 
 
 if __name__ == "__main__":
-    # generate_observables()
+    generate_observables()
     sand = Sandbox()
-    #sand.run_LO()
+    sand.run_LO()
     sand.run_NLO()

benchmarks/yadmark/benchmark/db_interface.py

@@ -69,6 +69,7 @@ def __init__(self, mode, external=None, assert_external=None):
             "FNS": self.theory_query.FNS == "FFNS",
             "DAMP": self.theory_query.DAMP == 0,
             "TMC": self.theory_query.TMC == 0,
+            "IC": self.theory_query.IC == 0,
         }
 
     def _load_input_from_queries(self, theory_query, obs_query):
@@ -219,7 +220,9 @@ def run_queries_external(self, theory_query, obs_query, pdfs):
                     from .external import (  # pylint:disable=import-error,import-outside-toplevel
                         apfel_utils,
                     )
-
+                    if theory["IC"] != 0 and theory["PTO"] > 0:
+                        print(yad_tab)
+                        raise ValueError("APFEL is currently not able to run")
                     ext_tab = external.get_external_data(
                         theory,
                         obs,

benchmarks/yadmark/data/theories.py

@@ -54,6 +54,7 @@ def get_matrix(self):
         # sandbox
         return {
             "PTO": [0, 1],
+            "IC": [0, 1],
             "XIR": [0.5, 0.7, 1.0, 2.0],
             "XIF": [0.5, 0.7, 1.0, 2.0],
             "TMC": [0, 1, 2, 3],

docs/sphinx/source/conf.py

@@ -62,6 +62,7 @@
     "sphinx.ext.napoleon",
     "sphinx.ext.graphviz",
     "sphinxcontrib.bibtex",
+    "sphinxcontrib.details.directive",
     "sphinx_rtd_theme",
 ]
 

docs/sphinx/source/dev-tools/benchmarks.rst

@@ -27,9 +27,6 @@ Different definition of |SF|
 Due to a different definition |SF| in |APFEL| (and so `yadism`) and |QCDNUM| is
 not possible to compare all the structure functions in all the schemes.
 
-In particular |APFEL|/`yadism` definition of |SF| is based on the coupling to
-the photon (or |DIS| boson in general), cf. :doc:`../theory/light-flavors` and
-:doc:`../theory/heavy-flavors`.
 
 On the other hand |QCDNUM| is using a different definition of the |SF| that is
 not matching the other one and from which it is not possible to recover the

docs/sphinx/source/index.rst

@@ -51,75 +51,44 @@ In particular:
 .. _PyPI: https://pypi.org/project/yadism/
 .. _pip: https://pip.pypa.io/en/stable/
 
-Features
-~~~~~~~~
-
-**Features** already implemented include:
-
-- |LO| structure functions
-- |NLO| structure functions
-
-  1. scale variations
-  2. target mass corrections
-  3. flavor number schemes
-
-The implemented structure functions are :math:`F_2,~F_L,~F_3` for the following
-|DIS| processes:
-
-- |EM|
-- |NC|
-- |CC|
-
 .. toctree::
-   :maxdepth: 1
-   :caption: Package Reference
+  :maxdepth: 1
+  :caption: Overview
+  :hidden:
 
-   yadism.rst
-   tutorials/index.rst
-
-In these :underlined:`documents` you can find the reference for all *yadism* API and internals.
+  overview/features.rst
+  overview/tutorials/index.rst
+  overview/indices.rst
 
 .. toctree::
-   :maxdepth: 1
-   :caption: Theory
-
-   theory/DIS-intro.rst
-   theory/structure-functions.rst
-   theory/coeff-funcs.rst 
-   theory/light-flavors.rst
-   theory/heavy-flavors.rst
-   theory/processes.rst
-   theory/scale-variations.rst
-   theory/TMC.rst
+  :maxdepth: 1
+  :caption: Theory
+  :hidden:
 
-   zzz-refs.rst
-
-And here there is some *physics documentation*.
-
-Some of these are trivial and will be replaced, they are here just to create a
-structure.
+  theory/DIS-intro.rst
+  theory/fact.rst 
+  theory/coeff-funcs.rst 
+  theory/fns.rst
+  theory/scale-variations.rst
 
+  zzz-refs.rst
 
 .. toctree::
-   :maxdepth: 1
-   :caption: Dev Tools
-
-   dev-tools/tests.rst
-   dev-tools/db-suite.rst
-   dev-tools/yadmark.rst
-   dev-tools/extras.rst
-   dev-tools/third-party.rst
-   dev-tools/code-todos.rst
-
-These tools have been developed alongside `yadism`, in order to automatize tests
-and improve the quality of code (less bugs, more readable... hopefully...).
-
----------
+  :maxdepth: 1
+  :caption: Implementation
+  :hidden:
 
+  Overview <yadism.rst>
+  API <modules/yadism.rst>
 
-Indices and tables
-==================
-
-* :ref:`genindex`
-* :ref:`modindex`
-* :ref:`search`
+.. toctree::
+  :maxdepth: 1
+  :caption: Dev Tools
+  :hidden:
+
+  dev-tools/tests.rst
+  dev-tools/db-suite.rst
+  dev-tools/yadmark.rst
+  dev-tools/extras.rst
+  dev-tools/third-party.rst
+  dev-tools/code-todos.rst

docs/sphinx/source/overview/features.rst

@@ -0,0 +1,16 @@
+Features
+=========
+
+- |LO| structure functions
+- |NLO| structure functions
+
+  1. scale variations
+  2. target mass corrections
+  3. flavor number schemes
+
+The implemented structure functions are :math:`F_2,~F_L,~F_3` for the following
+|DIS| processes:
+
+- |EM|
+- |NC|
+- |CC|

docs/sphinx/source/overview/indices.rst

@@ -0,0 +1,6 @@
+Indices and tables
+==================
+
+* :ref:`genindex`
+* :ref:`modindex`
+* :ref:`search`

docs/sphinx/source/package-structure/TMC.rst

@@ -1,3 +1,37 @@
 Target Mass Corrections
 =======================
 
+
+
+Work in Progress notes
+----------------------
+
+
+Runner: holds a list of structure functions instances, that will act as our
+"soft singletons"
+StructureFunction: holds a reference to his parent runner, in this way its also
+able to access its StructureFunctions siblings (like F2 from FL, in order to
+compute TMCs, e.g.)
+
+Caching: it is managed at the level of StructureFunction, while the values are
+kept at the level of ESF (so the StructureFunction is just routing the caller to
+the correct instance of ESF to ask for values)
+
+Note 2 (caching)
+----------------
+Since the responsibility of caching is of SF as written above we decided the
+following layout:
+
+- SF instantiate ESF or ESFTMC according to TMC flag in theory dictionary,
+  and append it to `self.__ESFs` at load time, i.e. in `self.load()` (these
+  are the observables to be computed)
+- when asked for output if noTMC a ESF is called and the instance is
+  registered
+
+    - `self.get_ouput()` is used for getting the result passing through:
+    - `self.get_ESF()` is used for getting the instance and register to the
+      cache
+
+- if TMC a ESFTMC is called, and whenever he needs an ESF instance to
+  compute a point it will ask its parent SF with `SF.get_ESF()` method, in
+  this way passing through the cache

docs/sphinx/source/refs/cf.bib

@@ -14,3 +14,17 @@ @article{vogt-sv
     archivePrefix={arXiv},
     primaryClass={hep-ph}
 }
+@article{Gnendiger:2017pys,
+    author = "Gnendiger, C. and others",
+    title = "{To ${d}$, or not to ${d}$: recent developments and comparisons of regularization schemes}",
+    eprint = "1705.01827",
+    archivePrefix = "arXiv",
+    primaryClass = "hep-ph",
+    reportNumber = "PSI-PR-17-06, ZU-TH-10-17, TUM-HEP-1081-17, IFIC-17-17, TIF-UNIMI-2017-4",
+    doi = "10.1140/epjc/s10052-017-5023-2",
+    journal = "Eur. Phys. J. C",
+    volume = "77",
+    number = "7",
+    pages = "471",
+    year = "2017"
+}

docs/sphinx/source/refs/gen.bib

@@ -22,4 +22,23 @@ @article{nnpdf-1.0
     year = "2009",
     note = "[Erratum: Nucl.Phys.B 816, 293 (2009)]"
 }
-
+@article{Zyla:2020zbs,
+    author = "Zyla, P.A. and others",
+    collaboration = "Particle Data Group",
+    title = "{Review of Particle Physics}",
+    doi = "10.1093/ptep/ptaa104",
+    journal = "PTEP",
+    volume = "2020",
+    number = "8",
+    pages = "083C01",
+    year = "2020"
+}
+@article{Callan:1969uq,
+    author = "Callan, Curtis G., Jr. and Gross, David J.",
+    title = "{High-energy electroproduction and the constitution of the electric current}",
+    doi = "10.1103/PhysRevLett.22.156",
+    journal = "Phys. Rev. Lett.",
+    volume = "22",
+    pages = "156--159",
+    year = "1969"
+}