### 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: stadard deviation (of average multiplicity) of "ALL" particles
- #4: average multiplicity of "ALL CHARGED" particles
- #5: stadard deviation (of average multiplicity) of "ALL CHARGED" particles
- #6: average multiplicity of "charge -1" particles
- #7: stadard deviation (of average multiplicity) of "charge -1" particles
- #8: average multiplicity of "charge 0" particles
- #9: stadard deviation (of average multiplicity) of "charge 0" particles
- #10: average multiplicity of "charge +1" particles
- #11: stadard deviation (of average multiplicity) of "charge +1" particles
- #12: average multiplicity of "ALL pions"
- #13: stadard deviation (of average multiplicity) of "ALL pions"
- #14: average multiplicity of "pion -1"
- #15: stadard deviation (of average multiplicity) of "pion -1"
- #16: average multiplicity of "pion 0"
- #17: stadard deviation (of average multiplicity) of "pion 0"
- #18: average multiplicity of "pion +1"
- #19: stadard deviation (of average multiplicity) of "pion +1"
- #20: average multiplicity of "ALL nucleons"
- #21: stadard deviation (of average multiplicity) of "ALL nucleons"
- #22: average multiplicity of "neutrons"
- #23: stadard deviation (of average multiplicity) of "neutrons"
- #24: average multiplicity of "protons"
- #25: stadard deviation (of average multiplicity) of "protons"
- #26: average multiplicity of "ALL kaons"
- #27: stadard deviation (of average multiplicity) of "ALL kaons"
- #28: average multiplicity of "Kbar -1"
- #29: stadard deviation (of average multiplicity) of "Kbar -1"
- #30: average multiplicity of "Kbar 0"
- #31: stadard deviation (of average multiplicity) of "Kbar 0"
- #32: average multiplicity of "K 0"
- #33: stadard deviation (of average multiplicity) of "K 0"
- #34: average multiplicity of "K +1"
- #35: stadard deviation (of average multiplicity) of "K +1"
- #36: average multiplicity of "Lambda"
- #37: stadard 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 each value of "running variable" you have several (="num_runs_sameEnergy") blocks of output separated by two empty lines

Columns:

- #1: number of outgoing particles (integer)
- #2: percentage of events with number of "ALL" particles equal to #1
- #3: percentage of events with number of "ALL CHARGED" particles equal to #1
- #4: percentage of events with number of "charge -1" particles equal to #1
- #5: percentage of events with number of "charge 0" particles equal to #1
- #6: percentage of events with number of "charge +1" particles equal to #1
- #7: percentage of events with number of "ALL pions" equal to #1
- #8: percentage of events with number of "pion -1" equal to #1
- #9: percentage of events with number of "pion 0" equal to #1
- #10: percentage of events with number of "pion +1" equal to #1
- #11: percentage of events with number of "ALL baryons" equal to #1
- #12: percentage of events with number of "neutrons" equal to #1
- #13: percentage of events with number of "protons" equal to #1
- #14: percentage of events with number of "ALL kaons" (and antikaons) equal to #1
- #15: percentage of events with number of "Kbar -1" equal to #1
- #16: percentage of events with number of "Kbar 0" equal to #1
- #17: percentage of events with number of "K0" equal to #1
- #18: percentage of events with number of "K+" equal to #1
- #19: percentage of events with number of "Lambda" equal to #1
- #20: average multiplicity of "ALL" (here column 1 makes not sense, so the values are the same)
- #21: average multiplicity of "ALL CHARGED" (here column 1 makes not sense, so the values are the same)
- #22: average multiplicity of "charge -1" (here column 1 makes not sense, so the values are the same)
- #23: average multiplicity of "charge 0" (here column 1 makes not sense, so the values are the same)
- #24: average multiplicity of "charge +1" (here column 1 makes not sense, so the values are the same)
- #25: average multiplicity of "ALL pions" (here column 1 makes not sense, so the values are the same)
- #26: average multiplicity of "pion -1" (here column 1 makes not sense, so the values are the same)
- #27: average multiplicity of "pion 0" (here column 1 makes not sense, so the values are the same)
- #28: average multiplicity of "pion +1" (here column 1 makes not sense, so the values are the same)
- #29: average multiplicity of "ALL baryons" (here column 1 makes not sense, so the values are the same)
- #30: average multiplicity of "neutrons" (here column 1 makes not sense, so the values are the same)
- #31: average multiplicity of "protons" (here column 1 makes not sense, so the values are the same)
- #32: average multiplicity of "ALL kaons" (here column 1 makes not sense, so the values are the same)
- #33: average multiplicity of "Kbar -1" (here column 1 makes not sense, so the values are the same)
- #34: average multiplicity of "Kbar 0" (here column 1 makes not sense, so the values are the same)
- #35: average multiplicity of "K0" (here column 1 makes not sense, so the values are the same)
- #36: average multiplicity of "K+" (here column 1 makes not sense, so the values are the same)
- #37: average multiplicity of "Lambda" (here column 1 makes not sense, so the values are the same)

NOTES

The numbers in one raw 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