HBTAN/ltran12.F

   1
   2       SUBROUTINE ltran12
   3 C==>calculating particle momentum in PRF {EE,PPX,PPY,PPZ} from
   4 C-  the momentum of the first particle {E1,P1X,P1Y,P1Z) in NRF
   5       IMPLICIT REAL*8 (A-H,O-Z)
   6       COMMON/FSI_MOM/P1X,P1Y,P1Z,E1,P1,  !part. momenta in NRF
   7      1               P2X,P2Y,P2Z,E2,P2
   8       COMMON/FSI_PRF/PPX,PPY,PPZ,AK,AKS,
   9      1               X,Y,Z,T,RP,RPS
  10       COMMON/FSI_POC/AMN,AM1,AM2,CN,C1,C2,AC1,AC2
  11       COMMON/FSI_P12/P12X,P12Y,P12Z,E12,P12,AM12,EPM
  12       COMMON/FSI_COOR/X1,Y1,Z1,T1,R1, !4-coord. of emis. points in NRF
  13      1                X2,Y2,Z2,T2,R2
  14       COMMON/FSI_CONS/PI,PI2,SPI,DR,W
  15       COMMON/LEDWEIGHT/WEIF,WEI,WEIN,ITEST,IRANPOS
  16
  17
  18 C   fm --> 1/GeV
  19 c      write(*,*)'in LTRAN12 W', W
  20 c      write(*,*)'in LTRAN12 p1x p1y p1z', P1X,P1Y,P1Z
  21 c      write(*,*)'in LTRAN12 p2x p2y p2z', P2X,P2Y,P2Z
  22 c       write(*,*)'in LTRAN12 x y z t', X,Y,Z,T
  23 c      write(*,*)'in LTRAN12 x2 y2 z2', X2,Y2,Z2,T2
  24
  25 CMLV
  26        IF(IRANPOS.EQ.0)THEN
  27       X1=X1*W
  28       Y1=Y1*W
  29       Z1=Z1*W
  30       T1=T1*W
  31       X2=X2*W
  32       Y2=Y2*W
  33       Z2=Z2*W
  34       T2=T2*W
  35 C   calculating Ri, Pi and Ei
  36       R1=DSQRT(X1*X1+Y1*Y1+Z1*Z1)
  37       R2=DSQRT(X2*X2+Y2*Y2+Z2*Z2)
  38       ENDIF
  39 CMLV
  40
  41       P1S=P1X*P1X+P1Y*P1Y+P1Z*P1Z
  42       P2S=P2X*P2X+P2Y*P2Y+P2Z*P2Z
  43       P1=DSQRT(P1S)
  44       P2=DSQRT(P2S)
  45       E1=DSQRT(AM1*AM1+P1S)
  46       E2=DSQRT(AM2*AM2+P2S)
  47 C-----------------------------------------------------------------------
  48       E12=E1+E2
  49       P12X=P1X+P2X
  50       P12Y=P1Y+P2Y
  51       P12Z=P1Z+P2Z
  52       P12S=P12X**2+P12Y**2+P12Z**2
  53       AM12=DSQRT(E12**2-P12S)
  54       EPM=E12+AM12
  55       P12=DSQRT(P12S)
  56       P112=P1X*P12X+P1Y*P12Y+P1Z*P12Z
  57       H1=(P112/EPM-E1)/AM12
  58       PPX=P1X+P12X*H1
  59       PPY=P1Y+P12Y*H1
  60       PPZ=P1Z+P12Z*H1
  61       EE=(E12*E1-P112)/AM12
  62       AKS=EE**2-AM1**2
  63       AK=DSQRT(AKS)
  64 CW      WRITE(6,38)'AK ',AK,'K ',PPX,PPY,PPZ,EE
  65 38    FORMAT(A7,E11.4,A7,4E11.4)
  66       RETURN
  67       END
  68