CombineHarvester
HttSystematics_MSSMUpdate.cc
Go to the documentation of this file.
1 #include "CombineHarvester/CombineTools/interface/HttSystematics.h"
2 #include <vector>
3 #include <string>
4 #include "CombineHarvester/CombineTools/interface/Systematics.h"
5 #include "CombineHarvester/CombineTools/interface/Process.h"
6 #include "CombineHarvester/CombineTools/interface/Utilities.h"
7 
8 namespace ch {
9 
10 using ch::syst::SystMap;
11 using ch::syst::SystMapAsymm;
12 using ch::syst::era;
13 using ch::syst::channel;
14 using ch::syst::bin_id;
15 using ch::syst::process;
16 using ch::JoinStr;
17 
18 void AddMSSMUpdateSystematics_et_mt(CombineHarvester & cb, CombineHarvester src) {
19  //CombineHarvester src = cb.cp();
20  src.channel({"et", "mt"});
21 
22  // This regular expression will match any Higgs signal process.
23  // It's useful for catching all these processes when we don't know
24  // which of them will be signal or background
25  //std::string bb_rgx = "(bb[hHA])*";
26  //std::string gg_rgx = "(gg[hHA])*";
27 
28  std::vector<std::string> ggH = {"ggH", "ggH_Htautau", "ggA_Atautau", "ggh_htautau"};
29  std::vector<std::string> bbH = {"bbH", "bbH_Htautau", "bbA_Atautau", "bbh_htautau"};
30 
31  auto signal = Set2Vec(src.cp().signals().SetFromProcs(
32  std::mem_fn(&Process::process)));
33 
34  src.cp().process(JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}))
35  .AddSyst(cb, "lumi_$ERA", "lnN", SystMap<era>::init
36  //({"7TeV"}, 1.026)
37  ({"8TeV"}, 1.026));
38 
39  src.cp().process({"ZTT"})
40  .AddSyst(cb, "CMS_eff_m", "lnN", SystMap<channel>::init({"mt"}, 1.010));
41  src.cp().process({"ZTT"})
42  .AddSyst(cb, "CMS_eff_e", "lnN", SystMap<channel>::init({"et"}, 1.020));
43 
44  src.cp()
45  .AddSyst(cb, "CMS_eff_t_$CHANNEL_$ERA", "lnN", SystMap<bin_id, process>::init
46  ({10, 11, 12, 13, 14}, JoinStr({signal,{"ZTT", "VV", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.08)
47  ({10, 11, 12}, {"TT"}, 1.080));
48 
49  src.cp()
50  .AddSyst(cb, "CMS_eff_t_mssmHigh_mutau_$ERA", "shape", SystMap<channel, process>::init
51  ({"mt"}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1));
52 
53  src.cp()
54  .AddSyst(cb, "CMS_eff_t_mssmHigh_etau_$ERA", "shape", SystMap<channel, process>::init
55  ({"et"}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1));
56 
57 // src.cp().process(JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}))
58 // .AddSyst(cb, "CMS_eff_t_mssmHigh_$CHANNEL_$ERA", "shape", SystMap<>::init(1));
59 
60  src.cp()
61  .AddSyst(cb, "CMS_scale_t_mutau_$ERA", "shape", SystMap<channel, process>::init
62  ({"mt"}, JoinStr({signal,{"ZTT", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1));
63 
64  src.cp()
65  .AddSyst(cb, "CMS_scale_t_etau_$ERA", "shape", SystMap<channel, process>::init
66  ({"et"}, JoinStr({signal,{"ZTT", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1));
67 
68 // src.cp().process(JoinStr({signal,{"ZTT", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
69 // .AddSyst(cb, "CMS_scale_t_$CHANNEL_$ERA", "shape", SystMap<channel>::init(1));
70 
72  src.cp().process({"W"})
73  .AddSyst(cb, "CMS_shift1_muTau_nobtag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {10}, 1.000));
75  src.cp().process({"W"})
76  .AddSyst(cb, "CMS_shift2_muTau_nobtag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {10}, 1.000));
77  src.cp().process({"QCD"})
78  .AddSyst(cb, "CMS_shift1_muTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {10}, 1.000));
79  src.cp().process({"QCD"})
80  .AddSyst(cb, "CMS_shift2_muTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {10}, 1.000));
81  src.cp().process({"W"})
82  .AddSyst(cb, "CMS_shift1_muTau_nobtag_medium_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {11}, 1.000));
83  src.cp().process({"W"})
84  .AddSyst(cb, "CMS_shift2_muTau_nobtag_medium_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {11}, 1.000));
85  src.cp().process({"QCD"})
86  .AddSyst(cb, "CMS_shift1_muTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {11}, 1.000));
87  src.cp().process({"QCD"})
88  .AddSyst(cb, "CMS_shift2_muTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {11}, 1.000));
89  src.cp().process({"W"})
90  .AddSyst(cb, "CMS_shift1_muTau_nobtag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {12}, 1.000));
91  src.cp().process({"W"})
92  .AddSyst(cb, "CMS_shift2_muTau_nobtag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {12}, 1.000));
93  src.cp().process({"QCD"})
94  .AddSyst(cb, "CMS_shift1_muTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {12}, 1.000));
95  src.cp().process({"QCD"})
96  .AddSyst(cb, "CMS_shift2_muTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {12}, 1.000));
97  src.cp().process({"W"})
98  .AddSyst(cb, "CMS_shift1_muTau_btag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {13}, 1.000));
99  src.cp().process({"W"})
100  .AddSyst(cb, "CMS_shift2_muTau_btag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {13}, 1.000));
101  src.cp().process({"QCD"})
102  .AddSyst(cb, "CMS_shift1_muTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {13}, 1.000));
103  src.cp().process({"QCD"})
104  .AddSyst(cb, "CMS_shift2_muTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {13}, 1.000));
105  src.cp().process({"W"})
106  .AddSyst(cb, "CMS_shift1_muTau_btag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {14}, 1.000));
107  src.cp().process({"W"})
108  .AddSyst(cb, "CMS_shift2_muTau_btag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {14}, 1.000));
109  src.cp().process({"QCD"})
110  .AddSyst(cb, "CMS_shift1_muTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {14}, 1.000));
111  src.cp().process({"QCD"})
112  .AddSyst(cb, "CMS_shift2_muTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"mt"}, {14}, 1.000));
113  src.cp().process({"W"})
114  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {10}, 1.000));
115  src.cp().process({"W"})
116  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {10}, 1.000));
117  src.cp().process({"QCD"})
118  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {10}, 1.000));
119  src.cp().process({"QCD"})
120  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {10}, 1.000));
121  src.cp().process({"W"})
122  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_medium_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {11}, 1.000));
123  src.cp().process({"W"})
124  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_medium_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {11}, 1.000));
125  src.cp().process({"QCD"})
126  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {11}, 1.000));
127  src.cp().process({"QCD"})
128  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {11}, 1.000));
129  src.cp().process({"W"})
130  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {12}, 1.000));
131  src.cp().process({"W"})
132  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {12}, 1.000));
133  src.cp().process({"QCD"})
134  .AddSyst(cb, "CMS_shift1_eleTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {12}, 1.000));
135  src.cp().process({"QCD"})
136  .AddSyst(cb, "CMS_shift2_eleTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {12}, 1.000));
137  src.cp().process({"W"})
138  .AddSyst(cb, "CMS_shift1_eleTau_btag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {13}, 1.000));
139  src.cp().process({"W"})
140  .AddSyst(cb, "CMS_shift2_eleTau_btag_low_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {13}, 1.000));
141  src.cp().process({"QCD"})
142  .AddSyst(cb, "CMS_shift1_eleTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {13}, 1.000));
143  src.cp().process({"QCD"})
144  .AddSyst(cb, "CMS_shift2_eleTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {13}, 1.000));
145  src.cp().process({"W"})
146  .AddSyst(cb, "CMS_shift1_eleTau_btag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {14}, 1.000));
147  src.cp().process({"W"})
148  .AddSyst(cb, "CMS_shift2_eleTau_btag_high_$ERA_W_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {14}, 1.000));
149  src.cp().process({"QCD"})
150  .AddSyst(cb, "CMS_shift1_eleTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {14}, 1.000));
151  src.cp().process({"QCD"})
152  .AddSyst(cb, "CMS_shift2_eleTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<channel, bin_id>::init({"et"}, {14}, 1.000));
153 
154  src.cp()
155  .AddSyst(cb, "CMS_scale_j_$ERA", "lnN", SystMap<channel, bin_id, process>::init
156  ({"mt", "et"}, {10, 11, 12}, JoinStr({bbH, {"TT"}}), 0.99)
157  ({"mt"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 0.99)
158  ({"mt"}, {13, 14}, {bbH}, 1.01)
159  ({"mt"}, {13, 14}, {"ZJ"}, 0.980)
160  ({"mt"}, {13, 14}, {"ZL"}, 0.980)
161  ({"mt"}, {13, 14}, {"TT"}, 0.920)
162  ({"mt"}, {13, 14}, {"VV"}, 0.980)
163  ({"et"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 1.03)
164  ({"et"}, {13, 14}, {bbH}, 1.01)
165  ({"et"}, {13, 14}, {"ZJ"}, 1.060)
166  ({"et"}, {13, 14}, {"ZL"}, 0.960)
167  ({"et"}, {13, 14}, {"TT"}, 0.900)
168  ({"et"}, {13, 14}, {"VV"}, 0.940));
169 
170  src.cp()
171  .AddSyst(cb, "CMS_htt_scale_met_$ERA", "lnN", SystMap<bin_id, process>::init
172  ({10, 11, 12, 13, 14}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.99)
173  ({10, 11, 12, 13, 14}, {"ZJ", "ZL", "TT"}, 1.010)
174  ({13, 14}, {"W"}, 1.010));
175 
176  src.cp()
177  .AddSyst(cb, "CMS_eff_b_$ERA", "lnN", SystMap<channel, bin_id, process>::init
178  ({"et", "mt"}, {10, 11, 12}, {bbH}, 0.99)
179  ({"et", "mt"}, {10, 11, 12}, {"TT"}, 0.950)
180  ({"mt"}, {10, 11, 12}, {"VV"}, 0.980)
181  ({"et"}, {10, 11, 12}, {"VV"}, 0.990)
182  ({"mt"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 1.01)
183  ({"mt"}, {13, 14}, {bbH}, 1.03)
184  ({"mt"}, {13, 14}, {"ZL"}, 1.040)
185  ({"mt"}, {13, 14}, {"TT"}, 1.020)
186  ({"mt"}, {13, 14}, {"VV"}, 1.040)
187  ({"et"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 1.01)
188  ({"et"}, {13, 14}, {bbH}, 1.02)
189  ({"et"}, {13, 14}, {"ZL"}, 1.020)
190  ({"et"}, {13, 14}, {"ZJ"}, 1.040)
191  ({"et"}, {13, 14}, {"TT"}, 1.040)
192  ({"et"}, {13, 14}, {"VV"}, 1.040));
193 
194  src.cp()
195  .AddSyst(cb, "CMS_fake_b_$ERA", "lnN", SystMap<channel, bin_id, process>::init
196  ({"mt", "et"}, {10, 11, 12}, {"TT"}, 0.970)
197  ({"mt", "et"}, {10, 11, 12}, {"VV"}, 0.990)
198  ({"mt"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 1.05)
199  ({"mt"}, {13, 14}, {"ZL"}, 1.050)
200  ({"mt"}, {13, 14}, {"ZJ"}, 1.090)
201  ({"mt"}, {13, 14}, {"TT", "VV"}, 1.010)
202  ({"et"}, {13, 14}, JoinStr({ggH, {"ggH_SM125"}}), 1.03)
203  ({"et"}, {13, 14}, {"ZL"}, 1.020)
204  ({"et"}, {13, 14}, {"ZJ"}, 1.090)
205  ({"et"}, {13, 14}, {"TT", "VV"}, 1.010));
206 
207  src.cp().process({"ZTT", "ZJ", "ZL"})
208  .AddSyst(cb, "CMS_htt_zttNorm_$ERA", "lnN", SystMap<>::init(1.030));
209 
210  src.cp().process({"ZTT"})
211  .AddSyst(cb, "CMS_htt_extrap_ztt_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
212  ({11, 12, 13, 14}, 1.050));
213 
214  src.cp().process({"TT"})
215  .AddSyst(cb, "CMS_htt_ttbarNorm_$ERA", "lnN", SystMap<>::init(1.100));
216 
217  src.cp().process({"TT"})
218  .AddSyst(cb, "CMS_htt_ttbarPtReweight_$ERA", "shape", SystMap<>::init(1.000));
219 
220  src.cp().process({"TT"})
221  .AddSyst(cb, "CMS_htt_ttbarJetFake_$ERA", "shape", SystMap<>::init(1.000));
222 
223  src.cp().process({"TT"})
224  .AddSyst(cb, "CMS_htt_ttbar_emb_$ERA", "lnN", SystMap<bin_id>::init({13, 14}, 1.140));
225 
226  src.cp().process({"VV"})
227  .AddSyst(cb, "CMS_htt_DiBosonNorm_$ERA", "lnN", SystMap<>::init(1.150));
228 
229  src.cp().process({"W"})
230  .AddSyst(cb, "CMS_htt_WNorm_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
231  ({10, 11, 12}, 1.200)
232  ({13}, 1.300)
233  ({14}, 1.500));
234 
235  src.cp().process({"QCD"})
236  .AddSyst(cb, "CMS_htt_QCDSyst_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
237  ({10, 11, 12}, 1.100)
238  ({13, 14}, 1.200));
239 
240  src.cp().process({"QCD"})
241  .AddSyst(cb, "CMS_htt_QCDfrShape_mutau_$ERA", "shape", SystMap<channel>::init({"mt"}, 1.000));
242  src.cp().process({"QCD"})
243  .AddSyst(cb, "CMS_htt_QCDfrShape_etau_$ERA", "shape", SystMap<channel>::init({"et"}, 1.000));
244 
245  src.cp().process({"ZJ"})
246  .AddSyst(cb, "CMS_htt_ZJetFakeTau_$CHANNEL_$ERA", "lnN", SystMap<>::init(1.200));
247 
248  src.cp().process({"ZL"})
249  .AddSyst(cb, "CMS_htt_ZLeptonFakeTau_$CHANNEL_low_pTtau_$ERA", "lnN", SystMap<channel, bin_id>::init
250  ({"mt"}, {10, 13}, 1.300)
251  ({"et"}, {10, 13}, 1.200));
252 
253  src.cp().process({"ZL"})
254  .AddSyst(cb, "CMS_htt_ZLeptonFakeTau_$CHANNEL_medium_pTtau_$ERA", "lnN", SystMap<channel, bin_id>::init
255  ({"mt"}, {11, 14}, 1.300)
256  ({"et"}, {11, 14}, 1.200));
257 
258  src.cp().process({"ZL"})
259  .AddSyst(cb, "CMS_htt_ZLeptonFakeTau_$CHANNEL_high_pTtau_$ERA", "lnN", SystMap<channel, bin_id>::init
260  ({"mt"}, {12}, 1.300)
261  ({"et"}, {12}, 1.200));
262 
263  src.cp().process({"ZL"}).channel({"mt"})
264  .AddSyst(cb, "CMS_htt_ZLScale_mutau_$ERA", "shape", SystMap<>::init(1.000));
265  src.cp().process({"ZL"}).channel({"et"})
266  .AddSyst(cb, "CMS_htt_ZLScale_etau_$ERA", "shape", SystMap<>::init(1.000));
267 
268  src.cp().process({"W"})
269  .AddSyst(cb, "CMS_htt_WShape_mutau_nobtag_low_$ERA", "shape", SystMap<channel, bin_id>::init({"mt"}, {10}, 1.000));
270  src.cp().process({"W"})
271  .AddSyst(cb, "CMS_htt_WShape_mutau_nobtag_medium_$ERA", "shape", SystMap<channel, bin_id>::init({"mt"}, {11}, 1.000));
272  src.cp().process({"W"})
273  .AddSyst(cb, "CMS_htt_WShape_mutau_nobtag_high_$ERA", "shape", SystMap<channel, bin_id>::init({"mt"}, {12}, 1.000));
274  src.cp().process({"W"})
275  .AddSyst(cb, "CMS_htt_WShape_mutau_btag_low_$ERA", "shape", SystMap<channel, bin_id>::init({"mt"}, {13}, 1.000));
276  src.cp().process({"W"})
277  .AddSyst(cb, "CMS_htt_WShape_mutau_btag_high_$ERA", "shape", SystMap<channel, bin_id>::init({"mt"}, {14}, 1.000));
278  src.cp().process({"W"})
279  .AddSyst(cb, "CMS_htt_WShape_etau_nobtag_low_$ERA", "shape", SystMap<channel, bin_id>::init({"et"}, {10}, 1.000));
280  src.cp().process({"W"})
281  .AddSyst(cb, "CMS_htt_WShape_etau_nobtag_medium_$ERA", "shape", SystMap<channel, bin_id>::init({"et"}, {11}, 1.000));
282  src.cp().process({"W"})
283  .AddSyst(cb, "CMS_htt_WShape_etau_nobtag_high_$ERA", "shape", SystMap<channel, bin_id>::init({"et"}, {12}, 1.000));
284  src.cp().process({"W"})
285  .AddSyst(cb, "CMS_htt_WShape_etau_btag_low_$ERA", "shape", SystMap<channel, bin_id>::init({"et"}, {13}, 1.000));
286  src.cp().process({"W"})
287  .AddSyst(cb, "CMS_htt_WShape_etau_btag_high_$ERA", "shape", SystMap<channel, bin_id>::init({"et"}, {14}, 1.000));
288 
289  src.cp().process({"ggH_SM125", "qqH_SM125", "VH_SM125"})
290  .AddSyst(cb, "CMS_htt_SM125_mu", "lnN", SystMap<>::init(1.300));
291 
292  src.cp()
293  .AddSyst(cb, "pdf_qqbar", "lnN", SystMap<process>::init
294  ({"qqH_SM125"}, 1.036)
295  ({"VH_SM125"}, 1.040));
296 
297  src.cp().process({"ggH_SM125"})
298  .AddSyst(cb, "pdf_gg", "lnN", SystMap<>::init(1.097));
299 
300  src.cp()
301  .AddSyst(cb, "UEPS", "lnN", SystMap<bin_id, process>::init
302  ({10, 11, 12}, {"ggH_SM125"}, 1.013)
303  ({10, 11, 12}, {"qqH_SM125", "VH_SM125"}, 1.050)
304  ({13, 14}, {"ggH_SM125"}, 0.946)
305  ({13, 14}, {"qqH_SM125", "VH_SM125"}, 1.007));
306 
307  src.cp().process({"ggH_SM125"})
308  .AddSyst(cb, "QCDscale_ggH", "lnN", SystMap<channel, bin_id>::init
309  ({"mt", "et"}, {10, 11, 12}, 1.080)
310  ({"mt"}, {13, 14}, 1.105)
311  ({"et"}, {13, 14}, 1.125));
312 
313  src.cp().process({"qqH_SM125"})
314  .AddSyst(cb, "QCDscale_qqH", "lnN", SystMap<bin_id>::init
315  ({10, 11, 12}, 1.022)
316  ({13, 14}, 1.015));
317 
318  src.cp().process({"VH_SM125"})
319  .AddSyst(cb, "QCDscale_VH", "lnN", SystMap<bin_id>::init
320  ({10, 11, 12}, 1.010)
321  ({13, 14}, 1.040));
322 
323  src.cp().process({ggH})
324  .AddSyst(cb, "CMS_htt_higgsPtReweight_$ERA", "shape", SystMap<>::init(1));
325 
326  src.cp().process({ggH})
327  .AddSyst(cb, "CMS_htt_higgsPtReweight_scale_$ERA", "shape", SystMap<channel, bin_id>::init
328  ({"mt", "et"}, {10, 11, 12, 13}, 1)
329  ({"et"}, {14}, 1));
330 
331  src.cp().process({"ggH_SM125"})
332  .AddSyst(cb, "CMS_htt_higgsPtReweightSM_$ERA", "shape", SystMap<>::init(1.000));
333 }
334 
335 void AddMSSMUpdateSystematics_et_mt(CombineHarvester & cb) {
336 
337  CombineHarvester src = cb.cp();
338  AddMSSMUpdateSystematics_et_mt(cb, src);
339 
340 }
341 
342 void AddMSSMUpdateSystematics_em(CombineHarvester & cb, CombineHarvester src) {
343  //CombineHarvester src = cb.cp();
344  src.channel({"em"});
345 
346  std::vector<std::string> ggH = {"ggH", "ggH_Htautau", "ggA_Atautau", "ggh_htautau"};
347  std::vector<std::string> bbH = {"bbH", "bbH_Htautau", "bbA_Atautau", "bbh_htautau"};
348 
349  auto signal = Set2Vec(src.cp().signals().SetFromProcs(
350  std::mem_fn(&Process::process)));
351 
352  src.cp().process(JoinStr({signal,{"EWK", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
353  .AddSyst(cb, "lumi_$ERA", "lnN", SystMap<era>::init({"8TeV"}, 1.026));
354 
355  src.cp().process(JoinStr({signal,{"Ztt", "ttbar", "EWK", "Fakes", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
356  .AddSyst(cb, "CMS_eff_e_$ERA", "lnN", SystMap<>::init(1.02));
357 
358  src.cp().process(JoinStr({signal,{"Ztt", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
359  .AddSyst(cb, "CMS_scale_e_$ERA", "shape", SystMap<>::init(1));
360 
361  src.cp().process({"EWK"})
362  .AddSyst(cb, "CMS_shift1_emu_nobtag_$ERA_EWK_fine_binning", "shape", SystMap<bin_id>::init({8}, 1.000));
363 
364  src.cp().process({"EWK"})
365  .AddSyst(cb, "CMS_shift1_emu_btag_$ERA_EWK_fine_binning", "shape", SystMap<bin_id>::init({9}, 1.000));
366 
367  src.cp().process({"EWK"})
368  .AddSyst(cb, "CMS_shift2_emu_nobtag_$ERA_EWK_fine_binning", "shape", SystMap<bin_id>::init({8}, 1.000));
369  src.cp().process({"EWK"})
370  .AddSyst(cb, "CMS_shift2_emu_btag_$ERA_EWK_fine_binning", "shape", SystMap<bin_id>::init({9}, 1.000));
371 
372  src.cp().process(JoinStr({signal,{"Ztt", "ttbar", "EWK", "Fakes", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
373  .AddSyst(cb, "CMS_eff_m_$ERA", "lnN", SystMap<>::init(1.02));
374 
375  src.cp()
376  .AddSyst(cb, "CMS_scale_j_$ERA", "lnN", SystMap<bin_id, process>::init
377  ({8}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.01)
378  ({8}, {"ttbar", "EWK"}, 0.990)
379  ({9}, JoinStr({ggH, {"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.04)
380  ({9}, {bbH}, 1.01)
381  ({9}, {"ttbar"}, 0.930)
382  ({9}, {"EWK"}, 0.970));
383 
384  src.cp()
385  .AddSyst(cb, "CMS_htt_scale_met_$ERA", "lnN", SystMap<bin_id, process>::init
386  ({8, 9}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
387  ({8}, {"ttbar"}, 0.990)
388  ({9}, {"ttbar"}, 1.010));
389 
390  src.cp()
391  .AddSyst(cb, "CMS_eff_b_$ERA", "lnN", SystMap<bin_id, process>::init
392  ({8}, JoinStr({signal,{"EWK", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
393  ({8}, {"ttbar"}, 0.950)
394  ({9}, JoinStr({ggH, {"ggH_SM125"}}), 1.01)
395  ({9}, {bbH}, 1.05)
396  ({9}, {"ttbar"}, 1.020)
397  ({9}, {"EWK"}, 1.030));
398 
399  src.cp()
400  .AddSyst(cb, "CMS_fake_b_$ERA", "lnN", SystMap<bin_id, process>::init
401  ({8}, JoinStr({signal,{"ttbar", "EWK", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
402  ({9}, JoinStr({ggH, {"ggH_SM125"}}), 1.05)
403  ({9}, {"EWK"}, 1.030));
404 
405  src.cp().process({"Ztt"})
406  .AddSyst(cb, "CMS_htt_zttNorm_$ERA", "lnN", SystMap<>::init(1.030));
407 
408  src.cp().process({"Ztt"})
409  .AddSyst(cb, "CMS_htt_extrap_ztt_$CHANNEL_btag_$ERA", "lnN", SystMap<bin_id>::init({9}, 1.010));
410 
411  src.cp().process({"ttbar"})
412  .AddSyst(cb, "CMS_htt_ttbarNorm_$ERA", "lnN", SystMap<>::init(1.100));
413 
414  src.cp().process({"EWK"})
415  .AddSyst(cb, "CMS_htt_DiBosonNorm_$ERA", "lnN", SystMap<>::init(1.150));
416 
417  src.cp().process({"Fakes"})
418  .AddSyst(cb, "CMS_htt_fakes_$CHANNEL_$ERA", "lnN", SystMap<>::init(1.300));
419 
420  src.cp().process({"Fakes"})
421  .AddSyst(cb, "CMS_htt_fakes_$CHANNEL_btag_$ERA", "lnN", SystMap<bin_id>::init({9}, 1.090));
422 
423  src.cp().process({"ttbar"})
424  .AddSyst(cb, "CMS_htt_ttbar_emb_$ERA", "lnN", SystMap<bin_id>::init({9}, 1.020));
425 
426  src.cp().process({"Fakes"})
427  .AddSyst(cb, "CMS_htt_FakeShape_$CHANNEL_nobtag_$ERA", "shape", SystMap<bin_id>::init({8}, 1.000));
428  src.cp().process({"Fakes"})
429  .AddSyst(cb, "CMS_htt_FakeShape_$CHANNEL_btag_$ERA", "shape", SystMap<bin_id>::init({9}, 1.000));
430 
431  src.cp().process({"ttbar"})
432  .AddSyst(cb, "CMS_htt_TTbarShape_$CHANNEL_nobtag_$ERA", "shape", SystMap<bin_id>::init({8}, 1.000));
433  src.cp().process({"ttbar"})
434  .AddSyst(cb, "CMS_htt_TTbarShape_$CHANNEL_btag_$ERA", "shape", SystMap<bin_id>::init({9}, 1.000));
435 
436  src.cp().process({"ggH_SM125", "qqH_SM125", "VH_SM125"})
437  .AddSyst(cb, "CMS_htt_SM125_mu", "lnN", SystMap<>::init(1.300));
438 
439  src.cp()
440  .AddSyst(cb, "pdf_qqbar", "lnN", SystMap<bin_id, process>::init
441  ({8, 9}, {"qqH_SM125"}, 1.036)
442  ({8}, {"VH_SM125"}, 1.010)
443  ({9}, {"VH_SM125"}, 1.020));
444 
445  src.cp().process({"ggH_SM125"})
446  .AddSyst(cb, "pdf_gg", "lnN", SystMap<>::init(1.097));
447 
448  src.cp().process({"ggH_SM125"})
449  .AddSyst(cb, "QCDscale_ggH", "lnN", SystMap<bin_id>::init
450  ({8}, 1.080)
451  ({9}, 1.105));
452 
453  src.cp().process({"qqH_SM125"})
454  .AddSyst(cb, "QCDscale_qqH", "lnN", SystMap<bin_id>::init
455  ({8}, 1.034)
456  ({9}, 1.008));
457 
458  src.cp().process({"VH_SM125"})
459  .AddSyst(cb, "QCDscale_VH", "lnN", SystMap<>::init(1.040));
460 
461  src.cp()
462  .AddSyst(cb, "UEPS", "lnN", SystMap<bin_id, process>::init
463  ({8}, {"ggH_SM125"}, 1.035)
464  ({9}, {"ggH_SM125"}, 0.984)
465  ({8}, {"qqH_SM125", "VH_SM125"}, 1.089));
466 }
467 
468 void AddMSSMUpdateSystematics_em(CombineHarvester & cb) {
469 
470  CombineHarvester src = cb.cp();
471  AddMSSMUpdateSystematics_em(cb, src);
472 
473 }
474 
475 void AddMSSMUpdateSystematics_mm(CombineHarvester & cb, CombineHarvester src) {
476  //CombineHarvester src = cb.cp();
477  src.channel({"mm"});
478 
479  std::vector<std::string> ggH = {"ggH", "ggH_Htautau", "ggA_Atautau", "ggh_htautau"};
480  std::vector<std::string> bbH = {"bbH", "bbH_Htautau", "bbA_Atautau", "bbh_htautau"};
481 
482  auto signal = Set2Vec(src.cp().signals().SetFromProcs(
483  std::mem_fn(&Process::process)));
484 
485  src.cp()
486  .AddSyst(cb, "lumi_$ERA", "lnN", SystMap<bin_id, process, era>::init
487  ({8, 9}, JoinStr({signal,{"Dibosons", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), {"8TeV"}, 1.026)
488  ({8}, {"WJets"}, {"8TeV"}, 1.026));
489 
490  src.cp()
491  .AddSyst(cb, "CMS_eff_m_$ERA", "lnN", SystMap<bin_id, process>::init
492  ({8, 9}, JoinStr({signal,{"ZTT", "TTJ", "Dibosons", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.06)
493  ({8}, {"WJets"}, 1.060));
494 
495  src.cp()
496  .AddSyst(cb, "CMS_scale_j_$ERA", "lnN", SystMap<bin_id, process>::init
497  ({8}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.02)
498  ({8}, {"TTJ", "WJets", "Dibosons"}, 0.990)
499  ({9}, JoinStr({ggH,{"ggH_SM125"}}), 0.98)
500  ({9}, {bbH}, 1.01)
501  ({9}, {"TTJ"}, 0.920)
502  ({9}, {"Dibosons"}, 0.950));
503 
504  src.cp()
505  .AddSyst(cb, "CMS_eff_b_$ERA", "lnN", SystMap<bin_id, process>::init
506  ({8}, JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
507  ({8}, {"TTJ"}, 0.940)
508  ({8}, {"WJets", "Dibosons"}, 0.970)
509  ({9}, JoinStr({ggH, {"ggH_SM125"}}), 1.01)
510  ({9}, {bbH}, 1.05)
511  ({9}, {"TTJ"}, 1.020)
512  ({9}, {"Dibosons"}, 1.030));
513 
514  src.cp()
515  .AddSyst(cb, "CMS_fake_b_$ERA", "lnN", SystMap<bin_id, process>::init
516  ({8}, JoinStr({signal,{"TqTJ", "WJets", "Dibosons", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
517  ({9}, JoinStr({ggH, {"ggH_SM125"}}), 1.06)
518  ({9}, {"Dibosons"}, 1.050));
519 
520  src.cp().process({"ZTT"})
521  .AddSyst(cb, "CMS_htt_zttNorm_$ERA", "lnN", SystMap<>::init(1.050));
522 
523  src.cp().process({"TTJ"})
524  .AddSyst(cb, "CMS_htt_ttbarNorm_$ERA", "lnN", SystMap<>::init(1.100));
525 
526  src.cp().process({"Dibosons"})
527  .AddSyst(cb, "CMS_htt_$CHANNEL_DiBosonNorm_$ERA", "lnN", SystMap<>::init(1.300));
528 
529  src.cp().process({"ZMM"})
530  .AddSyst(cb, "CMS_htt_$CHANNEL_zmmNorm_$ERA", "lnN", SystMap<>::init(1.050));
531 
532  src.cp().process({"QCD"})
533  .AddSyst(cb, "CMS_htt_$CHANNEL_QCDNorm_$BIN_$ERA", "lnN", SystMap<bin_id>::init
534  ({8}, 1.100)
535  ({9}, 1.250));
536 
537  src.cp().process({"WJets"})
538  .AddSyst(cb, "CMS_htt_$CHANNEL_WJetsNorm_nobtag_$ERA", "lnN", SystMap<bin_id>::init({8}, 1.220));
539 
540  src.cp().process({"ZMM"})
541  .AddSyst(cb, "CMS_htt_$CHANNEL_res_met_$ERA", "shape", SystMap<bin_id>::init({8}, 1.000));
542 
543  src.cp().process({"ZTT"})
544  .AddSyst(cb, "CMS_htt_$CHANNEL_ztt_extrap_btag_$ERA", "lnN", SystMap<bin_id>::init({9}, 1.050));
545 
546  src.cp().process({"ZMM"})
547  .AddSyst(cb, "CMS_htt_$CHANNEL_zmm_extrap_btag_$ERA", "lnN", SystMap<bin_id>::init({9}, 1.050));
548 
549  src.cp().process({"ggH_SM125", "qqH_SM125", "VH_SM125"})
550  .AddSyst(cb, "CMS_htt_SM125_mu", "lnN", SystMap<>::init(1.300));
551 
552  src.cp().process({"qqH_SM125"})
553  .AddSyst(cb, "pdf_qqbar", "lnN", SystMap<>::init(1.036));
554 
555  src.cp().process({"ggH_SM125"})
556  .AddSyst(cb, "pdf_gg", "lnN", SystMap<>::init(1.097));
557 
558  src.cp().process({"VH_SM125"})
559  .AddSyst(cb, "pdf_vh", "lnN", SystMap<>::init(1.010));
560 
561  src.cp().process({"ggH_SM125"})
562  .AddSyst(cb, "QCDscale_ggH", "lnN", SystMap<bin_id>::init
563  ({8}, 1.103)
564  ({9}, 1.109));
565 
566  src.cp().process({"qqH_SM125"})
567  .AddSyst(cb, "QCDscale_qqH", "lnN", SystMap<bin_id>::init
568  ({8}, 1.034)
569  ({9}, 1.008));
570 
571  src.cp().process({"VH_SM125"})
572  .AddSyst(cb, "QCDscale_VH", "lnN", SystMap<>::init(1.040));
573 
574  src.cp()
575  .AddSyst(cb, "UEPS", "lnN", SystMap<bin_id, process>::init
576  ({8}, {"ggH_SM125"}, 1.035)
577  ({8}, {"qqH_SM125", "VH_SM125"}, 1.089)
578  ({9}, {"ggH_SM125"}, 0.984)
579  ({9}, {"qqH_SM125", "VH_SM125"}, 1.000));
580 }
581 
582 void AddMSSMUpdateSystematics_mm(CombineHarvester & cb) {
583 
584  CombineHarvester src = cb.cp();
585  AddMSSMUpdateSystematics_mm(cb, src);
586 
587 }
588 
589 void AddMSSMUpdateSystematics_tt(CombineHarvester & cb, CombineHarvester src) {
590  //CombineHarvester src = cb.cp();
591  src.channel({"tt"});
592 
593  std::vector<std::string> ggH = {"ggH", "ggH_Htautau", "ggA_Atautau", "ggh_htautau"};
594  std::vector<std::string> bbH = {"bbH", "bbH_Htautau", "bbA_Atautau", "bbh_htautau"};
595 
596  auto signal = Set2Vec(src.cp().signals().SetFromProcs(
597  std::mem_fn(&Process::process)));
598 
599  src.cp().process(JoinStr({signal,{"ggH_SM125", "qqH_SM125", "VH_SM125"}}))
600  .AddSyst(cb, "lumi_$ERA", "lnN", SystMap<era>::init({"8TeV"}, 1.026));
601 
602  src.cp().process(JoinStr({signal,{"ZTT", "TT", "VV", "ggH_SM125", "qqH_SM125", "VH_SM125"}}))
603  .AddSyst(cb, "CMS_eff_t_tautau_$ERA", "lnN", SystMap<>::init(1.19));
604 
605  src.cp().process(JoinStr({signal}))
606  .AddSyst(cb, "CMS_eff_t_mssmHigh_tautau_$ERA", "shape", SystMap<>::init(1));
607 
608  src.cp().process(JoinStr({signal,{"ZTT"}}))
609  .AddSyst(cb, "CMS_scale_t_tautau_$ERA", "shape", SystMap<>::init(1));
610 
611  src.cp().process({"QCD"})
612  .AddSyst(cb, "CMS_shift1_tauTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({10}, 1.000));
613  src.cp().process({"QCD"})
614  .AddSyst(cb, "CMS_shift1_tauTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({11}, 1.000));
615  src.cp().process({"QCD"})
616  .AddSyst(cb, "CMS_shift1_tauTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({12}, 1.000));
617  src.cp().process({"QCD"})
618  .AddSyst(cb, "CMS_shift1_tauTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({13}, 1.000));
619  src.cp().process({"QCD"})
620  .AddSyst(cb, "CMS_shift1_tauTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({14}, 1.000));
621 
622  src.cp().process({"QCD"})
623  .AddSyst(cb, "CMS_shift2_tauTau_nobtag_low_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({10}, 1.000));
624  src.cp().process({"QCD"})
625  .AddSyst(cb, "CMS_shift2_tauTau_nobtag_medium_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({11}, 1.000));
626  src.cp().process({"QCD"})
627  .AddSyst(cb, "CMS_shift2_tauTau_nobtag_high_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({12}, 1.000));
628  src.cp().process({"QCD"})
629  .AddSyst(cb, "CMS_shift2_tauTau_btag_low_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({13}, 1.000));
630  src.cp().process({"QCD"})
631  .AddSyst(cb, "CMS_shift2_tauTau_btag_high_$ERA_QCD_fine_binning", "shape", SystMap<bin_id>::init({14}, 1.000));
632 
633  src.cp()
634  .AddSyst(cb, "CMS_scale_j_$ERA", "lnN", SystMap<process>::init
635  ({ggH}, 1.05)
636  ({bbH}, 1.06)
637  ({"TT"}, 1.010)
638  ({"VV"}, 1.030)
639  ({"ggH_SM125"}, 1.050));
640 
641  src.cp()
642  .AddSyst(cb, "CMS_eff_b_$ERA", "lnN", SystMap<bin_id, process>::init
643  ({10, 11, 12}, JoinStr({signal,{"ZJ", "TT", "VV", "ZL", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
644  ({13, 14}, JoinStr({signal,{"TT", "VV", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.09));
645 
646  src.cp()
647  .AddSyst(cb, "CMS_fake_b_$ERA", "lnN", SystMap<bin_id, process>::init
648  ({10, 11, 12}, JoinStr({signal,{"ZJ", "TT", "VV", "ZL", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 0.98)
649  ({13, 14}, JoinStr({signal,{"TT", "VV", "ggH_SM125", "qqH_SM125", "VH_SM125"}}), 1.02));
650 
651  src.cp().process({"W"})
652  .AddSyst(cb, "CMS_htt_WNorm_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
653  ({10, 11, 12, 13}, 1.300)
654  ({14}, 1.500));
655 
656  src.cp().process({"QCD"})
657  .AddSyst(cb, "CMS_htt_QCDSyst_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
658  ({10, 11, 12}, 1.100)
659  ({13, 14}, 1.200));
660 
661  src.cp().process({"QCD"})
662  .AddSyst(cb, "CMS_htt_QCDfrShape_tautau_$ERA", "shape", SystMap<>::init(1.000));
663 
664  src.cp().process({"ZTT", "ZJ", "ZL"})
665  .AddSyst(cb, "CMS_htt_zttNorm_$ERA", "lnN", SystMap<>::init(1.033));
666 
667  src.cp().process({"ZTT"})
668  .AddSyst(cb, "CMS_htt_extrap_ztt_$CHANNEL_$BIN_$ERA", "lnN", SystMap<bin_id>::init
669  ({10, 11, 12}, 1.050)
670  ({13}, 1.090)
671  ({14}, 1.180));
672 
673  src.cp().process({"TT"})
674  .AddSyst(cb, "CMS_htt_ttbarNorm_$ERA", "lnN", SystMap<>::init(1.100));
675 
676  src.cp().process({"TT"})
677  .AddSyst(cb, "CMS_htt_ttbarNorm_$CHANNEL_$bin_$ERA", "lnN", SystMap<>::init(1.050));
678 
679  src.cp().process({"TT"})
680  .AddSyst(cb, "CMS_htt_ttbarPtReweight_$ERA", "shape", SystMap<>::init(1.000));
681 
682  src.cp().process({"VV"})
683  .AddSyst(cb, "CMS_htt_DiBosonNorm_$ERA", "lnN", SystMap<>::init(1.150));
684 
685  src.cp().process({"VV"})
686  .AddSyst(cb, "CMS_htt_DiBosonNorm_$CHANNEL_$BIN_$ERA", "lnN", SystMap<>::init(1.150));
687 
688  src.cp().process({"ZJ"})
689  .AddSyst(cb, "CMS_htt_ZJetFakeTau_$CHANNEL_$ERA", "lnN", SystMap<>::init(1.200));
690 
691  src.cp().process({"ZL"})
692  .AddSyst(cb, "CMS_htt_ZLeptonFakeTau_$CHANNEL_$ERA", "lnN", SystMap<>::init(1.200));
693 
694  src.cp().process({"W"})
695  .AddSyst(cb, "CMS_htt_WShape_tautau_$ERA", "shape", SystMap<>::init(1.000));
696 
697  src.cp().process({"ggH_SM125", "qqH_SM125", "VH_SM125"})
698  .AddSyst(cb, "CMS_htt_SM125_mu", "lnN", SystMap<>::init(1.300));
699 
700  src.cp().process({"ggH_SM125"})
701  .AddSyst(cb, "pdf_gg", "lnN", SystMap<>::init(1.097));
702 
703  src.cp()
704  .AddSyst(cb, "pdf_qqbar", "lnN", SystMap<process>::init
705  ({"qqH_SM125"}, 1.036)
706  ({"VH_SM125"}, 1.020));
707 
708  src.cp().process({"ggH_SM125"})
709  .AddSyst(cb, "QCDscale_ggH", "lnN", SystMap<>::init(1.205));
710 
711  src.cp().process({"qqH_SM125"})
712  .AddSyst(cb, "QCDscale_qqH", "lnN", SystMap<>::init(1.012));
713 
714  src.cp().process({"VH_SM125"})
715  .AddSyst(cb, "QCDscale_VH", "lnN", SystMap<>::init(1.040));
716 
717  src.cp()
718  .AddSyst(cb, "UEPS", "lnN", SystMap<process>::init
719  ({"ggH_SM125"}, 0.975)
720  ({"qqH_SM125", "VH_SM125"}, 1.025));
721 
722  src.cp().process({ggH})
723  .AddSyst(cb, "CMS_htt_higgsPtReweight_$ERA", "shape", SystMap<>::init(1));
724 
725  src.cp().process({"ggH_SM125"})
726  .AddSyst(cb, "CMS_htt_higgsPtReweightSM_$ERA", "shape", SystMap<>::init(1.000));
727 }
728 
729 void AddMSSMUpdateSystematics_tt(CombineHarvester & cb) {
730 
731  CombineHarvester src = cb.cp();
732  AddMSSMUpdateSystematics_tt(cb, src);
733 
734 }
735 
736 }
ch::CombineHarvester::signals
CombineHarvester & signals()
Definition: CombineHarvester_Filters.cc:146
ch::CombineHarvester::process
CombineHarvester & process(std::vector< std::string > const &vec, bool cond=true)
Definition: CombineHarvester_Filters.cc:35
ch::CombineHarvester::SetFromProcs
std::set< R > SetFromProcs(T func)
Fill an std::set using only the Process entries.
Definition: CombineHarvester.h:615
ch::CombineHarvester::AddSyst
void AddSyst(CombineHarvester &target, std::string const &name, std::string const &type, Map const &valmap)
Definition: CombineHarvester.h:682
ch::CombineHarvester::cp
CombineHarvester cp()
Creates and returns a shallow copy of the CombineHarvester instance.
Definition: CombineHarvester.cc:220
ch::CombineHarvester::channel
CombineHarvester & channel(std::vector< std::string > const &vec, bool cond=true)
Definition: CombineHarvester_Filters.cc:83
ch::Object::process
virtual std::string const & process() const
Definition: Object.h:20
ch::syst::bin_id
Definition: Systematics.h:41
ch
Definition: Algorithm.h:10
ch::AddMSSMUpdateSystematics_et_mt
void AddMSSMUpdateSystematics_et_mt(CombineHarvester &cb, CombineHarvester src)
Definition: HttSystematics_MSSMUpdate.cc:18
ch::AddMSSMUpdateSystematics_mm
void AddMSSMUpdateSystematics_mm(CombineHarvester &cb, CombineHarvester src)
Definition: HttSystematics_MSSMUpdate.cc:475
ch::Set2Vec
std::vector< T > Set2Vec(std::set< T > const &in)
Definition: Utilities.h:101
ch::JoinStr
std::vector< std::string > JoinStr(std::vector< std::vector< std::string >> const &in)
Definition: Utilities.cc:153
ch::AddMSSMUpdateSystematics_tt
void AddMSSMUpdateSystematics_tt(CombineHarvester &cb, CombineHarvester src)
Definition: HttSystematics_MSSMUpdate.cc:589
ch::AddMSSMUpdateSystematics_em
void AddMSSMUpdateSystematics_em(CombineHarvester &cb, CombineHarvester src)
Definition: HttSystematics_MSSMUpdate.cc:342