C.....MINIMIZE PT1-SQUARED PLUS PT3-SQUARED.....
IF(P1.GE.P3) THEN
- PSI=.5*ULANGL_HIJING(P1**2+P3**2*COS(2.*BET),-P3**2*SIN(2.*BET))
+ PSI=.5*ULANGL_HIJING(P1**2+P3**2*COS(2.*BET),-P3**2*SIN(2.
+ 1 *BET))
PT1=P1*SIN(PSI)
PZ1=P1*COS(PSI)
PT3=P3*SIN(PSI+BET)
PZ3=P3*COS(PSI+BET)
ELSE IF(P3.GT.P1) THEN
- PSI=.5*ULANGL_HIJING(P3**2+P1**2*COS(2.*BET),-P1**2*SIN(2.*BET))
+ PSI=.5*ULANGL_HIJING(P3**2+P1**2*COS(2.*BET),-P1**2*SIN(2.
+ 1 *BET))
PT1=P1*SIN(BET+PSI)
PZ1=-P1*COS(BET+PSI)
PT3=P3*SIN(PSI)
PZ3=-P3*COS(PSI)
ENDIF
-
+
DEL=2.0*HIPR1(40)*RLU_HIJING(0)
P(JL,4)=E1
P(JL,1)=PT1*SIN(DEL)