Mocadi home  |   Manual  |   Download area  |   Program  


MOCADI is a Monte Carlo simulation program to calculate the transport of primary beams, projectile fragments, and fission fragments or fusion products [4] through ion optical systems described by third or fifth order transfer matrices and through layers of matter. The ion optical matrices can be calculated with the program GICOSY. Currently, input data describing the FRS with all branches, other fragment separators and storage rings are available for different ion optical settings. Atomic interactions [1] and nuclear interactions of relativistic heavy ions with matter are calculated according to [2]. Charge state distributions can be calculated with the help of special cross-section tables [3], otherwise they have to be be given explicitly as an input parameter. They can be obtained from the programs CHARGE or GLOBAL. In the case of projectile fragments the calculation of cross sections via the EPAX formula is included. ROOT trees, CERN-standard histograms (HBOOK) and list-mode output (Ntuple) for particle mass, charge, rigidity, energy, positions, angles, TOF, and more are available at any position along the optical system. They can conveniently be analyzed with PAW or ROOT.

How to use MOCADI at GSI

MOCADI at GSI is available on the Linux system . The input parameters for MOCADI like beam properties, magnet settings and the geometric information of the experimental setup (magnets, collimators, wedges etc.) are supplied by an ASCII input file. This input file also includes references to the data files containing the matrices for the particle transfer through the various magnets. You can start with an example input file from the download area. To call the program:

/u/weick/mocadi/exe/mocadi   < filename >    (full version but without Root tree output, linked static for 32 and 64bit)
/u/weick/mocadi/exe/mocadiR < filename >    (full version with Root tree output, linked dynamically)
/u/weick/mocadi/exe/gmocadi                         (GUI for mocadi)

where < filename > is the name of the MOCADI input file. The default extension of the input file is ".in" and can be omitted. MOCADI uses a precalculated table for energy-loss and angular stragglig in matter. Two versions are available ATIMA 1.3 as default and since 2018 also the improved ATIMA 1.4. To select ATIMA 1.4 please modify the path name set in an environmant variable:
In a bash or korn shell: >export MOCADI_SPLINES_GZ='/u/weick/mocadi/splines_gz_1.4'

To obtain your own version and more details please check MOCADI program.

Distribution policy and  copyright

H. Geissel, N. Iwasa, and H. Weick hold the copyright for MOCADI, except for the subroutine for the calculation of the atomic interactions, which was made by C. Scheidenberger and P. Malzacher.
If you are interested to run MOCADI on LINUX in your institute you are welcome to contact H. Weick to obtain a copy of the executable. However, in order to prevent unauthorized modifications we decided not to distribute the source code.


As of November 1999 Helmut Weick takes care of MOCADI at GSI. Please send comments, suggestions, or bug reports to him.

History and Acknowledgement

The program was developed in PL/I by T. Schwab, H. Geissel, and A. Magel at Giessen University in Germany. N. Iwasa translated MOCADI 1.34 into C on a DEC VMS operating system and developed it further.
N. Iwasa thanks K. Sümmerer and T. Baumann for their help on writing the MOCADI manual. N. Iwasa is grateful to M. Dahlinger for his help on programming (in particular, CERN-library and C-fortran). Many thanks are due to P. Malzacher and C. Scheidenberger for providing the splines calculated by ATIMA-1.2
Since then it has been developed further on LINUX by N. Iwasa and H. Weick with many extensions and new features.


[1]  C. Scheidenberger et al., Phys. Rev. Lett. 73, 50 (1994); C. Scheidenberger et al., Phys. Rev. Lett. 77, 3987 (1996).
[2]  N. Iwasa et al., Nucl. Instr. Meth. B126 (1997) 284.
[3]  N. Iwasa, H. Weick, H. Geissel, Nucl. Instr. Meth. B269 (2011) 752.
[4]  M. Mazzocco et al., Eur. Phys. J. Special Topics 150 (2007) 363.

This page was written by Dr. Helmut Weick, 25th May 2018, contact h.weick(at),  Imprint (Impressum), Privacy Policy (Datenschutzerklärung)