c Version 1.38
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Jun 2010 14:02:13 +0000 (14:02 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Jun 2010 14:02:13 +0000 (14:02 +0000)
c
c     The following common block is added to record the number of elastic
c     (NELT, NELP) and inelastic (NINT, NINP) participants
c
c        COMMON/HIJGLBR/NELT,NINT,NELP,NINP
c        SAVE  /HIJGLBR/

HIJING/hijing/hijglbr.inc [new file with mode: 0644]
HIJING/hijing1_36/hijing.F

diff --git a/HIJING/hijing/hijglbr.inc b/HIJING/hijing/hijglbr.inc
new file mode 100644 (file)
index 0000000..752816c
--- /dev/null
@@ -0,0 +1,2 @@
+        COMMON/HIJGLBR/NELT,NINT,NELP,NINP
+        SAVE  /HIJGLBR/
\ No newline at end of file
index fb077dc..d8f3c38 100644 (file)
@@ -159,6 +159,7 @@ C****************************************************************
 #include "hiparnt.inc"
 C
 #include "hijcrdn.inc"
+#include "hijglbr.inv"
 #include "himain1.inc"
 #include "himain2.inc"
 #include "histrng.inc"
@@ -282,6 +283,10 @@ C                  ********Initialize for a new event
        N01=0
        N10=0
        N11=0
+        NELT=0
+        NINT=0
+        NELP=0
+        NINP=0
        NSG=0
        NCOLT=0
 
@@ -501,6 +506,26 @@ C          ********conduct soft scattering between JP and JT
 
 200    CONTINUE
 
+c
+c**************************
+c
+       DO 201 JP=1,IHNT2(1)
+           IF(NFP(JP,5).GT.2) THEN
+              NINP=NINP+1
+           ELSE IF(NFP(JP,5).EQ.2.OR.NFP(JP,5).EQ.1) THEN
+              NELP=NELP+1
+           ENDIF
+ 201    continue
+       DO 202 JT=1,IHNT2(3)
+           IF(NFT(JT,5).GT.2) THEN
+              NINT=NINT+1
+           ELSE IF(NFT(JT,5).EQ.2.OR.NFT(JT,5).EQ.1) THEN
+              NELT=NELT+1
+           ENDIF
+ 202    continue
+c     
+c*******************************
+
 C********perform jet quenching for jets with PT>HIPR1(11)**********
 
        IF((IHPR2(8).NE.0.OR.IHPR2(3).NE.0).AND.IHPR2(4).GT.0.AND.