station,
det,
number_concidences=2,
threshold=100 * units.mV)
if station.get_trigger('default_simple_threshold').has_triggered():
channelBandPassFilter.run(
evt,
station,
det,
passband=[80 * units.MHz, 500 * units.MHz],
filter_type='butter',
order=10)
eventTypeIdentifier.run(evt, station, "forced", 'cosmic_ray')
channelStopFilter.run(evt, station, det)
channelBandPassFilter.run(
evt,
station,
det,
passband=[60 * units.MHz, 600 * units.MHz],
filter_type='rectangular')
channelSignalReconstructor.run(evt, station, det)
hardwareResponseIncorporator.run(evt, station, det)
correlationDirectionFitter.run(evt,
'If set, the run number and event id will be set to an increasing value.')
args = parser.parse_args()
# initialize modules
det = detector.GenericDetector(json_filename=args.detectordescription,
default_station=102,
default_channel=0)
readCoREASShower = NuRadioReco.modules.io.coreas.readCoREASShower.readCoREASShower(
)
readCoREASShower.begin(args.inputfilename,
det,
set_ascending_run_and_event_number=args.set_run_number)
efieldToVoltageConverter = NuRadioReco.modules.efieldToVoltageConverter.efieldToVoltageConverter(
)
efieldToVoltageConverter.begin()
eventTypeIdentifier = NuRadioReco.modules.eventTypeIdentifier.eventTypeIdentifier(
)
eventWriter = NuRadioReco.modules.io.eventWriter.eventWriter()
eventWriter.begin(args.output_filename)
for event, gen_det in readCoREASShower.run():
print('Event {} {}'.format(event.get_run_number(), event.get_id()))
for station in event.get_stations():
eventTypeIdentifier.run(event, station, 'forced', 'cosmic_ray')
efieldToVoltageConverter.run(event, station, gen_det)
eventWriter.run(event, gen_det)
nevents = eventWriter.end()
print("Finished processing, {} events".format(nevents))