171 |
|
// evaluate MVA *before* regression correction (but cut on mvaVal *after*) |
172 |
|
double mvaVal = getElectronIDMVAval(ctrl, ele, vtx); |
173 |
|
|
174 |
< |
electron_momentum_correction.correct_electron_momentum(ctrl, ele, info, puEnergyDensity->At(0)->RhoKt6PFJets(), vtxArr->GetEntries()); |
174 |
> |
float combination_perr = electron_momentum_correction.correct_electron_momentum(ctrl, ele, info, puEnergyDensity->At(0)->RhoKt6PFJets(), vtxArr->GetEntries()); |
175 |
|
|
176 |
|
if(ctrl.debug) cout << setprecision(8) << " corr el: scEt " << setw(12) << ele->SCluster()->Et() << " P: " |
177 |
|
<< setw(12) << ele->P() << " pT: " << setw(12) << ele->Pt() << " scEta " << setw(12) << ele->SCluster()->Eta() << setprecision(5) << endl; |
200 |
|
tmplep.vec.SetPtEtaPhiM( ele->Pt(), ele->Eta(), ele->Phi(), ELECTRON_MASS ); |
201 |
|
tmplep.type = 11; |
202 |
|
tmplep.index = i; |
203 |
+ |
tmplep.perr = combination_perr; |
204 |
|
tmplep.charge = ele->Charge(); |
205 |
|
tmplep.ip3dSig = ele->Ip3dPVSignificance(); |
206 |
|
tmplep.is4l = false; |