From 444b31f97498ce3b49f029d7d8b6841f2f202888 Mon Sep 17 00:00:00 2001
From: Marshall Scott <mbscott@anl.gov>
Date: Thu, 28 Oct 2021 00:37:02 +0000
Subject: [PATCH] Resolve "Job Failed #408711: imaging_ecal_pion_rejection has
ZeroDivisionError: Weights sum to zero, can't be normalized"
---
benchmarks/imaging_ecal/scripts/draw_cluster.py | 0
benchmarks/imaging_ecal/scripts/draw_cluster_layers.py | 2 +-
2 files changed, 1 insertion(+), 1 deletion(-)
mode change 100644 => 100755 benchmarks/imaging_ecal/scripts/draw_cluster.py
diff --git a/benchmarks/imaging_ecal/scripts/draw_cluster.py b/benchmarks/imaging_ecal/scripts/draw_cluster.py
old mode 100644
new mode 100755
diff --git a/benchmarks/imaging_ecal/scripts/draw_cluster_layers.py b/benchmarks/imaging_ecal/scripts/draw_cluster_layers.py
index 89578104..42b80d43 100644
--- a/benchmarks/imaging_ecal/scripts/draw_cluster_layers.py
+++ b/benchmarks/imaging_ecal/scripts/draw_cluster_layers.py
@@ -110,7 +110,7 @@ if __name__ == '__main__':
else:
dfc = df[(df['eta'] <= df['eta'].quantile(0.95)) & (df['eta'] >= df['eta'].quantile(0.05)) &
(df['phi'] <= df['phi'].quantile(0.95)) & (df['phi'] >= df['phi'].quantile(0.05))]
- vec = np.average(dfc[['x', 'y', 'z']].values, axis=0, weights=dfc['energy'].values)
+ vec = np.average(dfc[['x', 'y', 'z']].values, axis=0, weights= dfc['energy'].values if np.sum(dfc['energy'].values) > 0 else None)
vec = vec/np.linalg.norm(vec)
# particle line from (0, 0, 0) to the inner Ecal surface
--
GitLab