next up previous index
Next: Y201 Find Histogram-Channel Up: CERNLIB Previous: W505 General Monte-Carlo

W515 N-Body Monte-Carlo Event Generator

Routine ID: W515
Author(s): F. JamesLibrary: POOL
Submitter: Submitted: 20.10.1975
Language: FortranRevised:

GENBOD generates a multi-particle weighted event according to Lorentz-invariant Fermi phase space. It is a modification of the routine GENEV (in FOWL (W505)) and uses the method of Raubold and Lynch (see Ref. 1). The total CM energy as well as the number and masses of the outgoing particles are specified by the user, but may be changed from event to event. GENBOD generates the CM vector momenta (and energies) of the outgoing particles and gives the weight which must be associated with each event. The weight may then be multiplied by any 'matrix element' or geometrical detection function calculated by the user.

Structure:

SUBROUTINE subprogram
User Entry Names: GENBOD
Files Referenced: Printer
External References: FLPSOR (M103), RNDM (V104), PDK, ROTES2 (W505)
COMMON Block Names and Lengths : /GENIN/ 21, /GENOUT/ 91

Usage:

    COMMON /GENIN /NP,TECM,AMASS(18),KGENEV
    COMMON /GENOUT/ PCM(5,18),WT
    CALL GENBOD
Input:
NP
( INTEGER) Number of outgoing particles (2 ≤NP ≤18) .
TECM
( REAL) Total CM energy.
AMASS
( REAL) Array where element I contains the mass of the I-th outgoing particle.
KGENEV
( INTEGER) = 1 for cross section constant with energy, = 2 for Fermi energy dependence.
Output:
PCM(1,I)
( REAL) Px of I-th particle.
PCM(2,I)
( REAL) Py of I-th particle.
PCM(3,I)
( REAL) Pz of I-th particle.
PCM(4,I)
( REAL) Energy of I-th particle.
PCM(5,I)
( REAL) P of I-th particle.
WT
( REAL) Weight of the event.
See also the Long Write-up for FOWL (W505).

References:

  1. F. James, Monte Carlo Phase Space, CERN 68-15 (1968)

Statistical Data Analysis and Presentation

Y201



next up previous index
Next: Y201 Find Histogram-Channel Up: CERNLIB Previous: W505 General Monte-Carlo


Janne Saarela
Mon Apr 3 15:06:23 METDST 1995