From 8b7ea311ef0421b45e767277e50312ca8bf1f524 Mon Sep 17 00:00:00 2001 From: Andreas Morsch Date: Mon, 16 Dec 2013 23:30:26 +0100 Subject: [PATCH] Do not shift position of nucleons --- HIJING/hijing1_36/hijing.F | 28 +++++++++++++--------------- 1 file changed, 13 insertions(+), 15 deletions(-) diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F index 23c319a89d9..75610273184 100644 --- a/HIJING/hijing1_36/hijing.F +++ b/HIJING/hijing1_36/hijing.F @@ -395,26 +395,24 @@ c *** cl glauber *** xm2=0D0 ym2=0D0 xym=0D0 - +c IF(NCOLT>0) THEN DO 1110 JP=1,IHNT2(1) - YP(1,JP)=YP(1,JP)+BBX - YP(2,JP)=YP(2,JP)+BBY - xmeana=xmeana+YP(1,JP) - ymeana=ymeana+YP(2,JP) + xmeana=xmeana+YP(1,JP)+BBX + ymeana=ymeana+YP(2,JP)+BBY DO 1120 JT=1,IHNT2(3) IF(SCIP2(JP,JT).LT.2.0D0) THEN npart=npart+1 - xmeanp=xmeanp+YP(1,JP) - ymeanp=ymeanp+YP(2,JP) - xm2=xm2+YP(1,JP)*YP(1,JP) - ym2=ym2+YP(2,JP)*YP(2,JP) - xym=xym+YP(1,JP)*YP(2,JP) + xmeanp=xmeanp+YP(1,JP)+BBX + ymeanp=ymeanp+YP(2,JP)+BBY + xm2=xm2+(YP(1,JP)+BBX)*(YP(1,JP)+BBX) + ym2=ym2+(YP(2,JP)+BBY)*(YP(2,JP)+BBY) + xym=xym+(YP(1,JP)+BBX)*(YP(2,JP)+BBY) goto 1110 end if 1120 continue 1110 continue - +c DO 1130 JT=1,IHNT2(3) xmeanb=xmeanb+YT(1,JT) ymeanb=ymeanb+YT(2,JT) @@ -441,11 +439,11 @@ c *** cl glauber *** xm2=xm2/npart ym2=ym2/npart xym=xym/npart - +c sx2=xm2-xmeanp*xmeanp sy2=ym2-ymeanp*ymeanp sxy=xym-xmeanp*ymeanp - +c delx=xmeanb-xmeana dely=ymeanb-ymeana dtmp=delx**2+dely**2 @@ -457,10 +455,10 @@ c *** cl glauber *** eccpart=sqrt(dtmp)/(sx2+sy2) eccmc=(sy2-sx2)/(sy2+sx2) write(*,*),'HOUT: ',bb,' ',bbtrue,' ',ncolt,' ',npart, - 1 ' ',eccrp,' ',eccpart + 1 ' ',eccrp,' ',eccpart, bbx, bby end if end if - +c C ********total number interactions proj and targ has C suffered -- 2.43.0