class PChannel: public TObject

 Channel constructor by pointer to array of pointers to particles
 (parent, decay products), number of decay particles (default 2).

 Flags: (obsolete)

 Channel constructor by decay-mode index

 Channel constructor by decay-mode index and parent reference

 Identify beam and target if the current channel is a beam + target reaction.

 Make quasi-channel for d+p reactions

 parent is midrapidity thermal source

 parent is dilepton source

On the 1st call, set the grandparent, etc...
look for tentative distributions

 The N-body phase-space decay function is based on the CERNLIB
 routine GENBOD.

 Return code:

 status = 0  all ok
        = 1  not used
        = 2  weighting negative
        = 3  not enough energy
        = 4  not enough energy to decay to final state
        = 5  genbod error (energy violation)
        = 6  genbod error (kinematics failed)
        = 7  genbod error (energy violation)
        = 8  EOF on input in PFileInput::readParticle()
        = 9  not used
        = 10 Abort distribution rejection after 10 times, do chain again
        = 11 Fermi sampling failed for deuteron wave function
 > 80 is for PReaction
        = 80 Empty event
   STATUS_NOT_DECAYED 99

 N-body phase-space decay of a particle in its rest frame,
 wrapper function to all the distribution methods

return -1 if distribution does not match the channel

returns the name from the data base

Add a bulk interface to the list
Each bulk object will be executed during the event loop
after the normal decay

Add a bulk interface to the list
Each bulk object will be executed during the event loop
before the normal decay

 Channel constructor by particle array

 Channel constructor by decay-mode index and parent reference

 N-body phase-space decay function based on the CERNLIB routine
 GENBOD. Resonance and dilepton masses (for Dalitz decays)
 are sampled in PData. Scattering angles are sampled here, for
 selected channels (see idChannel).

 address of array of pointers to the channel particles

 address of array of channel-particle ids

 number of decay products

******for quasi-elastic handling:

 1 or 2 for elementary or quasi particles respectively

returns the decay mode index of this channel

  double ds_dt(double);
 ...> moved to PParticle (IF)
 ds/dt(cos_th_cm) in the cm for N+N->N+Delta. With the mass of Delta
 sampled independently in PData, the simulated events are distributed
 as ds/dOmega (normalized). Ref: NPA459 (1986) 503

this is for envelope distributions:

 parent is midrapidity thermal source