| 1 |
dkralph |
1.1 |
#include "FR_struct.h"
|
| 2 |
|
|
//----------------------------------------------------------------------------------------
|
| 3 |
|
|
FR_struct::FR_struct(TString mufile, TString elefile)
|
| 4 |
|
|
{
|
| 5 |
|
|
cout << "loading fake file: " << mufile << endl;
|
| 6 |
|
|
// kevin's format
|
| 7 |
|
|
mufrfile = TFile::Open(mufile); assert(mufrfile->IsOpen());
|
| 8 |
|
|
hFRmu = (TH2D*)mufrfile->Get("frEtaPt"); hFRmu->SetDirectory(0);
|
| 9 |
|
|
hFRmuErrl = (TH2D*)mufrfile->Get("errlEtaPt"); hFRmuErrl->SetDirectory(0);
|
| 10 |
|
|
hFRmuErrh = (TH2D*)mufrfile->Get("errhEtaPt"); hFRmuErrh->SetDirectory(0);
|
| 11 |
|
|
mufrfile->Close();
|
| 12 |
|
|
cout << "loading fake file: " << elefile << endl;
|
| 13 |
|
|
elefrfile = TFile::Open(elefile); assert(elefrfile->IsOpen());
|
| 14 |
|
|
hFRele = (TH2D*)elefrfile->Get("frEtaPt"); hFRele->SetDirectory(0);
|
| 15 |
|
|
hFReleErrl = (TH2D*)elefrfile->Get("errlEtaPt"); hFReleErrl->SetDirectory(0);
|
| 16 |
|
|
hFReleErrh = (TH2D*)elefrfile->Get("errhEtaPt"); hFReleErrh->SetDirectory(0);
|
| 17 |
|
|
elefrfile->Close();
|
| 18 |
|
|
// si's format
|
| 19 |
|
|
dynamicMuErrHists = dynamicEleErrHists = false;
|
| 20 |
|
|
if(!hFRmu) {
|
| 21 |
|
|
TH2DAsymErr *hAsymMu = (TH2DAsymErr*)mufrfile->Get("efficiency_denominator_v1_eta_pt");
|
| 22 |
|
|
hFRmu = (TH2D*)hAsymMu;
|
| 23 |
|
|
dynamicMuErrHists = true;
|
| 24 |
|
|
hFRmuErrl = new TH2D(*hFRmu); hFRmuErrl->Reset();
|
| 25 |
|
|
hFRmuErrh = new TH2D(*hFRmu); hFRmuErrh->Reset();
|
| 26 |
|
|
mufr.loadEff(hFRmu,hFRmuErrl,hFRmuErrh);
|
| 27 |
|
|
setErrors(hAsymMu,&mufr);
|
| 28 |
|
|
}
|
| 29 |
|
|
if(!hFRele) {
|
| 30 |
|
|
TH2DAsymErr *hAsymEle = (TH2DAsymErr*)elefrfile->Get("efficiency_denominator_v3_eta_pt");
|
| 31 |
|
|
hFRele = (TH2D*)hAsymEle;
|
| 32 |
|
|
dynamicEleErrHists = true;
|
| 33 |
|
|
hFReleErrl = new TH2D(*hFRele); hFReleErrl->Reset();
|
| 34 |
|
|
hFReleErrh = new TH2D(*hFRele); hFReleErrh->Reset();
|
| 35 |
|
|
elefr.loadEff(hFRele,hFReleErrl,hFReleErrh);
|
| 36 |
|
|
setErrors(hAsymEle,&elefr);
|
| 37 |
|
|
}
|
| 38 |
|
|
mufr.loadEff(hFRmu,hFRmuErrl,hFRmuErrh);
|
| 39 |
|
|
elefr.loadEff(hFRele,hFReleErrl,hFReleErrh);
|
| 40 |
|
|
|
| 41 |
|
|
// figure out whether to use abs eta
|
| 42 |
|
|
if(hFRmu->GetXaxis()->GetBinLowEdge(1) < 0) absEtaMu = false;
|
| 43 |
|
|
else absEtaMu = true;
|
| 44 |
|
|
if(absEtaMu) cout << "\nUsing abs eta in muon fake rates: " << setprecision(3) << hFRmu->GetXaxis()->GetBinLowEdge(1) << endl;
|
| 45 |
|
|
else cout << "\nNot using abs eta in muon fake rates: " << setprecision(3) << hFRmu->GetXaxis()->GetBinLowEdge(1) << endl;
|
| 46 |
|
|
if(hFRele->GetXaxis()->GetBinLowEdge(1) < 0) absEtaEle = false;
|
| 47 |
|
|
else absEtaEle = true;
|
| 48 |
|
|
if(absEtaEle) cout << "\nUsing abs eta in electron fake rates: " << setprecision(3) << hFRele->GetXaxis()->GetBinLowEdge(1) << endl;
|
| 49 |
|
|
else cout << "\nNot using abs eta in electron fake rates: " << setprecision(3) << hFRele->GetXaxis()->GetBinLowEdge(1) << endl;
|
| 50 |
|
|
|
| 51 |
|
|
// mufr.printEff(cout);
|
| 52 |
|
|
// elefr.printEff(cout);
|
| 53 |
|
|
// mufr.printErrLow(cout);
|
| 54 |
|
|
// elefr.printErrLow(cout);
|
| 55 |
|
|
// mufr.printErrHigh(cout);
|
| 56 |
|
|
// elefr.printErrHigh(cout);
|
| 57 |
|
|
}
|
| 58 |
|
|
//----------------------------------------------------------------------------------------
|
| 59 |
|
|
FR_struct::~FR_struct()
|
| 60 |
|
|
{
|
| 61 |
|
|
if(dynamicMuErrHists) {
|
| 62 |
|
|
delete hFRmuErrl;
|
| 63 |
|
|
delete hFRmuErrh;
|
| 64 |
|
|
}
|
| 65 |
|
|
if(dynamicEleErrHists) {
|
| 66 |
|
|
delete hFReleErrl;
|
| 67 |
|
|
delete hFReleErrh;
|
| 68 |
|
|
}
|
| 69 |
|
|
if(mufrfile && mufrfile->IsOpen()) mufrfile->Close();
|
| 70 |
|
|
if(elefrfile && elefrfile->IsOpen()) elefrfile->Close();
|
| 71 |
|
|
}
|
| 72 |
|
|
//----------------------------------------------------------------------------------------
|
| 73 |
|
|
double FR_struct::get_fake_weight(TString uncert, SimpleLepton fake_lep)
|
| 74 |
|
|
{
|
| 75 |
|
|
double rate=0;
|
| 76 |
|
|
assert(uncert=="" || uncert=="lo" || uncert=="hi");
|
| 77 |
|
|
double x = fake_lep.vec.Eta();
|
| 78 |
|
|
double y = fake_lep.vec.Pt();
|
| 79 |
|
|
if(abs(fake_lep.type) == 13) { // fake muons:
|
| 80 |
|
|
if(absEtaMu) x = fabs(x);
|
| 81 |
|
|
rate = mufr.getEff(x,y,true);
|
| 82 |
|
|
if(uncert == "lo") rate -= mufr.getErrLow(x,y,true);
|
| 83 |
|
|
if(uncert == "hi") rate += mufr.getErrHigh(x,y,true);
|
| 84 |
|
|
} else if(abs(fake_lep.type) == 11) {
|
| 85 |
|
|
if(absEtaEle) x = fabs(x);
|
| 86 |
|
|
rate = elefr.getEff(x,y,true);
|
| 87 |
|
|
if(uncert == "lo") rate -= elefr.getErrLow(x,y,true);
|
| 88 |
|
|
if(uncert == "hi") rate += elefr.getErrHigh(x,y,true);
|
| 89 |
|
|
} else { cout << "Bad fake type: " << fake_lep.type << endl; assert(0); }
|
| 90 |
|
|
|
| 91 |
|
|
if(rate==1) cout << "ERROR: rate = 1! " << endl;
|
| 92 |
|
|
|
| 93 |
|
|
double weight = rate/(1-rate);
|
| 94 |
|
|
assert(weight == weight);
|
| 95 |
|
|
return weight;
|
| 96 |
|
|
}
|
| 97 |
|
|
//----------------------------------------------------------------------------------------
|
| 98 |
|
|
void FR_struct::setErrors(TH2DAsymErr *h, CEffUser2D *fr)
|
| 99 |
|
|
{
|
| 100 |
|
|
assert(fr->getHerrl() && fr->getHerrl());
|
| 101 |
|
|
for(int xbin=0; xbin<h->GetXaxis()->GetNbins()+2; xbin++) {
|
| 102 |
|
|
for(int ybin=0; ybin<h->GetYaxis()->GetNbins()+2; ybin++) {
|
| 103 |
|
|
int ncell = h->GetBin(xbin,ybin);
|
| 104 |
|
|
int ncellChk = fr->getHerrl()->GetBin(xbin,ybin);
|
| 105 |
|
|
assert(ncell==ncellChk);
|
| 106 |
|
|
|
| 107 |
|
|
double errl = h->GetBinError(xbin,ybin,TH2DAsymErr::kStatErrLow);
|
| 108 |
|
|
fr->getHerrl()->SetBinContent(ncell,errl);
|
| 109 |
|
|
double errh = h->GetBinError(xbin,ybin,TH2DAsymErr::kStatErrHigh);
|
| 110 |
|
|
fr->getHerrh()->SetBinContent(ncell,errh);
|
| 111 |
|
|
}
|
| 112 |
|
|
}
|
| 113 |
|
|
}
|
| 114 |
|
|
//----------------------------------------------------------------------------------------
|
| 115 |
|
|
void FR_struct::printSiEff(TH2DAsymErr *h, TH2DAsymErr::EErrType errType)
|
| 116 |
|
|
{
|
| 117 |
|
|
const Int_t nx = h->GetNbinsX();
|
| 118 |
|
|
const Int_t ny = h->GetNbinsY();
|
| 119 |
|
|
|
| 120 |
|
|
for(Int_t iy=0; iy<=ny; iy++) {
|
| 121 |
|
|
for(Int_t ix=0; ix<=nx; ix++) {
|
| 122 |
|
|
if(ix==0 && iy==0) {
|
| 123 |
|
|
cout << setw(11) << "";
|
| 124 |
|
|
} else if(ix==0) {
|
| 125 |
|
|
cout << "[" << setw(4) << h->GetYaxis()->GetBinLowEdge(iy) << "," << setw(4) << h->GetYaxis()->GetBinLowEdge(iy+1) << "]";
|
| 126 |
|
|
} else if(iy==0) {
|
| 127 |
|
|
cout << "[" << setw(4) << h->GetXaxis()->GetBinLowEdge(ix) << "," << setw(4) << h->GetXaxis()->GetBinLowEdge(ix+1) << "]";
|
| 128 |
|
|
} else {
|
| 129 |
|
|
ios_base::fmtflags flags = cout.flags();
|
| 130 |
|
|
cout.precision(7);
|
| 131 |
|
|
cout << " " << setw(9) << fixed << h->GetBinError(ix,iy,errType) << " ";
|
| 132 |
|
|
cout.flags(flags);
|
| 133 |
|
|
}
|
| 134 |
|
|
}
|
| 135 |
|
|
cout << endl;
|
| 136 |
|
|
}
|
| 137 |
|
|
}
|
