diff --git a/resources/images/README/title.key b/resources/images/README/title.key index eed519bf6..98739c8c5 100755 Binary files a/resources/images/README/title.key and b/resources/images/README/title.key differ diff --git a/resources/images/README/title/title.001.png b/resources/images/README/title/title.001.png index 52de44a17..dcf93d807 100644 Binary files a/resources/images/README/title/title.001.png and b/resources/images/README/title/title.001.png differ diff --git a/resources/images/organisation/logo.key b/resources/images/organisation/logo.key index 8a34ebc93..3561bb5ef 100755 Binary files a/resources/images/organisation/logo.key and b/resources/images/organisation/logo.key differ diff --git a/resources/images/organisation/logo/logo.001.png b/resources/images/organisation/logo/logo.001.png index 38b0b6540..2e2c3402b 100644 Binary files a/resources/images/organisation/logo/logo.001.png and b/resources/images/organisation/logo/logo.001.png differ diff --git a/src/StatisticalInference/JointProbability/include/BarlowBeestonBanff2022.h b/src/StatisticalInference/JointProbability/include/BarlowBeestonBanff2022.h index 00c6c963e..1d4cfea38 100644 --- a/src/StatisticalInference/JointProbability/include/BarlowBeestonBanff2022.h +++ b/src/StatisticalInference/JointProbability/include/BarlowBeestonBanff2022.h @@ -26,7 +26,7 @@ namespace JointProbability{ void printConfiguration() const; mutable int verboseLevel{0}; - bool throwIfInfLlh{false}; + bool throwIfInfLlh{true}; bool allowZeroMcWhenZeroData{true}; bool usePoissonLikelihood{false}; bool BBNoUpdateWeights{false}; // OA 2021 bug reimplementation @@ -72,6 +72,11 @@ namespace JointProbability{ double dataVal = samplePair_.data->getHistogram().getBinContentList()[bin_].sumWeights; double predVal = samplePair_.model->getHistogram().getBinContentList()[bin_].sumWeights; + if( predVal == 0 ){ + LogError << bin_ << " bin predicting 0 rate -> llh not defined / inf" << std::endl; + LogError << samplePair_.model->getSummary() << std::endl; + } + { /// the first time we reach this point, we assume the predMC is at its nominal value std::lock_guard g(_mutex_);