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;*/
}
