Tree [r195] /

File Date Author Commit
lib 2012-02-25 elseifkk [r185] bug fix
massdata 2011-10-10 elseifkk [r116] bug fixes
tools 2011-10-12 elseifkk [r118] source header changed
ChangeLog 2013-09-04 elseifkk [r195] ChangeLog format updated
README 2013-08-29 elseifkk [r194] release 0.5.4 closed
about.ui 2013-08-29 elseifkk [r194] release 0.5.4 closed
aboutcls.cpp 2012-02-22 elseifkk [r184] bug fixes
aboutcls.h 2012-01-12 elseifkk [r164] bug fix
kinemaPlot.ui 2011-12-18 elseifkk [r158] icon changed
kinemaplotcls.cpp 2012-02-22 elseifkk [r184] bug fixes
kinemaplotcls.h 2012-02-22 elseifkk [r184] bug fixes
lineread.ui 2011-12-18 elseifkk [r158] icon changed
linereadcls.cpp 2012-01-23 elseifkk [r170] bug fix
linereadcls.h 2012-01-12 elseifkk [r164] bug fix
main.cpp 2013-08-18 elseifkk [r192] release v0.5.4
mboxcls.cpp 2012-01-12 elseifkk [r164] bug fix
mboxcls.h 2012-01-12 elseifkk [r164] bug fix
mdmDlg.ui 2013-08-25 elseifkk [r193] v0.5.4 RC2
mdmcls.cpp 2013-08-25 elseifkk [r193] v0.5.4 RC2
mdmcls.h 2013-08-25 elseifkk [r193] v0.5.4 RC2
message.ui 2011-12-18 elseifkk [r158] icon changed
misc.h 2012-02-22 elseifkk [r184] bug fixes
pix.h 2012-01-12 elseifkk [r164] bug fix
rclist.ui 2011-12-18 elseifkk [r158] icon changed
rclistcls.cpp 2012-01-12 elseifkk [r164] bug fix
rclistcls.h 2012-01-12 elseifkk [r164] bug fix
relkinemaCls.cpp 2013-08-25 elseifkk [r193] v0.5.4 RC2
relkinemaCls.h 2013-08-25 elseifkk [r193] v0.5.4 RC2
relkinemaDlg.ui 2013-08-25 elseifkk [r193] v0.5.4 RC2
resultWindow.ui 2012-02-25 elseifkk [r186] bug fix
resultwindowcls.cpp 2013-08-25 elseifkk [r193] v0.5.4 RC2
resultwindowcls.h 2012-02-22 elseifkk [r184] bug fixes
rkCalcDlg.ui 2012-01-25 elseifkk [r174] layout change in mdm and clearHist button added...
rkcalccls.cpp 2013-08-25 elseifkk [r193] v0.5.4 RC2
rkcalccls.h 2012-01-25 elseifkk [r174] layout change in mdm and clearHist button added...
rwthreadcls.cpp 2013-08-18 elseifkk [r192] release v0.5.4
rwthreadcls.h 2012-02-18 elseifkk [r177] bug fix

Read Me

RelKinema version 0.5.4
========= ======= =====

Relativistic two-body kinematics calculator for nuclear reaction,
including mass data manager, table browser and plotter.


RelKinema mainly depends on the following packages:

gfortran-4.6.1 or higher

Note that it cannot build with qt4/kde4 or above.


RelKinema uses fortran codes for calculations of nuclear physics and
the configure script should know what to use as a fortran compiler.
The calculation is performed with quad-precision, which is supported
by gfortran since version 4.6, and truncated to double precision on
display. The compiler name and flags can be passed to configure script
by FC and FCFLAGS variable, respectively. For some fortran compilers
that does not support quad-precision, the following parameter might be

setenv FCFLAGS "-D_NO_REAL10_ -D_NO_REAL16_"

RelKinema depends on libfzcal for evaluating user-defined expressions
with kinematic parameters, available at:

When compiled by gfortran libfzcal will depend on libgfortran, and
path for the library might be necessary for LDFLAGS.

