Commits

Anonymous committed 9aab670

fix acerdet-herwig converter

Comments (0)

Files changed (1)

acerdet_src/acerdet.f

  
       INTEGER I
 
-        DO 180 I=1,NHEP
-           K(I,1)=0
-           IF(ISTHEP(I).EQ.1) K(I,1)=1
-           IF(ISTHEP(I).EQ.2) K(I,1)=11      !parton before hadronisation
-           IF(ISTHEP(I).EQ.3) K(I,1)=21
-           IF(ISTHEP(I).GT.10.AND.ISTHEP(I).LT.100) K(I,1)=ISTHEP(I)+30
-           IF(ISTHEP(I).GE.101.AND.ISTHEP(I).LE.103) K(I,1)=21
-           IF(ISTHEP(I).GE.120.AND.ISTHEP(I).LE.124) K(I,1)=21
-           K(I,2)=IDHEP(I)
-           K(I,3)=JMOHEP(1,I)
-           K(I,4)=JDAHEP(1,I)
-           K(I,5)=JDAHEP(2,I)
+      N=NHEP
+      DO 180 I=1,NHEP
+!print *, "==",I,NHEP
+         K(I,1)=0
+         IF(ISTHEP(I).EQ.1) K(I,1)=1
+         IF(ISTHEP(I).EQ.2) K(I,1)=11 !parton before hadronisation
+         IF(ISTHEP(I).EQ.3) K(I,1)=21
+         IF(ISTHEP(I).GT.10.AND.ISTHEP(I).LT.100) K(I,1)=ISTHEP(I)+30
+         IF(ISTHEP(I).GE.101.AND.ISTHEP(I).LE.103) K(I,1)=21
+         IF(ISTHEP(I).GE.120.AND.ISTHEP(I).LE.124) K(I,1)=21
+         K(I,2)=IDHEP(I)
+         K(I,3)=JMOHEP(1,I)
+         K(I,4)=JDAHEP(1,I)
+         K(I,5)=JDAHEP(2,I)
 
-           DO 160 J=1,5
+         DO 160 J=1,5
 
-              P(I,J)=REAL(PHEP(J,I))
- 160       CONTINUE
-           DO 170 J=1,4
-              V(I,J)=REAL(VHEP(J,I))
- 170       CONTINUE
-           V(I,5)=0.
-           IF(ISTHEP(I).EQ.2.AND.PHEP(4,I).GT.PHEP(5,I).AND.
-     +          PHEP(4,I).NE.0.) THEN
-              I1=JDAHEP(1,I)
-              IF(I1.GT.0.AND.I1.LE.NHEP) V(I,5)=(VHEP(4,I1)-VHEP(4,I))
-     &             *PHEP(5,I)/PHEP(4,I)
-           ENDIF
-           
- 180    CONTINUE
+            P(I,J)=REAL(PHEP(J,I))
+ 160     CONTINUE
+         DO 170 J=1,4
+            V(I,J)=REAL(VHEP(J,I))
+ 170     CONTINUE
+         V(I,5)=0.
+         IF(ISTHEP(I).EQ.2.AND.PHEP(4,I).GT.PHEP(5,I).AND.
+     +        PHEP(4,I).NE.0.) THEN
+            I1=JDAHEP(1,I)
+            IF(I1.GT.0.AND.I1.LE.NHEP) V(I,5)=(VHEP(4,I1)-VHEP(4,I))
+     &           *PHEP(5,I)/PHEP(4,I)
+         ENDIF
+         
+ 180  CONTINUE
 
       END
 C