import NuRadioReco.detector.antennapattern import matplotlib.pyplot as plt from NuRadioReco.utilities import units import logging import numpy as np from NuRadioReco.modules.base import module logger = module.setup_logger(level=logging.DEBUG) provider = NuRadioReco.detector.antennapattern.AntennaPatternProvider() theta = 90 * units.deg ff = np.linspace(50 * units.MHz, 1 * units.GHz, 1000) bicone = provider.load_antenna_pattern("bicone_v8_inf_n1.78") bicone_n14 = provider.load_antenna_pattern("bicone_v8_inf_n1.4") bicone_air = provider.load_antenna_pattern("bicone_v8_InfAir") bicone_XFDTD = provider.load_antenna_pattern( "XFDTD_Vpol_CrossFeed_150mmHole_n1.78") VELs = bicone.get_antenna_response_vectorized(ff, 90 * units.deg, np.deg2rad(0), np.deg2rad(180), 0, np.deg2rad(90), np.deg2rad(0)) fig, (ax) = plt.subplots(1, 1, sharey=True) ax.plot(ff / units.MHz, np.abs(VELs['theta']), label='eTheta bicone n=1.78 (WIPL-D)') # ax.plot(ff / units.MHz, np.abs(VELs['phi']), label='ePhi bicone (WIPL-D)') VELs = bicone_XFDTD.get_antenna_response_vectorized(ff, 90 * units.deg,
import numpy as np import uuid import glob from NuRadioReco.eventbrowser.app import app from NuRadioReco.eventbrowser.apps import overview from NuRadioReco.eventbrowser.apps import traces from NuRadioReco.eventbrowser.apps import cosmic_rays from NuRadioReco.eventbrowser.apps import simulation import os import argparse import NuRadioReco.eventbrowser.dataprovider import logging import webbrowser from NuRadioReco.modules.base import module logger = module.setup_logger(level=logging.INFO) argparser = argparse.ArgumentParser( description= "Starts the Event Display, which then can be accessed via a webbrowser") argparser.add_argument('file_location', type=str, help="Path of folder or filename.") argparser.add_argument( '--open-window', const=True, default=False, action='store_const', help="Open the event display in a new browser tab on startup") argparser.add_argument('--port', default=8080,
import NuRadioReco.modules.channelGenericNoiseAdder import NuRadioReco.modules.trigger.simpleThreshold import NuRadioReco.modules.channelBandPassFilter import NuRadioReco.modules.eventTypeIdentifier import NuRadioReco.modules.channelStopFilter import NuRadioReco.modules.channelSignalReconstructor import NuRadioReco.modules.correlationDirectionFitter import NuRadioReco.modules.voltageToEfieldConverter import NuRadioReco.modules.electricFieldSignalReconstructor import NuRadioReco.modules.voltageToAnalyticEfieldConverter import NuRadioReco.modules.channelResampler import NuRadioReco.modules.electricFieldResampler import NuRadioReco.modules.io.eventWriter # Logging level import logging logger = module.setup_logger(name='NuRadioReco', level=logging.WARNING) plt.switch_backend('agg') """ Here, we show an example reconstruction of CoREAS data. A variety of modules are being used. Please refer to details in the modules themselves. Input parameters (all with a default provided) --------------------- Command line input: python FullReconstruction.py station_id input_file detector_file station_id: int station id to be used, default 32 input_file: str
type=str, help='outputfilename for the snr files') parser.add_argument('--outputfilenameNuRadioReco', type=str, nargs='?', default=None, help='outputfilename of NuRadioReco detector sim file') parser.add_argument('--nchannels', type=int, help='number of channels to phase', default=4) args = parser.parse_args() n_channels = args.nchannels logger = module.setup_logger(level=logging.WARNING) # initialize detector sim modules efieldToVoltageConverter = NuRadioReco.modules.efieldToVoltageConverter.efieldToVoltageConverter( ) efieldToVoltageConverter.begin(debug=False) triggerSimulator = NuRadioReco.modules.phasedarray.triggerSimulator.triggerSimulator( ) channelResampler = NuRadioReco.modules.channelResampler.channelResampler() channelBandPassFilter = NuRadioReco.modules.channelBandPassFilter.channelBandPassFilter( ) channelGenericNoiseAdder = NuRadioReco.modules.channelGenericNoiseAdder.channelGenericNoiseAdder( ) thresholdSimulator = NuRadioReco.modules.trigger.simpleThreshold.triggerSimulator( )