BuildRooMorphing now checks if the full set of shape systematics is p…

…resent for each mass points, and throws an exception if not - otherwise would lead to a segfault when trying to access the missing Systematic objects

CombinePdfs/src/MorphFunctions.cc

@@ -72,6 +72,16 @@ void BuildRooMorphing(RooWorkspace& ws, CombineHarvester& cb,
   // Make a list of the names of shape systematics affecting this process
   vector<string> ss_vec =
       Set2Vec(cb_bp.cp().syst_type({"shape"}).syst_name_set());
+  // Now check if all shape systematics are present for all mass points
+  for (auto const& s : m_str_vec) {
+    if (cb_bp.cp().syst_type({"shape"}).mass({s}).syst_name_set().size() !=
+        ss_vec.size()) {
+      throw std::runtime_error(FNERROR(
+          "Some mass points do not have the full set of shape systematics, "
+          "this is currently unsupported"));
+    }
+  }
+
   unsigned ss = ss_vec.size();  // number of shape systematics
 
   // ls = "lnN systematic"