Note that for some systems, qt3/kde3 environment might work
unexpectedly if qt4/kde4 coexists. Coexistence of different versions
of gcc (gfortran) can also cause problem at run-time, to which we
suggest use of LD_RUN_PATH at linking-time as a workaround.

Here's an example of build steps:

setenv FC /usr/local/bin/gfortran
setenv FCFLAGS "-cpp -fcray-pointer -O2 -g0 -fdefault-integer-8"
setenv CXXFLAGS "-O2 -g0"
./configure --prefix=/usr/local --without-arts
gmake install

Atomic mass data

RelKinema depends on atomic mass data files located in a local directory.
The current version uses files extracted from a source data file available
at (as of Aug 29, 2013):


* In main window, the reaction can be set when four of three nuclei
  are specified; hit enter for completion.
* massdata/A.dat contains particle masses in MeV, which includes
  photon(g), electron(e), proton(p), neutron(n) deuteron(d), triton(t)
  and alpha(a). 
* For a pure decay process, i.e., no target, you can set "g" as target.
* In plot window, hit H key to show key-bindings.
* Some configurations are saved in ~/.relkinemarc.

Key bindings
Escape	discard keyboard input focus
A	show popup menu
M	invoke mass data manager
T	open table&plot window
K	call calculator
R	refresh
S	turn to settings page
C	clear reaction
Ctrl+O	load saved reaction
Ctrl+S	save current reaction
1,2,3	show energy, momentum and custom page in the emission box
Q	quit


We assume a nuclear reaction of type m2(m1,m3)m4', where
  m1:    projectile,
  m2:    target,
  m3:    ejectile,
  m4':   residual including excitation energy (m4+Ex).

(Lab. frame)
                  / ) th3
m1 -----------> m2 -----------------
                  \ ) th4

Each particle has four-momentum whose magnitude is denoted as
    m^2 = E^2 - p^2,

where E and p are energy and magnitude of spatial momentum,
respectively. The kinetic energy of the particle is represented as K
and the total energy of the projectile is written as,

    E = K + m.

The four-momentum in center of mass (CM) frame can be written as, for

    E3c = gamma ( E3 - beta * p3 * cos(th3)),

where gamma=1/sqrt(1-beta^2), and beta is relative speed of Lab. frame
to CM.


The on-line calculator and the custom evaluation box in the main window
accept mathematical expressions with kinematic parameters, which include:

m1,m2,m3,m4;      mass,
p1,p2,p3,p4;      spatial momentum,
E1,E2,E3,E4;      total energy,
K1,K3,K4;         kinetic energy,
p1c,p2c,p3c,p4c;  spatial momentum in CM frame,
E1c,E2c,E3c,E4c;  total energy in CM frame,
K1c,K2c,K3c,K4c;  kinetic energy in CM frame,
QValue;           Q-value of the reaction,
Ex;               excitation energy of residual,
K1Min;            minimum K1,
p1Min;            minimum p1,
K1cMin;           minimum K1 in CM frame,
p1cMin;           minimum p1 in CM frame,
theM;             the invariant mass,
gamma, beta;      CM to Lab. conversion factor,
th3, th4;         theta3 and theta4,
th3c, th4c;       theta3 and theta4 in CM frame,
q;                momentum transfer in 1/fm,
qMin;             minimum momentum transfer,
qMax;             maximum momentum transfer,
J3,J4;            solid angle Jacobian,
KS;               kinematic shift (dK3/dth3), 
KF;               kinematic factor (dp3/dth3/p3).

These parameters are case-sensitive and read-only, and passed by
reference to libfzcal in quad-precsion, if possible.

The results are always in MeV for energy and the angles in radian. 

In the calculator you can define new parameters such as, v3=p3/E3,
which can also be used in the custom evaluation box. Note that the
newly defined parameters are temporary, i.e., parameters in RHS are
passed by value. If defined as v3="p3/E3", v3 becomes a macro para-
meter and will reflect changes (parameters in RHS will be passed by 


1) L.D. Landau, E.M. Lifshitz (1975), "The Classical Theory of
Fields", Vol. 2 (4th ed.). Butterworth-Heinemann. ISBN


All reports and requests are welcome to