do 1001 i = mul + 1, mul + incmul
ityp0(i) = 21
xmass0(i) = xmp
- call energy(e, temp)
+ call energya(e, temp)
call momntm(px, py, pz, e)
c 7/20/01:
c e = sqrt(e ** 2 + xmp ** 2)
bex = 0d0
bey = 0d0
bez = -tanh(eta(i))
- call lorenz(e, px, py, pz, bex, bey, bez)
+ call lorenza(e, px, py, pz, bex, bey, bez)
px0(i) = pxnew
py0(i) = pynew
pz0(i) = pznew
return
end
- subroutine energy(e, temp)
+ subroutine energya(e, temp)
c to generate the magnitude of the momentum e,
c knowing the temperature of the local thermal distribution temp
cc SAVE /lora/
SAVE
- external lorenz
+ external lorenza
bex = 0d0
bey = 0d0
py1 = gy0(i)
pz1 = gz0(i)
e1 = ft0(i)
- call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+ call lorenza(e1, px1, py1, pz1, bex, bey, bez)
gx0(i) = pxnew
gy0(i) = pynew
gz0(i) = pznew
py1 = py0(i)
pz1 = pz0(i)
e1 = e0(i)
- call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+ call lorenza(e1, px1, py1, pz1, bex, bey, bez)
px0(i) = pxnew
py0(i) = pynew
pz0(i) = pznew
bey = (py1 + py2) / (e1 + e2)
bez = (pz1 + pz2) / (e1 + e2)
- call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+ call lorenza(e1, px1, py1, pz1, bex, bey, bez)
cc SAVE pxnew, ..., values for later use.
px1 = pxnew
py1 = pynew
c we boost to the cm frame, get rotation axis, and rotate 1 particle
c momentum
- call lorenz(t1, x1, y1, z1, bex, bey, bez)
+ call lorenza(t1, x1, y1, z1, bex, bey, bez)
x1 = pxnew
y1 = pynew
z1 = pznew
- call lorenz(t2, x2, y2, z2, bex, bey, bez)
+ call lorenza(t2, x2, y2, z2, bex, bey, bez)
x2 = pxnew
y2 = pynew
e2 = dsqrt(px2**2+py2**2+pz2**2+xmass(jscat)**2)
c boost the 2 particle 4 momentum back to lab frame
- call lorenz(e1, px1, py1, pz1, -bex, -bey, -bez)
+ call lorenza(e1, px1, py1, pz1, -bex, -bey, -bez)
px(iscat) = pxnew
py(iscat) = pynew
pz(iscat) = pznew
e(iscat) = enenew
- call lorenz(e2, px2, py2, pz2, -bex, -bey, -bez)
+ call lorenza(e2, px2, py2, pz2, -bex, -bey, -bez)
px(jscat) = pxnew
py(jscat) = pynew
pz(jscat) = pznew
******************************************************************************
- subroutine lorenz(energy, px, py, pz, bex, bey, bez)
+ subroutine lorenza(energy, px, py, pz, bex, bey, bez)
c add in a cut for beta2 to prevent gam to be nan (infinity)