Retrofeed from the Release developement
[u/mrichter/AliRoot.git] / HBTAN / ltran12.F
1        SUBROUTINE ltran12
2 C
3
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
10      1               P2X,P2Y,P2Z,E2,P2
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
14      1               X,Y,Z,T,RP,RPS
15       COMMON/FSI_POC/AMN,AM1,AM2,CN,C1,C2,AC1,AC2
16       COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
17       COMMON/FSI_CVK/V,CVK
18 C-FSI ***************************************************
19       COMMON /PAIR/P12T,V12Z,GAMZ,V12T,CPHI,SPHI
20  
21  
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
27       P1=DSQRT(P1S)
28       P2=DSQRT(P2S)
29       E1=DSQRT(AM1*AM1+P1S)
30       E2=DSQRT(AM2*AM2+P2S)
31 C-----------------------------------------------------------------------
32       E12=E1+E2
33       P12X=P1X+P2X
34       P12Y=P1Y+P2Y
35       P12Z=P1Z+P2Z
36       P12S=P12X**2+P12Y**2+P12Z**2
37       AM12=DSQRT(E12**2-P12S)
38       EPM=E12+AM12
39       P12=DSQRT(P12S)
40       P112=P1X*P12X+P1Y*P12Y+P1Z*P12Z
41       H1=(P112/EPM-E1)/AM12
42       PPX=P1X+P12X*H1
43       PPY=P1Y+P12Y*H1
44       PPZ=P1Z+P12Z*H1
45       EE=(E12*E1-P112)/AM12
46       AKS=EE**2-AM1**2
47       AK=DSQRT(AKS)
48
49       RETURN
50       END
51 C====       
52 C==== ===============================================================
53 C==== ===============================================================
54 C==== ===============================================================
55
56       subroutine BoostToPrf()
57       IMPLICIT REAL*8 (A-H,O-Z)
58       COMMON/FSI_CVK/V,CVK
59       COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1,  !part. momenta in NRF
60      1               P2X,P2Y,P2Z,E2,P2
61       COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS,
62      1               X,Y,Z,T,RP,RPS
63       COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, ! 4-coord. of emis. points in NRF
64      1                X2,Y2,Z2,T2,R2
65       COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
66
67       XS=X1-X2
68       YS=Y1-Y2
69       ZS=Z1-Z2
70       TS=T1-T2
71       RS12=XS*P12X+YS*P12Y+ZS*P12Z
72       H1=(RS12/EPM-TS)/AM12
73       X=XS+P12X*H1
74       Y=YS+P12Y*H1
75       Z=ZS+P12Z*H1
76       T=(E12*TS-RS12)/AM12
77       RPS=X*X+Y*Y+Z*Z
78       RP=DSQRT(RPS)
79 CW      WRITE(6,38)'RP ',RP,'X ',X,Y,Z,T
80 38    FORMAT(A7,E11.4,A7,4E11.4)
81  
82       CVK=(P12X*PPX+P12Y*PPY+P12Z*PPZ)/(P12*AK)
83       V=P12/E12
84       return 
85       end
86 C==== ===============================================================
87 C==== ===============================================================
88
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)
95       ZH=GAMMA*(Z-BETA*T)     
96       TT=GAMMA*(T-BETA*Z)
97       ZT=ZH
98       RETURN
99       END
100 C====       
101       
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==== ===============================================================
107       ZH=GAMMA*(Z-BETA*T)     
108       TT=GAMMA*(T-BETA*Z)
109       ZT=ZH
110       RETURN
111       END
112 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)
119       XH=X*CF+Y*SF !Y              
120       YR=Y*CF-X*SF !   _-X'   
121       XR=XH        ! _- f      
122       RETURN       !------>
123       END          !     X
124
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
129      1               X,Y,Z,T,RP,RPS
130       RP=R
131       RPS=R*R
132       RETURN
133       END