NtupleAnalysis: MyHtautauAnalysis/src/HtautauHistos.cxx Source File

00001 //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
00002 // 06.12.2006, AUTHOR: MANFRED GROH
00003 //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
00004 
00006 // BOOK THE HISTOGRAM //
00008 
00009     char char_help_array[50];
00010 
00011     // for normal histos
00012     double ptmin(0.0);
00013     double ptmax(250.0);
00014     double etamin(-5.0);
00015     double etamax(5.0);
00016     double phimin(-5.0);
00017     double phimax(5.0);
00018     // for delta histos, e.g. delta eta of two particles
00019     double d_ptmax(50.0);
00020     double d_etamax(0.2);
00021     double d_phimax(0.2);
00022 
00023     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00024         dir_truth_decay = new TDirectoryFile("truth_decay","truth_decay");
00025     #else
00026         dir_truth_decay = new TDirectory("truth_decay","truth_decay");
00027     #endif
00028     //dir_truth_decay = new TDirectoryFile("truth_decay","truth_decay");
00029     dir_truth_decay->cd();
00030 
00031     // truth leptons (check if lepton filter was applied)
00032     h_emu_truth_filter_nb = new TH1F("emu_truth_filter_nb","truth e+mu, pt>5 GeV, eta<2.7",19, -1, 18);
00033     h_TTbarLeptonFilter = new TH1F("TTbarLeptonFilter","0:rej, 1: acc, 2:jec*weight, 3:acc*weight",10, 0., 10.);
00034 
00035     // truth decay mode histos
00036     h_primary_particle_mode = new
00037         TH1F("primary_particle_mode","Primary Particle Mode", 20, -1.5, 18.5);
00038     h_primary_particle_decay_mode = new
00039         TH1F("primary_particle_decay_mode","Primary Particle Decay Mode", 20, -1.5, 18.5);
00040     h_tau_0_decay_mode = new
00041         TH1F("tau_0_decay_mode","Tau 0 Decay Mode", 20, -1.5, 18.5);
00042     h_tau_1_decay_mode = new
00043         TH1F("tau_1_decay_mode","Tau 1 Decay Mode", 20, -1.5, 18.5);
00044     h_higgs_decay_mode = new
00045             TH1F("higgs_decay_mode","Higgs Decay Mode", 20, -1.5, 18.5);
00046     h_truth_taudecay_dr = new TH1F("truth_taudecay_dr", "truth_taudecay_dr",100, 0., 0.5);
00047     h_truth_antitaudecay_dr = new TH1F("truth_antitaudecay_dr", "truth_antitaudecay_dr",100, 0., 0.5);
00048     h_truth_2l_rf_dr = new TH1F("truth_2l_rf_dr", "truth_2l_rf_dr",50, 0., 6.);
00049     dir_truth_decay->cd("../");
00050 
00051 // raw aod content histos
00052     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00053         dir_aod_content = new TDirectoryFile("aod_content","aod_content");
00054     #else
00055         dir_aod_content = new TDirectory("aod_content","aod_content");
00056     #endif
00057     //dir_aod_content = new TDirectory("aod_content","aod_content");
00058     dir_aod_content->cd();
00059     h_etmiss_aod = new TH1F("etmiss_aod","etmiss_aod",100, ptmin, ptmax);
00060     h_e_reco_aod_nb = new TH1F("e_reco_aod_nb","e_reco_aod_nb",19, -1, 18);
00061     h_e_reco_aod_pt = new TH1F("e_reco_aod_pt","e_reco_aod_pt",100, ptmin, ptmax);
00062     h_e_reco_aod_eta = new TH1F("e_reco_aod_eta","e_reco_aod_eta",100, etamin, etamax);
00063     h_e_reco_aod_phi = new TH1F("e_reco_aod_phi","e_reco_aod_phi",100, phimin, phimax);
00064     h_e_reco_aod_isolation02 = new TH2D("h_e_reco_aod_isolation02","h_e_reco_aod_isolation02",104,-2.6,2.6,45,-5.,40.);
00065     h_mu_reco_aod_nb = new TH1F("mu_reco_aod_nb","mu_reco_aod_nb",19, -1, 18);
00066     h_mu_reco_aod_pt = new TH1F("mu_reco_aod_pt","mu_reco_aod_pt",100, ptmin, ptmax);
00067     h_mu_reco_aod_eta = new TH1F("mu_reco_aod_eta","mu_reco_aod_eta",100, etamin, etamax);
00068     h_mu_reco_aod_phi = new TH1F("mu_reco_aod_phi","mu_reco_aod_phi",100, phimin, phimax);
00069     h_mu_reco_aod_fitChi2OverDoF = new TH1F("mu_reco_aod_fitChi2OverDoF","mu_reco_aod_fitChi2OverDoF",100, -1, 9);
00070     h_taujet_reco_aod_nb = new TH1F("taujet_reco_aod_nb","taujet_reco_aod_nb",19, -1, 18);
00071     h_taujet_reco_aod_pt = new TH1F("taujet_reco_aod_pt","taujet_reco_aod_pt",100, ptmin, ptmax);
00072     h_taujet_reco_aod_eta = new TH1F("taujet_reco_aod_eta","taujet_reco_aod_eta",100, etamin, etamax);
00073     h_taujet_reco_aod_phi = new TH1F("taujet_reco_aod_phi","taujet_reco_aod_phi",100, phimin, phimax);
00074     h_jet_reco_aod_nb = new TH1F("jet_reco_aod_nb","jet_reco_aod_nb",19, -1, 18);
00075     h_jet_reco_aod_pt = new TH1F("jet_reco_aod_pt","jet_reco_aod_pt",100, ptmin, ptmax);
00076     h_jet_reco_aod_eta = new TH1F("jet_reco_aod_eta","jet_reco_aod_eta",100, etamin, etamax);
00077     h_jet_reco_aod_phi = new TH1F("jet_reco_aod_phi","jet_reco_aod_phi",100, phimin, phimax);
00078     h_e_truth_aod_nb = new TH1F("e_truth_aod_nb","e_truth_aod_nb",19, -1, 18);
00079     h_e_truth_aod_pt = new TH1F("e_truth_aod_pt","e_truth_aod_pt",100, ptmin, ptmax);
00080     h_e_truth_aod_eta = new TH1F("e_truth_aod_eta","e_truth_aod_eta",100, etamin, etamax);
00081     h_e_truth_aod_phi = new TH1F("e_truth_aod_phi","e_truth_aod_phi",100, phimin, phimax);
00082     h_mu_truth_aod_nb = new TH1F("mu_truth_aod_nb","mu_truth_aod_nb",19, -1, 18);
00083     h_mu_truth_aod_pt = new TH1F("mu_truth_aod_pt","mu_truth_aod_pt",100, ptmin, ptmax);
00084     h_mu_truth_aod_eta = new TH1F("mu_truth_aod_eta","mu_truth_aod_eta",100, etamin, etamax);
00085     h_mu_truth_aod_phi = new TH1F("mu_truth_aod_phi","mu_truth_aod_phi",100, phimin, phimax);
00086     h_tau_truth_aod_nb = new TH1F("tau_truth_aod_nb","tau_truth_aod_nb",19, -1, 18);
00087     h_tau_truth_aod_pt = new TH1F("tau_truth_aod_pt","tau_truth_aod_pt",100, ptmin, ptmax);
00088     h_tau_truth_aod_eta = new TH1F("tau_truth_aod_eta","tau_truth_aod_eta",100, etamin, etamax);
00089     h_tau_truth_aod_phi = new TH1F("tau_truth_aod_phi","tau_truth_aod_phi",100, phimin, phimax);
00090     h_tau_truth_aod_status = new TH1F("tau_truth_aod_status","tau_truth_aod_status",200, -1, 198);
00091     h_jet_truth_aod_nb = new TH1F("jet_truth_aod_nb","jet_truth_aod_nb",19, -1, 18);
00092     h_jet_truth_aod_pt = new TH1F("jet_truth_aod_pt","jet_truth_aod_pt",100, ptmin, ptmax);
00093     h_jet_truth_aod_eta = new TH1F("jet_truth_aod_eta","jet_truth_aod_eta",100, etamin, etamax);
00094     h_jet_truth_aod_phi = new TH1F("jet_truth_aod_phi","jet_truth_aod_phi",100, phimin, phimax);
00095     dir_aod_content->cd("../");
00096 
00097 
00098 // matching, efficiency histos
00099     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00100         dir_presel_matching = new TDirectoryFile("presel_matching","presel_matching");
00101     #else
00102         dir_presel_matching = new TDirectory("presel_matching","presel_matching");
00103     #endif
00104     //dir_presel_matching = new TDirectory("presel_matching","presel_matching");
00105     dir_presel_matching->cd();
00106     // electrons
00107     h_e_reco_presel_pt = new TH1F("e_reco_presel_pt","e_reco_presel_pt",100, ptmin, ptmax);
00108     h_e_reco_presel_eta = new TH1F("e_reco_presel_eta","e_reco_presel_eta",100, etamin, etamax);
00109     h_e_reco_presel_phi = new TH1F("e_reco_presel_phi","e_reco_presel_phi",100, phimin, phimax);
00110     h_e_reco_presel_nb = new TH1F("e_reco_presel_nb","e_reco_presel_nb",19, -1, 18);
00111     h_e_reco_presel_isem = new TH1F("e_reco_presel_isem","e_reco_presel_isem",4001,0.0,4000.0);
00112 
00113     h_e_truth_presel_pt = new TH1F("e_truth_presel_pt","e_truth_presel_pt",100, ptmin, ptmax);
00114     h_e_truth_presel_eta = new TH1F("e_truth_presel_eta","e_truth_presel_eta",100, etamin, etamax);
00115     h_e_truth_presel_phi = new TH1F("e_truth_presel_phi","e_truth_presel_phi",100, phimin, phimax);
00116     h_e_truth_presel_nb = new TH1F("e_truth_presel_nb","e_truth_presel_nb",19, -1, 18);
00117     h_e_truth_presel_status = new TH1F("e_truth_presel_status",
00118         "e_truth_presel_status", 200, -1, 198);
00119 
00120     h_e_reco_matched_pt = new TH1F("e_reco_matched_pt","e_reco_matched_pt",100, ptmin, ptmax);
00121     h_e_reco_matched_eta = new TH1F("e_reco_matched_eta","e_reco_matched_eta",100, etamin, etamax);
00122     h_e_reco_matched_phi = new TH1F("e_reco_matched_phi","e_reco_matched_phi",100, phimin, phimax);
00123     h_e_reco_matched_nb = new TH1F("e_reco_matched_nb","e_reco_matched_nb",19, -1, 18);
00124     h_e_reco_matched_isem = new TH1F("e_reco_matched_isem","e_reco_matched_isem",4001,0.0,4000.0);
00125     h_e_reco_matched_NN = new TH1F("e_reco_matched_NN","e_reco_matched_NN",100,-0.1,1.1);
00126 
00127     h_e_truth_matched_pt = new TH1F("e_truth_matched_pt","e_truth_matched_pt",100, ptmin, ptmax);
00128     h_e_truth_matched_eta = new TH1F("e_truth_matched_eta","e_truth_matched_eta",100, etamin, etamax);
00129     h_e_truth_matched_phi = new TH1F("e_truth_matched_phi","e_truth_matched_phi",100, phimin, phimax);
00130 
00131     h_e_reco_fake_pt = new TH1F("e_reco_fake_pt","e_reco_fake_pt",100, ptmin, ptmax);
00132     h_e_reco_fake_eta = new TH1F("e_reco_fake_eta","e_reco_fake_eta",100, etamin, etamax);
00133     h_e_reco_fake_phi = new TH1F("e_reco_fake_phi","e_reco_fake_phi",100, phimin, phimax);
00134     h_e_reco_fake_isem = new TH1F("e_reco_fake_isem","e_reco_fake_isem",4001,0.0,4000.0);
00135     h_e_reco_fake_nb = new TH1F("e_reco_fake_nb","e_reco_fake_nb",19, -1, 18);
00136 
00137     h_e_truth_notmatched_pt = new TH1F("e_truth_notmatched_pt","e_truth_notmatched_pt", 100, ptmin, ptmax);
00138     h_e_truth_notmatched_eta = new TH1F("e_truth_notmatched_eta","e_truth_notmatched_eta", 100, etamin, etamax);
00139     h_e_truth_notmatched_phi = new TH1F("e_truth_notmatched_phi","e_truth_notmatched_phi", 100, phimin, phimax);
00140     h_e_truth_notmatched_nb = new TH1F("e_truth_notmatched_nb","e_truth_notmatched_nb",19, -1, 18);
00141 
00142     h_e_matched_dpt = new TH1F("e_matched_dpt","e_matched_dpt", 200, -d_ptmax, d_ptmax);
00143     h_e_matched_deta = new TH1F("e_matched_deta","e_matched_deta", 200, -d_etamax, d_etamax);
00144     h_e_matched_dphi = new TH1F("e_matched_dphi","e_matched_dphi", 200, -d_phimax, d_phimax);
00145     h_e_matched_dr = new TH1F("e_matched_dr","e_matched_dr", 200, 0.0, 0.1);
00146 
00147     h_e_matched_resolution_pt = new TH2F("e_matched_resolution_pt","e_matched_resolution_pt",
00148                      20, ptmin, ptmax, 100, -1.0, 1.0);
00149     h_e_matched_resolution_eta = new TH2F("e_matched_resolution_eta","e_matched_resolution_eta",
00150                      20, etamin, etamax, 100, -1.0, 1.0);
00151     h_e_matched_resolution_phi = new TH2F("e_matched_resolution_phi","e_matched_resolution_phi",
00152                      20, phimin, phimax, 100, -1.0, 1.0);
00153 
00154     // muons
00155     h_mu_reco_presel_pt = new TH1F("mu_reco_presel_pt","mu_reco_presel_pt",100, ptmin, ptmax);
00156     h_mu_reco_presel_eta = new TH1F("mu_reco_presel_eta","mu_reco_presel_eta",100, etamin, etamax);
00157     h_mu_reco_presel_phi = new TH1F("mu_reco_presel_phi","mu_reco_presel_phi",100, phimin, phimax);
00158     h_mu_reco_presel_nb = new TH1F("mu_reco_presel_nb","mu_reco_presel_nb",19, -1, 18);
00159     h_mu_reco_presel_fitChi2OverDoF = new TH1F("mu_reco_presel_fitChi2OverDoF","mu_reco_presel_fitChi2OverDoF",100, -1, 9);
00160     h_mu_reco_presel_fitChi2OverDoFvsphi = new TH2F("mu_reco_presel_fitChi2OverDoFvsphi","mu_reco_presel_fitChi2OverDoFvsphi",
00161                      40, phimin, phimax, 100, -1.0, 9.0);
00162 
00163     h_mu_truth_presel_pt = new TH1F("mu_truth_presel_pt","mu_truth_presel_pt",100, ptmin, ptmax);
00164     h_mu_truth_presel_eta = new TH1F("mu_truth_presel_eta","mu_truth_presel_eta",100, etamin, etamax);
00165     h_mu_truth_presel_phi = new TH1F("mu_truth_presel_phi","mu_truth_presel_phi",100, phimin, phimax);
00166     h_mu_truth_presel_nb = new TH1F("mu_truth_presel_nb","mu_truth_presel_nb",19, -1, 18);
00167     h_mu_truth_presel_status = new TH1F("mu_truth_presel_status",
00168         "mu_truth_presel_status", 100, -5, 95);
00169 
00170     h_mu_reco_matched_pt = new TH1F("mu_reco_matched_pt","mu_reco_matched_pt",100, ptmin, ptmax);
00171     h_mu_reco_matched_eta = new TH1F("mu_reco_matched_eta","mu_reco_matched_eta",100, etamin, etamax);
00172     h_mu_reco_matched_phi = new TH1F("mu_reco_matched_phi","mu_reco_matched_phi",100, phimin, phimax);
00173     h_mu_reco_matched_nb = new TH1F("mu_reco_matched_nb","mu_reco_matched_nb",19, -1, 18);
00174     h_mu_reco_matched_fitChi2OverDoF = new TH1F("mu_reco_matched_fitChi2OverDoF","mu_reco_matched_fitChi2OverDoF",100, -1, 9);
00175 
00176     h_mu_truth_matched_pt = new TH1F("mu_truth_matched_pt","mu_truth_matched_pt",100, ptmin, ptmax);
00177     h_mu_truth_matched_eta = new TH1F("mu_truth_matched_eta","mu_truth_matched_eta",100, etamin, etamax);
00178     h_mu_truth_matched_phi = new TH1F("mu_truth_matched_phi","mu_truth_matched_phi",100, phimin, phimax);
00179 
00180     h_mu_reco_fake_pt = new TH1F("mu_reco_fake_pt","mu_reco_fake_pt",100, ptmin, ptmax);
00181     h_mu_reco_fake_eta = new TH1F("mu_reco_fake_eta","mu_reco_fake_eta",100, etamin, etamax);
00182     h_mu_reco_fake_phi = new TH1F("mu_reco_fake_phi","mu_reco_fake_phi",100, phimin, phimax);
00183     h_mu_reco_fake_nb = new TH1F("mu_reco_fake_nb","mu_reco_fake_nb",19, -1, 18);
00184     h_mu_reco_fake_fitChi2OverDoF = new TH1F("mu_reco_fake_fitChi2OverDoF","mu_reco_fake_fitChi2OverDoF",100, -1, 9);
00185 
00186     h_mu_truth_notmatched_pt = new TH1F("mu_truth_notmatched_pt","mu_truth_notmatched_pt", 100, ptmin, ptmax);
00187     h_mu_truth_notmatched_eta = new TH1F("mu_truth_notmatched_eta","mu_truth_notmatched_eta", 100, etamin, etamax);
00188     h_mu_truth_notmatched_phi = new TH1F("mu_truth_notmatched_phi","mu_truth_notmatched_phi", 100, phimin, phimax);
00189     h_mu_truth_notmatched_nb = new TH1F("mu_truth_notmatched_nb","mu_truth_notmatched_nb",19, -1, 18);
00190 
00191     h_mu_matched_dpt = new TH1F("mu_matched_dpt","mu_matched_dpt", 200, -d_ptmax, d_ptmax);
00192     h_mu_matched_deta = new TH1F("mu_matched_deta","mu_matched_deta", 200, -d_etamax, d_etamax);
00193     h_mu_matched_dphi = new TH1F("mu_matched_dphi","mu_matched_dphi", 200, -d_phimax, d_phimax);
00194     h_mu_matched_dr = new TH1F("mu_matched_dr","mu_matched_dr", 200, 0.0, 0.1);
00195 
00196     h_mu_matched_resolution_pt = new TH2F("mu_matched_resolution_pt","mu_matched_resolution_pt",
00197                      20, ptmin, ptmax, 100, -1.0, 1.0);
00198     h_mu_matched_resolution_ptvsphi = new TH2F("mu_matched_resolution_ptvsphi","mu_matched_resolution_ptvsphi",
00199                      40, phimin, phimax, 100, -1.0, 1.0);
00200     h_mu_matched_resolution_ptvseta = new TH2F("mu_matched_resolution_ptvseta","mu_matched_resolution_ptvseta",
00201                      20, etamin, etamax, 100, -1.0, 1.0);
00202     h_mu_matched_resolution_eta = new TH2F("mu_matched_resolution_eta","mu_matched_resolution_eta",
00203                      20, etamin, etamax, 100, -1.0, 1.0);
00204     h_mu_matched_resolution_phi = new TH2F("mu_matched_resolution_phi","mu_matched_resolution_phi",
00205                      20, phimin, phimax, 100, -1.0, 1.0);
00206 
00207     // taujets
00208     h_taujet_reco_presel_llh = new TH1F("taujet_reco_presel_llh","taujet_reco_presel_llh",101, -50.0, 50.0);
00209     h_taujet_reco_presel_pt = new TH1F("taujet_reco_presel_pt","taujet_reco_presel_pt",100, ptmin, ptmax);
00210     h_taujet_reco_presel_eta = new TH1F("taujet_reco_presel_eta","taujet_reco_presel_eta",100, etamin, etamax);
00211     h_taujet_reco_presel_phi = new TH1F("taujet_reco_presel_phi","taujet_reco_presel_phi",100, phimin, phimax);
00212     h_taujet_reco_presel_nb = new TH1F("taujet_reco_presel_nb","taujet_reco_presel_nb",19, -1, 18);
00213 
00214     h_taujet_truth_presel_pt = new TH1F("taujet_truth_presel_pt","taujet_truth_presel_pt",100, ptmin, ptmax);
00215     h_taujet_truth_presel_eta = new TH1F("taujet_truth_presel_eta","taujet_truth_presel_eta",100, etamin, etamax);
00216     h_taujet_truth_presel_phi = new TH1F("taujet_truth_presel_phi","taujet_truth_presel_phi",100, phimin, phimax);
00217     h_taujet_truth_presel_nb = new TH1F("taujet_truth_presel_nb","taujet_truth_presel_nb",19, -1, 18);
00218     h_taujet_truth_presel_status = new TH1F("taujet_truth_presel_status",
00219         "taujet_truth_presel_status", 200, -1, 198);
00220 
00221     h_taujet_reco_matched_pt = new TH1F("taujet_reco_matched_pt","taujet_reco_matched_pt",100, ptmin, ptmax);
00222     h_taujet_reco_matched_eta = new TH1F("taujet_reco_matched_eta","taujet_reco_matched_eta",100, etamin, etamax);
00223     h_taujet_reco_matched_phi = new TH1F("taujet_reco_matched_phi","taujet_reco_matched_phi",100, phimin, phimax);
00224     h_taujet_reco_matched_nb = new TH1F("taujet_reco_matched_nb","taujet_reco_matched_nb",19, -1, 18);
00225 
00226     h_taujet_truth_matched_pt = new TH1F("taujet_truth_matched_pt","taujet_truth_matched_pt",100, ptmin, ptmax);
00227     h_taujet_truth_matched_eta = new TH1F("taujet_truth_matched_eta","taujet_truth_matched_eta",100, etamin, etamax);
00228     h_taujet_truth_matched_phi = new TH1F("taujet_truth_matched_phi","taujet_truth_matched_phi",100, phimin, phimax);
00229 
00230     h_taujet_reco_fake_pt = new TH1F("taujet_reco_fake_pt","taujet_reco_fake_pt",100, ptmin, ptmax);
00231     h_taujet_reco_fake_eta = new TH1F("taujet_reco_fake_eta","taujet_reco_fake_eta",100, etamin, etamax);
00232     h_taujet_reco_fake_phi = new TH1F("taujet_reco_fake_phi","taujet_reco_fake_phi",100, phimin, phimax);
00233     h_taujet_reco_fake_nb = new TH1F("taujet_reco_fake_nb","taujet_reco_fake_nb",19, -1, 18);
00234 
00235     h_taujet_truth_notmatched_pt = new TH1F("taujet_truth_notmatched_pt","taujet_truth_notmatched_pt", 100, ptmin, ptmax);
00236     h_taujet_truth_notmatched_eta = new TH1F("taujet_truth_notmatched_eta","taujet_truth_notmatched_eta", 100, etamin, etamax);
00237     h_taujet_truth_notmatched_phi = new TH1F("taujet_truth_notmatched_phi","taujet_truth_notmatched_phi", 100, phimin, phimax);
00238     h_taujet_truth_notmatched_nb = new TH1F("taujet_truth_notmatched_nb","taujet_truth_notmatched_nb",19, -1, 18);
00239 
00240     h_taujet_matched_dpt = new TH1F("taujet_matched_dpt","taujet_matched_dpt", 200, -d_ptmax, d_ptmax);
00241     h_taujet_matched_deta = new TH1F("taujet_matched_deta","taujet_matched_deta", 200, -d_etamax, d_etamax);
00242     h_taujet_matched_dphi = new TH1F("taujet_matched_dphi","taujet_matched_dphi", 200, -d_phimax, d_phimax);
00243     h_taujet_matched_dr = new TH1F("taujet_matched_dr","taujet_matched_dr", 200, 0.0, 0.1);
00244 
00245     h_taujet_matched_resolution_pt = new TH2F("taujet_matched_resolution_pt","taujet_matched_resolution_pt",
00246                      20, ptmin, ptmax, 100, -1.0, 1.0);
00247     h_taujet_matched_resolution_eta = new TH2F("taujet_matched_resolution_eta","taujet_matched_resolution_eta",
00248                      20, etamin, etamax, 100, -1.0, 1.0);
00249     h_taujet_matched_resolution_phi = new TH2F("taujet_matched_resolution_phi","taujet_matched_resolution_phi",
00250                      20, phimin, phimax, 100, -1.0, 1.0);
00251     h_taujet_fake_match = new TH1F("taujet_fake_match","taujet_fake_match",10, 0, 10);
00252     
00253 
00254     h_jet_faking_tau_pt = new TH1F("jet_faking_tau_pt","jet_faking_tau_pt",100, ptmin, ptmax);
00255     h_jet_faking_tau_eta = new TH1F("jet_faking_tau_eta","jet_faking_tau_eta",100, etamin, etamax);
00256     h_jet_faking_tau_phi = new TH1F("jet_faking_tau_phi","jet_faking_tau_phi",100, phimin, phimax);
00257     h_jet_faking_tau_nb = new TH1F("jet_faking_tau_nb","jet_faking_tau_nb",19, -1, 18);
00258 
00259 
00260     h_tau_all_truth_nb = new TH1F("tau_all_truth_nb","tau_all_truth_nb",19, -1, 18);
00261     h_tau_cleaned_truth_nb = new TH1F("tau_cleaned_truth_nb","tau_cleaned_truth_nb",19, -1, 18);
00262     h_taujet_cleaned_truth_nb = new TH1F("taujet_cleaned_truth_nb","taujet_cleaned_truth_nb",19, -1, 18);
00263     h_tau_all_truth_status = new TH1F("tau_all_truth_status","tau_all_truth_status", 200, -1, 198);
00264     h_tau_cleaned_truth_status = new TH1F("tau_cleaned_truth_status","tau_cleaned_truth_status", 200, -1, 198);
00265     h_taujet_cleaned_truth_status = new TH1F("taujet_cleaned_truth_status","taujet_cleaned_truth_status", 200, -1, 198);
00266 
00267 
00268     // bjets
00269     h_bjet_reco_presel_pt = new TH1F("bjet_reco_presel_pt","bjet_reco_presel_pt",100, ptmin, ptmax);
00270     h_bjet_reco_presel_eta = new TH1F("bjet_reco_presel_eta","bjet_reco_presel_eta",100, etamin, etamax);
00271     h_bjet_reco_presel_phi = new TH1F("bjet_reco_presel_phi","bjet_reco_presel_phi",100, phimin, phimax);
00272     h_bjet_reco_presel_nb = new TH1F("bjet_reco_presel_nb","bjet_reco_presel_nb",19, -1, 18);
00273 
00274     h_bjet_truth_presel_pt = new TH1F("bjet_truth_presel_pt","bjet_truth_presel_pt",100, ptmin, ptmax);
00275     h_bjet_truth_presel_eta = new TH1F("bjet_truth_presel_eta","bjet_truth_presel_eta",100, etamin, etamax);
00276     h_bjet_truth_presel_phi = new TH1F("bjet_truth_presel_phi","bjet_truth_presel_phi",100, phimin, phimax);
00277     h_bjet_truth_presel_nb = new TH1F("bjet_truth_presel_nb","bjet_truth_presel_nb",19, -1, 18);
00278     
00279     h_bjet_reco_matched_pt = new TH1F("bjet_reco_matched_pt","bjet_reco_matched_pt",100, ptmin, ptmax);
00280     h_bjet_reco_matched_eta = new TH1F("bjet_reco_matched_eta","bjet_reco_matched_eta",100, etamin, etamax);
00281     h_bjet_reco_matched_phi = new TH1F("bjet_reco_matched_phi","bjet_reco_matched_phi",100, phimin, phimax);
00282     h_bjet_reco_matched_nb = new TH1F("bjet_reco_matched_nb","bjet_reco_matched_nb",19, -1, 18);
00283 
00284     h_bjet_truth_matched_pt = new TH1F("bjet_truth_matched_pt","bjet_truth_matched_pt",100, ptmin, ptmax);
00285     h_bjet_truth_matched_eta = new TH1F("bjet_truth_matched_eta","bjet_truth_matched_eta",100, etamin, etamax);
00286     h_bjet_truth_matched_phi = new TH1F("bjet_truth_matched_phi","bjet_truth_matched_phi",100, phimin, phimax);
00287 
00288     h_bjet_reco_fake_pt = new TH1F("bjet_reco_fake_pt","bjet_reco_fake_pt",100, ptmin, ptmax);
00289     h_bjet_reco_fake_eta = new TH1F("bjet_reco_fake_eta","bjet_reco_fake_eta",100, etamin, etamax);
00290     h_bjet_reco_fake_phi = new TH1F("bjet_reco_fake_phi","bjet_reco_fake_phi",100, phimin, phimax);
00291     h_bjet_reco_fake_nb = new TH1F("bjet_reco_fake_nb","bjet_reco_fake_nb",19, -1, 18);
00292 
00293     h_bjet_truth_notmatched_pt = new TH1F("bjet_truth_notmatched_pt","bjet_truth_notmatched_pt", 100, ptmin, ptmax);
00294     h_bjet_truth_notmatched_eta = new TH1F("bjet_truth_notmatched_eta","bjet_truth_notmatched_eta", 100, etamin, etamax);
00295     h_bjet_truth_notmatched_phi = new TH1F("bjet_truth_notmatched_phi","bjet_truth_notmatched_phi", 100, phimin, phimax);
00296     h_bjet_truth_notmatched_nb = new TH1F("bjet_truth_notmatched_nb","bjet_truth_notmatched_nb",19, -1, 18);
00297 
00298     h_bjet_matched_dpt = new TH1F("bjet_matched_dpt","bjet_matched_dpt", 200, -d_ptmax, d_ptmax);
00299     h_bjet_matched_deta = new TH1F("bjet_matched_deta","bjet_matched_deta", 200, -d_etamax, d_etamax);
00300     h_bjet_matched_dphi = new TH1F("bjet_matched_dphi","bjet_matched_dphi", 200, -d_phimax, d_phimax);
00301     h_bjet_matched_dr = new TH1F("bjet_matched_dr","bjet_matched_dr", 200, 0.0, 0.1);
00302 
00303     h_bjet_matched_resolution_pt = new TH2F("bjet_matched_resolution_pt","bjet_matched_resolution_pt",
00304                                             20, ptmin, ptmax, 100, -1.0, 1.0);
00305     h_bjet_matched_resolution_eta = new TH2F("bjet_matched_resolution_eta","bjet_matched_resolution_eta",
00306                                              20, etamin, etamax, 100, -1.0, 1.0);
00307     h_bjet_matched_resolution_phi = new TH2F("bjet_matched_resolution_phi","bjet_matched_resolution_phi",
00308                                              20, phimin, phimax, 100, -1.0, 1.0);
00309 
00310     
00311 
00312     // jets
00313     // after overlap removal with electrons
00314     h_jet_reco_orm_pt = new TH1F("jet_reco_orm_pt","jet_reco_orm_pt",100, ptmin, ptmax);
00315     h_jet_reco_orm_eta = new TH1F("jet_reco_orm_eta","jet_reco_orm_eta",100, etamin, etamax);
00316     h_jet_reco_orm_phi = new TH1F("jet_reco_orm_phi","jet_reco_orm_phi",100, phimin, phimax);
00317     h_jet_reco_orm_nb = new TH1F("jet_reco_orm_nb","jet_reco_orm_nb",19, -1, 18);
00318 
00319     h_jet_reco_presel_pt = new TH1F("jet_reco_presel_pt","jet_reco_presel_pt",100, ptmin, ptmax);
00320     h_jet_reco_presel_eta = new TH1F("jet_reco_presel_eta","jet_reco_presel_eta",100, etamin, etamax);
00321     h_jet_reco_presel_phi = new TH1F("jet_reco_presel_phi","jet_reco_presel_phi",100, phimin, phimax);
00322     h_jet_reco_presel_nb = new TH1F("jet_reco_presel_nb","jet_reco_presel_nb",19, -1, 18);
00323 
00324     h_jet_truth_presel_pt = new TH1F("jet_truth_presel_pt","jet_truth_presel_pt",100, ptmin, ptmax);
00325     h_jet_truth_presel_eta = new TH1F("jet_truth_presel_eta","jet_truth_presel_eta",100, etamin, etamax);
00326     h_jet_truth_presel_phi = new TH1F("jet_truth_presel_phi","jet_truth_presel_phi",100, phimin, phimax);
00327     h_jet_truth_presel_nb = new TH1F("jet_truth_presel_nb","jet_truth_presel_nb",19, -1, 18);
00328 
00329     h_jet_reco_matched_pt = new TH1F("jet_reco_matched_pt","jet_reco_matched_pt",100, ptmin, ptmax);
00330     h_jet_reco_matched_eta = new TH1F("jet_reco_matched_eta","jet_reco_matched_eta",100, etamin, etamax);
00331     h_jet_reco_matched_phi = new TH1F("jet_reco_matched_phi","jet_reco_matched_phi",100, phimin, phimax);
00332     h_jet_reco_matched_nb = new TH1F("jet_reco_matched_nb","jet_reco_matched_nb",19, -1, 18);
00333 
00334     h_jet_truth_matched_pt = new TH1F("jet_truth_matched_pt","jet_truth_matched_pt",100, ptmin, ptmax);
00335     h_jet_truth_matched_eta = new TH1F("jet_truth_matched_eta","jet_truth_matched_eta",100, etamin, etamax);
00336     h_jet_truth_matched_phi = new TH1F("jet_truth_matched_phi","jet_truth_matched_phi",100, phimin, phimax);
00337 
00338     h_jet_reco_fake_pt = new TH1F("jet_reco_fake_pt","jet_reco_fake_pt",100, ptmin, ptmax);
00339     h_jet_reco_fake_eta = new TH1F("jet_reco_fake_eta","jet_reco_fake_eta",100, etamin, etamax);
00340     h_jet_reco_fake_phi = new TH1F("jet_reco_fake_phi","jet_reco_fake_phi",100, phimin, phimax);
00341     h_jet_reco_fake_nb = new TH1F("jet_reco_fake_nb","jet_reco_fake_nb",19, -1, 18);
00342 
00343     h_jet_truth_notmatched_pt = new TH1F("jet_truth_notmatched_pt","jet_truth_notmatched_pt", 100, ptmin, ptmax);
00344     h_jet_truth_notmatched_eta = new TH1F("jet_truth_notmatched_eta","jet_truth_notmatched_eta", 100, etamin, etamax);
00345     h_jet_truth_notmatched_phi = new TH1F("jet_truth_notmatched_phi","jet_truth_notmatched_phi", 100, phimin, phimax);
00346     h_jet_truth_notmatched_nb = new TH1F("jet_truth_notmatched_nb","jet_truth_notmatched_nb",19, -1, 18);
00347 
00348     h_jet_matched_dpt = new TH1F("jet_matched_dpt","jet_matched_dpt", 200, -d_ptmax, d_ptmax);
00349     h_jet_matched_deta = new TH1F("jet_matched_deta","jet_matched_deta", 200, -d_etamax, d_etamax);
00350     h_jet_matched_dphi = new TH1F("jet_matched_dphi","jet_matched_dphi", 200, -d_phimax, d_phimax);
00351     h_jet_matched_dr = new TH1F("jet_matched_dr","jet_matched_dr", 200, 0.0, 0.1);
00352 
00353     h_jet_matched_resolution_pt = new TH2F("jet_matched_resolution_pt","jet_matched_resolution_pt",
00354                      20, ptmin, ptmax, 100, -1.0, 1.0);
00355     h_jet_matched_resolution_eta = new TH2F("jet_matched_resolution_eta","jet_matched_resolution_eta",
00356                      20, etamin, etamax, 100, -1.0, 1.0);
00357     h_jet_matched_resolution_phi = new TH2F("jet_matched_resolution_phi","jet_matched_resolution_phi",
00358                      20, phimin, phimax, 100, -1.0, 1.0);
00359         h_remove_double_dpt = new TH1F("remove_double_dpt","remove_double_dpt", 200, -d_ptmax, d_ptmax);
00360     h_remove_double_dr = new TH1F("remove_double_dr","remove_double_dr", 200, 0.0, 0.5);
00361     h_remove_double_nb_before = new TH1F("remove_double_nb_before","remove_double_nb_before",19, -1, 18);
00362     h_remove_double_nb_after = new TH1F("remove_double_nb_after","remove_double_nb_after",19, -1, 18);
00363     h_remove_double_nb_doubles = new TH1F("remove_double_nb_doubles","remove_double_nb_doubles",19, -1, 18);
00364 
00365     dir_presel_matching->cd("../");
00366 
00367 
00368 
00369 // MET histos
00370     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00371         dir_MET = new TDirectoryFile("MET","MET");
00372     #else
00373         dir_MET = new TDirectory("MET","MET");
00374     #endif
00375     //dir_MET = new TDirectory("MET","MET");
00376     dir_MET->cd();
00377     // recalculation of missing et (muon part)
00378     h_etmiss_muons_reco_truth_dx = new TH1F("etmiss_muons_reco_truth_dx","etmiss_muons_reco_truth_dx",100, -0.2*ptmax, 0.2*ptmax);
00379     h_etmiss_muons_reco_truth_dy = new TH1F("etmiss_muons_reco_truth_dy","etmiss_muons_reco_truth_dy",100, -0.2*ptmax, 0.2*ptmax);
00380     h_etmiss_muons_reco_truth_dpt = new TH1F("etmiss_muons_reco_truth_dpt","etmiss_muons_reco_truth_dpt",100, -0.2*ptmax, 0.2*ptmax);
00381     h_etmiss_muons_reco_x = new TH1F("etmiss_muons_reco_x","etmiss_muons_reco_x",100, -0.2*ptmax, 0.2*ptmax);
00382     h_etmiss_muons_reco_y = new TH1F("etmiss_muons_reco_y","etmiss_muons_reco_y",100, -0.2*ptmax, 0.2*ptmax);
00383     h_etmiss_muons_reco_pt = new TH1F("etmiss_muons_reco_pt","etmiss_muons_reco_pt",100, -0.2*ptmax, 0.2*ptmax);
00384     h_etmiss_muons_truth_x = new TH1F("etmiss_muons_truth_x","etmiss_muons_truth_x",100, -0.2*ptmax, 0.2*ptmax);
00385     h_etmiss_muons_truth_y = new TH1F("etmiss_muons_truth_y","etmiss_muons_truth_y",100, -0.2*ptmax, 0.2*ptmax);
00386     h_etmiss_muons_truth_pt = new TH1F("etmiss_muons_truth_pt","etmiss_muons_truth_pt",100, -0.2*ptmax, 0.2*ptmax);
00387 
00388     // complete MET:
00389     h_etmiss_reco_truth_dx = new TH1F("etmiss_reco_truth_dx","etmiss_reco_truth_dx",100, -0.2*ptmax, 0.2*ptmax);
00390     h_etmiss_reco_truth_dy = new TH1F("etmiss_reco_truth_dy","etmiss_reco_truth_dy",100, -0.2*ptmax, 0.2*ptmax);
00391     h_etmiss_reco_truth_dpt = new TH1F("etmiss_reco_truth_dpt","etmiss_reco_truth_dpt",100, -0.2*ptmax, 0.2*ptmax);
00392     h_etmiss_reco_x = new TH1F("etmiss_reco_x","etmiss_reco_x",100, -0.2*ptmax, 0.2*ptmax);
00393     h_etmiss_reco_y = new TH1F("etmiss_reco_y","etmiss_reco_y",100, -0.2*ptmax, 0.2*ptmax);
00394     h_etmiss_reco_pt = new TH1F("etmiss_reco_pt","etmiss_reco_pt",100, -0.2*ptmax, 0.2*ptmax);
00395     h_etmiss_truth_x = new TH1F("etmiss_truth_x","etmiss_truth_x",100, -0.2*ptmax, 0.2*ptmax);
00396     h_etmiss_truth_y = new TH1F("etmiss_truth_y","etmiss_truth_y",100, -0.2*ptmax, 0.2*ptmax);
00397     h_etmiss_truth_pt = new TH1F("etmiss_truth_pt","etmiss_truth_pt",100, -0.2*ptmax, 0.2*ptmax);
00398 
00399        h_etmiss_reco_truth_dptvscut = new TH2F("etmiss_reco_truth_dptvscut","etmiss_reco_truth_dptvscut",
00400                      22,-2.,20., 100, -0.2*ptmax, 0.2*ptmax);
00401     h_etmiss_reco_truth_dptvsmettruth = new TH2F("etmiss_reco_truth_dptvsmettruth","etmiss_reco_truth_dptvsmettruth",
00402                      100, ptmin, ptmax, 100, -0.2*ptmax, 0.2*ptmax);
00403     h_etmiss_reco_truth_dpt_vs_muondpt = new TH2F("etmiss_reco_truth_dpt_vs_muondpt",
00404         "MET reco(muon staco)-truth(muon) vs. reco(met_in_ana)-truth(nonint)", 100,-1*ptmax,ptmax, 100, -1*ptmax, 1*ptmax);
00405 
00406 
00407     h_etmiss_invest_reco_truth_dpx = new TH1F("etmiss_invest_reco_truth_dpx","etmiss_invest_reco_truth_dpx",100, -100, 100);
00408     h_etmiss_invest_reco_truth_dpy = new TH1F("etmiss_invest_reco_truth_dpy","etmiss_invest_reco_truth_dpy",100, -100, 100);
00409     h_etmiss_invest_reco_truth_dpt = new TH1F("etmiss_invest_reco_truth_dpt","etmiss_invest_reco_truth_dpt",100, 0, 200);
00410     h_etmiss_invest_dpt = new TH2F("etmiss_invest_dpt","etmiss_invest_dpt", 200, -200, 200, 200, -200, 200);
00411     h_etmiss_MB_dpt = new TH1F("etmiss_MB_dpt","etmiss_MB_dpt", 1000, -200, 200);
00412     h_etmiss_staco_dpt = new TH1F("etmiss_staco_dpt","etmiss_staco_dpt", 1000, -200, 200);
00413 
00414     dir_MET->cd("../");
00415 
00416 
00417 // histos for the analysis
00418     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00419         dir_analysis = new TDirectoryFile("analysis","analysis");
00420     #else
00421         dir_analysis = new TDirectory("analysis","analysis");
00422     #endif
00423     //dir_analysis = new TDirectory("analysis","analysis");
00424     dir_analysis->cd();
00425     // central jets
00426     h_centraljet_ana_nb = new TH1F("centraljet_ana_nb","centraljet_ana_nb",19, -1, 18);
00427     h_centraljet_ana_pt = new TH1F("centraljet_ana_pt","centraljet_ana_pt", 100, ptmin, ptmax);
00428     h_centraljet_ana_eta = new TH1F("centraljet_ana_eta","centraljet_ana_eta", 100, etamin, etamax);
00429     h_centraljet_ana_phi = new TH1F("centraljet_ana_phi","centraljet_ana_phi", 100, phimin, phimax);
00430     h_centraljet_ana_zeppvar = new TH1F("centraljet_ana_zeppvar","centraljet_ana_zeppvar", 100, etamin, etamax);
00431 
00432 
00433     h_fwj_cj_dr = new TH1F("fwj_cj_dr","Delta R between central jet and forward jets", 200, 0.0, 6.);
00434 
00435        // cut evolution
00436     h_cut_evolution = new TH1D("cut_evolution","cut_evolution",27,-2.,25.);
00437     h_cut_evolution_candidates = new TH1D("cut_evolution_candidates","cut_evolution_candidates",27,-2.,25.);
00438 
00439     h_trigger_online = new TH1F("trigger_online","trigger_online",19, -1, 18);
00440     h_trigger_offline = new TH1F("trigger_offline","trigger_offline",19, -1, 18);
00441     h_trigger_combined = new TH1F("trigger_combined","trigger_combined",19, -1, 18);
00442 
00443 
00444 
00445    for (int j=0; j<26; j++){
00446     string nb;
00447      sprintf(char_help_array,"%02d",j);
00448     nb = char_help_array;
00449     string s1  = "";
00450 
00451     s1 = "etmiss_ana_reco_pt_"+nb;
00452     v_h_etmiss_ana_reco_pt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00453     s1 = "etmiss_ana_reco_x_"+nb;
00454     v_h_etmiss_ana_reco_x.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00455     s1 = "etmiss_ana_reco_y_"+nb;
00456     v_h_etmiss_ana_reco_y.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00457 
00458     s1 = "etmiss_ana_truth_pt_"+nb;
00459     v_h_etmiss_ana_truth_pt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00460     s1 = "etmiss_ana_truth_x_"+nb;
00461     v_h_etmiss_ana_truth_x.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00462     s1 = "etmiss_ana_truth_y_"+nb;
00463     v_h_etmiss_ana_truth_y.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00464 
00465     s1 = "etmiss_ana_reco_truth_dpt_"+nb;
00466     v_h_etmiss_ana_reco_truth_dpt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 0., 0.5*ptmax));
00467     s1 = "etmiss_ana_reco_truth_dpt_scalar_"+nb;
00468     v_h_etmiss_ana_reco_truth_dpt_scalar.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.5*ptmax, 0.5*ptmax));
00469     s1 = "etmiss_ana_reco_truth_dx_"+nb;
00470     v_h_etmiss_ana_reco_truth_dx.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.3*ptmax, 0.3*ptmax));
00471     s1 = "etmiss_ana_reco_truth_dy_"+nb;
00472     v_h_etmiss_ana_reco_truth_dy.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.3*ptmax, 0.3*ptmax));
00473     
00474     s1 = "etmiss_muon_ana_reco_pt_"+nb;
00475     v_h_etmiss_muon_ana_reco_pt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00476     s1 = "etmiss_muon_ana_reco_x_"+nb;
00477     v_h_etmiss_muon_ana_reco_x.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00478     s1 = "etmiss_muon_ana_reco_y_"+nb;
00479     v_h_etmiss_muon_ana_reco_y.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00480     
00481     s1 = "etmiss_muon_ana_truth_pt_"+nb;
00482     v_h_etmiss_muon_ana_truth_pt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00483     s1 = "etmiss_muon_ana_truth_x_"+nb;
00484     v_h_etmiss_muon_ana_truth_x.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00485     s1 = "etmiss_muon_ana_truth_y_"+nb;
00486     v_h_etmiss_muon_ana_truth_y.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 1.5*ptmin, 1.5*ptmax));
00487     
00488     s1 = "etmiss_muon_ana_reco_truth_dpt_"+nb;
00489     v_h_etmiss_muon_ana_reco_truth_dpt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, 0., 0.5*ptmax));
00490     s1 = "etmiss_muon_ana_reco_truth_dpt_scalar_"+nb;
00491     v_h_etmiss_muon_ana_reco_truth_dpt_scalar.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.5*ptmax, 0.5*ptmax));
00492     s1 = "etmiss_muon_ana_reco_truth_dx_"+nb;
00493     v_h_etmiss_muon_ana_reco_truth_dx.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.3*ptmax, 0.3*ptmax));
00494     s1 = "etmiss_muon_ana_reco_truth_dy_"+nb;
00495     v_h_etmiss_muon_ana_reco_truth_dy.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.3*ptmax, 0.3*ptmax));
00496     
00497     s1 = "e_ana_nb_"+nb;
00498        v_h_e_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00499       
00500     s1 = "e_ana_pt_"+nb;
00501        v_h_e_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, ptmin, ptmax));
00502       
00503     s1 = "e_ana_eta_"+nb;
00504        v_h_e_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, etamin, etamax));
00505       
00506     s1 = "e_ana_phi_"+nb;
00507        v_h_e_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, phimin, phimax));
00508       
00509     s1 = "mu_ana_nb_"+nb;
00510        v_h_mu_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00511       
00512     s1 = "mu_ana_pt_"+nb;
00513        v_h_mu_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, ptmin, ptmax));
00514       
00515     s1 = "mu_ana_eta_"+nb;
00516        v_h_mu_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, etamin, etamax));
00517       
00518     s1 = "mu_ana_phi_"+nb;
00519        v_h_mu_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, phimin, phimax));
00520       
00521     s1 = "taujet_ana_nb_"+nb;
00522        v_h_taujet_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00523       
00524     s1 = "taujet_ana_pt_"+nb;
00525        v_h_taujet_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, ptmin, ptmax));
00526       
00527     s1 = "taujet_ana_eta_"+nb;
00528        v_h_taujet_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, etamin, etamax));
00529       
00530     s1 = "taujet_ana_phi_"+nb;
00531        v_h_taujet_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, phimin, phimax));
00532       
00533     s1 = "lepton_ana_nb_"+nb;
00534        v_h_lepton_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00535       
00536     s1 = "lepton_ana_pt_"+nb;
00537        v_h_lepton_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, ptmin, ptmax));
00538       
00539     s1 = "lepton_ana_eta_"+nb;
00540        v_h_lepton_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, etamin, etamax));
00541       
00542     s1 = "lepton_ana_phi_"+nb;
00543        v_h_lepton_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, phimin, phimax));
00544       
00545       s1 = "2l_ana_nb_"+nb;
00546     v_h_2l_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str(),11, -1, 10));
00547 
00548     s1 = "2l_ana_deta_"+nb;
00549     v_h_2l_ana_deta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -0.2, 6.));
00550     
00551     s1 = "2l_ana_dphi_"+nb;
00552     v_h_2l_ana_dphi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -0.1, 4.));
00553     
00554     s1 = "2l_ana_dpt_"+nb;
00555     v_h_2l_ana_dpt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -10., 200.));
00556     
00557     s1 = "2l_ana_dr_"+nb;
00558     v_h_2l_ana_dr.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 0., 6.));
00559     
00560     s1 = "2l_ana_m_"+nb;
00561     v_h_2l_ana_m.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 0., 400.));
00562     
00563     s1 = "jet_ana_nb_"+nb;
00564        v_h_jet_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00565       
00566     s1 = "jet_ana_pt_"+nb;
00567        v_h_jet_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 1.5*ptmin, 1.5*ptmax));
00568       
00569     s1 = "jet_ana_eta_"+nb;
00570        v_h_jet_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, etamin, etamax));
00571       
00572     s1 = "jet_ana_phi_"+nb;
00573        v_h_jet_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, phimin, phimax));
00574     
00575       s1 = "2j_ana_nb_"+nb;
00576        v_h_2j_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,11, -1, 10));
00577     
00578     s1 = "2j_ana_deta_"+nb;
00579     v_h_2j_ana_deta.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -0.2, 12.));
00580     
00581     s1 = "2j_ana_dphi_"+nb;
00582     v_h_2j_ana_dphi.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -0.1, 4.));
00583     
00584     s1 = "2j_ana_dpt_"+nb;
00585     v_h_2j_ana_dpt.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, -10., 300.));
00586     
00587     s1 = "2j_ana_dr_"+nb;
00588     v_h_2j_ana_dr.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 0.0, 12.));
00589     
00590     s1 = "2j_ana_m_"+nb;
00591     v_h_2j_ana_m.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 0., 5000.));
00592     
00593     s1 = "2l_extau1_"+nb;
00594     v_h_2l_extau1.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
00595     
00596     s1 = "2l_extau2_"+nb;
00597     v_h_2l_extau2.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
00598     
00599     s1 = "2l_extau12_"+nb;
00600     v_h_2l_extau12.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
00601     
00602     s1 = "2l_extau1vs2_"+nb;
00603     v_h_2l_extau1vs2.push_back(new TH2F(s1.c_str(),s1.c_str() ,100,-3.,3., 100,-3.,3.));
00604     
00605     s1 = "higgs_mass_"+nb;
00606     v_h_higgs_mass.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,0.,300.));
00607 
00608     s1 = "lj_deta_min_"+nb;
00609     v_h_lj_deta_min.push_back(new TH1F(s1.c_str(),s1.c_str() ,50,-5.,5.));
00610     
00611     s1 = "lj_deta_max_"+nb;
00612     v_h_lj_deta_max.push_back(new TH1F(s1.c_str(),s1.c_str() ,50,-5.,5.));
00613 
00614    }
00615    
00616    
00617 //    h_2l_pz1_9       = new TH1F("h_2l_pz1_9"    ,"h_2l_pz1_9"        , 500 , -500000. , 500000.);
00618 //    h_2l_pztotal_9  = new TH1F("h_2l_pztotal_9"  ,"h_2l_pztotal_9"   , 500 , -1000000. , 1000000.);
00619 //    h_2l_phi2_9     = new TH1F("h_2l_phi2_9"    ,"h_2l_phi2_9"      , 60 , -4. ,4.);
00620 //    h_2l_theta_9    = new TH1F("h_2l_theta_9"    ,"h_2l_theta_9"     , 60 , 0. ,4.);
00621 //    h_2l_projmag_9  = new TH1F("h_2l_projmag_9"  ,"h_2l_projmag_9"   , 100 , 0., 10.);
00622 //    h_2l_projtheta_9= new TH1F("h_2l_projtheta_9","h_2l_projtheta_9" , 600 , -4. ,4.);
00623 //    h_2l_projphi_9  = new TH1F("h_2l_projphi_9"  ,"h_2l_projphi_9"   , 60 , -4. ,4.);
00624 //    h_2l_projdist_9 = new TH1F("h_2l_projdist_9" ,"h_2l_projdist_9"  , 60 , 0. ,15.);
00625 //    h_2l_projdist2_9= new TH1F("h_2l_projdist2_9","h_2l_projdist2_9" , 60 , 0. ,15.);
00626 // // CUT 10 
00627 //    h_2l_norig     =new TH1F("h_2l_norig"   ,"h_2l_norig",20,0,20);
00628 //    h_2l_nbefore  =new TH1F("h_2l_nbefore" ,"h_2l_nbefore",20,0,20);
00629 //    h_2l_nbefore0 =new TH1F("h_2l_nbefore0","h_2l_nbefore0",20,0,20);
00630 //    h_2l_nmiddle1 =new TH1F("h_2l_nmiddle1","h_2l_nmiddle1",20,0,20);
00631 //    h_2l_nmiddle2 =new TH1F("h_2l_nmiddle2","h_2l_nmiddle2",20,0,20);
00632 //    h_2l_nafter   =new TH1F("h_2l_nafter"  ,"h_2l_nafter",20,0,20);
00633 //    h_2l_correct  =new TH1F("h_2l_correct" ,"h_2l_correct",20,0,20);
00634    
00635 // // CUT 11: all jet pairs 
00636 //    h_2j_leadingpt    = new TH1F("h_2j_leadingpt"    ,"h_2j_leadingpt"      ,200,0.,200000.);
00637 //    h_2j_subleadingpt = new TH1F("h_2j_subleadingpt" ,"h_2j_subleadingpt"  ,200,0.,200000.);
00638 // // CUT 12: jet pt cut
00639 //    h_2j_phi2_12        = new TH1F("h_2j_phi2_12"      ,"h_2j_phi2_12"      , 60 , -4. ,4.);
00640 //    h_2j_theta_12       = new TH1F("h_2j_theta_12"     ,"h_2j_theta_12"    , 60 , 0. ,4.);
00641 //    h_2j_projmag_12     = new TH1F("h_2j_projmag_12"   ,"h_2j_projmag_12"  , 100 , 0., 10.);
00642 //    h_2j_projtheta_12   = new TH1F("h_2j_projtheta_12" ,"h_2j_projtheta_12", 600 , -4. ,4.);
00643 //    h_2j_projphi_12     = new TH1F("h_2j_projphi_12"   ,"h_2j_projphi_12"  , 60 , -4. ,4.);
00644 //    h_2j_projdist_12    = new TH1F("h_2j_projdist_12"  ,"h_2j_projdist_12" , 60 , 0. ,15.);
00645 //    h_2j_projdist2_12   = new TH1F("h_2j_projdist2_12" ,"h_2j_projdist2_12", 60 , 0. ,15.);
00646 // // CUT 13: jet eta cut
00647 // // CUT 14: forward-backward jets
00648 // // CUT 15: jet eta separation 
00649 // // CUT 16: invariant mass of two jets
00650 //   h_2j_ncentral        = new TH1F("2j_ncentral"        ,"2j_ncentral"         , 15, 0 ,15);
00651 // // CUT 17: central jet veto
00652 // // CUT 18: delta phi between the jets
00653 //    h_etmiss_18         = new TH1F("h_etmiss_18"      ,"h_etmiss_18"       ,300,0.,300000.);
00654 //    h_2j_phi2_18        = new TH1F("h_2j_phi2_18"      ,"h_2j_phi2_18"      , 60 , -4. ,4.);
00655 //    h_2j_theta_18       = new TH1F("h_2j_theta_18"     ,"h_2j_theta_18"    , 60 , 0. ,4.);
00656 //    h_2j_projmag_18     = new TH1F("h_2j_projmag_18"   ,"h_2j_projmag_18"  , 100 , 0., 10.);
00657 //    h_2j_projtheta_18   = new TH1F("h_2j_projtheta_18" ,"h_2j_projtheta_18", 600 , -4. ,4.);
00658 //    h_2j_projphi_18     = new TH1F("h_2j_projphi_18"   ,"h_2j_projphi_18"  , 60 , -4. ,4.);
00659 //    h_2j_projdist_18    = new TH1F("h_2j_projdist_18"  ,"h_2j_projdist_18" , 60 , 0. ,15.);
00660 //    h_2j_projdist2_18   = new TH1F("h_2j_projdist2_18" ,"h_2j_projdist2_18", 60 , 0. ,15.);
00661 // CUT 19 : at least one correct jet pair, at least one correct lepton pair
00662 //    h_2j_norig     = new TH1F("2j_norig"   ,"2j_norig",20,0,20);
00663 //    h_2j_nbefore  = new TH1F("2j_nbefore" ,"2j_nbefore",20,0,20);
00664 //    h_2j_nmiddle1 = new TH1F("2j_nmiddle1","2j_nmiddle1",20,0,20);
00665 //    h_2j_nmiddle2 = new TH1F("2j_nmiddle2","2j_nmiddle2",20,0,20);
00666 //    h_2j_nmiddle3 = new TH1F("2j_nmiddle3","2j_nmiddle3",20,0,20);
00667 //    h_2j_nmiddle4 = new TH1F("2j_nmiddle4","2j_nmiddle4",20,0,20);
00668 //    h_2j_nmiddle5 = new TH1F("2j_nmiddle5","2j_nmiddle5",20,0,20);
00669 //    h_2j_nafter   = new TH1F("2j_nafter"  ,"2j_nafter",20,0,20);
00670  
00671 //    h_2j_pt_f           = new TH1F("2j_pt_f"          ,"2j_pt_f"        ,300,-100000.,200000.);
00672 //    h_2j_eta_f           = new TH1F("2j_eta_f"          ,"2j_eta_f"        ,120,-6.,6.);
00673 //    h_2j_phi_f           = new TH1F("2j_phi_f"          ,"2j_phi_f"        ,40,-4.,4.);
00674 //    h_2j_dphi_f         = new TH1F("2j_dphi_f"          ,"2j_dphi_f"        ,40,-4.,4.);
00675 //    h_2j_deta_f         = new TH1F("2j_deta_f"          ,"2j_deta_f"        ,100,-12.,12.);
00676 //    h_2j_m_f           = new TH1F("2j_m_f"          ,"2j_m_f"        ,100,0.,5000000.);
00677 //    h_etmiss_f           = new TH1F("etmiss_f"          ,"etmiss_f"       ,300,0.,300000.);
00678 //    h_2j_ncentral_f     = new TH1F("2j_ncentral_f"     ,"2j_ncentral_f"          , 15, 0 ,15);
00679 
00680     h_lmissinget_ana_mt = new TH1F("lmissinget_ana_mt","lmissinget_ana_mt",100, ptmin, ptmax);
00681 
00682 
00683 //   h_higgs_mass      = new TH1F("higgs_mass"   ,"higgs_mass"   ,100,0.,300.);
00684    h_higgs_pt         = new TH1F("higgs_pt"     ,"higgs_pt"     ,100,0.,500.);
00685    h_higgs_eta       = new TH1F("higgs_eta"    ,"higgs_eta"    ,120,-6.,6.);
00686    h_higgs_phi       = new TH1F("higgs_phi"    ,"higgs_phi"    ,40,-4.,4.);
00687    h_2l_phi2_final     = new TH1F("2l_phi2_final"      ,"2l_phi2_final"    , 60 , -4. ,4.);
00688    h_2l_theta_final    = new TH1F("2l_theta_final"      ,"2l_theta_final"   ,  60 , 0. ,4.);
00689    h_2l_projmag_final  = new TH1F("2l_projmag_final"      ,"2l_projmag_final" , 100 , 0., 10.);
00690    h_2l_projtheta_final= new TH1F("2l_projtheta_final"    ,"2l_projtheta_final", 600 , -4. ,4.);
00691    h_2l_projphi_final  = new TH1F("2l_projphi_final"      ,"2l_projphi_final"  , 60 , -4. ,4.);
00692    h_2l_projdist_final = new TH1F("2l_projdist_final"      ,"2l_projdist_final" , 60 , 0. ,15.);
00693    h_2l_projdist2_final= new TH1F("2l_projdist2_final"    ,"2l_projdist2_final", 60 , 0. ,15.);
00694 
00695 //    // truth info, VBF cut
00696 //    h_DecayMode    = new TH1F("DecayMode"    ,"DecayMode"     ,10, 0 , 10. );
00697 //    h_DecayMode_VBF= new TH1F("DecayMode_VBF","DecayMode_VBF" ,10, 0 , 10. );
00698 //    h_DecayMode_f  = new TH1F("DecayMode_f  ","DecayMode_f"   ,10, 0 , 10. );
00699 //    h_NTau      = new TH1F("NTau"        ,"NTau"         ,10 , 0 , 10  );
00700 //    h_NTau_VBF      = new TH1F("NTau_VBF"     ,"NTau_VBF"      ,10 , 0 , 10  );
00701 //    h_MInv      = new TH1F("MInv"        ,"MInv"         ,200 ,0 , 200000  );
00702 //    h_MInv_VBF      = new TH1F("MInv_VBF"     ,"MInv_VBF"      ,200 ,0 , 200000  );
00703 //    h_PtTauVis      = new TH1F("PtTauVis"     ,"PtTauVis"      ,200 ,0 , 200000  );
00704 //    h_PtTauVis_VBF = new TH1F("PtTauVis_VBF" ,"PtTauVis_VBF"  ,200 ,0 , 200000  );
00705 //    h_NJetOrig      = new TH1F("NJetOrig"     ,"NJetOrig"      ,30 , 0 , 30  );
00706 //    h_NJet      = new TH1F("NJet"        ,"NJet"         ,30 , 0 , 30  );
00707     dir_analysis->cd("../");
00708 
00709     
00710     
00711 // Folders and other Info ntuples
00712     #if ROOT_VERSION_CODE >= ROOT_VERSION(5,15,0)
00713         dir_cut_info = new TDirectoryFile("cut_info","cut_info");
00714     #else
00715         dir_cut_info = new TDirectory("cut_info","cut_info");
00716     #endif
00717     //dir_cut_info = new TDirectory("cut_info","cut_info");
00718     dir_cut_info->cd();
00719     tree_cut_order = new TTree("cut_order","cut_order");
00720     tree_cut_values = new TTree("cut_values","cut_values");
00721     dir_cut_info->cd("../");
00722 
00723     tree_settings = new TTree("settings","settings");
00724     tree_MVA_variables = new TTree("MVA_variables","MVA_variables");
00725     tree_cand_data = new TTree("cand_data","cand_data");
00726 
00727