pyhyperiso 1.0.0__tar.gz → 1.0.2__tar.gz

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyhyperiso
-Version: 1.0.0
+Version: 1.0.2
 Summary: Hyperiso Python binding (C++ backend via pybind11)
 Author: Niels Fardeau
 Author-Email: Théo Reymermier <t.reymermier@ip2i.in2p3.fr>
@@ -157,9 +157,9 @@ cmake --build build-examples -j
 The C++ project, Python package and release tag should use the same semantic version. For example:
 
 ```text
-CMake project version: 1.0.0
-Python package version: 1.0.0
-Git tag: v1.0.0
+CMake project version: 1.0.2
+Python package version: 1.0.2
+Git tag: v1.0.2
 ```
 
 Update all version locations in the same pull request.

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/README.md

@@ -116,9 +116,9 @@ cmake --build build-examples -j
 The C++ project, Python package and release tag should use the same semantic version. For example:
 
 ```text
-CMake project version: 1.0.0
-Python package version: 1.0.0
-Git tag: v1.0.0
+CMake project version: 1.0.2
+Python package version: 1.0.2
+Git tag: v1.0.2
 ```
 
 Update all version locations in the same pull request.

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/CMakeLists.txt

@@ -1,5 +1,5 @@
 cmake_minimum_required(VERSION 3.20)
-project(Hyperiso VERSION 1.0.0 LANGUAGES CXX)
+project(Hyperiso VERSION 1.0.2 LANGUAGES CXX)
 
 set(CMAKE_CXX_STANDARD 20)
 set(CMAKE_CXX_STANDARD_REQUIRED True)

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Common/Mapper/Map.cpp

@@ -2,11 +2,22 @@
 
 #include "Map.h"
 
+#include "ObservableTypeId.h"
+
 using enum Observables;
 using enum Decays;
 using enum WCoef;
 using enum WGroup;
 
+namespace {
+constexpr long kFermionPolDaughter2 =
+    make_polarization_type_id(PolarizationKind::Fermion, 2);
+constexpr long kLongitudinalVectorPolDaughter1 =
+    make_polarization_type_id(PolarizationKind::LongitudinalVector, 1);
+constexpr long kTransverseFraction = kTransverseFractionTypeId;
+constexpr long kAlphaK = kAlphaKTypeId;
+}
+
 const std::map<Observables, std::string>& observable_mapping() {
     static const std::map<Observables, std::string> m = {
         {BR_BS_MUMU, "BR_Bs__mu_mu"},
@@ -569,21 +580,21 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {BR_B_XS_GAMMA,                    {5,     1,      2,  3,     22     }},
         {BR_B__D0_TAU_NU,                  {521,   1,      3,  421,  -15,  16}},
         {A_FB_B__D0_TAU_NU,                {521,   5,      3,  421,  -15,  16}},
-        {P_TAU_B__D0_TAU_NU,               {521,   92015,  3,  421,  -15,  16}},
+        {P_TAU_B__D0_TAU_NU,               {521,   kFermionPolDaughter2, 3,  421,  -15,  16}},
         {R_D0,                             {521,   9132,   3,  421,  -15,  16}},
         {BR_B0__D_TAU_NU,                  {511,   1,      3,  411,  -15,  16}},
         {A_FB_B0__D_TAU_NU,                {511,   5,      3,  411,  -15,  16}},
-        {P_TAU_B0__D_TAU_NU,               {511,   92015,  3,  411,  -15,  16}},
+        {P_TAU_B0__D_TAU_NU,               {511,   kFermionPolDaughter2, 3,  411,  -15,  16}},
         {R_D,                              {511,   9132,   3,  411,  -15,  16}},
         {BR_B__DSTAR0_TAU_NU,              {521,   1,      3,  423,  -15,  16}},
         {A_FB_B__DSTAR0_TAU_NU,            {521,   5,      3,  423,  -15,  16}},
-        {P_TAU_B__DSTAR0_TAU_NU,           {521,   92015,  3,  423,  -15,  16}},
-        {P_D_B__DSTAR0_TAU_NU,             {521,   921423, 3,  423,  -15,  16}},
+        {P_TAU_B__DSTAR0_TAU_NU,           {521,   kFermionPolDaughter2, 3,  423,  -15,  16}},
+        {P_D_B__DSTAR0_TAU_NU,             {521,   kLongitudinalVectorPolDaughter1, 3,  423,  -15,  16}},
         {R_DSTAR0,                         {521,   9132,   3,  423,  -15,  16}},
         {BR_B0__DSTAR_TAU_NU,              {511,   1,      3,  413,  -15,  16}},
         {A_FB_B0__DSTAR_TAU_NU,            {511,   5,      3,  413,  -15,  16}},
-        {P_TAU_B0__DSTAR_TAU_NU,           {511,   92015,  3,  413,  -15,  16}},
-        {P_D_B0__DSTAR_TAU_NU,             {511,   921423, 3,  413,  -15,  16}},
+        {P_TAU_B0__DSTAR_TAU_NU,           {511,   kFermionPolDaughter2, 3,  413,  -15,  16}},
+        {P_D_B0__DSTAR_TAU_NU,             {511,   kLongitudinalVectorPolDaughter1, 3,  413,  -15,  16}},
         {R_DSTAR,                          {511,   9132,   3,  413,  -15,  16}},
         {BR_B__Xs_e_e,                     {5,     1,      3,  3,     11, -11}},
         {BR_B__Xs_mu_mu,                   {5,     1,      3,  3,     13, -13}},
@@ -599,7 +610,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B__KSTAR_E_E,             {521,   50,     3,  323,   11, -11}},
         {A_CP_B__KSTAR_E_E,                {521,   3,      3,  323,   11, -11}},      
         {F_L_B__KSTAR_E_E,                 {521,   931,    3,  323,   11, -11}},
-        {F_T_B__KSTAR_E_E,                 {521,   932,    3,  323,   11, -11}},
+        {F_T_B__KSTAR_E_E,                 {521,   kTransverseFraction,    3,  323,   11, -11}},
         {A_T_1_B__KSTAR_E_E,               {521,   9411,   3,  323,   11, -11}},
         {A_T_2_B__KSTAR_E_E,               {521,   9412,   3,  323,   11, -11}},
         {A_T_3_B__KSTAR_E_E,               {521,   9413,   3,  323,   11, -11}},
@@ -608,7 +619,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B__KSTAR_E_E,              {521,   9416,   3,  323,   11, -11}},
         {A_T_RE_CPV_B__KSTAR_E_E,          {521,  -9416,   3,  323,   11, -11}},
         {A_IM_B__KSTAR_E_E,                {521,   942,    3,  323,   11, -11}},
