'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))
channel_pairs=channel_pairs)
voltageToEfieldConverter.run(evt,
station,
det,
use_channels=used_channels_efield)
electricFieldSignalReconstructor.run(evt, station, det)
voltageToAnalyticEfieldConverter.run(
evt,
station,
det,
use_channels=used_channels_efield,
bandpass=[80 * units.MHz, 500 * units.MHz],
use_MC_direction=False)
channelResampler.run(evt,
station,
det,
sampling_rate=1 * units.GHz)
electricFieldResampler.run(evt,
station,
det,
sampling_rate=1 * units.GHz)
eventWriter.run(evt, det)
nevents = eventWriter.end()
print("Finished processing, {} events".format(nevents))
# traditional
voltageToEfieldConverter.run(evt,
station,
det,
use_channels=used_channels_efield,
use_MC_direction=True)
electricFieldBandPassFilter.run(
evt,
station,
det,
passband=[30 * units.MHz, 80 * units.MHz],
filter_type='butter',
order=10)
electricFieldSignalReconstructor.run(evt, station, det)
# channelSignalReconstructor.run(evt, station, det)
# new analytic approach
voltageToAnalyticEfieldConverter.run(
evt,
station,
det,
use_channels=used_channels_efield,
bandpass=[30 * units.MHz, 80 * units.MHz],
use_MC_direction=True)
eventWriter.run(evt)
nevents = eventWriter.end()
print("Finished processing, {} events".format(nevents))