NuRadioReco.modules.cosmicRayEnergyReconstructor module

Module to reconstruct the energy of an air shower from its radio signal

To reconstruct the air-shower energy, the radio signals should be bandpass filtered between 80 to 300 MHz with a 10th order Butterworth filter, and both the cosmic ray direction and the unfolded electric field should be present. The latter can be obtained (after the signal direction has been reconstructed) by running the NuRadioReco.modules.voltageToAnalyticEfieldConverter module.

For more details on this algorithm, see https://dx.doi.org/10.1088/1475-7516/2019/10/075

class NuRadioReco.modules.cosmicRayEnergyReconstructor.cosmicRayEnergyReconstructor

Bases: object

Reconstructs the energy of an air shower from its radio signal

Requires the following modules to be run beforehand:

• a 10th order Butterworth bandpass filter with passband 80-300 MHz

• a direction reconstruction

• the voltageToAnalyticEfieldConverter

Methods

begin([site])	Initialize the cosmicRayEnergyReconstructor (optional)
run(event, station, detector[, electric_field])	Determine the cosmic ray energy from the electric field fluence.

begin(site=None)

Initialize the cosmicRayEnergyReconstructor (optional)

Parameters:

sitestring | None (default: None)

Specifies the site of the station. The parameterization of the cosmic ray energy depends on the site of the detector.

If None, the site will be determined from the detector passed to the run function.

run(event, station, detector, electric_field=None)

Determine the cosmic ray energy from the electric field fluence.

The reconstructed cosmic ray energy will be stored in the station in the cr_energy_em parameter.

Parameters:

eventEvent

stationStation

The station containing the reconstructed electric field. If it contains multiple electric fields, only the last electric field will be used, unless another electric field is passed as a keyword-argument

detectorDetector

electric_fieldElectricField | None (default: None)

If not None, reconstruct the energy for this electric field. Otherwise, reconstruct the last electric field in the station. Useful if a station contains multiple reconstructed electric fields.

Returns:

rec_energyfloat

The reconstructed cosmic ray energy