[u/mrichter/AliRoot.git] / STARLIGHT / starlight / .svn / text-base / Readme.svn-base

The STARlight Monte Carlo models 2-photon and photon-Pomeron interactions\r
in ultra-peripheral heavy ion collisions. The physics approach for the photon-Pomeron \r
interactions is described in Klein and Nystrand, Phys. Rev. C60, 014903 (1999),\r
with the p_t spectrum (including vector meson interference) discussed in Phys. Rev. \r
Lett. 84, 2330 (2000). The 2-photon interactions are described in Baltz, Gorbunov, \r
Klein, Nystrand, Phys.Rev. C80 (2009) 044902. \r
\r
STARlight has several input files, all of which are expected to be in the\r
same directory as the starlight code.  User-specified input parameters are\r
read from a file named "slight.in"; these parameters are described\r
below.\r
\r
The output of the STARlight Monte Carlo will be written to an ascii file\r
named "slight.out".\r
\r
----------------------------------------------------------------------------\r
\r
The contents of a sample slight.in file are listed below, with \r
explanation.  The user should modify this file according to his or her\r
needs.\r
BEAM_1_Z = 82	  #Z of projectile\r
BEAM_1_A = 208    #A of projectile\r
BEAM_2_Z = 82     #Z of target\r
BEAM_2_A = 208 	  #A of target\r
BEAM_GAMMA = 1470 #Gamma of the colliding ions\r
\r
W_MAX = 12.0   #Max value of w\r
W_MIN = -1     #Min value of w\r
W_N_BINS = 40  #Bins i w\r
      	# maximum and minimum values for w (the gamma-gamma center of mass\r
	# energy, w = 4(E1)(E2), (Setting W_MIN = -1 tells STARlight to use the default values\r
	# specified in inputParameters.cpp; otherwise, specify wmin here, and the number\r
        # of w bins in the lookup tables\r
\r
RAP_MAX = 8.    # max rapidity\r
RAP_N_BINS = 80 # Number of y bins used in the cross section calculation\r
\r
CUT_PT = 0   	 #Cut in pT? 0 = (no, 1 = yes)\r
PT_MIN = 1.0 	 #Minimum pT in GeV\r
PT_MAX = 3.0 	 #Maximum pT in GeV\r
CUT_ETA = 0  	 #Cut in pseudorapidity? (0 = no, 1 = yes)\r
ETA_MIN = -10 	 #Minimum pseudorapidity\r
ETA_MAX = 10 	 #Maximum pseudorapidity\r
\r
PROD_MODE = 2    #gg or gP switch (1 = two-photon, 2 = coherent vector meson, 4 = incoherent vector meson )\r
\r
N_EVENTS = 10    #Number of events produced\r
\r
PROD_PID = 443013   #Channel of interest (in PDG notation); currently supported	options listed below\r
RND_SEED = 34533    #Random number seed\r
\r
OUTPUT_FORMAT = 2    #Form of the output (not considered at this state)\r
\r
BREAKUP_MODE = 5     #Controls the nuclear breakup\r
	             # Note that this option only works for lead or gold; it should work at any energy\r
		     # 1 = hard sphere nuclei (b>2R)\r
		     # 2 = require that both nuclei break up (Xn,Xn)\r
		     # 3 = a single neutron from each nucleus (1n,1n)\r
		     # 4 = require that neither nucleus breaks up (0n,0n) \r
		     # 5 = require that there be no hadronic break up \r
		     # 	   (This is similar to option 1, but with the actual hadronic interaction\r
		     # 	   probability)\r
		     # 6 = require break up of one or both nuclei (Xn)\r
                     # 7 = require break up of one nucleus (Xn,0n)\r
\r
INTERFERENCE = 0     #Interference (0 = off, 1 = on)\r
IF_STRENGTH = 1.     #% of intefernce (0.0 - 0.1)\r
COHERENT = 1         #Coherent=1,Incoherent=0 production for assym collisions\r
INCO_FACTOR = 1.     #percentage of incoherence\r
BFORD = 9.5          #Deuteron slope parameter (dsig/dt = exp(BFORD*t))\r
INT_PT_MAX = 0.24    #Maximum pt considered, when interference is turned on\r
INT_PT_N_BINS =120   #Number of pt bins when interference is turned on\r
\r
\r
----------------------------------------------------------------------------\r
\r
Currently supported 2-photon channel options:\r
      jetset id		particle\r
   ---------------------------------\r
	221		eta			\r
	331		eta-prime 		\r
	441		eta-c			\r
	9010221		f0(975)		\r
	225		f2(1270)\r
	115		a2(1320)\r
	335		f2(1525)\r
	33		rho0 pair\r
	11		e+/e- pair\r
	13 		mu+/mu- pair\r
	15		tau+/tau- pair\r
\r
Currently supported vector meson options:\r
      jetset id		particle\r
   ---------------------------------\r
\r
	113		rho0\r
	223		omega\r
	333		phi\r
	443011		J/psi --> e+e-\r
        443013          J/Psi --> mu+mu- \r
	444011		Psi(2S) --> e+e-\r
	444013		Psi(2S) --> mu+mu- \r
	553011 		Upsilon(1S) --> e+e- \r
	553013 		Upsilon(1S) --> mu+mu-  \r
	554011 		Upsilon(2S) --> e+e- \r
	554013 		Upsilon(2S) --> mu+mu-  \r
	555011 		Upsilon(3S) --> e+e- \r
	555013 		Upsilon(3S) --> mu+mu-  \r
	913		rho0 + direct pi+pi- (with interference)\r
The direct pi+pi- fraction is from the ZEUS results\r
\r
\r
\r
Logs:\r
\r
fixed types for class Particle \r
- Int -> Double \r
- fixed Makefile .. \r
   autoconfig will appear in the future .. make sure paths are fixed before compiling \r
   this is for includes and libraries in Makefile\r
- don't forget about link to pythia xml doc \r
  This is located in main.cpp\r
\r
- changed PID of the charged KAON to 321 and neutral kaon to 310 .. it should work now for phi -> kk \r
\r
- incoherent formfactor was hardcoded for dAu events, it should now accept coherent too.\r
- Had hardcoded narrow and wide resonances for VMs. The user can again select gg_or_gp (2 or 3)\r
except for the psi family.  This is still always narrow.\r
- Changed the spin dependence for VMs going to ee/mumu(gammaavectormeson::gettheta).  It was only for psi and now upsilon will follow too\r
inretrospect, i should probably have done this first and merge psifamily back into gammaavm.cpp\r
Commit	Line	Data
da32329d AM	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