| 3820ca8e |
1 | |
| 2 | CDECK ID>, HWUMBW. |
| 3 | |
| 4 | *CMZ :- -21/02/98 11.11.56 by Bryan Webber |
| 5 | |
| 6 | *-- Author : Bryan Webber |
| 7 | |
| 8 | C----------------------------------------------------------------------- |
| 9 | |
| 10 | FUNCTION HWUMBW(ID) |
| 11 | |
| 12 | C----------------------------------------------------------------------- |
| 13 | |
| 14 | C CHOOSES MASS ACCORDING TO BREIT-WIGNER DISTRIBUTION |
| 15 | |
| 16 | C----------------------------------------------------------------------- |
| 17 | |
| 18 | INCLUDE 'HERWIG61.INC' |
| 19 | |
| 20 | DOUBLE PRECISION HWUMBW,HWR,WMX,TAU,T |
| 21 | |
| 22 | INTEGER ID |
| 23 | |
| 24 | C--WMX IS MAX NUMBER OF WIDTHS FROM NOMINAL MASS |
| 25 | |
| 26 | DATA WMX/10D0/ |
| 27 | |
| 28 | SAVE WMX |
| 29 | |
| 30 | HWUMBW=RMASS(ID) |
| 31 | |
| 32 | TAU=RLTIM(ID) |
| 33 | |
| 34 | IF (TAU.EQ.ZERO.OR.TAU.GT.1D-18) RETURN |
| 35 | |
| 36 | 1 T=HALF*TAN(PIFAC*(HWR()-HALF)) |
| 37 | |
| 38 | IF (ABS(T).GT.WMX) GO TO 1 |
| 39 | |
| 40 | HWUMBW=HWUMBW+HBAR*T/TAU |
| 41 | |
| 42 | IF(HWUMBW.LT.ZERO) GOTO 1 |
| 43 | |
| 44 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame