Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] adlb_v1_0.tar.gz(137 Kbytes)
Manuscript Title: SOPHIA, Monte Carlo simulations of photohadronic processes in astrophysics.
Authors: A. Muecke, R. Engel, J.P. Rachen, R.J. Protheroe, T. Stanev
Program title: SOPHIA 2.0
Catalogue identifier: ADLB_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 124(2000)290
Programming language: Fortran.
Computer: DEC-Alpha workstation.
Operating system: UNIX, Linux, Open-VMS.
RAM: 1M words
Word size: 64
Keywords: Astrophysics, Nuclear processes, Photon-hadron, Interactions, Resonance production, Resonance decay, Photoproduction cross, Section, Monte Carlo event, Generator, Multiparticle production.
Classification: 1.2.

Nature of problem:
Simulation of minimum bias photoproduction for astrophysical applications.

Solution method:
Monte Carlo simulation of individual events for given nucleon and photon energies, photon energies are sampled from various distributions.

Restrictions:
Incident ambient photon fields limited to power law and blackbody spectra so far. No tests were done for center-of-mass energies sqrt(s) >= 1000 GeV.

Unusual features:
Other programs used in SOPHIA in modified form: RNDM, processor independent random number generator based on Ref [1]; JETSET 7.4, Lund Monte Carlo for jet fragmentation [2]; DECSIB, SIBYLL routine which decays unstable particles [3].

Running time:
10000 events at a center-of-mass energy of 1.5 GeV require a typical CPU time of about 75 seconds.

References:
[1] G. Marsaglia, A. Zaman, Toward a universal random number generator, Florida State Univ., Preprint FSU-SCRI-87-70, 1987.
[2] T. Sjostrand, Comp. Phys. Commun. 82 (1994) 74.
[3] R.S. Fletcher, T.K. Gaisser, P.Lipari, T. Stanev, Phys. Rev. D 50 (1994) 5710.