RP Package for Sunsort

Relies on: Coord, Vector, Reaction, Event-struct
C include file: evsubs.h
FORTRAN include files: subs.i and evsubs.i

This package handles the creation and manipulation of resonant particles. A resonant particle is formed by combining two particles to determine the properties of a potential parent.

Routines

particle *rp(particle *vrp, particle *hi, particle *li)
subroutine rp(integer vrp, integer hi, integer li)

This routine combines the two particles specified in hi and li and places the data into the particle slot specified by vrp. The C routine returns vrp.

The following parameters of the parent are calculated:

• The fold (the sum of the folds of the children).
• The mass (the sum of the masses of the children).
• The momentum (the vector sum of the momenta of the children).
• The total energy (etot) which is defined to be the sum of the particle energies of the singles particles that were combined to make this particle.
• The Q-value for the breakup of this particle to its final state (the difference of etot and the energy inferred from the particle's momentum).

particle *new_rp(particle *hi, particle *li)
integer function new_rp(integer vrp, integer hi, integer li)

This routine is the same as rp except that the parent particle is placed in the next free slot in the particle array that the particles are in (they must be in the same array).

The routine returns the ID of the new particle or NULL or -1 if there was no more room in the array.

void nrhit(double newmass, particle *p);
subroutine nrhit(double precision newmass, integer p)

This routine recalculates all resonant particles to take into account the change of mass of a specified singles particle. Note that calling evs_m directly for a particle does not do this.