Header keys missing in read_corsika_headers.py?

[jebuss]

Is it possible that the following header keys are missing in read_corsika_headers.py?

• n_showers
• x_scatter
• scatter_radius
• scatter_radius_unit

irf/irf/scripts/read_corsika_headers.py

Lines 33 to 39 in bc36aed

	run_header_columns = [
	'run_number',
	'date',
	'energy_spectrum_slope',
	'energy_min',
	'energy_max',
	]

I saw in pyfact that they are needed:

https://github.com/fact-project/pyfact/blob/91273b419e296e1d94e5b85f1d257872c5d5ac1b/fact/io.py#L484

[maxnoe]

Ah yes, I did not add them in the read_corsika_headers script, only created the new eventio script

[maxnoe]

But: x_scatter is note saved into the old headers. So you need to add it manually using h5py.

with h5py.File(file, mode='r+') as f:
    f['corsika_headers'].attrs['scatter_radius'] = 270
    f['corsika_headers'].attrs['scatter_radius_unit'] = 'm'

[jebuss]

I see. THX

[maxnoe]

Also in the fact io function, it looks for either n_showers or n_events, the meaning is the same.