void compare_iso()
{
Sample ttbar, DY;
DY.SetParameters("Drell-Yan");
DY.SetFileNames("/afs/cern.ch/work/i/ishvetso/EgammaWork/my_puppi_test/CMSSW_7_3_3/src/EgammaWork/ElectronNtupler/test/crab_projects/crab_Electron-Isolation_CITK_validation_DY_miniAOD_PUPPI_with_NoLeptons_update11May2015/results/DY.root");
ttbar.SetParameters("ttbar");
ttbar.SetFileNames("/afs/cern.ch/work/i/ishvetso/EgammaWork/my_puppi_test/CMSSW_7_3_3/src/EgammaWork/ElectronNtupler/test/crab_projects/crab_Electron-Isolation_CITK_validation_ttbar_miniAOD_PUPPI_with_NoLeptons_update11May2015/results/ttbar.root");
vector <Var> Vars_CH, Vars_NH, Vars_Gamma;
Var var;
//charged hadrons
var.VarName = "sumChargedHadronPt_CITK";
var.color = kBlue;
var.SetRange(0., 1.5, 30);
Vars_CH.push_back(var);
var.VarName = "sumChargedHadronPt_PUPPI";
var.color = kGreen;
var.SetRange(0., 1.5, 30);
Vars_CH.push_back(var);
var.VarName = "sumChargedHadronPt_PUPPI_NoLeptons";
var.color = kRed;
var.SetRange(0., 1.5, 30);
Vars_CH.push_back(var);
//neutral hadrons
var.VarName = "sumNeutralHadronPt_CITK";
var.color = kBlue;
var.SetRange(0., 1.5, 30);
Vars_NH.push_back(var);
var.VarName = "sumNeutralHadronPt_PUPPI";
var.color = kGreen;
var.SetRange(0., 1.5, 30);
Vars_NH.push_back(var);
var.VarName = "sumNeutralHadronPt_PUPPI_NoLeptons";
var.color = kRed;
var.SetRange(0., 1.5, 30);
Vars_NH.push_back(var);
//photons
var.VarName = "sumPhotonPt_CITK";
var.color = kBlue;
var.SetRange(0., 1.5, 30);
Vars_Gamma.push_back(var);
var.VarName = "sumPhotonPt_PUPPI";
var.color = kGreen;
var.SetRange(0., 1.5, 30);
Vars_Gamma.push_back(var);
var.VarName = "sumPhotonPt_PUPPI_NoLeptons";
var.color = kRed;
var.SetRange(0., 1.5, 30);
Vars_Gamma.push_back(var);
setTDRStyle();
draw(Vars_CH, ttbar, "ch", "barrel", "isEB == 1", "iso_matching");
draw(Vars_CH, ttbar, "ch", "endcap", "isEB != 1", "iso_matching");
draw(Vars_CH, DY, "ch", "endcap", "isEB != 1", "iso_matching");
draw(Vars_CH, DY, "ch", "barrel", "isEB == 1", "iso_matching");
draw(Vars_NH, ttbar, "nh", "barrel", "isEB == 1", "iso_matching");
draw(Vars_NH, ttbar, "nh", "endcap", "isEB != 1", "iso_matching");
draw(Vars_NH, DY, "nh", "endcap", "isEB != 1", "iso_matching");
draw(Vars_NH, DY, "nh", "barrel", "isEB == 1", "iso_matching");
draw(Vars_Gamma , ttbar, "gamma", "barrel", "isEB == 1", "iso_matching");
draw(Vars_Gamma, ttbar, "gamma", "endcap", "isEB != 1", "iso_matching");
draw(Vars_Gamma, DY, "gamma", "endcap", "isEB != 1", "iso_matching");
draw(Vars_Gamma, DY, "gamma", "barrel", "isEB == 1", "iso_matching");
}
