[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / jetini.F

* $Id$
C
C
C
C
C
	SUBROUTINE JETINI(JP,JT,I_TRIG)
C*******Initialize PYTHIA for jet production**********************
C	I_TRIG=0: for normal processes
C	I_TRIG=1: for triggered processes
C       JP: sequence number of the projectile
C       JT: sequence number of the target
C     For A+A collisions, one has to initilize pythia
C     separately for each type of collisions, pp, pn,np and nn,
C     or hp and hn for hA collisions. In this subroutine we use the following
C     catalogue for different type of collisions:
C     h+h: h+h (I_TYPE=1)
C     h+A: h+p (I_TYPE=1), h+n (I_TYPE=2)
C     A+h: p+h (I_TYPE=1), n+h (I_TYPE=2)
C     A+A: p+p (I_TYPE=1), p+n (I_TYPE=2), n+p (I_TYPE=3), n+n (I_TYPE=4)
C*****************************************************************
	CHARACTER BEAM*16,TARG*16
	DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
     &		MINT44(8),MINT45(8)
#include "hijcrdn.inc"
#include "hiparnt.inc"
#include "histrng.inc"
#include "hipyint.inc"
C
#include "ludat1_hijing.inc"
#include "ludat3_hijing.inc"
#include "pysubs_hijing.inc"
#include "pypars_hijing.inc"
#include "pyint1_hijing.inc"
#include "pyint2_hijing.inc"
#include "pyint5_hijing.inc"
	DATA INI/8*0/I_LAST/-1/
	SAVE
C
        IHNT2(11)=JP
        IHNT2(12)=JT
        IF(IHNT2(5).NE.0 .AND. IHNT2(6).NE.0) THEN
           I_TYPE=1
        ELSE IF(IHNT2(5).NE.0 .AND. IHNT2(6).EQ.0) THEN
           I_TYPE=1
           IF(NFT(JT,4).EQ.2112) I_TYPE=2
        ELSE IF(IHNT2(5).EQ.0 .AND. IHNT2(6).NE.0) THEN
           I_TYPE=1
           IF(NFP(JP,4).EQ.2112) I_TYPE=2
        ELSE
           IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2212) THEN
              I_TYPE=1
           ELSE IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2112) THEN
              I_TYPE=2
           ELSE IF(NFP(JP,4).EQ.2112 .AND. NFT(JT,4).EQ.2212) THEN
              I_TYPE=3
           ELSE
              I_TYPE=4
           ENDIF
        ENDIF
c
	IF(I_TRIG.NE.0) GO TO 160
        IF(I_TRIG.EQ.I_LAST) GO TO 150
        MSTP(2)=2
c			********second order running alpha_strong
        MSTP(33)=1
        PARP(31)=HIPR1(17)
C			********inclusion of K factor
        MSTP(51)=3
C			********Duke-Owens set 1 structure functions
        MSTP(61)=1
C			********INITIAL STATE RADIATION
        MSTP(71)=1
C			********FINAL STATE RADIATION
        IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.2) MSTP(61)=0
        IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.1) MSTP(71)=0
	MSTP(91) = 1
c
        MSTP(81)=0
C			******** NO MULTIPLE INTERACTION
        MSTP(82)=1
C			*******STRUCTURE OF MUTLIPLE INTERACTION
        MSTP(111)=0
C		********frag off(have to be done by local call)
        IF(IHPR2(10).EQ.0) MSTP(122)=0
