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