Alpgen Studies for W+HF

To do

• automatically obtain xsec (metadata?), either in Showering or Analysis
• delete tmp files of Showering
• create random number files for Production
• make analysis jobOptions flexible
• root script to analyse and scale output

Documentation

Scripts taken from https://twiki.cern.ch/twiki/bin/view/AtlasProtected/AlpgenProductionScripts

more info: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/AlpGenForAtlas

FAQ: https://twiki.cern.ch/twiki/bin/view/Main/AlpgenFAQ

Systematic studies: https://twiki.cern.ch/twiki/bin/view/Main/AtlasAlpgenWJets

Event generation

package: https://svnweb.cern.ch/trac/atlasusr/browser/clange/WplusHF/AlpgenProduction

check-out:

svn co svn+ssh://svn.cern.ch/reps/atlasusr/clange/WplusHF/AlpgenProduction/trunk AlpgenProduction

install Alpgen:

source installAlpgen.sh

setup for NAF batch system

source setup.sh

create a few soft-links for large file output:

ln -s /scratch/zn/current/atlas/clange/WplusHFOutput OutputFiles
ln -s /scratch/zn/current/atlas2/monthly/clange/WplusHFOutput_Consolidated OutputFiles_Consolidated

running Alpgen

Scripts are set up in submit, e.g. do:

source submit/submitNominal.sh

This will submit a lot of batch jobs and return tar balls of those.

merging of Alpgen jobs

./consolidate WJet0

Note: check that this gawk error is not a problem:

./consolidate: INFO: Found 100 tar files matching base name
Now processing tar files....................................................................................................Done
./consolidate: INFO: Unweighted events consolidated in to file OutputFiles_Consolidated//WJet1.unw
./consolidate: INFO: Adding following MLM matching settings based on alpgen default choice for the chosen generator settings:
./consolidate: INFO: ETCLUS=20. RCLUS=0.7 ETAMAX=6. IEXCL=1
gawk: (FILENAME=- FNR=16) fatal: division by zero attempted
./consolidate: INFO: Parameter file consolidated in to file OutputFiles_Consolidated//WJet1_unw.par
./consolidate: INFO: Consolidated file has 10355997 events, representing a total luminosity of 1143.64 pb-1 (avg. cross section = 9055.89 pb)

packaging files for Evgen_trf

Determine from .unw file how many lines are taken up by each Alpgen job

./package -n 5 -out OutputFiles_Consolidated/packaged/user.clange.alpgen.907710.WlnuNp0_pt20_7TeV.TXT.v1/user.clange.alpgen.907710.WlnuNp0_pt20_7TeV.TXT.v1. -evtsize 5 -evtcount 1000 OutputFiles_Consolidated/WJet0

Showering

package: https://svnweb.cern.ch/trac/atlasusr/browser/clange/WplusHF/Showering

check-out:

svn co svn+ssh://svn.cern.ch/reps/atlasusr/clange/WplusHF/Showering/trunk Showering

put jobOptions into jobOptions dir using subdirectory for dataset number, e.g.

jobOptions/907710/MC10.907710.AlpgenJimmyWlnuNp0_pt20.py

link input files:

ln -s /scratch/zn/current/atlas2/monthly/clange/WplusHFOutput_Consolidated/packaged input/MC10

link output files and logs:

ln -s /scratch/zn/current/atlas2/monthly/clange/WplusHFOutputPool output
ln -s /scratch/zn/current/atlas2/monthly/clange/WplusHFOutputPoolLogs log

The log files are needed to obtain sample cross-sections.

Analysis

package: https://svnweb.cern.ch/trac/atlasusr/browser/clange/WplusHF/AlpgenWJetsStudy

check-out:

svn co svn+ssh://svn.cern.ch/reps/atlasusr/clange/WplusHF/AlpgenWJetsStudy/trunk AlpgenWJetsStudy

Old Framework (deprecated)

export CVSROOT=:gserver:isscvs.cern.ch:/local/reps/qmulatlas
kinit clange@CERN.CH
klog clange@IFH.DE
cd /afs/ifh.de/group/atlas/users/scratch/clange/alpgenstudies/
cvs co vbjet_project
export ALPGEN=/afs/ifh.de/group/atlas/users/scratch/clange/alpgenstudies/alpgen
cd vbjet_project
./install.sh $ALPGEN
ini root528_64
./buildWorkDir.sh $ALPGEN workarea wjet

/afs/cern.ch/sw/lcg/external/cernlib/2005/x86_64-slc5-gcc43-opt/lib