C		********No printout of initialization information
        PARP(81)=HIPR1(8)
        CKIN(5)=HIPR1(8)
        CKIN(3)=HIPR1(8)
        CKIN(4)=HIPR1(9)
        IF(HIPR1(9).LE.HIPR1(8)) CKIN(4)=-1.0
        CKIN(9)=-10.0
        CKIN(10)=10.0
        MSEL=0
        DO 100 ISUB=1,200
           MSUB(ISUB)=0
 100    CONTINUE
        MSUB(11)=1
        MSUB(12)=1
        MSUB(13)=1
        MSUB(28)=1
        MSUB(53)=1
        MSUB(68)=1
        MSUB(81)=1
        MSUB(82)=1
        DO 110 J=1,MIN(8,MDCY(21,3))
 110    MDME(MDCY(21,2)+J-1,1)=0
        ISEL=4
        IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
        MDME(MDCY(21,2)+ISEL-1,1)=1
C			********QCD subprocesses
        MSUB(14)=1
        MSUB(18)=1
        MSUB(29)=1
C                       ******* direct photon production
 150    IF(INI(I_TYPE).NE.0) GO TO 800
	GO TO 400
C
C	*****triggered subprocesses, jet, photon, heavy quark and DY
C
 160    I_TYPE=4+I_TYPE
        IF(I_TRIG.EQ.I_LAST) GO TO 260
        PARP(81)=ABS(HIPR1(10))-0.25
        CKIN(5)=ABS(HIPR1(10))-0.25
        CKIN(3)=ABS(HIPR1(10))-0.25
        CKIN(4)=ABS(HIPR1(10))+0.25
        IF(HIPR1(10).LT.HIPR1(8)) CKIN(4)=-1.0
C	WRITE(6,*) "TRIGGERED", CKIN(3), CKIN(4), CKIN(5), IHPR2(3)
C
C ALICE
C
	IF (HIPR1(50) .EQ. 1) THEN
	   MSTP(61) = 0
	   MSTP(71) = 0
	   MSTP(91) = 0
	ENDIF

c
        MSEL=0
        DO 101 ISUB=1,200
           MSUB(ISUB)=0
 101    CONTINUE
        IF(IHPR2(3).EQ.1) THEN
           MSUB(11)=1
           MSUB(12)=1
           MSUB(13)=1
           MSUB(28)=1
           MSUB(53)=1
           MSUB(68)=1
           MSUB(81)=1
           MSUB(82)=1
           MSUB(14)=1
           MSUB(18)=1
           MSUB(29)=1
           DO 102 J=1,MIN(8,MDCY(21,3))
 102	   MDME(MDCY(21,2)+J-1,1)=0
           ISEL=4
           IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
           MDME(MDCY(21,2)+ISEL-1,1)=1
C			********QCD subprocesses
        ELSE IF(IHPR2(3).EQ.2) THEN
C	   WRITE(6,*) "DIRECT PHOTON PRODUCTION"
           MSUB(14)=1
           MSUB(18)=1
           MSUB(29)=1
C		********Direct photon production
c		q+qbar->g+gamma,q+qbar->gamma+gamma, q+g->q+gamma
        ELSE IF(IHPR2(3).EQ.3) THEN
           CKIN(3)=MAX(0.0,HIPR1(10))
           CKIN(5)=HIPR1(8)
           PARP(81)=HIPR1(8)
           MSUB(81)=1
           MSUB(82)=1
           DO 105 J=1,MIN(8,MDCY(21,3))
 105	   MDME(MDCY(21,2)+J-1,1)=0
           ISEL=4
           IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
           MDME(MDCY(21,2)+ISEL-1,1)=1
C             **********Heavy quark production
        ENDIF
