class AliPythia6: public TPythia6, public AliPythiaBase

 Initialise the process to generate

 Check Lund particle code (for debugging)

 Treat protons as inside nuclei with mass numbers a1 and a2
    The MSTP array in the PYPARS common block is used to enable and
    select the nuclear structure functions.
    MSTP(52)  : (D=1) choice of proton and nuclear structure-function library
            =1: internal PYTHIA acording to MSTP(51)
            =2: PDFLIB proton  s.f., with MSTP(51)  = 1000xNGROUP+NSET
    If the following mass number both not equal zero, nuclear corrections of the stf are used.
    MSTP(192) : Mass number of nucleus side 1
    MSTP(193) : Mass number of nucleus side 2

 Set random number generator

 Print list of particl properties

  Set default values for pythia decay switches

  Set default values for pythia decay switches

  Call Pythia join alogorithm to set up a string between
  npart partons, given by indices in array ipart[npart]

  Call qPythia showering

position of the hard scattering in the nuclear overlapping area.
just for qpythia.

  Call Pythia clustering algorithm

  Call Pythia jet reconstruction algorithm

 Get jet number i

  Call Pythia showering

 Initializes
 (1) The quenching model using quenching weights according to C. Salgado and U. Wiedemann
 (2) The nuclear geometry using the Glauber Model

  Simple Jet Quenching routine:

  The jet formed by all final state partons radiated by the parton created
  in the hard collisions is quenched by a factor (1-z) using light cone variables in
  the initial parton reference frame:
  (E + p_z)new = (1-z) (E + p_z)old




   The lost momentum is first balanced by one gluon with virtuality > 0.
   Subsequently the gluon splits to yield two gluons with E = p.

 Igor Lokthine's quenching routine
 http://lokhtin.web.cern.ch/lokhtin/pyquen/pyquen.txt

 Set the parameters for the PYQUEN package.
 See comments in PyquenCommon.h

 Load event into Pythia Common Block

 New multiple interaction scenario

 Return event specific quenching parameters

 Default configuration for Heavy Flavor production

 All QCD processes

 Configuration for the ATLAS tuning

 Configuration for the ATLAS tuning

 Set the pt hard range

 Set the y hard range

 Switch fragmentation on/off

  initial state radiation

 Set the inreinsic kt

 Switch off heavy flavor
 Maximum number of quark flavours used in pdf

 Set pycell parameters

 Modified splitting probability as a model for quenching

 Switch off hadronisarion

 Switch on hadronisarion

 Get x1, x2 and Q for this event

 Get the total cross-section

 Get the pT hard for this event

 Get the subprocess code

 End of run statistics

 End of run statistics

 Assignment operator

 Copy

 convert to compressed code and print result (for debugging only)

 Treat protons as inside nuclei

 Set colliding nuclei ("p","n",...)