From ee101ca347cd584f9bc878d64bd8989c5ec959ad Mon Sep 17 00:00:00 2001
From: Jelena Petereit <108803123+jelenakhlr@users.noreply.github.com>
Date: Mon, 4 Mar 2024 10:08:22 +0100
Subject: [PATCH] update zenith & azimuth conventions zen_GRAND = zen_Corsika
 az_GRAND = 180 - az_Corsika

---
 sim2root/CoREASRawRoot/CoreasToRawROOT.py | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/sim2root/CoREASRawRoot/CoreasToRawROOT.py b/sim2root/CoREASRawRoot/CoreasToRawROOT.py
index f8ad8c30..add180c6 100644
--- a/sim2root/CoREASRawRoot/CoreasToRawROOT.py
+++ b/sim2root/CoREASRawRoot/CoreasToRawROOT.py
@@ -145,9 +145,10 @@ def CoreasToRawRoot(file, simID=None):
     GeomagneticAngle = read_params(reas_input, "GeomagneticAngle") # in degrees
 
   else:
-    #theta_GRAND = 180-theta_Corsika
-    zenith = 180 - read_params(inp_input, "THETAP")
-    azimuth = read_params(inp_input, "PHIP")
+    #theta_GRAND = theta_Corsika
+    zenith = read_params(inp_input, "THETAP")
+    #azimuth_GRAND = 180 - azimuth_Corsika
+    azimuth = 180 - read_params(inp_input, "PHIP")
 
     Energy = read_params(inp_input, "ERANGE") # in GeV
     Primary = read_params(inp_input, "PRMPAR") # as defined in CORSIKA