Small modif to switch of gluon radiation and pT kick for triggered jets.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 31 Jan 2002 20:19:51 +0000 (20:19 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 31 Jan 2002 20:19:51 +0000 (20:19 +0000)
HIJING/hijing1_36/jetini.F

index bc31593..d2afed3 100644 (file)
@@ -37,7 +37,6 @@ C
        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
@@ -75,6 +74,7 @@ C                     ********INITIAL STATE RADIATION
 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
@@ -125,6 +125,16 @@ C
         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
@@ -149,6 +159,7 @@ c
            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
@@ -215,6 +226,7 @@ C
 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)
@@ -256,3 +268,12 @@ C
         MINT(12)=NFT(JT,4)
        RETURN
        END
+
+
+
+
+
+
+
+
+