Fixing bug in error calculation.

rdstar/offline_analysis/fitting/dedicated_fit_approach/dedicated_fit_routine.py

@@ -594,7 +594,7 @@ class SysResultItem:
                 tuple(base_result.param_values[:2]),
             )
 
-        errors: np.ndarray = base_result.errors[:2] - this_result.errors[:2]
+        errors: np.ndarray = np.sqrt(np.square(np.abs(base_result.errors[:2]) - np.square(this_result.errors[:2])))
         relative_errors: np.ndarray = errors / base_result.param_values[:2]
 
         rd_error: float = errors[0]
@@ -1000,8 +1000,8 @@ class RDStarFitEvaluator:
             logger.info(
                 f"\tResults for sys_name = {sys_name}\n"
                 f"\t\t the results are\n"
-                f"\t\t\t sigma R(D) = {_r.rd_error} (absolute) or {_r.rd_relative_error} (relative)"
-                f"\t\t\t sigma R(D*) = {_r.rd_star_error} (absolute) or {_r.rd_star_relative_error} (relative)"
+                f"\t\t\t sigma R(D) = {_r.rd_error} (absolute) or {_r.rd_relative_error * 100} % (relative)\n"
+                f"\t\t\t sigma R(D*) = {_r.rd_star_error} (absolute) or {_r.rd_star_relative_error * 100} % (relative)"
             )
 
         # TODO: Loop also over background setups!