1 FUNCTION PHOCOR(MPASQR,MCHREN,ME)
2 C.----------------------------------------------------------------------
4 C. PHOTOS: PHOton radiation in decays CORrection weight from
7 C. Purpose: Calculate photon angle. The reshaping functions will
8 C. have to depend on the spin S of the charged particle.
9 C. We define: ME = 2 * S + 1 !
11 C. Input Parameters: MPASQR: Parent mass squared,
12 C. MCHREN: Renormalised mass of charged system,
13 C. ME: 2 * spin + 1 determines matrix element
15 C. Output Parameter: Function value.
17 C. Author(s): Z. Was, B. van Eijk Created at: 26/11/89
18 C. Last Update: 21/03/93
20 C.----------------------------------------------------------------------
22 DOUBLE PRECISION MPASQR,MCHREN,BETA,XX,YY,DATA
24 REAL*8 PHOCOR,PHOFAC,WT1,WT2,WT3
25 DOUBLE PRECISION MCHSQR,MNESQR
27 COMMON/PHOMOM/MCHSQR,MNESQR,PNEUTR(5)
28 DOUBLE PRECISION COSTHG,SINTHG
30 COMMON/PHOPHS/XPHMAX,XPHOTO,COSTHG,SINTHG
33 COMMON/PHOPRO/PROBH,CORWT,XF,IREP
35 C-- Shaping (modified by ZW)...
36 XX=4.D0*MCHSQR/MPASQR*(1.D0-XPHOTO)/(1.D0-XPHOTO+(MCHSQR-MNESQR)/
40 WT3=(1.D0-XPHOTO/XPHMAX)/((1.D0+(1.D0-XPHOTO/XPHMAX)**2)/2.D0)
42 YY=0.5D0*(1.D0-XPHOTO/XPHMAX+1.D0/(1.D0-XPHOTO/XPHMAX))
44 ELSEIF ((ME.EQ.3).OR.(ME.EQ.4).OR.(ME.EQ.5)) THEN
46 WT3=(1.D0+(1.D0-XPHOTO/XPHMAX)**2-(XPHOTO/XPHMAX)**3)/
47 & (1.D0+(1.D0-XPHOTO/XPHMAX)** 2)
50 CALL PHOERR(6,'PHOCOR',DATA)
55 WT1=(1.D0-COSTHG*SQRT(1.D0-MCHREN))/(1.D0-COSTHG*BETA)
56 WT2=(1.D0-XX/YY/(1.D0-BETA**2*COSTHG**2))*(1.D0+COSTHG*BETA)/2.D0
60 IF (PHOCOR.GT.1.D0) THEN
62 CALL PHOERR(3,'PHOCOR',DATA)