CommonTools/AnalysisMuon.cc at master · pandolf/CommonTools · GitHub

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
#include "AnalysisMuon.h"


bool AnalysisMuon::passedMuonID() {

  bool passed = ( isGlobalMuonPromptTight && isAllTrackerMuon
                && pixelHits>0 && trackerHits>10
                && nMatchedStations>=2
                && (fabs(dxy)<0.02) && (fabs(dz)<1.) );

  return passed;

}


bool AnalysisMuon::passedMuonID_tight2012() {

  bool passed = ( isGlobalMuon && isPFMuon
                && normChiSquare < 10.
                && pixelHits>0 && trackerHits>5
                && nValidMuonHits > 0
                && nMatchedStations>=2
                && (fabs(dxy)<0.02) && (fabs(dz)<0.5) );

  return passed;

}


bool AnalysisMuon::isIsolated() {

  bool isIsolated = this->combinedIsoRel() < 0.15;

  return isIsolated;

}


bool AnalysisMuon::isIsolated2012_CutsDeltaBeta() {

  bool isIsolated = (this->pfIsolationSumPUPtR04/this->Pt()) < 0.15;

  return isIsolated;

}


bool AnalysisMuon::isIsolated2012_MVA() {

  float abseta = fabs(this->Eta());
  float pt = this->Pt();
  float bdtiso = mvaisoMuon;

  bool isIsolated = (
                             ( pt <= 20 && abseta >= 0.000 && abseta < 0.479 && bdtiso > 0.86 ) ||
                             ( pt <= 20 && abseta >= 0.479 && abseta < 2.400 && bdtiso > 0.82 ) ||
                             ( pt >  20 && abseta >= 0.000 && abseta < 0.479 && bdtiso > 0.82 ) ||
                             ( pt >  20 && abseta >= 0.479 && abseta < 2.400 && bdtiso > 0.86 )
                             );


  return isIsolated;

}


float AnalysisMuon::combinedIsoRel() {

 return (sumPt03 + emEt03 + hadEt03)/this->Pt();

}


bool AnalysisMuon::passedVBTF() {

  bool isIsolated = this->isIsolated();
  bool passedMuonID = this->passedMuonID();

  return ( isIsolated && passedMuonID );

}


bool AnalysisMuon::isTightMuon2012() {

  bool isIsolated = this->isIsolated2012_CutsDeltaBeta();
  bool passedMuonID = this->passedMuonID_tight2012();

  return ( isIsolated && passedMuonID );

}