[u/mrichter/AliRoot.git] / HERWIG / src / hwuedt.f


CDECK  ID>, HWUEDT.

*CMZ :-        -09/12/91  12.07.08  by  Mike Seymour

*-- Author :    Mike Seymour

C-----------------------------------------------------------------------

      SUBROUTINE HWUEDT(N,IEDT)

C-----------------------------------------------------------------------

C     EDIT THE EVENT RECORD

C     IF N>0 DELETE THE N ENTRIES IN IEDT FROM EVENT RECORD

C     IF N<0 INSERT LINES AFTER THE -N ENTRIES IN IEDT

C-----------------------------------------------------------------------

      INCLUDE 'HERWIG61.INC'

      INTEGER N,IEDT(*),IMAP(0:NMXHEP),IHEP,I,J,I1,I2

      COMMON /HWUMAP/IMAP

C---MOVE ENTRIES AND CALCULATE MAPPING OF POINTERS

      IF (N.EQ.0) THEN

        RETURN

      ELSEIF (N.GT.0) THEN

        I=1

        I1=1

        I2=NHEP

      ELSE

        I=NHEP-N

        I1=NHEP

        I2=1

      ENDIF

      DO 110 IHEP=I1,I2,SIGN(1,I2-I1)

        IMAP(IHEP)=I

        DO 100 J=1,ABS(N)

          IF (IHEP.EQ.IEDT(J)) THEN

            IF (N.GT.0) IMAP(IHEP)=0

            I=I-1

            IF (N.LT.0) IMAP(IHEP)=I

          ENDIF

 100    CONTINUE

        IF (IMAP(IHEP).EQ.I .AND. IHEP.NE.I) THEN

          ISTHEP(I)=ISTHEP(IHEP)

          IDHW(I)=IDHW(IHEP)

          IDHEP(I)=IDHEP(IHEP)

          JMOHEP(1,I)=JMOHEP(1,IHEP)

          JMOHEP(2,I)=JMOHEP(2,IHEP)

          JDAHEP(1,I)=JDAHEP(1,IHEP)

          JDAHEP(2,I)=JDAHEP(2,IHEP)

          CALL HWVEQU(5,PHEP(1,IHEP),PHEP(1,I))

          CALL HWVEQU(4,VHEP(1,IHEP),VHEP(1,I))

          ISTHEP(IHEP)=0

          IDHW(IHEP)=20

          IDHEP(IHEP)=0

          JMOHEP(1,IHEP)=0

          JMOHEP(2,IHEP)=0

          JDAHEP(1,IHEP)=0

          JDAHEP(2,IHEP)=0

          CALL HWVZRO(5,PHEP(1,IHEP))

          CALL HWVZRO(4,VHEP(1,IHEP))

        ENDIF

        I=I+SIGN(1,N)

 110  CONTINUE

      NHEP=NHEP-N

C---RELABEL POINTERS, SETTING ANY WHICH WERE TO DELETED ENTRIES TO ZERO

      IMAP(0)=0

      DO 200 IHEP=1,NHEP

        JMOHEP(1,IHEP)=IMAP(JMOHEP(1,IHEP))

        JMOHEP(2,IHEP)=IMAP(JMOHEP(2,IHEP))

        JDAHEP(1,IHEP)=IMAP(JDAHEP(1,IHEP))

        JDAHEP(2,IHEP)=IMAP(JDAHEP(2,IHEP))

 200  CONTINUE

      END
Commit	Line	Data
3820ca8e	1
	2	CDECK ID>, HWUEDT.
	3
	4	*CMZ :- -09/12/91 12.07.08 by Mike Seymour
	5
	6	*-- Author : Mike Seymour
	7
	8	C-----------------------------------------------------------------------
	9
	10	SUBROUTINE HWUEDT(N,IEDT)
	11
	12	C-----------------------------------------------------------------------
	13
	14	C EDIT THE EVENT RECORD
	15
	16	C IF N>0 DELETE THE N ENTRIES IN IEDT FROM EVENT RECORD
	17
	18	C IF N<0 INSERT LINES AFTER THE -N ENTRIES IN IEDT
	19
	20	C-----------------------------------------------------------------------
	21
	22	INCLUDE 'HERWIG61.INC'
	23
	24	INTEGER N,IEDT(*),IMAP(0:NMXHEP),IHEP,I,J,I1,I2
	25
	26	COMMON /HWUMAP/IMAP
	27
	28	C---MOVE ENTRIES AND CALCULATE MAPPING OF POINTERS
	29
	30	IF (N.EQ.0) THEN
	31
	32	RETURN
	33
	34	ELSEIF (N.GT.0) THEN
	35
	36	I=1
	37
	38	I1=1
	39
	40	I2=NHEP
	41
	42	ELSE
	43
	44	I=NHEP-N
	45
	46	I1=NHEP
	47
	48	I2=1
	49
	50	ENDIF
	51
	52	DO 110 IHEP=I1,I2,SIGN(1,I2-I1)
	53
	54	IMAP(IHEP)=I
	55
	56	DO 100 J=1,ABS(N)
	57
	58	IF (IHEP.EQ.IEDT(J)) THEN
	59
	60	IF (N.GT.0) IMAP(IHEP)=0
	61
	62	I=I-1
	63
	64	IF (N.LT.0) IMAP(IHEP)=I
65
66	ENDIF
67
68	100 CONTINUE
69
70	IF (IMAP(IHEP).EQ.I .AND. IHEP.NE.I) THEN
71
72	ISTHEP(I)=ISTHEP(IHEP)
73
74	IDHW(I)=IDHW(IHEP)
75
76	IDHEP(I)=IDHEP(IHEP)
77
78	JMOHEP(1,I)=JMOHEP(1,IHEP)
79
80	JMOHEP(2,I)=JMOHEP(2,IHEP)
81
82	JDAHEP(1,I)=JDAHEP(1,IHEP)
83
84	JDAHEP(2,I)=JDAHEP(2,IHEP)
85
86	CALL HWVEQU(5,PHEP(1,IHEP),PHEP(1,I))
87
88	CALL HWVEQU(4,VHEP(1,IHEP),VHEP(1,I))
89
90	ISTHEP(IHEP)=0
91
92	IDHW(IHEP)=20
93
94	IDHEP(IHEP)=0
95
96	JMOHEP(1,IHEP)=0
97
98	JMOHEP(2,IHEP)=0
99
100	JDAHEP(1,IHEP)=0
101
102	JDAHEP(2,IHEP)=0
103
104	CALL HWVZRO(5,PHEP(1,IHEP))
105
106	CALL HWVZRO(4,VHEP(1,IHEP))
107
108	ENDIF
109
110	I=I+SIGN(1,N)
111
112	110 CONTINUE
113
114	NHEP=NHEP-N
115
116	C---RELABEL POINTERS, SETTING ANY WHICH WERE TO DELETED ENTRIES TO ZERO
117
118	IMAP(0)=0
119
120	DO 200 IHEP=1,NHEP
121
122	JMOHEP(1,IHEP)=IMAP(JMOHEP(1,IHEP))
123
124	JMOHEP(2,IHEP)=IMAP(JMOHEP(2,IHEP))
125
126	JDAHEP(1,IHEP)=IMAP(JDAHEP(1,IHEP))
127
128	JDAHEP(2,IHEP)=IMAP(JDAHEP(2,IHEP))
129
130	200 CONTINUE
131
132	END