'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))