Difference: Rivet1Howto (2 vs. 3)

Revision 32014-04-23 - CarlosEscobar

Line: 1 to 1
 
META TOPICPARENT name="RivetHowTo"
Line: 10 to 10
 

Producing EVGEN files

Added:
>
>
Information not yet ready frown
 

Running Rivet 1.X

Changed:
<
<
First, log on lxplus5 (I haven't had luck with SLC6 yet) and create a temporal directory to work there:
>
>
First, log on lxplus5 and create a directory to work, let's say a temporal one:
 
tmpdir=$(mktemp -d)
mkdir ${tmpdir}/evgen
Changed:
<
<
Rivet will run over EVGEN files, so one can produce these files or they can be retrieved from the GRID if they already exists. Let's do the latter:
>
>
Rivet will run over EVGEN files, so one can produce these files, they can be retrieved from the GRID or read from any location if they already exists. Let's do the latter and read EVGEN files from EOS here /eos/atlas/user/c/cescobar/top/validation/evgen/mc12_8TeV.117050.PowhegPythia_P2011C_ttbar.evgen.EVNT.v2/
 
Changed:
<
<
cd ${tmpdir}/evgen
export ATLAS_LOCAL_ROOT_BASE=/cvmfs/atlas.cern.ch/repo/ATLASLocalRootBase
alias setupATLAS='source ${ATLAS_LOCAL_ROOT_BASE}/user/atlasLocalSetup.sh'
export ALRB_localConfigDir=$HOME/localConfig
setupATLAS --quiet
localSetupDQ2Client --quiet --skipConfirm
voms-proxy-init -voms atlas
>
>
Now, one can produce the rivet output file (i.e. AIDA file) just running the produce_rivet_validation.py script which is available. It can be retrieved doing this:
 
Deleted:
<
<
Take for example:
 

Changed:
<
<
dq2-get -f "*.root" -n 2 mc12_8TeV.110101.AcerMCPythia_P2011CCTEQ6L1_singletop_tchan_l.evgen.EVNT.e1731/
>
>
cd ${tmpdir} svn export svn+ssh://svn.cern.ch/reps/atlasphys/Physics/Top/Software/MCvalidation/Rivet/Rivet1.X/trunk/scripts/produce_rivet_validation.py
 
Changed:
<
<
Now, one can easily produce the AIDA files just running the runRivetValidation.sh script which is available here. It can be retrieved simply doing:
>
>
Now, just run this script:
 

Changed:
<
<
cd ${tmpdir} svn export svn+ssh://svn.cern.ch/reps/atlas-cescobar/cescobar/TopMC/validation/scripts/runRivetValidation.sh
>
>
python produce_rivet_validation.py -i /eos/atlas/user/c/cescobar/top/validation/evgen/mc12_8TeV.117050.PowhegPythia_P2011C_ttbar.evgen.EVNT.v2/
 
Changed:
<
<
Just run this script in the temporal directory that you just produce:
>
>
This script will run over 1 input file (i.e. 5000 events) in the input directory and will produce an AIDA file with the plots. You can explore the options of the previous script doing:
 
Changed:
<
<
./runRivetValidation.sh
>
>
python produce_rivet_validation.py -h
 
Changed:
<
<
This script will run over all the directories you have in ${tmpdir}/evgen and it will run locally one athena job per directory in parallel. Once all the jobs have finished you will find a AIDA file per DSID in ${tmpdir} directory and also a copy in your EOS area. Now let's produce some plots.
>
>
Assuming that you have produced the AIDA file with the statistics you wanted for all the samples you want to compare, let's now produce some plots.
 

Producing plots using Rivet

Changed:
<
<
Setup rivet:
>
>
Setup Rivet in either lxplus or lxplus5:
 
asetup 17.2.10,here

Line: 64 to 55
 
mkdir plots

Changed:
<
<
compare-histos mc12_8TeV.110101.aida:'Title=AcerMC+Pythia,P2011C,110101' mc12_8TeV.110105.aida:'Title=AcerMC+Pythia,P2011C,MorePS,110105' mc12_8TeV.110106.aida:'Title=AcerMC+Pythia,P2011C,LessPS,110106'
>
>
compare-histos myfile1.aida:'Title=myTitle1' myfile2.aida:'Title=myTitle2'
 make-plots --png *.dat mv MC_* plots
Line: 72 to 63
 Or one can use the rivet-mkhtml utility to produce the plots together with a (ugly) webpage.


Changed:
<
<
rivet-mkhtml mc12_8TeV.110101.aida:'Title=AcerMC+Pythia,P2011C,110101' mc12_8TeV.110105.aida:'Title=AcerMC+Pythia,P2011C,MorePS,110105' mc12_8TeV.110106.aida:'Title=AcerMC+Pythia,P2011C,LessPS,110106' --mc-errs -s
>
>
rivet-mkhtml myfile1.aida:'Title=myTitle1':'LineColor=black':'LineStyle=solid' myfile2.aida:'Title=myTitle2':'LineColor=red':'LineStyle=dashed':LineDash'=1pt 1pt 1pt 1pt' --mc-errs -s
 
Added:
>
>
You can also add options such as --mc-errs, -s, -n 20, --times, etc...
 -- CarlosEscobar - 14-May-2013 \ No newline at end of file
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback