STARLIGHT/starlight/.svn/text-base/Readme.svn-base

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