260	IF(INI(I_TYPE).NE.0) GO TO 800
C
C
400	INI(I_TYPE)=1
	IF(IHPR2(10).EQ.0) MSTP(122)=0
	IF(NFP(JP,4).EQ.2212) THEN
		BEAM='P'
	ELSE IF(NFP(JP,4).EQ.-2212) THEN
		BEAM='P~'
	ELSE IF(NFP(JP,4).EQ.2112) THEN
		BEAM='N'
	ELSE IF(NFP(JP,4).EQ.-2112) THEN
		BEAM='N~'
	ELSE IF(NFP(JP,4).EQ.211) THEN
		BEAM='PI+'
	ELSE IF(NFP(JP,4).EQ.-211) THEN
		BEAM='PI-'
	ELSE IF(NFP(JP,4).EQ.321) THEN
		BEAM='PI+'
	ELSE IF(NFP(JP,4).EQ.-321) THEN
		BEAM='PI-'
	ELSE
		WRITE(6,*) 'unavailable beam type', NFP(JP,4)
	ENDIF
	IF(NFT(JT,4).EQ.2212) THEN
		TARG='P'
	ELSE IF(NFT(JT,4).EQ.-2212) THEN
		TARG='P~'
	ELSE IF(NFT(JT,4).EQ.2112) THEN
		TARG='N'
	ELSE IF(NFT(JT,4).EQ.-2112) THEN
		TARG='N~'
	ELSE IF(NFT(JT,4).EQ.211) THEN
		TARG='PI+'
	ELSE IF(NFT(JT,4).EQ.-211) THEN
		TARG='PI-'
	ELSE IF(NFT(JT,4).EQ.321) THEN
		TARG='PI+'
	ELSE IF(NFT(JT,4).EQ.-321) THEN
		TARG='PI-'
	ELSE
		WRITE(6,*) 'unavailable target type', NFT(JT,4)
	ENDIF
C
	IHNT2(16)=1
C       ******************indicate for initialization use when
C                         structure functions are called in PYTHIA
C
C	WRITE(6,*) "CALL PYINIT WITH", CKIN(3), CKIN(4), MSTP(61), MSTP(71), MSTP(91)
	CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
	MINT4=MINT(44)
	MINT5=MINT(45)
	MINT44(I_TYPE)=MINT(44)
	MINT45(I_TYPE)=MINT(45)
	ATXS(0)=XSEC(0,1)
	XSEC0(I_TYPE,0)=XSEC(0,1)
	DO 500 I=1,200
		ATXS(I)=XSEC(I,1)
		XSEC0(I_TYPE,I)=XSEC(I,1)
		DO 500 J=1,20
			ATCO(I,J)=COEF(I,J)
			COEF0(I_TYPE,I,J)=COEF(I,J)
500	CONTINUE
C
	IHNT2(16)=0
C
        I_LAST=I_TRIG
	RETURN
C		********Store the initialization information for
C				late use
C
C
800	MINT(44)=MINT44(I_TYPE)
	MINT(45)=MINT45(I_TYPE)
	MINT4=MINT(44)
	MINT5=MINT(45)
	XSEC(0,1)=XSEC0(I_TYPE,0)
	ATXS(0)=XSEC(0,1)
	DO 900 I=1,200
		XSEC(I,1)=XSEC0(I_TYPE,I)
		ATXS(I)=XSEC(I,1)
	DO 900 J=1,20
		COEF(I,J)=COEF0(I_TYPE,I,J)
		ATCO(I,J)=COEF(I,J)
900	CONTINUE
        I_LAST=I_TRIG
        MINT(11)=NFP(JP,4)
        MINT(12)=NFT(JT,4)
	RETURN
	END
