Merge branch 'master' of https://git.cern.ch/reps/AliRoot

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

diff --git a/HIJING/hijing1_36/jetini.F b/HIJING/hijing1_36/jetini.F
index b7e277016c04dbb579c0529659383ce3ce909582..3307b7e5cde59d7dc288d82a6b18068e8e823538 100644 (file)
--- a/HIJING/hijing1_36/jetini.F
+++ b/HIJING/hijing1_36/jetini.F
@@ -22,6 +22,7 @@ C*****************************************************************
        CHARACTER BEAM*16,TARG*16
        DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
      &         MINT44(8),MINT45(8)
+#define BLANKET_SAVE
 #include "hijcrdn.inc"
 #include "hiparnt.inc"
 #include "histrng.inc"
@@ -59,6 +60,7 @@ C
            ENDIF
         ENDIF
 c
+
        IF(I_TRIG.NE.0) GO TO 160
         IF(I_TRIG.EQ.I_LAST) GO TO 150
         MSTP(2)=2
@@ -74,6 +76,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
@@ -100,8 +103,8 @@ C           ********No printout of initialization information
         MSUB(28)=1
         MSUB(53)=1
         MSUB(68)=1
-        MSUB(81)=1
-        MSUB(82)=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
@@ -124,6 +127,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
@@ -148,6 +161,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
@@ -214,6 +228,33 @@ 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)
+C ALICE ->
+       IF (IHPR2(49) .EQ. 1) THEN
+C SWITCH OFF MASSIVE FINAL QUARKS 
+          MSEL     = 0
+C QCD on
+          MSUB(11) = 1
+          MSUB(12) = 1
+          MSUB(13) = 1
+          MSUB(28) = 1
+          MSUB(53) = 1
+          MSUB(68) = 1
+C HF OFF
+          MSUB(81) = 0 
+          MSUB(82) = 0
+          MDME(127, 1) = 1
+          MDME(128, 1) = 1
+          MDME(129, 1) = 1
+C g->QQbar off
+          MDME(130, 1) = 0
+          MDME(131, 1) = 0
+C highest flavor in final state shower
+          MSTJ(45) = 3
+C highest flavor in initial state shower and pdf
+          MSTP(54) = 3
+       ENDIF
+C <- ALICE
        CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
        MINT4=MINT(44)
        MINT5=MINT(45)
@@ -231,6 +272,7 @@ C
 C
        IHNT2(16)=0
 C
+        I_LAST=I_TRIG
        RETURN
 C              ********Store the initialization information for
 C                              late use
@@ -254,3 +296,12 @@ C
         MINT(12)=NFT(JT,4)
        RETURN
        END
+
+
+
+
+
+
+
+
+