+++ /dev/null
- SUBROUTINE ltran12
-C
-
-c SUBROUTINE TRANS(icrf,irot)
-C==> TRANSformation to the Co-moving frame (icrf>0) and
-C Rotation to the system where (Pt || X),(irot=1).
-C-FSI ***************************************************
- IMPLICIT REAL*8 (A-H,O-Z)
- COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1, ! momenta in NRF
- 1 P2X,P2Y,P2Z,E2,P2
- COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis.
- 1 X2,Y2,Z2,T2,R2 ! points in NRF
- COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS, ! momenta in PRF
- 1 X,Y,Z,T,RP,RPS
- COMMON/FSI_POC/AMN,AM1,AM2,CN,C1,C2,AC1,AC2
- COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
- COMMON/FSI_CVK/V,CVK
-C-FSI ***************************************************
- COMMON /PAIR/P12T,V12Z,GAMZ,V12T,CPHI,SPHI
-
-
-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
- 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
- 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
- PPY=P1Y+P12Y*H1
- PPZ=P1Z+P12Z*H1
- EE=(E12*E1-P112)/AM12
- AKS=EE**2-AM1**2
- AK=DSQRT(AKS)
-
- 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
- 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==== ===============================================================
-
- 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.
-COut: ZT,TT- " " after transformation.
-C==== ===============================================================
- IMPLICIT REAL*8 (A-H,O-Z)
- ZH=GAMMA*(Z-BETA*T)
- TT=GAMMA*(T-BETA*Z)
- ZT=ZH
- RETURN
- END
-C====
-
- SUBROUTINE LTR4(Z,T,BETA,GAMMA,ZT,TT)
-C===> Lorentz Transf. of Z(Pz) and T(E) to moving ref. frame.(real*4)
-CInp: Z,T-Zcoord,Time before tr., BETA,GAMMA- velocity, Lor.fact.
-COut: ZT,TT- " " after transformation.
-C==== ===============================================================
- ZH=GAMMA*(Z-BETA*T)
- TT=GAMMA*(T-BETA*Z)
- ZT=ZH
- RETURN
- END
-C====
- SUBROUTINE ROT8(X,Y,SF,CF,XR,YR)
-C===> Rotation with the angle f. (REAL*8)
-CInp: X,Y-coord. before rotation; SF=sin(f), CF=cos(f),
-COut: XR,YR - coordinates after rotation.
-C==== =================================================
- IMPLICIT REAL*8 (A-H,O-Z)
- XH=X*CF+Y*SF !Y
- YR=Y*CF-X*SF ! _-X'
- XR=XH ! _- f
- RETURN !------>
- END ! X
-
- SUBROUTINE SETPDIST(R)
-C=====Just sets distance between particles
- 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
- RP=R
- RPS=R*R
- RETURN
- END
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff