Data implementation for Lambda production

.gitignore

@@ -55,3 +55,6 @@ inc/MontBlanc/version.h
 # Intermediary results
 testruns/
 fits/
+
+# macOS auxiliary folder
+.DS_Store

CMakeLists.txt

@@ -1,31 +1,58 @@
-cmake_minimum_required (VERSION 3.5)
-# Set a default build type if none was specified
-if(NOT CMAKE_BUILD_TYPE AND NOT CMAKE_CONFIGURATION_TYPES)
-  message(STATUS "Setting build type to 'Release' as none was specified.")
-  set(CMAKE_BUILD_TYPE Release CACHE STRING "Choose the type of build." FORCE)
-  # Set the possible values of build type for cmake-gui
-  set_property(CACHE CMAKE_BUILD_TYPE PROPERTY STRINGS "Debug" "Release"
-    "MinSizeRel" "RelWithDebInfo")
-endif()
+# ============================
+# BASIC DEFINITIONS
+# ============================
 
-if (NOT DEFINED CMAKE_MACOSX_RPATH)
-   set(CMAKE_MACOSX_RPATH 0)
-endif()
+# Define code version
+set(MontBlanc_VERSION 4.8.0)
 
-# the project name
-project(MontBlanc)
+# Define minimum version of cmake
+cmake_minimum_required(VERSION 3.10)
 
-# activating some global properties for the project
-set(MontBlanc_VERSION 1.0.0)
+# Define project name and its language
+project(MontBlanc CXX)
+
+# Define c++ standard and issue all the warning demanded by this standard
+set(CMAKE_BUILD_TYPE Release)
 set(CMAKE_CXX_STANDARD 17)
 set(CMAKE_CXX_STANDARD_REQUIRED ON)
-set(CMAKE_CXX_EXTENSIONS OFF)
+set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wall -Wextra -fPIC -O3")
 