-        {ALPHA_K_B__KSTAR_E_E,             {521,   932,    3,  323,   11, -11}},
+        {ALPHA_K_B__KSTAR_E_E,             {521,   kAlphaK,    3,  323,   11, -11}},
         {H_T_1_B__KSTAR_E_E,               {521,   9431,   3,  323,   11, -11}},
         {H_T_2_B__KSTAR_E_E,               {521,   9432,   3,  323,   11, -11}},
         {H_T_3_B__KSTAR_E_E,               {521,   9433,   3,  323,   11, -11}},
@@ -659,7 +670,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B__KSTAR_MU_MU,           {521,   50,     3,  323,   13, -13}},
         {A_CP_B__KSTAR_MU_MU,              {521,   3,      3,  323,   13, -13}},      
         {F_L_B__KSTAR_MU_MU,               {521,   931,    3,  323,   13, -13}},
-        {F_T_B__KSTAR_MU_MU,               {521,   932,    3,  323,   13, -13}},
+        {F_T_B__KSTAR_MU_MU,               {521,   kTransverseFraction,    3,  323,   13, -13}},
         {A_T_1_B__KSTAR_MU_MU,             {521,   9411,   3,  323,   13, -13}},
         {A_T_2_B__KSTAR_MU_MU,             {521,   9412,   3,  323,   13, -13}},
         {A_T_3_B__KSTAR_MU_MU,             {521,   9413,   3,  323,   13, -13}},
@@ -668,7 +679,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B__KSTAR_MU_MU,            {521,   9416,   3,  323,   13, -13}},
         {A_T_RE_CPV_B__KSTAR_MU_MU,        {521,  -9416,   3,  323,   13, -13}},
         {A_IM_B__KSTAR_MU_MU,              {521,   942,    3,  323,   13, -13}},
-        {ALPHA_K_B__KSTAR_MU_MU,           {521,   932,    3,  323,   13, -13}},
+        {ALPHA_K_B__KSTAR_MU_MU,           {521,   kAlphaK,    3,  323,   13, -13}},
         {H_T_1_B__KSTAR_MU_MU,             {521,   9431,   3,  323,   13, -13}},
         {H_T_2_B__KSTAR_MU_MU,             {521,   9432,   3,  323,   13, -13}},
         {H_T_3_B__KSTAR_MU_MU,             {521,   9433,   3,  323,   13, -13}},
