### TABLE OF CONTENTS

## /MultiplicityAnalysis [ Modules ]

NAME

module **MultiplicityAnalysis**

PURPOSE

This module contains collision independent analysis routines concerning multiplicity of particles

INPUTS

(none)

## MultiplicityAnalysis/Multiplicity_Reset [ Subroutines ]

[ Top ] [ MultiplicityAnalysis ] [ Subroutines ]

NAME

subroutine **Multiplicity_Reset**

PURPOSE

## MultiplicityAnalysis/Multiplicity_AddEvent [ Subroutines ]

[ Top ] [ MultiplicityAnalysis ] [ Subroutines ]

NAME

subroutine **Multiplicity_AddEvent**(E)

PURPOSE

## MultiplicityAnalysis/Multiplicity_Write [ Subroutines ]

[ Top ] [ MultiplicityAnalysis ] [ Subroutines ]

NAME

subroutine **Multiplicity_Write**(Prefix,iFile)

PURPOSE

## MultiplicityAnalysis/Multiplicity_Ave.dat [ Output files ]

[ Top ] [ MultiplicityAnalysis ] [ Output files ]

NAME

file **Multiplicity_Ave.dat**

PURPOSE

The file shows the average multiplicity and its standard deviation for the various outging hadrons

Columns:

- #1: running variable (see Multiplicity_Distr.dat for explanation)
- #2: average multiplicity of "ALL" particles
- #3: standard deviation (of average multiplicity) of "ALL" particles
- #4: average multiplicity of "ALL CHARGED" particles
- #5: standard deviation (of average multiplicity) of "ALL CHARGED" particles
- #6: average multiplicity of "charge -1" particles
- #7: standard deviation (of average multiplicity) of "charge -1" particles
- #8: average multiplicity of "charge 0" particles
- #9: standard deviation (of average multiplicity) of "charge 0" particles
- #10: average multiplicity of "charge +1" particles
- #11: standard deviation (of average multiplicity) of "charge +1" particles
- #12: average multiplicity of "ALL pions"
- #13: standard deviation (of average multiplicity) of "ALL pions"
- #14: average multiplicity of "pion -1"
- #15: standard deviation (of average multiplicity) of "pion -1"
- #16: average multiplicity of "pion 0"
- #17: standard deviation (of average multiplicity) of "pion 0"
- #18: average multiplicity of "pion +1"
- #19: standard deviation (of average multiplicity) of "pion +1"
- #20: average multiplicity of "ALL nucleons"
- #21: standard deviation (of average multiplicity) of "ALL nucleons"
- #22: average multiplicity of "neutrons"
- #23: standard deviation (of average multiplicity) of "neutrons"
- #24: average multiplicity of "protons"
- #25: standard deviation (of average multiplicity) of "protons"
- #26: average multiplicity of "ALL kaons"
- #27: standard deviation (of average multiplicity) of "ALL kaons"
- #28: average multiplicity of "Kbar -1"
- #29: standard deviation (of average multiplicity) of "Kbar -1"
- #30: average multiplicity of "Kbar 0"
- #31: standard deviation (of average multiplicity) of "Kbar 0"
- #32: average multiplicity of "K 0"
- #33: standard deviation (of average multiplicity) of "K 0"
- #34: average multiplicity of "K +1"
- #35: standard deviation (of average multiplicity) of "K +1"
- #36: average multiplicity of "Lambda"
- #37: standard deviation (of average multiplicity) of "Lambda"

## MultiplicityAnalysis/Multiplicity_Distr.dat [ Output files ]

[ Top ] [ MultiplicityAnalysis ] [ Output files ]

NAME

file **Multiplicity_Distr.dat**

PURPOSE

The file shows the multiplicity distributions of various hadrons in the final state

The commented line with one number only (looks like # 1.000) is the "running variable" from your jobcard, For example, for neutrino mode runs versus neutrino energy (mode 0 or mode 6) this is neutrino energy for neutrino mode runs versus Q2 (mode 3) this is Q2 For neutrino model calculations with a given flux this 'running variable' is a dummy, its the last energy in the flux sampling For each value of "running variable" you have several (="num_runs_sameEnergy") blocks of output separated by two empty lines For flux-averaged results these blocks should be (within statistics) the same. This can be used to check if statistics are sufficient

Columns:

- #1: number of outgoing particles (integer)
- #2: percentage of events with #1 particles of type "ALL"
- #3: percentage of events with #1 particles of type "ALL CHARGED"
- #4: percentage of events with #1 particles of type "charge -1"
- #5: percentage of events with #1 particles of type "charge 0"
- #6: percentage of events with #1 particles of type "charge +1"
- #7: percentage of events with #1 particles of type "ALL pions"
- #8: percentage of events with #1 particles of type "pion -1"
- #9: percentage of events with #1 particles of type "pion 0"
- #10: percentage of events with #1 particles of type "pion +1"
- #11: percentage of events with #1 particles of type "ALL baryons"
- #12: percentage of events with #1 particles of type "neutrons"
- #13: percentage of events with #1 particles of type "protons"
- #14: percentage of events with #1 particles of type "ALL kaons" (incl antikaons)
- #15: percentage of events with #1 particles of type "Kbar -1"
- #16: percentage of events with #1 particles of type "Kbar 0"
- #17: percentage of events with #1 particles of type "K0"
- #18: percentage of events with #1 particles of type "K+"
- #19: percentage of events with #1 particles of type "Lambda"
- In the following columns all values in one column are idential, since col. 1 is irrelevant for multiplicities
- #20: average multiplicity of "ALL"
- #21: average multiplicity of "ALL CHARGED"
- #22: average multiplicity of "charge -1"
- #23: average multiplicity of "charge 0"
- #24: average multiplicity of "charge +1"
- #25: average multiplicity of "ALL pions"
- #26: average multiplicity of "pion -1"
- #27: average multiplicity of "pion 0"
- #28: average multiplicity of "pion +1"
- #29: average multiplicity of "ALL baryons"
- #30: average multiplicity of "neutrons"
- #31: average multiplicity of "protons"
- #32: average multiplicity of "ALL kaons"
- #33: average multiplicity of "Kbar -1"
- #34: average multiplicity of "Kbar 0"
- #35: average multiplicity of "K0"
- #36: average multiplicity of "K+"
- #37: average multiplicity of "Lambda"

NOTES

The numbers in one row should NOT add-up in any way

Example: consider an event with 1 neutron, 1 pi+ and 1pi0 in the final state; it will contribute to "3 ALL", "1 ALL CHARGED", "0 charge -1", "2 charge 0", "1 charge +1",

"2 ALL pions", "0 pion -1", "1 pion 0", "1 pion+", "1 ALL baryons", "1 neutrons", "0 protons"

In each column (in each block) the probabilities should sum up to 1