Skip to content

Commit 8292cfe

Browse files
committed
up templates/v19
1 parent 32e4df2 commit 8292cfe

File tree

3 files changed

+68
-16
lines changed

3 files changed

+68
-16
lines changed

binder/combine.ipynb

Lines changed: 65 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -160,38 +160,38 @@
160160
},
161161
{
162162
"cell_type": "code",
163-
"execution_count": 36,
163+
"execution_count": 79,
164164
"id": "c6707244",
165165
"metadata": {},
166166
"outputs": [],
167167
"source": [
168168
"# get lumi\n",
169169
"import pickle as pkl\n",
170170
"\n",
171-
"version = \"vtest\"\n",
171+
"version = \"v19\"\n",
172172
"year = \"Run2\"\n",
173173
"with open(f\"../combine/templates/{version}/hists_templates_{year}.pkl\", \"rb\") as f:\n",
174174
" h = pkl.load(f)"
175175
]
176176
},
177177
{
178178
"cell_type": "code",
179-
"execution_count": 37,
179+
"execution_count": 76,
180180
"id": "ac045a42",
181181
"metadata": {},
182182
"outputs": [
183183
{
184184
"data": {
185185
"text/plain": [
186186
"Hist(\n",
187-
" StrCategory(['ggFpt200to300', 'ggFpt300to450', 'ggFpt450toInf', 'ggF', 'Fake'], growth=True, name='Sample'),\n",
187+
" StrCategory(['Data', 'EWKvjets', 'WH', 'TTbar', 'SingleTop', 'DYJets', 'Diboson', 'ttH', 'WJetsLNu', 'ZH', 'mjj1000toInf', 'VBF', 'ggFpt200to300', 'ggFpt300to450', 'ggFpt450toInf', 'ggF', 'Fake'], growth=True, name='Sample'),\n",
188188
" StrCategory(['nominal', 'trigger_ele_SF_up', 'trigger_ele_SF_down', 'EW_up', 'EW_down', 'weight_pdf_acceptance_up', 'weight_pdf_acceptance_down', 'weight_qcd_scale_up', 'weight_qcd_scale_down', 'weight_btagSFlightCorrelated_up', 'weight_btagSFlightCorrelated_down', 'weight_btagSFbcCorrelated_up', 'weight_btagSFbcCorrelated_down', 'weight_btagSFlight_2016_up', 'weight_btagSFlight_2016_down', 'weight_btagSFbc_2016_up', 'weight_btagSFbc_2016_down', 'weight_btagSFlight_2016APV_up', 'weight_btagSFlight_2016APV_down', 'weight_btagSFbc_2016APV_up', 'weight_btagSFbc_2016APV_down', 'weight_btagSFlight_2017_up', 'weight_btagSFlight_2017_down', 'weight_btagSFbc_2017_up', 'weight_btagSFbc_2017_down', 'weight_btagSFlight_2018_up', 'weight_btagSFlight_2018_down', 'weight_btagSFbc_2018_up', 'weight_btagSFbc_2018_down', 'weight_pileup_id_up', 'weight_pileup_id_down', 'weight_PSFSR_up', 'weight_PSFSR_down', 'weight_PSISR_up', 'weight_PSISR_down', 'weight_d1kappa_EW_up', 'weight_d1kappa_EW_down', 'weight_Z_d2kappa_EW_up', 'weight_Z_d2kappa_EW_down', 'weight_Z_d3kappa_EW_up', 'weight_Z_d3kappa_EW_down', 'weight_ele_id_up', 'weight_ele_id_down', 'weight_ele_reco_up', 'weight_ele_reco_down', 'weight_mu_isolation_up', 'weight_mu_isolation_down', 'weight_mu_id_stat_up', 'weight_mu_id_stat_down', 'weight_mu_id_syst_up', 'weight_mu_id_syst_down', 'weight_mu_trigger_iso_up', 'weight_mu_trigger_iso_down', 'weight_mu_trigger_noniso_up', 'weight_mu_trigger_noniso_down', 'weight_TopPtReweight_up', 'weight_TopPtReweight_down', 'weight_pileup_2016_up', 'weight_pileup_2016_down', 'weight_PSFSR_wjets_2016_up', 'weight_PSFSR_wjets_2016_down', 'weight_L1Prefiring_2016_up', 'weight_L1Prefiring_2016_down', 'weight_pileup_2016APV_up', 'weight_pileup_2016APV_down', 'weight_PSFSR_wjets_2016APV_up', 'weight_PSFSR_wjets_2016APV_down', 'weight_L1Prefiring_2016APV_up', 'weight_L1Prefiring_2016APV_down', 'weight_pileup_2017_up', 'weight_pileup_2017_down', 'weight_PSFSR_wjets_2017_up', 'weight_PSFSR_wjets_2017_down', 'weight_L1Prefiring_2017_up', 'weight_L1Prefiring_2017_down', 'weight_pileup_2018_up', 'weight_pileup_2018_down', 'weight_PSFSR_wjets_2018_up', 'weight_PSFSR_wjets_2018_down', 'JES_FlavorQCD_up', 'JES_FlavorQCD_down', 'JES_RelativeBal_up', 'JES_RelativeBal_down', 'JES_HF_up', 'JES_HF_down', 'JES_BBEC1_up', 'JES_BBEC1_down', 'JES_EC2_up', 'JES_EC2_down', 'JES_Absolute_up', 'JES_Absolute_down', 'UES_up', 'UES_down', 'JES_BBEC1_2016_up', 'JES_BBEC1_2016_down', 'JES_RelativeSample_2016_up', 'JES_RelativeSample_2016_down', 'JES_EC2_2016_up', 'JES_EC2_2016_down', 'JES_HF_2016_up', 'JES_HF_2016_down', 'JES_Absolute_2016_up', 'JES_Absolute_2016_down', 'JER_2016_up', 'JER_2016_down', 'JMR_2016_up', 'JMR_2016_down', 'JMS_2016_up', 'JMS_2016_down', 'JES_BBEC1_2016APV_up', 'JES_BBEC1_2016APV_down', 'JES_RelativeSample_2016APV_up', 'JES_RelativeSample_2016APV_down', 'JES_EC2_2016APV_up', 'JES_EC2_2016APV_down', 'JES_HF_2016APV_up', 'JES_HF_2016APV_down', 'JES_Absolute_2016APV_up', 'JES_Absolute_2016APV_down', 'JER_2016APV_up', 'JER_2016APV_down', 'JMR_2016APV_up', 'JMR_2016APV_down', 'JMS_2016APV_up', 'JMS_2016APV_down', 'JES_BBEC1_2017_up', 'JES_BBEC1_2017_down', 'JES_RelativeSample_2017_up', 'JES_RelativeSample_2017_down', 'JES_EC2_2017_up', 'JES_EC2_2017_down', 'JES_HF_2017_up', 'JES_HF_2017_down', 'JES_Absolute_2017_up', 'JES_Absolute_2017_down', 'JER_2017_up', 'JER_2017_down', 'JMR_2017_up', 'JMR_2017_down', 'JMS_2017_up', 'JMS_2017_down', 'JES_BBEC1_2018_up', 'JES_BBEC1_2018_down', 'JES_RelativeSample_2018_up', 'JES_RelativeSample_2018_down', 'JES_EC2_2018_up', 'JES_EC2_2018_down', 'JES_HF_2018_up', 'JES_HF_2018_down', 'JES_Absolute_2018_up', 'JES_Absolute_2018_down', 'JER_2018_up', 'JER_2018_down', 'JMR_2018_up', 'JMR_2018_down', 'JMS_2018_up', 'JMS_2018_down', 'FR_stat_Up', 'FR_stat_Down', 'EWK_SF_Up', 'EWK_SF_Down'], growth=True, name='Systematic'),\n",
189-
" StrCategory(['VBF'], growth=True, name='Region'),\n",
189+
" StrCategory(['VBF', 'ggFpt250to350', 'ggFpt350to500', 'ggFpt500toInf', 'WJetsCR', 'TopCR'], growth=True, name='Region'),\n",
190190
" Variable([75, 95, 115, 135, 155, 175, 195, 215, 235], name='mass_observable', label='Higgs reconstructed mass [GeV]'),\n",
191-
" storage=Weight()) # Sum: WeightedSum(value=479.993, variance=6.57083) (WeightedSum(value=487.132, variance=6.74091) with flow)"
191+
" storage=Weight()) # Sum: WeightedSum(value=2.30656e+06, variance=1.60248e+06) (WeightedSum(value=2.59663e+06, variance=1.80048e+06) with flow)"
192192
]
193193
},
194-
"execution_count": 37,
194+
"execution_count": 76,
195195
"metadata": {},
196196
"output_type": "execute_result"
197197
}
@@ -202,26 +202,78 @@
202202
},
203203
{
204204
"cell_type": "code",
205-
"execution_count": 39,
206-
"id": "0fd202da",
205+
"execution_count": 123,
206+
"id": "42745855",
207+
"metadata": {},
208+
"outputs": [
209+
{
210+
"data": {
211+
"text/plain": [
212+
"array([ 43.32175251, 83.59237614, 142.94556479, 179.7033345 ,\n",
213+
" 168.2463245 , 142.88461987, 111.70961338, 68.35236569])"
214+
]
215+
},
216+
"execution_count": 123,
217+
"metadata": {},
218+
"output_type": "execute_result"
219+
}
220+
],
221+
"source": [
222+
"h17[{\"Sample\": \"TTbar\", \"Region\": \"WJetsCR\", \"Systematic\": \"JER_2018_down\"}].values()"
223+
]
224+
},
225+
{
226+
"cell_type": "code",
227+
"execution_count": 124,
228+
"id": "a6bb4809",
229+
"metadata": {},
230+
"outputs": [
231+
{
232+
"data": {
233+
"text/plain": [
234+
"array([ 86.64350501, 167.18475228, 285.89112958, 359.40666901,\n",
235+
" 336.49264899, 285.76923973, 223.41922675, 136.70473138])"
236+
]
237+
},
238+
"execution_count": 124,
239+
"metadata": {},
240+
"output_type": "execute_result"
241+
}
242+
],
243+
"source": [
244+
"h19[{\"Sample\": \"TTbar\", \"Region\": \"WJetsCR\", \"Systematic\": \"JER_2018_down\"}].values()"
245+
]
246+
},
247+
{
248+
"cell_type": "code",
249+
"execution_count": 99,
250+
"id": "76b1e85e",
207251
"metadata": {},
208252
"outputs": [
209253
{
210254
"data": {
211255
"text/plain": [
212-
"array([0.06911826, 0.001 , 0.001 , 0.00278477, 0.001 ,\n",
213-
" 0.0490055 , 0.001 , 0.001 ])"
256+
"array([ 86.09231681, 167.56327396, 286.76986975, 360.7549742 ,\n",
257+
" 336.1691917 , 284.3532072 , 220.93714581, 139.37088821])"
214258
]
215259
},
216-
"execution_count": 39,
260+
"execution_count": 99,
217261
"metadata": {},
218262
"output_type": "execute_result"
219263
}
220264
],
221265
"source": [
222-
"h[{\"Sample\": \"Fake\", \"Region\": \"VBF\", \"Systematic\": \"nominal\"}].values()"
266+
"h20[{\"Sample\": \"TTbar\", \"Region\": \"WJetsCR\", \"Systematic\": \"JER_2018_down\"}].values()"
223267
]
224268
},
269+
{
270+
"cell_type": "code",
271+
"execution_count": null,
272+
"id": "9a3af5b2",
273+
"metadata": {},
274+
"outputs": [],
275+
"source": []
276+
},
225277
{
226278
"cell_type": "code",
227279
"execution_count": 40,

combine/config_make_templates.yaml

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -20,9 +20,9 @@ regions_sel:
2020
"ggFpt350to500": (n_bjets_T==0) & (THWW>0.93) & ( (mjj<1000) | (deta<3.5) | (NumOtherJets<2) ) & (rec_higgs_pt>350) & (rec_higgs_pt<500) & (VH_fj_VScore<0.9)
2121
"ggFpt500toInf": (n_bjets_T==0) & (THWW>0.93) & ( (mjj<1000) | (deta<3.5) | (NumOtherJets<2) ) & (rec_higgs_pt>500) & (VH_fj_VScore<0.9)
2222
"WJetsCR": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>250) & (NumFatjets>=2)
23-
"WJetsCRpt250to350": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>250) & (NumFatjets>=2) & (rec_higgs_pt<350)
24-
"WJetsCRpt350to500": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>350) & (NumFatjets>=2) & (rec_higgs_pt<500)
25-
"WJetsCRpt500toInf": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>500) & (NumFatjets>=2)
23+
# "WJetsCRpt250to350": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>250) & (NumFatjets>=2) & (rec_higgs_pt<350)
24+
# "WJetsCRpt350to500": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>350) & (NumFatjets>=2) & (rec_higgs_pt<500)
25+
# "WJetsCRpt500toInf": (n_bjets_T==0) & (THWW<0.905) & (rec_higgs_pt>500) & (NumFatjets>=2)
2626
"TopCR": (n_bjets_T>0) & (rec_higgs_pt>250) & (THWW>0.90)
2727

2828
samples_dir:
-1.25 MB
Binary file not shown.

0 commit comments

Comments
 (0)