More decay options
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 31 Oct 2010 18:38:57 +0000 (18:38 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 31 Oct 2010 18:38:57 +0000 (18:38 +0000)
IHPR2(12)
1: old behavior
2: old behavior + heavy baryons
3: as AliPythia default

HIJING/hijing1_36/hijset.F

index b92cb55..0c2d73a 100644 (file)
@@ -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