ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/Shapes/run.sh
(Generate patch)

Comparing UserCode/MitHzz4l/Shapes/run.sh (file contents):
Revision 1.1 by dkralph, Tue Nov 22 20:15:06 2011 UTC vs.
Revision 1.2 by dkralph, Wed Nov 23 15:40:07 2011 UTC

# Line 1 | Line 1
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

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines