void CutFlow::allPlots(AllSamples samples){ // Variable::Variable(TString name_temp, TString xTitle_temp, double minX_temp, double maxX_temp, int rebinFact_temp) Variable ttbar("TTbarMuPlusJetsRefSelection", "selection step", -1.5, 9.5, 1); saveCutFlowPlot(samples, ttbar); }
void manualFRTest() { TFile* f = TFile::Open("/eos/user/v/vischia/www/taufakes_new/plotter_wjet.root", "READ"); TString data("data (single mu)"), ttbar("TTbar"), wjets("W#rightarrow l#nu inclusive"), qcd("QCD_EMEnr_Incl"); TString numerator("/wjet_step5byLooseCombinedIsolationDeltaBetaCorr3Hitseta_numerator"), denominator("/wjet_step5eta_denominator"); TH1F* data_num = (TH1F*) f->Get(data+numerator); TH1F* data_den = (TH1F*) f->Get(data+denominator); TH1F* ttbar_num = (TH1F*) f->Get(ttbar+numerator); TH1F* ttbar_den = (TH1F*) f->Get(ttbar+denominator); TH1F* wjets_num = (TH1F*) f->Get(wjets+numerator); TH1F* wjets_den = (TH1F*) f->Get(wjets+denominator); TH1F* qcd_num = (TH1F*) f->Get(qcd+numerator); TH1F* qcd_den = (TH1F*) f->Get(qcd+denominator); TH1F* mc_num = (TH1F*) ttbar_num->Clone("mc_num"); mc_num->Add(wjets_num); //mc_num->Add(qcd_num); TH1F* mc_den = (TH1F*) ttbar_den->Clone("mc_den"); mc_den->Add(wjets_den); //mc_den->Add(qcd_den); data_num->SetMarkerColor(kBlack); mc_num->SetMarkerColor(kRed); data_den->SetMarkerColor(kBlack); mc_den->SetMarkerColor(kRed); TCanvas* d = new TCanvas("d", "d", 600, 1200); d->Divide(1,2); d->cd(1); gPad->SetTitle("numerator"); data_num->DrawCopy("pe"); mc_num->DrawCopy("pesame"); d->cd(2); gPad->SetTitle("denominator"); data_den->DrawCopy("pe"); mc_den->DrawCopy("pesame"); TH1F* data_fr = (TH1F*) data_num->Clone("data_fr"); data_fr->Divide(data_den); TH1F* mc_fr = (TH1F*) mc_num->Clone("mc_fr"); mc_fr->Divide(mc_den); data_fr->SetMarkerColor(kBlack); mc_fr->SetMarkerColor(kRed); data_fr->SetMarkerSize(2); mc_fr->SetMarkerSize(2); TCanvas* c = new TCanvas("c", "c", 600, 600); c->cd(); data_fr->Draw("pe"); mc_fr->Draw("pesame"); }