[u/mrichter/AliRoot.git] / TDPMjet / DPMcommon.h

#ifndef ROOT_DPMCommon
#define ROOT_DPMCommon

#ifndef __CFORTRAN_LOADED
//*KEEP,cfortran.
#include "cfortran.h"
//*KEND.
#endif

extern "C" {
/*========================================================*/
/* COMMON/DTEVNO/NEVENT,ICASCA				  */
/*--------------------------------------------------------*/
typedef struct {
   Int_t    nevent;
   Int_t    icasca;
} DtevnoCommon;

#define DTEVNO COMMON_BLOCK(DTEVNO,dtevno)
COMMON_BLOCK_DEF(DtevnoCommon,DTEVNO);

/**********************************************************/
/*           D E S C R I P T I O N :                      */
/*--------------------------------------------------------*/
/* 			Event flag		          */
/*========================================================*/

/*========================================================*/
/* COMMON/DTEVT1/NHKK,NEVHKK,ISTHKK(NMXHKK),IDHKK(NMXHKK),*/
/*	JMOHKK(2,NMXHKK),JDAHKK(2,NMXHKK),		  */
/*	PHKK(5,NMXHKK),VHKK(4,NMXHKK),WHKK(4,NMXHKK)	  */
/*--------------------------------------------------------*/
typedef struct {
   Int_t    nhkk;
   Int_t    nevhkk;
   Int_t    isthkk[200000];
   Int_t    idhkk[200000];
   Int_t    jmohkk[200000][2];
   Int_t    jdahkk[200000][2];
   Double_t  phkk[200000][5];
   Double_t  vhkk[200000][4];
   Double_t  whkk[200000][4];
} Dtevt1Common;

#define DTEVT1 COMMON_BLOCK(DTEVT1,dtevt1)
COMMON_BLOCK_DEF(Dtevt1Common,DTEVT1);

/**********************************************************/
/*           D E S C R I P T I O N :                      */
/*--------------------------------------------------------*/
/* 			Event history		          */
/* 						          */
/* NHKK - number of entries in common block		  */
/* NEVHKK - number of the event				  */
/* ISTHKK(i) - status code for entry i with following     */
/*	       meanings:			  	  */
/*	     = 1 final state particle produced in         */
/*	     	photon-/hadron-/nucleon-nucleon collisions*/
/*	     	 or in intranuclear cascade proc.         */
/*	     =-1 nucleons, deuterons, H3, He3, He4        */
/*	     	 evaporated from excited nucleus and      */
/*	     	 photons produced in nuclear deexcitation */
/*	     	 processes;			          */
/*	     = 1001  residual nucleus (ground state).     */
/* IDHKK(i) - particle identity according to PDG code;    */
/*	      for nuclei (evaporation products and     	  */
/* 	      residual nucleus): IDHKK(IHKK)=80000     	  */
/* JMOHKK(1,i) - pointer to the position where the mother */
/*		 is stored; the initial value is 0  	  */
/* JMOHKK(2,i) - pointer to the position of the 2nd mother*/
/* 		 Normally only 1 mother exists, in which  */
/* 		case the value 0 is used. In cluster      */
/*		fragmentation models, the 2 mothers would */
/*		correspond to the q and qbar which join to*/
/*		form a cluster. In string fragmentation,  */
/*		the two mothers of a particle produced in */
/*		the fragmentation would be the 2 endpoints*/
/*		of the string.				  */
/* JDAHKK(1,i) - pointer to the position of the 1st       */
/*		daughter; if an entry has not decayed =0. */
/* JDAHKK(1,i) - pointer to the position of the last      */
/*		daughter; if an entry has not decayed =0. */
/* PHKK(1,i) - momentum in x direction in GeV/c           */
/* PHKK(2,i) - momentum in y direction in GeV/c           */
/* PHKK(3,i) - momentum in z direction in GeV/c           */
/* PHKK(4,i) - energy in GeV			          */
/* PHKK(5,i) - mass in GeV/c^2; for space-like partons    */
/*	    it is allowed to use a mass<0, according      */
/*	    PHKK(5,IHKK) = -sqrt(-m^2)  	          */ 
/* VHKK(1,i) - production vertex in x position in mm      */
/* VHKK(2,i) - production vertex in y position in mm      */
/* VHKK(3,i) - production vertex in z position in mm      */
/* VHKK(4,i) - production time in mm/c (=3.33*10^(-12)s   */
/* WHKK(I,i) - gives positions and times in projectile    */ 
/*		frame, the chains are created on the posi-*/
/*		tions of the projectile nucleons in the   */
/*		projectile frame (target nucleons in target*/
/*		frame); both positions are therefore not  */
/*		completely consistent. The times in the   */
/*		projectile frame are obtained by a Lorentz*/
/*		transformation from the LAB system.	  */
/*========================================================*/

/*========================================================*/
/* COMMON/DTEVT2/IDRES(NMXHKK),IDXRES(NMXHKK),NOBAM(NMXHKK),*/
/*               IDBAM(NMXHKK),IDCH(NMXHKK),NPOINT(10),   */
/*               IHIST(2,NMXHKK)			  */
/*--------------------------------------------------------*/
typedef struct {
   Int_t    idres[200000];
   Int_t    idxres[200000];
   Int_t    nobam[200000];
   Int_t    idbam[200000];
   Int_t    idch[200000];
   Int_t    npoint[10];
   Int_t    ihist[200000][2];
} Dtevt2Common;

#define DTEVT2 COMMON_BLOCK(DTEVT2,dtevt2)
COMMON_BLOCK_DEF(Dtevt2Common,DTEVT2);

/**********************************************************/
/*           D E S C R I P T I O N :                      */
/*--------------------------------------------------------*/
/* 		Extended event history		          */
/* 						          */
/* NMXHKK - max. num. of entries (partons/particles) that */
/* 	  	can be stored in the common block	  */
/* IDRES(IHKK) - mass num. A in case of nuclear fragments */
/*		or residual nuclei (IDHKK(IHKK)=80000).   */
/* IDXRES(IHKK) - charge Zin case of nuclear fragments    */
/*		or residual nuclei (IDHKK(IHKK)=80000).   */
/* NOBAM(IHKK) =1 for particles from proj. fragmentation  */
/* 	       =2 for particles from target fragmentation.*/
/* IDBAM(IHKK) - internal dpmjet particle code(BAMJET code)*/
/*========================================================*/

}

