| 1 |
dkralph |
1.2 |
# gg2zz central value is really this:
|
| 2 |
|
|
#pdfset=cteq6ll.LHpdf
|
| 3 |
dkralph |
1.1 |
|
| 4 |
dkralph |
1.2 |
#----------------------------------------------------------------------------------------
|
| 5 |
|
|
# make dummy flat ntuple with information needed to get pdf weights
|
| 6 |
|
|
#----------------------------------------------------------------------------------------
|
| 7 |
|
|
read_bacon() {
|
| 8 |
|
|
exe=./bin/readBacon.exe
|
| 9 |
|
|
|
| 10 |
|
|
in=./data/datafiles-gg2zz.txt
|
| 11 |
|
|
label=gg2zz
|
| 12 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label.root &
|
| 13 |
|
|
|
| 14 |
|
|
in=./data/datafiles-ZZ-powheg.txt
|
| 15 |
|
|
label=ZZ-powheg
|
| 16 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label.root &
|
| 17 |
|
|
}
|
| 18 |
|
|
#----------------------------------------------------------------------------------------
|
| 19 |
|
|
# add weight to ntuples for each pdf set
|
| 20 |
|
|
#----------------------------------------------------------------------------------------
|
| 21 |
|
|
gg2zz_Wgts() {
|
| 22 |
|
|
exe=./bin/pdfreweighter.exe
|
| 23 |
|
|
label=gg2zz
|
| 24 |
|
|
echo ./root/ntuples/$label.root > ./data/$label.txt
|
| 25 |
|
|
in=./data/$label.txt
|
| 26 |
|
|
for imem in {0..40}; do
|
| 27 |
|
|
if [ $imem == 0 ]; then pdfset=cteq6m.LHpdf; else pdfset=cteq6.LHpdf; fi
|
| 28 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label-$pdfset-$imem.root --pdfset $pdfset --imem $imem &
|
| 29 |
|
|
done
|
| 30 |
|
|
}
|
| 31 |
dkralph |
1.1 |
|
| 32 |
dkralph |
1.2 |
ZZ_powheg_Wgts() {
|
| 33 |
|
|
exe=./bin/pdfreweighter.exe
|
| 34 |
|
|
label=ZZ-powheg
|
| 35 |
|
|
echo ./root/ntuples/$label.root > ./data/$label.txt
|
| 36 |
|
|
in=./data/$label.txt
|
| 37 |
|
|
pdfset=CT10.LHgrid
|
| 38 |
|
|
for imem in {0..52}; do
|
| 39 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label-$pdfset-$imem.root --pdfset $pdfset --imem $imem &
|
| 40 |
|
|
done
|
| 41 |
|
|
}
|
| 42 |
|
|
|
| 43 |
|
|
#----------------------------------------------------------------------------------------
|
| 44 |
|
|
# run over all the ntuples generated with previous step, sum up the up and down errors
|
| 45 |
|
|
#----------------------------------------------------------------------------------------
|
| 46 |
|
|
gg2zz_UpDown() {
|
| 47 |
|
|
label=gg2zz
|
| 48 |
|
|
in=./data/pdfset-files.txt; rm -f $in
|
| 49 |
|
|
# note: assumes central member is at top of file
|
| 50 |
|
|
pdfset=cteq6m.LHpdf
|
| 51 |
|
|
for file in `ls ./root/ntuples/* | grep $label-$pdfset-[0-9][0-9]*.root`; do echo $file >> $in; done
|
| 52 |
|
|
pdfset=cteq6.LHpdf
|
| 53 |
|
|
for file in `ls ./root/ntuples/* | grep $label-$pdfset-[0-9][0-9]*.root`; do echo $file >> $in; done
|
| 54 |
|
|
exe=./bin/getUpDowns.exe
|
| 55 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label-out.root
|
| 56 |
|
|
}
|
| 57 |
|
|
|
| 58 |
|
|
ZZ_powheg_UpDown() {
|
| 59 |
|
|
label=ZZ-powheg
|
| 60 |
|
|
pdfset=CT10.LHgrid
|
| 61 |
|
|
in=./data/pdfset-files.txt; rm -f $in
|
| 62 |
|
|
# note: assumes central member is at top of file
|
| 63 |
|
|
for file in `ls ./root/ntuples/* | grep $label-$pdfset-[0-9][0-9]*.root`; do echo $file >> $in; done
|
| 64 |
|
|
exe=./bin/getUpDowns.exe
|
| 65 |
|
|
$exe --inputfile $in --output ./root/ntuples/$label-out.root
|
| 66 |
|
|
}
|
| 67 |
|
|
|
| 68 |
|
|
gg2zz_UpDown |
