4 #ifndef __CFORTRAN_LOADED
12 /*=========================================================*/
13 /* COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)*/
14 /*---------------------------------------------------------*/
22 #define HIPARNT COMMON_BLOCK(HIPARNT,hiparnt)
23 COMMON_BLOCK_DEF(HiparntCommon,HIPARNT);
25 /**************************************************************************/
26 /* D E S C R I P T I O N : */
27 /*------------------------------------------------------------------------*/
28 /*COMMON/HIPARNT/HIPR1(100), IHPR2(50), HINT1(100), IHNT2(50) */
29 /*Purpose: contains input parameters (HIPR1, IHPR2) for event options */
30 /* and some extra information (HINT1, IHNT2) of current event. */
31 /*HIPR1(1): (D=1.5 GeV/$c^2$) minimum value for the invariant mass of */
32 /* the excited string system in a hadron-hadron interaction. */
33 /*HIPR1(2): (D=0.35 GeV) width of the Gaussian $P_T$ distribution of */
34 /* produced hadron in Lund string fragmentation */
35 /* (PARJ(21) in JETSET 7.2). */
36 /*HIPR1(3), HIPR1(4): (D=0.5, 0.9 GeV$^{-2}$) give the $a$ and $b$ */
37 /* parameters of the symmetric Lund fragmentation function */
38 /* (PARJ(41), PARJ(42) in JETSET 7.2). */
39 /*HIPR1(5): (D=2.0 GeV/$c^2$) invariant mass cut-off for the dipole */
40 /* radiation of a string system below which soft gluon */
41 /* radiations are terminated. */
42 /*HIPR1(6): (D=0.1) the depth of shadowing of structure functions */
43 /* at $x=0$:$\alpha_A=\mbox{HIPR1(6)}\times(A^{1/3}-1)$. */
44 /*HIPR1(7): not used */
45 /*HIPR1(8): (D=2.0 GeV/$c$) minimum $P_T$ transfer in hard or */
46 /* semihard scatterings. */
47 /*HIPR1(9): (D=$-1.0$ GeV/$c$) maximum $P_T$ transfer in hard or */
48 /* semihard scatterings. If negative, the limit is set */
49 /* by the colliding energy. */
50 /*HIPR1(10): (D=$-2.25$ GeV/$c$) specifies the value of $P_T$ for */
51 /* each triggered hard scattering generated per event */
52 /* (see Section \ref{sec:jet2}). If HIPR1(10) is negative, */
53 /* its absolute value gives the low limit of the */
54 /* $P_T$ of the triggered jets. */
55 /*HIPR1(11): (D=2.0 GeV/$c$) minimum $P_T$ of a jet which will interact */
56 /* with excited nuclear matter. When the $P_T$ of a jet */
57 /* is smaller than HIPR1(11) it will stop interacting further. */
58 /*HIPR1(12): (D=1.0 fm) transverse distance between a traversing jet */
59 /* and an excited nucleon (string system) below which they */
60 /* will interact and the jet will lose energy and momentum */
61 /* to that string system. */
62 /*HIPR1(13): (D=1.0 fm) the mean free path of a jet when it goes */
63 /* through the excited nuclear matter. */
64 /*HIPR1(14): (D=2.0 GeV/fm) the energy loss $dE/dz$ of a gluon */
65 /* jet inside the excited nuclear matter. The energy loss */
66 /* for a quark jet is half of the energy loss of a gluon. */
67 /*HIPR1(15): (D=0.2 GeV/$c$) the scale $\Lambda$ in the */
68 /* calculation of $\alpha_s$. */
69 /*HIPR1(16): (D=2.0 GeV/$c$) the initial scale $Q_0$ for the */
70 /* evolution of the structure functions. */
71 /*HIPR1(17): (D=2.0) $K$ factor for the differential jet cross */
72 /* sections in the lowest order pQCD calculation. */
73 /*HIPR1(18): not used */
74 /*HIPR1(19), HIPR1(20): (D=0.1, 1.4 GeV/$c$) parameters in the */
75 /* distribution for the $P_T$ kick from soft interactions , */
76 /* $1/[(\mbox{HIPR1(19)}^2+P_T^2)(\mbox{HIPR1(20)}^2+P_T^2)]$. */
77 /*HIPR1(21): (D=1.6 GeV/$c$) the maximum $P_T$ for soft interactions, */
78 /* beyond which a Gaussian distribution as specified by */
79 /* HIPR1(2) will be used. */
80 /*HIPR1(22): (D=2.0 GeV/$c$) the scale in the form factor to suppress */
81 /* the $P_T$ transfer to diquarks in hard scatterings, */
82 /*HIPR1(23)--HIPR1(28): not used. */
83 /*HIPR1(29): (D=0.4 fm) the minimum distance between two nucleons */
84 /* inside a nucleus when the coordinates of the nucleons */
85 /* inside a nucleus are initialized. */
86 /*HIPR1(30): (D=2$\times$HIPR1(31)=57.0 mb) the inclusive cross */
87 /* section $\sigma_{soft}$ for soft interactions. The default */
88 /* value $\sigma_{soft}=2\sigma_0$ is used to ensure the */
89 /* geometrical scaling of $pp$ interaction cross sections */
90 /* at low energies. */
91 /*HIPR1(31): (D=28.5 mb) the cross section $\sigma_0$ which */
92 /* characterizes the geometrical size of a nucleon */
93 /* ($\pi b_0^2=\sigma_0$, see Eq.~\ref{eq:over2}). */
94 /* The default value is only for high-energy */
95 /* limit ($\sqrt{s}>200$ GeV). At lower energies, a slight */
96 /* decrease which depends on energy is parametrized in the */
97 /* program. The default values of the two parameters */
98 /* HIPR1(30), HIPR1(31) are only for $NN$ type interactions. */
99 /* For other kinds of projectile or target hadrons, users */
100 /* should change these values so that correct inelastic */
101 /* and total cross sections (HINT1(12), HINT1(13)) are */
102 /* obtained by the program. */
103 /*HIPR1(32): (D=3.90) parameter $\mu_0$ in Eq.~\ref{eq:over2} for */
104 /* the scaled eikonal function. */
105 /*HIPR1(33): fractional cross section of single-diffractive */
106 /* interaction as parametrized in Ref.~\cite{goulianos}. */
107 /*HIPR1(34): maximum radial coordinate for projectile nucleons */
108 /* to be given by the initialization program HIJSET. */
109 /*HIPR1(35): maximum radial coordinate for target nucleons */
110 /* to be given by the initialization program HIJSET. */
111 /*HIPR1(36)-HIPR1(39): not used. */
112 /*HIPR1(40): (D=3.141592654) value of $\pi$. */
113 /*HIPR1(41)--HIPR1(42): not used. */
114 /*HIPR1(43): (D=0.01) fractional energy error relative to the */
115 /* colliding energy permitted per nucleon-nucleon collision. */
116 /*HIPR1(44), HIPR1(45), HIPR1(46): (D=1.5, 0.1 GeV, 0.25) parameters */
117 /* $\alpha$, $c$ and $\beta$ in the valence quark */
118 /* distributions for soft string excitation, */
119 /* $(1-x)^{\alpha}/(x^2+c^2/s)^{\beta}$ for baryons, */
120 /* $1/{(x^2+c^2/s)[(1-x)^2+c^2/s)]}^{\beta}$ for mesons. */
121 /*HIPR1(47), HIPR1(48): (D=0.0, 0.5) parameters $\alpha$ and $\beta$ */
122 /* in valence quark distribution, */
123 /* $(1-x)^{\alpha}/(x^2+c^2/s)^{\beta}$, for the */
124 /* disassociated excitation in a single diffractive collision. */
125 /*HIPR1(49)--HIPR1(100): not used. */
126 /*IHPR2(1): (D=1) switch for dipole-approximated QCD radiation */
127 /* of the string system in soft interactions. */
128 /*IHPR2(2): (D=3) option for initial and final state radiation in */
129 /* the hard scattering. */
130 /* =0: both initial and final radiation are off. */
131 /* =1: initial radiation on and final radiation off. */
132 /* =2: initial radiation off and final radiation on. */
133 /* =3: both initial and final radiation are on. */
134 /*IHPR2(3): (D=0) switch for triggered hard scattering with specified */
135 /* $P_T\geq$HIPR1(10). */
136 /* =0: no triggered jet production. */
137 /* =1: ordinary hard processes. */
138 /* =2: only direct photon production. */
139 /*IHPR2(4): (D=1) switch for jet quenching in the excited */
140 /* nuclear matter. */
141 /*IHPR2(5): (D=1) switch for the $P_T$ kick due to soft interactions. */
142 /*IHPR2(6): (D=1) switch for the nuclear effect on the parton */
143 /* distribution function such as shadowing. */
144 /*IHPR2(7): (D=1) selection of Duke-Owens set (1 or 2) of parametrization */
145 /* of nucleon structure functions. */
146 /*IHPR2(8): (D=10) maximum number of hard scatterings per */
147 /* nucleon-nucleon interaction. When IHPR2(8)=0, jet */
148 /* production will be turned off. When IHPR2(8)$<0$, the */
149 /* number of jet production will be fixed at its absolute */
150 /* value for each NN collision. */
151 /*IHPR2(9): (D=0) switch to guarantee at least one pair of minijets */
152 /* production per event ($pp$, $pA$ or $AB$). */
153 /*IHPR2(10): (D=0) option to print warning messages about errors that */
154 /* might happen. When a fatal error happens the current event */
155 /* will be abandoned and a new one is generated. */
156 /*IHPR2(11): (D=1) choice of baryon production model. */
157 /* =0: no baryon-antibaryon pair production, initial */
158 /* diquark treated as a unit. */
159 /* =1: diquark-antidiquark pair production allowed, */
160 /* initial diquark treated as a unit. */
161 /* =2: diquark-antidiquark pair production allowed, */
162 /* with the possibility for diquark to split */
163 /* according to the ``popcorn'' scheme (see the */
164 /* documentation of JETSET 7.2). */
165 /*IHPR2(12): (D=1) option to turn off the automatic decay of the */
166 /* following particles: */
167 /* $\pi^0$, $K^0_S$, $D^{\pm}$, $\Lambda$, $\Sigma^{\pm}$. */
168 /*IHPR2(13): (D=1) option to turn on single diffractive reactions. */
169 /*IHPR2(14): (D=1) option to turn on elastic scattering. */
170 /*IHPR2(15)--IHPR2(18): not used. */
171 /*IHPR2(19): (D=1) option to turn on initial state soft interaction. */
172 /*IHPR2(20): (D=1) switch for the final fragmentation. */
173 /*IHPR2(21): (D=0) option to keep the information of all particles */
174 /* including those which have decayed and the decay history */
175 /* in the common block HIMAIN2. The line number of the parent */
176 /* particle is KATT(I,3). The status of a partcile, */
177 /* whether it is a finally produced particle (KATT(I,4)=1) */
178 /* or a decayed particle (KATT(I,4)=11) is also kept. */
179 /*IHPR2(22)-IHPR2(50): not used. */
180 /*HINT1(1): (GeV) colliding energy in the c.m. frame of nucleon-nucleon */
182 /*HINT1(2): Lorentz transformation variable $\beta$ from laboratory */
183 /* to c.m. frame of nucleon nucleon collisions. */
184 /*HINT1(3): rapidity $y_{cm}$ of the c.m. frame */
185 /* $\beta=\tanh y_{cm}$. */
186 /*HINT1(4): rapidity of projectile nucleons (hadron) $y_{proj}$. */
187 /*HINT1(5): rapidity of target nucleons (hadron) $y_{targ}$. */
188 /*HINT1(6): (GeV) energy of the projectile nucleons (hadron) in the */
190 /*HINT1(7): (GeV) energy of the target nucleons (hadron) in the */
192 /*HINT1(8): (GeV) the rest mass of projectile particles. */
193 /*HINT1(9): (GeV) the rest mass of target particles. */
194 /*HINT1(10): (mb) the averaged cross section for jet production */
195 /* per nucleon-nucleon collisions, */
196 /* $\int d^2b\{1-\exp[-\sigma_{jet}T_N(b)]\}$. */
197 /*HINT1(11): (mb) the averaged inclusive cross section $\sigma_{jet}$ */
198 /* for jet production per nucleon-nucleon collisions. */
199 /*HINT1(12): (mb) the averaged inelastic cross section of */
200 /* nucleon-nucleon collisions. */
201 /*HINT1(13): (mb) the averaged total cross section of nucleon-nucleon */
203 /*HINT1(14): (mb) the jet production cross section without nuclear */
204 /* shadowing effect $\sigma_{jet}^0$ (see Eq.~\ref{eq:sjetab}). */
205 /*HINT1(15): (mb) the cross section $\sigma_{jet}^A$ to account for */
206 /* the projectile shadowing correction term in the jet cross */
207 /* section (see Eq.~\ref{eq:sjetab}). */
208 /*HINT1(16): (mb) the cross section $\sigma_{jet}^B$ to account for */
209 /* the target shadowing correction term in the jet cross */
210 /* section (see Eq.~\ref{eq:sjetab}). */
211 /*HINT1(17): (mb) the cross section $\sigma_{jet}^{AB}$ to account */
212 /* for the cross term of shadowing correction in the jet */
214 /*HINT1(18): (mb) the effective cross section */
215 /* $\sigma_{jet}^{eff}(r_A,r_B)$ for jet production */
216 /* of the latest nucleon-nucleon collision which depends */
217 /* on the transverse coordinates of the colliding */
219 /*HINT1(19): (fm) the (absolute value of) impact parameter of the */
221 /*HINT1(20): (radians) the azimuthal angle $\phi$ of the impact */
222 /* parameter vector in the transverse plane of the latest */
224 /*HINT1(21)--HINT1(25): the four momentum and mass ($p_x,p_y,p_z,E,M$) */
225 /* (GeV/$c$, GeV, GeV/$c^2$) of the first scattered parton */
226 /* in the triggered hard scattering. This is before the final */
227 /* state radiation but after the initial state radiation. */
228 /*HINT1(26)--HINT1(30): not used. */
229 /*HINT1(31)--HINT1(35): the four momentum and mass ($p_x,p_y,p_z,E,M$) */
230 /* (GeV/$c$, GeV, GeV/$c^2$) of the second scattered parton */
231 /* in the triggered hard scattering. This is before the final */
232 /* state radiation but after the initial state radiation. */
233 /*HINT1(46)--HINT1(40): not used. */
234 /*HINT1(41)--HINT1(45): the four momentum and mass ($p_x,p_y,p_z,E,M$) */
235 /* (GeV/$c$, GeV, GeV/$c^2$) of the first scattered parton */
236 /* in the latest hard scattering of the latest event. */
237 /*HINT1(46): $P_T$ (GeV/$c$) of the first scattered parton in the */
238 /* latest hard scattering of the latest event. */
239 /*HINT1(47)--HINT1(50): not used. */
240 /*HINT1(51)--HINT1(55): the four momentum and mass ($p_x,p_y,p_z,E,M$) */
241 /* (GeV/$c$, GeV, GeV/$c^2$) of the second scattered parton */
242 /* in the latest hard scattering of the latest event. */
243 /*HINT1(56): $P_T$ (GeV/$c$) of the second scattered parton in the */
244 /* latest hard scattering of the latest event. */
245 /*HINT1(57)--HINT1(58): not used. */
246 /*HINT1(59): (mb) the averaged cross section of the */
247 /* triggered jet production (with $P_T$ specified by HIPR1(10) */
248 /* and with switch by IHPR2(3)) per nucleon-nucleon */
250 /* $\int d^2b\{1-\exp[-\sigma_{jet}^{trig}T_N(b)]\}$ */
251 /*HINT1(60): (mb) the averaged inclusive cross section of the */
252 /* triggered jet production $\sigma_{jet}^{trig}$ */
253 /* (with $P_T$ specified by */
254 /* HIPR1(10) and with switch by IHPR2(3)) per */
255 /* nucleon-nucleon collision. */
256 /*HINT1(61): (mb) the triggered jet production cross section without */
257 /* nuclear shadowing effect (similar to HINT1(14)). */
258 /*HINT1(62): (mb) the cross section to account for the projectile */
259 /* shadowing correction term in the triggered jet cross */
260 /* section (similar to HINT1(15)). */
261 /*HINT1(63): (mb) the cross section to account for the target */
262 /* shadowing correction term in the triggered jet cross */
263 /* section (similar to HINT1(16)). */
264 /*HINT1(64): (mb) the cross section to account for the cross */
265 /* term of shadowing correction in the triggered jet */
266 /* cross section (similar to HINT1(17). */
267 /*HINT1(65): (mb) the inclusive cross section for latest triggered */
268 /* jet production which depends on the transverse coordinates */
269 /* of the colliding nucleons (similar to HINT1(18)). */
270 /*HINT1(67)--HINT1(71): not used. */
271 /*HINT1(72)--HINT1(75): three parameters for the Wood-Saxon */
272 /* projectile nuclear distribution and the normalization */
273 /* read from a table inside the program, */
274 /* $\rho(r)=C[1+W(r/R_A)^2]/\{1+\exp[(r-R_A)/D]\}$, */
275 /* $R_A$=HINT1(72), $D$=HINT1(73), $W$=HINT1(74), $C$=HINT1(75). */
276 /*HINT1(76)--HINT1(79): three parameters for the Wood-Saxon */
277 /* projectile nuclear distribution and the normalization */
278 /* read from a table inside the program, */
279 /* $\rho(r)=C[1+W(r/R_A)^2]/\{1+\exp[(r-R_A)/D]\}$, */
280 /* $R_A$=HINT1(76), $D$=HINT1(77), $W$=HINT1(78), $C$=HINT1(79). */
281 /*HINT1(80)--HINT1(100): the probability of $j=0-20$ number of hard */
282 /* scatterings per nucleon-nucleon collisions. */
283 /*IHNT2(1): the mass number of the projectile nucleus (1 for a hadron). */
284 /*IHNT2(2): the charge number of the projectile nucleus. If the */
285 /* projectile is a hadron, it gives the charge of the hadron. */
286 /*IHNT2(3): the mass number of the target nucleus (1 for a hadron). */
287 /*IHNT2(4): the charge number of the target nucleus. If the target */
288 /* is a hadron, it gives the charge of the hadron. */
289 /*IHNT2(5): the flavor code of the projectile hadron (0 for nucleus). */
290 /*IHNT2(6): the flavor code of the target hadron (0 for nucleus). */
291 /*IHNT2(7)--IHNT2(8): not used. */
292 /*IHNT2(9): the flavor code of the first scattered parton in the */
293 /* triggered hard scattering. */
294 /*IHNT2(10): the flavor code of the second scattered parton in the */
295 /* triggered hard scattering. */
296 /*IHNT2(11): the sequence number of the projectile nucleon in the */
297 /* latest nucleon-nucleon interaction of the latest event. */
298 /*IHNT2(12): the sequence number of the target nucleon in the latest */
299 /* nucleon-nucleon interaction of the latest event. */
300 /*IHNT2(13): status of the latest soft string excitation. */
301 /* =1: double diffractive. */
302 /* =2: single diffractive. */
303 /* =3: non-single diffractive. */
304 /*IHNT2(14): the flavor code of the first scattered parton in the */
305 /* latest hard scattering of the latest event. */
306 /*IHNT2(15): the flavor code of the second scattered parton in the */
307 /* latest hard scattering of the latest event. */
308 /*IHNT2(16)--IHNT2(50): not used. */
309 /*========================================================================*/
311 /*========================================================================*/
312 /* COMMON/HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11,BB */
313 /*------------------------------------------------------------------------*/
327 #define HIMAIN1 COMMON_BLOCK(HIMAIN1,himain1)
328 COMMON_BLOCK_DEF(Himain1Common,HIMAIN1);
329 /*************************************************************************/
330 /* D E S C R I P T I O N : */
331 /*-----------------------------------------------------------------------*/
332 /*COMMON/HIMAIN1/NATT, EATT, JATT, NT, NP, N0, N01, N10, N11 */
333 /*Purpose: to give the overall information of the generated event. */
334 /*NATT: total number of produced stable and undecayed particles of */
335 /* the current event. */
336 /*EATT: the total energy of the produced particles in c.m. frame */
337 /* of the collision to check energy conservation. */
338 /*JATT: the total number of hard scatterings in the current event. */
339 /*NP, NT: the number of participant projectile and target nucleons */
340 /* in the current event. */
341 /*N0, N01, N10, N11: number of $N$-$N$, $N$-$N_{wounded}$, */
342 /* $N_{wounded}$-$N$, and */
343 /* $N_{wounded}$-$N_{wounded}$ collisions in */
344 /* the current event ($N$, $N_{wounded}$ stand */
345 /* for nucleon and wounded nucleon respectively). */
346 /*=======================================================================*/
349 /*========================================================*/
350 /* COMMON/HIMAIN2/KATT(200000,4),PATT(200000,4),VATT(200000,4)*/
351 /*--------------------------------------------------------*/
353 Int_t katt[4][200000];
354 Float_t patt[4][200000];
355 Float_t vatt[4][200000];
358 #define HIMAIN2 COMMON_BLOCK(HIMAIN2,himain2)
359 COMMON_BLOCK_DEF(Himain2Common,HIMAIN2);
360 /*************************************************************************/
361 /* D E S C R I P T I O N : */
362 /*-----------------------------------------------------------------------*/
363 /*Purpose: to give information of produced stable and undecayed */
364 /* particles. Parent particles which decayed are not included */
366 /*KATT(I, 1): (I=1,$\cdots$,NATT) flavor codes (see appendix) of */
367 /* the produced particles. */
368 /*KATT(I, 2): (I=1,$\cdots$,NATT) status codes to identify the */
369 /* sources from which the particles come. */
370 /* =0: projectile nucleon (or hadron) which has */
371 /* not interacted at all. */
372 /* =1: projectile nucleon (or hadron) which */
373 /* only suffers an elastic collision. */
374 /* =2: from a diffractive projectile nucleon (or hadron) */
375 /* in a single diffractive interaction. */
376 /* =3: from the fragmentation of a projectile string */
377 /* system (including gluon jets). */
378 /* =10 target nucleon (or hadron) which has not */
379 /* interacted at all. */
380 /* =11: target nucleon (or hadron) which only */
381 /* suffers an elastic collision. */
382 /* =12: from a diffractive target nucleon (or hadron) */
383 /* in a single diffractive interaction. */
384 /* =13: from the fragmentation of a target string */
385 /* system (including gluon jets). */
386 /* =20: from scattered partons which form string */
387 /* systems themselves. */
388 /* =40: from direct production in the hard processes */
389 /* ( currently, only direct photons are included). */
390 /*KATT(I,3): (I=1,$\cdots$,NATT) line number of the parent particle. */
391 /* For finally produced or directly produced (not from */
392 /* the decay of another particle) particles, it is set */
393 /* to 0 (The option to keep the information of all */
394 /* particles including the decayed ones is IHPR2(21)=1). */
395 /*KATT(I,4): (I=1,$\cdots$,NATT) status number of the particle. */
396 /* =1: finally or directly produced particles. */
397 /* =11: particles which has already decayed. */
398 /*PATT(I, 1-4): (I=1,$\cdots$,NATT) four-momentum ($p_x,p_y,p_z,E$) */
399 /* (GeV/$c$, GeV) of the produced particles. */
401 /*=======================================================================*/
403 /*=======================================================================*/
404 /* COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),PJPY(300,500) */
405 /* & ,PJPZ(300,500),PJPE(300,500),PJPM(300,500) */
406 /* & ,NTJ(300),KFTJ(300,500),PJTX(300,500),PJTY(300,500) */
407 /* & ,PJTZ(300,500),PJTE(300,500),PJTM(300,500) */
408 /*-----------------------------------------------------------------------*/
411 Int_t kfpj[500][300];
412 Float_t pjpx[500][300];
413 Float_t pjpy[500][300];
414 Float_t pjpz[500][300];
415 Float_t pjpe[500][300];
416 Float_t pjpm[500][300];
418 Int_t kftj[500][300];
419 Float_t pjtx[500][300];
420 Float_t pjty[500][300];
421 Float_t pjtz[500][300];
422 Float_t pjte[500][300];
423 Float_t pjtm[500][300];
426 #define HIJJET1 COMMON_BLOCK(HIJJET1,hijjet1)
427 COMMON_BLOCK_DEF(Hijjet1Common,HIJJET1);
428 /*************************************************************************/
429 /* D E S C R I P T I O N : */
430 /*-----------------------------------------------------------------------*/
431 /*Purpose: contains information about produced partons which are */
432 /* connected with the valence quarks and diquarks of */
433 /* projectile or target nucleons (or hadron) to form */
434 /* string systems for fragmentation. The momentum and */
435 /* energy of all produced partons are calculated in */
436 /* the c.m. frame of the collision. IAP, IAT are the */
437 /* numbers of nucleons in projectile and target nucleus */
438 /* respectively (IAP, IAT=1 for hadron projectile or target). */
439 /*NPJ(I): (I=1,$\cdots$,IAP) number of partons associated with projectile*/
441 /*KFPJ(I, J): (I=1,$\cdots$,IAP, J=1,$\cdots$,NPJ(I)) parton */
442 /* flavor code of the */
443 /* parton J associated with projectile nucleon I. */
444 /*PJPX(I, J), PJPY(I, J), PJPZ(I, J), PJPE(I, J), PJPM(I, J): the four */
445 /* momentum and mass ($p_x,p_y,p_z,E,M$) */
446 /* (GeV/$c$, GeV, GeV/$c^2$) of parton J associated with */
447 /* the projectile nucleon I. */
448 /*NTJ(I): (I=1,$\cdots$,IAT) number of partons associated with */
449 /* target nucleon I. */
450 /*KFTJ(I, J): (I=1,$\cdots$,IAT, J=1,$\cdots$,NTJ(I)): parton */
451 /* flavor code of the parton J associated with */
452 /* target nucleon I. */
453 /*PJTX(I, J), PJTY(I, J), PJTZ(I, J), PJTE(I, J), PJTM(I, J): the four */
454 /* momentum and mass ($p_x,p_y,p_z,E,M$) */
455 /* (GeV/$c$, GeV, GeV/$c^2$) of parton J associated with */
456 /* target nucleon I. */
458 /*=======================================================================*/
460 /*=======================================================================*/
461 /* COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100), */
462 /* & K2SG(900,100),PXSG(900,100),PYSG(900,100), */
463 /* & PZSG(900,100),PESG(900,100),PMSG(900,100) */
464 /*-----------------------------------------------------------------------*/
469 Int_t k1sg[100][900];
470 Int_t k2sg[100][900];
471 Float_t pxsg[100][900];
472 Float_t pysg[100][900];
473 Float_t pzsg[100][900];
474 Float_t pesg[100][900];
475 Float_t pmsg[100][900];
478 #define HIJJET2 COMMON_BLOCK(HIJJET2,hijjet2)
479 COMMON_BLOCK_DEF(Hijjet2Common,HIJJET2);
480 /*************************************************************************/
481 /* D E S C R I P T I O N : */
482 /*-----------------------------------------------------------------------*/
483 /*Purpose: contains information about the produced partons which */
484 /* will form string systems themselves without being */
485 /* connected to valence quarks and diquarks. */
486 /*NSG: the total number of such string systems. */
487 /*NJSG(I): (I=1,$\cdots$,NSG) number of partons in the string system I. */
488 /*IASG(I, 1), IASG(I, 2): to specify which projectile and target */
489 /* nucleons produce string system I. */
490 /*IASG(I, 3): to indicate whether the jets will be quenched (0) */
491 /* or will not be quenched (1). */
492 /*K1SG(I, J): (J=1,$\cdots$,NJSG(I)) color flow information of parton J */
493 /* in string system I (see JETSET 7.2 for detailed */
495 /*K2SG(I, J): (J=1,$\cdots$,NJSG(I)) flavor code of parton J in string */
497 /*PXSG(I, J), PYSG(I, J), PZSG(I, J), PESG(I, J), PMSG(I, J): four */
498 /* momentum and mass ($p_x,p_y,p_z,E,M$) */
499 /* ( GeV/$c$, GeV, GeV/$c^2$) of parton J in string system I. */
500 /*=======================================================================*/
502 /*=======================================================================*/
503 /* COMMON/HISTRNG/NFP(300,15),PP(300,15),NFT(300,15),PT(300,15) */
504 /*-----------------------------------------------------------------------*/
512 #define HISTRNG COMMON_BLOCK(HISTRNG,histrng)
513 COMMON_BLOCK_DEF(HistrngCommon,HISTRNG);
514 /*************************************************************************/
515 /* D E S C R I P T I O N : */
516 /*-----------------------------------------------------------------------*/
517 /*Purpose: contains information about the projectile and */
518 /* target nucleons (hadron) and the corresponding constituent */
519 /* quarks, diquarks. IAP, IAT are the numbers of nucleons in */
520 /* projectile and target nucleus respectively (IAP, IAT=1 */
521 /* for hadron projectile or target). */
522 /*NFP(I, 1): (I=1,$\cdots$,IAP) flavor code of the valence quark in */
523 /* projectile nucleon (hadron) I. */
524 /*NFP(I, 2): flavor code of diquark in projectile nucleon (anti-quark */
525 /* in projectile meson) I. */
526 /*NFP(I, 3): present flavor code of the projectile nucleon (hadron) I */
527 /* ( a nucleon or meson can be excited to its vector resonance). */
528 /*NFP(I, 4): original flavor code of projectile nucleon (hadron) I. */
529 /*NFP(I, 5): collision status of projectile nucleon (hadron) I. */
530 /* =0: suffered no collision. */
531 /* =1: suffered an elastic collision. */
532 /* =2: being the diffractive one in a single-diffractive */
534 /* =3: became an excited string after an inelastic */
536 /*NFP(I, 6): the total number of hard scatterings associated with */
537 /* projectile nucleon (hadron) I. If NFP(I,6)$<0$, it can not */
538 /* produce jets any more due to energy conservation. */
539 /*NFP(I, 10): to indicate whether the valence quarks or diquarks */
540 /* (anti-quarks) in projectile nucleon (hadron) I */
541 /* suffered a hard scattering, */
542 /* =0: has not suffered a hard scattering. */
543 /* =1: suffered one or more hard scatterings in */
544 /* current binary nucleon-nucleon collision. */
545 /* =-1: suffered one or more hard scatterings in */
546 /* previous binary nucleon-nucleon collisions. */
547 /*NFP(I, 11): total number of interactions projectile nucleon (hadron) */
548 /* I has suffered so far. */
549 /*PP(I, 1), PP(I, 2), PP(I, 3), PP(I, 4), PP(I, 5): four momentum and */
550 /* the invariant mass ($p_x,p_y,p_z,E,M$) */
551 /* (GeV/$c$, GeV, GeV/$c^2$) of projectile nucleon (hadron) I. */
552 /*PP(I, 6), PP(I, 7): transverse momentum ($p_x,p_y$) (GeV/$c$) of the */
553 /* valence quark in projectile nucleon (hadron) I. */
554 /*PP(I, 8), PP(I, 9): transverse momentum ($p_x,p_y$) (GeV/$c$) of the */
555 /* diquark (anti-quark) in projectile nucleon (hadron) I. */
556 /*PP(I, 10), PP(I, 11), PP(I, 12): three momentum ($p_x,p_y,p_z$) */
557 /* (GeV/$c$) transferred to the quark or diquark (anti-quark) */
558 /* in projectile nucleon (hadron) I from the last hard */
560 /*PP(I, 14): mass (GeV/$c^2$) of the quark in projectile nucleon */
562 /*PP(I, 15): mass of the diquark (anti-quark) in projectile */
563 /* nucleon (hadron) I. */
564 /*NFT(I, 1--15), PT(I,1--15): give the same */
565 /* information for the target nucleons (hadron) and the */
566 /* corresponding quarks and diquarks (anti-quarks) as for */
567 /* the projectile nucleons. */
569 /*=======================================================================*/
571 // COMMON/LUDAT1_HIJING/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
572 // SAVE /LUDAT1_HIJING/
578 } Ludat1_HijingCommon;
580 #define LUDAT1_HIJING COMMON_BLOCK(LUDAT1_HIJING,ludat1_hijing)
581 COMMON_BLOCK_DEF(Ludat1_HijingCommon,LUDAT1_HIJING);