/*========================================================*/
/* COMMON/DTPRTA/IT,ITZ,IP,IPZ,IJPROJ,IBPROJ,IJTARG,IBTARG*/
/*--------------------------------------------------------*/
typedef struct {
   Int_t    it;
   Int_t    itz;
   Int_t    ip;
   Int_t    ipz;
   Int_t    ijproj;
   Int_t    ibproj;
   Int_t    ijtarg;
   Int_t    ibtarg;
} DtprtaCommon;

#define DTPRTA COMMON_BLOCK(DTPRTA,dtprta)
COMMON_BLOCK_DEF(DtprtaCommon,DTPRTA);
/**********************************************************/
/*	     D E S C R I P T I O N :			  */
/*--------------------------------------------------------*/
/* IT, ITZ - nucleon/atomic number of target nucleus	  */
/* IP, IPZ - nucleon/atomic number of projectile nucleus  */
/*	   for incident hadrons IP=IPZ=1		  */
/*========================================================*/

/*========================================================*/
/* COMMON /DTIMPA/ BIMIN,BIMAX,XSFRAC,ICENTR		  */
/*--------------------------------------------------------*/
typedef struct {
   Double_t    bimin;
   Double_t    bimax;
   Double_t    xsfrac;
   Double_t    icent;
} DtimpaCommon;

