1 |
< |
in=./data/flatfiles.txt |
2 |
< |
mc="--mc" |
3 |
< |
exe=./bin/pdfreweighter.exe |
4 |
< |
label=foo |
1 |
> |
# gg2zz central value is really this: |
2 |
> |
#pdfset=cteq6ll.LHpdf |
3 |
|
|
4 |
< |
for imem in {11..52}; do |
5 |
< |
break |
6 |
< |
$exe $mc --inputfile $in --output ./root/ntuples/$label-$imem.root --imem $imem & |
7 |
< |
done |
4 |
> |
#---------------------------------------------------------------------------------------- |
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 |
|
|
32 |
< |
in=./data/pdfset-files.txt |
33 |
< |
rm -f $in |
34 |
< |
for file in `ls ./root/ntuples/* | grep $label-[0-9][0-9]*.root`; do echo $file >> $in; done |
35 |
< |
exe=./bin/getUpDowns.exe |
36 |
< |
$exe $mc --inputfile $in --output ./root/ntuples/$label-out.root |
32 |
> |
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 |