4 c SUBROUTINE TRANS(icrf,irot)
5 C==> TRANSformation to the Co-moving frame (icrf>0) and
6 C Rotation to the system where (Pt || X),(irot=1).
7 C-FSI ***************************************************
8 IMPLICIT REAL*8 (A-H,O-Z)
9 COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1, ! momenta in NRF
11 COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis.
12 1 X2,Y2,Z2,T2,R2 ! points in NRF
13 COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF
15 COMMON/FSI_POC/AMN,AM1,AM2,CN,C1,C2,AC1,AC2
16 COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
18 C-FSI ***************************************************
19 COMMON /PAIR/P12T,V12Z,GAMZ,V12T,CPHI,SPHI
22 C calculating Ri, Pi and Ei
23 R1=DSQRT(X1*X1+Y1*Y1+Z1*Z1)
24 R2=DSQRT(X2*X2+Y2*Y2+Z2*Z2)
25 P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
26 P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
31 C-----------------------------------------------------------------------
36 P12S=P12X**2+P12Y**2+P12Z**2
37 AM12=DSQRT(E12**2-P12S)
40 P112=P1X*P12X+P1Y*P12Y+P1Z*P12Z
52 C==== ===============================================================
53 C==== ===============================================================
54 C==== ===============================================================
56 subroutine BoostToPrf()
57 IMPLICIT REAL*8 (A-H,O-Z)
59 COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1, !part. momenta in NRF
61 COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS,
63 COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis. points in NRF
65 COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
71 RS12=XS*P12X+YS*P12Y+ZS*P12Z
79 CW WRITE(6,38)'RP ',RP,'X ',X,Y,Z,T
80 38 FORMAT(A7,E11.4,A7,4E11.4)
82 CVK=(P12X*PPX+P12Y*PPY+P12Z*PPZ)/(P12*AK)
86 C==== ===============================================================
87 C==== ===============================================================
89 SUBROUTINE LTR8(Z,T,BETA,GAMMA,ZT,TT)
90 C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(REAL*8)
91 CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
92 COut: ZT,TT- " " after transformation.
93 C==== ===============================================================
94 IMPLICIT REAL*8 (A-H,O-Z)
102 SUBROUTINE LTR4(Z,T,BETA,GAMMA,ZT,TT)
103 C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(real*4)
104 CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
105 COut: ZT,TT- " " after transformation.
106 C==== ===============================================================
113 SUBROUTINE ROT8(X,Y,SF,CF,XR,YR)
114 C===> Rotation with the angle f. (REAL*8)
115 CInp: X,Y-coord. before rotation; SF=sin(f), CF=cos(f),
116 COut: XR,YR - coordinates after rotation.
117 C==== =================================================
118 IMPLICIT REAL*8 (A-H,O-Z)
125 SUBROUTINE SETPDIST(R)
126 C=====Just sets distance between particles
127 IMPLICIT REAL*8 (A-H,O-Z)
128 COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF