forked from piedraj/AnalysisCMS
-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathAnalysisCMS.h
More file actions
601 lines (511 loc) · 22.2 KB
/
Copy pathAnalysisCMS.h
File metadata and controls
601 lines (511 loc) · 22.2 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
#ifndef AnalysisCMS_h
#define AnalysisCMS_h
#include "AnalysisBase.h"
#include "Constants.h"
#include <fstream>
#include <iostream>
#include <TH1.h>
#include <TH2.h>
#include <TLorentzVector.h>
#include <TMatrixDSym.h>
#include <TMatrixDSymEigen.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TTree.h>
#include <TSpline.h>
#include "TRandom3.h"
enum {Loose, Tight};
struct Lepton
{
int index;
int type;
int flavour;
float idisoW;
float iso;
float motherPID;
TLorentzVector v;
};
struct Jet
{
int index;
float cmvav2;
float csvv2ivf;
float mass;
TLorentzVector v;
};
class AnalysisCMS : public AnalysisBase
{
public :
AnalysisCMS(TTree* tree = 0, TString systematic = "nominal");
void ApplyWeights ();
void DefineHistograms (int ichannel,
int icut,
int ijet,
TString suffix);
float ElectronIsolation (int k);
void EndJob ();
void EventDump (Bool_t leptonInfo = false);
void EventSetup (float jet_eta_max = -1,
float jet_pt_min = 30);
void FillHistograms (int ichannel,
int icut,
int ijet);
void GetFakeWeights ();
void GetGenPtllWeight ();
void GetGenWeightsLHE ();
void GetJets (float jet_eta_max,
float jet_pt_min);
void GetLeptons ();
void GetTops ();
void GetGenLeptonsAndNeutrinos();
void GetDark ();
void GetMlb ();
void GetTopReco ();
void GetMET (float module,
float phi);
void GetTrkMET (float module,
float phi);
float MuonIsolation (int k);
void OpenMinitree ();
bool PassTrigger ();
bool ApplyMETFilters (bool ApplyGiovanniFilters = false,
bool ApplyICHEPAdditionalFilters = true);
void PrintProgress (Long64_t counter,
Long64_t total);
void SetSaveMinitree (Bool_t saveminitree) {_saveminitree = saveminitree;}
void SetMinitreePath (TString minitreepath) {_minitreepath = minitreepath;}
void Setup (TString analysis,
TString filename,
float luminosity,
TString suffix = "");
//void ShapeWZ ();
//void ShapeZZ ();
//void ShapeWZforZZ ();
void Summary (TString analysis,
TString precision,
TString title);
// New analysis variables
//----------------------------------------------------------------------------
void GetDeltaPhi ();
void GetDeltaPhiVeto ();
void GetDeltaR ();
void GetHt ();
void GetJetPtSum ();
void GetMc ();
void GetMetVar ();
void GetMpMet ();
void GetPtWW ();
void GetSoftMuon ();
void GetStarVar ();
void GetZHCRVar ();
void GetZZGenVar ();
void GetHZZTo2L2Nu ();
void GetStopVar ();
void GetRazor ();
void GetMt (Lepton lep,
float& transverse_mass);
double ComputeMT2 (TLorentzVector VisibleA,
TLorentzVector VisibleB,
TLorentzVector Invisible,
int MT2Type = 0,
double MT2Precision = 0.);
int GetMotherPID (int index);
void GetScaleAndResolution();
TRandom3* _rand;
TRandom3* _rand2;
// Data members
//----------------------------------------------------------------------------
std::vector<Jet> AnalysisJets;
std::vector<Lepton> AnalysisLeptons;
Lepton Lepton1;
Lepton Lepton2;
TLorentzVector MET;
TLorentzVector trkMET;
std::vector<float> _jet_eta;
std::vector<float> _jet_phi;
std::vector<float> _jet_pt;
std::vector<float> _bjet30csvv2m_eta;
std::vector<float> _bjet30csvv2m_phi;
std::vector<float> _bjet30csvv2m_pt;
// to testing for the fakes
//std::vector<float> _lep1GenDeltaR;
//std::vector<float> _lep2GenDeltaR;
//std::vector<float> _testMotherID;
TH1F* _shapemlb;
bool _saveminitree;
bool _eventdump;
bool _foundsoftmuon;
bool _ismc;
bool _isminitree;
bool _isfastsim;
bool _passdphiveto;
bool _systematic_btag_do;
bool _systematic_btag_up;
bool _systematic_idiso_do;
bool _systematic_idiso_up;
bool _systematic_trigger_do;
bool _systematic_trigger_up;
bool _systematic_reco_do;
bool _systematic_reco_up;
bool _systematic_fastsim_do;
bool _systematic_fastsim_up;
bool _systematic_toppt;
bool _applytopptreweighting;
bool _applytriggeremulation;
TString _analysis;
TString _dataperiod;
TString _isdatadriven;
TString _filename;
TString _suffix;
TString _sample;
TString _systematic;
TString _longname;
TString _minitreepath;
float _channel;
float _darkpt;
float _darkpt_gen;
float _deltarjet1met;
float _deltarjet2met;
float _deltarjj;
float _deltarjjmet;
float _deltarlep1jet1;
float _deltarlep1jet2;
float _deltarlep2jet1;
float _deltarlep2jet2;
float _deltarllmet;
float _deltarl1met;
float _deltarl2met;
float _dphijet1met;
float _dphijet2met;
float _dphijj;
float _dphijjmet;
float _dphilep1jet1;
float _dphilep1jet2;
float _dphilep2jet1;
float _dphilep2jet2;
float _dphillmet;
float _dphillstar;
float _dphitt_gen;
float _detatt_gen;
float _event_weight;
float _event_weight_Btagup;
float _event_weight_Btagdo;
float _event_weight_Btaglightup;
float _event_weight_Btaglightdo;
float _event_weight_BtagFSup;
float _event_weight_BtagFSdo;
float _event_weight_Idisoup;
float _event_weight_Idisodo;
float _event_weight_Idisoeleup;
float _event_weight_Idisoeledo;
float _event_weight_Idisomuup;
float _event_weight_Idisomudo;
float _event_weight_Triggerup;
float _event_weight_Triggerdo;
float _event_weight_Recoup;
float _event_weight_Recodo;
float _event_weight_Fastsimup;
float _event_weight_Fastsimdo;
float _event_weight_Toppt;
float _event_weight_genmatched;
float _fake_weight;
float _fake_weight_elUp;
float _fake_weight_elDown;
float _fake_weight_elStatUp;
float _fake_weight_elStatDown;
float _fake_weight_muUp;
float _fake_weight_muDown;
float _fake_weight_muStatUp;
float _fake_weight_muStatDown;
float _gen_ptll_weight;
float _gen_ptll;
float _fullpmet;
float _ht;
float _htvisible;
float _htjets;
float _htnojets;
float _htgen;
float _isrpt;
float _detall;
float _lep1id;
float _lep1eta;
float _lep1phi;
float _lep1pt;
float _lep1mass;
float _lep1isfake;
float _lep1eta_gen;
float _lep1phi_gen;
float _lep1pt_gen;
float _lep2id;
float _lep2eta;
float _lep2phi;
float _lep2pt;
float _lep2mass;
float _lep2isfake;
float _lep2eta_gen;
float _lep2phi_gen;
float _lep2pt_gen;
float _lep1tau_gen;
float _lep1id_gen;
float _lep1motherid_gen;
float _lep2tau_gen;
float _lep2id_gen;
float _lep2motherid_gen;
float _lep1mid;
float _lep2mid;
float _lep3mid;
float _lep4mid;
float _lep3id;
float _lep3eta;
float _lep3phi;
float _lep3pt;
float _lep3mass;
float _lep3type;
float _lep3idisoW;
float _lep4id;
float _lep4eta;
float _lep4phi;
float _lep4pt;
float _lep4mass;
float _lep4type;
float _lep4idisoW;
float _nu1pt_gen;
float _nu1tau_gen;
float _nu2pt_gen;
float _nu2tau_gen;
float _luminosity;
float _mc;
float _metvar;
float _mllstar;
float _mpmet;
float _mtw;
float _m2l;
float _m3l;
float _m2t_gen;
float _leadingPtCSVv2L;
float _leadingPtCSVv2M;
float _leadingPtCSVv2T;
float _trailingPtCSVv2L;
float _trailingPtCSVv2M;
float _trailingPtCSVv2T;
float _nbjet30csvv2l;
float _nbjet30csvv2m;
float _nbjet30csvv2t;
float _nbjet20cmvav2l;
float _nbjet20cmvav2m;
float _nbjet20cmvav2t;
float _nbjet30cmvav2l;
float _nbjet30cmvav2m;
float _nbjet30cmvav2t;
float _njet;
float _nisrjet;
float _ptww;
float _pt2l;
float _sumjpt12;
float _trkpmet;
float _top1eta_gen;
float _top1phi_gen;
float _top1pt_gen;
float _top2eta_gen;
float _top2phi_gen;
float _top2pt_gen;
float _topReco;
float _topRecoW;
float _mll13;
float _mll23;
float _mll14;
float _mll24;
float _mll34;
float _dyll;
float _ptbll;
float _mt2ll;
float _mt2llgen;
float _dphimetptbll;
float _dphimetjet;
float _dphimetbbll;
float _mllbb;
float _meff;
float _mt2bb;
float _mt2lblb;
float _mlb1;
float _mlb2;
float _mt2lblbcomb;
float _mt2bbtrue;
float _mt2lblbtrue;
float _mt2lblbmatch;
float _mlb1comb;
float _mlb2comb;
float _mlb1true;
float _mlb2true;
float _bjet1pt;
float _bjet1eta;
float _bjet1phi;
float _bjet1mass;
float _bjet1csvv2ivf;
float _bjet2pt;
float _bjet2eta;
float _bjet2phi;
float _bjet2mass;
float _bjet2csvv2ivf;
float _tjet1pt;
float _tjet1eta;
float _tjet1phi;
float _tjet1mass;
float _tjet1csvv2ivf;
float _tjet1assignment;
float _tjet2pt;
float _tjet2eta;
float _tjet2phi;
float _tjet2mass;
float _tjet2csvv2ivf;
float _tjet2assignment;
float _ZZdphi;
float _ZZpt;
float _ZZmass;
float _MR;
float _R2;
float _Rpt;
float _invGamma;
float _Mdr;
float _DeltaPhiRll;
Long64_t _nentries;
unsigned int _jetbin;
unsigned int _nelectron;
//unsigned int _nlepton;
//unsigned int _ntightlepton;
int _nlepton;
int _ntightlepton;
int _verbosity;
ofstream txt_eventdump;
ofstream txt_summary;
TFile* root_output;
TFile* root_minitree;
TTree* minitree;
float _scale;
float _uPara;
float _uPerp;
// ZZ k-factors
// https://twiki.cern.ch/twiki/bin/viewauth/CMS/HiggsZZ4l2015#K_factorshttps://twiki.cern.ch/twiki/bin/viewauth/CMS/HiggsZZ4l2015#K_factors
TSpline3* ggZZ_kf[9];
void ReadGluGluHZZkfactors();
float kfactor_ggHZZ_qcd (float GENmassZZ, TString SystTitle);
float kfactor_qqZZ_qcd_dPhi(float GENdPhiZZ, int finalState);
float kfactor_qqZZ_qcd_M (float GENmassZZ, int finalState);
float kfactor_qqZZ_qcd_Pt (float GENpTZZ, int finalState);
// TH1 histograms
//----------------------------------------------------------------------------
TH1D* h_counterRaw [nchannel][ncut][njetbin+1];
TH1D* h_counterLum [nchannel][ncut][njetbin+1];
TH1D* h_fakes [nchannel][ncut][njetbin+1];
TH1D* h_nvtx [nchannel][ncut][njetbin+1];
TH1D* h_ptww [nchannel][ncut][njetbin+1];
TH1D* h_pt2l [nchannel][ncut][njetbin+1];
TH1D* h_sumjpt12 [nchannel][ncut][njetbin+1];
TH1D* h_sumpt12 [nchannel][ncut][njetbin+1];
TH1D* h_dyll [nchannel][ncut][njetbin+1];
TH1D* h_mllbb [nchannel][ncut][njetbin+1];
TH1D* h_dphimetjet [nchannel][ncut][njetbin+1];
TH1D* h_meff [nchannel][ncut][njetbin+1];
TH1D* h_ptbll [nchannel][ncut][njetbin+1];
TH1D* h_dphimetptbll [nchannel][ncut][njetbin+1];
TH1D* h_dphimetbbll [nchannel][ncut][njetbin+1];
TH1D* h_mt2ll [nchannel][ncut][njetbin+1];
TH1D* h_mt2bb [nchannel][ncut][njetbin+1];
TH1D* h_mt2lblb [nchannel][ncut][njetbin+1];
TH1D* h_mlb1 [nchannel][ncut][njetbin+1];
TH1D* h_mlb2 [nchannel][ncut][njetbin+1];
// TH1 histograms with minitree variables
//----------------------------------------------------------------------------
TH1D* h_channel [nchannel][ncut][njetbin+1];
TH1D* h_dphijet1met [nchannel][ncut][njetbin+1];
TH1D* h_dphijet2met [nchannel][ncut][njetbin+1];
TH1D* h_dphijj [nchannel][ncut][njetbin+1];
TH1D* h_dphijjmet [nchannel][ncut][njetbin+1];
TH1D* h_dphilep1jet1 [nchannel][ncut][njetbin+1];
TH1D* h_dphilep1jet2 [nchannel][ncut][njetbin+1];
TH1D* h_dphilep2jet1 [nchannel][ncut][njetbin+1];
TH1D* h_dphilep2jet2 [nchannel][ncut][njetbin+1];
TH1D* h_dphill [nchannel][ncut][njetbin+1];
TH1D* h_dphillmet [nchannel][ncut][njetbin+1];
TH1D* h_dphillstar [nchannel][ncut][njetbin+1];
TH1D* h_dphilmet1 [nchannel][ncut][njetbin+1];
TH1D* h_dphilmet2 [nchannel][ncut][njetbin+1];
TH1D* h_dphitt_gen [nchannel][ncut][njetbin+1];
TH1D* h_detatt_gen [nchannel][ncut][njetbin+1];
TH1D* h_drll [nchannel][ncut][njetbin+1];
TH1D* h_ht [nchannel][ncut][njetbin+1];
TH1D* h_htvisible [nchannel][ncut][njetbin+1];
TH1D* h_htgen [nchannel][ncut][njetbin+1];
TH1D* h_htjets [nchannel][ncut][njetbin+1];
TH1D* h_htnojets [nchannel][ncut][njetbin+1];
TH1D* h_jet1eta [nchannel][ncut][njetbin+1];
TH1D* h_jet1mass [nchannel][ncut][njetbin+1];
TH1D* h_jet1phi [nchannel][ncut][njetbin+1];
TH1D* h_jet1pt [nchannel][ncut][njetbin+1];
TH1D* h_jet2eta [nchannel][ncut][njetbin+1];
TH1D* h_jet2mass [nchannel][ncut][njetbin+1];
TH1D* h_jet2phi [nchannel][ncut][njetbin+1];
TH1D* h_jet2pt [nchannel][ncut][njetbin+1];
TH1D* h_detall [nchannel][ncut][njetbin+1];
TH1D* h_lep1eta [nchannel][ncut][njetbin+1];
TH1D* h_lep1phi [nchannel][ncut][njetbin+1];
TH1D* h_lep1pt [nchannel][ncut][njetbin+1];
TH1D* h_lep1eta_gen [nchannel][ncut][njetbin+1];
TH1D* h_lep1phi_gen [nchannel][ncut][njetbin+1];
TH1D* h_lep1pt_gen [nchannel][ncut][njetbin+1];
TH1D* h_lep2eta [nchannel][ncut][njetbin+1];
TH1D* h_lep2phi [nchannel][ncut][njetbin+1];
TH1D* h_lep2pt [nchannel][ncut][njetbin+1];
TH1D* h_lep2eta_gen [nchannel][ncut][njetbin+1];
TH1D* h_lep2phi_gen [nchannel][ncut][njetbin+1];
TH1D* h_lep2pt_gen [nchannel][ncut][njetbin+1];
TH1D* h_mc [nchannel][ncut][njetbin+1];
TH1D* h_metPfType1 [nchannel][ncut][njetbin+1];
TH1D* h_metPfType1Phi [nchannel][ncut][njetbin+1];
TH1D* h_metTtrk [nchannel][ncut][njetbin+1];
TH1D* h_metTtrkPhi [nchannel][ncut][njetbin+1];
TH1D* h_mpmet [nchannel][ncut][njetbin+1];
TH1D* h_metPuppi [nchannel][ncut][njetbin+1];
TH1D* h_mth [nchannel][ncut][njetbin+1];
TH1D* h_mtw1 [nchannel][ncut][njetbin+1];
TH1D* h_mtw2 [nchannel][ncut][njetbin+1];
TH1D* h_m2l [nchannel][ncut][njetbin+1];
TH1D* h_m2t_gen [nchannel][ncut][njetbin+1];
TH1D* h_nbjet15csvv2l [nchannel][ncut][njetbin+1];
TH1D* h_nbjet15csvv2m [nchannel][ncut][njetbin+1];
TH1D* h_nbjet15csvv2t [nchannel][ncut][njetbin+1];
TH1D* h_nbjet30csvv2l [nchannel][ncut][njetbin+1];
TH1D* h_nbjet30csvv2m [nchannel][ncut][njetbin+1];
TH1D* h_nbjet30csvv2t [nchannel][ncut][njetbin+1];
TH1D* h_nbjet20cmvav2l[nchannel][ncut][njetbin+1];
TH1D* h_nbjet20cmvav2m[nchannel][ncut][njetbin+1];
TH1D* h_nbjet20cmvav2t[nchannel][ncut][njetbin+1];
TH1D* h_nbjet30cmvav2l[nchannel][ncut][njetbin+1];
TH1D* h_nbjet30cmvav2m[nchannel][ncut][njetbin+1];
TH1D* h_nbjet30cmvav2t[nchannel][ncut][njetbin+1];
TH1D* h_njet [nchannel][ncut][njetbin+1];
TH1D* h_top1eta_gen [nchannel][ncut][njetbin+1];
TH1D* h_top1phi_gen [nchannel][ncut][njetbin+1];
TH1D* h_top1pt_gen [nchannel][ncut][njetbin+1];
TH1D* h_top2eta_gen [nchannel][ncut][njetbin+1];
TH1D* h_top2phi_gen [nchannel][ncut][njetbin+1];
TH1D* h_top2pt_gen [nchannel][ncut][njetbin+1];
TH1D* h_topReco [nchannel][ncut][njetbin+1];
TH1D* h_met_over_pt2l [nchannel][ncut][njetbin+1];
// Razor
//----------------------------------------------------------------------------
TH1D* h_MR [nchannel][ncut][njetbin+1];
TH1D* h_R2 [nchannel][ncut][njetbin+1];
TH1D* h_Rpt [nchannel][ncut][njetbin+1];
TH1D* h_invGamma [nchannel][ncut][njetbin+1];
TH1D* h_Mdr [nchannel][ncut][njetbin+1];
TH1D* h_DeltaPhiRll [nchannel][ncut][njetbin+1];
// TH2 histograms
//----------------------------------------------------------------------------
TH2D* h_metPfType1_m2l[nchannel][ncut][njetbin+1];
TH2D* h_mpmet_m2l [nchannel][ncut][njetbin+1];
TH2D* h_mt2ll_m2l [nchannel][ncut][njetbin+1];
TH2D* h_2ht [nchannel][ncut][njetbin+1];
TH2D* h_dym [nchannel][ncut][njetbin+1];
};
#endif