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.

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.

[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.