5 c SUBROUTINE TRANS(icrf,irot)
6 C==> TRANSformation to the Co-moving frame (icrf>0) and
7 C Rotation to the system where (Pt || X),(irot=1).
8 C-FSI ***************************************************
9 IMPLICIT REAL*8 (A-H,O-Z)
10 COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1, ! momenta in NRF
12 COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis.
13 1 X2,Y2,Z2,T2,R2 ! points in NRF
14 COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF
16 COMMON/FSI_POC/AMN,AM1,AM2,CN,C1,C2,AC1,AC2
17 COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
19 C-FSI ***************************************************
20 COMMON /PAIR/P12T,V12Z,GAMZ,V12T,CPHI,SPHI
25 C---> Particle energies ---------
26 P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
27 P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
30 C---> Pair parameters -----------
35 P12S=P12X**2+P12Y**2+P12Z**2
36 P12 =DSQRT(P12S)! Momentum
37 V12 =P12/E12 ! Velocity
38 CTH =P12Z/P12 ! cos(theta)
39 STH =DSQRT(1.D0-CTH**2) !sin
40 V12Z=V12*CTH ! Longit. V
41 GAMZ=1.D0/DSQRT(1.D0-V12Z**2)
42 C-- V12T=V12*STH ! Transv. V in CMS (not needed)
43 P12TS=P12X*P12X+P12Y*P12Y
45 C===> Azimuthal rotation (Pt||X) ============
47 CPHI=P12X/P12T ! cos(phi)
48 SPHI=P12Y/P12T ! sin(phi)
50 CALL ROT8(P1X,P1Y,SPHI,CPHI,P1X,P1Y)
51 CALL ROT8(P2X,P2Y,SPHI,CPHI,P2X,P2Y)
52 CALL ROT8(X1,Y1,SPHI,CPHI,X1,Y1)
53 CALL ROT8(X2,Y2,SPHI,CPHI,X2,Y2)
55 ELSE ! Rotation impossible
57 SPHI=2.D0 ! using it !
59 C===> Co-moving ref. frame ============
61 CALL LTR8(P1Z,E1,V12Z,GAMZ,P1Z,E1a)
62 CALL LTR8(P2Z,E2,V12Z,GAMZ,P2Z,E2a)
63 P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
64 P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
67 CALL LTR8(Z1,T1,V12Z,GAMZ,Z1,T1)
68 CALL LTR8(Z2,T2,V12Z,GAMZ,Z2,T2)
70 C===> Pair reference frame ============
77 P12S=P12X**2+P12Y**2+P12Z**2
82 P112=P1X*P12X+P1Y*P12Y+P1Z*P12Z
90 CVK=(P12X*PPX+P12Y*PPY+P12Z*PPZ)/(P12*AK)
92 V12T=P12T/SQRT(AM12S+P12TS) ! transverse velocity in LCMS
93 C---> Coordinates -----------------------------
98 RS12=XS*P12X+YS*P12Y+ZS*P12Z
112 SUBROUTINE LTR8(Z,T,BETA,GAMMA,ZT,TT)
113 C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(REAL*8)
114 CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
115 COut: ZT,TT- " " after transformation.
116 C==== ===============================================================
117 IMPLICIT REAL*8 (A-H,O-Z)
125 SUBROUTINE LTR4(Z,T,BETA,GAMMA,ZT,TT)
126 C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(real*4)
127 CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
128 COut: ZT,TT- " " after transformation.
129 C==== ===============================================================
136 SUBROUTINE ROT8(X,Y,SF,CF,XR,YR)
137 C===> Rotation with the angle f. (REAL*8)
138 CInp: X,Y-coord. before rotation; SF=sin(f), CF=cos(f),
139 COut: XR,YR - coordinates after rotation.
140 C==== =================================================
141 IMPLICIT REAL*8 (A-H,O-Z)
148 SUBROUTINE SETPDIST(R)
149 C=====Just sets distance between particles
150 IMPLICIT REAL*8 (A-H,O-Z)
151 COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF