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 OUTPUT_FORMAT = 2    #Form of the output (not considered at this state)
  53
  54 BREAKUP_MODE = 5     #Controls the nuclear breakup
  55                      # Note that this option only works for lead or gold; it should work at any energy
  56                      # 1 = hard sphere nuclei (b>2R)
  57                      # 2 = require that both nuclei break up (Xn,Xn)
  58                      # 3 = a single neutron from each nucleus (1n,1n)
  59                      # 4 = require that neither nucleus breaks up (0n,0n)
  60                      # 5 = require that there be no hadronic break up
  61                      #     (This is similar to option 1, but with the actual hadronic interaction
  62                      #     probability)
  63                      # 6 = require break up of one or both nuclei (Xn)
  64                      # 7 = require break up of one nucleus (Xn,0n)
  65
  66 INTERFERENCE = 0     #Interference (0 = off, 1 = on)
  67 IF_STRENGTH = 1.     #% of intefernce (0.0 - 0.1)
  68 COHERENT = 1         #Coherent=1,Incoherent=0 production for assym collisions
  69 INCO_FACTOR = 1.     #percentage of incoherence
  70 BFORD = 9.5          #Deuteron slope parameter (dsig/dt = exp(BFORD*t))
  71 INT_PT_MAX = 0.24    #Maximum pt considered, when interference is turned on
  72 INT_PT_N_BINS =120   #Number of pt bins when interference is turned on
  73
  74
  75 ----------------------------------------------------------------------------
  76
  77 Currently supported 2-photon channel options:
  78       jetset id         particle
  79    ---------------------------------
  80         221             eta
  81         331             eta-prime
  82         441             eta-c
  83         9010221         f0(975)
  84         225             f2(1270)
  85         115             a2(1320)
  86         335             f2(1525)
  87         33              rho0 pair
  88         11              e+/e- pair
  89         13              mu+/mu- pair
  90         15              tau+/tau- pair
  91
  92 Currently supported vector meson options:
  93       jetset id         particle
  94    ---------------------------------
  95
  96         113             rho0
  97         223             omega
  98         333             phi
  99         443011          J/psi --> e+e-
 100         443013          J/Psi --> mu+mu-
 101         444011          Psi(2S) --> e+e-
 102         444013          Psi(2S) --> mu+mu-
 103         553011          Upsilon(1S) --> e+e-
 104         553013          Upsilon(1S) --> mu+mu-
 105         554011          Upsilon(2S) --> e+e-
 106         554013          Upsilon(2S) --> mu+mu-
 107         555011          Upsilon(3S) --> e+e-
 108         555013          Upsilon(3S) --> mu+mu-
 109         913             rho0 + direct pi+pi- (with interference)
 110 The direct pi+pi- fraction is from the ZEUS results
 111
 112
 113
 114 Logs:
 115
 116 fixed types for class Particle
 117 - Int -> Double
 118 - fixed Makefile ..
 119    autoconfig will appear in the future .. make sure paths are fixed before compiling
 120    this is for includes and libraries in Makefile
 121 - don't forget about link to pythia xml doc
 122   This is located in main.cpp
 123
 124 - changed PID of the charged KAON to 321 and neutral kaon to 310 .. it should work now for phi -> kk
 125
 126 - incoherent formfactor was hardcoded for dAu events, it should now accept coherent too.
 127 - Had hardcoded narrow and wide resonances for VMs. The user can again select gg_or_gp (2 or 3)
 128 except for the psi family.  This is still always narrow.
 129 - Changed the spin dependence for VMs going to ee/mumu(gammaavectormeson::gettheta).  It was only for psi and now upsilon will follow too
 130 inretrospect, i should probably have done this first and merge psifamily back into gammaavm.cpp