@@ -719,7 +730,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B__KSTAR_TAU_TAU,         {521,   50,     3,  323,   15, -15}},
         {A_CP_B__KSTAR_TAU_TAU,            {521,   3,      3,  323,   15, -15}},      
         {F_L_B__KSTAR_TAU_TAU,             {521,   931,    3,  323,   15, -15}},
-        {F_T_B__KSTAR_TAU_TAU,             {521,   932,    3,  323,   15, -15}},
+        {F_T_B__KSTAR_TAU_TAU,             {521,   kTransverseFraction,    3,  323,   15, -15}},
         {A_T_1_B__KSTAR_TAU_TAU,           {521,   9411,   3,  323,   15, -15}},
         {A_T_2_B__KSTAR_TAU_TAU,           {521,   9412,   3,  323,   15, -15}},
         {A_T_3_B__KSTAR_TAU_TAU,           {521,   9413,   3,  323,   15, -15}},
@@ -728,7 +739,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B__KSTAR_TAU_TAU,          {521,   9416,   3,  323,   15, -15}},
         {A_T_RE_CPV_B__KSTAR_TAU_TAU,      {521,  -9416,   3,  323,   15, -15}},
         {A_IM_B__KSTAR_TAU_TAU,            {521,   942,    3,  323,   15, -15}},
-        {ALPHA_K_B__KSTAR_TAU_TAU,         {521,   932,    3,  323,   15, -15}},
+        {ALPHA_K_B__KSTAR_TAU_TAU,         {521,   kAlphaK,    3,  323,   15, -15}},
         {H_T_1_B__KSTAR_TAU_TAU,           {521,   9431,   3,  323,   15, -15}},
         {H_T_2_B__KSTAR_TAU_TAU,           {521,   9432,   3,  323,   15, -15}},
         {H_T_3_B__KSTAR_TAU_TAU,           {521,   9433,   3,  323,   15, -15}},
@@ -782,7 +793,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B0__KSTAR0_E_E,              {511,   50,     3,  313,   11, -11}},
         {A_CP_B0__KSTAR0_E_E,                 {511,   3,      3,  313,   11, -11}},      
         {F_L_B0__KSTAR0_E_E,                  {511,   931,    3,  313,   11, -11}},
-        {F_T_B0__KSTAR0_E_E,                  {511,   932,    3,  313,   11, -11}},
+        {F_T_B0__KSTAR0_E_E,                  {511,   kTransverseFraction,    3,  313,   11, -11}},
         {A_T_1_B0__KSTAR0_E_E,                {511,   9411,   3,  313,   11, -11}},
         {A_T_2_B0__KSTAR0_E_E,                {511,   9412,   3,  313,   11, -11}},
         {A_T_3_B0__KSTAR0_E_E,                {511,   9413,   3,  313,   11, -11}},
@@ -791,7 +802,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B0__KSTAR0_E_E,               {511,   9416,   3,  313,   11, -11}},
         {A_T_RE_CPV_B0__KSTAR0_E_E,           {511,  -9416,   3,  313,   11, -11}},
         {A_IM_B0__KSTAR0_E_E,                 {511,   942,    3,  313,   11, -11}},
-        {ALPHA_K_B0__KSTAR0_E_E,              {511,   932,    3,  313,   11, -11}},
+        {ALPHA_K_B0__KSTAR0_E_E,              {511,   kAlphaK,    3,  313,   11, -11}},
         {H_T_1_B0__KSTAR0_E_E,                {511,   9431,   3,  313,   11, -11}},
         {H_T_2_B0__KSTAR0_E_E,                {511,   9432,   3,  313,   11, -11}},
         {H_T_3_B0__KSTAR0_E_E,                {511,   9433,   3,  313,   11, -11}},
@@ -842,7 +853,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B0__KSTAR0_MU_MU,            {511,   50,     3,  313,   13, -13}},
         {A_CP_B0__KSTAR0_MU_MU,               {511,   3,      3,  313,   13, -13}},      
         {F_L_B0__KSTAR0_MU_MU,                {511,   931,    3,  313,   13, -13}},
-        {F_T_B0__KSTAR0_MU_MU,                {511,   932,    3,  313,   13, -13}},
+        {F_T_B0__KSTAR0_MU_MU,                {511,   kTransverseFraction,    3,  313,   13, -13}},
         {A_T_1_B0__KSTAR0_MU_MU,              {511,   9411,   3,  313,   13, -13}},
         {A_T_2_B0__KSTAR0_MU_MU,              {511,   9412,   3,  313,   13, -13}},
         {A_T_3_B0__KSTAR0_MU_MU,              {511,   9413,   3,  313,   13, -13}},
