Do not shift position of nucleons
authorAndreas Morsch <morsch@pcalice102b.(none)>
Mon, 16 Dec 2013 22:30:26 +0000 (23:30 +0100)
committerAndreas Morsch <morsch@pcalice102b.(none)>
Mon, 16 Dec 2013 22:30:26 +0000 (23:30 +0100)
HIJING/hijing1_36/hijing.F

index 23c319a..7561027 100644 (file)
@@ -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