Commit	Line	Data
e74335a4	1	* $Id$
	2	C
	3	C
	4	C
	5	C
	6	C
	7	SUBROUTINE JETINI(JP,JT,I_TRIG)
	8	C*****Initialize PYTHIA for jet production********************
	9	C I_TRIG=0: for normal processes
	10	C I_TRIG=1: for triggered processes
	11	C JP: sequence number of the projectile
	12	C JT: sequence number of the target
	13	C For A+A collisions, one has to initilize pythia
	14	C separately for each type of collisions, pp, pn,np and nn,
	15	C or hp and hn for hA collisions. In this subroutine we use the following
	16	C catalogue for different type of collisions:
	17	C h+h: h+h (I_TYPE=1)
	18	C h+A: h+p (I_TYPE=1), h+n (I_TYPE=2)
	19	C A+h: p+h (I_TYPE=1), n+h (I_TYPE=2)
	20	C A+A: p+p (I_TYPE=1), p+n (I_TYPE=2), n+p (I_TYPE=3), n+n (I_TYPE=4)
	21	C*****************************************************************
	22	CHARACTER BEAM16,TARG16
	23	DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
	24	& MINT44(8),MINT45(8)
	25	#include "hijcrdn.inc"
	26	#include "hiparnt.inc"
	27	#include "histrng.inc"
	28	#include "hipyint.inc"
	29	C
	30	#include "ludat1_hijing.inc"
	31	#include "ludat3_hijing.inc"
	32	#include "pysubs_hijing.inc"
	33	#include "pypars_hijing.inc"
	34	#include "pyint1_hijing.inc"
	35	#include "pyint2_hijing.inc"
	36	#include "pyint5_hijing.inc"
	37	DATA INI/8*0/I_LAST/-1/
	38	SAVE
	39	C
	40	IHNT2(11)=JP
	41	IHNT2(12)=JT
	42	IF(IHNT2(5).NE.0 .AND. IHNT2(6).NE.0) THEN
	43	I_TYPE=1
	44	ELSE IF(IHNT2(5).NE.0 .AND. IHNT2(6).EQ.0) THEN
	45	I_TYPE=1
	46	IF(NFT(JT,4).EQ.2112) I_TYPE=2
	47	ELSE IF(IHNT2(5).EQ.0 .AND. IHNT2(6).NE.0) THEN
	48	I_TYPE=1
	49	IF(NFP(JP,4).EQ.2112) I_TYPE=2
	50	ELSE
	51	IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2212) THEN
	52	I_TYPE=1
	53	ELSE IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2112) THEN
	54	I_TYPE=2
	55	ELSE IF(NFP(JP,4).EQ.2112 .AND. NFT(JT,4).EQ.2212) THEN
	56	I_TYPE=3
	57	ELSE
	58	I_TYPE=4
	59	ENDIF
	60	ENDIF
	61	c
	62	IF(I_TRIG.NE.0) GO TO 160
	63	IF(I_TRIG.EQ.I_LAST) GO TO 150
	64	MSTP(2)=2
65	c ********second order running alpha_strong
66	MSTP(33)=1
67	PARP(31)=HIPR1(17)
68	C ********inclusion of K factor
69	MSTP(51)=3
70	C ********Duke-Owens set 1 structure functions
71	MSTP(61)=1
72	C ********INITIAL STATE RADIATION
73	MSTP(71)=1
74	C ********FINAL STATE RADIATION
75	IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.2) MSTP(61)=0
76	IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.1) MSTP(71)=0
25d2dbca	77	MSTP(91) = 1
e74335a4	78	c
	79	MSTP(81)=0
	80	C ******** NO MULTIPLE INTERACTION
	81	MSTP(82)=1
	82	C *******STRUCTURE OF MUTLIPLE INTERACTION
	83	MSTP(111)=0
	84	C ********frag off(have to be done by local call)
	85	IF(IHPR2(10).EQ.0) MSTP(122)=0
	86	C ********No printout of initialization information
	87	PARP(81)=HIPR1(8)
	88	CKIN(5)=HIPR1(8)
	89	CKIN(3)=HIPR1(8)
	90	CKIN(4)=HIPR1(9)
	91	IF(HIPR1(9).LE.HIPR1(8)) CKIN(4)=-1.0
	92	CKIN(9)=-10.0
	93	CKIN(10)=10.0
	94	MSEL=0
	95	DO 100 ISUB=1,200
	96	MSUB(ISUB)=0
	97	100 CONTINUE
	98	MSUB(11)=1
	99	MSUB(12)=1
	100	MSUB(13)=1
	101	MSUB(28)=1
	102	MSUB(53)=1
	103	MSUB(68)=1
	104	MSUB(81)=1
	105	MSUB(82)=1
	106	DO 110 J=1,MIN(8,MDCY(21,3))
	107	110 MDME(MDCY(21,2)+J-1,1)=0
	108	ISEL=4
	109	IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
	110	MDME(MDCY(21,2)+ISEL-1,1)=1
	111	C ********QCD subprocesses
	112	MSUB(14)=1
	113	MSUB(18)=1
	114	MSUB(29)=1
	115	C ******* direct photon production
	116	150 IF(INI(I_TYPE).NE.0) GO TO 800
	117	GO TO 400
	118	C
	119	C *****triggered subprocesses, jet, photon, heavy quark and DY
	120	C
	121	160 I_TYPE=4+I_TYPE
	122	IF(I_TRIG.EQ.I_LAST) GO TO 260
	123	PARP(81)=ABS(HIPR1(10))-0.25
	124	CKIN(5)=ABS(HIPR1(10))-0.25
	125	CKIN(3)=ABS(HIPR1(10))-0.25
	126	CKIN(4)=ABS(HIPR1(10))+0.25
	127	IF(HIPR1(10).LT.HIPR1(8)) CKIN(4)=-1.0
