From: Giuseppe Battistoni (Giuseppe.Battistoni@mi.infn.it)
Date: Wed Jul 18 2007 - 17:01:14 CEST
Yes: for a pure phase-space decay there is SUBROUTINE EXPLOD(...)
Its description and instructions for use are the following:
SUBROUTINE EXPLOD ( NPEXPL, AMEXPL, ETOTEX, ETEXPL, PXEXPL,
& PYEXPL, PZEXPL )
*----------------------------------------------------------------------*
* *
* Copyright (C) 1995-2006 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* EXPLOD: makes a phase space like explosion out of Npexpl parti- *
* cles of masses Amexpl and total energy Etotex. The *
* resulting total energies and momenta are stored into *
* etexpl and px,y,zexpl. The "explosion" is supposed to *
* occur in the CMS of the particles *
* *
* Created on 08 february 1995 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 28-may-06 by Alfredo Ferrari *
* *
* Input variables: *
* *
* Npexpl = number of particles *
* Amexpl(i) = mass of the ith particle, i=1,npexpl *
* Etotex = invariant mass of npexpl particle system *
* *
* Output variables: *
* *
* Etexpl(i) = (total) energy of the ith particle in the *
* CMS of the npexpl particle system *
* Px,y,zexpl(i) = momentum components of the ith particle in *
* the CMS of the npexpl particle system *
* *
* The unit of energy and momentum is immaterial provided Etotex *
* and Amexpl are given in the same units. Of course momenta come *
* out in the same unit too *
* *
* Limitations: the number of particles must not exceed Ncycmx+2 *
* The efficiency becomes smaller the larger the par- *
* ticle number and the larger the system energy with *
* respect to the total particle mass *
* *
Best regards
Giuseppe Battistoni
On Wed, 18 Jul 2007, Biagio Di Micco wrote:
> Dear all,
>
> there is a multiparticle phase-space decay routine in fluka?
> I need to implement a decay in the source routine and I would not like to
> use external generator to not lose random number seed sequence.
>
> thanks in advance,
>
> Biagio
>
>
This archive was generated by hypermail 2.1.6 : Wed Jul 18 2007 - 17:19:54 CEST