Module documentation
Fitting LORA data for new shower parameters. Currently implemented for different directions and core positions. Defaults are for set-up pf 20 LORA stations.
Run as (all parameters are of shape as default in database):
cr.trun(“LoraAnalysis”,particles_lora=particles_lora, azimuth=azimuth, elevation=elevation,core=core, core_type=core_type,use_title=True,save_plots=True,plot_prefix=’./’,plotlist=[])
Will return fit parameters as well as plots when they are stored by save_plots
Module author: Anna Nelles <a.nelles@astro.ru.nl>
-
class pycrtools.tasks.loraanalysis.LoraAnalysis(ws=None, parfile=None, **kwargs)
Refitting particle data using different shower parameters
Input parameters
- azimuth [default value: None ]
- Azimuth of arriving cosmic ray (LOFAR convention)
- core [default value: None ]
- New core position that should be fitted
- core_type [default value: “LORA” ]
- Type of core that is fitted, label on axis
- elevation [default value: None ]
- Elevation of arriving cosmic ray (LOFAR convention)
- lora_x [default value: [11.21, -29.79, -57.79, -3.79, -120.79, -82.79, -162.79, -134.79, 78.21, 155.21, 112.21, 133.21, 74.21, 118.21, 41.21, 80.21, -53.79, 3.21, -2.79, -48.79] ]
- Database parameter lora_posx, Set to default to 20 stations set-up
- lora_y [default value: [-94.07, -83.07, -125.07, -158.07, 3.93, -40.07, -26.07, -74.07, 75.93, 50.93, 24.93, 100.93, -121.07, -93.07, -76.07, -40.07, 129.93, 111.93, 61.93, 56.93] ]
- Database parameter lora_posy, Set to default to 20 stations set-up
- lora_z [default value: [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0] ]
- Database parameter lora_posz, Set to default to 20 stations set-up
- particles_lora [default value: None ]
- Database parameter lora_particle_density__m2
- plot_prefix [default value: ”./” ]
- Where to store the figures
- plot_type
- Plot type (e.g. png, jpeg, pdf)
- plotlist [default value: [] ]
- Plot list for website
- save_plots [default value: False ]
- Save the new particle LDF
- use_title [default value: “True” ]
- Plot with title
Output parameters
- lora_chi [default value: None ]
- Chi square of fit with new core
- lora_fit_parameters [default value: None ]
- Best fitting parameters of fit, N_charged, R_moliere
- lora_ndof [default value: None ]
- Number of degrees of freedom of fit
- lora_red_chi [default value: None ]
- Reduced chi square of fit with new core
-
run()
Run the task.
-
pycrtools.tasks.loraanalysis.get_distance(core, direction, positions)
Calculate distance in shower axis