void k_Vs_relIso() { TFile* file1 = TFile::Open("outputFile_Mu_Run2010A_Sep17thReReco.root"); TH2F* obs = static_cast<TH2F*>(file1->Get("relIso_Vs_d0_tagRateObs")); TH2F* pred = static_cast<TH2F*>(file1->Get("relIso_Vs_d0_tagRatePred")); TFile* file2 = TFile::Open("outputFile_Mu_Run2010B_PromptReco_Part1.root"); obs->Add(static_cast<TH2F*>(file2->Get("relIso_Vs_d0_tagRateObs"))); pred->Add(static_cast<TH2F*>(file2->Get("relIso_Vs_d0_tagRatePred"))); TFile* file3 = TFile::Open("outputFile_Mu_Run2010B_PromptReco_Part2.root"); obs->Add(static_cast<TH2F*>(file3->Get("relIso_Vs_d0_tagRateObs"))); pred->Add(static_cast<TH2F*>(file3->Get("relIso_Vs_d0_tagRatePred"))); TFile* file4 = TFile::Open("outputFile_Mu_Run2010B_PromptReco_Part3.root"); obs->Add(static_cast<TH2F*>(file4->Get("relIso_Vs_d0_tagRateObs"))); pred->Add(static_cast<TH2F*>(file4->Get("relIso_Vs_d0_tagRatePred"))); TH2F* k = new TH2F("k", "k Factor vs relIso", 100, 0.0, 0.25, 100, 0.0, 5.0); k->Divide(obs, pred); k->ProfileY()->Draw(); /*file1->Close();*/ /*delete file1;*/ /*file2->Close();*/ /*delete file2;*/ /*file3->Close();*/ /*delete file3;*/ /*file4->Close();*/ /*delete file4;*/ }