@@ -851,7 +862,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B0__KSTAR0_MU_MU,             {511,   9416,   3,  313,   13, -13}},
         {A_T_RE_CPV_B0__KSTAR0_MU_MU,         {511,  -9416,   3,  313,   13, -13}},
         {A_IM_B0__KSTAR0_MU_MU,               {511,   942,    3,  313,   13, -13}},
-        {ALPHA_K_B0__KSTAR0_MU_MU,            {511,   932,    3,  313,   13, -13}},
+        {ALPHA_K_B0__KSTAR0_MU_MU,            {511,   kAlphaK,    3,  313,   13, -13}},
         {H_T_1_B0__KSTAR0_MU_MU,              {511,   9431,   3,  313,   13, -13}},
         {H_T_2_B0__KSTAR0_MU_MU,              {511,   9432,   3,  313,   13, -13}},
         {H_T_3_B0__KSTAR0_MU_MU,              {511,   9433,   3,  313,   13, -13}},
@@ -902,7 +913,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {Q0_A_FB_B0__KSTAR0_TAU_TAU,          {511,   50,     3,  313,   15, -15}},
         {A_CP_B0__KSTAR0_TAU_TAU,             {511,   3,      3,  313,   15, -15}},      
         {F_L_B0__KSTAR0_TAU_TAU,              {511,   931,    3,  313,   15, -15}},
-        {F_T_B0__KSTAR0_TAU_TAU,              {511,   932,    3,  313,   15, -15}},
+        {F_T_B0__KSTAR0_TAU_TAU,              {511,   kTransverseFraction,    3,  313,   15, -15}},
         {A_T_1_B0__KSTAR0_TAU_TAU,            {511,   9411,   3,  313,   15, -15}},
         {A_T_2_B0__KSTAR0_TAU_TAU,            {511,   9412,   3,  313,   15, -15}},
         {A_T_3_B0__KSTAR0_TAU_TAU,            {511,   9413,   3,  313,   15, -15}},
@@ -911,7 +922,7 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_T_RE_B0__KSTAR0_TAU_TAU,           {511,   9416,   3,  313,   15, -15}},
         {A_T_RE_CPV_B0__KSTAR0_TAU_TAU,       {511,  -9416,   3,  313,   15, -15}},
         {A_IM_B0__KSTAR0_TAU_TAU,             {511,   942,    3,  313,   15, -15}},
-        {ALPHA_K_B0__KSTAR0_TAU_TAU,          {511,   932,    3,  313,   15, -15}},
+        {ALPHA_K_B0__KSTAR0_TAU_TAU,          {511,   kAlphaK,    3,  313,   15, -15}},
         {H_T_1_B0__KSTAR0_TAU_TAU,            {511,   9431,   3,  313,   15, -15}},
         {H_T_2_B0__KSTAR0_TAU_TAU,            {511,   9432,   3,  313,   15, -15}},
         {H_T_3_B0__KSTAR0_TAU_TAU,            {511,   9433,   3,  313,   15, -15}},
