Greg Archbold's Shower Library Format

This is taken from an email that Greg sent Brian Connolly.

I (AON) am using the output from gmf with SENECA to produce full shower profiles. This is accomplished using scripts/gmf2seneca, outputting a SENECA -l.dat file (reasonably documented within the actual file and on HaJo's site.

gmf2seneca also calls scripts/seneca2greg which stores the output from SENECA in a file which (should) correspond to Greg's shower library. The format of that file is as follows (taken from an email from Greg)

p q 20.0 588
-1 42.43 813.17 59517399040.0
1 0.00 1
1 7.07 1
1 14.14 1
1 21.21 1
1 28.28 1
1 35.36 1
1 42.43 2042
1 49.50 1292390
.
.
.
1 1442.50 6643330048
1 1449.57 6440869888
-2 7.07 813.17 62267899904.0
2 0.00 1
2 7.07 1442
2 14.14 52109

The first line of the code tells you, in order, the species, model, energy, and number of showers in the file.

A line starting with a negative number tells you that you are starting a new shower, with -N being the start of shower N. The rest of the line is summary info for that shower: first interaction depth (depth where there are more than one charged particle, so it must have interacted), Xmax, and Nmax. The following lines--all the way until the next -N line--give the bin-by-bin number of charge particles in the format shower number, X (gm/cm2), number of charged particles.

NOTE: I am using bins as produced by SENECA, don't know if that matters yet

NOTE: I have more info in the filename

      Greg: mc20.0pq.txt proton QGSjet
            mc20.0ps.txt proton SIBYLL
            mc20.0fq.txt iron QGSJet

(unfinished)