class PReaction: public TObject

 Reaction constructor by pointer to array of pointers to Channels,
 output file name (see below), options (see below), pointer to external event tree
 (if called from the decay manager).
 Options:
   f0: 0=Tracked, 1=All particles in ROOT output file_name.root
   f1: unused
   f2: 0=Do not calculate, 1=Calculate production vertices for product particles
   f3: 0=Do not produce, 1=Produce GEANT ascii output file_name.evt (tracked particles always)

   Format of the ASCII output in file_name.evt:
   EvSeqNb NumPart Ebeam bPar phiEv flag <== number of particles in event
   E1 px1  py1  pz1  ID1  weight1      <-- 3-momentum components in GeV/c,
   E2 px2  py2  pz2  ID2  weight2        <-- GEANT particle ID,
   E3 px3  py3  pz3  ID3  weight3          <-- weight
   .    .    .    .   .    .   .   .   .   .  .   .

 See also PProjector/PBatch for additional conventions if filters are used.

 same as above, but passing one unsigned int as flag
 position of arguments shuffled to break ambiguity with previous constructor

 build reaction completely from descriptor string

 build reaction completely from descriptor string

 Build "reaction" completely from descriptor string  (author: V. Hejny)
 For "command, the batch script syntax is supported
 to set the beam parameters. The following variables
 can be used:
 _T1, _T2, _P1, _P2 as kinetic beam energy/momentum
 for the beam/target (in GeV), and
 _theta1, _theta2, _phi as beam inclination (in rad).
 It should be mentioned, that for the target, the
 notation is opposite, i.e. theta2 has the opposite
 rotation as theta1, and T2 goes into the opposite
 direction as T1

 If no variable is used, only T1 as the beam energy is
 assumed

 Examples:
 "2.2"               -> Beam energy is 2.2 GeV
 "_T1=1.5; _T2=1.2"  -> Collider experiment
 "_T1=1.5; _theta=2*TMath::DegToRad();"
                     -> Beam inclination of 2 deg

 build reaction completely from descriptor string  (author: V. Hejny)

 Reaction destructor

    cout << name << endl;

 parse reaction descriptor string  (author: V. Hejny)

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

 compute the reaction id from the list of product ids of 1st channel
 (to be used in SetUp() as sourceId)

 get the channels and particles, identify the physics, set up the defaults.

 Sets up output file names

 resets the dynamical objects of the reaction, creates tree and opens output files

 Simulate "nevents" events

outdated

 close root output file

kept for backward compatibility

 number of events to simulate

 returns number of all the particles in the reaction

 returns product channel index

 turn off vertex calculation after initialization

 put all particles on root file after initialization

 put tracked particles on root file file after initialization

 enable ascii output for GEANT after initialization

 disable ascii output for GEANT after initialization

 sets up output file names

 This function is invoked by the Filter class. It directs loop to
 impose a filter at address (*filter) following the decay represented
 by the reaction Channel number ch_id. See also PFilter.

 close root output file

 from the HGeant prompt

pure inline PReaction, e.g. FAIRROOT

 resets the dynamical objects of the reaction and opens files

 ascii file