@@ -1063,21 +1074,21 @@ const std::map<Observables, LhaID>& observable_flha_mapping() {
         {A_FB_H_LAMBDA_B__LAMBDA_E_E,         {5122,  52,     3,  3122,  11, -11}},
         {A_FB_LH_LAMBDA_B__LAMBDA_E_E,        {5122,  53,     3,  3122,  11, -11}},
         {F_L_LAMBDA_B__LAMBDA_E_E,            {5122,  931,    3,  3122,  11, -11}},
-        {F_T_LAMBDA_B__LAMBDA_E_E,            {5122,  932,    3,  3122,  11, -11}},
+        {F_T_LAMBDA_B__LAMBDA_E_E,            {5122,  kTransverseFraction,    3,  3122,  11, -11}},
         {DGAMMA_DQ2_LAMBDA_B__LAMBDA_MU_MU,   {5122,  10,     3,  3122,  13, -13}}, 
         {DBR_DQ2_LAMBDA_B__LAMBDA_MU_MU,      {5122,  1,      3,  3122,  13, -13}}, 
         {A_FB_L_LAMBDA_B__LAMBDA_MU_MU,       {5122,  51,     3,  3122,  13, -13}},
         {A_FB_H_LAMBDA_B__LAMBDA_MU_MU,       {5122,  52,     3,  3122,  13, -13}},
         {A_FB_LH_LAMBDA_B__LAMBDA_MU_MU,      {5122,  53,     3,  3122,  13, -13}},
         {F_L_LAMBDA_B__LAMBDA_MU_MU,          {5122,  931,    3,  3122,  13, -13}},
-        {F_T_LAMBDA_B__LAMBDA_MU_MU,          {5122,  932,    3,  3122,  13, -13}}, 
+        {F_T_LAMBDA_B__LAMBDA_MU_MU,          {5122,  kTransverseFraction,    3,  3122,  13, -13}},
         {DGAMMA_DQ2_LAMBDA_B__LAMBDA_TAU_TAU, {5122,  10,      3,  3122,  15, -15}}, 
         {DBR_DQ2_LAMBDA_B__LAMBDA_TAU_TAU,    {5122,  1,      3,  3122,  15, -15}}, 
         {A_FB_L_LAMBDA_B__LAMBDA_TAU_TAU,     {5122,  51,     3,  3122,  15, -15}},
         {A_FB_H_LAMBDA_B__LAMBDA_TAU_TAU,     {5122,  52,     3,  3122,  15, -15}},
         {A_FB_LH_LAMBDA_B__LAMBDA_TAU_TAU,    {5122,  53,     3,  3122,  15, -15}},
         {F_L_LAMBDA_B__LAMBDA_TAU_TAU,        {5122,  931,    3,  3122,  15, -15}},
-        {F_T_LAMBDA_B__LAMBDA_TAU_TAU,        {5122,  932,    3,  3122,  15, -15}},
+        {F_T_LAMBDA_B__LAMBDA_TAU_TAU,        {5122,  kTransverseFraction,    3,  3122,  15, -15}},
         {PHI_D,                               {511,   71,     1, -511}},
         {DELTA_M_BD,                          {511,   7,      1, -511}},
         {PHI_S,                               {531,   71,     1, -531}},

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Common/Mapper/observable_ids.cpp

@@ -2,10 +2,36 @@
 
 #include "BinnedObservableId.h"
 #include "decay_ids.hpp"
+#include "ObservableTypeId.h"
 
+#include <optional>
 #include <stdexcept>
 #include <utility>
 
+namespace {
+
+std::optional<LhaID> canonicalize_legacy_polarization_id(const LhaID& ext) {
+    auto parts = ext.get_parts();
+    if (parts.size() < 2) {
+        return std::nullopt;
+    }
+
+    if (parts[1] == kLegacyFermionPolarizationTypeId) {
+        parts[1] = make_polarization_type_id(PolarizationKind::Fermion, 2);
+    } else if (parts[1] == kLegacyLongitudinalVectorPolarizationTypeId) {
+        parts[1] = make_polarization_type_id(
+            PolarizationKind::LongitudinalVector,
+            1
+        );
+    } else {
+        return std::nullopt;
+    }
+
+    return LhaID(std::move(parts));
+}
+
+} // namespace
+
 Observables ObservableMapper::enum_elt_legacy(std::string_view s) {
     const auto key = nk(s);
 
@@ -19,7 +45,7 @@ Observables ObservableMapper::enum_elt_legacy(std::string_view s) {
 }
 
 Observables ObservableMapper::enum_elt(const LhaID& ext) {
-    auto maybeId = from_external(ext);
+    auto maybeId = from_flha(ext);
     if (!maybeId) {
         throw std::out_of_range("Unknown observable (FLHA): " + ext.to_string());
     }
@@ -33,7 +59,15 @@ Observables ObservableMapper::enum_elt(const LhaID& ext) {
 }
 
 std::optional<ObservableId> ObservableMapper::from_flha(const LhaID& ext) {
-    return from_external(ext);
+    if (auto direct = from_external(ext)) {
+        return direct;
+    }
+
+    const auto canonical = canonicalize_legacy_polarization_id(ext);
+    if (!canonical) {
+        return std::nullopt;
+    }
+    return from_external(*canonical);
 }
 
 std::optional<LhaID> ObservableMapper::flha_of(const ObservableId& id) {

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Common/Mapper/observable_ids.hpp

@@ -110,6 +110,11 @@ public:
     /**
      * @brief Resolve a FLHA id to a dynamic observable id.
      *
+     * The canonical v1.0.2 polarization pattern is @c 92ij. For input
+     * compatibility, the unambiguous v1.0.0--v1.0.1 polarization identifiers
+     * @c 92015 and @c 921423 are accepted and translated to @c 9212 and
+     * @c 9221, respectively. Output always uses the canonical identifiers.
+     *
      * @param ext FLHA identifier.
      * @return Observable id if the FLHA id is known.
      */

pyhyperiso-1.0.2/core/src/Common/ObservableTypeId.h

@@ -0,0 +1,73 @@
+#ifndef OBSERVABLE_TYPE_ID_H
+#define OBSERVABLE_TYPE_ID_H
+
+#include <stdexcept>
+
+/**
+ * @file ObservableTypeId.h
+ * @brief Helpers for project-defined FLHA observable-type identifiers.
+ */
+
+/**
+ * @brief Polarization category encoded by the project-specific FLHA pattern
+ *        @c 92ij.
+ *
+ * The decimal digit @c i identifies the polarization category:
+ *
+ * - 1: fermion polarization;
+ * - 2: longitudinal vector polarization;
+ * - 3: transverse vector polarization.
+ */
+enum class PolarizationKind : long {
+    Fermion = 1,
+    LongitudinalVector = 2,
+    TransverseVector = 3,
+};
+
+/** @brief FLHA observable type for the transverse fraction @f$F_T@f$. */
+inline constexpr long kTransverseFractionTypeId = 932;
+
+/** @brief FLHA observable type for the angular coefficient @f$\alpha_K@f$. */
+inline constexpr long kAlphaKTypeId = 933;
+
+/** @brief Legacy v1.0.0--v1.0.1 type used for charged-lepton polarization. */
+inline constexpr long kLegacyFermionPolarizationTypeId = 92015;
+
+/** @brief Legacy v1.0.0--v1.0.1 type used for longitudinal vector polarization. */
+inline constexpr long kLegacyLongitudinalVectorPolarizationTypeId = 921423;
+
+/**
+ * @brief Build a project-defined FLHA polarization type id using @c 92ij.
+ *
+ * The decimal digit @c j is the one-based position of the particle of interest
+ * in the ordered list of daughter particles stored in the FLHA observable id.
+ * For example, in @f$B\to D\tau\nu@f$, the tau is daughter 2, hence its
+ * fermion-polarization type id is @c 9212.  In
+ * @f$B\to D^*\tau\nu@f$, the @f$D^*@f$ is daughter 1, hence its longitudinal
+ * vector-polarization type id is @c 9221.
+ *
+ * @param kind Polarization category encoded in digit @c i.
+ * @param daughter_position One-based daughter position encoded in digit @c j.
+ * @return The four-digit type id @c 92ij.
+ * @throws std::invalid_argument if @p kind is outside 1--3 or if
+ *         @p daughter_position is outside 1--9.
+ */
+constexpr long make_polarization_type_id(
+    PolarizationKind kind,
+    unsigned int daughter_position
+) {
+    const long kind_digit = static_cast<long>(kind);
+    if (kind_digit < 1 || kind_digit > 3) {
+        throw std::invalid_argument("polarization kind must be encoded by i=1,2,3");
+    }
+    if (daughter_position < 1 || daughter_position > 9) {
+        throw std::invalid_argument("polarized daughter position must be in [1,9]");
+    }
+    return 9200L + 10L * kind_digit + static_cast<long>(daughter_position);
+}
+
+static_assert(make_polarization_type_id(PolarizationKind::Fermion, 2) == 9212);
+static_assert(make_polarization_type_id(PolarizationKind::LongitudinalVector, 1) == 9221);
+static_assert(make_polarization_type_id(PolarizationKind::TransverseVector, 1) == 9231);
+
+#endif // OBSERVABLE_TYPE_ID_H

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Core/adapters/APIAdapter.h

@@ -49,7 +49,9 @@
  * @see IMonitor
  * @see IPathProvider
  */
-class APIAdapter : public IMonitor<ExternalFlag>, IPathProvider<APIPath>, IDataMonitor {
+class APIAdapter : public IMonitor<ExternalFlag>,
+                   public IPathProvider<APIPath>,
+                   public IDataMonitor {
 public:
     /**
      * @brief Checks the status of a specific external flag.
@@ -71,7 +73,7 @@ public:
      * @param path_name The requested path.
      * @return The filesystem path corresponding to the enum.
      */
-    fs::path get_path(APIPath path_name);
+    fs::path get_path(APIPath path_name) override;
 
     /**
      * @brief Retrieves all block names across all available model types.

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Core/domain/ParamOptimizer.cpp

@@ -20,54 +20,66 @@ void ParamOptimizer::commit(bool coalesce) {
 
     freeze_all_();
 
-    const auto plan = coalesce ? coalesce_ops_() : ops_;
-
-    std::unordered_map<std::pair<std::string,std::string>, bool, BAKeyHash> existed;
-    existed.reserve(plan.size());
-    for (const auto& v : plan) {
-        std::visit([&](auto&& op){
-            auto blk = find_block_(op.block);
-            existed[{op.block, op.id.to_string()}] = blk->contains(op.id);
-        }, v);
-    }
-
-    std::unordered_set<std::shared_ptr<Block>> blocks_needing_notify;
-
-    for (const auto& v : plan) {
-        std::visit([&](auto&& op) {
-            using T = std::decay_t<decltype(op)>;
-            auto blk = find_block_(op.block);
-            const bool had = existed.at({op.block, op.id.to_string()});
+    try {
+        const auto plan = coalesce ? coalesce_ops_() : ops_;
+
+        std::unordered_map<std::pair<std::string,std::string>, bool, BAKeyHash> existed;
+        existed.reserve(plan.size());
+        for (const auto& v : plan) {
+            std::visit([&](auto&& op){
+                auto blk = find_block_(op.block);
+                existed[{op.block, op.id.to_string()}] = blk->contains(op.id);
+            }, v);
+        }
 
-            if constexpr (std::is_same_v<T, OpSetValue>) {
-                if (had) {
-                    blk->assign(op.id, op.value); 
-                } else {
-                    blk->store(op.id, std::make_shared<Parameter>(ParamId(blk->get_name(), op.id), op.value, 0., 0.));
-                    blocks_needing_notify.insert(blk);
-                }
-            } else if constexpr (std::is_same_v<T, OpSetParam>) {
-                if (had) {
-                    blk->assign(op.id, op.param); 
-                } else {
-                    blk->store(op.id, op.param);
-                    blocks_needing_notify.insert(blk);
+        std::unordered_set<std::shared_ptr<Block>> blocks_needing_notify;
+
+        for (const auto& v : plan) {
+            std::visit([&](auto&& op) {
+                using T = std::decay_t<decltype(op)>;
+                auto blk = find_block_(op.block);
+                const bool had = existed.at({op.block, op.id.to_string()});
+
+                if constexpr (std::is_same_v<T, OpSetValue>) {
+                    if (had) {
+                        blk->assign(op.id, op.value);
+                    } else {
+                        blk->store(op.id, std::make_shared<Parameter>(ParamId(blk->get_name(), op.id), op.value, 0., 0.));
+                        blocks_needing_notify.insert(blk);
+                    }
+                } else if constexpr (std::is_same_v<T, OpSetParam>) {
+                    if (had) {
+                        blk->assign(op.id, op.param);
+                    } else {
+                        blk->store(op.id, op.param);
+                        blocks_needing_notify.insert(blk);
+                    }
+                } else if constexpr (std::is_same_v<T, OpRemove>) {
+                    if (!had) {
+                        return;
+                    }
+                    blk->remove(op.id);
                 }
-            } else if constexpr (std::is_same_v<T, OpRemove>) {
-                if (!had) {
-                    return;
-                }
-                blk->remove(op.id);
-            }
-        }, v);
-    }
+            }, v);
+        }
 
-    for (const auto& blk : blocks_needing_notify) {
-        if (blk) blk->notifyObservers();
-    }
+        for (const auto& blk : blocks_needing_notify) {
+            if (blk) blk->notifyObservers();
+        }
 
-    unfreeze_all_();
-    ops_.clear();
+        unfreeze_all_();
+        ops_.clear();
+    } catch (...) {
+        // A failed lookup or assignment must not leave the global parameter
+        // graph frozen, nor replay stale operations on the next likelihood call.
+        try {
+            unfreeze_all_();
+        } catch (...) {
+            // Preserve the original failure.
+        }
+        ops_.clear();
+        throw;
+    }
 }
 
 void ParamOptimizer::clear() { ops_.clear(); }

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Python/CMakeLists.txt

@@ -2,7 +2,7 @@ cmake_minimum_required(VERSION 3.20)
 if(DEFINED SKBUILD_PROJECT_VERSION)
     project(phyperiso VERSION "${SKBUILD_PROJECT_VERSION}" LANGUAGES CXX)
 else()
-    project(phyperiso VERSION 1.0.0 LANGUAGES CXX)
+    project(phyperiso VERSION 1.0.2 LANGUAGES CXX)
 endif()
 
 set(CMAKE_CXX_STANDARD 20)

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Python/src/common/binding_common.cpp

@@ -11,6 +11,7 @@
 #include "observable_ids.hpp"
 #include "decay_ids.hpp"
 #include "BinnedObservableId.h"
+#include "ExperimentObs.h"
 #include "SourcesView.h"
 #include <optional>
 #include <stdexcept>
@@ -572,6 +573,7 @@ void init_common(py::module &m) {
         .value("A_FB_B__K_TAU_TAU", Observables::A_FB_B__K_TAU_TAU)
         .value("F_H_B__K_TAU_TAU", Observables::F_H_B__K_TAU_TAU)
         .value("R_1_B__K_L_L", Observables::R_1_B__K_L_L)
+        .value("R_1_B0__K0_L_L", Observables::R_1_B0__K0_L_L)
         .value("DGAMMA_DQ2_LAMBDA_B__LAMBDA_E_E", Observables::DGAMMA_DQ2_LAMBDA_B__LAMBDA_E_E)
         .value("DBR_DQ2_LAMBDA_B__LAMBDA_E_E", Observables::DBR_DQ2_LAMBDA_B__LAMBDA_E_E)
         .value("A_FB_L_LAMBDA_B__LAMBDA_E_E", Observables::A_FB_L_LAMBDA_B__LAMBDA_E_E)
@@ -1274,4 +1276,20 @@ The object is non-owning and should only be used during the callback call.
         [](BinnedObservableId const& x) -> std::size_t { return std::hash<BinnedObservableId>{}(x); },
         py::is_method(m.attr("BinnedObservableId"))
     );
+
+    py::class_<ExperimentObs>(m, "ExperimentObs")
+        .def(py::init<std::string, BinnedObservableId>(),
+             py::arg("experiment"), py::arg("obs"))
+        .def_readwrite("experiment", &ExperimentObs::experiment)
+        .def_readwrite("obs", &ExperimentObs::obs)
+        .def("str", &ExperimentObs::str)
+        .def("__str__", &ExperimentObs::str)
+        .def("__repr__", [](ExperimentObs const& x) {
+            return std::string("ExperimentObs(") + x.str() + ")";
+        })
+        .def("__hash__", [](ExperimentObs const& x) {
+            return std::hash<ExperimentObs>{}(x);
+        })
+        .def(py::self == py::self)
+        .def(py::self < py::self);
 }

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Python/src/statistic/binding_statistic.cpp

@@ -747,6 +747,15 @@ void init_statistic_interface(py::module_& m) {
         .def("select_experiments_all", &StatisticInterface::select_experiments_all)
         .def("has_experiment_selection", &StatisticInterface::has_experiment_selection)
         .def("selected_experiments", &StatisticInterface::selected_experiments)
+        .def("select_experiment_observables",
+             &StatisticInterface::select_experiment_observables,
+             py::arg("observables"))
+        .def("select_experiment_observables_all",
+             &StatisticInterface::select_experiment_observables_all)
+        .def("has_experiment_observable_selection",
+             &StatisticInterface::has_experiment_observable_selection)
+        .def("selected_experiment_observables",
+             &StatisticInterface::selected_experiment_observables)
         .def("compute_uncertainties", &StatisticInterface::compute_uncertainties,
              py::call_guard<py::gil_scoped_release>())
         .def("compute_uncertainties_and_sampling", &StatisticInterface::compute_uncertainties_and_sampling,

{pyhyperiso-1.0.0 → pyhyperiso-1.0.2}/core/src/Statistic/CMakeLists.txt

@@ -56,6 +56,7 @@ add_library(StatisticLib STATIC
     adapters/StatParamOptimizerProxy.cpp
     adapters/StatDependencyPruner.cpp
     adapters/StatisticInterface.cpp
+    adapters/RuntimeNuisancePathsAdapter.cpp
     adapters/ObservableInterfaceProxy.cpp
     domain/NuisanceReader.cpp
     # domain/RNGHelper.cpp

pyhyperiso-1.0.2/core/src/Statistic/adapters/RuntimeNuisancePathsAdapter.cpp

@@ -0,0 +1,33 @@
+#include "RuntimeNuisancePathsAdapter.h"
+
+#include <stdexcept>
+#include <utility>
+
+#include "APIAdapter.h"
+
+RuntimeNuisancePathsAdapter::RuntimeNuisancePathsAdapter()
+    : RuntimeNuisancePathsAdapter(std::make_shared<APIAdapter>())
+{
+}
+
+RuntimeNuisancePathsAdapter::RuntimeNuisancePathsAdapter(
+    std::shared_ptr<IPathProvider<APIPath>> path_provider
+)
+    : path_provider_(std::move(path_provider))
+{
+    if (!path_provider_) {
+        throw std::invalid_argument(
+            "RuntimeNuisancePathsAdapter requires a non-null path provider"
+        );
+    }
+}
+
+fs::path RuntimeNuisancePathsAdapter::default_nuisances_path() const
+{
+    return path_provider_->get_path(APIPath::DEFAULT_NUISANCES);
+}
+
+fs::path RuntimeNuisancePathsAdapter::user_nuisances_path() const
+{
+    return path_provider_->get_path(APIPath::USER_NUISANCES);
+}