class AliFastGlauber: public TObject

  Default Constructor
  Defaults for Pb

 Copy constructor

 Set random number generator

 Destructor

Set all parameters for RHIC

Set all parameters for LHC

 Initialisation
    mode = 0; all functions are calculated
    mode = 1; overlap function is read from file (for Pb-Pb only)
    mode = 2; interaction almond functions are read from file
              USE THIS FOR PATH LENGTH CALC.!

 Reset dynamic allocated formulas
 in case init is called twice

  Draw Wood-Saxon Nuclear Density Function

  Draw Overlap Function

  Draw Number of Participants Npart

  Draw Thickness Function

  Draw Geometrical Cross-Section

  Draw Binary Cross-Section

  Draw Binaries per event (Ncoll)

  Draw Kernel

  Draw Interaction Almond

  Draw energy density

  Draw Path Length

  Draw Path Length

  Draw Interaction Radius

  Woods-Saxon Parameterisation
  as a function of radius (xx)

  Woods-Saxon Parameterisation
  as a function of radius (xx)
  times r**2

  Wood Saxon Parameterisation
  as a function of z and  b

  Wood Saxon Parameterisation
  as a function of z for fixed b

  Thickness function T_A
  as a function of b

  Kernel for overlap function: T_A(s)*T_A(s-b)
  as a function of r and phi

  Overlap function
  T_{AB}=Int d2s T_A(s)*T_B(s-b)
  as a function of b
 (normalized to fA*fB)

  Kernel for number of participants
  as a function of r and phi

  Number of Participants as
  a function of b

  Geometrical Cross-Section
  as a function of b

  Number of binary hard collisions
  as a function of b

  Number of hard processes per event
  as a function of b

  Initial energy density
  as a function of the impact parameter

  Almond shaped interaction region
  as a function of cartesian x,y.

  Average interaction density over radius
  at which interaction takes place
  as a function of radius

  Path Length as a function of phi
  for interaction point fixed at (0,0)
  as a function of phi-direction

  Phi direction in Almond

  Path Length as a function of phi
  Interaction point from random distribution
  as a function of the phi-direction

 Return the geometrical cross-section integrated from b1 to b2

 Return the hard cross-section integrated from b1 to b2

 Return fraction of hard cross-section integrated from b1 to b2

  Number of binary hard collisions
  as a function of b (nucl/ex/0302016 eq. 19)

 Return number of binary hard collisions normalized to 1 at b=0

 Calculate the mean overlap for impact parameter range b1 .. b2

 Calculate the mean number of collisions per event for impact parameter range b1 .. b2

 Return number of binary hard collisions at b

 Return the number of participants normalized to 1 at b=0

 Return the number of participants for impact parameter b

 Return the number of collisions for impact parameter b

 Return the number of collisions per event (at least one collision)
 for impact parameter b

  Simulates Trigger

 Gives back a random impact parameter, hard trigger probability and multiplicity

 Gives back a random impact parameter bin, and hard trigger decission

 Gives back a random impact parameter in the range bmin .. bmax

 Store in file functions

 Store in file
 40 almonds for b = (0.25+k*0.5) fm (k=0->39)

 Set limits of centrality class as fractions
 of the geomtrical cross section

 Get random impact parameter according to distribution of
 hard (binary) cross-section, in the range defined by the centrality class

 Get random position of parton production point according to
 product of thickness functions

 Get random parton azimuthal propagation direction

 Calculate path length for a parton with production point (x0,y0)
 and propagation direction (ux=cos(phi0),uy=sin(phi0))
 in a collision with impact parameter b

 Return length from random b, x0, y0, phi0
 Return also phi0

 Return length from random b, x0, y0, phi0

 Return 2 lengths back to back from random b, x0, y0, phi0
 Return also phi0

 Return 2 lengths back to back from random b, x0, y0, phi0

 Returns lenghts for n partons with azimuthal angles phi[n]
 from random b, x0, y0

 Plot distribution of n impact parameters

 Plot length distribution

 Plot lengths back-to-back distributions

 Plot almonds for some impact parameters

 Calculate integrals:
  integral0 = \int_0^ellCut dl*(T_A*T_B)(x0+l*ux,y0+l*uy)
  integral1 = \int_0^ellCut dl*l*(T_A*T_B)(x0+l*ux,y0+l*uy)

 for a parton with production point (x0,y0)
 and propagation direction (ux=cos(phi0),uy=sin(phi0))
 in a collision with impact parameter b

 Return I0 and I1 from random b, x0, y0, phi0
 Return also phi

 Return I0 and I1 from random b, x0, y0, phi0

 Return 2 pairs of I0 and I1 back to back from random b, x0, y0, phi0
 Return also phi0

 Return 2 pairs of I0 and I1 back to back from random b, x0, y0, phi0
 Return also phi0

 Return 2 pairs of I0 and I1 back to back from random b, x0, y0, phi0

 Returns I0 and I1 pairs for n partons with azimuthal angles phi[n]
 from random b, x0, y0

 Returns I0 and I1 pairs for n partons with azimuthal angles phi[n]
 from random b, x0, y0 and return x0,y0

 Plot I0-I1 distribution

 Plot I0-I1 back-to-back distributions

 Assignment operator

 Copy