Skip to content

Commit

Permalink
Removing 7TeV electron scale factors so that 8TeV factors are always …
Browse files Browse the repository at this point in the history
…used (until we find/get the correct 7TeV factors). Also corrected a mistake in one of the existing 8TeV corrections.
  • Loading branch information
jjacob committed May 28, 2014
1 parent 64f013b commit e295efa
Showing 1 changed file with 22 additions and 22 deletions.
44 changes: 22 additions & 22 deletions src/RecoObjects/Electron.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -438,31 +438,31 @@ double Electron::getEfficiencyCorrection(bool qcd, int electron_scale_factor_sys
double correction(1.);
double electronEta(eta());

if(Globals::energyInTeV == 7){
// hard coded values since pickle file not available...not updated because I couldn't find any values on twikis!
if (electronEta < -1.5)
correction = 1.003;
else if (electronEta >= -1.5 && electronEta < -1.2)
correction = 0.980;
else if (electronEta >= -1.2 && electronEta < -0.9)
correction = 0.941;
else if (electronEta >= -0.9 && electronEta < 0)
correction = 0.974;
else if (electronEta >= 0 && electronEta < 0.9)
correction = 0.977;
else if (electronEta >= 0.9 && electronEta < 1.2)
correction = 0.939;
else if (electronEta >= 1.2 && electronEta < 1.5)
correction = 0.967;
else if (electronEta >= 1.5)
correction = 1.023;
}
// if(Globals::energyInTeV == 7){
// // hard coded values since pickle file not available...not updated because I couldn't find any values on twikis!
// if (electronEta < -1.5)
// correction = 1.003;
// else if (electronEta >= -1.5 && electronEta < -1.2)
// correction = 0.980;
// else if (electronEta >= -1.2 && electronEta < -0.9)
// correction = 0.941;
// else if (electronEta >= -0.9 && electronEta < 0)
// correction = 0.974;
// else if (electronEta >= 0 && electronEta < 0.9)
// correction = 0.977;
// else if (electronEta >= 0.9 && electronEta < 1.2)
// correction = 0.939;
// else if (electronEta >= 1.2 && electronEta < 1.5)
// correction = 0.967;
// else if (electronEta >= 1.5)
// correction = 1.023;
// }
//8TeV scale factors from https://twiki.cern.ch/twiki/bin/viewauth/CMS/KoPFAElectronTagAndProbe
//Only factors from PromptReco available (in the "Efficiency for e+jet channel (promptreco)" section)
//Specifically: ID & Iso: "ID/Isolation efficiency" sub-section
//Specifically: Trigger: "Trigger efficiency" sub-section
//These values are hard coded because, unlike for Muons, there is no pickle file provided.
else if(qcd == false){ //corrections for (ID & Iso) and Trigger respectively
// else if(qcd == false){ //corrections for (ID & Iso) and Trigger respectively
if(fabs(electronEta)<0.8) {
if(20<=pt() && pt()<30) { //Note: Trigger scale factors only provided down to electron pt of 30GeV in the link above, so I have used the same as for the 30GeV-40GeV range.
switch (electron_scale_factor_systematic) {
Expand Down Expand Up @@ -559,7 +559,7 @@ double Electron::getEfficiencyCorrection(bool qcd, int electron_scale_factor_sys
correction = (0.959+0.003)*(0.988+0.002);
break;
default:
correction = 0.954*0.988;
correction = 0.959*0.988;
}
}
}
Expand Down Expand Up @@ -613,7 +613,7 @@ double Electron::getEfficiencyCorrection(bool qcd, int electron_scale_factor_sys
}
}
}
}
// }
//ID scale factor for QCD is technically irrelevant since we use conversion veto for QCD, but ID and isolation scale factors are provided as one so unable to split them to remove the ID scale factor.
// else if(qcd == true){
// }
Expand Down

0 comments on commit e295efa

Please sign in to comment.