From 8e529af3a5451eab66f4415cf3f43aeedf7bec1a Mon Sep 17 00:00:00 2001 From: morsch Date: Wed, 9 Jun 2010 14:02:13 +0000 Subject: [PATCH] c Version 1.38 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 | 2 ++ HIJING/hijing1_36/hijing.F | 25 +++++++++++++++++++++++++ 2 files changed, 27 insertions(+) create mode 100644 HIJING/hijing/hijglbr.inc diff --git a/HIJING/hijing/hijglbr.inc b/HIJING/hijing/hijglbr.inc new file mode 100644 index 00000000000..752816cf221 --- /dev/null +++ b/HIJING/hijing/hijglbr.inc @@ -0,0 +1,2 @@ + COMMON/HIJGLBR/NELT,NINT,NELP,NINP + SAVE /HIJGLBR/ \ No newline at end of file diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F index fb077dcbcc4..d8f3c3823f3 100644 --- a/HIJING/hijing1_36/hijing.F +++ b/HIJING/hijing1_36/hijing.F @@ -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. -- 2.39.3