helper

helper.gps_to_datetime(gps)

conversion between GPS seconds and a python datetime object (taking into account leap seconds)

helper.datetime_to_gps(date)

helper.GPS_to_UTC(gps)

helper.UTC_to_GPS(utc)

helper.datetime_to_UTC(dt)

helper.spherical_to_cartesian(zenith, azimuth)

helper.cartesian_to_spherical(x, y, z)

helper.get_angle(v1, v2)

helper.get_normalized_angle(angle, degree=False)

helper.get_declination(magnetic_field_vector)

helper.get_magnetic_field_vector(site=None)

get the geomagnetic field vector in Gauss. x points to geographic East and y towards geographic North

helper.get_angle_to_magnetic_field_vector(zenith, azimuth)

returns the angle between shower axis and magnetic field

helper.get_magneticfield_azimuth(magnetic_field_declination)

helper.get_magneticfield_zenith(magnetic_field_inclination)

helper.get_lorentzforce_vector(zenith, azimuth, magnetic_field_vector=None)

helper.get_sine_angle_to_lorentzforce(zenith, azimuth, magnetic_field_vector=None)

helper.get_chargeexcess_vector(core, zenith, azimuth, stationPosition)

helper.get_chargeexcess_correction_factor(core, zenithSd, azimuthSd, stationPositions, a=0.14, magnetic_field_vector=None)

helper.get_polarization_vector_max(trace)

calculates polarization vector of efield trace (vector at maximum pulse position)

helper.get_interval_hilbert(trace, scale=0.5)

helper.get_FWHM_hilbert(trace)

helper.get_polarization_vector_FWHM(trace)

calculates polarization vector of efield trace, all vectors in the FWHM interval are averaged, the amplitude is set to the maximum of the hilbert envelope

helper.get_expected_efield_vector(core, zenith, azimuth, stationPositions, a=0.14, magnetic_field_vector=None)

helper.get_expected_efield_vector_vxB_vxvxB(station_positions, zenith, azimuth, a=0.14)

also returns the expected electric field vector, but station positions and the returned field vectors are in the vxB-vxvxB coordinate system

helper.get_angle_to_efieldexpectation_in_showerplane(efield, core, zenith, azimuth, stationPositions, a=0.14, magnetic_field_vector=None)

calculated the angle between a measured efield vector and the expectation from the geomagnetic and chargeexcess emission model. Thereby, the angular difference is evaluated in the showerfront, components not perpendicular to the shower axis are thus neglected.

helper.get_distance_to_showeraxis(core, zenith, azimuth, antennaPosition)

helper.get_position_at_height(pos, height, zenith, azimuth)

helper.get_2d_probability(x, y, xx, yy, xx_error, yy_error, xy_correlation, sigma=False)

helper.is_equal(a, b, rel_precision=1e-05)

helper.has_same_direction(zenith1, azimuth1, zenith2, azimuth2, distancecut=20)

helper.get_cherenkov_angle(h, model=1)

returns the cherenkov angle for the density at height above ground assuming that the particle speed is the speed of light

helper.get_cherenkov_ellipse(zenith, xmax, model=1)

returns the major and minor axis of the cherenkov cone projected on the ground plane reference: 10.1016/j.astropartphys.2014.04.004

helper.gaisser_hillas1_parametric(x, xmax, nmax=1)

return one parametric form of Gaisser-Hillers function Reference: http://en.wikipedia.org/wiki/Gaisser%E2%80%93Hillas_function and Darko Veberic (2012). “Lambert W Function for Applications in Physics”. Computer Physics Communications 183 (12): 2622–2628. arXiv:1209.0735. doi:10.1016/j.cpc.2012.07.008.

helper.gaisser_hillas(X, xmax, X0, lam, nmax=1)

returns the Gaisser-Hillers function Reference: http://en.wikipedia.org/wiki/Gaisser%E2%80%93Hillas_function

helper.is_confined(x, y, station_positions, delta_confinement=0)

returns True if core (x, y coordinate) is confined within stations given by ‘station_positions’. If the ‘delta_confinement’ parameter is given, the stations need this minimum distance to the core for the core to be confined.

helper.is_confined_weak(x, y, station_positions, delta_confinement=0)

returns True if core (x, y coordinate) is confined within stations given by ‘station_positions’. If the ‘delta_confinement’ parameter is given, the stations need this minimum distance to the core for the core to be confined. The criterion is weaker than in isConfined(). Here, at least one station has to be above, below, left or right the core position.

helper.in_hull(p, hull)

helper.is_confined2(x, y, station_positions, delta_confinement=0)

helper.get_efield_in_shower_plane(ex, ey, ez, zenith, azimuth)

helper.get_dirac_pulse(samples, binning=1.0, low_freq=30.0, up_freq=80.0)

generate dirac pulse

helper.rotate_vector_in_2d(v, angle)

helper.get_sd_core_error_ellipse(easting_error, northing_error, error_correlation, p)

returns semi major and semi minor axis of the confidence region with p-value p

helper.transform_error_ellipse_into_vxB_vxvxB(semi_major_axis, semi_minor_axis, zenith, azimuth)

accepts semi major and semi minor axis of an ellipse in standard auger coordinates. transforms the ellipse into vxB-vxvxB system and projects it onto the vxB-vxvxB_plane. returns the semi major and semi minor axis of the result.

helper.is_in_quantile(center, station_position, easting_error, northing_error, error_correlation, p)

returns true if station_position is within the p-quantile around center, false otherwise

helper.get_ellipse_tangents_through_point(point, semi_major_axis, semi_minor_axis)

determines the points where the tangents to an ellipse with given semi major and semi minor axis that go through point touch the ellipse. returns none if point is inside the ellipse

helper.covariance_to_correlation(M)

converts covariance matrix into correlation matrix