]>
Commit | Line | Data |
---|---|---|
1d319207 | 1 | SUBROUTINE SETPOWWGHT(POWER) |
2 | IMPLICIT DOUBLE PRECISION(A-H, O-Z) | |
3 | COMMON/PYWGHT/WGHTPOW | |
4 | SAVE /PYWGHT/ | |
5 | C.. Divide power by 2 because we get p_T^2 from the Pythia | |
6 | WGHTPOW = POWER/2.D0 | |
7 | RETURN | |
8 | END | |
9 | ||
58f5aaa6 | 10 | C...PYEVWT |
1d319207 | 11 | C...Implementation routine for event weighting |
12 | C...Uses a pt-hard based weighting with w ~ (p_T,hard/p0)^n | |
13 | C...Where n can be set using a common block variable | |
14 | C...p0 = 5 GeV/c | |
15 | C...Multiplies the | |
58f5aaa6 | 16 | C...standard PYTHIA differential cross-section by a process- and |
17 | C...kinematics-dependent factor WTXS. For MSTP(142)=1 this corresponds | |
18 | C...to generation of weighted events, with weight 1/WTXS, while for | |
19 | C...MSTP(142)=2 it corresponds to a modification of the underlying | |
20 | C...physics. | |
21 | ||
22 | SUBROUTINE PYEVWT(WTXS) | |
23 | ||
24 | C...Double precision and integer declarations. | |
25 | IMPLICIT DOUBLE PRECISION(A-H, O-Z) | |
26 | IMPLICIT INTEGER(I-N) | |
27 | INTEGER PYK,PYCHGE,PYCOMP | |
28 | C...Commonblocks. | |
29 | COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) | |
30 | COMMON/PYINT1/MINT(400),VINT(400) | |
31 | COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2) | |
1d319207 | 32 | COMMON/PYWGHT/WGHTPOW |
33 | SAVE /PYDAT1/,/PYINT1/,/PYINT2/,/PYWGHT/ | |
34 | DATA P02 /25.D0/ | |
35 | ||
58f5aaa6 | 36 | PT2=VINT(48) |
37 | ||
1d319207 | 38 | WTXS=(PT2/P02)**WGHTPOW |
58f5aaa6 | 39 | |
40 | RETURN | |
41 | END | |
1d319207 | 42 | |
43 | C********************************************************************* |