| 1 | #include "SampleAnalyzer/User/Analyzer/sfs_cms_sus_16_048.h"
|
|---|
| 2 | using namespace MA5;
|
|---|
| 3 | using namespace std;
|
|---|
| 4 |
|
|---|
| 5 | // -----------------------------------------------------------------------------
|
|---|
| 6 | // Initialize
|
|---|
| 7 | // function called one time at the beginning of the analysis
|
|---|
| 8 | // -----------------------------------------------------------------------------
|
|---|
| 9 | bool sfs_cms_sus_16_048::Initialize(const MA5::Configuration& cfg,
|
|---|
| 10 | const std::map<std::string,std::string>& parameters)
|
|---|
| 11 | {
|
|---|
| 12 | INFO << " <><><><><><><><><><><><><><><><><><><><><><><><><>" << endmsg;
|
|---|
| 13 | INFO << " <> Analysis: CMS_SUS_16_048, arXiv:1801.01846 <>" << endmsg;
|
|---|
| 14 | INFO << " <> (soft dilepton LHC13, 35.9 fb^-1)<>" << endmsg;
|
|---|
| 15 | INFO << " <> Recasted by: B. Fuks <>" << endmsg;
|
|---|
| 16 | INFO << " <> Contacts: fuks@lpthe.jussieu.fr <>" << endmsg;
|
|---|
| 17 | INFO << " <> Based on MadAnalysis 5 v1.8 <>" << endmsg;
|
|---|
| 18 | INFO << " <><><><><><><><><><><><><><><><><><><><><><><><><>" << endmsg;
|
|---|
| 19 |
|
|---|
| 20 | // Initializing PhysicsService for MC
|
|---|
| 21 | PHYSICS->mcConfig().Reset();
|
|---|
| 22 |
|
|---|
| 23 |
|
|---|
| 24 | // definition of the multiparticle "hadronic"
|
|---|
| 25 | PHYSICS->mcConfig().AddHadronicId(-20543);
|
|---|
| 26 | PHYSICS->mcConfig().AddHadronicId(-20533);
|
|---|
| 27 | PHYSICS->mcConfig().AddHadronicId(-20523);
|
|---|
| 28 | PHYSICS->mcConfig().AddHadronicId(-20513);
|
|---|
| 29 | PHYSICS->mcConfig().AddHadronicId(-20433);
|
|---|
| 30 | PHYSICS->mcConfig().AddHadronicId(-20423);
|
|---|
| 31 | PHYSICS->mcConfig().AddHadronicId(-20413);
|
|---|
| 32 | PHYSICS->mcConfig().AddHadronicId(-20323);
|
|---|
| 33 | PHYSICS->mcConfig().AddHadronicId(-20313);
|
|---|
| 34 | PHYSICS->mcConfig().AddHadronicId(-20213);
|
|---|
| 35 | PHYSICS->mcConfig().AddHadronicId(-10543);
|
|---|
| 36 | PHYSICS->mcConfig().AddHadronicId(-10541);
|
|---|
| 37 | PHYSICS->mcConfig().AddHadronicId(-10533);
|
|---|
| 38 | PHYSICS->mcConfig().AddHadronicId(-10531);
|
|---|
| 39 | PHYSICS->mcConfig().AddHadronicId(-10523);
|
|---|
| 40 | PHYSICS->mcConfig().AddHadronicId(-10521);
|
|---|
| 41 | PHYSICS->mcConfig().AddHadronicId(-10513);
|
|---|
| 42 | PHYSICS->mcConfig().AddHadronicId(-10511);
|
|---|
| 43 | PHYSICS->mcConfig().AddHadronicId(-10433);
|
|---|
| 44 | PHYSICS->mcConfig().AddHadronicId(-10431);
|
|---|
| 45 | PHYSICS->mcConfig().AddHadronicId(-10423);
|
|---|
| 46 | PHYSICS->mcConfig().AddHadronicId(-10421);
|
|---|
| 47 | PHYSICS->mcConfig().AddHadronicId(-10413);
|
|---|
| 48 | PHYSICS->mcConfig().AddHadronicId(-10411);
|
|---|
| 49 | PHYSICS->mcConfig().AddHadronicId(-10323);
|
|---|
| 50 | PHYSICS->mcConfig().AddHadronicId(-10321);
|
|---|
| 51 | PHYSICS->mcConfig().AddHadronicId(-10313);
|
|---|
| 52 | PHYSICS->mcConfig().AddHadronicId(-10311);
|
|---|
| 53 | PHYSICS->mcConfig().AddHadronicId(-10213);
|
|---|
| 54 | PHYSICS->mcConfig().AddHadronicId(-10211);
|
|---|
| 55 | PHYSICS->mcConfig().AddHadronicId(-5554);
|
|---|
| 56 | PHYSICS->mcConfig().AddHadronicId(-5544);
|
|---|
| 57 | PHYSICS->mcConfig().AddHadronicId(-5542);
|
|---|
| 58 | PHYSICS->mcConfig().AddHadronicId(-5534);
|
|---|
| 59 | PHYSICS->mcConfig().AddHadronicId(-5532);
|
|---|
| 60 | PHYSICS->mcConfig().AddHadronicId(-5524);
|
|---|
| 61 | PHYSICS->mcConfig().AddHadronicId(-5522);
|
|---|
| 62 | PHYSICS->mcConfig().AddHadronicId(-5514);
|
|---|
| 63 | PHYSICS->mcConfig().AddHadronicId(-5512);
|
|---|
| 64 | PHYSICS->mcConfig().AddHadronicId(-5503);
|
|---|
| 65 | PHYSICS->mcConfig().AddHadronicId(-5444);
|
|---|
| 66 | PHYSICS->mcConfig().AddHadronicId(-5442);
|
|---|
| 67 | PHYSICS->mcConfig().AddHadronicId(-5434);
|
|---|
| 68 | PHYSICS->mcConfig().AddHadronicId(-5432);
|
|---|
| 69 | PHYSICS->mcConfig().AddHadronicId(-5424);
|
|---|
| 70 | PHYSICS->mcConfig().AddHadronicId(-5422);
|
|---|
| 71 | PHYSICS->mcConfig().AddHadronicId(-5414);
|
|---|
| 72 | PHYSICS->mcConfig().AddHadronicId(-5412);
|
|---|
| 73 | PHYSICS->mcConfig().AddHadronicId(-5403);
|
|---|
| 74 | PHYSICS->mcConfig().AddHadronicId(-5401);
|
|---|
| 75 | PHYSICS->mcConfig().AddHadronicId(-5342);
|
|---|
| 76 | PHYSICS->mcConfig().AddHadronicId(-5334);
|
|---|
| 77 | PHYSICS->mcConfig().AddHadronicId(-5332);
|
|---|
| 78 | PHYSICS->mcConfig().AddHadronicId(-5324);
|
|---|
| 79 | PHYSICS->mcConfig().AddHadronicId(-5322);
|
|---|
| 80 | PHYSICS->mcConfig().AddHadronicId(-5314);
|
|---|
| 81 | PHYSICS->mcConfig().AddHadronicId(-5312);
|
|---|
| 82 | PHYSICS->mcConfig().AddHadronicId(-5303);
|
|---|
| 83 | PHYSICS->mcConfig().AddHadronicId(-5301);
|
|---|
| 84 | PHYSICS->mcConfig().AddHadronicId(-5242);
|
|---|
| 85 | PHYSICS->mcConfig().AddHadronicId(-5232);
|
|---|
| 86 | PHYSICS->mcConfig().AddHadronicId(-5224);
|
|---|
| 87 | PHYSICS->mcConfig().AddHadronicId(-5222);
|
|---|
| 88 | PHYSICS->mcConfig().AddHadronicId(-5214);
|
|---|
| 89 | PHYSICS->mcConfig().AddHadronicId(-5212);
|
|---|
| 90 | PHYSICS->mcConfig().AddHadronicId(-5203);
|
|---|
| 91 | PHYSICS->mcConfig().AddHadronicId(-5201);
|
|---|
| 92 | PHYSICS->mcConfig().AddHadronicId(-5142);
|
|---|
| 93 | PHYSICS->mcConfig().AddHadronicId(-5132);
|
|---|
| 94 | PHYSICS->mcConfig().AddHadronicId(-5122);
|
|---|
| 95 | PHYSICS->mcConfig().AddHadronicId(-5114);
|
|---|
| 96 | PHYSICS->mcConfig().AddHadronicId(-5112);
|
|---|
| 97 | PHYSICS->mcConfig().AddHadronicId(-5103);
|
|---|
| 98 | PHYSICS->mcConfig().AddHadronicId(-5101);
|
|---|
| 99 | PHYSICS->mcConfig().AddHadronicId(-4444);
|
|---|
| 100 | PHYSICS->mcConfig().AddHadronicId(-4434);
|
|---|
| 101 | PHYSICS->mcConfig().AddHadronicId(-4432);
|
|---|
| 102 | PHYSICS->mcConfig().AddHadronicId(-4424);
|
|---|
| 103 | PHYSICS->mcConfig().AddHadronicId(-4422);
|
|---|
| 104 | PHYSICS->mcConfig().AddHadronicId(-4414);
|
|---|
| 105 | PHYSICS->mcConfig().AddHadronicId(-4412);
|
|---|
| 106 | PHYSICS->mcConfig().AddHadronicId(-4403);
|
|---|
| 107 | PHYSICS->mcConfig().AddHadronicId(-4334);
|
|---|
| 108 | PHYSICS->mcConfig().AddHadronicId(-4332);
|
|---|
| 109 | PHYSICS->mcConfig().AddHadronicId(-4324);
|
|---|
| 110 | PHYSICS->mcConfig().AddHadronicId(-4322);
|
|---|
| 111 | PHYSICS->mcConfig().AddHadronicId(-4314);
|
|---|
| 112 | PHYSICS->mcConfig().AddHadronicId(-4312);
|
|---|
| 113 | PHYSICS->mcConfig().AddHadronicId(-4303);
|
|---|
| 114 | PHYSICS->mcConfig().AddHadronicId(-4301);
|
|---|
| 115 | PHYSICS->mcConfig().AddHadronicId(-4232);
|
|---|
| 116 | PHYSICS->mcConfig().AddHadronicId(-4224);
|
|---|
| 117 | PHYSICS->mcConfig().AddHadronicId(-4222);
|
|---|
| 118 | PHYSICS->mcConfig().AddHadronicId(-4214);
|
|---|
| 119 | PHYSICS->mcConfig().AddHadronicId(-4212);
|
|---|
| 120 | PHYSICS->mcConfig().AddHadronicId(-4203);
|
|---|
| 121 | PHYSICS->mcConfig().AddHadronicId(-4201);
|
|---|
| 122 | PHYSICS->mcConfig().AddHadronicId(-4132);
|
|---|
| 123 | PHYSICS->mcConfig().AddHadronicId(-4122);
|
|---|
| 124 | PHYSICS->mcConfig().AddHadronicId(-4114);
|
|---|
| 125 | PHYSICS->mcConfig().AddHadronicId(-4112);
|
|---|
| 126 | PHYSICS->mcConfig().AddHadronicId(-4103);
|
|---|
| 127 | PHYSICS->mcConfig().AddHadronicId(-4101);
|
|---|
| 128 | PHYSICS->mcConfig().AddHadronicId(-3334);
|
|---|
| 129 | PHYSICS->mcConfig().AddHadronicId(-3324);
|
|---|
| 130 | PHYSICS->mcConfig().AddHadronicId(-3322);
|
|---|
| 131 | PHYSICS->mcConfig().AddHadronicId(-3314);
|
|---|
| 132 | PHYSICS->mcConfig().AddHadronicId(-3312);
|
|---|
| 133 | PHYSICS->mcConfig().AddHadronicId(-3303);
|
|---|
| 134 | PHYSICS->mcConfig().AddHadronicId(-3224);
|
|---|
| 135 | PHYSICS->mcConfig().AddHadronicId(-3222);
|
|---|
| 136 | PHYSICS->mcConfig().AddHadronicId(-3214);
|
|---|
| 137 | PHYSICS->mcConfig().AddHadronicId(-3212);
|
|---|
| 138 | PHYSICS->mcConfig().AddHadronicId(-3203);
|
|---|
| 139 | PHYSICS->mcConfig().AddHadronicId(-3201);
|
|---|
| 140 | PHYSICS->mcConfig().AddHadronicId(-3122);
|
|---|
| 141 | PHYSICS->mcConfig().AddHadronicId(-3114);
|
|---|
| 142 | PHYSICS->mcConfig().AddHadronicId(-3112);
|
|---|
| 143 | PHYSICS->mcConfig().AddHadronicId(-3103);
|
|---|
| 144 | PHYSICS->mcConfig().AddHadronicId(-3101);
|
|---|
| 145 | PHYSICS->mcConfig().AddHadronicId(-2224);
|
|---|
| 146 | PHYSICS->mcConfig().AddHadronicId(-2214);
|
|---|
| 147 | PHYSICS->mcConfig().AddHadronicId(-2212);
|
|---|
| 148 | PHYSICS->mcConfig().AddHadronicId(-2203);
|
|---|
| 149 | PHYSICS->mcConfig().AddHadronicId(-2114);
|
|---|
| 150 | PHYSICS->mcConfig().AddHadronicId(-2112);
|
|---|
| 151 | PHYSICS->mcConfig().AddHadronicId(-2103);
|
|---|
| 152 | PHYSICS->mcConfig().AddHadronicId(-2101);
|
|---|
| 153 | PHYSICS->mcConfig().AddHadronicId(-1114);
|
|---|
| 154 | PHYSICS->mcConfig().AddHadronicId(-1103);
|
|---|
| 155 | PHYSICS->mcConfig().AddHadronicId(-545);
|
|---|
| 156 | PHYSICS->mcConfig().AddHadronicId(-543);
|
|---|
| 157 | PHYSICS->mcConfig().AddHadronicId(-541);
|
|---|
| 158 | PHYSICS->mcConfig().AddHadronicId(-535);
|
|---|
| 159 | PHYSICS->mcConfig().AddHadronicId(-533);
|
|---|
| 160 | PHYSICS->mcConfig().AddHadronicId(-531);
|
|---|
| 161 | PHYSICS->mcConfig().AddHadronicId(-525);
|
|---|
| 162 | PHYSICS->mcConfig().AddHadronicId(-523);
|
|---|
| 163 | PHYSICS->mcConfig().AddHadronicId(-521);
|
|---|
| 164 | PHYSICS->mcConfig().AddHadronicId(-515);
|
|---|
| 165 | PHYSICS->mcConfig().AddHadronicId(-513);
|
|---|
| 166 | PHYSICS->mcConfig().AddHadronicId(-511);
|
|---|
| 167 | PHYSICS->mcConfig().AddHadronicId(-435);
|
|---|
| 168 | PHYSICS->mcConfig().AddHadronicId(-433);
|
|---|
| 169 | PHYSICS->mcConfig().AddHadronicId(-431);
|
|---|
| 170 | PHYSICS->mcConfig().AddHadronicId(-425);
|
|---|
| 171 | PHYSICS->mcConfig().AddHadronicId(-423);
|
|---|
| 172 | PHYSICS->mcConfig().AddHadronicId(-421);
|
|---|
| 173 | PHYSICS->mcConfig().AddHadronicId(-415);
|
|---|
| 174 | PHYSICS->mcConfig().AddHadronicId(-413);
|
|---|
| 175 | PHYSICS->mcConfig().AddHadronicId(-411);
|
|---|
| 176 | PHYSICS->mcConfig().AddHadronicId(-325);
|
|---|
| 177 | PHYSICS->mcConfig().AddHadronicId(-323);
|
|---|
| 178 | PHYSICS->mcConfig().AddHadronicId(-321);
|
|---|
| 179 | PHYSICS->mcConfig().AddHadronicId(-315);
|
|---|
| 180 | PHYSICS->mcConfig().AddHadronicId(-313);
|
|---|
| 181 | PHYSICS->mcConfig().AddHadronicId(-311);
|
|---|
| 182 | PHYSICS->mcConfig().AddHadronicId(-215);
|
|---|
| 183 | PHYSICS->mcConfig().AddHadronicId(-213);
|
|---|
| 184 | PHYSICS->mcConfig().AddHadronicId(-211);
|
|---|
| 185 | PHYSICS->mcConfig().AddHadronicId(111);
|
|---|
| 186 | PHYSICS->mcConfig().AddHadronicId(113);
|
|---|
| 187 | PHYSICS->mcConfig().AddHadronicId(115);
|
|---|
| 188 | PHYSICS->mcConfig().AddHadronicId(130);
|
|---|
| 189 | PHYSICS->mcConfig().AddHadronicId(211);
|
|---|
| 190 | PHYSICS->mcConfig().AddHadronicId(213);
|
|---|
| 191 | PHYSICS->mcConfig().AddHadronicId(215);
|
|---|
| 192 | PHYSICS->mcConfig().AddHadronicId(221);
|
|---|
| 193 | PHYSICS->mcConfig().AddHadronicId(223);
|
|---|
| 194 | PHYSICS->mcConfig().AddHadronicId(225);
|
|---|
| 195 | PHYSICS->mcConfig().AddHadronicId(310);
|
|---|
| 196 | PHYSICS->mcConfig().AddHadronicId(311);
|
|---|
| 197 | PHYSICS->mcConfig().AddHadronicId(313);
|
|---|
| 198 | PHYSICS->mcConfig().AddHadronicId(315);
|
|---|
| 199 | PHYSICS->mcConfig().AddHadronicId(321);
|
|---|
| 200 | PHYSICS->mcConfig().AddHadronicId(323);
|
|---|
| 201 | PHYSICS->mcConfig().AddHadronicId(325);
|
|---|
| 202 | PHYSICS->mcConfig().AddHadronicId(331);
|
|---|
| 203 | PHYSICS->mcConfig().AddHadronicId(333);
|
|---|
| 204 | PHYSICS->mcConfig().AddHadronicId(335);
|
|---|
| 205 | PHYSICS->mcConfig().AddHadronicId(411);
|
|---|
| 206 | PHYSICS->mcConfig().AddHadronicId(413);
|
|---|
| 207 | PHYSICS->mcConfig().AddHadronicId(415);
|
|---|
| 208 | PHYSICS->mcConfig().AddHadronicId(421);
|
|---|
| 209 | PHYSICS->mcConfig().AddHadronicId(423);
|
|---|
| 210 | PHYSICS->mcConfig().AddHadronicId(425);
|
|---|
| 211 | PHYSICS->mcConfig().AddHadronicId(431);
|
|---|
| 212 | PHYSICS->mcConfig().AddHadronicId(433);
|
|---|
| 213 | PHYSICS->mcConfig().AddHadronicId(435);
|
|---|
| 214 | PHYSICS->mcConfig().AddHadronicId(441);
|
|---|
| 215 | PHYSICS->mcConfig().AddHadronicId(443);
|
|---|
| 216 | PHYSICS->mcConfig().AddHadronicId(445);
|
|---|
| 217 | PHYSICS->mcConfig().AddHadronicId(511);
|
|---|
| 218 | PHYSICS->mcConfig().AddHadronicId(513);
|
|---|
| 219 | PHYSICS->mcConfig().AddHadronicId(515);
|
|---|
| 220 | PHYSICS->mcConfig().AddHadronicId(521);
|
|---|
| 221 | PHYSICS->mcConfig().AddHadronicId(523);
|
|---|
| 222 | PHYSICS->mcConfig().AddHadronicId(525);
|
|---|
| 223 | PHYSICS->mcConfig().AddHadronicId(531);
|
|---|
| 224 | PHYSICS->mcConfig().AddHadronicId(533);
|
|---|
| 225 | PHYSICS->mcConfig().AddHadronicId(535);
|
|---|
| 226 | PHYSICS->mcConfig().AddHadronicId(541);
|
|---|
| 227 | PHYSICS->mcConfig().AddHadronicId(543);
|
|---|
| 228 | PHYSICS->mcConfig().AddHadronicId(545);
|
|---|
| 229 | PHYSICS->mcConfig().AddHadronicId(551);
|
|---|
| 230 | PHYSICS->mcConfig().AddHadronicId(553);
|
|---|
| 231 | PHYSICS->mcConfig().AddHadronicId(555);
|
|---|
| 232 | PHYSICS->mcConfig().AddHadronicId(1103);
|
|---|
| 233 | PHYSICS->mcConfig().AddHadronicId(1114);
|
|---|
| 234 | PHYSICS->mcConfig().AddHadronicId(2101);
|
|---|
| 235 | PHYSICS->mcConfig().AddHadronicId(2103);
|
|---|
| 236 | PHYSICS->mcConfig().AddHadronicId(2112);
|
|---|
| 237 | PHYSICS->mcConfig().AddHadronicId(2114);
|
|---|
| 238 | PHYSICS->mcConfig().AddHadronicId(2203);
|
|---|
| 239 | PHYSICS->mcConfig().AddHadronicId(2212);
|
|---|
| 240 | PHYSICS->mcConfig().AddHadronicId(2214);
|
|---|
| 241 | PHYSICS->mcConfig().AddHadronicId(2224);
|
|---|
| 242 | PHYSICS->mcConfig().AddHadronicId(3101);
|
|---|
| 243 | PHYSICS->mcConfig().AddHadronicId(3103);
|
|---|
| 244 | PHYSICS->mcConfig().AddHadronicId(3112);
|
|---|
| 245 | PHYSICS->mcConfig().AddHadronicId(3114);
|
|---|
| 246 | PHYSICS->mcConfig().AddHadronicId(3122);
|
|---|
| 247 | PHYSICS->mcConfig().AddHadronicId(3201);
|
|---|
| 248 | PHYSICS->mcConfig().AddHadronicId(3203);
|
|---|
| 249 | PHYSICS->mcConfig().AddHadronicId(3212);
|
|---|
| 250 | PHYSICS->mcConfig().AddHadronicId(3214);
|
|---|
| 251 | PHYSICS->mcConfig().AddHadronicId(3222);
|
|---|
| 252 | PHYSICS->mcConfig().AddHadronicId(3224);
|
|---|
| 253 | PHYSICS->mcConfig().AddHadronicId(3303);
|
|---|
| 254 | PHYSICS->mcConfig().AddHadronicId(3312);
|
|---|
| 255 | PHYSICS->mcConfig().AddHadronicId(3314);
|
|---|
| 256 | PHYSICS->mcConfig().AddHadronicId(3322);
|
|---|
| 257 | PHYSICS->mcConfig().AddHadronicId(3324);
|
|---|
| 258 | PHYSICS->mcConfig().AddHadronicId(3334);
|
|---|
| 259 | PHYSICS->mcConfig().AddHadronicId(4101);
|
|---|
| 260 | PHYSICS->mcConfig().AddHadronicId(4103);
|
|---|
| 261 | PHYSICS->mcConfig().AddHadronicId(4112);
|
|---|
| 262 | PHYSICS->mcConfig().AddHadronicId(4114);
|
|---|
| 263 | PHYSICS->mcConfig().AddHadronicId(4122);
|
|---|
| 264 | PHYSICS->mcConfig().AddHadronicId(4132);
|
|---|
| 265 | PHYSICS->mcConfig().AddHadronicId(4201);
|
|---|
| 266 | PHYSICS->mcConfig().AddHadronicId(4203);
|
|---|
| 267 | PHYSICS->mcConfig().AddHadronicId(4212);
|
|---|
| 268 | PHYSICS->mcConfig().AddHadronicId(4214);
|
|---|
| 269 | PHYSICS->mcConfig().AddHadronicId(4222);
|
|---|
| 270 | PHYSICS->mcConfig().AddHadronicId(4224);
|
|---|
| 271 | PHYSICS->mcConfig().AddHadronicId(4232);
|
|---|
| 272 | PHYSICS->mcConfig().AddHadronicId(4301);
|
|---|
| 273 | PHYSICS->mcConfig().AddHadronicId(4303);
|
|---|
| 274 | PHYSICS->mcConfig().AddHadronicId(4312);
|
|---|
| 275 | PHYSICS->mcConfig().AddHadronicId(4314);
|
|---|
| 276 | PHYSICS->mcConfig().AddHadronicId(4322);
|
|---|
| 277 | PHYSICS->mcConfig().AddHadronicId(4324);
|
|---|
| 278 | PHYSICS->mcConfig().AddHadronicId(4332);
|
|---|
| 279 | PHYSICS->mcConfig().AddHadronicId(4334);
|
|---|
| 280 | PHYSICS->mcConfig().AddHadronicId(4403);
|
|---|
| 281 | PHYSICS->mcConfig().AddHadronicId(4412);
|
|---|
| 282 | PHYSICS->mcConfig().AddHadronicId(4414);
|
|---|
| 283 | PHYSICS->mcConfig().AddHadronicId(4422);
|
|---|
| 284 | PHYSICS->mcConfig().AddHadronicId(4424);
|
|---|
| 285 | PHYSICS->mcConfig().AddHadronicId(4432);
|
|---|
| 286 | PHYSICS->mcConfig().AddHadronicId(4434);
|
|---|
| 287 | PHYSICS->mcConfig().AddHadronicId(4444);
|
|---|
| 288 | PHYSICS->mcConfig().AddHadronicId(5101);
|
|---|
| 289 | PHYSICS->mcConfig().AddHadronicId(5103);
|
|---|
| 290 | PHYSICS->mcConfig().AddHadronicId(5112);
|
|---|
| 291 | PHYSICS->mcConfig().AddHadronicId(5114);
|
|---|
| 292 | PHYSICS->mcConfig().AddHadronicId(5122);
|
|---|
| 293 | PHYSICS->mcConfig().AddHadronicId(5132);
|
|---|
| 294 | PHYSICS->mcConfig().AddHadronicId(5142);
|
|---|
| 295 | PHYSICS->mcConfig().AddHadronicId(5201);
|
|---|
| 296 | PHYSICS->mcConfig().AddHadronicId(5203);
|
|---|
| 297 | PHYSICS->mcConfig().AddHadronicId(5212);
|
|---|
| 298 | PHYSICS->mcConfig().AddHadronicId(5214);
|
|---|
| 299 | PHYSICS->mcConfig().AddHadronicId(5222);
|
|---|
| 300 | PHYSICS->mcConfig().AddHadronicId(5224);
|
|---|
| 301 | PHYSICS->mcConfig().AddHadronicId(5232);
|
|---|
| 302 | PHYSICS->mcConfig().AddHadronicId(5242);
|
|---|
| 303 | PHYSICS->mcConfig().AddHadronicId(5301);
|
|---|
| 304 | PHYSICS->mcConfig().AddHadronicId(5303);
|
|---|
| 305 | PHYSICS->mcConfig().AddHadronicId(5312);
|
|---|
| 306 | PHYSICS->mcConfig().AddHadronicId(5314);
|
|---|
| 307 | PHYSICS->mcConfig().AddHadronicId(5322);
|
|---|
| 308 | PHYSICS->mcConfig().AddHadronicId(5324);
|
|---|
| 309 | PHYSICS->mcConfig().AddHadronicId(5332);
|
|---|
| 310 | PHYSICS->mcConfig().AddHadronicId(5334);
|
|---|
| 311 | PHYSICS->mcConfig().AddHadronicId(5342);
|
|---|
| 312 | PHYSICS->mcConfig().AddHadronicId(5401);
|
|---|
| 313 | PHYSICS->mcConfig().AddHadronicId(5403);
|
|---|
| 314 | PHYSICS->mcConfig().AddHadronicId(5412);
|
|---|
| 315 | PHYSICS->mcConfig().AddHadronicId(5414);
|
|---|
| 316 | PHYSICS->mcConfig().AddHadronicId(5422);
|
|---|
| 317 | PHYSICS->mcConfig().AddHadronicId(5424);
|
|---|
| 318 | PHYSICS->mcConfig().AddHadronicId(5432);
|
|---|
| 319 | PHYSICS->mcConfig().AddHadronicId(5434);
|
|---|
| 320 | PHYSICS->mcConfig().AddHadronicId(5442);
|
|---|
| 321 | PHYSICS->mcConfig().AddHadronicId(5444);
|
|---|
| 322 | PHYSICS->mcConfig().AddHadronicId(5503);
|
|---|
| 323 | PHYSICS->mcConfig().AddHadronicId(5512);
|
|---|
| 324 | PHYSICS->mcConfig().AddHadronicId(5514);
|
|---|
| 325 | PHYSICS->mcConfig().AddHadronicId(5522);
|
|---|
| 326 | PHYSICS->mcConfig().AddHadronicId(5524);
|
|---|
| 327 | PHYSICS->mcConfig().AddHadronicId(5532);
|
|---|
| 328 | PHYSICS->mcConfig().AddHadronicId(5534);
|
|---|
| 329 | PHYSICS->mcConfig().AddHadronicId(5542);
|
|---|
| 330 | PHYSICS->mcConfig().AddHadronicId(5544);
|
|---|
| 331 | PHYSICS->mcConfig().AddHadronicId(5554);
|
|---|
| 332 | PHYSICS->mcConfig().AddHadronicId(10111);
|
|---|
| 333 | PHYSICS->mcConfig().AddHadronicId(10113);
|
|---|
| 334 | PHYSICS->mcConfig().AddHadronicId(10211);
|
|---|
| 335 | PHYSICS->mcConfig().AddHadronicId(10213);
|
|---|
| 336 | PHYSICS->mcConfig().AddHadronicId(10221);
|
|---|
| 337 | PHYSICS->mcConfig().AddHadronicId(10223);
|
|---|
| 338 | PHYSICS->mcConfig().AddHadronicId(10311);
|
|---|
| 339 | PHYSICS->mcConfig().AddHadronicId(10313);
|
|---|
| 340 | PHYSICS->mcConfig().AddHadronicId(10321);
|
|---|
| 341 | PHYSICS->mcConfig().AddHadronicId(10323);
|
|---|
| 342 | PHYSICS->mcConfig().AddHadronicId(10331);
|
|---|
| 343 | PHYSICS->mcConfig().AddHadronicId(10333);
|
|---|
| 344 | PHYSICS->mcConfig().AddHadronicId(10411);
|
|---|
| 345 | PHYSICS->mcConfig().AddHadronicId(10413);
|
|---|
| 346 | PHYSICS->mcConfig().AddHadronicId(10421);
|
|---|
| 347 | PHYSICS->mcConfig().AddHadronicId(10423);
|
|---|
| 348 | PHYSICS->mcConfig().AddHadronicId(10431);
|
|---|
| 349 | PHYSICS->mcConfig().AddHadronicId(10433);
|
|---|
| 350 | PHYSICS->mcConfig().AddHadronicId(10441);
|
|---|
| 351 | PHYSICS->mcConfig().AddHadronicId(10443);
|
|---|
| 352 | PHYSICS->mcConfig().AddHadronicId(10511);
|
|---|
| 353 | PHYSICS->mcConfig().AddHadronicId(10513);
|
|---|
| 354 | PHYSICS->mcConfig().AddHadronicId(10521);
|
|---|
| 355 | PHYSICS->mcConfig().AddHadronicId(10523);
|
|---|
| 356 | PHYSICS->mcConfig().AddHadronicId(10531);
|
|---|
| 357 | PHYSICS->mcConfig().AddHadronicId(10533);
|
|---|
| 358 | PHYSICS->mcConfig().AddHadronicId(10541);
|
|---|
| 359 | PHYSICS->mcConfig().AddHadronicId(10543);
|
|---|
| 360 | PHYSICS->mcConfig().AddHadronicId(10551);
|
|---|
| 361 | PHYSICS->mcConfig().AddHadronicId(10553);
|
|---|
| 362 | PHYSICS->mcConfig().AddHadronicId(20113);
|
|---|
| 363 | PHYSICS->mcConfig().AddHadronicId(20213);
|
|---|
| 364 | PHYSICS->mcConfig().AddHadronicId(20223);
|
|---|
| 365 | PHYSICS->mcConfig().AddHadronicId(20313);
|
|---|
| 366 | PHYSICS->mcConfig().AddHadronicId(20323);
|
|---|
| 367 | PHYSICS->mcConfig().AddHadronicId(20333);
|
|---|
| 368 | PHYSICS->mcConfig().AddHadronicId(20413);
|
|---|
| 369 | PHYSICS->mcConfig().AddHadronicId(20423);
|
|---|
| 370 | PHYSICS->mcConfig().AddHadronicId(20433);
|
|---|
| 371 | PHYSICS->mcConfig().AddHadronicId(20443);
|
|---|
| 372 | PHYSICS->mcConfig().AddHadronicId(20513);
|
|---|
| 373 | PHYSICS->mcConfig().AddHadronicId(20523);
|
|---|
| 374 | PHYSICS->mcConfig().AddHadronicId(20533);
|
|---|
| 375 | PHYSICS->mcConfig().AddHadronicId(20543);
|
|---|
| 376 | PHYSICS->mcConfig().AddHadronicId(20553);
|
|---|
| 377 | PHYSICS->mcConfig().AddHadronicId(100443);
|
|---|
| 378 | PHYSICS->mcConfig().AddHadronicId(100553);
|
|---|
| 379 | PHYSICS->mcConfig().AddHadronicId(9900440);
|
|---|
| 380 | PHYSICS->mcConfig().AddHadronicId(9900441);
|
|---|
| 381 | PHYSICS->mcConfig().AddHadronicId(9900443);
|
|---|
| 382 | PHYSICS->mcConfig().AddHadronicId(9900551);
|
|---|
| 383 | PHYSICS->mcConfig().AddHadronicId(9900553);
|
|---|
| 384 | PHYSICS->mcConfig().AddHadronicId(9910441);
|
|---|
| 385 | PHYSICS->mcConfig().AddHadronicId(9910551);
|
|---|
| 386 |
|
|---|
| 387 | // definition of the multiparticle "invisible"
|
|---|
| 388 | PHYSICS->mcConfig().AddInvisibleId(-16);
|
|---|
| 389 | PHYSICS->mcConfig().AddInvisibleId(-14);
|
|---|
| 390 | PHYSICS->mcConfig().AddInvisibleId(-12);
|
|---|
| 391 | PHYSICS->mcConfig().AddInvisibleId(12);
|
|---|
| 392 | PHYSICS->mcConfig().AddInvisibleId(14);
|
|---|
| 393 | PHYSICS->mcConfig().AddInvisibleId(16);
|
|---|
| 394 | PHYSICS->mcConfig().AddInvisibleId(1000022);
|
|---|
| 395 | PHYSICS->mcConfig().AddInvisibleId(1000039);
|
|---|
| 396 |
|
|---|
| 397 | // Initializing PhysicsService for RECO
|
|---|
| 398 | PHYSICS->recConfig().Reset();
|
|---|
| 399 |
|
|---|
| 400 | // Signal regions
|
|---|
| 401 | Manager()->AddRegionSelection("Ewkino_lowMET_M_4to9");
|
|---|
| 402 | Manager()->AddRegionSelection("Ewkino_lowMET_M_10to20");
|
|---|
| 403 | Manager()->AddRegionSelection("Ewkino_lowMET_M_20to30");
|
|---|
| 404 | Manager()->AddRegionSelection("Ewkino_lowMET_M_30to50");
|
|---|
| 405 | Manager()->AddRegionSelection("Ewkino_medMET_M_4to9");
|
|---|
| 406 | Manager()->AddRegionSelection("Ewkino_medMET_M_10to20");
|
|---|
| 407 | Manager()->AddRegionSelection("Ewkino_medMET_M_20to30");
|
|---|
| 408 | Manager()->AddRegionSelection("Ewkino_medMET_M_30to50");
|
|---|
| 409 | Manager()->AddRegionSelection("Ewkino_highMET_M_4to9");
|
|---|
| 410 | Manager()->AddRegionSelection("Ewkino_highMET_M_10to20");
|
|---|
| 411 | Manager()->AddRegionSelection("Ewkino_highMET_M_20to30");
|
|---|
| 412 | Manager()->AddRegionSelection("Ewkino_highMET_M_30to50");
|
|---|
| 413 |
|
|---|
| 414 | Manager()->AddRegionSelection("stop_lowMET_PT_5to12");
|
|---|
| 415 | Manager()->AddRegionSelection("stop_lowMET_PT_12to20");
|
|---|
| 416 | Manager()->AddRegionSelection("stop_lowMET_PT_20to30");
|
|---|
| 417 | Manager()->AddRegionSelection("stop_medMET_PT_5to12");
|
|---|
| 418 | Manager()->AddRegionSelection("stop_medMET_PT_12to20");
|
|---|
| 419 | Manager()->AddRegionSelection("stop_medMET_PT_20to30");
|
|---|
| 420 | Manager()->AddRegionSelection("stop_highMET_PT_5to12");
|
|---|
| 421 | Manager()->AddRegionSelection("stop_highMET_PT_12to20");
|
|---|
| 422 | Manager()->AddRegionSelection("stop_highMET_PT_20to30");
|
|---|
| 423 |
|
|---|
| 424 | // Preselection cuts
|
|---|
| 425 | std::string ewkinos[] = {"Ewkino_lowMET_M_4to9", "Ewkino_lowMET_M_10to20",
|
|---|
| 426 | "Ewkino_lowMET_M_20to30", "Ewkino_lowMET_M_30to50", "Ewkino_medMET_M_4to9",
|
|---|
| 427 | "Ewkino_medMET_M_10to20", "Ewkino_medMET_M_20to30", "Ewkino_medMET_M_30to50",
|
|---|
| 428 | "Ewkino_highMET_M_4to9", "Ewkino_highMET_M_10to20","Ewkino_highMET_M_20to30",
|
|---|
| 429 | "Ewkino_highMET_M_30to50"};
|
|---|
| 430 | std::string non_stpHM[] = {"Ewkino_lowMET_M_4to9", "Ewkino_lowMET_M_10to20",
|
|---|
| 431 | "Ewkino_lowMET_M_20to30", "Ewkino_lowMET_M_30to50", "Ewkino_medMET_M_4to9",
|
|---|
| 432 | "Ewkino_medMET_M_10to20", "Ewkino_medMET_M_20to30", "Ewkino_medMET_M_30to50",
|
|---|
| 433 | "Ewkino_highMET_M_4to9", "Ewkino_highMET_M_10to20","Ewkino_highMET_M_20to30",
|
|---|
| 434 | "Ewkino_highMET_M_30to50", "stop_lowMET_PT_5to12","stop_lowMET_PT_12to20",
|
|---|
| 435 | "stop_lowMET_PT_20to30", "stop_medMET_PT_5to12", "stop_medMET_PT_12to20",
|
|---|
| 436 | "stop_medMET_PT_20to30"};
|
|---|
| 437 | std::string lowMET[] = {"Ewkino_lowMET_M_4to9", "Ewkino_lowMET_M_10to20",
|
|---|
| 438 | "Ewkino_lowMET_M_20to30", "Ewkino_lowMET_M_30to50", "stop_lowMET_PT_5to12",
|
|---|
| 439 | "stop_lowMET_PT_12to20", "stop_lowMET_PT_20to30"};
|
|---|
| 440 | Manager()->AddCut("Lepton_1_init");
|
|---|
| 441 | Manager()->AddCut("Lepton_1");
|
|---|
| 442 | Manager()->AddCut("Lepton_2");
|
|---|
| 443 | Manager()->AddCut("Lepton_2_nonStopHM", non_stpHM);
|
|---|
| 444 | Manager()->AddCut("Dimuon",lowMET);
|
|---|
| 445 | Manager()->AddCut("SameFlavour",ewkinos);
|
|---|
| 446 | Manager()->AddCut("OSdilepton");
|
|---|
| 447 | Manager()->AddCut("pt(ll)");
|
|---|
| 448 | Manager()->AddCut("M(ll)_bounds", ewkinos);
|
|---|
| 449 | Manager()->AddCut("M(ll)_hadronic_veto", ewkinos);
|
|---|
| 450 |
|
|---|
| 451 | // MET cuts
|
|---|
| 452 | std::string ewk_medMET[] = {"Ewkino_medMET_M_4to9", "Ewkino_medMET_M_10to20",
|
|---|
| 453 | "Ewkino_medMET_M_20to30", "Ewkino_medMET_M_30to50"};
|
|---|
| 454 | std::string ewk_highMET[] = {"Ewkino_highMET_M_4to9",
|
|---|
| 455 | "Ewkino_highMET_M_10to20","Ewkino_highMET_M_20to30",
|
|---|
| 456 | "Ewkino_highMET_M_30to50"};
|
|---|
| 457 | std::string stp_medMET[] = {"stop_medMET_PT_5to12", "stop_medMET_PT_12to20",
|
|---|
| 458 | "stop_medMET_PT_20to30"};
|
|---|
| 459 | std::string stp_highMET[] = {"stop_highMET_PT_5to12", "stop_highMET_PT_12to20",
|
|---|
| 460 | "stop_highMET_PT_20to30"};
|
|---|
| 461 | Manager()->AddCut("MET");
|
|---|
| 462 | Manager()->AddCut("lowMET",lowMET);
|
|---|
| 463 | Manager()->AddCut("MET_corrected");
|
|---|
| 464 | Manager()->AddCut("lowMET_corrected",lowMET);
|
|---|
| 465 | Manager()->AddCut("ewk_medMET", ewk_medMET);
|
|---|
| 466 | Manager()->AddCut("ewk_highMET",ewk_highMET);
|
|---|
| 467 | Manager()->AddCut("stp_medMET", stp_medMET);
|
|---|
| 468 | Manager()->AddCut("stp_highMET",stp_highMET);
|
|---|
| 469 |
|
|---|
| 470 | // Follow up preselection
|
|---|
| 471 | Manager()->AddCut("ISR_jet");
|
|---|
| 472 | Manager()->AddCut("HT");
|
|---|
| 473 | Manager()->AddCut("MET/HT");
|
|---|
| 474 | Manager()->AddCut("b-veto");
|
|---|
| 475 | Manager()->AddCut("Mtautau");
|
|---|
| 476 | Manager()->AddCut("MT_MET",ewkinos);
|
|---|
| 477 |
|
|---|
| 478 | // Mll cuts
|
|---|
| 479 | std::string MLL4_9[] = {"Ewkino_lowMET_M_4to9", "Ewkino_medMET_M_4to9",
|
|---|
| 480 | "Ewkino_highMET_M_4to9"};
|
|---|
| 481 | std::string MLL10_20[] = {"Ewkino_lowMET_M_10to20", "Ewkino_medMET_M_10to20",
|
|---|
| 482 | "Ewkino_highMET_M_10to20"};
|
|---|
| 483 | std::string MLL20_30[] = {"Ewkino_lowMET_M_20to30", "Ewkino_medMET_M_20to30",
|
|---|
| 484 | "Ewkino_highMET_M_20to30"};
|
|---|
| 485 | std::string MLL30_50[] = {"Ewkino_lowMET_M_30to50", "Ewkino_medMET_M_30to50",
|
|---|
| 486 | "Ewkino_highMET_M_30to50"};
|
|---|
| 487 | Manager()->AddCut("MLL4_9",MLL4_9);
|
|---|
| 488 | Manager()->AddCut("MLL10_20",MLL10_20);
|
|---|
| 489 | Manager()->AddCut("MLL20_30",MLL20_30);
|
|---|
| 490 | Manager()->AddCut("MLL30_50",MLL30_50);
|
|---|
| 491 |
|
|---|
| 492 | // PT(lepton1)
|
|---|
| 493 | std::string PTL5_12[] = {"stop_lowMET_PT_5to12", "stop_medMET_PT_5to12",
|
|---|
| 494 | "stop_highMET_PT_5to12"};
|
|---|
| 495 | std::string PTL12_20[] = {"stop_lowMET_PT_12to20", "stop_medMET_PT_12to20",
|
|---|
| 496 | "stop_highMET_PT_12to20"};
|
|---|
| 497 | std::string PTL20_30[] = {"stop_lowMET_PT_20to30", "stop_medMET_PT_20to30",
|
|---|
| 498 | "stop_highMET_PT_20to30"};
|
|---|
| 499 | Manager()->AddCut("PTL5_12", PTL5_12);
|
|---|
| 500 | Manager()->AddCut("PTL12_20",PTL12_20);
|
|---|
| 501 | Manager()->AddCut("PTL20_30",PTL20_30);
|
|---|
| 502 |
|
|---|
| 503 |
|
|---|
| 504 | return true;
|
|---|
| 505 | }
|
|---|
| 506 |
|
|---|
| 507 | // -----------------------------------------------------------------------------
|
|---|
| 508 | // Finalize
|
|---|
| 509 | // function called one time at the end of the analysis
|
|---|
| 510 | // -----------------------------------------------------------------------------
|
|---|
| 511 | void sfs_cms_sus_16_048::Finalize(const SampleFormat& summary,
|
|---|
| 512 | const std::vector<SampleFormat>& files)
|
|---|
| 513 | {
|
|---|
| 514 | }
|
|---|
| 515 |
|
|---|
| 516 | // -----------------------------------------------------------------------------
|
|---|
| 517 | // Execute
|
|---|
| 518 | // function called each time one event is read
|
|---|
| 519 | // -----------------------------------------------------------------------------
|
|---|
| 520 | bool sfs_cms_sus_16_048::Execute(SampleFormat& sample, const EventFormat& event)
|
|---|
| 521 | {
|
|---|
| 522 | // Security for empty events
|
|---|
| 523 | if (event.rec()==0) return true;
|
|---|
| 524 |
|
|---|
| 525 | // Event weight
|
|---|
| 526 | double myEventWeight;
|
|---|
| 527 | if(Configuration().IsNoEventWeight()) myEventWeight=1.;
|
|---|
| 528 | else if(event.mc()->weight()!=0.) myEventWeight = event.mc()->weight();
|
|---|
| 529 | else return false;
|
|---|
| 530 | Manager()->InitializeForNewEvent(myEventWeight);
|
|---|
| 531 |
|
|---|
| 532 | // Leptons
|
|---|
| 533 | std::vector<const RecLeptonFormat*> Leptons;
|
|---|
| 534 | for(unsigned int ii=0; ii<event.rec()->electrons().size(); ii++)
|
|---|
| 535 | {
|
|---|
| 536 | const RecLeptonFormat *myElec = &(event.rec()->electrons()[ii]);
|
|---|
| 537 | if(myElec->pt()<5. || myElec->abseta()>2.5) continue;
|
|---|
| 538 | Leptons.push_back(myElec);
|
|---|
| 539 | }
|
|---|
| 540 | for(unsigned int ii=0; ii<event.rec()->muons().size(); ii++)
|
|---|
| 541 | {
|
|---|
| 542 | const RecLeptonFormat *myMuon = &(event.rec()->muons()[ii]);
|
|---|
| 543 | if(myMuon->pt()<3.5 || myMuon->abseta()>2.4) continue;
|
|---|
| 544 | Leptons.push_back(myMuon);
|
|---|
| 545 | }
|
|---|
| 546 | SORTER->sort(Leptons);
|
|---|
| 547 | unsigned int Nl = Leptons.size();
|
|---|
| 548 | if(!Manager()->ApplyCut((event.rec()->electrons().size()+event.rec()->muons().size())>0, "Lepton_1_init")) return true;
|
|---|
| 549 |
|
|---|
| 550 | // Jets
|
|---|
| 551 | std::vector<const RecJetFormat*> SignalJets;
|
|---|
| 552 | for(unsigned int ii=0; ii<event.rec()->jets().size(); ii++)
|
|---|
| 553 | {
|
|---|
| 554 | const RecJetFormat * myJet = &(event.rec()->jets()[ii]);
|
|---|
| 555 | if(myJet->pt()<25. || myJet->abseta()>2.4) continue;
|
|---|
| 556 | SignalJets.push_back(myJet);
|
|---|
| 557 | }
|
|---|
| 558 | SignalJets = PHYSICS->Isol->JetCleaning(SignalJets, Leptons, 0.2);
|
|---|
| 559 | double HT=0.;
|
|---|
| 560 | unsigned int Nb=0, Nj = SignalJets.size();
|
|---|
| 561 | for(unsigned int ii=0; ii<SignalJets.size(); ii++)
|
|---|
| 562 | {
|
|---|
| 563 | HT+=SignalJets[ii]->pt();
|
|---|
| 564 | if(SignalJets[ii]->btag()) Nb++;
|
|---|
| 565 | }
|
|---|
| 566 |
|
|---|
| 567 | // Requirement on the first lepton
|
|---|
| 568 | bool firstlpt = (Nl>=1) && (Leptons[0]->pt()>5.) && (Leptons[0]->pt()<30.);
|
|---|
| 569 | if(!Manager()->ApplyCut(firstlpt, "Lepton_1")) return true;
|
|---|
| 570 |
|
|---|
| 571 | // Requirement on the second lepton
|
|---|
| 572 | bool seclpt_general = (Nl>=2) && (Leptons[1]->pt()<30.);
|
|---|
| 573 | if(!Manager()->ApplyCut(seclpt_general, "Lepton_2")) return true;
|
|---|
| 574 | bool seclpt_non_stpHM = (Leptons[1]->pt()>5.);
|
|---|
| 575 | if(!Manager()->ApplyCut(seclpt_non_stpHM, "Lepton_2_nonStopHM")) return true;
|
|---|
| 576 |
|
|---|
| 577 | // 1 OS dilepton / flavour requirements
|
|---|
| 578 | bool osdilep = (Nl==2) && (Leptons[0]->charge()*Leptons[1]->charge()<0);
|
|---|
| 579 | bool is_SF = (Leptons[0]->isMuon()==Leptons[1]->isMuon());
|
|---|
| 580 | bool is_dimuon = is_SF && Leptons[0]->isMuon();
|
|---|
| 581 | if(!Manager()->ApplyCut(is_dimuon, "Dimuon")) return true;
|
|---|
| 582 | if(!Manager()->ApplyCut(is_SF, "SameFlavour")) return true;
|
|---|
| 583 | if(!Manager()->ApplyCut(osdilep, "OSdilepton")) return true;
|
|---|
| 584 |
|
|---|
| 585 | // Missing Transverse Energy
|
|---|
| 586 | MALorentzVector ptmiss = event.rec()->MET().momentum();
|
|---|
| 587 | double MET = ptmiss.Pt();
|
|---|
| 588 | double MET_cor = (ptmiss + Leptons[0]->momentum()+Leptons[1]->momentum()).Pt();
|
|---|
| 589 |
|
|---|
| 590 | // PT(ll) and Mll
|
|---|
| 591 | double ptll = (Leptons[0]->momentum()+Leptons[1]->momentum()).Pt();
|
|---|
| 592 | double mll = (Leptons[0]->momentum()+Leptons[1]->momentum()).M();
|
|---|
| 593 | bool mll_lims = is_SF && (mll>4.) && (mll<50.);
|
|---|
| 594 | bool mll_veto = is_SF && ((mll<9.) || (mll>10.5));
|
|---|
| 595 | if(!Manager()->ApplyCut(ptll>3., "pt(ll)")) return true;
|
|---|
| 596 | if(!Manager()->ApplyCut(mll_lims, "M(ll)_bounds")) return true;
|
|---|
| 597 | if(!Manager()->ApplyCut(mll_veto, "M(ll)_hadronic_veto")) return true;
|
|---|
| 598 |
|
|---|
| 599 | // MET cut
|
|---|
| 600 | if(!Manager()->ApplyCut(MET>125., "MET")) return true;
|
|---|
| 601 | if(!Manager()->ApplyCut(is_dimuon && MET<200., "lowMET")) return true;
|
|---|
| 602 | // Trigger efficiency
|
|---|
| 603 | if(is_dimuon && MET < 200.)
|
|---|
| 604 | { Manager()->SetCurrentEventWeight(0.65*myEventWeight); }
|
|---|
| 605 | // Corrected MET cuts
|
|---|
| 606 | if(!Manager()->ApplyCut(is_dimuon && MET_cor>125., "MET_corrected")) return true;
|
|---|
| 607 | if(!Manager()->ApplyCut(is_dimuon && MET_cor<200., "lowMET_corrected"))
|
|---|
| 608 | return true;
|
|---|
| 609 | if(!Manager()->ApplyCut(MET>200. && MET<250., "ewk_medMET")) return true;
|
|---|
| 610 | if(!Manager()->ApplyCut(MET>200. && MET<300., "stp_medMET")) return true;
|
|---|
| 611 | if(!Manager()->ApplyCut(MET>250., "ewk_highMET")) return true;
|
|---|
| 612 | if(!Manager()->ApplyCut(MET>300., "stp_highMET")) return true;
|
|---|
| 613 |
|
|---|
| 614 | // ISR jet
|
|---|
| 615 | if(!Manager()->ApplyCut(Nj>0, "ISR_jet")) return true;
|
|---|
| 616 | if(!Manager()->ApplyCut(HT>100., "HT")) return true;
|
|---|
| 617 |
|
|---|
| 618 | // MET to HT
|
|---|
| 619 | if(!Manager()->ApplyCut(MET/HT<1.4 && MET/HT>0.6, "MET/HT")) return true;
|
|---|
| 620 |
|
|---|
| 621 | // b-veto
|
|---|
| 622 | if(!Manager()->ApplyCut(Nb==0, "b-veto")) return true;
|
|---|
| 623 |
|
|---|
| 624 | // M(tau tau)
|
|---|
| 625 | double AA = -(ptmiss.Px()*Leptons[1]->py() - ptmiss.Py()*Leptons[1]->px())/
|
|---|
| 626 | (Leptons[0]->py()*Leptons[1]->px() - Leptons[0]->px()*Leptons[1]->py());
|
|---|
| 627 | double BB = (ptmiss.Px()*Leptons[0]->py() - ptmiss.Py()*Leptons[0]->px())/
|
|---|
| 628 | (Leptons[0]->py()*Leptons[1]->px() - Leptons[0]->px()*Leptons[1]->py());
|
|---|
| 629 | double mtautau = ((AA+1.)*Leptons[0]->momentum() + (BB+1.)*Leptons[1]->momentum()).M();
|
|---|
| 630 | if(!Manager()->ApplyCut(mtautau<0. || mtautau>160., "Mtautau")) return true;
|
|---|
| 631 |
|
|---|
| 632 | // MT(li, met)
|
|---|
| 633 | double mt1 = sqrt(2.*Leptons[0]->pt()*MET*(1.-cos(Leptons[0]->dphi_0_pi(ptmiss))));
|
|---|
| 634 | double mt2 = sqrt(2.*Leptons[1]->pt()*MET*(1.-cos(Leptons[1]->dphi_0_pi(ptmiss))));
|
|---|
| 635 | if(!Manager()->ApplyCut(std::max(mt1,mt2)<70., "MT_MET")) return true;
|
|---|
| 636 |
|
|---|
| 637 | // Mll
|
|---|
| 638 | if(!Manager()->ApplyCut(mll<9., "MLL4_9")) return true;
|
|---|
| 639 | if(!Manager()->ApplyCut(mll>10.5 && mll<20., "MLL10_20")) return true;
|
|---|
| 640 | if(!Manager()->ApplyCut(mll>20. && mll<30., "MLL20_30")) return true;
|
|---|
| 641 | if(!Manager()->ApplyCut(mll>30., "MLL30_50")) return true;
|
|---|
| 642 |
|
|---|
| 643 | // pt(l1)
|
|---|
| 644 | double pt1=Leptons[0]->pt();
|
|---|
| 645 | if(!Manager()->ApplyCut(pt1<12., "PTL5_12")) return true;
|
|---|
| 646 | if(!Manager()->ApplyCut(pt1>12. && pt1<20., "PTL12_20")) return true;
|
|---|
| 647 | if(!Manager()->ApplyCut(pt1>20. && pt1<30., "PTL20_30")) return true;
|
|---|
| 648 |
|
|---|
| 649 | // exit
|
|---|
| 650 | return true;
|
|---|
| 651 | }
|
|---|
| 652 |
|
|---|
| 653 | // Overlap Removal
|
|---|
| 654 | template<typename T1, typename T2> std::vector<const T1*> sfs_cms_sus_16_048::Removal(std::vector<const T1*> &v1,
|
|---|
| 655 | std::vector<const T2*> &v2, const double &drmin)
|
|---|
| 656 | {
|
|---|
| 657 | // Determining with objects should be removed
|
|---|
| 658 | std::vector<bool> mask(v1.size(),false);
|
|---|
| 659 | for (unsigned int j=0;j<v1.size();j++)
|
|---|
| 660 | for (unsigned int i=0;i<v2.size();i++)
|
|---|
| 661 | if (v2[i]->dr(v1[j]) < drmin)
|
|---|
| 662 | {
|
|---|
| 663 | mask[j]=true;
|
|---|
| 664 | break;
|
|---|
| 665 | }
|
|---|
| 666 |
|
|---|
| 667 | // Building the cleaned container
|
|---|
| 668 | std::vector<const T1*> cleaned_v1;
|
|---|
| 669 | for (unsigned int i=0;i<v1.size();i++)
|
|---|
| 670 | if (!mask[i]) cleaned_v1.push_back(v1[i]);
|
|---|
| 671 |
|
|---|
| 672 | return cleaned_v1;
|
|---|
| 673 | }
|
|---|
| 674 |
|
|---|
