]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PYTHIA6/pyevwt.f
Script to extract ROOT file from ZIP archive in case we got that from the Grid
[u/mrichter/AliRoot.git] / PYTHIA6 / pyevwt.f
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       
10 C...PYEVWT
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
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)
32       COMMON/PYWGHT/WGHTPOW
33       SAVE /PYDAT1/,/PYINT1/,/PYINT2/,/PYWGHT/
34       DATA P02 /25.D0/
35
36       PT2=VINT(48)
37  
38       WTXS=(PT2/P02)**WGHTPOW
39  
40       RETURN
41       END
42  
43 C*********************************************************************