Hello @sdw21 ,
To answer to your questions:
- I believe this is a result of the way you have now defined and use the BXDRAW entry for writing the data in your updated mgdraw.f routine. Perhaps it is more productive for you to stick to USDRAW as you did before and include the WRITE statement there, i.e. only write the data when an inelastic event occurs. This will also make your life easier for the next questions:
- As mentioned, you can take inspiration from this post to further detail the parent particle information. You can use the variables NCASE (the index of the primary source particle) to know which event history it belongs to and LTRACK indicating the particle generation number, i.e. a primary particle has LTRACK=1 and any secondary, tertiary, … particle has LTRACK>1.
- Hence, if in your mgdraw.f you select particles only secondary particles (LTRACK.EQ.1) originating from a primary neutrino parent (JTRACK.EQ.5).
Hopefully this helps, let me know in case you run into trouble.
Cheers,
Andreas