Skip to content

Commit 535da60

Browse files
andrzejnovakandrzejnovak
committed
fix: update scale/smear fix (wip)
1 parent f1a536a commit 535da60

File tree

5 files changed

+299
-176
lines changed

5 files changed

+299
-176
lines changed

‎README.md‎

Lines changed: 32 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -61,7 +61,7 @@ combine -M FitDiagnostics --expectSignal 1 -d model_combined.root --cminDefaultM
6161
combine -M FitDiagnostics --expectSignal 1 -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 --saveShapes --saveWithUncertainties -t -1 --toysFrequentist --setParameters z=1
6262
combine -M Significance model_combined.root --expectSignal 1 --redefineSignalPOIs z -t -1 --toysFrequentist
6363
combineTool.py -M AsymptoticLimits -m 125 -d model_combined.root --expectSignal 1 --setParameters z=1 --redefineSignalPOIs r -t -1 --toysFrequentist
64-
python ../HiggsAnalysis/CombinedLimit/test/diffNuisances.py tempModel/fitDiagnostics.root
64+
python ../../../HiggsAnalysis/CombinedLimit/test/diffNuisances.py fitDiagnostics.root
6565
6666
6767
python ../plot.py --data
@@ -98,26 +98,49 @@ plotImpacts.py -i impacts.json -o plots/impacts_out_H --blind
9898
Ensure signal min/max are sufficiently large
9999
```
100100
export BIAS=bias0
101-
combineTool.py -M FitDiagnostics --expectSignal 0 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
101+
combineTool.py -M FitDiagnostics --expectSignal 0 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 20 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
102102
export BIAS=bias1
103-
combineTool.py -M FitDiagnostics --expectSignal 1 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
103+
combineTool.py -M FitDiagnostics --expectSignal 1 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
104104
export BIAS=bias5
105-
combineTool.py -M FitDiagnostics --expectSignal 5 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
105+
combineTool.py -M FitDiagnostics --expectSignal 5 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
106106
export BIAS=bias10
107-
combineTool.py -M FitDiagnostics --expectSignal 10 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
107+
combineTool.py -M FitDiagnostics --expectSignal 10 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
108+
export BIAS=bias30
109+
combineTool.py -M FitDiagnostics --expectSignal 30 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
108110
export BIAS=bias50
109-
combineTool.py -M FitDiagnostics --expectSignal 50 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
111+
combineTool.py -M FitDiagnostics --expectSignal 50 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
112+
export BIAS=bias70
113+
combineTool.py -M FitDiagnostics --expectSignal 70 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
110114
export BIAS=bias100
111-
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:40:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
115+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1:50:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS --random
112116
```
113117
```
114-
for BIAS in bias0 bias1 bias5 bias10 bias50 bias100
118+
for BIAS in bias0 bias1 bias5 bias10 bias30 bias50 bias70 bias100
115119
do
116-
hadd $BIAS.root *Combine$BIAS.*
120+
hadd -f $BIAS.root *Combine$BIAS.*
117121
done
118122
119123
```
120124

125+
```
126+
export BIAS=bias1k
127+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 1000:1020:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
128+
export BIAS=bias5k
129+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 5000:5020:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
130+
export BIAS=bias10k
131+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 10000:10020:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
132+
export BIAS=bias100k
133+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 100000:100020:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
134+
export BIAS=bias500k
135+
combineTool.py -M FitDiagnostics --expectSignal 100 -n $BIAS -d model_combined.root --cminDefaultMinimizerStrategy 0 --robustFit=1 -t 25 -s 500000:500020:1 --toysFrequentist --job-mode condor --sub-opts='+JobFlavour = "workday"' --task-name ggHcc$BIAS
136+
```
137+
```
138+
for BIAS in bias1k bias5k bias10k bias100k bias500k
139+
do
140+
hadd -f $BIAS.root *Combine$BIAS.*
141+
done
142+
143+
```
121144

122145

123146
### Running likelihood scan

‎config_Hxx.py‎

Lines changed: 22 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -57,3 +57,25 @@ def _get(s):
5757
TF_cf_map = np.array(xparlist).reshape(rhodeg + 1, ptdeg + 1)
5858

5959
return TF_cf_map, rhodeg, ptdeg
60+
61+
62+
# SFs
63+
def ddxSF(pbin, flav):
64+
# ptbins = np.array([450, 500, 550, 600, 675, 800, 1200])
65+
_SFdict = {
66+
'cc': {
67+
"SF": [0.899, 1.152, 0.692],
68+
"up": [0.254, 0.428, 0.309],
69+
"down": [0.254, 0.426, 0.278],
70+
}
71+
}
72+
if flav not in ['bb', 'cc', 'qq']:
73+
raise ValueError("``flav`` has be one of ['bb', 'cc', 'qq'].")
74+
if pbin in [0, 1, 2]:
75+
return _SFdict[flav]['SF'][0], _SFdict[flav]['up'][0], _SFdict[flav]['down'][0]
76+
elif pbin in [3, 4]:
77+
return _SFdict[flav]['SF'][1], _SFdict[flav]['up'][1], _SFdict[flav]['down'][1]
78+
elif pbin in [5]:
79+
return _SFdict[flav]['SF'][2], _SFdict[flav]['up'][2], _SFdict[flav]['down'][2]
80+
else:
81+
raise RuntimeError()

0 commit comments

Comments
 (0)