IAEA Benchmark of Spallation Models - Tools

You will find in this page some tools used within this benchmark. The goal is to give the possibility for those of you who want to get more figures (for example comparison of two , three, four models). The tools are given with explanation, and people can use them as they want (any question and/or remarks, contact jean-christophe.david@cea.fr). In this first version, only tools for particles are given. Tools for residue production will be provided in a second step. The tools have been made on unix (linux). Don't forget to read carefully the last section "More explanations and warnings" just after downloading the file.

Download

The first step is to download the file benchspal.tar.gz. Then, unzip (ex.: gzip -d benchspal.tar.gz) and untar (tar xvf benchspal.tar). You will get the directory benchspal.

Content of benchspal

2 directories: tools and results:

• tools contains all files to build figures and figures of merit (scripts, programs, etc, and also experimental data (expdata) and calculation results (codes)).
• results is the directory where you will work.
• In tools and results you will find several directories:
  • neutron,
  • multiplicity,
  • lcp,
  • pion.
• In each directory of results dealing with particles (so all, but latter residue and excfn should be also available), you find 2 files: a text file (input) and an executable (runbench).
• input is an input file where you decide which reaction and models you want to study.
• The first time you execute runbench (which downloads all needed files from tools and launches neutron_all or multip or lcp_all or pion_all). The next times you have just to run neutron_all or multip or lcp_all or pion_all.

The files you get are, for instance: (always: proj_target_Energy_observable_author(exp.data)_model.*; example below is for the 17th model, i.e. mcnpx-bert-dres)

p_Pb_800_ddxsn_leray_0000000000000000100000000.plt

This file is read by gnuplot to draw the figures (postscript - see file below)

p_Pb_800_ddxsn_leray_0000000000000000100000000.ps

and this ps is converted in pdf via ps2pdf

p_Pb_800_ddxsn_leray_0000000000000000100000000.pdf

Then 4 other files (tables) for the Figures of Merit

p_Pb_800_ddxsn_leray_0000000000000000100000000_fom_high.txt

p_Pb_800_ddxsn_leray_0000000000000000100000000_fom_low.txt

p_Pb_800_ddxsn_leray_0000000000000000100000000_fom_medium.txt

p_Pb_800_ddxsn_leray_0000000000000000100000000_fom_tot.txt

These 4 files can be plotted latter with the perl executable plot_fom.pl (originally in tools, but copied in results by runbench)

To use plot_fom.pl:

1. open it,
2. write on the second line the reaction to be studied ($head="reaction"; ex. reaction=p_Pb_3000_ddxsn_ishibashi),
3. put the right path for gnuplot,
4. run plot_fom.pl,
5. and follow instructions.

More explanations and warnings

• You need: fortran77, unix, perl, gnuplot, xpdf and ps2pdf.
• Use your fortran77 compiler in neutron_all or multip or lcp_all or pion_all (in tools/*).
• Use the right path in perlfix1 (in tools) for gnuplot, ps2pdf and xpdf (beginning of the file).
• Use the right path in plot_fom.pl (in tools) for gnuplot (line 8) and put the right reaction (line 1).
• In input, the numbers start column 19. The number of each reaction is given on the left column in number.html (in tools). Put 1 behind each model you want to study and 0 behind the others.
• Sometimes the scales are not the best ones, so you have to open the *.plt obtained, replace

      set autoscale 
  

  by

      set xrange [minX:maxX] 
      set yrange [minY:maxY] 
  

  and re-run "gnuplot *.plt" (you will get the new postscript file)
• We tested f77, g77, gfortran (on linux), g95 (on mac).