class AliPythia: public TPythia6, public AliRndm

 Default Constructor

  Set random number

 Copy Constructor

 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
    MSTP(194) : Nuclear structure function: 0: EKS98 8:EPS08 9:EPS09LO 19:EPS09NLO

 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 clustering algorithm

  Call Pythia jet reconstruction algorithm

  Call Pythia jet reconstruction algorithm

C
C ITUNE    NAME (detailed descriptions below)
C     0 Default : No settings changed => linked Pythia version's defaults.
C ====== Old UE, Q2-ordered showers ==========================================
C   100       A : Rick Field's CDF Tune A
C   101      AW : Rick Field's CDF Tune AW
C   102      BW : Rick Field's CDF Tune BW
C   103      DW : Rick Field's CDF Tune DW
C   104     DWT : Rick Field's CDF Tune DW with slower UE energy scaling
C   105      QW : Rick Field's CDF Tune QW (NB: needs CTEQ6.1M pdfs externally)
C   106 ATLAS-DC2: Arthur Moraes' (old) ATLAS tune (ATLAS DC2 / Rome)
C   107     ACR : Tune A modified with annealing CR
C   108      D6 : Rick Field's CDF Tune D6 (NB: needs CTEQ6L pdfs externally)
C   109     D6T : Rick Field's CDF Tune D6T (NB: needs CTEQ6L pdfs externally)
C ====== Intermediate Models =================================================
C   200    IM 1 : Intermediate model: new UE, Q2-ordered showers, annealing CR
C   201     APT : Tune A modified to use pT-ordered final-state showers
C ====== New UE, interleaved pT-ordered showers, annealing CR ================
C   300      S0 : Sandhoff-Skands Tune 0
C   301      S1 : Sandhoff-Skands Tune 1
C   302      S2 : Sandhoff-Skands Tune 2
C   303     S0A : S0 with "Tune A" UE energy scaling
C   304    NOCR : New UE "best try" without colour reconnections
C   305     Old : New UE, original (primitive) colour reconnections
C   306 ATLAS-CSC: Arthur Moraes' (new) ATLAS tune (needs CTEQ6L externally)
C ======= The Uppsala models =================================================
C   ( NB! must be run with special modified Pythia 6.215 version )
C   ( available from http://www.isv.uu.se/thep/MC/scigal/        )
C   400   GAL 0 : Generalized area-law model. Old parameters
C   401   SCI 0 : Soft-Colour-Interaction model. Old parameters
C   402   GAL 1 : Generalized area-law model. Tevatron MB retuned (Skands)

 Inset 2-parton system at line idx

 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

 New multiple interaction scenario

  Call medium-modified Pythia jet reconstruction algorithm

 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

 Assignment operator

 Copy

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

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