-
SUBROUTINE ltran12
C
COMMON/FSI_CVK/V,CVK
C-FSI ***************************************************
COMMON /PAIR/P12T,V12Z,GAMZ,V12T,CPHI,SPHI
-
- icrf=1
- irot=1
-
-C---> Particle energies ---------
+
+
+C calculating Ri, Pi and Ei
+ R1=DSQRT(X1*X1+Y1*Y1+Z1*Z1)
+ R2=DSQRT(X2*X2+Y2*Y2+Z2*Z2)
P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
- E1=DSQRT(AM1*AM1+P1S)
- E2=DSQRT(AM2*AM2+P2S)
-C---> Pair parameters -----------
- E12=E1+E2 ! Energy
- P12X=P1X+P2X ! Px
- P12Y=P1Y+P2Y ! Py
- P12Z=P1Z+P2Z ! Pz
- P12S=P12X**2+P12Y**2+P12Z**2
- P12 =DSQRT(P12S)! Momentum
- V12 =P12/E12 ! Velocity
- CTH =P12Z/P12 ! cos(theta)
- STH =DSQRT(1.D0-CTH**2) !sin
- V12Z=V12*CTH ! Longit. V
- GAMZ=1.D0/DSQRT(1.D0-V12Z**2)
-C-- V12T=V12*STH ! Transv. V in CMS (not needed)
- P12TS=P12X*P12X+P12Y*P12Y
- P12T=DSQRT(P12TS) !Pt
-C===> Azimuthal rotation (Pt||X) ============
- IF(V12T.NE.0.D0) THEN
- CPHI=P12X/P12T ! cos(phi)
- SPHI=P12Y/P12T ! sin(phi)
- IF((irot.eq.1)) THEN
- CALL ROT8(P1X,P1Y,SPHI,CPHI,P1X,P1Y)
- CALL ROT8(P2X,P2Y,SPHI,CPHI,P2X,P2Y)
- CALL ROT8(X1,Y1,SPHI,CPHI,X1,Y1)
- CALL ROT8(X2,Y2,SPHI,CPHI,X2,Y2)
- END IF
- ELSE ! Rotation impossible
- CPHI=2.D0 ! to avoid
- SPHI=2.D0 ! using it !
- END IF
-C===> Co-moving ref. frame ============
- IF(icrf.gt.0) THEN
- CALL LTR8(P1Z,E1,V12Z,GAMZ,P1Z,E1a)
- CALL LTR8(P2Z,E2,V12Z,GAMZ,P2Z,E2a)
- P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
- P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
- E1=DSQRT(AM1*AM1+P1S)
- E2=DSQRT(AM2*AM2+P2S)
- CALL LTR8(Z1,T1,V12Z,GAMZ,Z1,T1)
- CALL LTR8(Z2,T2,V12Z,GAMZ,Z2,T2)
- END IF
-C===> Pair reference frame ============
P1=DSQRT(P1S)
P2=DSQRT(P2S)
+ E1=DSQRT(AM1*AM1+P1S)
+ E2=DSQRT(AM2*AM2+P2S)
+C-----------------------------------------------------------------------
E12=E1+E2
P12X=P1X+P2X
P12Y=P1Y+P2Y
P12Z=P1Z+P2Z
P12S=P12X**2+P12Y**2+P12Z**2
- P12 =DSQRT(P12S)
- AM12S=E12*E12-P12S
- AM12=DSQRT(AM12S)
+ AM12=DSQRT(E12**2-P12S)
EPM=E12+AM12
+ P12=DSQRT(P12S)
P112=P1X*P12X+P1Y*P12Y+P1Z*P12Z
H1=(P112/EPM-E1)/AM12
PPX=P1X+P12X*H1
EE=(E12*E1-P112)/AM12
AKS=EE**2-AM1**2
AK=DSQRT(AKS)
- CVK=(P12X*PPX+P12Y*PPY+P12Z*PPZ)/(P12*AK)
- V=P12/E12
- V12T=P12T/SQRT(AM12S+P12TS) ! transverse velocity in LCMS
-C---> Coordinates -----------------------------
+
+ RETURN
+ END
+C====
+C==== ===============================================================
+C==== ===============================================================
+C==== ===============================================================
+
+ subroutine BoostToPrf()
+ IMPLICIT REAL*8 (A-H,O-Z)
+ COMMON/FSI_CVK/V,CVK
+ COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1, !part. momenta in NRF
+ 1 P2X,P2Y,P2Z,E2,P2
+ COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS,
+ 1 X,Y,Z,T,RP,RPS
+ COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis. points in NRF
+ 1 X2,Y2,Z2,T2,R2
+ COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
+
XS=X1-X2
YS=Y1-Y2
ZS=Z1-Z2
- TS=T1-T2
+ TS=T1-T2
RS12=XS*P12X+YS*P12Y+ZS*P12Z
H1=(RS12/EPM-TS)/AM12
-
X=XS+P12X*H1
Y=YS+P12Y*H1
Z=ZS+P12Z*H1
T=(E12*TS-RS12)/AM12
RPS=X*X+Y*Y+Z*Z
RP=DSQRT(RPS)
+CW WRITE(6,38)'RP ',RP,'X ',X,Y,Z,T
+38 FORMAT(A7,E11.4,A7,4E11.4)
+
+ CVK=(P12X*PPX+P12Y*PPY+P12Z*PPZ)/(P12*AK)
+ V=P12/E12
+ return
+ end
+C==== ===============================================================
+C==== ===============================================================
-
- RETURN
- END
-C====
SUBROUTINE LTR8(Z,T,BETA,GAMMA,ZT,TT)
C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(REAL*8)
CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
IMPLICIT REAL*8 (A-H,O-Z)
COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF
1 X,Y,Z,T,RP,RPS
- RPS=R
- RP=R*R
+ RP=R
+ RPS=R*R
RETURN
END
-
-
-
-
-
-
-
-
-
-
-
-