STARLIGHT/starlight/Readme

   1 The STARlight Monte Carlo models 2-photon and photon-Pomeron interactions
   2 in ultra-peripheral heavy ion collisions. The physics approach for the photon-Pomeron
   3 interactions is described in Klein and Nystrand, Phys. Rev. C60, 014903 (1999),
   4 with the p_t spectrum (including vector meson interference) discussed in Phys. Rev.
   5 Lett. 84, 2330 (2000). The 2-photon interactions are described in Baltz, Gorbunov,
   6 Klein, Nystrand, Phys.Rev. C80 (2009) 044902.
   7
   8 STARlight has several input files, all of which are expected to be in the
   9 same directory as the starlight code.  User-specified input parameters are
  10 read from a file named "slight.in"; these parameters are described
  11 below.
  12
  13 The output of the STARlight Monte Carlo will be written to an ascii file
  14 named "slight.out".
  15
  16 ----------------------------------------------------------------------------
  17
  18 The contents of a sample slight.in file are listed below, with
  19 explanation.  The user should modify this file according to his or her
  20 needs.
  21 BEAM_1_Z = 82     #Z of projectile
  22 BEAM_1_A = 208    #A of projectile
  23 BEAM_2_Z = 82     #Z of target
  24 BEAM_2_A = 208    #A of target
  25 BEAM_GAMMA = 1470 #Gamma of the colliding ions
  26
  27 W_MAX = 12.0   #Max value of w
  28 W_MIN = -1     #Min value of w
  29 W_N_BINS = 40  #Bins i w
  30         # maximum and minimum values for w (the gamma-gamma center of mass
  31         # energy, w = 4(E1)(E2), (Setting W_MIN = -1 tells STARlight to use the default values
  32         # specified in inputParameters.cpp; otherwise, specify wmin here, and the number
  33         # of w bins in the lookup tables
  34
  35 RAP_MAX = 8.    # max rapidity
  36 RAP_N_BINS = 80 # Number of y bins used in the cross section calculation
  37
  38 CUT_PT = 0       #Cut in pT? 0 = (no, 1 = yes)
  39 PT_MIN = 1.0     #Minimum pT in GeV
  40 PT_MAX = 3.0     #Maximum pT in GeV
  41 CUT_ETA = 0      #Cut in pseudorapidity? (0 = no, 1 = yes)
  42 ETA_MIN = -10    #Minimum pseudorapidity
  43 ETA_MAX = 10     #Maximum pseudorapidity
  44
  45 PROD_MODE = 2    #gg or gP switch (1 = two-photon, 2 = coherent vector meson, 4 = incoherent vector meson )
  46
  47 N_EVENTS = 10    #Number of events produced
  48
  49 PROD_PID = 443013   #Channel of interest (in PDG notation); currently supported options listed below
  50 RND_SEED = 34533    #Random number seed
  51
  52 BREAKUP_MODE = 5     #Controls the nuclear breakup
  53                      # Note that this option only works for lead or gold; it should work at any energy
  54                      # 1 = hard sphere nuclei (b>2R)
  55                      # 2 = require that both nuclei break up (Xn,Xn)
  56                      # 3 = a single neutron from each nucleus (1n,1n)
  57                      # 4 = require that neither nucleus breaks up (0n,0n)
  58                      # 5 = require that there be no hadronic break up
  59                      #     (This is similar to option 1, but with the actual hadronic interaction
  60                      #     probability)
  61                      # 6 = require break up of one or both nuclei (Xn)
  62                      # 7 = require break up of one nucleus (Xn,0n)
  63
  64 INTERFERENCE = 0     #Interference (0 = off, 1 = on)
  65 IF_STRENGTH = 1.     #% of intefernce (0.0 - 0.1)
  66 INT_PT_MAX = 0.24    #Maximum pt considered, when interference is turned on
  67 INT_PT_N_BINS =120   #Number of pt bins when interference is turned on
  68 COHERENT = 1         #Coherent=1,Incoherent=0 production for assym collisions
  69 INCO_FACTOR = 1.     #percentage of incoherence
  70
  71
  72 ----------------------------------------------------------------------------
  73
  74 Currently supported 2-photon channel options:
  75       jetset id         particle
  76    ---------------------------------
  77         221             eta
  78         331             eta-prime
  79         441             eta-c
  80         9010221         f0(975)
  81         225             f2(1270)
  82         115             a2(1320)
  83         335             f2(1525)
  84         33              rho0 pair
  85         11              e+/e- pair
  86         13              mu+/mu- pair
  87         15              tau+/tau- pair
  88
  89 Currently supported vector meson options:
  90       jetset id         particle
  91    ---------------------------------
  92
  93         113             rho0
  94         223             omega
  95         333             phi
  96         443011          J/psi --> e+e-
  97         443013          J/Psi --> mu+mu-
  98         444011          Psi(2S) --> e+e-
  99         444013          Psi(2S) --> mu+mu-
 100         553011          Upsilon(1S) --> e+e-
 101         553013          Upsilon(1S) --> mu+mu-
 102         554011          Upsilon(2S) --> e+e-
 103         554013          Upsilon(2S) --> mu+mu-
 104         555011          Upsilon(3S) --> e+e-
 105         555013          Upsilon(3S) --> mu+mu-
 106         913             rho0 + direct pi+pi- (with interference)
 107 The direct pi+pi- fraction is from the ZEUS results
 108
 109
 110
 111 Logs:
 112
 113 fixed types for class Particle
 114 - Int -> Double
 115 - fixed Makefile ..
 116    autoconfig will appear in the future .. make sure paths are fixed before compiling
 117    this is for includes and libraries in Makefile
 118 - don't forget about link to pythia xml doc
 119   This is located in main.cpp
 120
 121 - changed PID of the charged KAON to 321 and neutral kaon to 310 .. it should work now for phi -> kk
 122