-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathrunFullChain.sh
executable file
·122 lines (113 loc) · 3.92 KB
/
runFullChain.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
#!/bin/bash
object=$1
year=$2
version=$3
declare -a WPs_TopvsQCD
declare -a WPs_WvsQCD
declare -a WPs_Wvs_QCD_MD
declare -A WPs_FullVer_vs_QCD
declare -A WP_MDVer_vs_QCD
declare -a mist_rates
if [ ${object} == "T" ];
then
category="top"
elif [ ${object} == "W" ];
then
category="w"
fi
if [ ${object} == "T" ];
then
mist_rates=("1p0" "0p5" "0p1")
#mist_rates=("1p0")
if [ ${year} == 2022 ];
then
if [ ${version} == "Nominal" ];
then
WPs_FullVer_vs_QCD=(["1p0"]="0.60" ["0p5"]="0.80" ["0p1"]="0.90") #0: 1p0, 1: 0p5, 2: 0p1 #2015
fi
elif [ ${year} == 2018 ];
then
if [ ${version} == "Nominal" ];
then
WPs_FullVer_vs_QCD=(["1p0"]="0.58" ["0p5"]="0.80" ["0p1"]="0.97") #0: 1p0, 1: 0p5, 2: 0p1
fi
elif [ ${year} == 2017 ];
then
if [ ${version} == "Nominal" ];
then
WPs_FullVer_vs_QCD=(["1p0"]="0.581" ["0p5"]="0.801" ["0p1"]="0.970") #0: 1p0, 1: 0p5, 2: 0p1
fi
elif [ ${year} == 2016 ];
then
if [ ${version} == "Nominal" ];
then
WPs_FullVer_vs_QCD=(["1p0"]="0.490" ["0p5"]="0.738" ["0p1"]="0.957") #0: 1p0, 1: 0p5, 2: 0p1 #2015
#WPs_FullVer_vs_QCD=(["1p0"]="0.495" ["0p5"]="0.733" ["0p1"]="0.958") #0: 1p0, 1: 0p5, 2: 0p1 #2016
fi
fi
elif [ ${object} == "W" ];
then
if [ ${year} == 2022 ];
then
if [ ${version} == "Nominal" ];
then
mist_rates=("5p0" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["5p0"]="0.60" ["1p0"]="0.80" ["0p5"]="0.90") #0: 5p0, 1: 1p0, 2: 0p5
elif [ ${version} == "MD" ];
then
mist_rates=("2p5" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["2p5"]="0.60" ["1p0"]="0.80" ["0p5"]="0.90") #0: 2p5, 1: 1p0, 2: 0p5
fi
elif [ ${year} == 2018 ];
then
if [ ${version} == "Nominal" ];
then
mist_rates=("5p0" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["5p0"]="0.70" ["1p0"]="0.94" ["0p5"]="0.98") #0: 5p0, 1: 1p0, 2: 0p5
elif [ ${version} == "MD" ];
then
mist_rates=("2p5" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["2p5"]="0.59" ["1p0"]="0.82" ["0p5"]="0.90") #0: 2p5, 1: 1p0, 2: 0p5
fi
elif [ ${year} == 2017 ];
then
if [ ${version} == "Nominal" ];
then
mist_rates=("5p0" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["5p0"]="0.709" ["1p0"]="0.944" ["0p5"]="0.978") #0: 5p0, 1: 1p0, 2: 0p5
elif [ ${version} == "MD" ];
then
mist_rates=("2p5" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["2p5"]="0.579" ["1p0"]="0.810" ["0p5"]="0.891") #0: 2p5, 1: 1p0, 2: 0p5
fi
elif [ ${year} == 2016 ];
then
if [ ${version} == "Nominal" ];
then
mist_rates=("5p0" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["5p0"]="0.677" ["1p0"]="0.935" ["0p5"]="0.974") #0: 5p0, 1: 1p0, 2: 0p5 #2015
#WPs_FullVer_vs_QCD=(["5p0"]="0.668" ["1p0"]="0.934" ["0p5"]="0.974") #0: 5p0, 1: 1p0, 2: 0p5 #2016
elif [ ${version} == "MD" ];
then
mist_rates=("2p5" "1p0" "0p5")
WPs_FullVer_vs_QCD=(["2p5"]="0.637" ["1p0"]="0.845" ["0p5"]="0.910") #0: 2p5, 1: 1p0, 2: 0p5 #2015
#WPs_FullVer_vs_QCD=(["2p5"]="0.642" ["1p0"]="0.842" ["0p5"]="0.907") #0: 2p5, 1: 1p0, 2: 0p5 #2016
fi
fi
fi
for era in ${year}
do
for mistRate in "${mist_rates[@]}";
do
unset wpmin
wpmin="${WPs_FullVer_vs_QCD[${mistRate}]}"
cmd_templates2d=$(echo 'make2DTemplates.C("tt1l","'${era}'","'${wpmin}'","1.00")')
cmd_templates1d=$(echo 'HeavyFlavourZCandleStudies.C("'${era}'","tt1l","'${category}'","'${wpmin}'","1.00",false,"pass")')
cmd_datacards=$(echo 'makeDatacards.C("'${era}'","tt1l","'${category}'","'${wpmin}'","1.00")')
cmd_makefits=$(echo 'makeFits.C("'${era}'","'${category}'","'${wpmin}'","1.00","tt1l")')
root -l -q ${cmd_templates2d}
root -l -q ${cmd_templates1d}
root -l -q ${cmd_datacards}
root -l -q ${cmd_makefits} | tee ${object}_${year}_${version}_${wpmin}.txt
done
done