|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
SUBROUTINE DBPART
C Print the HEPEUP.
IMPLICIT NONE
integer iup
INTEGER MAXNUP
PARAMETER (MAXNUP=500)
INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP,SPINUP
COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP,
& IDUP(MAXNUP),ISTUP(MAXNUP),MOTHUP(2,MAXNUP),
& ICOLUP(2,MAXNUP),PUP(5,MAXNUP),VTIMUP(MAXNUP),
& SPINUP(MAXNUP)
WRITE(42,*) '++++++++++++++++'
WRITE(42,*) NUP
DO iup=1,NUP
WRITE(42,*) 'Particle #',iup,': IDUP =',IDUP(iup),
& ' ISTUP =',ISTUP(iup),
& ' MOTHUP(1) =',MOTHUP(1,iup),
& ' MOTHUP(2) =',MOTHUP(2,iup),
& ' ICOLUP(1) =',ICOLUP(1,iup),
& ' ICOLUP(2) =',ICOLUP(2,iup)
ENDDO
END
SUBROUTINE DBINRD
C Print the list of INORAD'd particles
IMPLICIT NONE
INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
DOUBLE PRECISION PHEP,VHEP
PARAMETER (NMXHEP=4000)
COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
&JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
INTEGER IUP,IEND,INORAD
INTEGER MAXNUP
PARAMETER (MAXNUP=500)
COMMON/SHVETO/IEND,INORAD(MAXNUP)
WRITE(43,*) '++++++++++++++++'
DO IUP=1,NHEP
WRITE(43,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
& ' ISTHEP =',ISTHEP(iup),
& ' INORAD =',INORAD(iup)
ENDDO
END
SUBROUTINE DBCAL
C Print the particles as they go into the calorimeter.
IMPLICIT NONE
INTEGER NMXHEP,NEVHEP,NHEP,ISTHEP,IDHEP,JMOHEP,JDAHEP
DOUBLE PRECISION PHEP,VHEP
PARAMETER (NMXHEP=4000)
COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
&JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
INTEGER IUP,IEND,INORAD
INTEGER MAXNUP
PARAMETER (MAXNUP=500)
COMMON/SHVETO/IEND,INORAD(MAXNUP)
WRITE(44,*) '++++++++++++++++'
DO IUP=1,NHEP
WRITE(44,*) 'Particle #',iup,': IDHEP =',IDHEP(iup),
& ' ISTHEP =',ISTHEP(iup)
ENDDO
END
|