27 |
|
#include "SelectionFuncs.h" |
28 |
|
#include "ElectronEnergyCorrection.h" |
29 |
|
#include "ElectronEnergyRegression.h" |
30 |
+ |
#include "MuonMomentumCorrection.h" |
31 |
|
|
32 |
|
|
33 |
|
extern vector<SimpleLepton> failingLeptons; |
110 |
|
<< endl; |
111 |
|
} |
112 |
|
|
113 |
< |
|
113 |
> |
//correct muon momentum |
114 |
> |
if(ctrl.correct_muon_momentum){ |
115 |
> |
for(int i=0; i<muonArr->GetEntries(); i++) { |
116 |
> |
const mithep::Muon *const_mu = (mithep::Muon*)((*muonArr)[i]); |
117 |
> |
mithep::Muon *mu = const_cast<mithep::Muon*>(const_mu); |
118 |
> |
correct_muon_momentum(mu,info->EvtNum()); |
119 |
> |
} |
120 |
> |
} |
121 |
|
|
122 |
|
|
123 |
|
//******************************************************** |
127 |
|
int nlep_above_20=0; |
128 |
|
for(int i=0; i<muonArr->GetEntries(); i++) |
129 |
|
{ |
130 |
< |
const mithep::Muon *mu = (mithep::Muon*)((*muonArr)[i]); |
130 |
> |
const mithep::Muon *mu = (mithep::Muon*)((*muonArr)[i]); |
131 |
> |
|
132 |
> |
|
133 |
|
if( !(mu->IsTrackerMuon() || mu->IsGlobalMuon()) ) continue; |
134 |
|
if( fabs(mu->Eta()) > 2.4 ) continue; |
135 |
|
if( mu->Pt() > 10 ) nlep_above_10++; |