Context Navigation

Back to MA5SandBox

MA5SandBox: cms_sus_19_006.cpp

File cms_sus_19_006.cpp, 77.2 KB (added by Benjamin Fuks, 5 years ago)

Line
1	#include "SampleAnalyzer/User/Analyzer/cms_sus_19_006.h"
2
3	using namespace MA5;
4	using namespace std;
5
6	// -----------------------------------------------------------------------------
7	// Initialize
8	// function called one time at the beginning of the analysis
9	// -----------------------------------------------------------------------------
10	bool cms_sus_19_006::Initialize(const MA5::Configuration& cfg, const std::map<std::string,std::string>& parameters)
11	{
12	cout << "BEGIN Initialization" << endl;
13	// initialize variables, histos
14	// Information on the analysis, authors, ...
15	// VERY IMPORTANT FOR DOCUMENTATION, TRACEABILITY, BUG REPORTS
16	INFO << " <><><><><><><><><><><><><><><><><><><><><><><><>" << endmsg;
17	INFO << " <> Analysis: cms_sus_19_006 <>" << endmsg;
18	INFO << " <> Multijet,HT,MHT = 13 TeV, 137 fb^-1 luminosity <>" << endmsg;
19	INFO << " <> Recast by: M.Mrowietz, S.Bein <>" << endmsg;
20	INFO << " <> Contact: m.mrowietz@gmail.com , samuel.bein@gmail.com <>" << endmsg;
21	INFO << " <> Based on MadAnalysis 5 v1.7 and above <>" << endmsg;
22	INFO << " <> For more information, see <>" << endmsg;
23	INFO << " <> http://madanalysis.irmp.ucl.ac.be/wiki/PublicAnalysisDatabase <>" << endmsg;
24	// Define the 174 signal regions
25	//2-3 jets, 0bjets
26	Manager()->AddRegionSelection("SR1_Njet23_Nb0_HT300600_MHT300350");
27	Manager()->AddRegionSelection("SR2_Njet23_Nb0_HT6001200_MHT300350");
28	Manager()->AddRegionSelection("SR3_Njet23_Nb0_HT1200_MHT300350");
29	Manager()->AddRegionSelection("SR4_Njet23_Nb0_HT350600_MHT350600");
30	Manager()->AddRegionSelection("SR5_Njet23_Nb0_HT6001200_MHT350600");
31	Manager()->AddRegionSelection("SR6_Njet23_Nb0_HT1200_MHT350600");
32	Manager()->AddRegionSelection("SR7_Njet23_Nb0_HT6001200_MHT600850");
33	Manager()->AddRegionSelection("SR8_Njet23_Nb0_HT1200_MHT600850");
34	Manager()->AddRegionSelection("SR9_Njet23_Nb0_HT8501700_MHT850");
35	Manager()->AddRegionSelection("SR10_Njet23_Nb0_HT1700_MHT850");
36	// 2-3 jets, 1 bjet
37	Manager()->AddRegionSelection("SR11_Njet23_Nb1_HT300600_MHT300350");
38	Manager()->AddRegionSelection("SR12_Njet23_Nb1_HT6001200_MHT300350");
39	Manager()->AddRegionSelection("SR13_Njet23_Nb1_HT1200_MHT300350");
40	Manager()->AddRegionSelection("SR14_Njet23_Nb1_HT350600_MHT350600");
41	Manager()->AddRegionSelection("SR15_Njet23_Nb1_HT6001200_MHT350600");
42	Manager()->AddRegionSelection("SR16_Njet23_Nb1_HT1200_MHT350600");
43	Manager()->AddRegionSelection("SR17_Njet23_Nb1_HT6001200_MHT600850");
44	Manager()->AddRegionSelection("SR18_Njet23_Nb1_HT1200_MHT600850");
45	Manager()->AddRegionSelection("SR19_Njet23_Nb1_HT8501700_MHT850");
46	Manager()->AddRegionSelection("SR20_Njet23_Nb1_HT1700_MHT850");
47	// 2-3 jets, 2 bjets
48	Manager()->AddRegionSelection("SR21_Njet23_Nb2_HT300600_MHT300350");
49	Manager()->AddRegionSelection("SR22_Njet23_Nb2_HT6001200_MHT300350");
50	Manager()->AddRegionSelection("SR23_Njet23_Nb2_HT1200_MHT300350");
51	Manager()->AddRegionSelection("SR24_Njet23_Nb2_HT350600_MHT350600");
52	Manager()->AddRegionSelection("SR25_Njet23_Nb2_HT6001200_MHT350600");
53	Manager()->AddRegionSelection("SR26_Njet23_Nb2_HT1200_MHT350600");
54	Manager()->AddRegionSelection("SR27_Njet23_Nb2_HT6001200_MHT600850");
55	Manager()->AddRegionSelection("SR28_Njet23_Nb2_HT1200_MHT600850");
56	Manager()->AddRegionSelection("SR29_Njet23_Nb2_HT8501700_MHT850");
57	Manager()->AddRegionSelection("SR30_Njet23_Nb2_HT1700_MHT850");
58
59	//4-5 jets, 0bjets
60	Manager()->AddRegionSelection("SR31_Njet45_Nb0_HT300600_MHT300350");
61	Manager()->AddRegionSelection("SR32_Njet45_Nb0_HT6001200_MHT300350");
62	Manager()->AddRegionSelection("SR33_Njet45_Nb0_HT1200_MHT300350");
63	Manager()->AddRegionSelection("SR34_Njet45_Nb0_HT350600_MHT350600");
64	Manager()->AddRegionSelection("SR35_Njet45_Nb0_HT6001200_MHT350600");
65	Manager()->AddRegionSelection("SR36_Njet45_Nb0_HT1200_MHT350600");
66	Manager()->AddRegionSelection("SR37_Njet45_Nb0_HT6001200_MHT600850");
67	Manager()->AddRegionSelection("SR38_Njet45_Nb0_HT1200_MHT600850");
68	Manager()->AddRegionSelection("SR39_Njet45_Nb0_HT8501700_MHT850");
69	Manager()->AddRegionSelection("SR40_Njet45_Nb0_HT1700_MHT850");
70	//4-5 jets, 1bjets
71	Manager()->AddRegionSelection("SR41_Njet45_Nb1_HT300600_MHT300350");
72	Manager()->AddRegionSelection("SR42_Njet45_Nb1_HT6001200_MHT300350");
73	Manager()->AddRegionSelection("SR43_Njet45_Nb1_HT1200_MHT300350");
74	Manager()->AddRegionSelection("SR44_Njet45_Nb1_HT350600_MHT350600");
75	Manager()->AddRegionSelection("SR45_Njet45_Nb1_HT6001200_MHT350600");
76	Manager()->AddRegionSelection("SR46_Njet45_Nb1_HT1200_MHT350600");
77	Manager()->AddRegionSelection("SR47_Njet45_Nb1_HT6001200_MHT600850");
78	Manager()->AddRegionSelection("SR48_Njet45_Nb1_HT1200_MHT600850");
79	Manager()->AddRegionSelection("SR49_Njet45_Nb1_HT8501700_MHT850");
80	Manager()->AddRegionSelection("SR50_Njet45_Nb1_HT1700_MHT850");
81	//4-5 jets, 2bjets
82	Manager()->AddRegionSelection("SR51_Njet45_Nb2_HT300600_MHT300350");
83	Manager()->AddRegionSelection("SR52_Njet45_Nb2_HT6001200_MHT300350");
84	Manager()->AddRegionSelection("SR53_Njet45_Nb2_HT1200_MHT300350");
85	Manager()->AddRegionSelection("SR54_Njet45_Nb2_HT350600_MHT350600");
86	Manager()->AddRegionSelection("SR55_Njet45_Nb2_HT6001200_MHT350600");
87	Manager()->AddRegionSelection("SR56_Njet45_Nb2_HT1200_MHT350600");
88	Manager()->AddRegionSelection("SR57_Njet45_Nb2_HT6001200_MHT600850");
89	Manager()->AddRegionSelection("SR58_Njet45_Nb2_HT1200_MHT600850");
90	Manager()->AddRegionSelection("SR59_Njet45_Nb2_HT8501700_MHT850");
91	Manager()->AddRegionSelection("SR60_Njet45_Nb2_HT1700_MHT850");
92	//4-5 jets, 3bjets
93	Manager()->AddRegionSelection("SR61_Njet45_Nb3_HT300600_MHT300350");
94	Manager()->AddRegionSelection("SR62_Njet45_Nb3_HT6001200_MHT300350");
95	Manager()->AddRegionSelection("SR63_Njet45_Nb3_HT1200_MHT300350");
96	Manager()->AddRegionSelection("SR64_Njet45_Nb3_HT350600_MHT350600");
97	Manager()->AddRegionSelection("SR65_Njet45_Nb3_HT6001200_MHT350600");
98	Manager()->AddRegionSelection("SR66_Njet45_Nb3_HT1200_MHT350600");
99	Manager()->AddRegionSelection("SR67_Njet45_Nb3_HT6001200_MHT600850");
100	Manager()->AddRegionSelection("SR68_Njet45_Nb3_HT1200_MHT600850");
101	Manager()->AddRegionSelection("SR69_Njet45_Nb3_HT8501700_MHT850");
102	Manager()->AddRegionSelection("SR70_Njet45_Nb3_HT1700_MHT850");
103
104
105	//6-7 jets, 0bjets
106	Manager()->AddRegionSelection("SR71_Njet67_Nb0_HT300600_MHT300350");
107	Manager()->AddRegionSelection("SR72_Njet67_Nb0_HT6001200_MHT300350");
108	Manager()->AddRegionSelection("SR73_Njet67_Nb0_HT1200_MHT300350");
109	Manager()->AddRegionSelection("SR74_Njet67_Nb0_HT350600_MHT350600");
110	Manager()->AddRegionSelection("SR75_Njet67_Nb0_HT6001200_MHT350600");
111	Manager()->AddRegionSelection("SR76_Njet67_Nb0_HT1200_MHT350600");
112	Manager()->AddRegionSelection("SR77_Njet67_Nb0_HT6001200_MHT600850");
113	Manager()->AddRegionSelection("SR78_Njet67_Nb0_HT1200_MHT600850");
114	Manager()->AddRegionSelection("SR79_Njet67_Nb0_HT8501700_MHT850");
115	Manager()->AddRegionSelection("SR80_Njet67_Nb0_HT1700_MHT850");
116	//6-7 jets, 1bjets
117	Manager()->AddRegionSelection("SR81_Njet67_Nb1_HT300600_MHT300350");
118	Manager()->AddRegionSelection("SR82_Njet67_Nb1_HT6001200_MHT300350");
119	Manager()->AddRegionSelection("SR83_Njet67_Nb1_HT1200_MHT300350");
120	Manager()->AddRegionSelection("SR84_Njet67_Nb1_HT350600_MHT350600");
121	Manager()->AddRegionSelection("SR85_Njet67_Nb1_HT6001200_MHT350600");
122	Manager()->AddRegionSelection("SR86_Njet67_Nb1_HT1200_MHT350600");
123	Manager()->AddRegionSelection("SR87_Njet67_Nb1_HT6001200_MHT600850");
124	Manager()->AddRegionSelection("SR88_Njet67_Nb1_HT1200_MHT600850");
125	Manager()->AddRegionSelection("SR89_Njet67_Nb1_HT8501700_MHT850");
126	Manager()->AddRegionSelection("SR90_Njet67_Nb1_HT1700_MHT850");
127	//6-7 jets, 2bjets
128	Manager()->AddRegionSelection("SR91_Njet67_Nb2_HT300600_MHT300350");
129	Manager()->AddRegionSelection("SR92_Njet67_Nb2_HT6001200_MHT300350");
130	Manager()->AddRegionSelection("SR93_Njet67_Nb2_HT1200_MHT300350");
131	Manager()->AddRegionSelection("SR94_Njet67_Nb2_HT350600_MHT350600");
132	Manager()->AddRegionSelection("SR95_Njet67_Nb2_HT6001200_MHT350600");
133	Manager()->AddRegionSelection("SR96_Njet67_Nb2_HT1200_MHT350600");
134	Manager()->AddRegionSelection("SR97_Njet67_Nb2_HT6001200_MHT600850");
135	Manager()->AddRegionSelection("SR98_Njet67_Nb2_HT1200_MHT600850");
136	Manager()->AddRegionSelection("SR99_Njet67_Nb2_HT8501700_MHT850");
137	Manager()->AddRegionSelection("SR100_Njet67_Nb2_HT1700_MHT850");
138	//6-7 jets, 3bjets
139	Manager()->AddRegionSelection("SR101_Njet67_Nb3_HT300600_MHT300350");
140	Manager()->AddRegionSelection("SR102_Njet67_Nb3_HT6001200_MHT300350");
141	Manager()->AddRegionSelection("SR103_Njet67_Nb3_HT1200_MHT300350");
142	Manager()->AddRegionSelection("SR104_Njet67_Nb3_HT350600_MHT350600");
143	Manager()->AddRegionSelection("SR105_Njet67_Nb3_HT6001200_MHT350600");
144	Manager()->AddRegionSelection("SR106_Njet67_Nb3_HT1200_MHT350600");
145	Manager()->AddRegionSelection("SR107_Njet67_Nb3_HT6001200_MHT600850");
146	Manager()->AddRegionSelection("SR108_Njet67_Nb3_HT1200_MHT600850");
147	Manager()->AddRegionSelection("SR109_Njet67_Nb3_HT8501700_MHT850");
148	Manager()->AddRegionSelection("SR110_Njet67_Nb3_HT1700_MHT850");
149
150	//8-9 jets, 0bjets
151	Manager()->AddRegionSelection("SR111_Njet89_Nb0_HT6001200_MHT300350");
152	Manager()->AddRegionSelection("SR112_Njet89_Nb0_HT1200_MHT300350");
153	Manager()->AddRegionSelection("SR113_Njet89_Nb0_HT6001200_MHT350600");
154	Manager()->AddRegionSelection("SR114_Njet89_Nb0_HT1200_MHT350600");
155	Manager()->AddRegionSelection("SR115_Njet89_Nb0_HT6001200_MHT600850");
156	Manager()->AddRegionSelection("SR116_Njet89_Nb0_HT1200_MHT600850");
157	Manager()->AddRegionSelection("SR117_Njet89_Nb0_HT8501700_MHT850");
158	Manager()->AddRegionSelection("SR118_Njet89_Nb0_HT1700_MHT850");
159	//8-9 jets, 1bjets
160	Manager()->AddRegionSelection("SR119_Njet89_Nb1_HT6001200_MHT300350");
161	Manager()->AddRegionSelection("SR120_Njet89_Nb1_HT1200_MHT300350");
162	Manager()->AddRegionSelection("SR121_Njet89_Nb1_HT6001200_MHT350600");
163	Manager()->AddRegionSelection("SR122_Njet89_Nb1_HT1200_MHT350600");
164	Manager()->AddRegionSelection("SR123_Njet89_Nb1_HT6001200_MHT600850");
165	Manager()->AddRegionSelection("SR124_Njet89_Nb1_HT1200_MHT600850");
166	Manager()->AddRegionSelection("SR125_Njet89_Nb1_HT8501700_MHT850");
167	Manager()->AddRegionSelection("SR126_Njet89_Nb1_HT1700_MHT850");
168	//8-9 jets, 2bjets
169	Manager()->AddRegionSelection("SR127_Njet89_Nb2_HT6001200_MHT300350");
170	Manager()->AddRegionSelection("SR128_Njet89_Nb2_HT1200_MHT300350");
171	Manager()->AddRegionSelection("SR129_Njet89_Nb2_HT6001200_MHT350600");
172	Manager()->AddRegionSelection("SR130_Njet89_Nb2_HT1200_MHT350600");
173	Manager()->AddRegionSelection("SR131_Njet89_Nb2_HT6001200_MHT600850");
174	Manager()->AddRegionSelection("SR132_Njet89_Nb2_HT1200_MHT600850");
175	Manager()->AddRegionSelection("SR133_Njet89_Nb2_HT8501700_MHT850");
176	Manager()->AddRegionSelection("SR134_Njet89_Nb2_HT1700_MHT850");
177	//8-9 jets, 3bjets
178	Manager()->AddRegionSelection("SR135_Njet89_Nb3_HT6001200_MHT300350");
179	Manager()->AddRegionSelection("SR136_Njet89_Nb3_HT1200_MHT300350");
180	Manager()->AddRegionSelection("SR137_Njet89_Nb3_HT6001200_MHT350600");
181	Manager()->AddRegionSelection("SR138_Njet89_Nb3_HT1200_MHT350600");
182	Manager()->AddRegionSelection("SR139_Njet89_Nb3_HT6001200_MHT600850");
183	Manager()->AddRegionSelection("SR140_Njet89_Nb3_HT1200_MHT600850");
184	Manager()->AddRegionSelection("SR141_Njet89_Nb3_HT8501700_MHT850");
185	Manager()->AddRegionSelection("SR142_Njet89_Nb3_HT1700_MHT850");
186
187	//>10 jets, 0bjets
188	Manager()->AddRegionSelection("SR143_Njet10_Nb0_HT6001200_MHT300350");
189	Manager()->AddRegionSelection("SR144_Njet10_Nb0_HT1200_MHT300350");
190	Manager()->AddRegionSelection("SR145_Njet10_Nb0_HT6001200_MHT350600");
191	Manager()->AddRegionSelection("SR146_Njet10_Nb0_HT1200_MHT350600");
192	Manager()->AddRegionSelection("SR147_Njet10_Nb0_HT6001200_MHT600850");
193	Manager()->AddRegionSelection("SR148_Njet10_Nb0_HT1200_MHT600850");
194	Manager()->AddRegionSelection("SR149_Njet10_Nb0_HT8501700_MHT850");
195	Manager()->AddRegionSelection("SR150_Njet10_Nb0_HT1700_MHT850");
196	//>10 jets, 1bjets
197	Manager()->AddRegionSelection("SR151_Njet10_Nb1_HT6001200_MHT300350");
198	Manager()->AddRegionSelection("SR152_Njet10_Nb1_HT1200_MHT300350");
199	Manager()->AddRegionSelection("SR153_Njet10_Nb1_HT6001200_MHT350600");
200	Manager()->AddRegionSelection("SR154_Njet10_Nb1_HT1200_MHT350600");
201	Manager()->AddRegionSelection("SR155_Njet10_Nb1_HT6001200_MHT600850");
202	Manager()->AddRegionSelection("SR156_Njet10_Nb1_HT1200_MHT600850");
203	Manager()->AddRegionSelection("SR157_Njet10_Nb1_HT8501700_MHT850");
204	Manager()->AddRegionSelection("SR158_Njet10_Nb1_HT1700_MHT850");
205	//>10 jets, 2bjets
206	Manager()->AddRegionSelection("SR159_Njet10_Nb2_HT6001200_MHT300350");
207	Manager()->AddRegionSelection("SR160_Njet10_Nb2_HT1200_MHT300350");
208	Manager()->AddRegionSelection("SR161_Njet10_Nb2_HT6001200_MHT350600");
209	Manager()->AddRegionSelection("SR162_Njet10_Nb2_HT1200_MHT350600");
210	Manager()->AddRegionSelection("SR163_Njet10_Nb2_HT6001200_MHT600850");
211	Manager()->AddRegionSelection("SR164_Njet10_Nb2_HT1200_MHT600850");
212	Manager()->AddRegionSelection("SR165_Njet10_Nb2_HT8501700_MHT850");
213	Manager()->AddRegionSelection("SR166_Njet10_Nb2_HT1700_MHT850");
214	//>10 jets, 3bjets
215	Manager()->AddRegionSelection("SR167_Njet10_Nb3_HT6001200_MHT300350");
216	Manager()->AddRegionSelection("SR168_Njet10_Nb3_HT1200_MHT300350");
217	Manager()->AddRegionSelection("SR169_Njet10_Nb3_HT6001200_MHT350600");
218	Manager()->AddRegionSelection("SR170_Njet10_Nb3_HT1200_MHT350600");
219	Manager()->AddRegionSelection("SR171_Njet10_Nb3_HT6001200_MHT600850");
220	Manager()->AddRegionSelection("SR172_Njet10_Nb3_HT1200_MHT600850");
221	Manager()->AddRegionSelection("SR173_Njet10_Nb3_HT8501700_MHT850");
222	Manager()->AddRegionSelection("SR174_Njet10_Nb3_HT1700_MHT850");
223
224
225
226	// Define the 12 aggregated signal region
227	Manager()->AddRegionSelection("AGGSR1_Njet2_Nb0_HT600_MHT600");
228	Manager()->AddRegionSelection("AGGSR2_Njet4_Nb0_HT1700_MHT850");
229	Manager()->AddRegionSelection("AGGSR3_Njet6_Nb0_HT600_MHT_600");
230	Manager()->AddRegionSelection("AGGSR4_Njet8_Nb01_HT600_MHT600");
231	Manager()->AddRegionSelection("AGGSR5_Njet10_Nb01_HT1700_MHT850");
232	Manager()->AddRegionSelection("AGGSR6_Njet4_Nb2_HT300_MHT300");
233	Manager()->AddRegionSelection("AGGSR7_Njet2_Nb2_HT600_MHT600");
234	Manager()->AddRegionSelection("AGGSR8_Njet6_Nb2_HT350_MHT350");
235	Manager()->AddRegionSelection("AGGSR9_NJet4_Nb2_HT600_MHT600");
236	Manager()->AddRegionSelection("AGGSR10_Njet8_Nb3_HT300_MHT300");
237	Manager()->AddRegionSelection("AGGSR11_Njet6_Nb1_HT600_MHT600");
238	Manager()->AddRegionSelection("AGGSR12_Njet10_Nb3_HT850_MHT850");
239
240	// Define the preselection cuts (not including Jet,Nb,HT,MHT binnings)
241	Manager()->AddCut("Njet>=2") ; // more or equal 2 jets with \|eta\|<2.4
242	Manager()->AddCut("HT>300") ; // HT = scalar sum of jets pt>30 GeV
243	Manager()->AddCut("MHT>300") ; // HTM = magnitude of the the vector HTMiss, i.e. negative of the vector sum of all jets pt within \|eta\|<5
244	Manager()->AddCut("HT>MHT") ; // HTM < HT is vetoed, as it is likely to arise from mismeasurement
245	Manager()->AddCut("NoIsoMuons") ; // No isolated muons with pt>10 GeV
246	Manager()->AddCut("NoMuonsTracks") ; // No isolated muons tracks
247	Manager()->AddCut("NoIsoElectrons") ; // No isolated electron with pt>10 GeV
248	Manager()->AddCut("NoElectronsTracks") ; // No isolated electron tracks
249	Manager()->AddCut("NoIsoTracks") ; // no isolated tracks with mT<100 GeV and pT>10 GeV
250	// (pT > 5 GeV if track is a PF electron or muon)
251	// mT = transverse mass formed with the pTmiss (all pT of the PF objects) and the iso track
252	Manager()->AddCut("NoIsoPhotons") ; // no isolated photon with pT> 100 GeV to retain orthogonality to other SUSY searches
253
254
255	Manager()->AddCut("DPhi(MHTj1)>0.5") ; // Azimutal angle between MHT and the pT of the 'i' jet
256	Manager()->AddCut("DPhi(MHTj2)>0.5") ;
257	Manager()->AddCut("DPhi(MHTj3)>0.3") ;
258	Manager()->AddCut("DPhi(MHTj4)>0.3") ;
259
260	// Jets multiplicity selection
261
262	// the 174 signal regions
263	//2-3 jets, 0bjets
264	Manager()->AddCut("SR1_Njet23_Nb0_HT300600_MHT300350","SR1_Njet23_Nb0_HT300600_MHT300350");
265	Manager()->AddCut("SR2_Njet23_Nb0_HT6001200_MHT300350","SR2_Njet23_Nb0_HT6001200_MHT300350");
266	Manager()->AddCut("SR3_Njet23_Nb0_HT1200_MHT300350","SR3_Njet23_Nb0_HT1200_MHT300350");
267	Manager()->AddCut("SR4_Njet23_Nb0_HT350600_MHT350600","SR4_Njet23_Nb0_HT350600_MHT350600");
268	Manager()->AddCut("SR5_Njet23_Nb0_HT6001200_MHT350600","SR5_Njet23_Nb0_HT6001200_MHT350600");
269	Manager()->AddCut("SR6_Njet23_Nb0_HT1200_MHT350600","SR6_Njet23_Nb0_HT1200_MHT350600");
270	Manager()->AddCut("SR7_Njet23_Nb0_HT6001200_MHT600850","SR7_Njet23_Nb0_HT6001200_MHT600850");
271	Manager()->AddCut("SR8_Njet23_Nb0_HT1200_MHT600850","SR8_Njet23_Nb0_HT1200_MHT600850");
272	Manager()->AddCut("SR9_Njet23_Nb0_HT8501700_MHT850","SR9_Njet23_Nb0_HT8501700_MHT850");
273	Manager()->AddCut("SR10_Njet23_Nb0_HT1700_MHT850","SR10_Njet23_Nb0_HT1700_MHT850");
274	// 2-3 jets, 1 bjet
275	Manager()->AddCut("SR11_Njet23_Nb1_HT300600_MHT300350","SR11_Njet23_Nb1_HT300600_MHT300350");
276	Manager()->AddCut("SR12_Njet23_Nb1_HT6001200_MHT300350","SR12_Njet23_Nb1_HT6001200_MHT300350");
277	Manager()->AddCut("SR13_Njet23_Nb1_HT1200_MHT300350","SR13_Njet23_Nb1_HT1200_MHT300350");
278	Manager()->AddCut("SR14_Njet23_Nb1_HT350600_MHT350600","SR14_Njet23_Nb1_HT350600_MHT350600");
279	Manager()->AddCut("SR15_Njet23_Nb1_HT6001200_MHT350600","SR15_Njet23_Nb1_HT6001200_MHT350600");
280	Manager()->AddCut("SR16_Njet23_Nb1_HT1200_MHT350600","SR16_Njet23_Nb1_HT1200_MHT350600");
281	Manager()->AddCut("SR17_Njet23_Nb1_HT6001200_MHT600850","SR17_Njet23_Nb1_HT6001200_MHT600850");
282	Manager()->AddCut("SR18_Njet23_Nb1_HT1200_MHT600850","SR18_Njet23_Nb1_HT1200_MHT600850");
283	Manager()->AddCut("SR19_Njet23_Nb1_HT8501700_MHT850","SR19_Njet23_Nb1_HT8501700_MHT850");
284	Manager()->AddCut("SR20_Njet23_Nb1_HT1700_MHT850","SR20_Njet23_Nb1_HT1700_MHT850");
285	// 2-3 jets, 2 bjets
286	Manager()->AddCut("SR21_Njet23_Nb2_HT300600_MHT300350","SR21_Njet23_Nb2_HT300600_MHT300350");
287	Manager()->AddCut("SR22_Njet23_Nb2_HT6001200_MHT300350","SR22_Njet23_Nb2_HT6001200_MHT300350");
288	Manager()->AddCut("SR23_Njet23_Nb2_HT1200_MHT300350","SR23_Njet23_Nb2_HT1200_MHT300350");
289	Manager()->AddCut("SR24_Njet23_Nb2_HT350600_MHT350600","SR24_Njet23_Nb2_HT350600_MHT350600");
290	Manager()->AddCut("SR25_Njet23_Nb2_HT6001200_MHT350600","SR25_Njet23_Nb2_HT6001200_MHT350600");
291	Manager()->AddCut("SR26_Njet23_Nb2_HT1200_MHT350600","SR26_Njet23_Nb2_HT1200_MHT350600");
292	Manager()->AddCut("SR27_Njet23_Nb2_HT6001200_MHT600850","SR27_Njet23_Nb2_HT6001200_MHT600850");
293	Manager()->AddCut("SR28_Njet23_Nb2_HT1200_MHT600850","SR28_Njet23_Nb2_HT1200_MHT600850");
294	Manager()->AddCut("SR29_Njet23_Nb2_HT8501700_MHT850","SR29_Njet23_Nb2_HT8501700_MHT850");
295	Manager()->AddCut("SR30_Njet23_Nb2_HT1700_MHT850","SR30_Njet23_Nb2_HT1700_MHT850");
296
297	//4-5 jets, 0bjets
298	Manager()->AddCut("SR31_Njet45_Nb0_HT300600_MHT300350","SR31_Njet45_Nb0_HT300600_MHT300350");
299	Manager()->AddCut("SR32_Njet45_Nb0_HT6001200_MHT300350","SR32_Njet45_Nb0_HT6001200_MHT300350");
300	Manager()->AddCut("SR33_Njet45_Nb0_HT1200_MHT300350","SR33_Njet45_Nb0_HT1200_MHT300350");
301	Manager()->AddCut("SR34_Njet45_Nb0_HT350600_MHT350600","SR34_Njet45_Nb0_HT350600_MHT350600");
302	Manager()->AddCut("SR35_Njet45_Nb0_HT6001200_MHT350600","SR35_Njet45_Nb0_HT6001200_MHT350600");
303	Manager()->AddCut("SR36_Njet45_Nb0_HT1200_MHT350600","SR36_Njet45_Nb0_HT1200_MHT350600");
304	Manager()->AddCut("SR37_Njet45_Nb0_HT6001200_MHT600850","SR37_Njet45_Nb0_HT6001200_MHT600850");
305	Manager()->AddCut("SR38_Njet45_Nb0_HT1200_MHT600850","SR38_Njet45_Nb0_HT1200_MHT600850");
306	Manager()->AddCut("SR39_Njet45_Nb0_HT8501700_MHT850","SR39_Njet45_Nb0_HT8501700_MHT850");
307	Manager()->AddCut("SR40_Njet45_Nb0_HT1700_MHT850","SR40_Njet45_Nb0_HT1700_MHT850");
308	//4-5 jets, 1bjets
309	Manager()->AddCut("SR41_Njet45_Nb1_HT300600_MHT300350","SR41_Njet45_Nb1_HT300600_MHT300350");
310	Manager()->AddCut("SR42_Njet45_Nb1_HT6001200_MHT300350","SR42_Njet45_Nb1_HT6001200_MHT300350");
311	Manager()->AddCut("SR43_Njet45_Nb1_HT1200_MHT300350","SR43_Njet45_Nb1_HT1200_MHT300350");
312	Manager()->AddCut("SR44_Njet45_Nb1_HT350600_MHT350600","SR44_Njet45_Nb1_HT350600_MHT350600");
313	Manager()->AddCut("SR45_Njet45_Nb1_HT6001200_MHT350600","SR45_Njet45_Nb1_HT6001200_MHT350600");
314	Manager()->AddCut("SR46_Njet45_Nb1_HT1200_MHT350600","SR46_Njet45_Nb1_HT1200_MHT350600");
315	Manager()->AddCut("SR47_Njet45_Nb1_HT6001200_MHT600850","SR47_Njet45_Nb1_HT6001200_MHT600850");
316	Manager()->AddCut("SR48_Njet45_Nb1_HT1200_MHT600850","SR48_Njet45_Nb1_HT1200_MHT600850");
317	Manager()->AddCut("SR49_Njet45_Nb1_HT8501700_MHT850","SR49_Njet45_Nb1_HT8501700_MHT850");
318	Manager()->AddCut("SR50_Njet45_Nb1_HT1700_MHT850","SR50_Njet45_Nb1_HT1700_MHT850");
319	//4-5 jets, 2bjets
320	Manager()->AddCut("SR51_Njet45_Nb2_HT300600_MHT300350","SR51_Njet45_Nb2_HT300600_MHT300350");
321	Manager()->AddCut("SR52_Njet45_Nb2_HT6001200_MHT300350","SR52_Njet45_Nb2_HT6001200_MHT300350");
322	Manager()->AddCut("SR53_Njet45_Nb2_HT1200_MHT300350","SR53_Njet45_Nb2_HT1200_MHT300350");
323	Manager()->AddCut("SR54_Njet45_Nb2_HT350600_MHT350600","SR54_Njet45_Nb2_HT350600_MHT350600");
324	Manager()->AddCut("SR55_Njet45_Nb2_HT6001200_MHT350600","SR55_Njet45_Nb2_HT6001200_MHT350600");
325	Manager()->AddCut("SR56_Njet45_Nb2_HT1200_MHT350600","SR56_Njet45_Nb2_HT1200_MHT350600");
326	Manager()->AddCut("SR57_Njet45_Nb2_HT6001200_MHT600850","SR57_Njet45_Nb2_HT6001200_MHT600850");
327	Manager()->AddCut("SR58_Njet45_Nb2_HT1200_MHT600850","SR58_Njet45_Nb2_HT1200_MHT600850");
328	Manager()->AddCut("SR59_Njet45_Nb2_HT8501700_MHT850","SR59_Njet45_Nb2_HT8501700_MHT850");
329	Manager()->AddCut("SR60_Njet45_Nb2_HT1700_MHT850","SR60_Njet45_Nb2_HT1700_MHT850");
330	//4-5 jets, 3bjets
331	Manager()->AddCut("SR61_Njet45_Nb3_HT300600_MHT300350","SR61_Njet45_Nb3_HT300600_MHT300350");
332	Manager()->AddCut("SR62_Njet45_Nb3_HT6001200_MHT300350","SR62_Njet45_Nb3_HT6001200_MHT300350");
333	Manager()->AddCut("SR63_Njet45_Nb3_HT1200_MHT300350","SR63_Njet45_Nb3_HT1200_MHT300350");
334	Manager()->AddCut("SR64_Njet45_Nb3_HT350600_MHT350600","SR64_Njet45_Nb3_HT350600_MHT350600");
335	Manager()->AddCut("SR65_Njet45_Nb3_HT6001200_MHT350600","SR65_Njet45_Nb3_HT6001200_MHT350600");
336	Manager()->AddCut("SR66_Njet45_Nb3_HT1200_MHT350600","SR66_Njet45_Nb3_HT1200_MHT350600");
337	Manager()->AddCut("SR67_Njet45_Nb3_HT6001200_MHT600850","SR67_Njet45_Nb3_HT6001200_MHT600850");
338	Manager()->AddCut("SR68_Njet45_Nb3_HT1200_MHT600850","SR68_Njet45_Nb3_HT1200_MHT600850");
339	Manager()->AddCut("SR69_Njet45_Nb3_HT8501700_MHT850","SR69_Njet45_Nb3_HT8501700_MHT850");
340	Manager()->AddCut("SR70_Njet45_Nb3_HT1700_MHT850","SR70_Njet45_Nb3_HT1700_MHT850");
341
342
343	//6-7 jets, 0bjets
344	Manager()->AddCut("SR71_Njet67_Nb0_HT300600_MHT300350","SR71_Njet67_Nb0_HT300600_MHT300350");
345	Manager()->AddCut("SR72_Njet67_Nb0_HT6001200_MHT300350","SR72_Njet67_Nb0_HT6001200_MHT300350");
346	Manager()->AddCut("SR73_Njet67_Nb0_HT1200_MHT300350","SR73_Njet67_Nb0_HT1200_MHT300350");
347	Manager()->AddCut("SR74_Njet67_Nb0_HT350600_MHT350600","SR74_Njet67_Nb0_HT350600_MHT350600");
348	Manager()->AddCut("SR75_Njet67_Nb0_HT6001200_MHT350600","SR75_Njet67_Nb0_HT6001200_MHT350600");
349	Manager()->AddCut("SR76_Njet67_Nb0_HT1200_MHT350600","SR76_Njet67_Nb0_HT1200_MHT350600");
350	Manager()->AddCut("SR77_Njet67_Nb0_HT6001200_MHT600850","SR77_Njet67_Nb0_HT6001200_MHT600850");
351	Manager()->AddCut("SR78_Njet67_Nb0_HT1200_MHT600850","SR78_Njet67_Nb0_HT1200_MHT600850");
352	Manager()->AddCut("SR79_Njet67_Nb0_HT8501700_MHT850","SR79_Njet67_Nb0_HT8501700_MHT850");
353	Manager()->AddCut("SR80_Njet67_Nb0_HT1700_MHT850","SR80_Njet67_Nb0_HT1700_MHT850");
354	//6-7 jets, 1bjets
355	Manager()->AddCut("SR81_Njet67_Nb1_HT300600_MHT300350","SR81_Njet67_Nb1_HT300600_MHT300350");
356	Manager()->AddCut("SR82_Njet67_Nb1_HT6001200_MHT300350","SR82_Njet67_Nb1_HT6001200_MHT300350");
357	Manager()->AddCut("SR83_Njet67_Nb1_HT1200_MHT300350","SR83_Njet67_Nb1_HT1200_MHT300350");
358	Manager()->AddCut("SR84_Njet67_Nb1_HT350600_MHT350600","SR84_Njet67_Nb1_HT350600_MHT350600");
359	Manager()->AddCut("SR85_Njet67_Nb1_HT6001200_MHT350600","SR85_Njet67_Nb1_HT6001200_MHT350600");
360	Manager()->AddCut("SR86_Njet67_Nb1_HT1200_MHT350600","SR86_Njet67_Nb1_HT1200_MHT350600");
361	Manager()->AddCut("SR87_Njet67_Nb1_HT6001200_MHT600850","SR87_Njet67_Nb1_HT6001200_MHT600850");
362	Manager()->AddCut("SR88_Njet67_Nb1_HT1200_MHT600850","SR88_Njet67_Nb1_HT1200_MHT600850");
363	Manager()->AddCut("SR89_Njet67_Nb1_HT8501700_MHT850","SR89_Njet67_Nb1_HT8501700_MHT850");
364	Manager()->AddCut("SR90_Njet67_Nb1_HT1700_MHT850","SR90_Njet67_Nb1_HT1700_MHT850");
365	//6-7 jets, 2bjets
366	Manager()->AddCut("SR91_Njet67_Nb2_HT300600_MHT300350","SR91_Njet67_Nb2_HT300600_MHT300350");
367	Manager()->AddCut("SR92_Njet67_Nb2_HT6001200_MHT300350","SR92_Njet67_Nb2_HT6001200_MHT300350");
368	Manager()->AddCut("SR93_Njet67_Nb2_HT1200_MHT300350","SR93_Njet67_Nb2_HT1200_MHT300350");
369	Manager()->AddCut("SR94_Njet67_Nb2_HT350600_MHT350600","SR94_Njet67_Nb2_HT350600_MHT350600");
370	Manager()->AddCut("SR95_Njet67_Nb2_HT6001200_MHT350600","SR95_Njet67_Nb2_HT6001200_MHT350600");
371	Manager()->AddCut("SR96_Njet67_Nb2_HT1200_MHT350600","SR96_Njet67_Nb2_HT1200_MHT350600");
372	Manager()->AddCut("SR97_Njet67_Nb2_HT6001200_MHT600850","SR97_Njet67_Nb2_HT6001200_MHT600850");
373	Manager()->AddCut("SR98_Njet67_Nb2_HT1200_MHT600850","SR98_Njet67_Nb2_HT1200_MHT600850");
374	Manager()->AddCut("SR99_Njet67_Nb2_HT8501700_MHT850","SR99_Njet67_Nb2_HT8501700_MHT850");
375	Manager()->AddCut("SR100_Njet67_Nb2_HT1700_MHT850","SR100_Njet67_Nb2_HT1700_MHT850");
376	//6-7 jets, 3bjets
377	Manager()->AddCut("SR101_Njet67_Nb3_HT300600_MHT300350","SR101_Njet67_Nb3_HT300600_MHT300350");
378	Manager()->AddCut("SR102_Njet67_Nb3_HT6001200_MHT300350","SR102_Njet67_Nb3_HT6001200_MHT300350");
379	Manager()->AddCut("SR103_Njet67_Nb3_HT1200_MHT300350","SR103_Njet67_Nb3_HT1200_MHT300350");
380	Manager()->AddCut("SR104_Njet67_Nb3_HT350600_MHT350600","SR104_Njet67_Nb3_HT350600_MHT350600");
381	Manager()->AddCut("SR105_Njet67_Nb3_HT6001200_MHT350600","SR105_Njet67_Nb3_HT6001200_MHT350600");
382	Manager()->AddCut("SR106_Njet67_Nb3_HT1200_MHT350600","SR106_Njet67_Nb3_HT1200_MHT350600");
383	Manager()->AddCut("SR107_Njet67_Nb3_HT6001200_MHT600850","SR107_Njet67_Nb3_HT6001200_MHT600850");
384	Manager()->AddCut("SR108_Njet67_Nb3_HT1200_MHT600850","SR108_Njet67_Nb3_HT1200_MHT600850");
385	Manager()->AddCut("SR109_Njet67_Nb3_HT8501700_MHT850","SR109_Njet67_Nb3_HT8501700_MHT850");
386	Manager()->AddCut("SR110_Njet67_Nb3_HT1700_MHT850","SR110_Njet67_Nb3_HT1700_MHT850");
387
388	//8-9 jets, 0bjets
389	Manager()->AddCut("SR111_Njet89_Nb0_HT6001200_MHT300350","SR111_Njet89_Nb0_HT6001200_MHT300350");
390	Manager()->AddCut("SR112_Njet89_Nb0_HT1200_MHT300350","SR112_Njet89_Nb0_HT1200_MHT300350");
391	Manager()->AddCut("SR113_Njet89_Nb0_HT6001200_MHT350600","SR113_Njet89_Nb0_HT6001200_MHT350600");
392	Manager()->AddCut("SR114_Njet89_Nb0_HT1200_MHT350600","SR114_Njet89_Nb0_HT1200_MHT350600");
393	Manager()->AddCut("SR115_Njet89_Nb0_HT6001200_MHT600850","SR115_Njet89_Nb0_HT6001200_MHT600850");
394	Manager()->AddCut("SR116_Njet89_Nb0_HT1200_MHT600850","SR116_Njet89_Nb0_HT1200_MHT600850");
395	Manager()->AddCut("SR117_Njet89_Nb0_HT8501700_MHT850","SR117_Njet89_Nb0_HT8501700_MHT850");
396	Manager()->AddCut("SR118_Njet89_Nb0_HT1700_MHT850","SR118_Njet89_Nb0_HT1700_MHT850");
397	//8-9 jets, 1bjets
398	Manager()->AddCut("SR119_Njet89_Nb1_HT6001200_MHT300350","SR119_Njet89_Nb1_HT6001200_MHT300350");
399	Manager()->AddCut("SR120_Njet89_Nb1_HT1200_MHT300350","SR120_Njet89_Nb1_HT1200_MHT300350");
400	Manager()->AddCut("SR121_Njet89_Nb1_HT6001200_MHT350600","SR121_Njet89_Nb1_HT6001200_MHT350600");
401	Manager()->AddCut("SR122_Njet89_Nb1_HT1200_MHT350600","SR122_Njet89_Nb1_HT1200_MHT350600");
402	Manager()->AddCut("SR123_Njet89_Nb1_HT6001200_MHT600850","SR123_Njet89_Nb1_HT6001200_MHT600850");
403	Manager()->AddCut("SR124_Njet89_Nb1_HT1200_MHT600850","SR124_Njet89_Nb1_HT1200_MHT600850");
404	Manager()->AddCut("SR125_Njet89_Nb1_HT8501700_MHT850","SR125_Njet89_Nb1_HT8501700_MHT850");
405	Manager()->AddCut("SR126_Njet89_Nb1_HT1700_MHT850","SR126_Njet89_Nb1_HT1700_MHT850");
406	//8-9 jets, 2bjets
407	Manager()->AddCut("SR127_Njet89_Nb2_HT6001200_MHT300350","SR127_Njet89_Nb2_HT6001200_MHT300350");
408	Manager()->AddCut("SR128_Njet89_Nb2_HT1200_MHT300350","SR128_Njet89_Nb2_HT1200_MHT300350");
409	Manager()->AddCut("SR129_Njet89_Nb2_HT6001200_MHT350600","SR129_Njet89_Nb2_HT6001200_MHT350600");
410	Manager()->AddCut("SR130_Njet89_Nb2_HT1200_MHT350600","SR130_Njet89_Nb2_HT1200_MHT350600");
411	Manager()->AddCut("SR131_Njet89_Nb2_HT6001200_MHT600850","SR131_Njet89_Nb2_HT6001200_MHT600850");
412	Manager()->AddCut("SR132_Njet89_Nb2_HT1200_MHT600850","SR132_Njet89_Nb2_HT1200_MHT600850");
413	Manager()->AddCut("SR133_Njet89_Nb2_HT8501700_MHT850","SR133_Njet89_Nb2_HT8501700_MHT850");
414	Manager()->AddCut("SR134_Njet89_Nb2_HT1700_MHT850","SR134_Njet89_Nb2_HT1700_MHT850");
415	//8-9 jets, 3bjets
416	Manager()->AddCut("SR135_Njet89_Nb3_HT6001200_MHT300350","SR135_Njet89_Nb3_HT6001200_MHT300350");
417	Manager()->AddCut("SR136_Njet89_Nb3_HT1200_MHT300350","SR136_Njet89_Nb3_HT1200_MHT300350");
418	Manager()->AddCut("SR137_Njet89_Nb3_HT6001200_MHT350600","SR137_Njet89_Nb3_HT6001200_MHT350600");
419	Manager()->AddCut("SR138_Njet89_Nb3_HT1200_MHT350600","SR138_Njet89_Nb3_HT1200_MHT350600");
420	Manager()->AddCut("SR139_Njet89_Nb3_HT6001200_MHT600850","SR139_Njet89_Nb3_HT6001200_MHT600850");
421	Manager()->AddCut("SR140_Njet89_Nb3_HT1200_MHT600850","SR140_Njet89_Nb3_HT1200_MHT600850");
422	Manager()->AddCut("SR141_Njet89_Nb3_HT8501700_MHT850","SR141_Njet89_Nb3_HT8501700_MHT850");
423	Manager()->AddCut("SR142_Njet89_Nb3_HT1700_MHT850","SR142_Njet89_Nb3_HT1700_MHT850");
424
425	//>10 jets, 0bjets
426	Manager()->AddCut("SR143_Njet10_Nb0_HT6001200_MHT300350","SR143_Njet10_Nb0_HT6001200_MHT300350");
427	Manager()->AddCut("SR144_Njet10_Nb0_HT1200_MHT300350","SR144_Njet10_Nb0_HT1200_MHT300350");
428	Manager()->AddCut("SR145_Njet10_Nb0_HT6001200_MHT350600","SR145_Njet10_Nb0_HT6001200_MHT350600");
429	Manager()->AddCut("SR146_Njet10_Nb0_HT1200_MHT350600","SR146_Njet10_Nb0_HT1200_MHT350600");
430	Manager()->AddCut("SR147_Njet10_Nb0_HT6001200_MHT600850","SR147_Njet10_Nb0_HT6001200_MHT600850");
431	Manager()->AddCut("SR148_Njet10_Nb0_HT1200_MHT600850","SR148_Njet10_Nb0_HT1200_MHT600850");
432	Manager()->AddCut("SR149_Njet10_Nb0_HT8501700_MHT850","SR149_Njet10_Nb0_HT8501700_MHT850");
433	Manager()->AddCut("SR150_Njet10_Nb0_HT1700_MHT850","SR150_Njet10_Nb0_HT1700_MHT850");
434	//>10 jets, 1bjets
435	Manager()->AddCut("SR151_Njet10_Nb1_HT6001200_MHT300350","SR151_Njet10_Nb1_HT6001200_MHT300350");
436	Manager()->AddCut("SR152_Njet10_Nb1_HT1200_MHT300350","SR152_Njet10_Nb1_HT1200_MHT300350");
437	Manager()->AddCut("SR153_Njet10_Nb1_HT6001200_MHT350600","SR153_Njet10_Nb1_HT6001200_MHT350600");
438	Manager()->AddCut("SR154_Njet10_Nb1_HT1200_MHT350600","SR154_Njet10_Nb1_HT1200_MHT350600");
439	Manager()->AddCut("SR155_Njet10_Nb1_HT6001200_MHT600850","SR155_Njet10_Nb1_HT6001200_MHT600850");
440	Manager()->AddCut("SR156_Njet10_Nb1_HT1200_MHT600850","SR156_Njet10_Nb1_HT1200_MHT600850");
441	Manager()->AddCut("SR157_Njet10_Nb1_HT8501700_MHT850","SR157_Njet10_Nb1_HT8501700_MHT850");
442	Manager()->AddCut("SR158_Njet10_Nb1_HT1700_MHT850","SR158_Njet10_Nb1_HT1700_MHT850");
443	//>10 jets, 2bjets
444	Manager()->AddCut("SR159_Njet10_Nb2_HT6001200_MHT300350","SR159_Njet10_Nb2_HT6001200_MHT300350");
445	Manager()->AddCut("SR160_Njet10_Nb2_HT1200_MHT300350","SR160_Njet10_Nb2_HT1200_MHT300350");
446	Manager()->AddCut("SR161_Njet10_Nb2_HT6001200_MHT350600","SR161_Njet10_Nb2_HT6001200_MHT350600");
447	Manager()->AddCut("SR162_Njet10_Nb2_HT1200_MHT350600","SR162_Njet10_Nb2_HT1200_MHT350600");
448	Manager()->AddCut("SR163_Njet10_Nb2_HT6001200_MHT600850","SR163_Njet10_Nb2_HT6001200_MHT600850");
449	Manager()->AddCut("SR164_Njet10_Nb2_HT1200_MHT600850","SR164_Njet10_Nb2_HT1200_MHT600850");
450	Manager()->AddCut("SR165_Njet10_Nb2_HT8501700_MHT850","SR165_Njet10_Nb2_HT8501700_MHT850");
451	Manager()->AddCut("SR166_Njet10_Nb2_HT1700_MHT850","SR166_Njet10_Nb2_HT1700_MHT850");
452	//>10 jets, 3bjets
453	Manager()->AddCut("SR167_Njet10_Nb3_HT6001200_MHT300350","SR167_Njet10_Nb3_HT6001200_MHT300350");
454	Manager()->AddCut("SR168_Njet10_Nb3_HT1200_MHT300350","SR168_Njet10_Nb3_HT1200_MHT300350");
455	Manager()->AddCut("SR169_Njet10_Nb3_HT6001200_MHT350600","SR169_Njet10_Nb3_HT6001200_MHT350600");
456	Manager()->AddCut("SR170_Njet10_Nb3_HT1200_MHT350600","SR170_Njet10_Nb3_HT1200_MHT350600");
457	Manager()->AddCut("SR171_Njet10_Nb3_HT6001200_MHT600850","SR171_Njet10_Nb3_HT6001200_MHT600850");
458	Manager()->AddCut("SR172_Njet10_Nb3_HT1200_MHT600850","SR172_Njet10_Nb3_HT1200_MHT600850");
459	Manager()->AddCut("SR173_Njet10_Nb3_HT8501700_MHT850","SR173_Njet10_Nb3_HT8501700_MHT850");
460	Manager()->AddCut("SR174_Njet10_Nb3_HT1700_MHT850","SR174_Njet10_Nb3_HT1700_MHT850");
461
462
463	//aggregated regions
464	Manager()->AddCut("AGGSR1_Njet2_Nb0_HT600_MHT600","AGGSR1_Njet2_Nb0_HT600_MHT600");
465	Manager()->AddCut("AGGSR2_Njet4_Nb0_HT1700_MHT850","AGGSR2_Njet4_Nb0_HT1700_MHT850");
466	Manager()->AddCut("AGGSR3_Njet6_Nb0_HT600_MHT_600","AGGSR3_Njet6_Nb0_HT600_MHT_600");
467	Manager()->AddCut("AGGSR4_Njet8_Nb01_HT600_MHT600","AGGSR4_Njet8_Nb01_HT600_MHT600");
468	Manager()->AddCut("AGGSR5_Njet10_Nb01_HT1700_MHT850","AGGSR5_Njet10_Nb01_HT1700_MHT850");
469	Manager()->AddCut("AGGSR6_Njet4_Nb2_HT300_MHT300","AGGSR6_Njet4_Nb2_HT300_MHT300");
470	Manager()->AddCut("AGGSR7_Njet2_Nb2_HT600_MHT600","AGGSR7_Njet2_Nb2_HT600_MHT600");
471	Manager()->AddCut("AGGSR8_Njet6_Nb2_HT350_MHT350","AGGSR8_Njet6_Nb2_HT350_MHT350");
472	Manager()->AddCut("AGGSR9_NJet4_Nb2_HT600_MHT600","AGGSR9_NJet4_Nb2_HT600_MHT600");
473	Manager()->AddCut("AGGSR10_Njet8_Nb3_HT300_MHT300","AGGSR10_Njet8_Nb3_HT300_MHT300");
474	Manager()->AddCut("AGGSR11_Njet6_Nb1_HT600_MHT600","AGGSR11_Njet6_Nb1_HT600_MHT600");
475	Manager()->AddCut("AGGSR12_Njet10_Nb3_HT850_MHT850","AGGSR12_Njet10_Nb3_HT850_MHT850");
476
477	// ********* Kinematic Distributions for the points in the supplementary material )
478	// still need to test this
479	Manager()->AddHisto("nJets_Pre", 10,0,9);
480	Manager()->AddHisto("bJets_Pre", 4,-0.5,3.5);
481	Manager()->AddHisto("HT_Pre" , 40,0,3000);
482	Manager()->AddHisto("MHT_Pre" , 30,0,1500);
483
484	Manager()->AddHisto("nJets_Post", 10,0,9);
485	Manager()->AddHisto("bJets_Post", 4,-0.5,3.5);
486	Manager()->AddHisto("HT_Post" , 40,0,3000);
487	Manager()->AddHisto("MHT_Post" , 30,0,1500);
488	// cout << "END Initialization" << endl;
489
490	return true;
491	}
492
493	// -----------------------------------------------------------------------------
494	// Finalize
495	// function called one time at the end of the analysis
496	// -----------------------------------------------------------------------------
497	void cms_sus_19_006::Finalize(const SampleFormat& summary, const std::vector<SampleFormat>& files)
498	{
499	cout << "BEGIN Finalization" << endl;
500	// saving histos
501	cout << "END Finalization" << endl;
502	}
503
504	// -----------------------------------------------------------------------------
505	// Execute
506	// function called each time one event is read
507	// -----------------------------------------------------------------------------
508	bool cms_sus_19_006::Execute(SampleFormat& sample, const EventFormat& event)
509	{
510
511
512	if (event.rec()!=0)
513	{
514	double myEventWeight;
515	if(Configuration().IsNoEventWeight()) myEventWeight=1.;
516	else if(event.mc()->weight()!=0.) myEventWeight=event.mc()->weight();
517	else
518	{WARNING << "Found one event with a zero weight. Skipping..." << endmsg;
519	return false; }
520	// apply a 95% efficiency on the weight to account for the event quality filter in the cutflow, which is not implemented otherwise
521	myEventWeight = 0.95*myEventWeight;
522	Manager()->InitializeForNewEvent(myEventWeight);
523
524	vector<const RecLeptonFormat*> electrons, muons, isoElectrons, isoMuons ; // isolated electrons and muons (probably non isolated are not necessary...)
525	// vector<const RecJetFormat*> signalJets, bJets, jetsLargeEta ; // jets, bjets, large eta jets
526	vector<MALorentzVector> signalJets, bJets, jetsLargeEta,unmatched ; // jets, bjets, large eta jets
527	vector<const RecTrackFormat*> electronIsoTracks, muonIsoTracks, isoTracks ; // tracks, tracks identified as electons, and as muons
528	vector<const RecPhotonFormat*> isoPhotons; // isolated photons
529	vector<const MCParticleFormat*> genHSCP,genMuonHSCP, hadrons,bhadrons; // genelectron and muons
530
531
532
533	// **********************************************
534	// ** Objects Selection and Reconstruction **
535	// **********************************************
536	// check if LL HSCP are reconstructed
537	// build HSCP collection
538	for(unsigned int i=0;i<event.mc()->particles().size(); i++)
539	{
540	const MCParticleFormat* prt =&event.mc()->particles()[i];
541	if(std::abs(prt->pdgid()) == 1000024 \|\| std::abs(prt->pdgid()) == 1000011\|\| std::abs(prt->pdgid()) == 1000013\|\| std::abs(prt->pdgid()) == 1000015\|\| std::abs(prt->pdgid()) == 2000011\|\| std::abs(prt->pdgid()) == 2000013\|\| std::abs(prt->pdgid()) == 2000015){
542	genHSCP.push_back(prt);
543	}
544	}
545
546	// Selecting JETS
547
548	for(unsigned int i=0; i<event.rec()->jets().size(); i++)
549	{ const RecJetFormat *thisJet= &(event.rec()->jets()[i]);
550	double pT = thisJet->pt();
551	double eta = thisJet->eta() ;
552	double phi = thisJet->phi();
553	double energy = thisJet->e();
554	MALorentzVector CurrentJet;
555	CurrentJet.SetPtEtaPhiE(pT, eta,phi,energy);
556	if (pT > 30 && fabs(eta) < 5 ) {jetsLargeEta.push_back(CurrentJet) ;} // this extended eta range will be used to calculate MHT
557	if (pT > 30 && fabs(eta) < 2.4 )
558	{
559	signalJets.push_back(CurrentJet);
560	if (thisJet->btag()){
561	bJets.push_back(CurrentJet);
562	}
563	}
564	}
565
566	// Adding Taus back in to jet collection
567	for(unsigned int i=0; i<event.rec()->taus().size(); i++)
568	{ const RecTauFormat *thisTau= &(event.rec()->taus()[i]);
569	double pT = thisTau->pt();
570	double eta = thisTau->eta() ;
571	double phi = thisTau->phi();
572	double energy = thisTau->e();
573	MALorentzVector CurrentTau;
574	CurrentTau.SetPtEtaPhiE(pT, eta,phi,energy);
575	if (pT > 30 && fabs(eta) < 5 ) {jetsLargeEta.push_back(CurrentTau) ;} // this extended eta range will be used to calculate MHT
576	if (pT > 30 && fabs(eta) < 2.4 ){
577	signalJets.push_back(CurrentTau);
578	}
579	}
580	// Selecting muons
581	for(unsigned int i=0; i<event.rec()->muons().size(); i++)
582	{ const RecLeptonFormat *thisMuon = &(event.rec()->muons()[i]);
583	double pt = thisMuon->pt();
584	double eta = thisMuon->eta();
585
586	if (!(pt>10 && fabs(eta) < 2.4)) continue;
587	muons.push_back(thisMuon);
588	double IsoCone = Determine_DeltaR(thisMuon); // selecting the proper Radius
589	double All = PHYSICS->Isol->eflow->sumIsolation(thisMuon,event.rec() , IsoCone ,0.,IsolationEFlow::ALL_COMPONENTS);
590	double Isolation = (All)/pt ;
591
592	if (Isolation < 0.2) {
593	isoMuons.push_back(thisMuon);
594	}
595	}
596	// Selecting long-lived HSCP
597	for(unsigned int i=0;i<genHSCP.size();i++)
598	{
599	const MCParticleFormat *thisHSCP = genHSCP[i];
600	double pt = thisHSCP->pt();
601	double eta = thisHSCP->eta();
602	double beta = thisHSCP->beta();
603	if (!(pt>10 && fabs(eta) < 2.4)) continue;
604	if (thisHSCP->daughters().size()==0){
605	genMuonHSCP.push_back(thisHSCP);
606	}
607	else {
608	double decaylength = beta * 300000000*thisHSCP->daughters()[0]->ctau();
609	if (fabs(decaylength * cos(thisHSCP->theta()))>6){
610	genMuonHSCP.push_back(thisHSCP);
611	}
612	else if (fabs(decaylength * sin(thisHSCP->theta()))>4){
613	genMuonHSCP.push_back(thisHSCP);
614	}
615	}
616	}
617
618	// Selecting electrons
619	for(unsigned int i=0; i<event.rec()->electrons().size(); i++)
620	{ const RecLeptonFormat *thisElectron = &(event.rec()->electrons()[i]);
621	double pt = thisElectron->pt();
622	double eta = thisElectron->eta();
623	if (!(pt>10 && fabs(eta) < 2.5)) continue;
624	electrons.push_back(thisElectron);
625	double IsoCone = Determine_DeltaR(thisElectron); // selecting the proper Radius
626	double All = PHYSICS->Isol->eflow->sumIsolation(thisElectron,event.rec() , IsoCone ,0.,IsolationEFlow::ALL_COMPONENTS);
627	double Isolation = All/pt;
628	if (Isolation < 0.1) { isoElectrons.push_back(thisElectron); }
629	}
630
631
632
633	// Selecting tracks
634	for(unsigned int i=0; i<event.rec()->tracks().size(); i++)
635	{ const RecTrackFormat *thisTrack = &(event.rec()->tracks()[i]);
636	double pt = thisTrack->pt();
637	double eta = thisTrack->eta();
638	double mT = thisTrack->mt_met(event.rec()->MET().momentum());
639	if(!(fabs(eta) < 2.4)) continue; // selecting tracks pT,eta
640
641	int particleId = thisTrack->pdgid() ; // here I identify if the track is an electron or a muon, since this gives different iso requirements
642	bool isElectron = ( particleId== 11 \|\| particleId == -11) ;
643	bool isMuon = ( particleId== 13 \|\| particleId == -13) ;
644	double IsoCone = 0.3 ;
645	double ChargedSum = PHYSICS->Isol->eflow->sumIsolation(thisTrack,event.rec() , IsoCone ,0.,IsolationEFlow::TRACK_COMPONENT) / pt;
646	if ( isElectron && pt > 5 && mT < 100 && ChargedSum < 0.2 ) // selecting isolated electrons tracks
647	{ electronIsoTracks.push_back(thisTrack); }
648	if ( isMuon && pt > 5 && mT < 100 && ChargedSum < 0.2 )
649	{ muonIsoTracks.push_back(thisTrack); }
650	if ( !(isMuon \|\| isElectron) && pt > 10 && mT < 100 && ChargedSum < 0.1)
651	{isoTracks.push_back(thisTrack); }
652
653	}
654
655	// Selecting photons
656	for(unsigned int i=0; i<event.rec()->photons().size(); i++)
657	{ const RecPhotonFormat *thisPhoton = &(event.rec()->photons()[i]);
658	double pt = thisPhoton->pt();
659	double eta = thisPhoton->eta();
660	if (!(pt>100 && fabs(eta) < 2.5)) continue;
661	double IsoCone = 0.3;
662	double All = PHYSICS->Isol->eflow->sumIsolation(thisPhoton,event.rec() , IsoCone ,0.,IsolationEFlow::ALL_COMPONENTS);
663	double Isolation = (All)/pt ;
664	if (Isolation < 0.2) { isoPhotons.push_back(thisPhoton);}
665	}
666
667	// resort the signalJets by pT, now that taus are included
668	sort(signalJets.begin(),signalJets.end(),[](const MALorentzVector A, const MALorentzVector B){return A.Pt()>B.Pt();} );
669
670	// SORTER->sort(signalJets)
671	// SORTER->sort(jetsLargeEta);
672	// SORTER->sort(bJets);
673	// SORTER->sort(isoMuons);
674	// SORTER->sort(isoElectrons);
675	// SORTER->sort(isoPhotons);
676
677
678	// signalJets = Removal(signalJets, isoMuons, 0.2);
679	// signalJets = Removal(signalJets, isoElectrons, 0.2);
680
681
682	// *******************************************************
683	// Preselection Cuts
684	// *******************************************************
685
686	// Defining the kinematic variables
687	int nb = bJets.size();
688	int nJets = signalJets.size();
689	// Calculate the HT
690	double HT = 0;
691
692	for(unsigned int j = 0; j < signalJets.size(); j++)
693	{
694	HT += signalJets[j].Pt();
695	}
696	// Calculate the MHT
697	double MHT = 0;
698	MALorentzVector VectorMHT = MALorentzVector();
699	// MALorentzVector CurrentJet;
700	for(unsigned int j = 0; j < jetsLargeEta.size(); j++)
701	{
702	// CurrentJet.SetPtEtaPhiE(jetsLargeEta[j].Pt(), jetsLargeEta[j].Eta(),jetsLargeEta[j]->Phi(),jetsLargeEta[j]->E());
703
704	// VectorMHT -= CurrentJet;
705	VectorMHT -= jetsLargeEta[j];
706	}
707
708	MHT = VectorMHT.Pt();
709
710
711
712	double DeltaPhi_j1MHT , DeltaPhi_j2MHT , DeltaPhi_j3MHT , DeltaPhi_j4MHT ;
713
714	// Calculating the DeltaPhi.
715	// Since the preselection requires certain vales of DPhi based on different number of jets,
716	// the DPhi for j3 and j4 assume the required values for passing the cut, so that if in the event there are less than 3 and 4 jets, the cut on the 3rd and 4th jets has no effect
717	if ( nJets >= 1 ) { DeltaPhi_j1MHT = fabs(VectorMHT.DeltaPhi(signalJets[0])); }
718	if ( nJets < 1 ) { DeltaPhi_j1MHT = 0.5 ; }
719
720	if ( nJets >= 2 ) { DeltaPhi_j2MHT = fabs(VectorMHT.DeltaPhi(signalJets[1])); }
721	if ( nJets < 2 ) { DeltaPhi_j2MHT = 0.5 ; }
722
723	if ( nJets >= 3 ) { DeltaPhi_j3MHT = fabs(VectorMHT.DeltaPhi(signalJets[2])); }
724	if ( nJets < 3 ) { DeltaPhi_j3MHT = 0.3 ; }
725
726	if ( nJets >= 4 ) { DeltaPhi_j4MHT = fabs(VectorMHT.DeltaPhi(signalJets[3])); }
727	if ( nJets < 4 ) { DeltaPhi_j4MHT = 0.3 ; }
728
729
730	bool DPhi_1 = DeltaPhi_j1MHT >= 0.5 ;
731	bool DPhi_2 = DeltaPhi_j2MHT >= 0.5 ;
732	bool DPhi_3 = DeltaPhi_j3MHT >= 0.3 ;
733	bool DPhi_4 = DeltaPhi_j4MHT >= 0.3 ;
734	bool DPhi = DPhi_1 && DPhi_2 && DPhi_3 && DPhi_4 ; // Total condition on DPhi angles
735
736	// ******************************************************
737	// Filling the histograms (for validation) *
738	// ******************************************************
739	// NEED TO ADAPT THIS TO OUR VALIDATION HISTOGRAMS
740	int isoE = isoElectrons.size(), isoM = isoMuons.size()+genMuonHSCP.size() , isoEtracks = electronIsoTracks.size(), isoMtracks = muonIsoTracks.size() , isoTr = isoTracks.size() , isoP = isoPhotons.size();
741	bool Isolation = (isoE == 0 && isoM==0 && isoEtracks == 0 && isoMtracks == 0 && isoTr == 0 && isoP ==0) ;
742	bool passes = (Isolation && HT>300 && MHT >300 && DPhi && HT>MHT && nJets>=2);
743	// histograms excluding the cut on the variable shown
744	if (Isolation && HT > 300 && MHT > 300 && DPhi && HT>MHT) {
745	Manager()->FillHisto("nJets_Pre", nJets ) ;}
746	if (Isolation && nJets >=2 && HT > 300 && MHT >300 && DPhi && HT>MHT) { Manager()->FillHisto("bJets_Pre", nb ) ;}
747	if (Isolation && nJets >=2 && MHT > 300 && DPhi) { Manager()->FillHisto("HT_Pre", HT ) ;}
748	if (Isolation && nJets >=2 && HT > 300 && DPhi) { Manager()->FillHisto("MHT_Pre", MHT ) ;}
749	// histograms post preselection cuts
750	if (passes) {
751	Manager()->FillHisto("nJets_Post", nJets );
752	Manager()->FillHisto("bJets_Post", nb );
753	Manager()->FillHisto("HT_Post", HT );
754	Manager()->FillHisto("MHT_Post", MHT );
755
756	}
757
758
759	// ******************************************************
760	// Applying the preselection cuts *
761	// ******************************************************
762
763	// 1) Cut on NJet >= 2
764	if(!Manager()->ApplyCut((nJets >= 2),"Njet>=2")) return true;
765	// 2) Cut on HT > 300
766	if(!Manager()->ApplyCut((HT>300), "HT>300")) return true;
767	// 3) Cut on MHT > 300
768	if(!Manager()->ApplyCut((MHT>300), "MHT>300")) return true;
769	// 4) Cut on HT > MHT
770	if(!Manager()->ApplyCut((HT>MHT), "HT>MHT")) return true;
771	// 5) No isolated muons
772	if(!Manager()->ApplyCut((isoM == 0),"NoIsoMuons")) return true;
773	// 6) No isolated muons tracks
774	if(!Manager()->ApplyCut((isoMtracks == 0),"NoMuonsTracks")) return true;
775	// 7) No isolated electrons
776	if(!Manager()->ApplyCut((isoE == 0),"NoIsoElectrons")) return true;
777	// 8) No electrons tracks
778	if(!Manager()->ApplyCut((isoEtracks == 0),"NoElectronsTracks")) return true;
779	// 9) No hadronic tracks
780	if(!Manager()->ApplyCut((isoTr == 0),"NoIsoTracks")) return true;
781	// 10) No isolated photons
782	if(!Manager()->ApplyCut((isoP == 0),"NoIsoPhotons")) return true;
783	// 11) DeltaPhi between jets and MHT
784
785	if(!Manager()->ApplyCut(DPhi_1 , "DPhi(MHTj1)>0.5") ) return true;
786	if(!Manager()->ApplyCut(DPhi_2 , "DPhi(MHTj2)>0.5") ) return true;
787	if(!Manager()->ApplyCut(DPhi_3 , "DPhi(MHTj3)>0.3") ) return true;
788	if(!Manager()->ApplyCut(DPhi_4 , "DPhi(MHTj4)>0.3") ) return true;
789
790
791	// *************************************************
792	// Specific cuts for the aggregated Signal Regions *
793	// *************************************************
794
795	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 2 && nb == 0, "AGGSR1_Njet2_Nb0_HT600_MHT600" ) ) return true;
796	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 3 && nb == 0, "AGGSR2_Njet4_Nb0_HT1700_MHT850" ) ) return true;
797	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 6 && nb == 0, "AGGSR3_Njet6_Nb0_HT600_MHT_600" ) ) return true;
798	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 8 && nb <=1, "AGGSR4_Njet8_Nb01_HT600_MHT600" ) ) return true;
799	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 10 && nb <= 1, "AGGSR5_Njet10_Nb01_HT1700_MHT850" ) ) return true;
800	if(!Manager()->ApplyCut(HT > 300 && MHT > 300 && nJets >= 4 && nb >= 2, "AGGSR6_Njet4_Nb2_HT300_MHT300" ) ) return true;
801	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 2 && nb >= 2, "AGGSR7_Njet2_Nb2_HT600_MHT600" ) ) return true;
802	if(!Manager()->ApplyCut(HT > 350 && MHT > 350 && nJets >= 6 && nb >= 2, "AGGSR8_Njet6_Nb2_HT350_MHT350" ) ) return true;
803	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 4 && nb >= 2, "AGGSR9_NJet4_Nb2_HT600_MHT600" ) ) return true;
804	if(!Manager()->ApplyCut(HT > 300 && MHT > 300 && nJets >= 8 && nb >= 3, "AGGSR10_Njet8_Nb3_HT300_MHT300" ) ) return true;
805	if(!Manager()->ApplyCut(HT > 600 && MHT > 600 && nJets >= 6 && nb >= 1, "AGGSR11_Njet6_Nb1_HT600_MHT600" ) ) return true;
806	if(!Manager()->ApplyCut(HT > 850 && MHT > 850 && nJets >= 10 && nb >= 3, "AGGSR12_Njet10_Nb3_HT850_MHT850" ) ) return true;
807
808
809	// *************************************************
810	// Specific cuts for the 174 Signal Regions *
811	// *************************************************
812
813	//2-3 jets, 0bjets
814	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR1_Njet23_Nb0_HT300600_MHT300350")) return true;
815	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR2_Njet23_Nb0_HT6001200_MHT300350")) return true;
816	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR3_Njet23_Nb0_HT1200_MHT300350")) return true;
817	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR4_Njet23_Nb0_HT350600_MHT350600")) return true;
818	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR5_Njet23_Nb0_HT6001200_MHT350600")) return true;
819	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR6_Njet23_Nb0_HT1200_MHT350600")) return true;
820	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR7_Njet23_Nb0_HT6001200_MHT600850")) return true;
821	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR8_Njet23_Nb0_HT1200_MHT600850")) return true;
822	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR9_Njet23_Nb0_HT8501700_MHT850")) return true;
823	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb == 0,"SR10_Njet23_Nb0_HT1700_MHT850")) return true;
824	// 2-3 jets, 1 bjet
825	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR11_Njet23_Nb1_HT300600_MHT300350")) return true;
826	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR12_Njet23_Nb1_HT6001200_MHT300350")) return true;
827	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR13_Njet23_Nb1_HT1200_MHT300350")) return true;
828	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR14_Njet23_Nb1_HT350600_MHT350600")) return true;
829	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR15_Njet23_Nb1_HT6001200_MHT350600")) return true;
830	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR16_Njet23_Nb1_HT1200_MHT350600")) return true;
831	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR17_Njet23_Nb1_HT6001200_MHT600850")) return true;
832	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR18_Njet23_Nb1_HT1200_MHT600850")) return true;
833	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR19_Njet23_Nb1_HT8501700_MHT850")) return true;
834	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb == 1,"SR20_Njet23_Nb1_HT1700_MHT850")) return true;
835	// 2-3 jets, >=2 bjets
836	if(!Manager()->ApplyCut(HT > 300 && HT <= 600&& MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR21_Njet23_Nb2_HT300600_MHT300350")) return true;
837	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR22_Njet23_Nb2_HT6001200_MHT300350")) return true;
838	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR23_Njet23_Nb2_HT1200_MHT300350")) return true;
839	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR24_Njet23_Nb2_HT350600_MHT350600")) return true;
840	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR25_Njet23_Nb2_HT6001200_MHT350600")) return true;
841	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR26_Njet23_Nb2_HT1200_MHT350600")) return true;
842	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR27_Njet23_Nb2_HT6001200_MHT600850")) return true;
843	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR28_Njet23_Nb2_HT1200_MHT600850")) return true;
844	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR29_Njet23_Nb2_HT8501700_MHT850")) return true;
845	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 2 \|\| nJets == 3) && nb >= 2,"SR30_Njet23_Nb2_HT1700_MHT850")) return true;
846
847	//4-5 jets, 0bjets
848	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR31_Njet45_Nb0_HT300600_MHT300350")) return true;
849	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR32_Njet45_Nb0_HT6001200_MHT300350")) return true;
850	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR33_Njet45_Nb0_HT1200_MHT300350")) return true;
851	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR34_Njet45_Nb0_HT350600_MHT350600")) return true;
852	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR35_Njet45_Nb0_HT6001200_MHT350600")) return true;
853	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR36_Njet45_Nb0_HT1200_MHT350600")) return true;
854	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR37_Njet45_Nb0_HT6001200_MHT600850")) return true;
855	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR38_Njet45_Nb0_HT1200_MHT600850")) return true;
856	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR39_Njet45_Nb0_HT8501700_MHT850")) return true;
857	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 0,"SR40_Njet45_Nb0_HT1700_MHT850")) return true;
858	//4-5 jets, 1bjets
859	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR41_Njet45_Nb1_HT300600_MHT300350")) return true;
860	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR42_Njet45_Nb1_HT6001200_MHT300350")) return true;
861	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR43_Njet45_Nb1_HT1200_MHT300350")) return true;
862	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR44_Njet45_Nb1_HT350600_MHT350600")) return true;
863	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR45_Njet45_Nb1_HT6001200_MHT350600")) return true;
864	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR46_Njet45_Nb1_HT1200_MHT350600")) return true;
865	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR47_Njet45_Nb1_HT6001200_MHT600850")) return true;
866	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR48_Njet45_Nb1_HT1200_MHT600850")) return true;
867	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR49_Njet45_Nb1_HT8501700_MHT850")) return true;
868	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 1,"SR50_Njet45_Nb1_HT1700_MHT850")) return true;
869	//4-5 jets, 2bjets
870	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR51_Njet45_Nb2_HT300600_MHT300350")) return true;
871	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR52_Njet45_Nb2_HT6001200_MHT300350")) return true;
872	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR53_Njet45_Nb2_HT1200_MHT300350")) return true;
873	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR54_Njet45_Nb2_HT350600_MHT350600")) return true;
874	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR55_Njet45_Nb2_HT6001200_MHT350600")) return true;
875	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR56_Njet45_Nb2_HT1200_MHT350600")) return true;
876	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR57_Njet45_Nb2_HT6001200_MHT600850")) return true;
877	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR58_Njet45_Nb2_HT1200_MHT600850")) return true;
878	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR59_Njet45_Nb2_HT8501700_MHT850")) return true;
879	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb == 2,"SR60_Njet45_Nb2_HT1700_MHT850")) return true;
880	//4-5 jets, 3bjets
881	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR61_Njet45_Nb3_HT300600_MHT300350")) return true;
882	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR62_Njet45_Nb3_HT6001200_MHT300350")) return true;
883	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR63_Njet45_Nb3_HT1200_MHT300350")) return true;
884	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR64_Njet45_Nb3_HT350600_MHT350600")) return true;
885	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR65_Njet45_Nb3_HT6001200_MHT350600")) return true;
886	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR66_Njet45_Nb3_HT1200_MHT350600")) return true;
887	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR67_Njet45_Nb3_HT6001200_MHT600850")) return true;
888	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR68_Njet45_Nb3_HT1200_MHT600850")) return true;
889	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR69_Njet45_Nb3_HT8501700_MHT850")) return true;
890	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 4 \|\| nJets == 5) && nb >= 3,"SR70_Njet45_Nb3_HT1700_MHT850")) return true;
891
892
893	//6-7 jets, 0bjets
894	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR71_Njet67_Nb0_HT300600_MHT300350")) return true;
895	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR72_Njet67_Nb0_HT6001200_MHT300350")) return true;
896	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR73_Njet67_Nb0_HT1200_MHT300350")) return true;
897	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR74_Njet67_Nb0_HT350600_MHT350600")) return true;
898	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR75_Njet67_Nb0_HT6001200_MHT350600")) return true;
899	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR76_Njet67_Nb0_HT1200_MHT350600")) return true;
900	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR77_Njet67_Nb0_HT6001200_MHT600850")) return true;
901	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR78_Njet67_Nb0_HT1200_MHT600850")) return true;
902	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR79_Njet67_Nb0_HT8501700_MHT850")) return true;
903	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 0,"SR80_Njet67_Nb0_HT1700_MHT850")) return true;
904	//6-7 jets, 1bjets
905	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR81_Njet67_Nb1_HT300600_MHT300350")) return true;
906	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR82_Njet67_Nb1_HT6001200_MHT300350")) return true;
907	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR83_Njet67_Nb1_HT1200_MHT300350")) return true;
908	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR84_Njet67_Nb1_HT350600_MHT350600")) return true;
909	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR85_Njet67_Nb1_HT6001200_MHT350600")) return true;
910	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR86_Njet67_Nb1_HT1200_MHT350600")) return true;
911	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR87_Njet67_Nb1_HT6001200_MHT600850")) return true;
912	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR88_Njet67_Nb1_HT1200_MHT600850")) return true;
913	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR89_Njet67_Nb1_HT8501700_MHT850")) return true;
914	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 1,"SR90_Njet67_Nb1_HT1700_MHT850")) return true;
915	//6-7 jets, 2bjets
916	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR91_Njet67_Nb2_HT300600_MHT300350")) return true;
917	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR92_Njet67_Nb2_HT6001200_MHT300350")) return true;
918	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR93_Njet67_Nb2_HT1200_MHT300350")) return true;
919	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR94_Njet67_Nb2_HT350600_MHT350600")) return true;
920	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR95_Njet67_Nb2_HT6001200_MHT350600")) return true;
921	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR96_Njet67_Nb2_HT1200_MHT350600")) return true;
922	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR97_Njet67_Nb2_HT6001200_MHT600850")) return true;
923	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR98_Njet67_Nb2_HT1200_MHT600850")) return true;
924	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR99_Njet67_Nb2_HT8501700_MHT850")) return true;
925	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb == 2,"SR100_Njet67_Nb2_HT1700_MHT850")) return true;
926	//6-7 jets, 3bjets
927	if(!Manager()->ApplyCut(HT > 300 && HT <= 600 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR101_Njet67_Nb3_HT300600_MHT300350")) return true;
928	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR102_Njet67_Nb3_HT6001200_MHT300350")) return true;
929	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR103_Njet67_Nb3_HT1200_MHT300350")) return true;
930	if(!Manager()->ApplyCut(HT > 350 && HT <= 600 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR104_Njet67_Nb3_HT350600_MHT350600")) return true;
931	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR105_Njet67_Nb3_HT6001200_MHT350600")) return true;
932	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR106_Njet67_Nb3_HT1200_MHT350600")) return true;
933	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR107_Njet67_Nb3_HT6001200_MHT600850")) return true;
934	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR108_Njet67_Nb3_HT1200_MHT600850")) return true;
935	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR109_Njet67_Nb3_HT8501700_MHT850")) return true;
936	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 6 \|\| nJets == 7) && nb >= 3,"SR110_Njet67_Nb3_HT1700_MHT850")) return true;
937
938	//8-9 jets, 0bjets
939	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR111_Njet89_Nb0_HT6001200_MHT300350")) return true;
940	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR112_Njet89_Nb0_HT1200_MHT300350")) return true;
941	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR113_Njet89_Nb0_HT6001200_MHT350600")) return true;
942	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR114_Njet89_Nb0_HT1200_MHT350600")) return true;
943	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR115_Njet89_Nb0_HT6001200_MHT600850")) return true;
944	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR116_Njet89_Nb0_HT1200_MHT600850")) return true;
945	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR117_Njet89_Nb0_HT8501700_MHT850")) return true;
946	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 0,"SR118_Njet89_Nb0_HT1700_MHT850")) return true;
947	//8-9 jets, 1bjets
948	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR119_Njet89_Nb1_HT6001200_MHT300350")) return true;
949	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR120_Njet89_Nb1_HT1200_MHT300350")) return true;
950	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR121_Njet89_Nb1_HT6001200_MHT350600")) return true;
951	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR122_Njet89_Nb1_HT1200_MHT350600")) return true;
952	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR123_Njet89_Nb1_HT6001200_MHT600850")) return true;
953	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR124_Njet89_Nb1_HT1200_MHT600850")) return true;
954	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR125_Njet89_Nb1_HT8501700_MHT850")) return true;
955	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 1,"SR126_Njet89_Nb1_HT1700_MHT850")) return true;
956	//8-9 jets, 2bjets
957	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR127_Njet89_Nb2_HT6001200_MHT300350")) return true;
958	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR128_Njet89_Nb2_HT1200_MHT300350")) return true;
959	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR129_Njet89_Nb2_HT6001200_MHT350600")) return true;
960	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR130_Njet89_Nb2_HT1200_MHT350600")) return true;
961	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR131_Njet89_Nb2_HT6001200_MHT600850")) return true;
962	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR132_Njet89_Nb2_HT1200_MHT600850")) return true;
963	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR133_Njet89_Nb2_HT8501700_MHT850")) return true;
964	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb == 2,"SR134_Njet89_Nb2_HT1700_MHT850")) return true;
965	//8-9 jets, 3bjets
966	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR135_Njet89_Nb3_HT6001200_MHT300350")) return true;
967	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR136_Njet89_Nb3_HT1200_MHT300350")) return true;
968	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR137_Njet89_Nb3_HT6001200_MHT350600")) return true;
969	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR138_Njet89_Nb3_HT1200_MHT350600")) return true;
970	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR139_Njet89_Nb3_HT6001200_MHT600850")) return true;
971	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR140_Njet89_Nb3_HT1200_MHT600850")) return true;
972	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR141_Njet89_Nb3_HT8501700_MHT850")) return true;
973	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && (nJets == 8 \|\| nJets == 9) && nb >= 3,"SR142_Njet89_Nb3_HT1700_MHT850")) return true;
974
975	//>10 jets, 0bjets
976	if(!Manager()->ApplyCut(HT > 600 && HT <=1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 0,"SR143_Njet10_Nb0_HT6001200_MHT300350")) return true;
977	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 0,"SR144_Njet10_Nb0_HT1200_MHT300350")) return true;
978	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 0,"SR145_Njet10_Nb0_HT6001200_MHT350600")) return true;
979	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 0,"SR146_Njet10_Nb0_HT1200_MHT350600")) return true;
980	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 0,"SR147_Njet10_Nb0_HT6001200_MHT600850")) return true;
981	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 0,"SR148_Njet10_Nb0_HT1200_MHT600850")) return true;
982	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && nJets >= 10 && nb == 0,"SR149_Njet10_Nb0_HT8501700_MHT850")) return true;
983	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 10 && nb == 0,"SR150_Njet10_Nb0_HT1700_MHT850")) return true;
984	//>10 jets, 1bjets
985	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 1,"SR151_Njet10_Nb1_HT6001200_MHT300350")) return true;
986	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 1,"SR152_Njet10_Nb1_HT1200_MHT300350")) return true;
987	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 1,"SR153_Njet10_Nb1_HT6001200_MHT350600")) return true;
988	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 1,"SR154_Njet10_Nb1_HT1200_MHT350600")) return true;
989	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 1,"SR155_Njet10_Nb1_HT6001200_MHT600850")) return true;
990	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 1,"SR156_Njet10_Nb1_HT1200_MHT600850")) return true;
991	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && nJets >= 10 && nb == 1,"SR157_Njet10_Nb1_HT8501700_MHT850")) return true;
992	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 10 && nb == 1,"SR158_Njet10_Nb1_HT1700_MHT850")) return true;
993	//>10 jets, 2bjets
994	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 2,"SR159_Njet10_Nb2_HT6001200_MHT300350")) return true;
995	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb == 2,"SR160_Njet10_Nb2_HT1200_MHT300350")) return true;
996	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 2,"SR161_Njet10_Nb2_HT6001200_MHT350600")) return true;
997	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb == 2,"SR162_Njet10_Nb2_HT1200_MHT350600")) return true;
998	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 2,"SR163_Njet10_Nb2_HT6001200_MHT600850")) return true;
999	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb == 2,"SR164_Njet10_Nb2_HT1200_MHT600850")) return true;
1000	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && nJets >= 10 && nb == 2,"SR165_Njet10_Nb2_HT8501700_MHT850")) return true;
1001	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 10 && nb == 2,"SR166_Njet10_Nb2_HT1700_MHT850")) return true;
1002	//>10 jets, 3bjets
1003	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb >= 3,"SR167_Njet10_Nb3_HT6001200_MHT300350")) return true;
1004	if(!Manager()->ApplyCut(HT > 1200 && MHT > 300 && MHT <= 350 && nJets >= 10 && nb >= 3,"SR168_Njet10_Nb3_HT1200_MHT300350")) return true;
1005	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb >= 3,"SR169_Njet10_Nb3_HT6001200_MHT350600")) return true;
1006	if(!Manager()->ApplyCut(HT > 1200 && MHT > 350 && MHT <= 600 && nJets >= 10 && nb >= 3,"SR170_Njet10_Nb3_HT1200_MHT350600")) return true;
1007	if(!Manager()->ApplyCut(HT > 600 && HT <= 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb >= 3,"SR171_Njet10_Nb3_HT6001200_MHT600850")) return true;
1008	if(!Manager()->ApplyCut(HT > 1200 && MHT > 600 && MHT <= 850 && nJets >= 10 && nb >= 3,"SR172_Njet10_Nb3_HT1200_MHT600850")) return true;
1009	if(!Manager()->ApplyCut(HT > 850 && HT <= 1700 && MHT > 850 && nJets >= 10 && nb >= 3,"SR173_Njet10_Nb3_HT8501700_MHT850")) return true;
1010	if(!Manager()->ApplyCut(HT > 1700 && MHT > 850 && nJets >= 10 && nb >= 3,"SR174_Njet10_Nb3_HT1700_MHT850")) return true;
1011
1012
1013	}
1014
1015	return true;
1016	}
1017
1018
1019
1020
1021	// This function returns the DeltaR depending on the pT of the lepton considered,
1022	// required for selecting isolated leptons
1023	double cms_sus_19_006::Determine_DeltaR(const RecLeptonFormat* thisLepton)
1024	{ double DeltaR=0.;
1025	double pT = thisLepton->pt();
1026	if (pT <= 50.0 ) { DeltaR = 0.2 ; }
1027	else if (pT > 50.0 && pT <= 200.0 ) { DeltaR = 10.0/pT ; }
1028	else if (pT > 200.0 ) { DeltaR = 0.05 ; }
1029	return DeltaR;
1030	}
1031

Download in other formats:

Original Format