icecube.NuRadioMCInterface.recover module

class icecube.NuRadioMCInterface.recover.ParticleTypeRecover(ctx)

Bases: I3Module

This module sets the particle type of particles (of events) for which the type was not stored and set to “UNKNOWN” by the reader [1]. The particle type for the primary neutrino and daughter lepton is only set if it can be determine unambiguously from “NuRadioMC_GeneratorInfo” (omitting the difference between (anti)neutrino as it does not matter at radio energies ~> 10 PeV).

If the primary particle is a neutrino, this module does nothing and passes the frame along. If the primary particle has an “UNKNOWN” particle type, the partile type is determined from the information stored in “NuRadioMC_GeneratorInfo”. If that is not possible the frame is not passed to the next!

[1]

Due to a bug in v2.2 hdf5 files the first interaction was only stored when the shower its produce triggered the radio detector (See also https://github.com/nu-radio/NuRadioMC/pull/486). If a shower of a secondary interaction (high energy loss of a muon produced in the first interaction) triggered the radio detector, the event appears in the HDF5 file but without the information of the neutrino flavor and inelasticity of the first interaction. The inelasticity can not be “recovered” and the HDF5 files have to be patched to contain this information. The flavor can be “recovered” when the radio simulations have been produced with one flavor per file.

DAQ((I3Module)arg1, (I3Frame)arg2) → None :

C++ signature :

void DAQ(PythonModule<I3Module> {lvalue},boost::shared_ptr<I3Frame>)

Finish((I3Module)arg1) → None :

C++ signature :

void Finish(PythonModule<I3Module> {lvalue})