25d2dbca	128	C WRITE(6,*) "TRIGGERED", CKIN(3), CKIN(4), CKIN(5), IHPR2(3)
	129	C
	130	C ALICE
	131	C
	132	IF (HIPR1(50) .EQ. 1) THEN
	133	MSTP(61) = 0
	134	MSTP(71) = 0
	135	MSTP(91) = 0
	136	ENDIF
	137
e74335a4	138	c
	139	MSEL=0
	140	DO 101 ISUB=1,200
	141	MSUB(ISUB)=0
	142	101 CONTINUE
	143	IF(IHPR2(3).EQ.1) THEN
	144	MSUB(11)=1
	145	MSUB(12)=1
	146	MSUB(13)=1
	147	MSUB(28)=1
	148	MSUB(53)=1
	149	MSUB(68)=1
	150	MSUB(81)=1
	151	MSUB(82)=1
	152	MSUB(14)=1
	153	MSUB(18)=1
	154	MSUB(29)=1
	155	DO 102 J=1,MIN(8,MDCY(21,3))
	156	102 MDME(MDCY(21,2)+J-1,1)=0
	157	ISEL=4
	158	IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
	159	MDME(MDCY(21,2)+ISEL-1,1)=1
	160	C ********QCD subprocesses
	161	ELSE IF(IHPR2(3).EQ.2) THEN
