Retrofeed from the Release developement
[u/mrichter/AliRoot.git] / HBTAN / ltran12.F
index c3fa195..caa8952 100644 (file)
@@ -1,4 +1,3 @@
-
        SUBROUTINE ltran12
 C
 
@@ -18,67 +17,26 @@ C-FSI ***************************************************
       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
@@ -87,28 +45,47 @@ C===> Pair reference frame       ============
       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.
@@ -150,19 +127,7 @@ 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
-      RPS=R
-      RP=R*R
+      RP=R
+      RPS=R*R
       RETURN
       END
-
-
-
-
-
-
-
-
-
-
-