RE: tab.lis file different to plot, what are the correct units

From: Joachim Vollaire <>
Date: Tue, 11 May 2010 23:14:23 +0200

Dear Tristan

I think that the conversion factor between what is being plotted by
flair (what you see on gnuplot) and the numerical value in the output
file is equal to (Emax-Emin)*(log (Emax-Emin) ), where Emax and Emin are
respectively the higher and lower boundaries of each bin. Please check
and let me know if this does not explain the difference you are
observing. Hope it helps.


From: on behalf of
Sent: Tue 4/27/2010 4:07 PM
Subject: tab.lis file different to plot, what are the correct units


I have used a usryield scoring card to try and calculate pion+/- yield
from a target. I used usryield because I want to make a cut in
transverse momentum (only score pions that have a transverse momentum
below 0.5GeV/c) with energies in a defined range. I have specified
13bins in the energy range between 1.5GeV and 8GeV, i.e. 0.5GeV bin

I looked at the tab.lis file that is created and I also plotted the
result using FLAIR and gnuplot. I found that the numbers in the plot
were higher than the values in the tab.lis file. I then found that the
numbers in the plot are a transformation of the numbers in the file.
Each number in the file has been multiplied by the energy in the centre
of the bin.

I have not multiplied by the surface area of the target or 2pi() at this
stage. Does that mean the numbers in the file have units of
pions/GeV/primary/cm2/sr and the numbers on the graph have units of
pions/primary/cm2/sr? I would appreciate confirmation on these units
from an expert as I was a little surprised to find the numbers in the
tab.lis file were different to the numbers on the plot.

I attach my input file if that is of any help.

Many Thanks


Tristan Davenne
STFC Technology
G19 R65
Rutherford Appleton Laboratory
Chilton, Didcot, Oxfordshire OX11 0QX, UK
Tel: +44 (0)1235 445485
Email: <>
Received on Wed May 12 2010 - 09:16:44 CEST

This archive was generated by hypermail 2.2.0 : Wed May 12 2010 - 09:16:46 CEST