#define DTIMPA COMMON_BLOCK(DTIMPA,dtimpa)
COMMON_BLOCK_DEF(DtimpaCommon,DTIMPA);
/**********************************************************/
/*	     D E S C R I P T I O N :			  */
/*--------------------------------------------------------*/
/* BIMIN, BIMAX - min., max. b values (default bmin = 0)  */
/* XSFRAC - fraction of x-section (default: 1)  	  */
/* ICENTR =1. central production forced (default: 0)	  */
/*	<0 && >-100 -> bmin = BIMIN, bmax = BIMAX	  */
/*	<-99 -> fraction of x-sec. = XSFRAC		  */
/*	=-1. -> evaporation/fzc suppressed		  */
/*	<-1. -> evaporation/fzc suppressed		  */
/*========================================================*/


/*========================================================*/
/* COMMON /DTGLCP/RPROJ,RTARG,BIMPAC,NWTSAM,NWASAM,NWBSAM,*/
/*		  NWTACC,NWAACC,NWBACC		  	  */
/*--------------------------------------------------------*/
typedef struct {
   Double_t	rproj;
   Double_t	rtarg;
   Double_t	bimpac;
   Int_t	nwtsam;
   Int_t	nwasam;
   Int_t	nwbsam;
   Int_t	nwtacc;
   Int_t	nwtaac;
   Int_t	nwtbac;
} DtglcpCommon;

#define DTGLCP COMMON_BLOCK(DTGLCP,dtglcp)
COMMON_BLOCK_DEF(DtglcpCommon,DTGLCP);
/**********************************************************/
/*	     D E S C R I P T I O N :			  */
/*--------------------------------------------------------*/
/* RPROJ  = radius of projectile nucleus		  */
/* RPROJ  = radius of target nucleus			  */
/* BIMPAC = impact parameter of the collision		  */
/* NWTSAM = total number of wounded nucleons		  */
/* NWASAM = number of wounded nucleons in projectile	  */
/* NWBSAM = number of wounded nucleons in target	  */
/* NWTACC = total number of interacting nucleons	  */
/* NWTAAC = total number of interacting nucleons in proj. */
/* NWTBAC = total number of interacting nucleons in target*/
/*========================================================*/
/*========================================================*/
/* COMMON /POPRCS/IPROCE,IDNODF,IDIFR1,IDIFR2,IDDPOM,     */
/*		  IPRON(15,4)   		  	          */
/*--------------------------------------------------------*/
typedef struct {
   Int_t	iproce;  
   Int_t        idnodf;
   Int_t        idifr1;
   Int_t        idifr2;
   Int_t        iddpom;
   Int_t        ipron[15][4];
} PoprcsCommon;

#define POPRCS COMMON_BLOCK(POPRCS,poprcs)
COMMON_BLOCK_DEF(PoprcsCommon,POPRCS);


/**********************************************************/
/*	     D E S C R I P T I O N :			  */
/*--------------------------------------------------------*/
/*

    COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
     &                JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
     &                VHEP(4,NMXHEP)
*/
const Int_t nmxhep = 4000;

typedef struct {
  Int_t    nevhep;
  Int_t    nhep;
  Int_t    isthep[nmxhep];
  Int_t    idhep [nmxhep];
  Int_t    jmohep[nmxhep][2];
  Int_t    jdahep[nmxhep][2];
  Double_t phep  [nmxhep][5];
  Double_t vhep  [nmxhep][4];
} Poevt1Common;

#define POEVT1 COMMON_BLOCK(POEVT1,poevt1)
COMMON_BLOCK_DEF(Poevt1Common,POEVT1);