-# export version to file
-configure_file(
-  "${PROJECT_SOURCE_DIR}/inc/MontBlanc/version.h.in"
-  "${PROJECT_SOURCE_DIR}/inc/MontBlanc/version.h"
-)
+if (NOT DEFINED CMAKE_MACOSX_RPATH)
+   set(CMAKE_MACOSX_RPATH 0)
+endif ()
+
+if (CMAKE_VERSION VERSION_GREATER_EQUAL "3.27.0")
+   cmake_policy(SET CMP0148 OLD)
+endif ()
+
+# ============================
+# SIDIS TABLES
+# ============================s
+# Set a default path (relative to the source directory)
+set(DEFAULT_SOURCE_DIR "${CMAKE_SOURCE_DIR}")
+
+# Allow the user to override the default using -DCMAKE_SOURCE_DIR=<custom_path>
+set(SOURCE_DIR ${DEFAULT_SOURCE_DIR} CACHE PATH "Path to the source directory where tables are stored")
+
+# Pass the path as a preprocessor definition
+add_compile_definitions(SOURCE_DIR="${SOURCE_DIR}")
+
+# ============================
+# MAIN
+# ============================
+
+# Define include directory
+include_directories(${PROJECT_SOURCE_DIR}/inc)
+
+# Generate list of source files
+file(GLOB_RECURSE  source_files src/*)
+
+# Define target library
+set(LIBRARY_OUTPUT_PATH ${PROJECT_BINARY_DIR}/lib)
+add_library(MontBlanc SHARED ${source_files})
+
+# Define libraries to be linked
+target_link_libraries(MontBlanc)
 
 find_package(PkgConfig REQUIRED)
 pkg_search_module(YAML REQUIRED yaml-cpp)
@@ -98,20 +125,32 @@ if (NNAD_CONFIG)
    set(NNAD_CFLAGS ${NNAD_CFLAGS} CACHE STRING INTERNAL)
 endif(NNAD_CONFIG)
 
-set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS}") 
+target_link_libraries(MontBlanc ${YAML_LDFLAGS} ${EIGEN3_LDFLAGS} ${CERES_LIBRARIES} ${GLOG_LDFLAGS} ${GFLAGS_LDFLAGS} ${APFELXX_LIBRARIES} ${NANGAPARBAT_LIBRARIES} ${LHAPDF_LIBRARIES} ${GSL_LIBRARIES})
 set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ${NNAD_CFLAGS} ${EIGEN3_CFLAGS} ${Ceres_CFLAGS} ${GLOG_CFLAGS} ${GFLAGS_CFLAGS} ${LHAPDF_CXX_FLAGS} ${APFELXX_CXX_FLAGS} ${NANGAPARBAT_CXX_FLAGS} ${GSL_CXX_FLAGS} ${YAML_CFLAGS} -Wall -fPIC -Wunused")
 
-# build target library
-add_subdirectory(src)
-add_subdirectory(run)
+# Enable testing and build test codes
 add_subdirectory(tests)
+add_subdirectory(run)
+
+# ============================
+# CONFIGURATION
+# ============================
+
+# Configure version to file
+configure_file("${PROJECT_SOURCE_DIR}/inc/MontBlanc/version.h.in" "${PROJECT_SOURCE_DIR}/inc/MontBlanc/version.h")
+
+# Configure configuration script
+set(prefix ${CMAKE_INSTALL_PREFIX})
+set(includedir "${prefix}/include")
+set(libdir "${prefix}/lib")
+set(cflags "-I${prefix}/include")
+set(ldflags "-L${prefix}/lib -lMontBlanc")
+configure_file("${PROJECT_SOURCE_DIR}/bin/MontBlanc-config.in" "${PROJECT_BINARY_DIR}/bin/MontBlanc-config")
 
-#option(BUILD_TESTS "Build tests code" ON)
-#if(BUILD_TESTS)
-#  enable_testing()
-#  add_subdirectory(tests)
-#endif(BUILD_TESTS)
+# ============================
+# INSTALLATION
+# ============================
 
-install(FILES ${PROJECT_SOURCE_DIR}/bin/MontBlanc-config
-	   DESTINATION bin
-	   PERMISSIONS OWNER_READ OWNER_WRITE OWNER_EXECUTE GROUP_READ GROUP_EXECUTE WORLD_READ WORLD_EXECUTE)
+# Installation of code
+install(DIRECTORY inc/MontBlanc DESTINATION include)
+install(FILES ${PROJECT_BINARY_DIR}/bin/MontBlanc-config DESTINATION bin PERMISSIONS OWNER_READ OWNER_WRITE OWNER_EXECUTE GROUP_READ GROUP_EXECUTE WORLD_READ WORLD_EXECUTE)

config/MAPLAMBDA/250107-ac-lambda.yaml

@@ -0,0 +1,111 @@
+#################################
+# Predictions
+#################################
+Predictions:
+  # Perturbative order, 0: LO, 1: NLO
+  perturbative order: 1
+
+  # Initial scale in GeV to be used for the DGLAP evolution of the FFs.
+  mu0: 5
+
+  # Quark thresholds
+  thresholds: [0, 0, 0, 1.51, 4.92]
+
+  # Strong coupling
+  alphas:
+    aref: 0.118
+    Qref: 91.1876
+
+  # Fine-structure constant
+  alphaem:
+    aref: 0.00776578395589
+    Qref: 91.1876
+
+  # APFEL++ grid
+  xgrid:
+    - [100, 1e-2, 3]
+    - [60,  2e-1, 3]
+    - [40,  8e-1, 3]
+
+  # PDF set
+  pdfset:
+    name: NNPDF31_nlo_pch_as_0118
+    member: -1 #N>=0 for a specific member (0 for central) ; N<0 for all members to be used randomly(flat dist)
+
+#################################
+# Optimiser
+#################################
+# Parameters of the optimiser managed by ceres-solver
+Optimizer:
+  max_num_iterations: 3000
+  chi2_tolerance: 3
+#################################
+# NNAD
+#################################
+NNAD:
+  # Initialisation seed
+  seed: 4
+  # Architecture
+  architecture: [1, 20, 7]
+  # The output function can be either the activation function of the
+  # NN (0), or linear (1), or quadratic (2). If this entry is absent a
+  # linear function is assumed.
+  output function: 2
+  # The flavour map gives the the specific combinations of
+  # distributions in the physical-basis (d, u, s, ..) to be fitted to
+  # the data. The number of combinations has to match the number of
+  # nodes of the output layer of the NN given in the
+  # architecture. When defining the flavour map, one should keep in
+  # mind that the code computes predictions using the QCD-evolution
+  # basis (Sigma, V, T3, ...). Therefore, in general,
+  # QCD-evolution-like combinations should be preferred. Moreover, the
+  # distributions are assumed to be for the positive charge of the a
+  # given hadronic species. The negative distributions are derived
+  # using charge conjugation (q->qbar).
+
+  #             tb bb cb sb ub db g  d  u  s  c  b  t  # Combinations to be fitted (i.e. output of the NN):
+  flavour map: [0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, # - sb
+                0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, # - g
+                0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, # - s = ubar
+                0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, # - u
+                0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, # - d+
+                0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, # - c+
+                0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0] # - b+
+               #tb bb cb sb ub db g  d  u  s  c  b  t
+  combine: true
+#################################
+# Data
+#################################
+Data:
+  # Seed used for the replica generation and the splitting between
+  # training and validation (do not use a too large number here in
+  # order to avoid correlation in the replica generation).
+  seed: 0
+  # Hadronic species being fitted. This is used to make sure that all
+  # data sets included in the fit actually refer to the correct
+  # species.
+  hadron: LA
+  # Datasets to be included in the fit along with specific cuts and
+  # traning fraction. Each single dataset can implement an arbitrary
+  # number of cuts determined by the name of the appropriate function
+  # (e.g. zcut) and the allowed range.
+  sets:
+    - {name: "ALEPH $\\Lambda \\bar{\\Lambda}$", file: "ALEPH_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "ARGUS $\\Lambda \\bar{\\Lambda}$", file: "ARGUS_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "BELLE $\\Lambda \\bar{\\Lambda}$", file: "BELLE_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "CELLO $\\Lambda \\bar{\\Lambda}$", file: "CELLO_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 91.2 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_91.2_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 183 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_183_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 189 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_189_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "HRS $\Lambda \\bar{\\Lambda}$", file: "HRS_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "MARK II $\\Lambda \\bar{\\Lambda}$", file: "MARKII_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "OPAL $\\Lambda \\bar{\\Lambda}$", file: "OPAL_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD B $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_B.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD C $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_C.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD UDS $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_UDS.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 14 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_14_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 22 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_22_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 33 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_33_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 34.8 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_34.8_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 42.1 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_42.1_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}

config/MAPLAMBDA/250120-ac-01-lambda.yaml

@@ -0,0 +1,112 @@
+#################################
+# Predictions
+#################################
+Predictions:
+  # Perturbative order, 0: LO, 1: NLO, 2: NNLO
+  perturbative order: 1
+
+  # Initial scale in GeV to be used for the DGLAP evolution of the FFs.
+  mu0: 5
+
+  # Quark thresholds
+  thresholds: [0, 0, 0, 1.51, 4.92]
+
+  # Strong coupling
+  alphas:
+    aref: 0.118
+    Qref: 91.1876
+
+  # Fine-structure constant
+  alphaem:
+    aref: 0.00776578395589
+    Qref: 91.1876
+
+  # APFEL++ grid
+  xgrid:
+    - [100, 1e-2, 3]
+    - [60,  2e-1, 3]
+    - [40,  8e-1, 3]
+
+  # PDF set
+  pdfset:
+    name: NNPDF31_nlo_pch_as_0118
+    member: -1 #N>=0 for a specific member (0 for central) ; N<0 for all members to be used randomly(flat dist)
+
+#################################
+# Optimiser
+#################################
+# Parameters of the optimiser managed by ceres-solver
+Optimizer:
+  max_num_iterations: 3000
+  chi2_tolerance: 3
+#################################
+# NNAD
+#################################
+NNAD:
+  # Initialisation seed
+  seed: 4
+  # Architecture
+  architecture: [1, 20, 6]
+  # The output function can be either the activation function of the
+  # NN (0), or linear (1), or quadratic (2). If this entry is absent a
+  # linear function is assumed.
+  output function: 2
+  # The flavour map gives the the specific combinations of
+  # distributions in the physical-basis (d, u, s, ..) to be fitted to
+  # the data. The number of combinations has to match the number of
+  # nodes of the output layer of the NN given in the
+  # architecture. When defining the flavour map, one should keep in
+  # mind that the code computes predictions using the QCD-evolution
+  # basis (Sigma, V, T3, ...). Therefore, in general,
+  # QCD-evolution-like combinations should be preferred. Moreover, the
+  # distributions are assumed to be for the positive charge of the a
+  # given hadronic species. The negative distributions are derived
+  # using charge conjugation (q->qbar).
+
+  #             tb bb cb sb ub db g  d  u  s  c  b  t  # Combinations to be fitted (i.e. output of the NN):
+  flavour map: [
+                0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, # - u+ 
+                0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, # - d+
+                0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, # - s+
+                0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, # - c+
+                0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, # - b+
+                0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, # - g
+               #tb bb cb sb ub db g  d  u  s  c  b  t
+               ]
+  combine: true
+#################################
+# Data
+#################################
+Data:
+  # Seed used for the replica generation and the splitting between
+  # training and validation (do not use a too large number here in
+  # order to avoid correlation in the replica generation).
+  seed: 0
+  # Hadronic species being fitted. This is used to make sure that all
+  # data sets included in the fit actually refer to the correct
+  # species.
+  hadron: LA
+  # Datasets to be included in the fit along with specific cuts and
+  # traning fraction. Each single dataset can implement an arbitrary
+  # number of cuts determined by the name of the appropriate function
+  # (e.g. zcut) and the allowed range.
+  sets:
+    - {name: "ALEPH $\\Lambda \\bar{\\Lambda}$", file: "ALEPH_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "ARGUS $\\Lambda \\bar{\\Lambda}$", file: "ARGUS_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.85}], training fraction: 0.5}
+    - {name: "BELLE $\\Lambda \\bar{\\Lambda}$", file: "BELLE_LA_LABAR_wrc.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "CELLO $\\Lambda \\bar{\\Lambda}$", file: "CELLO_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 91.2 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_91.2_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 183 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_183_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "Delphi 189 GeV $\\Lambda \\bar{\\Lambda}$", file: "DELPHI_189_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "HRS $\\Lambda \\bar{\\Lambda}$", file: "HRS_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "MARK II $\\Lambda \\bar{\\Lambda}$", file: "MARKII_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "OPAL $\\Lambda \\bar{\\Lambda}$", file: "OPAL_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD B $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_B.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD C $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_C.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD UDS $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR_UDS.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "SLD $\\Lambda \\bar{\\Lambda}$", file: "SLD_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.02, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 14 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_14_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 22 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_22_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 33 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_33_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 34.8 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_34.8_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}
+    - {name: "TASSO 42.1 GeV $\\Lambda \\bar{\\Lambda}$", file: "TASSO_42.1_LA_LABAR.yaml", cuts: [{name: zcut, min: 0.075, max: 0.9}], training fraction: 0.5}