gibuu is hosted by Hepforge, IPPP Durham
GiBUU

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:

/chrome/site/files/CollRep1.png /chrome/site/files/CollRep2.png

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

Last modified 6 years ago Last modified on Jul 19, 2013, 4:49:36 PM