SLANT option makes Long files unreadable

[sultamehr]

Hello,

This might be a user issue, rather than an issue with the code. I am hoping I can get some feed back.

I am using CORSIKA v 77500. With the Slant option on, I get the following error

"ValueError: Some errors were detected !
Line #679 (got 9 columns instead of 10)"

Is there something required for the function "read_longitudinal_distributions(path)" for the slant option?

Thank you so much. Sorry for a random question. I am new to CORSIKA and found your package helpful to learn.

Best,

sultamehr

[maxnoe]

Dear Sultamehr,

Could you provide a test file, inputcard and compilation options that exhibit this issue?

[maxnoe]

There is a unit test, testing reading of longitudinal files with the SLANT option here:

https://github.com/cta-observatory/pycorsikaio/blob/main/tests/test_longitudinal.py

this is the test file used in that test:

https://github.com/cta-observatory/pycorsikaio/blob/main/tests/resources/corsika_77500_particle_slant.long

[sultamehr]

Hello!

Thank you for responding. Here is the data card:

RUNNR 757 run number
EVTNR 1 number of first shower event
NSHOW 1 number of showers to generate
PRMPAR 14 particle type of prim. particle
ESLOPE -2.7 slope of primary energy spectrum
ERANGE 1.E8 1.E8 energy range of primary particle
THETAP 0. 45. range of zenith angle (degree)
PHIP -180. 180. range of azimuth angle (degree)
SEED 1 0 0 seed for 1. random number sequence
SEED 2 0 0 seed for 2. random number sequence
OBSLEV 5.E5 observation level (in cm)
MAGNET 20.0 42.8 magnetic field centr. Europe
HADFLG 0 0 0 0 0 2 flags hadr.interact.&fragmentation
ECUTS .3 0.3 1e-4 1e-4 energy cuts for particles
MUADDI F additional info for muons
MUMULT F muon multiple scattering angle
ELMFLG F T em. interaction flags (NKG,EGS)
RADNKG 200.E2 outer radius for NKG lat.dens.distr.
LONGI T 1. F T longit.distr. & step size & fit & out
MAXPRT 100 max. number of printed events
DIRECT ./dataoutput/ output directory
USER you user
DEBUG F 6 F 1000000 debug flag and log.unit for out
EXIT terminates input

I had QGSII + GHEISHA with CERENKOV, CURVED, SLANT options. (attaching the output text file for reference).

output_QGSII_gheisha_curved_SLANT_test_with_cerenkov.txt

DATA000757_LONG.zip

Any help you can provide would be great! I will try to see if I can run the test from the package to compare with.