Comparison with IceProd Level3
-------------------------------

Mass-produced 5-component Corsika with primary energy > 30 TeV/nucleon

.. plot::
	:scale: 50
	:caption: Level 2 reconstructions
	
	from plot_observable import *
	plot_observable('spezenith', 'corsika_iceprod.hdf5', True)
	plot_observable('cscdllh_rlogl', 'corsika_iceprod.hdf5')
	plot_observable('spe_rlogl', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Level 3 cut parameters
	
	from plot_observable import *
	plot_observable('llhratio', 'corsika_iceprod.hdf5')
	plot_observable('acer_energy', 'corsika_iceprod.hdf5')
	plot_observable('fill_ratio', 'corsika_iceprod.hdf5')

.. plot::
	:scale: 50
	:caption: Distance between core-split CscdLlh vertices
	
	from plot_observable import *
	plot_observable('coresplit_distance', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Core-split LineFit speeds
	
	from plot_observable import *
	plot_observable('coresplit_linefit_speed_0', 'corsika_iceprod.hdf5')
	plot_observable('coresplit_linefit_speed_1', 'corsika_iceprod.hdf5')
	plot_observable('coresplit_linefit_speed_diff', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Core-split SPE zenith
	
	from plot_observable import *
	plot_observable('coresplit_spezenith0', 'corsika_iceprod.hdf5')
	plot_observable('coresplit_spezenith1', 'corsika_iceprod.hdf5')
	plot_observable('coresplit_spediff', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Distance between time-split CscdLlh vertices
	
	from plot_observable import *
	plot_observable('timesplit_distance', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Time-split LineFit speeds
	
	from plot_observable import *
	plot_observable('timesplit_linefit_speed_0', 'corsika_iceprod.hdf5')
	plot_observable('timesplit_linefit_speed_1', 'corsika_iceprod.hdf5')
	plot_observable('timesplit_linefit_speed_diff', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Time-split SPE zenith
	
	from plot_observable import *
	plot_observable('timesplit_spezenith0', 'corsika_iceprod.hdf5')
	plot_observable('timesplit_spezenith1', 'corsika_iceprod.hdf5')
	plot_observable('timesplit_spediff', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Reconstructions with TWSRT-cleaned pulses
	
	from plot_observable import *
	plot_observable('dipole_speed', 'corsika_iceprod.hdf5')
	plot_observable('dipole_speedz', 'corsika_iceprod.hdf5')
	plot_observable('srtlinefit_speed', 'corsika_iceprod.hdf5')
	
.. plot::
	:scale: 50
	:caption: Credo
	
	from plot_observable import *
	plot_observable('credo_zenith', 'corsika_iceprod.hdf5')
	plot_observable('credo_azimuth', 'corsika_iceprod.hdf5')
	plot_observable('credo_z', 'corsika_iceprod.hdf5')
	plot_observable('credo_energy', 'corsika_iceprod.hdf5')
	plot_observable('credo_logl', 'corsika_iceprod.hdf5')
	plot_observable('credo_rlogl', 'corsika_iceprod.hdf5')	
	
.. plot::
	:scale: 50
	:caption: Monopod
	
	from plot_observable import *
	plot_observable('monopod_zenith', 'corsika_iceprod.hdf5')
	plot_observable('monopod_azimuth', 'corsika_iceprod.hdf5')
	plot_observable('monopod_z', 'corsika_iceprod.hdf5')
	plot_observable('monopod_energy', 'corsika_iceprod.hdf5')
	plot_observable('monopod_logl', 'corsika_iceprod.hdf5')
	plot_observable('monopod_rlogl', 'corsika_iceprod.hdf5')
	plot_observable('monopod_qratio', 'corsika_iceprod.hdf5')