How to run generator-level analysis using Pythia with ROOT

setup procedure:

newrel -t 5.1.0 ewk
cd ewk
setup cdfsoft2 5.1.0   ; # can also be 4.10.4
addpkg Stntuple
addpkg Electroweak
make Stntuple.nobin USESHLIBS=1
make Stntuple._fit USESHLIBS=1
make Electroweak.nobin USESHLIBS=1
setenv CTEQ6_DATA    $PWD/Electroweak/cteq6/data
setenv MRST2003_DATA $PWD/Electroweak/mrst2003
cp Electroweak/scripts/.rootrc .

we also need PDFLIB and some of the CERN shared libs: grab libminicern.so and libpdf.so from fcdflnx3:/cdf/scratch/cdfopr/test_ewk/shlibs/Linux2-KCC_4_0 and copy them to shlibs/Linux2-KCC_4_0 in your test release (in this example - ewk) or create a symbolic link

now everything is ready:

stnfit.exe
root [0] .L Electroweak/scripts/pythia_ana.C
root [1] pythia_ana("wenu",30200,1000)

specifying PDF set to pythia_ana: PDF code = 10000*MSTP(52)+MSTP(51)

CTEQ6 : 30200-30241
MRST2003: 40001-40002
PDFLIB codes start from 20000, for example, use 24046 for CTEQ5L

quick validation: see and try to run scripts

Electroweak/scripts/test_cteq6.C
Electroweak/scripts/test_mrst2003.C
Electroweak/scripts/test_pdflib.C

last updated: Dec 15 2003 by Pasha Murat