def _detector_simulation_filter_amp(self, evt, station, det):
hardwareResponseIncorporator.run(evt, station, det, sim_to_data=True)
channelResampler.begin()
eventWriter = NuRadioReco.modules.io.eventWriter.eventWriter()
output_filename = "MC_example_station_{}.nur".format(station_id)
eventWriter.begin(output_filename)
# Loop over all events in file as initialized in readCoRREAS and perform analysis
for iE, evt in enumerate(readCoREAS.run(detector=det)):
logger.info("processing event {:d} with id {:d}".format(iE, evt.get_id()))
station = evt.get_station(station_id)
station.set_station_time(astropy.time.Time('2019-01-01T00:00:00'))
det.update(station.get_station_time())
if simulationSelector.run(evt, station.get_sim_station(), det):
efieldToVoltageConverter.run(evt, station, det)
hardwareResponseIncorporator.run(evt, station, det, sim_to_data=True)
channelGenericNoiseAdder.run(evt,
station,
det,
type="rayleigh",
amplitude=20 * units.mV)
triggerSimulator.run(evt,
station,
det,
number_concidences=2,
threshold=100 * units.mV)
if station.get_trigger('default_simple_threshold').has_triggered():
print("{} depth = {:.1f}m -> {:.2f} {:.2f} (solution type {})".format(
t, d, zen / units.deg, az / units.deg, r.get_solution_type(0)))
channelResampler.run(evt, station, det, 50 * units.GHz)
channelBandPassFilter.run(evt,
station,
det,
passband=[120 * units.MHz, 300 * units.MHz],
filter_type='butterabs',
order=10)
channelBandPassFilter.run(evt,
station,
det,
passband=[10 * units.MHz, 1000 * units.MHz],
filter_type='rectangular')
hardwareResponseIncorporator.run(evt, station, det)
channelSignalReconstructor.run(evt, station, det)
# channelTimeWindow.run(evt, station, det, window_function='hanning', around_pulse=True, window_width=20*units.ns,
# window_rise_time=20*units.ns)
correlationDirectionFitter.run(
evt,
station,
det,
n_index=1.353,
ZenLim=[90 * units.deg, 180 * units.deg],
AziLim=[300 * units.deg, 330 * units.deg],
channel_pairs=((0, 2), (1, 3)))
print(
"reco correlation LPDAs = {:.1f} ({:.1f}) {:.1f} ({:.1f})".format(