]>
Commit | Line | Data |
---|---|---|
f5ab1a71 | 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 |