Newer
Older
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
c $Id: newpart.f,v 1.6 1999/01/18 09:57:09 ernst Exp $
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
c include-file newpart
c
cdes this file contains arrays for scattered or new created particles
cdes and is used to communicate between the different routines which
cdes are involved in hadling the kinematics
c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
integer mprt,oprt
c maximum number of new particles:
parameter(mprt=200) ! maximum number of produced particles
parameter(oprt=2) ! maximum number of incoming particles
c pslot : slots of incoming particles
c itypnew: ityps of new particles
c i3new: $2*I_3$ of new particles
c inew: array-indices of new particles (will be assigned in scatter-routines)
c nexit: number of particles in exit-channel
c iline: tag for out-channel process
c nstring1 : number of particles in string 1
c nstring2 : number of particles in string 2
c strcount : ??
c sidnew : stringID for produced particles
c pslot : slots of incoming particles
c itot: $2*I_{tot}$ of new particles (will be assigned in scatter-routines)
c pnew(5,mprt) : momenta, energy and mass of produced particles
c pnew(1,*) = px
c pnew(2,*) = py
c pnew(3,*) = pz
c pnew(4,*) = e
c pnew(5,*) = mass
c xnew(4,mprt) : locations and time of produced particles
c xnew(1,*) = rx
c xnew(2,*) = ry
c xnew(3,*) = rz
c xnew(4,*) = r0
c pold(5,oprt) : momenta, defined like pnew
c itypold : incoming itypes
c iso3old : incoming iso3's
c xtotfacold: xtotfacs of incoming particles
c mstring() masses of strings 1 and 2 (or particles 1 and 2)
c leadfac(mprt): (1-leadfac) is the factor, by which the total cross
c section of a hadron is multiplied within its formation
c time (.ne.1 only for leading hadrons)
integer itypnew(mprt),i3new(mprt),itot(mprt),inew(mprt),nexit
integer nstring1, nstring2,iline,itypold(oprt),iso3old(oprt)
integer pslot(oprt),strcount
real*8 pnew(5,mprt),xnew(4,mprt),mstring(2),leadfac(mprt)
real*8 pold(5,oprt),xtotfacold(oprt)
integer sidnew(mprt)
c relative velocity/between comp. frame and two particle rest frame
c is betax, betay, betaz (needed for lotrans)
c momentum vector in two particle restframe is p0nn,pxnn,pynn,pznn
real*8 betax,betay,betaz,p0nn,pxnn,pynn,pznn,pnn,pnnout
common /inewpart/ itypnew,i3new,itot,inew,nexit,iline,strcount,
& pslot,nstring1,nstring2,sidnew,itypold,iso3old
common /rnewpart/ pnew,xnew,betax,betay,betaz,pold,
& p0nn,pxnn,pynn,pznn,pnn,mstring,pnnout,
& xtotfacold
common /fnewpart/ leadfac