Commit 05468ad7 authored by Alessio Berti's avatar Alessio Berti
Browse files

Do not use MAGIC numbers.

parent 8ac2ab96
Pipeline #116343 failed with stage
in 7 seconds
...@@ -722,7 +722,7 @@ class MAGICEventSource(EventSource): ...@@ -722,7 +722,7 @@ class MAGICEventSource(EventSource):
data.simulation.shower = SimulatedShowerContainer( data.simulation.shower = SimulatedShowerContainer(
energy = u.Quantity(event_data['true_energy'], u.GeV), energy = u.Quantity(event_data['true_energy'], u.GeV),
alt = Angle((np.pi/2 - event_data['true_zd']), u.rad), alt = Angle((np.pi/2 - event_data['true_zd']), u.rad),
az = Angle(-1 * (event_data['true_az'] - np.deg2rad(180 - 7)), u.rad), az = Angle(-1 * (event_data['true_az'] - np.deg2rad(180 - rot_corsika.value)), u.rad),
shower_primary_id = 1 - event_data['true_shower_primary_id'], shower_primary_id = 1 - event_data['true_shower_primary_id'],
h_first_int = u.Quantity(event_data['true_h_first_int'], u.cm), h_first_int = u.Quantity(event_data['true_h_first_int'], u.cm),
core_x = u.Quantity((event_data['true_core_x']*np.cos(rot_corsika) - event_data['true_core_y']*np.sin(rot_corsika)), u.cm), core_x = u.Quantity((event_data['true_core_x']*np.cos(rot_corsika) - event_data['true_core_y']*np.sin(rot_corsika)), u.cm),
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment