Do not include from subdirectories
[u/mrichter/AliRoot.git] / THijing / Hcommon.h
CommitLineData
7dffa423 1#ifndef ROOT_HCommon
2#define ROOT_HCommon
3
4#ifndef __CFORTRAN_LOADED
5//*KEEP,cfortran.
6#include "cfortran.h"
7//*KEND.
8#endif
9
10extern "C" {
11
12/*=========================================================*/
13/* COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)*/
14/*---------------------------------------------------------*/
15typedef struct {
16 Float_t hipr1[100];
17 Int_t ihpr2[50];
18 Float_t hint1[100];
19 Int_t ihnt2[50];
20} HiparntCommon;
21
22#define HIPARNT COMMON_BLOCK(HIPARNT,hiparnt)
23COMMON_BLOCK_DEF(HiparntCommon,HIPARNT);
24
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 */
181/* collisions. */
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 */
189/* given frame. */
190/*HINT1(7): (GeV) energy of the target nucleons (hadron) in the */
191/* given frame. */
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 */
202/* collisions. */
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 */
213/* cross section. */
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 */
218/* nucleons. */
219/*HINT1(19): (fm) the (absolute value of) impact parameter of the */
220/* latest event. */
221/*HINT1(20): (radians) the azimuthal angle $\phi$ of the impact */
222/* parameter vector in the transverse plane of the latest */
223/* event. */
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 */
249/* collision, */
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/*========================================================================*/
310
311/*========================================================================*/
312/* COMMON/HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11,BB */
313/*------------------------------------------------------------------------*/
314typedef struct {
315 Int_t natt;
316 Float_t eatt;
317 Int_t jatt;
318 Int_t nt;
319 Int_t np;
320 Int_t n0;
321 Int_t n01;
322 Int_t n10;
323 Int_t n11;
369eeebb 324 Float_t bb;
325 Int_t npart;
7dffa423 326 } Himain1Common;
327
328#define HIMAIN1 COMMON_BLOCK(HIMAIN1,himain1)
329COMMON_BLOCK_DEF(Himain1Common,HIMAIN1);
330/*************************************************************************/
331/* D E S C R I P T I O N : */
332/*-----------------------------------------------------------------------*/
333/*COMMON/HIMAIN1/NATT, EATT, JATT, NT, NP, N0, N01, N10, N11 */
334/*Purpose: to give the overall information of the generated event. */
335/*NATT: total number of produced stable and undecayed particles of */
336/* the current event. */
337/*EATT: the total energy of the produced particles in c.m. frame */
338/* of the collision to check energy conservation. */
339/*JATT: the total number of hard scatterings in the current event. */
340/*NP, NT: the number of participant projectile and target nucleons */
341/* in the current event. */
342/*N0, N01, N10, N11: number of $N$-$N$, $N$-$N_{wounded}$, */
343/* $N_{wounded}$-$N$, and */
344/* $N_{wounded}$-$N_{wounded}$ collisions in */
345/* the current event ($N$, $N_{wounded}$ stand */
346/* for nucleon and wounded nucleon respectively). */
347/*=======================================================================*/
348
349
350/*========================================================*/
34c6d0b9 351/* COMMON/HIMAIN2/KATT(200000,4),PATT(200000,4),VATT(200000,4)*/
7dffa423 352/*--------------------------------------------------------*/
353typedef struct {
c5af11ff 354 Int_t katt[4][200000];
355 Float_t patt[4][200000];
c9080d26 356 Float_t vatt[4][200000];
7dffa423 357} Himain2Common;
358
359#define HIMAIN2 COMMON_BLOCK(HIMAIN2,himain2)
360COMMON_BLOCK_DEF(Himain2Common,HIMAIN2);
361/*************************************************************************/
362/* D E S C R I P T I O N : */
363/*-----------------------------------------------------------------------*/
364/*Purpose: to give information of produced stable and undecayed */
365/* particles. Parent particles which decayed are not included */
366/* here. */
367/*KATT(I, 1): (I=1,$\cdots$,NATT) flavor codes (see appendix) of */
368/* the produced particles. */
369/*KATT(I, 2): (I=1,$\cdots$,NATT) status codes to identify the */
370/* sources from which the particles come. */
371/* =0: projectile nucleon (or hadron) which has */
372/* not interacted at all. */
373/* =1: projectile nucleon (or hadron) which */
374/* only suffers an elastic collision. */
375/* =2: from a diffractive projectile nucleon (or hadron) */
376/* in a single diffractive interaction. */
377/* =3: from the fragmentation of a projectile string */
378/* system (including gluon jets). */
379/* =10 target nucleon (or hadron) which has not */
380/* interacted at all. */
381/* =11: target nucleon (or hadron) which only */
382/* suffers an elastic collision. */
383/* =12: from a diffractive target nucleon (or hadron) */
384/* in a single diffractive interaction. */
385/* =13: from the fragmentation of a target string */
386/* system (including gluon jets). */
387/* =20: from scattered partons which form string */
388/* systems themselves. */
389/* =40: from direct production in the hard processes */
390/* ( currently, only direct photons are included). */
391/*KATT(I,3): (I=1,$\cdots$,NATT) line number of the parent particle. */
392/* For finally produced or directly produced (not from */
393/* the decay of another particle) particles, it is set */
394/* to 0 (The option to keep the information of all */
395/* particles including the decayed ones is IHPR2(21)=1). */
396/*KATT(I,4): (I=1,$\cdots$,NATT) status number of the particle. */
397/* =1: finally or directly produced particles. */
398/* =11: particles which has already decayed. */
399/*PATT(I, 1-4): (I=1,$\cdots$,NATT) four-momentum ($p_x,p_y,p_z,E$) */
400/* (GeV/$c$, GeV) of the produced particles. */
401/* */
402/*=======================================================================*/
403
404/*=======================================================================*/
405/* COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),PJPY(300,500) */
406/* & ,PJPZ(300,500),PJPE(300,500),PJPM(300,500) */
407/* & ,NTJ(300),KFTJ(300,500),PJTX(300,500),PJTY(300,500) */
408/* & ,PJTZ(300,500),PJTE(300,500),PJTM(300,500) */
409/*-----------------------------------------------------------------------*/
410typedef struct {
411 Int_t npj[300];
412 Int_t kfpj[500][300];
413 Float_t pjpx[500][300];
414 Float_t pjpy[500][300];
415 Float_t pjpz[500][300];
416 Float_t pjpe[500][300];
417 Float_t pjpm[500][300];
418 Int_t ntj[300];
419 Int_t kftj[500][300];
420 Float_t pjtx[500][300];
421 Float_t pjty[500][300];
422 Float_t pjtz[500][300];
423 Float_t pjte[500][300];
424 Float_t pjtm[500][300];
425} Hijjet1Common;
426
427#define HIJJET1 COMMON_BLOCK(HIJJET1,hijjet1)
428COMMON_BLOCK_DEF(Hijjet1Common,HIJJET1);
429/*************************************************************************/
430/* D E S C R I P T I O N : */
431/*-----------------------------------------------------------------------*/
432/*Purpose: contains information about produced partons which are */
433/* connected with the valence quarks and diquarks of */
434/* projectile or target nucleons (or hadron) to form */
435/* string systems for fragmentation. The momentum and */
436/* energy of all produced partons are calculated in */
437/* the c.m. frame of the collision. IAP, IAT are the */
438/* numbers of nucleons in projectile and target nucleus */
439/* respectively (IAP, IAT=1 for hadron projectile or target). */
440/*NPJ(I): (I=1,$\cdots$,IAP) number of partons associated with projectile*/
441/* nucleon I. */
442/*KFPJ(I, J): (I=1,$\cdots$,IAP, J=1,$\cdots$,NPJ(I)) parton */
443/* flavor code of the */
444/* parton J associated with projectile nucleon I. */
445/*PJPX(I, J), PJPY(I, J), PJPZ(I, J), PJPE(I, J), PJPM(I, J): the four */
446/* momentum and mass ($p_x,p_y,p_z,E,M$) */
447/* (GeV/$c$, GeV, GeV/$c^2$) of parton J associated with */
448/* the projectile nucleon I. */
449/*NTJ(I): (I=1,$\cdots$,IAT) number of partons associated with */
450/* target nucleon I. */
451/*KFTJ(I, J): (I=1,$\cdots$,IAT, J=1,$\cdots$,NTJ(I)): parton */
452/* flavor code of the parton J associated with */
453/* target nucleon I. */
454/*PJTX(I, J), PJTY(I, J), PJTZ(I, J), PJTE(I, J), PJTM(I, J): the four */
455/* momentum and mass ($p_x,p_y,p_z,E,M$) */
456/* (GeV/$c$, GeV, GeV/$c^2$) of parton J associated with */
457/* target nucleon I. */
458/* */
459/*=======================================================================*/
460
461/*=======================================================================*/
462/* COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100), */
463/* & K2SG(900,100),PXSG(900,100),PYSG(900,100), */
464/* & PZSG(900,100),PESG(900,100),PMSG(900,100) */
465/*-----------------------------------------------------------------------*/
466typedef struct {
467 Int_t nsg;
468 Int_t njsg[900];
469 Int_t iasg[3][900];
470 Int_t k1sg[100][900];
471 Int_t k2sg[100][900];
472 Float_t pxsg[100][900];
473 Float_t pysg[100][900];
474 Float_t pzsg[100][900];
475 Float_t pesg[100][900];
476 Float_t pmsg[100][900];
477} Hijjet2Common;
478
479#define HIJJET2 COMMON_BLOCK(HIJJET2,hijjet2)
480COMMON_BLOCK_DEF(Hijjet2Common,HIJJET2);
481/*************************************************************************/
482/* D E S C R I P T I O N : */
483/*-----------------------------------------------------------------------*/
484/*Purpose: contains information about the produced partons which */
485/* will form string systems themselves without being */
486/* connected to valence quarks and diquarks. */
487/*NSG: the total number of such string systems. */
488/*NJSG(I): (I=1,$\cdots$,NSG) number of partons in the string system I. */
489/*IASG(I, 1), IASG(I, 2): to specify which projectile and target */
490/* nucleons produce string system I. */
491/*IASG(I, 3): to indicate whether the jets will be quenched (0) */
492/* or will not be quenched (1). */
493/*K1SG(I, J): (J=1,$\cdots$,NJSG(I)) color flow information of parton J */
494/* in string system I (see JETSET 7.2 for detailed */
495/* explanation). */
496/*K2SG(I, J): (J=1,$\cdots$,NJSG(I)) flavor code of parton J in string */
497/* system I. */
498/*PXSG(I, J), PYSG(I, J), PZSG(I, J), PESG(I, J), PMSG(I, J): four */
499/* momentum and mass ($p_x,p_y,p_z,E,M$) */
500/* ( GeV/$c$, GeV, GeV/$c^2$) of parton J in string system I. */
501/*=======================================================================*/
502
503/*=======================================================================*/
504/* COMMON/HISTRNG/NFP(300,15),PP(300,15),NFT(300,15),PT(300,15) */
505/*-----------------------------------------------------------------------*/
506typedef struct {
507 Int_t nfp[15][300];
508 Float_t pp[15][300];
509 Int_t nft[15][300];
510 Float_t pt[15][300];
511} HistrngCommon;
512
513#define HISTRNG COMMON_BLOCK(HISTRNG,histrng)
514COMMON_BLOCK_DEF(HistrngCommon,HISTRNG);
65a17253 515
516
517/*************************************************************************/
518/* D E S C R I P T I O N : */
519/*-----------------------------------------------------------------------*/
520/* The following common block is added to record the number of elastic*/
521/* (NELT, NELP) and inelastic (NINT, NINP) participants */
522/* */
523/*=======================================================================*/
524/* COMMON/HIJGLBR/NELT,NINT,NELP,NINP */
525/* SAVE /HIJGLBR/ */
526/*=======================================================================*/
527typedef struct {
528 Int_t nelt;
529 Int_t nint;
530 Int_t nelp;
531 Int_t ninp;
532 Int_t npspecp;
533 Int_t nnspecp;
534 Int_t npspect;
535 Int_t nnspect;
536} HijglbrCommon;
537
538#define HIJGLBR COMMON_BLOCK(HIJGLBR,hijglbr)
539COMMON_BLOCK_DEF(HijglbrCommon,HIJGLBR);
540
541
7dffa423 542/*************************************************************************/
543/* D E S C R I P T I O N : */
544/*-----------------------------------------------------------------------*/
545/*Purpose: contains information about the projectile and */
546/* target nucleons (hadron) and the corresponding constituent */
547/* quarks, diquarks. IAP, IAT are the numbers of nucleons in */
548/* projectile and target nucleus respectively (IAP, IAT=1 */
549/* for hadron projectile or target). */
550/*NFP(I, 1): (I=1,$\cdots$,IAP) flavor code of the valence quark in */
551/* projectile nucleon (hadron) I. */
552/*NFP(I, 2): flavor code of diquark in projectile nucleon (anti-quark */
553/* in projectile meson) I. */
554/*NFP(I, 3): present flavor code of the projectile nucleon (hadron) I */
555/* ( a nucleon or meson can be excited to its vector resonance). */
556/*NFP(I, 4): original flavor code of projectile nucleon (hadron) I. */
557/*NFP(I, 5): collision status of projectile nucleon (hadron) I. */
558/* =0: suffered no collision. */
559/* =1: suffered an elastic collision. */
560/* =2: being the diffractive one in a single-diffractive */
561/* collision. */
562/* =3: became an excited string after an inelastic */
563/* collision. */
564/*NFP(I, 6): the total number of hard scatterings associated with */
565/* projectile nucleon (hadron) I. If NFP(I,6)$<0$, it can not */
566/* produce jets any more due to energy conservation. */
567/*NFP(I, 10): to indicate whether the valence quarks or diquarks */
568/* (anti-quarks) in projectile nucleon (hadron) I */
569/* suffered a hard scattering, */
570/* =0: has not suffered a hard scattering. */
571/* =1: suffered one or more hard scatterings in */
572/* current binary nucleon-nucleon collision. */
573/* =-1: suffered one or more hard scatterings in */
574/* previous binary nucleon-nucleon collisions. */
575/*NFP(I, 11): total number of interactions projectile nucleon (hadron) */
576/* I has suffered so far. */
577/*PP(I, 1), PP(I, 2), PP(I, 3), PP(I, 4), PP(I, 5): four momentum and */
578/* the invariant mass ($p_x,p_y,p_z,E,M$) */
579/* (GeV/$c$, GeV, GeV/$c^2$) of projectile nucleon (hadron) I. */
580/*PP(I, 6), PP(I, 7): transverse momentum ($p_x,p_y$) (GeV/$c$) of the */
581/* valence quark in projectile nucleon (hadron) I. */
582/*PP(I, 8), PP(I, 9): transverse momentum ($p_x,p_y$) (GeV/$c$) of the */
583/* diquark (anti-quark) in projectile nucleon (hadron) I. */
584/*PP(I, 10), PP(I, 11), PP(I, 12): three momentum ($p_x,p_y,p_z$) */
585/* (GeV/$c$) transferred to the quark or diquark (anti-quark) */
586/* in projectile nucleon (hadron) I from the last hard */
587/* scattering. */
588/*PP(I, 14): mass (GeV/$c^2$) of the quark in projectile nucleon */
589/* (hadron) I. */
590/*PP(I, 15): mass of the diquark (anti-quark) in projectile */
591/* nucleon (hadron) I. */
592/*NFT(I, 1--15), PT(I,1--15): give the same */
593/* information for the target nucleons (hadron) and the */
594/* corresponding quarks and diquarks (anti-quarks) as for */
595/* the projectile nucleons. */
596/* */
597/*=======================================================================*/
c9080d26 598
599// COMMON/LUDAT1_HIJING/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
600// SAVE /LUDAT1_HIJING/
601 typedef struct {
602 Int_t mstu[200];
603 Float_t paru[200];
604 Int_t mstj[200];
605 Float_t parj[200];
606} Ludat1_HijingCommon;
607
608#define LUDAT1_HIJING COMMON_BLOCK(LUDAT1_HIJING,ludat1_hijing)
609COMMON_BLOCK_DEF(Ludat1_HijingCommon,LUDAT1_HIJING);
7dffa423 610}
611
0f07b55b 612//COMMON/LUDAT3_HIJING/MDCY(500, 3),MDME(2000, 2),BRAT(2000),KFDP(2000, 5)
613 typedef struct {
614 Int_t mdcy[3][500];
615 Int_t mdme[2][2000];
616 Float_t brat[2000];
617 Float_t kfdp[5][2000];
618} Ludat3_HijingCommon;
619
620#define LUDAT3_HIJING COMMON_BLOCK(LUDAT3_HIJING,ludat3_hijing)
621COMMON_BLOCK_DEF(Ludat3_HijingCommon,LUDAT3_HIJING);
65a17253 622
7dffa423 623#endif