HOWTO: Report sqrt(s), etc of collisions per time step

The module collisionReporter is prepared to collect data about (two-body-) collisions.

source:/workingCode/code/collisions/collisionReporter.f90

It provides routines to store and to print the number of colisions as function of time and of sqrt(s). All information is split up into several chanels, cf. fritiofType:

• 07: BaB (Baryon-Antibaryon-Annihilation)
• 08: Manni (Meson-Baryon-Annihilation)
• 09: Elastic
• 10: == Low Energy ==
• 11: FRITIOF
• 12: PYTHIA

With gnuplot-scripts like the following...

reset
set xlabel "t [fm]
set ylabel "sqrt(s) [GeV]"
set xrange [0:]
set yrange [1:]
set log zcb
#splot  "CollRep.007.dat" u 1:2:($3>1e-10?$3:1/0) w l lt 2 t "BaB"
splot "CollRep.008.dat" u 1:2:($3>1e-10?$3:1/0) w l lt 1 t "Manni"
#replot "CollRep.009.dat" u 1:2:($3>1e-10?$3:1/0) w l lt 5 t "Elast"
replot "CollRep.010.dat" u 1:2:($3>1e-10?$3:1/0) w l lt 3 t "low energy"
replot "CollRep.012.dat" u 1:2:($3>1e-10?$3:1/0) w l lt 4 t "Pythia"

...you can generate figures like this:

These are examples taken from electron induced reactions according Hermes conditions (E_e=27.3GeV)