25d2dbca	162	C WRITE(6,*) "DIRECT PHOTON PRODUCTION"
e74335a4	163	MSUB(14)=1
	164	MSUB(18)=1
	165	MSUB(29)=1
	166	C ********Direct photon production
	167	c q+qbar->g+gamma,q+qbar->gamma+gamma, q+g->q+gamma
	168	ELSE IF(IHPR2(3).EQ.3) THEN
	169	CKIN(3)=MAX(0.0,HIPR1(10))
	170	CKIN(5)=HIPR1(8)
	171	PARP(81)=HIPR1(8)
	172	MSUB(81)=1
	173	MSUB(82)=1
	174	DO 105 J=1,MIN(8,MDCY(21,3))
	175	105 MDME(MDCY(21,2)+J-1,1)=0
	176	ISEL=4
	177	IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
	178	MDME(MDCY(21,2)+ISEL-1,1)=1
	179	C **********Heavy quark production
	180	ENDIF
	181	260 IF(INI(I_TYPE).NE.0) GO TO 800
	182	C
	183	C
	184	400 INI(I_TYPE)=1
	185	IF(IHPR2(10).EQ.0) MSTP(122)=0
	186	IF(NFP(JP,4).EQ.2212) THEN
	187	BEAM='P'
	188	ELSE IF(NFP(JP,4).EQ.-2212) THEN
	189	BEAM='P~'
	190	ELSE IF(NFP(JP,4).EQ.2112) THEN
	191	BEAM='N'
	192	ELSE IF(NFP(JP,4).EQ.-2112) THEN
	193	BEAM='N~'
	194	ELSE IF(NFP(JP,4).EQ.211) THEN
	195	BEAM='PI+'
	196	ELSE IF(NFP(JP,4).EQ.-211) THEN
	197	BEAM='PI-'
	198	ELSE IF(NFP(JP,4).EQ.321) THEN
	199	BEAM='PI+'
	200	ELSE IF(NFP(JP,4).EQ.-321) THEN
	201	BEAM='PI-'
	202	ELSE
	203	WRITE(6,*) 'unavailable beam type', NFP(JP,4)
	204	ENDIF
	205	IF(NFT(JT,4).EQ.2212) THEN
	206	TARG='P'
	207	ELSE IF(NFT(JT,4).EQ.-2212) THEN
	208	TARG='P~'
	209	ELSE IF(NFT(JT,4).EQ.2112) THEN
	210	TARG='N'
	211	ELSE IF(NFT(JT,4).EQ.-2112) THEN
	212	TARG='N~'
	213	ELSE IF(NFT(JT,4).EQ.211) THEN
	214	TARG='PI+'
	215	ELSE IF(NFT(JT,4).EQ.-211) THEN
	216	TARG='PI-'
	217	ELSE IF(NFT(JT,4).EQ.321) THEN
	218	TARG='PI+'
	219	ELSE IF(NFT(JT,4).EQ.-321) THEN
	220	TARG='PI-'
	221	ELSE
	222	WRITE(6,*) 'unavailable target type', NFT(JT,4)
	223	ENDIF
	224	C
	225	IHNT2(16)=1
	226	C ******************indicate for initialization use when
227	C structure functions are called in PYTHIA
228	C
25d2dbca	229	C WRITE(6,*) "CALL PYINIT WITH", CKIN(3), CKIN(4), MSTP(61), MSTP(71), MSTP(91)
e74335a4	230	CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
	231	MINT4=MINT(44)
	232	MINT5=MINT(45)
	233	MINT44(I_TYPE)=MINT(44)
	234	MINT45(I_TYPE)=MINT(45)
	235	ATXS(0)=XSEC(0,1)
	236	XSEC0(I_TYPE,0)=XSEC(0,1)
	237	DO 500 I=1,200
	238	ATXS(I)=XSEC(I,1)
	239	XSEC0(I_TYPE,I)=XSEC(I,1)
	240	DO 500 J=1,20
	241	ATCO(I,J)=COEF(I,J)
	242	COEF0(I_TYPE,I,J)=COEF(I,J)
	243	500 CONTINUE
	244	C
	245	IHNT2(16)=0
	246	C
57367a3d	247	I_LAST=I_TRIG
e74335a4	248	RETURN
	249	C ********Store the initialization information for
	250	C late use
	251	C
	252	C
	253	800 MINT(44)=MINT44(I_TYPE)
	254	MINT(45)=MINT45(I_TYPE)
	255	MINT4=MINT(44)
	256	MINT5=MINT(45)
	257	XSEC(0,1)=XSEC0(I_TYPE,0)
	258	ATXS(0)=XSEC(0,1)
	259	DO 900 I=1,200
	260	XSEC(I,1)=XSEC0(I_TYPE,I)
	261	ATXS(I)=XSEC(I,1)
	262	DO 900 J=1,20
	263	COEF(I,J)=COEF0(I_TYPE,I,J)
	264	ATCO(I,J)=COEF(I,J)
	265	900 CONTINUE
	266	I_LAST=I_TRIG
	267	MINT(11)=NFP(JP,4)
	268	MINT(12)=NFT(JT,4)
	269	RETURN
	270	END
25d2dbca	271
	272
	273
	274
	275
	276
	277
	278
	279