-
Notifications
You must be signed in to change notification settings - Fork 7
/
Plotting_Janos.h
879 lines (826 loc) · 88.7 KB
/
Plotting_Janos.h
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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
// Implementation of Plotting class defined in PlottingBase.h
#ifndef PLOTTING_H
#define PLOTTING_H
#include "include/Analysis.h"
//_______________________________________________________
// Calculate variables
//_______________________________________________________
// Define Histo options: Filling/Postfixes
void
Plotting::define_additional_histo_settings(const Weighting& w, EventSelections& evt_sel,
const unsigned int& syst_nSyst, const unsigned int& syst_index)
{
}
//_______________________________________________________
// Define Histograms here
void
Plotting::define_additional_histos(const Weighting& w, const unsigned int& syst_nSyst, const unsigned int& syst_index) {
//__________________________________
// Define Smarthistos
//__________________________________
// Usage of AddHistos:
// sh.AddHistos(type, { .fill=fill parameter(s), .pfs={list of posfix names}, .cuts={list of cuts},.draw="ROOT Draw() options", .opt="SmartHistos options", .ranges={plot ranges and overridden legend location}});
// 1st param: string, the object/event type defined by the first parameter in AddHistoType()
// 2nd param (.fill): string, either the name of a FillParam using the first parameter of AddNeFillParams() for a 1D histo, or two fillparams separated by "_vs_" for a 2D one
// 3rd param (.pfs): vector of strings for each Postfix name defined by the first parameter in AddNewPostfix(), supports up to 5 Postfix
// N.B: for 1D plots, the first postfix is used to draw different plots on the same canvas, while creating multiple plots with the rest
// 4th param (.cuts): vector of cuts before filling (I do not encourage the usage, because usually postfixes work nicer by attaching a string to the plot name)
// 5th param (.draw): The option to use when using the TH1D/TH2D.Draw() method
// 6th param (.opt): SmartHistos specific options, some common examples:
// "Stack"+"N" : Make a stack histo for all postfix, except the first N (can be data, signals, etc.)
// "Log" : Log Y or Z scale
// "LogX" : Log X scale
// "TwoCol"+"XY" : Split legend to two columns, and specify how many rows should each hav maximum
// "AddInt" : Add integral of the histogram to the legend
// "AddRatio" : Create a ratio of either 1st pf / 2nd pf, or 1st and pf/stacked histo
// There are many more, please peek into source code for complete list or see other examples
// 7th .ranges: List of axis ranges in order: xlow,xhigh, ylow,yhigh, then zlow,zhigh for 2D plot
// OR the location of the top left edge of the legend (x1,y2), the legend size is auto adjusted
// (by setting x2 and y1 depending on the number of entries)
// Systematics:
// The systematics of the plot are hidden in an additional dimension added to the histogram
// eg. a TH1D and it's systematic variations are stored in a TH2D
// Each extra dimensional bin contain a variation of the underlying distribution
// To enable this, a special dummy axis needs to be added which is called "Counts"
// The string "Counts_vs_" is added to the plot name
// and the "syst " is attached to the histogram type so the code knows to fill this
// additional axis
bool doSyst = (syst_nSyst>0);
bool systematics = 1;
std::string s = "";
std::string c = "";
if (systematics) {
s = "syst ";
c = "Counts_vs_";
}
std::vector<std::string> standard_plots;
if (doSyst) standard_plots = {"HT", "MET", "HTMET", "MRFine", "MR", "R2Fine", "R2", "MRR2"};
else standard_plots = {"HT", "HTFine", "METFine", "MET", "HTMET", "MRFine", "MR", "R2Fine", "R2", "MRR2"};
// -------------------------------------------------------------------------
// Trigger
sh.SetHistoWeights({ [] { return 1; } });
// Hadronic triggers
//for (auto plot : {"HadronicMeasurements", "HadronicMeasurementsNoR2", "NonIsoLepMeasurements"}) {
// std::string cut = "Baseline";
// if (std::string(plot)=="HadronicMeasurementsNoR2") cut = "BaselineNoR2";
// for (auto trigger : std::vector<std::string>({
// "HLTEff_PFHT1050",
// "HLTEff_PFMET120_PFMHT120",
// "HLTEff_PFMETNoMu120_PFMHTNoMu120",
// "HLTEff_PFHT500_PFMET100_PFMHT100",
// "HLTEff_PFHT700_PFMET85_PFMHT85",
// "HLTEff_PFHT800_PFMET75_PFMHT75",
// "HLTEff_HT_MET_or_HT",
// "HLTEff_AllMET_or_HT_MET_or_HT"
// })) {
// for (auto std_plot : standard_plots)
// sh.AddHistos("evt", { .fill=trigger+"_vs_"+std_plot, .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_Bin", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HT2DBins", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_MET2DBins", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HTMET", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HT1MET", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HT2MET", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HT3MET", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// sh.AddHistos("evt", { .fill=trigger+"_vs_HTMETPrev", .pfs={ plot,"Year",cut}, .cuts={}, .draw="PE1", .opt=o_1or2d_d+"TwoCol22", .ranges={0,0, 0,1, 0.5,0.35} });
// }
//}
// ----------------------------------------------------------------------------------------------
// New preselection segions
//-----------------------------------------------------------------------------------------------
if (doSyst) standard_plots = {"HT", "MET", "MRFine", "MR", "R2Fine", "R2"};
else standard_plots = {"HT", "HTFine", "METFine", "MET", "MRFine", "MR", "R2Fine", "R2"};
//standard_plots.push_back("HTFine");
//standard_plots.push_back("HTMET");
for (auto region : {Region::Pre_Had}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
// Signal discriminators
sh.AddHistos("evt", { .fill="MRFine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MRFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="METFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_R2Fine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Boosted jet Index
sh.AddHistos("hadw", { .fill="JetAK8Index_vs_HadWIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadW"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("hadz", { .fill="JetAK8Index_vs_HadZIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadZ"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("hadtop", { .fill="JetAK8Index_vs_HadTopIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadTop"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("leptop", { .fill="JetAK8Index_vs_LepTopIndex", .pfs={"Signals_Background","Year",cut,"MatchGenLepTop"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadw", { .fill="GenMatchedAK8JetIndex_vs_NGenHadW", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadz", { .fill="GenMatchedAK8JetIndex_vs_NGenHadZ", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadtop", { .fill="GenMatchedAK8JetIndex_vs_NGenHadTop", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen leptop", { .fill="GenMatchedAK8JetIndex_vs_NGenLepTop", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// MGlunio/MSquark plots
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBins", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBins", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBins", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
// Stack plots
std::string opt = o_stk_d_S;
for (auto nobj_bin : {"NJet45",
"NW1", "NW2", "NW3",
//"NTop1", "NTop2", "NTop3", "NTop4", "NTop5",
//"NWDeepMD1", "NWDeepMD2", "NWDeepMD3",
//"NZDeepMD1", "NZDeepMD2", "NZDeepMD3",
"NTopDeepMD1", "NTopDeepMD2", "NTopDeepMD3", "NTopDeepMD4",
"NWDeep1", "NWDeep2", "NWDeep3",
"NZDeep1", "NZDeep2", "NZDeep3",
"NVDeep1", "NVDeep2", "NVDeep3",
"NTopDeep1", "NTopDeep2", "NTopDeep3", "NTopDeep4"}) {
std::string signal = "StackPlotTopSignal";
opt = o_stk_d_T;
std::string bveto = "Nisob";
if (TString(nobj_bin).Contains("NZ")||TString(nobj_bin).Contains("NV")) {
signal = "StackPlotVSignal";
opt = o_stk_d_V;
bveto = "Nisob_NJet24";
}
for (auto std_plot : standard_plots) {
sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,bveto}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"MTBoost800"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
sh.AddHistos(s+"evt", { .fill=c+"MRR2", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
if (nobj_bin != std::string("NJet45")) {
sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,bveto}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"MTBoost800"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
sh.AddHistos( "evt", { .fill="DeltaPhi", .pfs={signal,"Year",cut+"_ExcldPhi",nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// MTBoost
sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,bveto}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
// Exclusive W combinations
//for (std::string nobj_bin : {"Nb", "Nisob", "2b", "NHadTop", "NTopDeepMD1", "NTopDeepMD2", "NTopDeep1", "NTopDeep2"}) {
// for (std::string nw_bin : {"NHadW"}) {
// opt = o_stk_d_T;
// for (auto std_plot : standard_plots)
// sh.AddHistos("evt", { .fill=std_plot, .pfs={"StackPlotTopSignal","Year",cut,nobj_bin,nw_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos(s+"evt", { .fill=c+"MRR2", .pfs={"StackPlotTopSignal","Year",cut,nobj_bin,nw_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="NJet", .pfs={"StackPlotTopSignal","Year",cut,nobj_bin,nw_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="MTBoost", .pfs={"StackPlotTopSignal","Year",cut,nobj_bin,nw_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="DeltaPhi", .pfs={"StackPlotTopSignal","Year",cut+"_ExcldPhi",nobj_bin,nw_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// }
//}
}
for (auto region : {Region::Pre_NonIso}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
// Leptons
sh.AddHistos("evt", { .fill="NLepNoIso_vs_NGenLepFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="NEleNoIso_vs_NGenEleFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="NMuNoIso_vs_NGenMuFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="NLepNonIso_vs_NGenLepFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="NEleNonIso_vs_NGenEleFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="NMuNonIso_vs_NGenMuFromTop", .pfs={"Signals_Background","Year",cut+"_Excl1LepMTdPhiJet"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Boosted jet Index
sh.AddHistos("hadw", { .fill="JetAK8Index_vs_HadWIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadW"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("hadz", { .fill="JetAK8Index_vs_HadZIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadZ"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("hadtop", { .fill="JetAK8Index_vs_HadTopIndex", .pfs={"Signals_Background","Year",cut,"MatchGenHadTop"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("leptop", { .fill="JetAK8Index_vs_LepTopIndex", .pfs={"Signals_Background","Year",cut,"MatchedGenLeptonMother"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadw", { .fill="GenMatchedAK8JetIndex_vs_NGenHadW", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadz", { .fill="GenMatchedAK8JetIndex_vs_NGenHadZ", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen hadtop", { .fill="GenMatchedAK8JetIndex_vs_NGenHadTop", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("gen leptop", { .fill="GenMatchedAK8JetIndex_vs_NGenLepTop", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Signal discriminators (2D)
sh.AddHistos("evt", { .fill="MRFine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MRFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="METFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MTNonIso", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MTNonIso", .pfs={"Signals_Background","Year",cut+"_ExclMT"},.cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_R2Fine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// MGlunio/MSquark plots
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBinsLep", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBinsLep", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRR2", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="HTFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="METFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="MRFine", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
sh.AddHistos("evt", { .fill="RazorBinsLep", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut}, .cuts={},.draw=d,.opt=o_norm_s,.ranges={}});
// Stack plots
std::string opt = o_stk_d_S;
// Boosted leptonic jets
sh.AddHistos("AK8", { .fill="LepJetPt", .pfs={"StackPlotSignal","Year",cut,"LepJetNoPt"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetPt", .pfs={"StackPlotSignal","Year",cut,"LepJetHighMass"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetEta", .pfs={"StackPlotSignal","Year",cut,"LepJet"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetPhi", .pfs={"StackPlotSignal","Year",cut,"LepJet"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetMass", .pfs={"StackPlotSignal","Year",cut,"LepJetCand"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetSubJetDeepB", .pfs={"StackPlotSignal","Year",cut,"LepTopNoSubJetB"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetTau21", .pfs={"StackPlotSignal","Year",cut,"LepJet"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetLSF", .pfs={"StackPlotSignal","Year",cut,"LepJet"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepJetNSubJet", .pfs={"StackPlotSignal","Year",cut,"LepJet"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepJetLSFNoIso", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepJetNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepJetDRmin", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepJetNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepJetPtrel", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepJetNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
// Boosted leptonic tops
sh.AddHistos("AK8", { .fill="LepTopPt", .pfs={"StackPlotSignal","Year",cut,"LepTopNoPt"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopPt", .pfs={"StackPlotSignal","Year",cut,"LepTopHighMass"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopEta", .pfs={"StackPlotSignal","Year",cut,"LepTop"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopPhi", .pfs={"StackPlotSignal","Year",cut,"LepTop"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopMass", .pfs={"StackPlotSignal","Year",cut,"LepTopNoMass"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopSubJetDeepB", .pfs={"StackPlotSignal","Year",cut,"LepTopNoSubJetB"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopTau21", .pfs={"StackPlotSignal","Year",cut,"LepTop"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopLSF", .pfs={"StackPlotSignal","Year",cut,"LepTop"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos("AK8", { .fill="LepTopNSubJet", .pfs={"StackPlotSignal","Year",cut,"LepTop"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepTopLSFNoIso", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepTopNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepJetDRmin", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepTopNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
//sh.AddHistos("AK8", { .fill="LepJetPtrel", .pfs={"StackPlotSignal","Year",cut+"_Excl1LepMTdPhiJet","LepTopNoIso"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
for (auto nobj_bin : {"NJet45",
//"NTop1", "NTop2", "NTop3", "NTop4", "NTop5",
//"NWDeepMD1", "NWDeepMD2", "NWDeepMD3",
//"NZDeepMD1", "NZDeepMD2", "NZDeepMD3",
"NTopDeepMD1", "NTopDeepMD2", "NTopDeepMD3", "NTopDeepMD4",
"NLepTop", "NLepTopHighMass",
"NTopDeep1", "NTopDeep2", "NTopDeep3", "NTopDeep4",
"NW1", "NW2", "NW3",
"NWDeep1", "NWDeep2", "NWDeep3",
"NLepJetCand", "NLepJet", "NLepJetHighMass",
"NZDeep1", "NZDeep2", "NZDeep3",
"NVDeep1", "NVDeep2", "NVDeep3"}) {
std::string signal = "StackPlotTopSignal";
opt = o_stk_d_T;
if (TString(nobj_bin).Contains("NZ")||TString(nobj_bin).Contains("NV")||TString(nobj_bin).Contains("LepJet")) {
signal = "StackPlotVSignal";
opt = o_stk_d_V;
}
for (auto std_plot : standard_plots) {
sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"Nisob"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill=std_plot, .pfs={signal,"Year",cut,nobj_bin,"MTBoost800"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
sh.AddHistos(s+"evt", { .fill=c+"MRR2", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
if (nobj_bin != std::string("NJet45")) {
sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"Nisob"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos("evt", { .fill="NJet", .pfs={signal,"Year",cut,nobj_bin,"MTBoost800"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
sh.AddHistos( "evt", { .fill="DeltaPhi", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos( "evt", { .fill="DeltaPhi", .pfs={signal,"Year",cut+"_ExclMT",nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// MTBoost
sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,"Nb"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,"Nisob"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
//sh.AddHistos( "evt", { .fill="MTBoost", .pfs={signal,"Year",cut,nobj_bin,"2b"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
// Exclusive Leptonic+Hadronic combinations
//for (std::string nhad_bin : {"Nb", "Nisob", "2b", "NHadTop", "NHadW", "NTopDeepMD1", "NTopDeepMD2", "NTopDeep1", "NTopDeep2"}) {
// for (std::string nlep_bin : {"NLepJetCand", "NLepJet", "NLepJetHighMass", "NLepTop", "NLepTopHighMass"}) {
// opt = o_stk_d_T;
// for (auto std_plot : standard_plots) {
// sh.AddHistos("evt", { .fill=std_plot, .pfs={"StackPlotTopSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill=std_plot, .pfs={"StackPlotTopSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// }
// sh.AddHistos(s+evt", { .fill=c+"MRR2", .pfs={"StackPlotTopSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos(s+evt", { .fill=c+"MRR2", .pfs={"StackPlotTopSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="NJet", .pfs={"StackPlotTopSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="NJet", .pfs={"StackPlotTopSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="MTBoost", .pfs={"StackPlotTopSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
// sh.AddHistos( "evt", { .fill="DeltaPhi", .pfs={"StackPlotTopSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// sh.AddHistos( "evt", { .fill="DeltaPhi", .pfs={"StackPlotTopSignal","Year",cut+"_ExclMT",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// }
//}
//for (std::string nhad_bin : {"Nb", "Nisob", "2b", "NHadW", "NZDeep1", "NZDeep2", "NZDeep3", "NVDeep1", "NVDeep2", "NVDeep3"}) {
// for (std::string nlep_bin : {"NLepJetCand", "NLepJet", "NLepJetHighMass"}) {
// opt = o_stk_d_V;
// for (auto std_plot : standard_plots) {
// sh.AddHistos("evt", { .fill=std_plot, .pfs={"StackPlotVSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill=std_plot, .pfs={"StackPlotVSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// }
// sh.AddHistos(s+"evt", { .fill=c+"MRR2", .pfs={"StackPlotVSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos(s+"evt", { .fill=c+"MRR2", .pfs={"StackPlotVSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="NJet", .pfs={"StackPlotVSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="NJet", .pfs={"StackPlotVSignal","Year",cut+"_Excl1LepMTdPhiJet",nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
// sh.AddHistos("evt", { .fill="MTBoost", .pfs={"StackPlotVSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
// sh.AddHistos("evt", { .fill="DeltaPhi", .pfs={"StackPlotVSignal","Year",cut, nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// sh.AddHistos("evt", { .fill="DeltaPhi", .pfs={"StackPlotVSignal","Year",cut+"_ExclMT", nhad_bin,nlep_bin}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
// }
//}
// Gen Truth plots
for (const auto& cut2 : { cut, cut+"_Excl1LepMTdPhiJet" }) {
// gen lep/top pt
sh.AddHistos("gen lep", { .fill="GenLepPtBins", .pfs={"Signals_Background","Year",cut2,"GenLepMother","GenLeptonFlavour"}, .cuts={}, .draw="HIST", .opt=o_1or2d_s+"Norm", .ranges={0,0, 0,1, 0.5,0.53} });
//sh.AddHistos("gen leptop", { .fill="GenLepTopMatchedGenLepPt", .pfs={"Signals_Background","Year",cut2,"GenTopLeptonFlavour"}, .cuts={}, .draw="HIST", .opt=o_1or2d_s+"Norm", .ranges={0,0, 0,1, 0.5,0.53} }); // Equivalent with above
sh.AddHistos("gen lep", { .fill="GenLepPtBins", .pfs={"MGluinoPoints", "MLSPPoints","GluinoSignalScans", cut2,"GenLepMother"}, .cuts={}, .draw="HIST", .opt=o_1or2d_s+"Norm", .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen lep", { .fill="GenLepPtBins", .pfs={"MSquarkPoints", "MLSPPoints","SquarkSignalScans", cut2,"GenLepMother"}, .cuts={}, .draw="HIST", .opt=o_1or2d_s+"Norm", .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen lep", { .fill="GenLepPtBins", .pfs={"MEWKPoints","MLSPPoints","EWKSignalScans",cut2,"GenLepMother"}, .cuts={}, .draw="HIST", .opt=o_1or2d_s+"Norm", .ranges={0,0, 0,1, 0.5,0.53} });
//sh.AddHistos("gen leptop", { .fill="GenLepTopMatchedGenLepPt_vs_GenLepTopPt", .pfs={"Signals_Background","Year",cut2,"GenTopLeptonFlavour"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
}
}
// Lepton efficiencies
//sh.SetHistoWeights({ [&w] { return isData ? 0 : w.sf_weight[Region::Pre]; } });
sh.SetHistoWeights({ [this] { return v.isData ? 0 : 1; } });
sh.AddHistos("gen ele", { .fill="EleCBVetoEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen ele", { .fill="EleCBVetoNoIsoEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen mu", { .fill="MuonCBLooseEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen mu", { .fill="MuonCBLooseNoIsoEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen mu", { .fill="MuonCBMediumEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
sh.AddHistos("gen mu", { .fill="MuonCBMediumNoIsoEfficiency_vs_GenLepPtBins", .pfs={"GenLepMother","Pre","Year","Signals_Background"}, .cuts={}, .draw="PE1",.opt=o_1or2d_s, .ranges={0,0, 0,1, 0.5,0.53} });
// ROC Curves
sh.SetHistoWeights({ [&w] { return w.sf_weight[Region::Pre]; } });
// ---------------------- Leptons ---------------------
for (std::string benchmark : {"Bkg"}) {
// ---------------------- Veto ----------------------
for (std::string mother : {"W", "HardProcess"}) {
// Electrons
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+benchmark, .pfs={"elept5_cut", "eleeta_veto_cut",
"eleip_medium", "eleid_mva_loose", "eleminiiso_loose"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+benchmark, .pfs={"elept5_cut", "eleeta_veto_cut", "eleip_medium_cut",
"eleid_mva_miniiso", "eleid_mva_2diso", "eleid_mva_iso", "eleid_cut_iso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+benchmark, .pfs={"elept5_cut", "eleeta_veto_cut", "eleip_medium_cut", "eleid_mva_loose_cut",
"eleminiiso", "elepfiso", "ele2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
// Muons
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+benchmark, .pfs={"mupt5_cut", "mueta_cut",
"muip_loose", "muid_soft", "muminiiso_loose"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+benchmark, .pfs={"mupt5_cut", "mueta_cut", "muip_loose_cut", "muminiiso_loose_cut",
"muid_mva", "muid_soft_comp", "muid_medium_comp", "muid_tight_comp"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+benchmark, .pfs={"mupt5_cut", "mueta_cut", "muip_loose_cut", "muid_soft_cut",
"muminiiso", "mupfiso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0.5,1, 0.9,1}});
}
}
// --------------------- Isolated ---------------------
for (const auto& bm : all_benchmarks) {
// Electrons
for (std::string mother : {"W", "HardProcess"}) {
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+bm.first, .pfs={"eleeta_cut",
"elept10", "eleip_medium", "eleid_mva_medium", "eleminiiso_tight"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+bm.first, .pfs={"elept10_cut", "eleeta_cut", "eleip_medium_cut",
"eleid_mva_miniiso", "eleid_mva_2diso", "eleid_mva_iso", "eleid_cut_iso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
sh.AddHistos("ele", { .fill="EleFrom"+mother+"_"+bm.first, .pfs={"elept10_cut", "eleeta_cut", "eleip_medium_cut", "eleid_mva_medium_cut",
"eleminiiso", "elepfiso", "ele2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
// Muons
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+bm.first, .pfs={"mueta_cut",
"mupt10", "muip_medium", "muid_mediumprompt", "muminiiso_medium"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+bm.first, .pfs={"mupt10_cut", "mueta_cut", "muip_medium_cut", "muminiiso_medium_cut",
"muid_mva", "muid_soft_comp", "muid_medium_comp", "muid_tight_comp"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
sh.AddHistos("mu", { .fill="MuFrom"+mother+"_"+bm.first, .pfs={"mupt10_cut", "mueta_cut", "muip_medium_cut", "muid_mediumprompt_cut",
"muminiiso", "mupfiso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.96,1.0}});
}
}
for (const auto& bm : top_benchmarks) {
// ------------------- Non-Isolated -------------------
for (std::string mother : {"Top"}) {
// Electrons and their neutrinos
for (std::string object : {"Ele","EleNu"}) {
sh.AddHistos("ele", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"eleeta_cut",
"elept30", "eleip_medium", "eleid_mva_loose", "ele2diso15"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
sh.AddHistos("ele", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"elept30_cut", "eleeta_cut", "eleip_medium_cut",
"eleid_mva_miniiso", "eleid_mva_2diso", "eleid_mva_iso", "eleid_cut_iso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
sh.AddHistos("ele", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"elept30_cut", "eleeta_cut", "eleip_medium_cut", "eleid_mva_loose_cut",
"eleminiiso", "elepfiso", "ele2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
}
// Muons and their neutrinos
for (std::string object : {"Mu","MuNu"}) {
sh.AddHistos("mu", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"mueta_cut",
"mupt30", "muip_tight", "muid_mva_loose", "mu2diso15"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
sh.AddHistos("mu", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"mupt30_cut", "mueta_cut", "muip_tight_cut", "mu2diso15_cut",
"muid_mva", "muid_soft_comp", "muid_medium_comp", "muid_tight_comp"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
sh.AddHistos("mu", { .fill=object+"From"+mother+"_"+bm.first, .pfs={"mupt30_cut", "mueta_cut", "muip_tight_cut", "muid_mva_loose_cut",
"muminiiso", "mupfiso", "mu2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
}
}
// ------------------- Leptonic top -------------------
// Top with W decaying to electron and the neutrino
for (std::string object : {"EleLepTop"}) {//,"EleNuFromLepTop"}) {
sh.AddHistos("ele", { .fill=object+"_"+bm.first, .pfs={"eleeta_cut", "ele_ak8eta_cut",
"ele_ak8pt200", "ele_msoftdrop0", "ele_sjbtag_loose", //"deepTag_top",
"ele_neutrinodr", //"ele_neutrinojetdr",
"elept30", "eleip_medium", "eleid_mva_loose", "ele2diso15"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
sh.AddHistos("ele", { .fill=object+"_"+bm.first, .pfs={"elept30_cut", "eleeta_cut", "eleip_medium_cut",
"ele_ak8eta_cut", "ele_ak8pt200_cut", "ele_sjbtag_loose_cut",
"eleid_mva_miniiso", "eleid_mva_2diso", "eleid_mva_iso", "eleid_cut_iso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
sh.AddHistos("ele", { .fill=object+"_"+bm.first, .pfs={"elept30_cut", "eleeta_cut", "eleip_medium_cut", "eleid_mva_loose_cut",
"ele_ak8eta_cut", "ele_ak8pt200_cut", "ele_sjbtag_loose_cut",
"eleminiiso", "elepfiso", "ele2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
}
// Top with W decaying to muon and the neutrino
for (std::string object : {"MuLepTop"}) {//,"MuNuFromLepTop"}) {
sh.AddHistos("mu", { .fill=object+"_"+bm.first, .pfs={"mueta_cut", "mu_ak8eta_cut",
"mu_ak8pt200", "mu_msoftdrop0", "mu_sjbtag_loose", //"deepTag_top",
"mu_neutrinodr", //"mu_neutrinojetdr",
"mupt30", "muip_tight", "muid_mva_loose", "mu2diso15"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
sh.AddHistos("mu", { .fill=object+"_"+bm.first, .pfs={"mupt30_cut", "mueta_cut", "muip_tight_cut", "mu2diso15_cut",
"mu_ak8eta_cut", "mu_ak8pt200_cut", "mu_sjbtag_loose_cut",
"muid_mva", "muid_soft_comp", "muid_medium_comp", "muid_tight_comp"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
sh.AddHistos("mu", { .fill=object+"_"+bm.first, .pfs={"mupt30_cut", "mueta_cut", "muip_tight_cut", "muid_mva_loose_cut",
"mu_ak8eta_cut", "mu_ak8pt200_cut", "mu_sjbtag_loose_cut",
"muminiiso", "mupfiso", "mu2diso"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.8, 0.97,1.0}});
}
// ------------------- Hadronic top -------------------
sh.AddHistos("AK8", { .fill="HadTop_"+bm.first, .pfs={"ak8pt500_cut", "tau32_softdrop_sjbtag", "deepTagMD_top", "deepTagMD_top_minsd", "deepTag_top"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.9,1.0}});
sh.AddHistos("AK8", { .fill="HadTop_"+bm.first, .pfs={"ak8pt500", "msoftdrop105", "deepTagMD_top_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.9,1.0}});
sh.AddHistos("AK8", { .fill="HadTop_"+bm.first, .pfs={"ak8pt500", "msoftdrop_min","deepTag_top_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.9,1.0}});
}
for (const auto& bm : all_benchmarks) {
// -------------------- Hadronic W --------------------
sh.AddHistos("AK8", { .fill="HadW_"+bm.first, .pfs={"ak8pt200_cut", "tau21_softdrop",
"deepTagMD_W", "deepTagMD_W_minsd", "deepTagMD_W_sd", "deepTag_W"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.6, 0.6,1.0}});
//sh.AddHistos("AK8", { .fill="HadW_"+bm.first, .pfs={"ak8pt200", "msoftdrop65", "msoftdrop_max", "deepTagMD_W_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
//sh.AddHistos("AK8", { .fill="HadW_"+bm.first, .pfs={"ak8pt200", "msoftdrop_min", "msoftdrop_max", "deepTag_W_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
//sh.AddHistos("AK8", { .fill="HadW_"+bm.first, .pfs={"ak8pt200", "msoftdrop65", "msoftdrop_max", "deepTagMD_W_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
sh.AddHistos("AK8", { .fill="HadW_"+bm.first, .pfs={"ak8pt200", "msoftdrop_min", "msoftdrop_max", "deepTag_W_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
}
for (const auto& bm : WZH_benchmarks) {
// -------------------- Hadronic Z --------------------
sh.AddHistos("AK8", { .fill="HadZ_"+bm.first, .pfs={"ak8pt200_cut", "tau21_softdrop",
"deepTagMD_Z", "deepTagMD_Z_minsd", "deepTagMD_Z_sd", "deepTag_Z"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
//sh.AddHistos("AK8", { .fill="HadZ_"+bm.first, .pfs={"ak8pt200", "msoftdrop80", "msoftdrop_max", "deepTagMD_Z_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
//sh.AddHistos("AK8", { .fill="HadZ_"+bm.first, .pfs={"ak8pt200", "msoftdrop_min", "msoftdrop_max", "deepTag_Z_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
sh.AddHistos("AK8", { .fill="HadZ_"+bm.first, .pfs={"ak8pt200", "msoftdrop80", "msoftdrop_max", "deepTagMD_Z_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
sh.AddHistos("AK8", { .fill="HadZ_"+bm.first, .pfs={"ak8pt200", "msoftdrop_min", "msoftdrop_max", "deepTag_Z_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0.5, 0.9,1.0}});
// -------------------- Hadronic H --------------------
sh.AddHistos("AK8", { .fill="HadH_"+bm.first, .pfs={"ak8pt300", "msoftdrop100", "msoftdrop_max", "deepTagMD_H_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
sh.AddHistos("AK8", { .fill="HadH_"+bm.first, .pfs={"ak8pt300", "deepTag_H_WP1"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
//sh.AddHistos("AK8", { .fill="HadH_"+bm.first, .pfs={"ak8pt300", "msoftdrop_min", "msoftdrop_max", "deepTag_H_WP2"}, .cuts={"P"},.draw="PLX",.opt="ROC", .ranges={0,0, 0.8,1.0}});
}
// ----------------------------------------------------------------------------------------------
// New signal segions
//-----------------------------------------------------------------------------------------------
// Fully hadronic Signal Regions
for (auto region : {Region::SR_Had_1htop, Region::SR_Had_2htop, Region::SR_Had_V_b_45j, Region::SR_Had_V_b_6j, Region::SR_Had_1V_0b_34j, Region::SR_Had_1V_0b_5j, Region::SR_Had_2V_0b_24j, Region::SR_Had_2V_0b_5j,
Region::SR_Had_H_b_45j, Region::SR_Had_H_b_6j, Region::SR_Had_H_b_6j, Region::SR_Had_HV_b_6j, Region::SR_Had_H_0b_34j, Region::SR_Had_H_0b_5j, Region::SR_Had_H_0b_34j, Region::SR_Had_H_0b_5j, Region::SR_Had_HV_0b_24j, Region::SR_Had_HV_0b_5j}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
std::vector<std::string> showdata = {"Blind"};
// Signal discriminators
sh.AddHistos("evt", { .fill="MRFine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MRFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="METFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_R2Fine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Stack plots
for (auto data : showdata ) {
std::string signal = "StackPlotSignal";
std::string opt = (data=="Blind") ? o_stk_s_S : o_stk_d_S;
if (region==Region::SR_Had_1htop||region==Region::SR_Had_2htop||region==Region::SR_Had_V_b_45j||region==Region::SR_Had_V_b_6j) {
opt = (data=="Blind") ? o_stk_s_T : o_stk_d_T;
signal = "StackPlotTopSignal";
} else if (region==Region::SR_Had_1V_0b_34j||region==Region::SR_Had_1V_0b_5j||region==Region::SR_Had_2V_0b_24j||region==Region::SR_Had_2V_0b_5j) {
opt = (data=="Blind") ? o_stk_s_V : o_stk_d_V;
signal = "StackPlotVSignal";
} else if (region==Region::SR_Had_H_b_45j||region==Region::SR_Had_H_b_6j||region==Region::SR_Had_H_b_6j||region==Region::SR_Had_HV_b_6j||
region==Region::SR_Had_H_0b_34j||region==Region::SR_Had_H_0b_5j||region==Region::SR_Had_H_0b_34j||region==Region::SR_Had_H_0b_5j||
region==Region::SR_Had_HV_0b_24j||region==Region::SR_Had_HV_0b_5j) {
opt = (data=="Blind") ? o_stk_s_H : o_stk_d_H;
signal = "StackPlotHSignal";
}
// Masses
sh.AddHistos( "hadtop", { .fill="HadTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos( "leptop", { .fill="LepTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos( "hadw", { .fill="HadWMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos( "hadz", { .fill="HadZMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
sh.AddHistos( "hadh", { .fill="HadHMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk5});
}
}
// Isolated lepton Signal Regions
for (auto region : {Region::SR_Lep_1htop, Region::SR_Lep_V_b, Region::SR_Lep_V_0b, Region::SR_Lep_H_b, Region::SR_Lep_H_0b}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
std::vector<std::string> showdata = {"Blind"};
// Signal discriminators
sh.AddHistos("evt", { .fill="MRFine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MRFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="METFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_R2Fine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Stack plots
for (auto data : showdata ) {
std::string signal = "StackPlotSignal";
std::string opt = (data=="Blind") ? o_stk_s_S : o_stk_d_S;
if (region==Region::SR_Lep_1htop) {
opt = (data=="Blind") ? o_stk_s_T : o_stk_d_T;
signal = "StackPlotTopSignal";
} else if (region==Region::SR_Lep_V_0b) {
opt = (data=="Blind") ? o_stk_s_V : o_stk_d_V;
signal = "StackPlotVSignal";
} else if (region==Region::SR_Lep_H_b||region==Region::SR_Lep_H_0b) {
opt = (data=="Blind") ? o_stk_s_H : o_stk_d_H;
signal = "StackPlotHSignal";
}
// Masses
sh.AddHistos( "hadtop", { .fill="HadTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "leptop", { .fill="LepTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadw", { .fill="HadWMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadz", { .fill="HadZMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadh", { .fill="HadHMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
}
}
// Non-isolated lepton Signal Regions
//for (auto region : {Region::SR_Leptop_0htop, Region::SR_Leptop_1htop, Region::SR_Lepjet_0V_24j, Region::SR_Lepjet_0V_5j, Region::SR_Lepjet_1V_24j, Region::SR_Lepjet_1V_5j}) {
for (auto region : {Region::SR_Lepjet_1htop, Region::SR_Lepjet_0V_24j, Region::SR_Lepjet_0V_5j, Region::SR_Lepjet_1V_24j, Region::SR_Lepjet_1V_5j}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
std::vector<std::string> showdata = {"Blind"};
// Signal discriminators
sh.AddHistos("evt", { .fill="MRFine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MRFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="METFine_vs_HTFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="R2Fine_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_METFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_R2Fine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("evt", { .fill="MTBoost_vs_MRFine", .pfs={"Signals_Background","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
// Stack plots
for (auto data : showdata ) {
std::string signal = "StackPlotSignal";
std::string opt = (data=="Blind") ? o_stk_s_S : o_stk_d_S;
// Masses
sh.AddHistos( "hadtop", { .fill="HadTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "leptop", { .fill="LepTopMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadw", { .fill="HadWMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadz", { .fill="HadZMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos( "hadh", { .fill="HadHMass", .pfs={signal,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
}
}
// -------------------------------------------------------------------------
// Selected AK4 Jets
for (const auto& regions : magic_enum::enum_entries<Region>()) {
Region region = regions.first;
std::string cut(regions.second);
if (TString(cut).BeginsWith("CR_")||TString(cut).BeginsWith("Val_")||TString(cut).BeginsWith("Pre_")) {
//if (region==Region::CR_Fake)
//sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_had_nor2; } });
//else
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string opt = o_stk_d;
sh.AddHistos("AK4", { .fill="JetPtBins", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhi", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetDeepB", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhotonDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhotonPtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhotonPtRatio", .pfs={Stack,"Year",cut,"JetPhotonDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhotonPtRatio", .pfs={Stack,"Year",cut,"JetPhotonDR0.4"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetPhotonDR_vs_JetPhotonPtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
sh.AddHistos("AK4", { .fill="JetEleDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetElePtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetElePtRatio", .pfs={Stack,"Year",cut,"JetEleDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetElePtRatio", .pfs={Stack,"Year",cut,"JetEleDR0.4"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetEleDR_vs_JetElePtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
sh.AddHistos("AK4", { .fill="JetMuonDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetMuonPtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetMuonPtRatio", .pfs={Stack,"Year",cut,"JetMuonDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetMuonPtRatio", .pfs={Stack,"Year",cut,"JetMuonDR0.4"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk9});
sh.AddHistos("AK4", { .fill="JetMuonDR_vs_JetMuonPtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
}
}
// -------------------------------------------------------------------------
// Leptons
// Leptons
for (auto region : {Region::Pre, Region::CR_QCD17_1Boost, Region::CR_Top17_1Boost, Region::CR_W17_1Boost, Region::CR_1LepInv, Region::CR_1LepInv_LepTrig, Region::CR_1PhoInv, Region::Val_Signal, Region::Val_QCD}) {
//for (auto region : {'S', Region::CR_QCD, Region::Val_Signal, Region::Val_QCD}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
std::string cut(magic_enum::enum_name(region));
if (region==Region::CR_2LepInv||region==Region::CR_1LepInv_LepTrig)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_lep; } });
if (region==Region::CR_1PhoInv&&v.year!=2016)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_pho; } });
if (region==Region::Pre) {
sh.AddHistos("evt", { .fill="NEleVeto", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoElePt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoEleEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NMuVeto", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NLepVeto", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NIsoTrk", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NTauVeto", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("tau veto", { .fill="VetoTauPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("tau veto", { .fill="VetoTauEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("ele tight noiso", { .fill="EleCleanJetDRmin", .pfs={"Signals_Background","Year",cut,"MatchGenEleFromTop"}, .cuts={},.draw="HIST",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("ele tight noiso", { .fill="EleCleanJetPtrel", .pfs={"Signals_Background","Year",cut,"MatchGenEleFromTop"}, .cuts={},.draw="HIST",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("ele tight noiso", { .fill="EleCleanJetPtrel_vs_EleCleanJetDRmin", .pfs={"Signals_Background","Year",cut,"MatchGenEleFromTop"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("mu tight noiso", { .fill="MuCleanJetDRmin", .pfs={"Signals_Background","Year",cut,"MatchGenMuFromTop"}, .cuts={},.draw="HIST",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("mu tight noiso", { .fill="MuCleanJetPtrel", .pfs={"Signals_Background","Year",cut,"MatchGenMuFromTop"}, .cuts={},.draw="HIST",.opt=o_1or2d_s+"Log",.ranges={}});
sh.AddHistos("mu tight noiso", { .fill="MuCleanJetPtrel_vs_MuCleanJetDRmin", .pfs={"Signals_Background","Year",cut,"MatchGenMuFromTop"}, .cuts={},.draw="COLZ",.opt=o_1or2d_s+"Log",.ranges={}});
} else if (region==Region::CR_Top17_1Boost||region==Region::CR_W17_1Boost||region==Region::CR_1LepInv||region==Region::CR_1LepInv_LepTrig||region==Region::Val_Signal_LostLep) {
sh.AddHistos("evt", { .fill="NEleVeto", .pfs={Stack,"Year",cut+"_Excl1LepMT"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoElePt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoEleEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NMuVeto", .pfs={Stack,"Year",cut+"_Excl1LepMT"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NLepVeto", .pfs={Stack,"Year",cut+"_Excl1LepMT"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
} else {
std::string opt = o_stk_d;
/*
sh.AddHistos("evt", { .fill="NEleVeto", .pfs={Stack,"Year",cut+"_Excl0Ele0Tau"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoElePt", .pfs={Stack,"Year",cut+"_Excl0Ele0Tau"}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("ele veto", { .fill="VetoEleEta", .pfs={Stack,"Year",cut+"_Excl0Ele0Tau"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NMuVeto", .pfs={Stack,"Year",cut+"_Excl0Mu0Tau"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuPt", .pfs={Stack,"Year",cut+"_Excl0Mu0Tau"}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("mu veto", { .fill="VetoMuEta", .pfs={Stack,"Year",cut+"_Excl0Mu0Tau"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NLepVeto", .pfs={Stack,"Year",cut+"_Excl0Ele0Mu0Tau"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NIsoTrk", .pfs={Stack,"Year",cut+"_Excl0Tau"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NTauVeto", .pfs={Stack,"Year",cut+"_Excl0Tau"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("tau veto", { .fill="VetoTauPt", .pfs={Stack,"Year",cut+"_Excl0Tau"}, .cuts={},.draw=d,.opt=o_stk_d+"logX",.ranges=r_Stk6});
sh.AddHistos("tau veto", { .fill="VetoTauEta", .pfs={Stack,"Year",cut+"_Excl0Tau"}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk6});
*/
}
}
// Selected Leptons
for (auto region : {Region::Pre, Region::CR_2LepInv}) {
//for (auto region : {Region::CR_2LepInv, Region::CR_Top, Region::CR_W}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
if (region==Region::CR_2LepInv)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_lep; } });
std::string cut(magic_enum::enum_name(region));
std::string cut2 = cut;
if (region==Region::CR_2LepInv) cut2 += "_Excl2LepR2";
//else cut2 += "_Excl1LepMT";
sh.AddHistos("evt", { .fill="ElePt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk7});
sh.AddHistos("evt", { .fill="EleEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk7});
sh.AddHistos("evt", { .fill="MuPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk7});
sh.AddHistos("evt", { .fill="MuEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk7});
sh.AddHistos("evt", { .fill="NLep", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NLep", .pfs={Stack,"run2",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NEle", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NMu", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NLepTight", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NLepNoIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NEleNoIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NMuNoIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NLepNonIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NLepNonIso", .pfs={Stack,"run2",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NEleNonIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NMuNonIso", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=o_stk_d,.ranges=r_Stk8});
}
// -------------------------------------------------------------------------
// bs
// Selected b-tags
for (auto region : {Region::Pre, Region::CR_Top17_1Boost, Region::CR_2LepInv, Region::CR_1PhoInv, Region::Val_Signal}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
if (region==Region::CR_2LepInv)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_lep; } });
if (region==Region::CR_1PhoInv&&v.year!=2016)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_pho; } });
std::string cut1(magic_enum::enum_name(region));
std::string cut2 = cut1;
//if (region==Region::Pre||region==Region::CR_2LepInv||region==Region::CR_1PhoInv)
// sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]*(*scale_factors['S'][2]); } });
//if (region=='S'||region==Region::Val_Signal||region==Region::CR_Top) cut2 += "_Excl1b";
std::string opt = o_stk_d;
sh.AddHistos("b", { .fill="BJetPtBins", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("b", { .fill="BJetPt", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("b", { .fill="BJetEta", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("b", { .fill="BJetPhi", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("b", { .fill="BJetDeepB", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("evt", { .fill="NBTag", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NBTag", .pfs={Stack,"run2",cut1}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
//if (region!=Region::Pre&®ion!=Region::CR_2LepInv&®ion!=Region::CR_1PhoInv)
// sh.AddHistos("evt", { .fill="NBTag", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
}
// Veto b-tags
for (auto region : {Region::Pre, Region::CR_QCD17_1Boost, Region::CR_W17_1Boost, Region::CR_1LepInv, Region::CR_1LepInv_LepTrig, Region::CR_2LepInv, Region::CR_1PhoInv, Region::Val_QCD}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
if (region==Region::CR_2LepInv||region==Region::CR_1LepInv_LepTrig)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_lep; } });
if (region==Region::CR_1PhoInv&&v.year!=2016)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_pho; } });
std::string cut1(magic_enum::enum_name(region));
//std::string cut2 = cut1;
//if (region==Region::CR_QCD||region==Region::Val_QCD||region==Region::CR_W||region==Region::CR_1LepInv) cut2 += "_Excl0b";
sh.AddHistos("evt", { .fill="NLooseBTag", .pfs={Stack,"Year",cut1}, .cuts={},.draw=d,.opt=O_stk_d+"NoComb5",.ranges={0,0, 1.01e-2,1e6, 0.5,0.86}});
//if (region!=Region::Pre&®ion!=Region::CR_2LepInv&®ion!=Region::CR_1PhoInv)
// sh.AddHistos("evt", { .fill="NLooseBTag", .pfs={Stack,"Year",cut2}, .cuts={},.draw=d,.opt=O_stk_d+"NoComb5",.ranges={0,0, 1.01e-2,1e6, 0.5,0.86}});
}
// -------------------------------------------------------------------------
// AK8/Ws Jets
for (auto region : {Region::Pre, Region::CR_Top17_1Boost,Region::CR_W17_1Boost,Region::CR_QCD17_1Boost, Region::CR_1LepInv, Region::CR_1LepInv_LepTrig, Region::CR_2LepInv, Region::CR_1PhoInv, Region::CR_Fake, Region::Val_Signal, Region::Val_QCD}) {
sh.SetHistoWeights({ [&w,region] { return w.sf_weight[region]; } });
if (region==Region::CR_2LepInv||region==Region::CR_1LepInv_LepTrig)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_lep; } });
if (region==Region::CR_1PhoInv&&v.year!=2016)
sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_pho; } });
//if (region==Region::CR_Fake||region==Region::CR_Fake_MET500)
// sh.SetHistoWeights({ [&w,region] { return w.w_nm1[region][9]*w.triggereff_had_nor2; } });
std::string cut1(magic_enum::enum_name(region));
//std::string cut2 = cut1;
//if (region==Region::Val_Signal||region==Region::CR_Top) cut2 += "_Excl1W";
//else if (region==Region::CR_QCD||region==Region::Val_QCD) cut2 += "_Excl1aW";
//else if (region==Region::CR_W) cut2 += "_Excl1mW";
//else if (region==Region::CR_1LepInv||region==Region::CR_2LepInv||region==Region::CR_1PhoInv) cut2 += "_Excl1mW";
std::vector<std::string> cuts;
cuts.push_back(cut1);
//if (cut2!=cut1) cuts.push_back(cut2);
for (auto cut : cuts) {
std::string opt = o_stk_d;
std::string Opt = O_stk_d;
// AK8 jets
sh.AddHistos("evt", { .fill="NJetAK8", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("evt", { .fill="NJetAK8", .pfs={Stack,"run2",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8PtBins", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Pt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Pt", .pfs={Stack,"run2",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Eta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Phi", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut,"Tau21Tagged"}, .cuts={},.draw=d,.opt=Opt,.ranges={0,0, 1.01e-2,1e8, 0.3,0.86, 0.21}});
sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut,"Tau32Tagged"}, .cuts={},.draw=d,.opt=Opt,.ranges={0,0, 1.01e-2,1e8, 0.3,0.86}});
sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut,"AK8_EB_EE"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
//sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut,"AK8_EB_EE","Tau21Tagged"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
//sh.AddHistos("AK8", { .fill="JetAK8Mass", .pfs={Stack,"Year",cut,"AK8_EB_EE","Tau32Tagged"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Tau21", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8Tau32", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="MaxAK8SubJetDeepB", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8PhotonDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8PhotonPtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8PhotonPtRatio", .pfs={Stack,"Year",cut,"JetAK8PhotonDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8PhotonPtRatio", .pfs={Stack,"Year",cut,"JetAK8PhotonDR0.8"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8PhotonDR_vs_JetAK8PhotonPtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
sh.AddHistos("AK8", { .fill="JetAK8EleDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8ElePtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8ElePtRatio", .pfs={Stack,"Year",cut,"JetAK8EleDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8ElePtRatio", .pfs={Stack,"Year",cut,"JetAK8EleDR0.8"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8EleDR_vs_JetAK8ElePtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
sh.AddHistos("AK8", { .fill="JetAK8MuonDR", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8MuonPtRatio", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8MuonPtRatio", .pfs={Stack,"Year",cut,"JetAK8MuonDR0.05"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8MuonPtRatio", .pfs={Stack,"Year",cut,"JetAK8MuonDR0.8"}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk6});
sh.AddHistos("AK8", { .fill="JetAK8MuonDR_vs_JetAK8MuonPtRatio", .pfs={"Data_MC","Year",cut}, .cuts={},.draw="COLZ",.opt=o_1or2d_d, .ranges={}});
// Deep AK8 tagger plots
sh.AddHistos("AK8", { .fill="JetAK8DeepTagTvsQCD", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagWvsQCD", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagZvsQCD", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagH", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagMDHbbvsQCD",.pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagMDWvsQCD", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
sh.AddHistos("AK8", { .fill="JetAK8DeepTagMDTvsQCD", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk8});
// Tagged Ws
// Add W scale factor
sh.AddHistos("hadw", { .fill="WPtBins", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadw", { .fill="WPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadw", { .fill="WEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadw", { .fill="WPhi", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadw", { .fill="WTau21", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=Opt,.ranges=r_Stk4});
sh.AddHistos("hadw", { .fill="WMass", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
// Tagged tops
// Add top scale factor
sh.AddHistos("hadtop", { .fill="TopPtBins", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadtop", { .fill="TopPt", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadtop", { .fill="TopEta", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadtop", { .fill="TopPhi", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadtop", { .fill="TopTau32", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
sh.AddHistos("hadtop", { .fill="TopMass", .pfs={Stack,"Year",cut}, .cuts={},.draw=d,.opt=opt,.ranges=r_Stk4});
}
}
//std::cout<<"-------------------------------------------------------"<<std::endl;
//std::cout<<" Total cell count for each SmartHisto:"<<std::endl;
//sh.GetTotalNCells();
//std::cout<<"-------------------------------------------------------"<<std::endl;
}
#endif // End header guard