Update PR task: drathee

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

diff --git a/HIJING/hijing1_36/hijset.F b/HIJING/hijing1_36/hijset.F
index b92cb55215ff9b6fe3e35a73e920e24a117ecc53..65f861d4698d4c5e671398223528af5efb8111d2 100644 (file)
--- a/HIJING/hijing1_36/hijset.F
+++ b/HIJING/hijing1_36/hijset.F
@@ -72,23 +72,58 @@ C
        ENDIF
 
 C...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-.
-       IF(IHPR2(12).GT.0) THEN
+       IF(IHPR2(12).EQ.1) THEN
        CALL LUGIVE_HIJING('MDCY(C111,1)=0')
        CALL LUGIVE_HIJING('MDCY(C310,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C411,1)=0;MDCY(C-411,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C421,1)=0;MDCY(C-421,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C431,1)=0;MDCY(C-431,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C511,1)=0;MDCY(C-511,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C521,1)=0;MDCY(C-521,1)=0')
-          CALL LUGIVE_HIJING('MDCY(C531,1)=0;MDCY(C-531,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C411,1)=0;MDCY(C-411,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C421,1)=0;MDCY(C-421,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C431,1)=0;MDCY(C-431,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C511,1)=0;MDCY(C-511,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C521,1)=0;MDCY(C-521,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C531,1)=0;MDCY(C-531,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3122,1)=0;MDCY(C-3122,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3112,1)=0;MDCY(C-3112,1)=0')
-        CALL LUGIVE_HIJING('MDCY(C3212,1)=0;MDCY(C-3212,1)=0')
+C        CALL LUGIVE_HIJING('MDCY(C3212,1)=0;MDCY(C-3212,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3222,1)=0;MDCY(C-3222,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3312,1)=0;MDCY(C-3312,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3322,1)=0;MDCY(C-3322,1)=0')
         CALL LUGIVE_HIJING('MDCY(C3334,1)=0;MDCY(C-3334,1)=0')
        ENDIF
+       IF(IHPR2(12).EQ.2) THEN
+       CALL LUGIVE_HIJING('MDCY(C111,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C310,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C411,1)=0;MDCY(C-411,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C421,1)=0;MDCY(C-421,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C431,1)=0;MDCY(C-431,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C4122,1)=0;MDCY(C-4122,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C4132,1)=0;MDCY(C-4132,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C4232,1)=0;MDCY(C-4232,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C4332,1)=0;MDCY(C-4332,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C511,1)=0;MDCY(C-511,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C521,1)=0;MDCY(C-521,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C531,1)=0;MDCY(C-531,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C5122,1)=0;MDCY(C-5122,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C5132,1)=0;MDCY(C-5132,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C5232,1)=0;MDCY(C-5232,1)=0')
+       CALL LUGIVE_HIJING('MDCY(C5332,1)=0;MDCY(C-5332,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3122,1)=0;MDCY(C-3122,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3112,1)=0;MDCY(C-3112,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3222,1)=0;MDCY(C-3222,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3312,1)=0;MDCY(C-3312,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3322,1)=0;MDCY(C-3322,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3334,1)=0;MDCY(C-3334,1)=0')
+       ENDIF
+       IF(IHPR2(12).EQ.3) THEN
+       CALL LUGIVE_HIJING('MDCY(C310,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3122,1)=0;MDCY(C-3122,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3112,1)=0;MDCY(C-3112,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3222,1)=0;MDCY(C-3222,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3312,1)=0;MDCY(C-3312,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3322,1)=0;MDCY(C-3322,1)=0')
+        CALL LUGIVE_HIJING('MDCY(C3334,1)=0;MDCY(C-3334,1)=0')
+       ENDIF
+
+
        MSTU(12)=0
        MSTU(21)=1
        IF(IHPR2(10).EQ.0) THEN
@@ -153,8 +188,8 @@ C
           HIDAT(J)=HIDAT0(J,I-1)+(HIDAT0(J,I)-HIDAT0(J,I-1))
      &    *(HINT1(1)-HIDAT0(10,I-1))/(HIDAT0(10,I)-HIDAT0(10,I-1))
 40     CONTINUE
-       HIPR1(31)=HIDAT(5)
-       HIPR1(30)=2.0*HIDAT(5)
+       IF(HIPR1(31).EQ.0) HIPR1(31)=HIDAT(5)
+       HIPR1(30)=2.0*HIPR1(31)
 C
 C
        CALL HIJCRS
@@ -180,5 +215,11 @@ C          ********booking for x distribution of valence quarks
      & 10X,'*',8X,'for ',
      & A4,'(',I3,',',I3,')',' + ',A4,'(',I3,',',I3,')',7X,'*'/
      & 10X,'**************************************************')
+
+       write (*,*),'Average jet cross section (mb): " ', HINT1(10)
+       write (*,*),'Jet cross section (mb): " ', HINT1(11)
+       write (*,*),'Geometrical N size (mb): ', HIPR1(31)
+       write (*,*),'Inelastic NN cross section (mb): ', HINT1(12)
+        write (*,*),'Total NN cross section (mb): ',HINT1(13) 
        RETURN
        END