Page MenuHomeHEPForge
Files F10881398

cmpplotv5v6
No OneTemporary
Actions

cmpplotv5v6
View Options

#! /usr/bin/env bash
SCRIPTDIR=$(dirname $0)
## The parton IDs to be compared (include t and tbar if --top flag used)
PIDS="0 1 -1 2 -2 3 -3 4 -4 5 -5"
(echo "$@" | egrep -- "--\<top\>" &> /dev/null) && PIDS="$PIDS 6 -6"
## PDF validation and plotting
for pid in $PIDS; do
## Run comparison calculations for various fixed Q values
for q in 10 50 100 200 500 1000 2000 5000; do
$SCRIPTDIR/cmpv5v6 xf_v*_scanx_q${q}_${pid}.dat -o xfacc_q${q}_${pid}.dat &
done
wait
## Run comparison calculations for various fixed x values
for x in 1e-08 1e-06 0.0001 0.01 0.1 0.2 0.5 0.8; do
$SCRIPTDIR/cmpv5v6 xf_v*_scanq_x${x}_${pid}.dat -o xfacc_x${x}_${pid}.dat &
done
wait
## Plot absolute PDF values and v5/v6 accuracies
$SCRIPTDIR/plotv5v6 xf_v*_scanx_q*_$pid.dat --type=xf --axis=x --log=x --ylim=0: -o cmpv5v6_xf_scanx_$pid.pdf &
$SCRIPTDIR/plotv5v6 xf_v*_scanq_x*_$pid.dat --type=xf --axis=q -o cmpv5v6_xf_scanq_$pid.pdf &
$SCRIPTDIR/plotv5v6 xfacc_q*_$pid.dat --axis=x -o cmpv5v6_xfacc_scanx_$pid.pdf &
$SCRIPTDIR/plotv5v6 xfacc_x*_$pid.dat --axis=q -o cmpv5v6_xfacc_scanq_$pid.pdf &
wait
done
## alpha_s comparisons and plotting
$SCRIPTDIR/cmpv5v6 --as as_v*.dat -o asacc.dat
$SCRIPTDIR/plotv5v6 as_v*.dat --type=as -o cmpv5v6_as.pdf
$SCRIPTDIR/plotv5v6 asacc.dat --type=asdiff -o cmpv5v6_asacc.pdf
## Make a combined PDF of validation
PLOTS="cmpv5v6_as.pdf cmpv5v6_asacc.pdf"
for pid in $PIDS; do
PLOTS="$PLOTS cmpv5v6_xf_scanx_$pid.pdf cmpv5v6_xfacc_scanx_$pid.pdf"
PLOTS="$PLOTS cmpv5v6_xf_scanq_$pid.pdf cmpv5v6_xfacc_scanq_$pid.pdf"
done
##pdftk $PLOTS cat output validation.pdf
gs -dBATCH -dNOPAUSE -q -sDEVICE=pdfwrite -sOutputFile=validation.pdf cmp*.pdf

File Metadata

Mime Type
text/plain
Expires
Sat, May 3, 6:19 AM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4982982
Default Alt Text
cmpplotv5v6 (1 KB)

Event Timeline