#endif
Commit	Line	Data
1a52e0ed	1	#ifndef ROOT_DPMCommon
	2	#define ROOT_DPMCommon
	3
	4	#ifndef __CFORTRAN_LOADED
	5	//*KEEP,cfortran.
	6	#include "cfortran.h"
	7	//*KEND.
	8	#endif
	9
	10	extern "C" {
	11	/========================================================/
	12	/* COMMON/DTEVNO/NEVENT,ICASCA */
	13	/--------------------------------------------------------/
	14	typedef struct {
	15	Int_t nevent;
	16	Int_t icasca;
	17	} DtevnoCommon;
	18
	19	#define DTEVNO COMMON_BLOCK(DTEVNO,dtevno)
	20	COMMON_BLOCK_DEF(DtevnoCommon,DTEVNO);
	21
	22	/**********************************************************/
	23	/* D E S C R I P T I O N : */
	24	/--------------------------------------------------------/
	25	/* Event flag */
	26	/========================================================/
	27
	28	/========================================================/
	29	/* COMMON/DTEVT1/NHKK,NEVHKK,ISTHKK(NMXHKK),IDHKK(NMXHKK),*/
	30	/* JMOHKK(2,NMXHKK),JDAHKK(2,NMXHKK), */
	31	/* PHKK(5,NMXHKK),VHKK(4,NMXHKK),WHKK(4,NMXHKK) */
	32	/--------------------------------------------------------/
	33	typedef struct {
	34	Int_t nhkk;
	35	Int_t nevhkk;
	36	Int_t isthkk[200000];
	37	Int_t idhkk[200000];
	38	Int_t jmohkk[200000][2];
	39	Int_t jdahkk[200000][2];
	40	Double_t phkk[200000][5];
	41	Double_t vhkk[200000][4];
	42	Double_t whkk[200000][4];
	43	} Dtevt1Common;
	44
	45	#define DTEVT1 COMMON_BLOCK(DTEVT1,dtevt1)
	46	COMMON_BLOCK_DEF(Dtevt1Common,DTEVT1);
	47
	48	/**********************************************************/
	49	/* D E S C R I P T I O N : */
	50	/--------------------------------------------------------/
	51	/* Event history */
	52	/* */
	53	/* NHKK - number of entries in common block */
	54	/* NEVHKK - number of the event */
	55	/* ISTHKK(i) - status code for entry i with following */
	56	/* meanings: */
	57	/* = 1 final state particle produced in */
	58	/* photon-/hadron-/nucleon-nucleon collisions*/
	59	/* or in intranuclear cascade proc. */
	60	/* =-1 nucleons, deuterons, H3, He3, He4 */
	61	/* evaporated from excited nucleus and */
	62	/* photons produced in nuclear deexcitation */
	63	/* processes; */
	64	/* = 1001 residual nucleus (ground state). */
65	/* IDHKK(i) - particle identity according to PDG code; */
66	/* for nuclei (evaporation products and */
67	/* residual nucleus): IDHKK(IHKK)=80000 */
68	/* JMOHKK(1,i) - pointer to the position where the mother */
69	/* is stored; the initial value is 0 */
70	/* JMOHKK(2,i) - pointer to the position of the 2nd mother*/
71	/* Normally only 1 mother exists, in which */
72	/* case the value 0 is used. In cluster */
73	/* fragmentation models, the 2 mothers would */
74	/* correspond to the q and qbar which join to*/
75	/* form a cluster. In string fragmentation, */
76	/* the two mothers of a particle produced in */
77	/* the fragmentation would be the 2 endpoints*/
78	/* of the string. */
79	/* JDAHKK(1,i) - pointer to the position of the 1st */
80	/* daughter; if an entry has not decayed =0. */
81	/* JDAHKK(1,i) - pointer to the position of the last */
82	/* daughter; if an entry has not decayed =0. */
83	/* PHKK(1,i) - momentum in x direction in GeV/c */
84	/* PHKK(2,i) - momentum in y direction in GeV/c */
85	/* PHKK(3,i) - momentum in z direction in GeV/c */
86	/* PHKK(4,i) - energy in GeV */
87	/* PHKK(5,i) - mass in GeV/c^2; for space-like partons */
88	/* it is allowed to use a mass<0, according */
89	/* PHKK(5,IHKK) = -sqrt(-m^2) */
90	/* VHKK(1,i) - production vertex in x position in mm */
91	/* VHKK(2,i) - production vertex in y position in mm */
92	/* VHKK(3,i) - production vertex in z position in mm */
93	/* VHKK(4,i) - production time in mm/c (=3.3310^(-12)s /
94	/* WHKK(I,i) - gives positions and times in projectile */
95	/* frame, the chains are created on the posi-*/
96	/* tions of the projectile nucleons in the */
97	/* projectile frame (target nucleons in target*/
98	/* frame); both positions are therefore not */
99	/* completely consistent. The times in the */
100	/* projectile frame are obtained by a Lorentz*/
101	/* transformation from the LAB system. */
102	/========================================================/
103
104	/========================================================/
105	/* COMMON/DTEVT2/IDRES(NMXHKK),IDXRES(NMXHKK),NOBAM(NMXHKK),*/
106	/* IDBAM(NMXHKK),IDCH(NMXHKK),NPOINT(10), */
107	/* IHIST(2,NMXHKK) */
108	/--------------------------------------------------------/
109	typedef struct {
110	Int_t idres[200000];
111	Int_t idxres[200000];
112	Int_t nobam[200000];
113	Int_t idbam[200000];
114	Int_t idch[200000];
115	Int_t npoint[10];
116	Int_t ihist[200000][2];
117	} Dtevt2Common;
118
119	#define DTEVT2 COMMON_BLOCK(DTEVT2,dtevt2)
120	COMMON_BLOCK_DEF(Dtevt2Common,DTEVT2);
121
122	/**********************************************************/
123	/* D E S C R I P T I O N : */
124	/--------------------------------------------------------/
125	/* Extended event history */
126	/* */
127	/* NMXHKK - max. num. of entries (partons/particles) that */
128	/* can be stored in the common block */
129	/* IDRES(IHKK) - mass num. A in case of nuclear fragments */
130	/* or residual nuclei (IDHKK(IHKK)=80000). */
131	/* IDXRES(IHKK) - charge Zin case of nuclear fragments */
132	/* or residual nuclei (IDHKK(IHKK)=80000). */
133	/* NOBAM(IHKK) =1 for particles from proj. fragmentation */
134	/* =2 for particles from target fragmentation.*/
135	/* IDBAM(IHKK) - internal dpmjet particle code(BAMJET code)*/
136	/========================================================/
137
138	}
139
140	/========================================================/
141	/* COMMON/DTPRTA/IT,ITZ,IP,IPZ,IJPROJ,IBPROJ,IJTARG,IBTARG*/
142	/--------------------------------------------------------/
143	typedef struct {
144	Int_t it;
145	Int_t itz;
146	Int_t ip;
147	Int_t ipz;
148	Int_t ijproj;
149	Int_t ibproj;
150	Int_t ijtarg;
151	Int_t ibtarg;
152	} DtprtaCommon;
153
154	#define DTPRTA COMMON_BLOCK(DTPRTA,dtprta)
155	COMMON_BLOCK_DEF(DtprtaCommon,DTPRTA);
156	/**********************************************************/
157	/* D E S C R I P T I O N : */
158	/--------------------------------------------------------/
159	/* IT, ITZ - nucleon/atomic number of target nucleus */
160	/* IP, IPZ - nucleon/atomic number of projectile nucleus */
161	/* for incident hadrons IP=IPZ=1 */
162	/========================================================/
163
164	/========================================================/
165	/* COMMON /DTIMPA/ BIMIN,BIMAX,XSFRAC,ICENTR */
166	/--------------------------------------------------------/
167	typedef struct {
168	Double_t bimin;
169	Double_t bimax;
170	Double_t xsfrac;
171	Double_t icent;
172	} DtimpaCommon;
173
174	#define DTIMPA COMMON_BLOCK(DTIMPA,dtimpa)
175	COMMON_BLOCK_DEF(DtimpaCommon,DTIMPA);
176	/**********************************************************/
177	/* D E S C R I P T I O N : */
178	/--------------------------------------------------------/
179	/* BIMIN, BIMAX - min., max. b values (default bmin = 0) */
180	/* XSFRAC - fraction of x-section (default: 1) */
181	/* ICENTR =1. central production forced (default: 0) */
182	/* <0 && >-100 -> bmin = BIMIN, bmax = BIMAX */
183	/* <-99 -> fraction of x-sec. = XSFRAC */
184	/* =-1. -> evaporation/fzc suppressed */
185	/* <-1. -> evaporation/fzc suppressed */
186	/========================================================/
187
188
189	/========================================================/
190	/* COMMON /DTGLCP/RPROJ,RTARG,BIMPAC,NWTSAM,NWASAM,NWBSAM,*/
191	/* NWTACC,NWAACC,NWBACC */
192	/--------------------------------------------------------/
193	typedef struct {
194	Double_t rproj;
195	Double_t rtarg;
196	Double_t bimpac;
197	Int_t nwtsam;
198	Int_t nwasam;
199	Int_t nwbsam;
200	Int_t nwtacc;
201	Int_t nwtaac;
202	Int_t nwtbac;
203	} DtglcpCommon;
204
205	#define DTGLCP COMMON_BLOCK(DTGLCP,dtglcp)
206	COMMON_BLOCK_DEF(DtglcpCommon,DTGLCP);
207	/**********************************************************/
208	/* D E S C R I P T I O N : */
209	/--------------------------------------------------------/
210	/* RPROJ = radius of projectile nucleus */
211	/* RPROJ = radius of target nucleus */
212	/* BIMPAC = impact parameter of the collision */
213	/* NWTSAM = total number of wounded nucleons */
214	/* NWASAM = number of wounded nucleons in projectile */
215	/* NWBSAM = number of wounded nucleons in target */
216	/* NWTACC = total number of interacting nucleons */
217	/* NWTAAC = total number of interacting nucleons in proj. */
218	/* NWTBAC = total number of interacting nucleons in target*/
219	/========================================================/
9dcc0beb	220	/========================================================/
	221	/* COMMON /POPRCS/IPROCE,IDNODF,IDIFR1,IDIFR2,IDDPOM, */
	222	/* IPRON(15,4) */
	223	/--------------------------------------------------------/
	224	typedef struct {
	225	Int_t iproce;
	226	Int_t idnodf;
	227	Int_t idifr1;
	228	Int_t idifr2;
	229	Int_t iddpom;
	230	Int_t ipron[15][4];
	231	} PoprcsCommon;
	232
	233	#define POPRCS COMMON_BLOCK(POPRCS,poprcs)
	234	COMMON_BLOCK_DEF(PoprcsCommon,POPRCS);
	235
	236
	237	/**********************************************************/
	238	/* D E S C R I P T I O N : */
	239	/--------------------------------------------------------/
fb941ba2	240	/*
	241
	242	COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
	243	& JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
	244	& VHEP(4,NMXHEP)
	245	*/
	246	const Int_t nmxhep = 4000;
	247
	248	typedef struct {
	249	Int_t nevhep;
	250	Int_t nhep;
	251	Int_t isthep[nmxhep];
	252	Int_t idhep [nmxhep];
	253	Int_t jmohep[nmxhep][2];
	254	Int_t jdahep[nmxhep][2];
	255	Double_t phep [nmxhep][5];
	256	Double_t vhep [nmxhep][4];
	257	} Poevt1Common;
	258
	259	#define POEVT1 COMMON_BLOCK(POEVT1,poevt1)
	260	COMMON_BLOCK_DEF(Poevt1Common,POEVT1);
1a52e0ed	261
1a52e0ed	262	#endif