MyHtoWWAnalysis/src/HtoWWHistos.cxx

// for normal histos
double ptmin(0.0);
double ptmax(250.0);
double etamin(-5.1);
double etamax(5.1);
double etamin_central(-2.5);
double etamax_central(2.5);
double phimin(-5.0);
double phimax(5.0);

int ptbins(60);  //50
int etabins(60); //50
int etabins_central(40); //50
int phibins(60); //60
int Ebins(60);   //60

double Emin(0.0);
double Emax(1000.5);
// for delta histos, e.g. delta eta of two particles
double d_ptmax(50.0);
double d_etamax(0.2);
double d_phimax(0.2);   
double d_Emax(50.0);


//cut opt
t_cutopt_std_jet = new TNtuple("cutopt_std_jet_ntuple", "cutopt_std_jet_ntuple",
                               "weight:ncjets:pt:eta:phi:ntrks:ntrk_primvert:ptratio:cmsptratio");
t_cutopt_trk_jet = new TNtuple("cutopt_trk_jet_ntuple", "cutopt_trk_jet_ntuple",
                               "weight:ncjets:pt:eta:phi");

t_cutopt_e_isol     = new TNtuple("cutopt_e_isol_ntuple", "cutopt_e_isol_ntuple", 
                                  "weight:etcone20:etcone30:etcone40:etcone45");
t_cutopt_e_nonisol  = new TNtuple("cutopt_e_nonisol_ntuple", "cutopt_e_nonisol_ntuple", 
                                  "weight:etcone20:etcone30:etcone40:etcone45");

t_cutopt_mu_isol    = new TNtuple("cutopt_mu_isol_ntuple", "cutopt_mu_isol_ntuple", 
                                  "weight:etcone10:etcone20:etcone30:etcone40:etcone45:etcone50:etcone60:etcone70");
t_cutopt_mu_nonisol = new TNtuple("cutopt_mu_nonisol_ntuple", "cutopt_mu_nonisol_ntuple", 
                                  "weight:etcone10:etcone20:etcone30:etcone40:etcone45:etcone50:etcone60:etcone70");

// vertex ntuple //
t_vertex_ntuple = new TNtuple("vertex_ntuple", "Vertex Ntuple",
                              "weight:rec_nb_vertices:chi2:rec_nb_tracks:rec_Et:rec_x:rec_y:rec_z:truth_nb_vertices:truth_nb_tracks:truth_Et:truth_x:truth_y:truth_z");

//dijet histos ->change
h_jets_dr = new TH1F("jets_dr","jets_dr",200,0.,6.);
h_jets_dr_cut = new TH1F("jets_dr_cut","jets_dr_cut",200,0.,6.);

// vertex histograms //
h_nb_vertices = new TH1F("nb_vertices",
                        "Number of Primary Vertices",
                        11, -0.5, 10.5);
        h_nb_vertices->SetXTitle("number of vertices");
h_Et_selected_vertex = new TH1F("Et_selected_vertex", 
                        "Total Transverse Track Energy for the Selected Vertex",
                        100, -0.5, 400.5);
        h_Et_selected_vertex->SetXTitle("E_T (GeV)");
        h_Et_selected_vertex->SetYTitle("dn/dE_T (1/(0.25 GeV))");
h_Et_v2_vs_Et_v1 = new TH2F("Et_v2_vs_Et_v1",
                        "Et(vertex 2) versus Et(vertex 1)",
                        401, -0.5, 400.5, 401, -0.5, 400.5);
        h_Et_v2_vs_Et_v1->SetXTitle("Et(vertex 1) (GeV)");
        h_Et_v2_vs_Et_v1->SetYTitle("Et(vertex 2) (GeV)");

h_delta_vertices = new TH1F("delta_vertices",
                        "Distance of Vertices",
                        210, -0.105, 0.105);
        h_delta_vertices->SetXTitle("#Delta(vertex) (mm)");
        h_delta_vertices->SetYTitle("dn/d#delta(vertex) (1/1 #mu m)");

h_rec_nb_tracks_from_vertex = new TH1F("rec_nb_tracks_from_vertex",
                        "Number of Tracks Emerging from the Reconstructed Vertex",
                        151, -0.5, 150.5);
        h_rec_nb_tracks_from_vertex->SetXTitle("number of tracks");

h_nb_truth_vertices = new TH1F("nb_truth_vertices",
                        "Number of Vertices in Truth",
                        11, -0.5, 10.5);
        h_nb_truth_vertices->SetXTitle("number of vertices");
h_Et_selected_truth_vertex = new TH1F("Et_selected_truth_vertex", 
                        "Total Transverse Track Energy for the Selected Vertex in Truth",
                        100, -0.5, 400.5);
        h_Et_selected_truth_vertex->SetXTitle("E_T (GeV)");
        h_Et_selected_truth_vertex->SetYTitle("dn/dE_T (1/(0.25 GeV))");
h_truth_Et_v2_vs_Et_v1 = new TH2F("truth_Et_v2_vs_Et_v1",
                        "Et(vertex 2) versus Et(vertex 1) in Truth",
                        401, -0.5, 400.5, 401, -0.5, 400.5);
        h_truth_Et_v2_vs_Et_v1->SetXTitle("Et(vertex 1) (GeV)");
        h_truth_Et_v2_vs_Et_v1->SetYTitle("Et(vertex 2) (GeV)");

h_delta_selected_vertices = new TH1F("delta_selected_vertices",
                        "Distrance of Reconstructed Vertex from Truth Vertex",
                        110, -0.055, 0.055);
        h_delta_selected_vertices->SetXTitle("#Delta z(vertex) (mm)");
        h_delta_selected_vertices->SetYTitle("dn/d#Delta z(vertex) (1/1 #mu m)");

h_truth_nb_tracks_from_vertex = new TH1F("truth_nb_tracks_from_vertex",
                        "Number of Tracks Emerging from the Truth Vertex",
                        151, -0.5, 150.5);
        h_truth_nb_tracks_from_vertex->SetXTitle("number of tracks");

h_vertex_ratio_number_of_tracks = new TH1F("vertex_ratio_number_of_tracks",
                        "Number of Tracks from Reconstructed Vertex/Number of Tracks from Truth Vertex",
                        87, -0.05, 2.55);
        h_vertex_ratio_number_of_tracks->SetXTitle("ratio");
        h_vertex_ratio_number_of_tracks->SetYTitle("dn/dratio (1/0.03)");

h_vertex_ratio_Et = new TH1F("vertex_ratio_Et",
                        "Et(reconstructed vertex)/Et(truth vertex)",
                        87, -0.05, 2.55);
        h_vertex_ratio_Et->SetXTitle("ratio");
        h_vertex_ratio_Et->SetYTitle("dn/dratio (1/0.03)");

// track-jet ntuples //
t_trj_ntuple = new TNtuple("trj_ntuple", "Track-Jet Ntuple",
                "event_weight:rec_nb_tracks:truth_nb_tracks:nb_trj_jets:nb_truth_trj_jets:Et:eta:phi:dR:truth_Et:truth_eta:truth_phi");
t_trj_ntuple_2 = new TNtuple("trj_ntuple_2", "Track-Jet Ntuple 2",
                "event_weight:min_Et:nb_jets");
t_std_ntuple_2 = new TNtuple("std_ntuple_2", "Standard-Jet Ntuple 2",
                "event_weight:min_Et:nb_jets");
t_std_ntuple_3 = new TNtuple("std_ntuple_3", "Standard-Jet Ntuple 3",
                "event_weight:min_Et:nb_jets");
t_truth_trj_ntuple = new TNtuple("truth_trj_ntuple", "Truth Track-Jet Ntuple",
                "event_weight:rec_nb_tracks:truth_nb_tracks:nb_trj_jets:nb_truth_trj_jets:Et:eta:phi:dR:truth_Et:truth_eta:truth_phi");

// track-jet histograms //
// h_trj_nor_nb_jets = new TH1F("trj_nor_nb_jets",
//              "Number of Track Jets (no overlap removal)",
//              16, -0.5, 15.5);
//      h_trj_nor_nb_jets->SetXTitle("number of jets");
// h_trj_nor_eta = new TH1F("trj_nor_eta",
//              "Pseudorapidity Distribution of Track Jets (no overlap removal)",
//              101, -5.05, 5.05);
//      h_trj_nor_eta->SetXTitle("#eta");
//      h_trj_nor_eta->SetYTitle("dn/d#eta (1/0.1)");
// h_trj_nor_phi = new TH1F("trj_nor_phi",
//              "Azimuth Distribution of Track Jets (no overlap removal)",
//              65, -3.2, 3.2);
//      h_trj_nor_phi->SetXTitle("#phi (rad)");
//      h_trj_nor_phi->SetYTitle("dn/d#phi (1/0.1)");
// h_trj_nor_Et = new TH1F("trj_nor_Et",
//              "Transverse-Energy Distribution of Track Jets (no overlap removal)",
//              301, -0.5, 300.5);
//      h_trj_nor_Et->SetXTitle("transverse energy (GeV)");
//      h_trj_nor_Et->SetYTitle("dn/dE_T (1/GeV)");
// h_trj_nor_Et_vs_eta = new TH2F("trj_nor_Et_vs_eta",
//              "Transverse-Energy Distribution of Track Jets versus Rapidity (no overlap removal)",
//              61, -3.05, 3.05, 301, -0.5, 300.5);
//      h_trj_nor_Et->SetZTitle("rapidity");
//      h_trj_nor_Et->SetYTitle("transverse energy (GeV)");
//      h_trj_nor_Et->SetZTitle("d^2n/dE_Td#eta (1/GeV0.1)");

// h_truth_trj_nor_nb_jets = new TH1F("truth_trj_nor_nb_jets",
//              "Number of Track Jets (no overlap removal)",
//              16, -0.5, 15.5);
//      h_truth_trj_nor_nb_jets->SetXTitle("number of jets");
// h_truth_trj_nor_eta = new TH1F("truth_trj_nor_eta",
//              "Pseudorapidity Distribution of Track Jets (no overlap removal)",
//              101, -5.05, 5.05);
//      h_truth_trj_nor_eta->SetXTitle("#eta");
//      h_truth_trj_nor_eta->SetYTitle("dn/d#eta (1/0.1)");
// h_truth_trj_nor_phi = new TH1F("truth_trj_nor_phi",
//              "Azimuth Distribution of Track Jets (no overlap removal)",
//              65, -3.2, 3.2);
//      h_truth_trj_nor_phi->SetXTitle("#phi (rad)");
//      h_truth_trj_nor_phi->SetYTitle("dn/d#phi (1/0.1)");
// h_truth_trj_nor_Et = new TH1F("truth_trj_nor_Et",
//              "Transverse-Energy Distribution of Track Jets (no overlap removal)",
//              301, -0.5, 300.5);
//      h_truth_trj_nor_Et->SetXTitle("transverse energy (GeV)");
//      h_truth_trj_nor_Et->SetYTitle("dn/dE_T (1/GeV)");
// h_truth_trj_nor_Et_vs_eta = new TH2F("truth_trj_nor_Et_vs_eta",
//              "Transverse-Energy Distribution of Track Jets versus Rapidity (no overlap removal)",
//              61, -3.05, 3.05, 301, -0.5, 300.5);
//      h_truth_trj_nor_Et->SetZTitle("rapidity");
//      h_truth_trj_nor_Et->SetYTitle("transverse energy (GeV)");
//      h_truth_trj_nor_Et->SetZTitle("d^2n/dE_Td#eta (1/GeV0.1)");


// h_stdj_nor_nb_jets = new TH1F("stdj_nor_nb_jets",
//              "Number of Standard Jets (no overlap removal)",
//              21, -0.5, 20.5);
//      h_stdj_nor_nb_jets->SetXTitle("number of jets");
// h_stdj_nor_eta = new TH1F("stdj_nor_eta",
//              "Pseudorapidity Distribution of Standard Jets (no overlap removal)",
//              101, -5.05, 5.05);
//      h_stdj_nor_eta->SetXTitle("#eta");
//      h_stdj_nor_eta->SetYTitle("dn/d#eta (1/0.1)");
// h_stdj_nor_phi = new TH1F("stdj_nor_phi",
//              "Azimuth Distribution of Standard Jets (no overlap removal)",
//              65, -3.2, 3.2);
//      h_stdj_nor_phi->SetXTitle("#phi (rad)");
//      h_stdj_nor_phi->SetYTitle("dn/d#phi (1/0.1)");
// h_stdj_nor_Et = new TH1F("stdj_nor_Et",
//              "Transverse-Energy Distribution of Standard Jets (no overlap removal)",
//              301, -0.5, 300.5);
//      h_stdj_nor_Et->SetXTitle("transverse energy (GeV)");
//      h_stdj_nor_Et->SetYTitle("dn/dE_T (1/GeV)");

// h_trj_or_nb_jets = new TH1F("trj_or_nb_jets",
//              "Number of Track Jets (overlap removal)",
//              16, -0.5, 15.5);
//      h_trj_or_nb_jets->SetXTitle("number of jets");
// h_trj_or_eta = new TH1F("trj_or_eta",
//              "Pseudorapidity Distribution of Track Jets (overlap removal)",
//              101, -5.05, 5.05);
//      h_trj_or_eta->SetXTitle("#eta");
//      h_trj_or_eta->SetYTitle("dn/d#eta (1/0.1)");
// h_trj_or_phi = new TH1F("trj_or_phi",
//              "Azimuth Distribution of Track Jets (overlap removal)",
//              65, -3.2, 3.2);
//      h_trj_or_phi->SetXTitle("#phi (rad)");
//      h_trj_or_phi->SetYTitle("dn/d#phi (1/0.1)");
// h_trj_or_Et = new TH1F("trj_or_Et",
//              "Transverse-Energy Distribution of Track Jets (overlap removal)",
//              301, -0.5, 300.5);
//      h_trj_or_Et->SetXTitle("transverse energy (GeV)");
//      h_trj_or_Et->SetYTitle("dn/dE_T (1/GeV)");
// h_trj_or_Et_vs_eta = new TH2F("trj_or_Et_vs_eta",
//              "Transverse-Energy Distribution of Track Jets versus Rapidity (overlap removal)",
//              61, -3.05, 3.05, 301, -0.5, 300.5);
//      h_trj_or_Et->SetZTitle("rapidity");
//      h_trj_or_Et->SetYTitle("transverse energy (GeV)");
//      h_trj_or_Et->SetZTitle("d^2n/dE_Td#eta (1/GeV0.1)");
 h_trj_nb_electron_overlaps = new TH1F("trj_nb_electron_overlaps",
                "Number of Overlaps of Track Jets with Electrons",
                11, -0.5, 10.5);
        h_trj_nb_electron_overlaps->SetXTitle("number of overlaps");
h_trj_nb_muon_overlaps = new TH1F("trj_nb_muon_overlaps",
                "Number of Overlaps of Track Jets with Muons",
                11, -0.5, 10.5);
        h_trj_nb_muon_overlaps->SetXTitle("number of overlaps");

// h_stdj_or_nb_jets = new TH1F("stdj_or_nb_jets",
//              "Number of Standard Jets (overlap removal)",
//              21, -0.5, 20.5);
//      h_stdj_nor_nb_jets->SetXTitle("number of jets");
// h_stdj_or_eta = new TH1F("stdj_or_eta",
//              "Pseudorapidity Distribution of Standard Jets (overlap removal)",
//              101, -5.05, 5.05);
//      h_stdj_nor_eta->SetXTitle("#eta");
//      h_stdj_nor_eta->SetYTitle("dn/d#eta (1/0.1)");
// h_stdj_or_phi = new TH1F("stdj_or_phi",
//              "Azimuth Distribution of Standard Jets (overlap removal)",
//              65, -3.2, 3.2);
//      h_stdj_nor_phi->SetXTitle("#phi (rad)");
//      h_stdj_nor_phi->SetYTitle("dn/d#phi (1/0.1)");
// h_stdj_or_Et = new TH1F("stdj_or_Et",
//              "Transverse-Energy Distribution of Standard Jets (overlap removal)",
//              301, -0.5, 300.5);
//      h_stdj_or_Et->SetXTitle("transverse energy (GeV)");
//      h_stdj_or_Et->SetYTitle("dn/dE_T (1/GeV)");

h_stdj_nb_electron_overlaps = new TH1F("stdj_nb_electron_overlaps",
                "Number of Overlaps of Standard Jets with Electrons",
                11, -0.5, 10.5);
        h_stdj_nb_electron_overlaps->SetXTitle("number of overlaps");
h_stdj_nb_muon_overlaps = new TH1F("stdj_nb_muon_overlaps",
                "Number of Overlaps of Standard Jets with Muons",
                11, -0.5, 10.5);
        h_stdj_nb_muon_overlaps->SetXTitle("number of overlaps");


// matching of track jets with truth and standard jets //
// h_trj_truth_jet_match = new TH1F("trj_truth_jet_match",
//                                  "Fraction of Track Jets Matched with Truth Jets",
//                                  2, -0.5, 1.5);

// h_trj_truth_jet_match_vs_eta = new TH2F("trj_truth_jet_match_vs_eta",
//                                         "Fraction of Track Jets Matched with Truth Jets vs eta",
//                                         65, -3.2, 3.2, 2, -0.5, 1.5);
        
// h_trj_truth_jet_match_vs_pt = new TH2F("trj_truth_jet_match_vs_pt",
//                                        "Fraction of Track Jets Matched with Truth Jets vs eta",
//                                        ptbins, ptmin, ptmax, 2, -0.5, 1.5);
        

// h_trj_std_jet_match = new TH1F("trj_std_jet_match",
//                                "Fraction of Track Jets Matched with Standard Jets",
//                                2, -0.5, 1.5);

// h_trj_std_jet_match_vs_eta = new TH2F("trj_std_jet_match_vs_eta",
//                                       "Fraction of Track Jets Matched with Standard Jets vs eta",
//                                       65, -3.2, 3.2, 2, -0.5, 1.5);

// h_trj_std_jet_match_vs_pt = new TH2F("trj_std_jet_match_vs_pt",
//                                      "Fraction of Track Jets Matched with Standard Jets vs pt",
//                                      ptbins, ptmin, ptmax, 2, -0.5, 1.5);



// h_std_trj_jet_match = new TH1F("std_trj_jet_match",
//                                "Fraction of Standard Jets Matched with Track Jets",
//                                2, -0.5, 1.5);
// h_std_trj_jet_match_vs_eta = new TH2F("std_trj_jet_match_vs_eta",
//                                       "Fraction of Standard Jets Matched with Track Jets vs eta",
//                                       65, -3.2, 3.2, 2, -0.5, 1.5);
// h_std_trj_jet_match_vs_pt = new TH2F("std_trj_jet_match_vs_pt",
//                                      "Fraction of Standard Jets Matched with Track Jets vs pt",
//                                      ptbins, ptmin, ptmax, 2, -0.5, 1.5);



// h_trj_std_jet_match_pt_ratio = new TH1F("trj_std_jet_match_pt_ratio",
//                                         "pt ratio", 60, 0, 2);

// h_trj_std_jet_match_2Dpt = new TH2F("trj_std_jet_match_2Dpt",
//                                     "pt reco vs pt track jet", 
//                                     ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);

// h_jet_reco_notmatched_pt = new TH1F("jet_reco_notmatched_pt","jet_reco_notmatched_pt", ptbins, ptmin, ptmax);
// h_jet_reco_notmatched_eta = new TH1F("jet_reco_notmatched_eta","jet_reco_notmatched_eta", etabins, etamin, etamax);
// h_jet_reco_notmatched_phi = new TH1F("jet_reco_notmatched_phi","jet_reco_notmatched_phi", phibins, phimin, phimax);
// h_jet_reco_notmatched_ntrks = new TH1F("jet_reco_notmatched_ntrks","jet_reco_notmatched_trks", 20, 0, 20);


//h_cut_evolution = new TH1F("cut_evolution","nevents after cuts",22,-2,20);

h_event_weight = new TH1F("event_weight","event weight",10,-5,5);
h_offline_trigger = new TH1F("offline_trigger","offline trigger",6,0,6);

h_cut_events = new TH1F("cut_events","nevents after cuts",nb_cuts_max+2,-2,nb_cuts_max);
h_cut_candidates = new TH1F("cut_candidates","ncandidates after cuts",nb_cuts_max+2,-2,nb_cuts_max);

h_multi_tagjet_pairs_nb = new TH1F("multi_tagjet_pairs_nb","nb multiple tagjet pairs",15,0,15);
h_multi_lepton_pairs_nb = new TH1F("multi_lepton_pairs_nb","nb multiple lepton pairs",15,0,15);

// raw aod content histos
h_etmiss_aod = new TH1F("etmiss_aod","etmiss_aod", 100, ptmin, ptmax);
        
h_e_reco_aod_nb = new TH1F("e_reco_aod_nb","e_reco_aod_nb",19, -1, 18);
h_e_reco_aod_pt = new TH1F("e_reco_aod_pt","e_reco_aod_pt", ptbins, ptmin, ptmax);
h_e_reco_aod_eta = new TH1F("e_reco_aod_eta","e_reco_aod_eta", etabins, etamin, etamax);
h_e_reco_aod_phi = new TH1F("e_reco_aod_phi","e_reco_aod_phi", phibins, phimin, phimax);
h_mu_reco_aod_nb = new TH1F("mu_reco_aod_nb","mu_reco_aod_nb",19, -1, 18);
h_mu_reco_aod_pt = new TH1F("mu_reco_aod_pt","mu_reco_aod_pt", ptbins, ptmin, ptmax);
h_mu_reco_aod_eta = new TH1F("mu_reco_aod_eta","mu_reco_aod_eta", etabins, etamin, etamax);
h_mu_reco_aod_phi = new TH1F("mu_reco_aod_phi","mu_reco_aod_phi", phibins, phimin, phimax);
h_mu_reco_aod_fitChi2          = new TH1F("mu_reco_aod_fitChi2","mu_reco_aod_fitChi2",100, -1, 50);
h_mu_reco_aod_fitChi2OverDoF   = new TH1F("mu_reco_aod_fitChi2OverDoF","mu_reco_aod_fitChi2OverDoF",100, -1, 9);
h_mu_reco_aod_matchChi2        = new TH1F("mu_reco_aod_matchChi2","mu_reco_aod_matchChi2",100, -1, 50);
h_mu_reco_aod_matchChi2OverDoF = new TH1F("mu_reco_aod_matchChi2OverDoF","mu_reco_aod_matchChi2OverDoF",100, -1, 9);
h_jet_reco_aod_nb = new TH1F("jet_reco_aod_nb","jet_reco_aod_nb",19, -1, 18);
h_jet_reco_aod_pt = new TH1F("jet_reco_aod_pt","jet_reco_aod_pt", ptbins, ptmin, ptmax);
h_jet_reco_aod_eta = new TH1F("jet_reco_aod_eta","jet_reco_aod_eta", etabins, etamin, etamax);
h_jet_reco_aod_phi = new TH1F("jet_reco_aod_phi","jet_reco_aod_phi", phibins, phimin, phimax);
h_e_truth_aod_nb = new TH1F("e_truth_aod_nb","e_truth_aod_nb",19, -1, 18);
h_e_truth_aod_pt = new TH1F("e_truth_aod_pt","e_truth_aod_pt", ptbins, ptmin, ptmax);
h_e_truth_aod_eta = new TH1F("e_truth_aod_eta","e_truth_aod_eta", etabins, etamin, etamax);
h_e_truth_aod_phi = new TH1F("e_truth_aod_phi","e_truth_aod_phi", phibins, phimin, phimax);
h_mu_truth_aod_nb = new TH1F("mu_truth_aod_nb","mu_truth_aod_nb",19, -1, 18);
h_mu_truth_aod_pt = new TH1F("mu_truth_aod_pt","mu_truth_aod_pt", ptbins, ptmin, ptmax);
h_mu_truth_aod_eta = new TH1F("mu_truth_aod_eta","mu_truth_aod_eta", etabins, etamin, etamax);
h_mu_truth_aod_phi = new TH1F("mu_truth_aod_phi","mu_truth_aod_phi", phibins, phimin, phimax);
h_jet_truth_aod_nb = new TH1F("jet_truth_aod_nb","jet_truth_aod_nb",19, -1, 18);
h_jet_truth_aod_pt = new TH1F("jet_truth_aod_pt","jet_truth_aod_pt", ptbins, ptmin, ptmax);
h_jet_truth_aod_eta = new TH1F("jet_truth_aod_eta","jet_truth_aod_eta", etabins, etamin, etamax);
h_jet_truth_aod_phi = new TH1F("jet_truth_aod_phi","jet_truth_aod_phi", phibins, phimin, phimax);

h_track_jet_reco_aod_nb = new TH1F("track_jet_reco_aod_nb","track_jet_reco_aod_nb",19, -1, 18);
h_track_jet_reco_aod_pt = new TH1F("track_jet_reco_aod_pt","track_jet_reco_aod_pt", ptbins, ptmin, ptmax);
h_track_jet_reco_aod_eta = new TH1F("track_jet_reco_aod_eta","track_jet_reco_aod_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_reco_aod_phi = new TH1F("track_jet_reco_aod_phi","track_jet_reco_aod_phi", phibins, phimin, phimax);
h_track_jet_reco_aod_pt_vs_eta = new TH2F("track_jet_reco_aod_pt_vs_eta", "track_jet_reco_aod_pt_vs_eta",
                                          etabins, etamin+2, etamax-2, ptbins, ptmin, ptmax);
h_track_jet_truth_aod_nb = new TH1F("track_jet_truth_aod_nb","track_jet_truth_aod_nb",19, -1, 18);
h_track_jet_truth_aod_pt = new TH1F("track_jet_truth_aod_pt","track_jet_truth_aod_pt", ptbins, ptmin, ptmax);
h_track_jet_truth_aod_eta = new TH1F("track_jet_truth_aod_eta","track_jet_truth_aod_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_truth_aod_phi = new TH1F("track_jet_truth_aod_phi","track_jet_truth_aod_phi", phibins, phimin, phimax);
h_track_jet_truth_aod_pt_vs_eta = new TH2F("track_jet_truth_aod_pt_vs_eta", "track_jet_truth_aod_pt_vs_eta",
                                           etabins, etamin+2, etamax-2, ptbins, ptmin, ptmax);

h_track_jet_reco_presel_nb = new TH1F("track_jet_reco_presel_nb","track_jet_reco_presel_nb",19, -1, 18);
h_track_jet_reco_presel_pt = new TH1F("track_jet_reco_presel_pt","track_jet_reco_presel_pt", ptbins, ptmin, ptmax);
h_track_jet_reco_presel_eta = new TH1F("track_jet_reco_presel_eta","track_jet_reco_presel_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_reco_presel_phi = new TH1F("track_jet_reco_presel_phi","track_jet_reco_presel_phi", phibins, phimin, phimax);
h_track_jet_reco_presel_pt_vs_eta = new TH2F("track_jet_reco_presel_pt_vs_eta", "track_jet_reco_presel_pt_vs_eta",
                                          etabins, etamin+2, etamax-2, ptbins, ptmin, ptmax);
h_track_jet_truth_presel_nb = new TH1F("track_jet_truth_presel_nb","track_jet_truth_presel_nb",19, -1, 18);
h_track_jet_truth_presel_pt = new TH1F("track_jet_truth_presel_pt","track_jet_truth_presel_pt", ptbins, ptmin, ptmax);
h_track_jet_truth_presel_eta = new TH1F("track_jet_truth_presel_eta","track_jet_truth_presel_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_truth_presel_phi = new TH1F("track_jet_truth_presel_phi","track_jet_truth_presel_phi", phibins, phimin, phimax);
h_track_jet_truth_presel_pt_vs_eta = new TH2F("track_jet_truth_presel_pt_vs_eta", "track_jet_truth_presel_pt_vs_eta",
                                           etabins, etamin+2, etamax-2, ptbins, ptmin, ptmax);
h_track_jet_reco_orm_nb = new TH1F("track_jet_reco_orm_nb","track_jet_reco_orm_nb",19, -1, 18);
h_track_jet_reco_orm_pt = new TH1F("track_jet_reco_orm_pt","track_jet_reco_orm_pt", ptbins, ptmin, ptmax);
h_track_jet_reco_orm_eta = new TH1F("track_jet_reco_orm_eta","track_jet_reco_orm_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_reco_orm_phi = new TH1F("track_jet_reco_orm_phi","track_jet_reco_orm_phi", phibins, phimin, phimax);
h_track_jet_reco_orm_pt_vs_eta = new TH2F("track_jet_reco_orm_pt_vs_eta", "track_jet_reco_orm_pt_vs_eta",
                                          etabins_central, etamin_central, etamax_central, ptbins, ptmin, ptmax);
h_track_jet_truth_orm_nb = new TH1F("track_jet_truth_orm_nb","track_jet_truth_orm_nb",19, -1, 18);
h_track_jet_truth_orm_pt = new TH1F("track_jet_truth_orm_pt","track_jet_truth_orm_pt", ptbins, ptmin, ptmax);
h_track_jet_truth_orm_eta = new TH1F("track_jet_truth_orm_eta","track_jet_truth_orm_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_truth_orm_phi = new TH1F("track_jet_truth_orm_phi","track_jet_truth_orm_phi", phibins, phimin, phimax);
h_track_jet_truth_orm_pt_vs_eta = new TH2F("track_jet_truth_orm_pt_vs_eta", "track_jet_truth_orm_pt_vs_eta",
                                           etabins_central, etamin_central, etamax_central, ptbins, ptmin, ptmax);
// matching, efficiency histos
// electrons
h_e_reco_presel_pt = new TH1F("e_reco_presel_pt","e_reco_presel_pt", ptbins, ptmin, ptmax);
h_e_reco_presel_eta = new TH1F("e_reco_presel_eta","e_reco_presel_eta", etabins, etamin+2.5, etamax-2.5);
h_e_reco_presel_phi = new TH1F("e_reco_presel_phi","e_reco_presel_phi", phibins, phimin, phimax);
h_e_reco_presel_nb = new TH1F("e_reco_presel_nb","e_reco_presel_nb",19, -1, 18);
h_e_reco_presel_isem = new TH1F("e_reco_presel_isem","e_reco_presel_isem",4001,0.0,4000.0);
        
h_e_truth_presel_pt = new TH1F("e_truth_presel_pt","e_truth_presel_pt", ptbins, ptmin, ptmax);
h_e_truth_presel_eta = new TH1F("e_truth_presel_eta","e_truth_presel_eta", etabins, etamin+2.5, etamax-2.5);
h_e_truth_presel_phi = new TH1F("e_truth_presel_phi","e_truth_presel_phi", phibins, phimin, phimax);
h_e_truth_presel_nb = new TH1F("e_truth_presel_nb","e_truth_presel_nb",19, -1, 18);
h_e_truth_presel_status = new TH1F("e_truth_presel_status",
                                   "e_truth_presel_status", 200, -1, 198);
        
h_e_reco_matched_pt = new TH1F("e_reco_matched_pt","e_reco_matched_pt", ptbins, ptmin, ptmax);
h_e_reco_matched_eta = new TH1F("e_reco_matched_eta","e_reco_matched_eta", etabins, etamin+2.5, etamax-2.5);
h_e_reco_matched_phi = new TH1F("e_reco_matched_phi","e_reco_matched_phi", phibins, phimin, phimax);
h_e_reco_matched_nb = new TH1F("e_reco_matched_nb","e_reco_matched_nb",19, -1, 18);
h_e_reco_matched_isem = new TH1F("e_reco_matched_isem","e_reco_matched_isem",4001,0.0,4000.0);
h_e_reco_matched_NN = new TH1F("e_reco_matched_NN","e_reco_matched_NN",100,-0.1,1.1);
        
h_e_truth_matched_pt = new TH1F("e_truth_matched_pt","e_truth_matched_pt", ptbins, ptmin, ptmax);
h_e_truth_matched_eta = new TH1F("e_truth_matched_eta","e_truth_matched_eta", etabins, etamin+2.5, etamax-2.5);
h_e_truth_matched_phi = new TH1F("e_truth_matched_phi","e_truth_matched_phi", phibins, phimin, phimax);
        
h_e_reco_fake_pt = new TH1F("e_reco_fake_pt","e_reco_fake_pt", ptbins, ptmin, ptmax);
h_e_reco_fake_eta = new TH1F("e_reco_fake_eta","e_reco_fake_eta", etabins, etamin+2.5, etamax-2.5);
h_e_reco_fake_phi = new TH1F("e_reco_fake_phi","e_reco_fake_phi", phibins, phimin, phimax);
h_e_reco_fake_isem = new TH1F("e_reco_fake_isem","e_reco_fake_isem",4001,0.0,4000.0);
h_e_reco_fake_nb = new TH1F("e_reco_fake_nb","e_reco_fake_nb",19, -1, 18);
        
h_e_truth_notmatched_pt = new TH1F("e_truth_notmatched_pt","e_truth_notmatched_pt",  ptbins, ptmin, ptmax);
h_e_truth_notmatched_eta = new TH1F("e_truth_notmatched_eta","e_truth_notmatched_eta",  etabins, etamin-2.5, etamax+2.5);
h_e_truth_notmatched_phi = new TH1F("e_truth_notmatched_phi","e_truth_notmatched_phi",  phibins, phimin, phimax);
h_e_truth_notmatched_nb = new TH1F("e_truth_notmatched_nb","e_truth_notmatched_nb",19, -1, 18);
        
h_e_matched_dpt = new TH1F("e_matched_dpt","e_matched_dpt", 200, -d_ptmax, d_ptmax);
h_e_matched_deta = new TH1F("e_matched_deta","e_matched_deta", 200, -d_etamax, d_etamax);
h_e_matched_dphi = new TH1F("e_matched_dphi","e_matched_dphi", 200, -d_phimax, d_phimax);
h_e_matched_dr = new TH1F("e_matched_dr","e_matched_dr", 200, 0.0, 0.1);
        
h_e_matched_resolution_pt = new TH2F("e_matched_resolution_pt","e_matched_resolution_pt",
                                     20, ptmin, ptmax, 100, -1.0, 1.0);
h_e_matched_resolution_eta = new TH2F("e_matched_resolution_eta","e_matched_resolution_eta",
                                      20, etamin-2.5, etamax+2.5, 100, -1.0, 1.0);
h_e_matched_resolution_phi = new TH2F("e_matched_resolution_phi","e_matched_resolution_phi",
                                      20, phimin, phimax, 100, -1.0, 1.0);

// muons
h_mu_reco_presel_pt = new TH1F("mu_reco_presel_pt","mu_reco_presel_pt", ptbins, ptmin, ptmax);
h_mu_reco_presel_eta = new TH1F("mu_reco_presel_eta","mu_reco_presel_eta", etabins, etamin+2.5, etamax-2.5);
h_mu_reco_presel_phi = new TH1F("mu_reco_presel_phi","mu_reco_presel_phi", phibins, phimin, phimax);
h_mu_reco_presel_nb = new TH1F("mu_reco_presel_nb","mu_reco_presel_nb",19, -1, 18);
h_mu_reco_presel_fitChi2OverDoF = new TH1F("mu_reco_presel_fitChi2OverDoF","mu_reco_presel_fitChi2OverDoF",100, -1, 9);
h_mu_reco_presel_fitChi2OverDoFvsphi = new TH2F("mu_reco_presel_fitChi2OverDoFvsphi","mu_reco_presel_fitChi2OverDoFvsphi",
                                                40, phimin, phimax, 100, -1.0, 9.0);
        
h_mu_truth_presel_pt = new TH1F("mu_truth_presel_pt","mu_truth_presel_pt", ptbins, ptmin, ptmax);
h_mu_truth_presel_eta = new TH1F("mu_truth_presel_eta","mu_truth_presel_eta", etabins, etamin+2.5, etamax-2.5);
h_mu_truth_presel_phi = new TH1F("mu_truth_presel_phi","mu_truth_presel_phi", phibins, phimin, phimax);
h_mu_truth_presel_nb = new TH1F("mu_truth_presel_nb","mu_truth_presel_nb",19, -1, 18);
h_mu_truth_presel_status = new TH1F("mu_truth_presel_status",
                                    "mu_truth_presel_status", 100, -5, 95);
        
h_mu_reco_matched_pt = new TH1F("mu_reco_matched_pt","mu_reco_matched_pt", ptbins, ptmin, ptmax);
h_mu_reco_matched_eta = new TH1F("mu_reco_matched_eta","mu_reco_matched_eta", etabins, etamin+2.5, etamax-2.5);
h_mu_reco_matched_phi = new TH1F("mu_reco_matched_phi","mu_reco_matched_phi", phibins, phimin, phimax);
h_mu_reco_matched_nb = new TH1F("mu_reco_matched_nb","mu_reco_matched_nb",19, -1, 18);
h_mu_reco_matched_fitChi2OverDoF = new TH1F("mu_reco_matched_fitChi2OverDoF","mu_reco_matched_fitChi2OverDoF",100, -1, 9);
        
h_mu_truth_matched_pt = new TH1F("mu_truth_matched_pt","mu_truth_matched_pt", ptbins, ptmin, ptmax);
h_mu_truth_matched_eta = new TH1F("mu_truth_matched_eta","mu_truth_matched_eta", etabins, etamin+2.5, etamax-2.5);
h_mu_truth_matched_phi = new TH1F("mu_truth_matched_phi","mu_truth_matched_phi", phibins, phimin, phimax);
        
h_mu_reco_fake_pt = new TH1F("mu_reco_fake_pt","mu_reco_fake_pt", ptbins, ptmin, ptmax);
h_mu_reco_fake_eta = new TH1F("mu_reco_fake_eta","mu_reco_fake_eta", etabins, etamin+2.5, etamax-2.5);
h_mu_reco_fake_phi = new TH1F("mu_reco_fake_phi","mu_reco_fake_phi", phibins, phimin, phimax);
h_mu_reco_fake_nb = new TH1F("mu_reco_fake_nb","mu_reco_fake_nb",19, -1, 18);
h_mu_reco_fake_fitChi2OverDoF = new TH1F("mu_reco_fake_fitChi2OverDoF","mu_reco_fake_fitChi2OverDoF",100, -1, 9);
        
h_mu_truth_notmatched_pt = new TH1F("mu_truth_notmatched_pt","mu_truth_notmatched_pt",  ptbins, ptmin, ptmax);
h_mu_truth_notmatched_eta = new TH1F("mu_truth_notmatched_eta","mu_truth_notmatched_eta",  etabins, etamin+2.5, etamax-2.5);
h_mu_truth_notmatched_phi = new TH1F("mu_truth_notmatched_phi","mu_truth_notmatched_phi",  phibins, phimin, phimax);
h_mu_truth_notmatched_nb = new TH1F("mu_truth_notmatched_nb","mu_truth_notmatched_nb",19, -1, 18);
        
h_mu_matched_dpt = new TH1F("mu_matched_dpt","mu_matched_dpt", 200, -d_ptmax, d_ptmax);
h_mu_matched_deta = new TH1F("mu_matched_deta","mu_matched_deta", 200, -d_etamax, d_etamax);
h_mu_matched_dphi = new TH1F("mu_matched_dphi","mu_matched_dphi", 200, -d_phimax, d_phimax);
h_mu_matched_dr = new TH1F("mu_matched_dr","mu_matched_dr", 200, 0.0, 0.1);
        
h_mu_matched_resolution_pt = new TH2F("mu_matched_resolution_pt","mu_matched_resolution_pt",
                                      20, ptmin, ptmax, 100, -1.0, 1.0);
h_mu_matched_resolution_ptvsphi = new TH2F("mu_matched_resolution_ptvsphi","mu_matched_resolution_ptvsphi",
                                           40, phimin, phimax, 100, -1.0, 1.0);
h_mu_matched_resolution_ptvseta = new TH2F("mu_matched_resolution_ptvseta","mu_matched_resolution_ptvseta",
                                           20, etamin+2.5, etamax-2.5, 100, -1.0, 1.0);
h_mu_matched_resolution_eta = new TH2F("mu_matched_resolution_eta","mu_matched_resolution_eta",
                                       20, etamin+2.5, etamax-2.5, 100, -1.0, 1.0);
h_mu_matched_resolution_phi = new TH2F("mu_matched_resolution_phi","mu_matched_resolution_phi",
                                       20, phimin, phimax, 100, -1.0, 1.0);
        

// jets
// after overlap removal with electrons
h_jet_reco_orm_pt = new TH1F("jet_reco_orm_pt","jet_reco_orm_pt", ptbins, ptmin, ptmax);
h_jet_reco_orm_eta = new TH1F("jet_reco_orm_eta","jet_reco_orm_eta", etabins, etamin, etamax);
h_jet_reco_orm_phi = new TH1F("jet_reco_orm_phi","jet_reco_orm_phi", phibins, phimin, phimax);
h_jet_reco_orm_nb = new TH1F("jet_reco_orm_nb","jet_reco_orm_nb",19, -1, 18);

h_jet_reco_central_pt = new TH1F("jet_reco_central_pt","jet_reco_central_pt", ptbins, ptmin, ptmax);
h_jet_reco_central_eta = new TH1F("jet_reco_central_eta","jet_reco_central_eta", etabins_central, etamin_central, etamax_central);
h_jet_reco_central_phi = new TH1F("jet_reco_central_phi","jet_reco_central_phi", phibins, phimin, phimax);
h_jet_reco_central_nb = new TH1F("jet_reco_central_nb","jet_reco_central_nb",19, -1, 18);
h_jet_reco_central_ptratio = new TH1F("jet_reco_central_ptratio","jet_reco_central_ptratio",100, 0,2);
h_jet_reco_central_cmsptratio = new TH1F("jet_reco_central_cmsptratio","jet_reco_central_cmsptratio",100, 0,2);
h_jet_reco_central_ptcorr  = new TH1F("jet_reco_central_ptcorr","jet_reco_central_ptcorr", ptbins, ptmin, ptmax);
h_jet_reco_central_pt_ptratiocut = new TH1F("jet_reco_central_pt_ptratiocut","jet_reco_central_pt_ptratiocut", ptbins, ptmin, ptmax);
h_jet_reco_central_pt_cmsptratiocut = new TH1F("jet_reco_central_pt_cmsptratiocut","jet_reco_central_pt_cmsptratiocut", ptbins, ptmin, ptmax);




h_jet_reco_presel_pt = new TH1F("jet_reco_presel_pt","jet_reco_presel_pt", ptbins, ptmin, ptmax);
h_jet_reco_presel_eta = new TH1F("jet_reco_presel_eta","jet_reco_presel_eta", etabins, etamin, etamax);
h_jet_reco_presel_phi = new TH1F("jet_reco_presel_phi","jet_reco_presel_phi", phibins, phimin, phimax);
h_jet_reco_presel_nb = new TH1F("jet_reco_presel_nb","jet_reco_presel_nb",19, -1, 18);
h_jet_reco_presel_ptratio = new TH1F("jet_reco_presel_ptratio","jet_reco_presel_ptratio",100, 0,2);
h_jet_reco_presel_cmsptratio = new TH1F("jet_reco_presel_cmsptratio","jet_reco_presel_cmsptratio",100, 0,2);
h_jet_reco_presel_ptcorr = new TH1F("jet_reco_presel_ptcorr","jet_reco_presel_ptcorr", ptbins, ptmin, ptmax);


h_jet_reco_presel_pt1 = new TH1F("jet_reco_presel_pt1","jet_reco_presel_pt1", ptbins, ptmin, ptmax);
h_jet_reco_presel_pt2 = new TH1F("jet_reco_presel_pt2","jet_reco_presel_pt2", ptbins, ptmin, ptmax);
h_jet_reco_presel_pt3 = new TH1F("jet_reco_presel_pt3","jet_reco_presel_pt3", ptbins, ptmin, ptmax);
h_jet_reco_presel_eta1 = new TH1F("jet_reco_presel_eta1","jet_reco_presel_eta1", etabins, etamin, etamax);
h_jet_reco_presel_eta2 = new TH1F("jet_reco_presel_eta2","jet_reco_presel_eta2", etabins, etamin, etamax);
h_jet_reco_presel_eta3 = new TH1F("jet_reco_presel_eta3","jet_reco_presel_eta3", etabins, etamin, etamax);

h_jet_truth_presel_pt = new TH1F("jet_truth_presel_pt","jet_truth_presel_pt", ptbins, ptmin, ptmax);
h_jet_truth_presel_eta = new TH1F("jet_truth_presel_eta","jet_truth_presel_eta", etabins, etamin, etamax);
h_jet_truth_presel_phi = new TH1F("jet_truth_presel_phi","jet_truth_presel_phi", phibins, phimin, phimax);
h_jet_truth_presel_nb = new TH1F("jet_truth_presel_nb","jet_truth_presel_nb",19, -1, 18);

h_jet_truth_central_pt = new TH1F("jet_truth_central_pt","jet_truth_central_pt", ptbins, ptmin, ptmax);
h_jet_truth_central_eta = new TH1F("jet_truth_central_eta","jet_truth_central_eta", etabins, etamin, etamax);
h_jet_truth_central_phi = new TH1F("jet_truth_central_phi","jet_truth_central_phi", phibins, phimin, phimax);
h_jet_truth_central_nb = new TH1F("jet_truth_central_nb","jet_truth_central_nb",19, -1, 18);
   


h_jet_reco_matched_pt = new TH1F("jet_reco_matched_pt","jet_reco_matched_pt", ptbins, ptmin, ptmax);
h_jet_reco_matched_eta = new TH1F("jet_reco_matched_eta","jet_reco_matched_eta", etabins, etamin, etamax);
h_jet_reco_matched_phi = new TH1F("jet_reco_matched_phi","jet_reco_matched_phi", phibins, phimin, phimax);
h_jet_reco_matched_nb = new TH1F("jet_reco_matched_nb","jet_reco_matched_nb",19, -1, 18);
h_jet_reco_matched_E = new TH1F("jet_reco_matched_E","jet_reco_matched_E",100, Emin, Emax);

h_jet_truth_matched_pt = new TH1F("jet_truth_matched_pt","jet_truth_matched_pt", ptbins, ptmin, ptmax);
h_jet_truth_matched_eta = new TH1F("jet_truth_matched_eta","jet_truth_matched_eta", etabins, etamin, etamax);
h_jet_truth_matched_phi = new TH1F("jet_truth_matched_phi","jet_truth_matched_phi", phibins, phimin, phimax);
h_jet_truth_matched_E = new TH1F("jet_truth_matched_E","jet_truth_matched_E",100, Emin, Emax);

h_jet_reco_fake_pt = new TH1F("jet_reco_fake_pt","jet_reco_fake_pt", ptbins, ptmin, ptmax);
h_jet_reco_fake_eta = new TH1F("jet_reco_fake_eta","jet_reco_fake_eta", etabins, etamin, etamax);
h_jet_reco_fake_phi = new TH1F("jet_reco_fake_phi","jet_reco_fake_phi", phibins, phimin, phimax);
h_jet_reco_fake_nb = new TH1F("jet_reco_fake_nb","jet_reco_fake_nb",19, -1, 18);
        
h_jet_truth_notmatched_pt = new TH1F("jet_truth_notmatched_pt","jet_truth_notmatched_pt",  ptbins, ptmin, ptmax);
h_jet_truth_notmatched_eta = new TH1F("jet_truth_notmatched_eta","jet_truth_notmatched_eta",  etabins, etamin, etamax);
h_jet_truth_notmatched_phi = new TH1F("jet_truth_notmatched_phi","jet_truth_notmatched_phi",  phibins, phimin, phimax);
h_jet_truth_notmatched_nb = new TH1F("jet_truth_notmatched_nb","jet_truth_notmatched_nb",19, -1, 18);
        
h_jet_matched_dpt = new TH1F("jet_matched_dpt","jet_matched_dpt", 200, -d_ptmax, d_ptmax);
h_jet_matched_deta = new TH1F("jet_matched_deta","jet_matched_deta", 200, -d_etamax, d_etamax);
h_jet_matched_dphi = new TH1F("jet_matched_dphi","jet_matched_dphi", 200, -d_phimax, d_phimax);
h_jet_matched_dr = new TH1F("jet_matched_dr","jet_matched_dr", 200, 0.0, 0.1);
h_jet_matched_dE = new TH1F("jet_matched_dE","jet_matched_dE", 200, -d_Emax, d_Emax);

h_jet_matched_ptratio_recoovertruth = new TH1F("jet_matched_ptratio_recoovertruth",
                                               "jet_matched_ptratio_recoovertruth", 
                                               60, 0, 2);
h_jet_matched_ptratio_recoovertruth_corr = new TH1F("jet_matched_ptratio_recoovertruth_corr",
                                                    "jet_matched_ptratio_recoovertruth_corr", 
                                                    60, 0, 2);
h_jet_matched_ptratio_truthoverreco = new TH1F("jet_matched_ptratio_truthoverreco",
                                               "jet_matched_ptratio_truthoverreco", 
                                               60, 0, 2);
h_jet_matched_ptratio_truthoverreco_corr = new TH1F("jet_matched_ptratio_truthoverreco_corr",
                                                    "jet_matched_ptratio_truthoverreco_corr", 
                                                    60, 0, 2);

h_jet_matched_bweightcut_ptratio_recoovertruth = new TH1F("jet_matched_bweightcut_ptratio_recoovertruth",
                                                          "jet_matched_bweightcut_ptratio_recoovertruth", 
                                                          60, 0, 2);
h_jet_matched_bweightcut_ptratio_recoovertruth_corr = new TH1F("jet_matched_bweightcut_ptratio_recoovertruth_corr",
                                                               "jet_matched_bweightcut_ptratio_recoovertruth_corr", 
                                                               60, 0, 2);
h_jet_matched_bweightcut_ptratio_truthoverreco = new TH1F("jet_matched_bweightcut_ptratio_truthoverreco",
                                                          "jet_matched_bweightcut_ptratio_truthoverreco", 
                                                          60, 0, 2);
h_jet_matched_bweightcut_ptratio_truthoverreco_corr = new TH1F("jet_matched_bweightcut_ptratio_truthoverreco_corr",
                                                               "jet_matched_bweightcut_ptratio_truthoverreco_corr", 
                                                               60, 0, 2);

h_jet_matched_resolution_pt = new TH2F("jet_matched_resolution_pt","jet_matched_resolution_pt",
                                       20, ptmin, ptmax, 100, -1.0, 1.0);
h_jet_matched_resolution_eta = new TH2F("jet_matched_resolution_eta","jet_matched_resolution_eta",
                                        20, etamin, etamax, 100, -1.0, 1.0);
h_jet_matched_resolution_phi = new TH2F("jet_matched_resolution_phi","jet_matched_resolution_phi",
                                        20, phimin, phimax, 100, -1.0, 1.0);

//=========
   
h_jet_matched_Eresolution = new TH1F("jet_matched_Eresolution", 
                                     "E RESOLUTION MATCHED JETS (GeV)", 
                                     80, -200, 200);
   
h_jet_matched_Eresolution_vs_pt = new TH2F("jet_matched_Eresolution_vs_pt", 
                                           "E RES MATCHED JETS (GeV) VS PT (GeV)", 
                                           ptbins, ptmin, ptmax, 80, -200, 200);
   
h_jet_matched_Eresolution_vs_eta = new TH2F("jet_matched_Eresolution_vs_eta", 
                                            "E RES MATCHED JETS (GeV) VS ETA", 
                                            20, etamin, etamax, 80, -200, 200);
   
h_jet_matched_Eresolution_vs_phi = new TH2F("jet_matched_Eresolution_vs_phi", 
                                            "E RES MATCHED JETS (GeV) VS PHI", 
                                            14, phimin, phimax, 80, -200, 200);
   
h_jet_matched_Eresolution_vs_E = new TH2F("jet_matched_Eresolution_vs_E", 
                                          "E RES MATCHED JETS (GeV) VS E (GeV)", 
                                          Ebins, Emin, Emax, 80, -200, 200);
//=========
   
h_jet_matched_relEresolution = new TH1F("jet_matched_relEresolution", 
                                        "Rel E RESOLUTION MATCHED JETS (GeV)", 
                                        80, -1, 1);
   
h_jet_matched_relEresolution_vs_pt = new TH2F("jet_matched_relEresolution_vs_pt", 
                                              "Rel E RES MATCHED JETS VS PT (GeV)", 
                                              ptbins, ptmin, ptmax, 80, -1, 1);
   
h_jet_matched_relEresolution_vs_eta = new TH2F("jet_matched_relEresolution_vs_eta", 
                                               "Rel E RES MATCHED JETS VS ETA", 
                                               20, etamin, etamax, 80, -1, 1);
   
h_jet_matched_relEresolution_vs_phi = new TH2F("jet_matched_relEresolution_vs_phi", 
                                               "Rel E RES MATCHED JETS VS PHI", 
                                               14, phimin, phimax, 80, -1, 1);

h_jet_matched_relEresolution_vs_E = new TH2F("jet_matched_relEresolution_vs_E", 
                                             "Rel E RES MATCHED JETS VS E (GeV)", 
                                             Ebins, Emin, Emax, 80, -1, 1);
   

h_jet_matched_invptresolution = new TH1F("jet_matched_invptresolution", 
                                         "jet_matched_invptresolution",
                                         80, -0.02, 0.02);
h_jet_matched_relinvptresolution = new TH1F("jet_matched_relinvptresolution", 
                                            "jet_matched_relinvptresolution",
                                            80, -1, 1);
h_jet_matched_invEresolution = new TH1F("jet_matched_invEresolution", 
                                        "jet_matched_invEresolution",
                                        80, -0.02, 0.02);
h_jet_matched_relinvEresolution = new TH1F("jet_matched_relinvEresolution", 
                                           "jet_matched_relinvEresolution",
                                           80, -1, 1);
h_jet_matched_invptresolution_vs_pt = new TH2F("jet_matched_invptresolution_vs_pt", 
                                               "jet_matched_invptresolution_vs_pt", 
                                               ptbins, ptmin, ptmax, 80, -0.02, 0.02);
h_jet_matched_relinvptresolution_vs_pt = new TH2F("jet_matched_relinvptresolution_vs_pt", 
                                                  "jet_matched_relinvptresolution_vs_pt", 
                                                  ptbins, ptmin, ptmax, 80, -1, 1);
h_jet_matched_invEresolution_vs_E = new TH2F("jet_matched_invEresolution_vs_E", 
                                             "jet_matched_invEresolution_vs_E", 
                                             Ebins, Emin, Emax, 80, -0.02, 0.02);
h_jet_matched_relinvEresolution_vs_E = new TH2F("jet_matched_relinvEresolution_vs_E", 
                                                "jet_matched_relinvEresolution_vs_E", 
                                                Ebins, Emin, Emax, 80, -1, 1);

//=========
// track jets
h_track_jet_fake_id  = new TH1F("track_jet_fake_id","track_jet_fake_id", 15, 0, 15);

h_track_jet_reco_matched_pt  = new TH1F("track_jet_reco_matched_pt","track_jet_reco_matched_pt", ptbins, ptmin, ptmax);
h_track_jet_reco_matched_eta = new TH1F("track_jet_reco_matched_eta","track_jet_reco_matched_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_reco_matched_phi = new TH1F("track_jet_reco_matched_phi","track_jet_reco_matched_phi", phibins, phimin, phimax);
h_track_jet_reco_matched_nb  = new TH1F("track_jet_reco_matched_nb","track_jet_reco_matched_nb",19, -1, 18);
h_track_jet_reco_matched_E   = new TH1F("track_jet_reco_matched_E","track_jet_reco_matched_E",100, Emin, Emax);

h_track_jet_truth_matched_pt  = new TH1F("track_jet_truth_matched_pt","track_jet_truth_matched_pt", ptbins, ptmin, ptmax);
h_track_jet_truth_matched_eta = new TH1F("track_jet_truth_matched_eta","track_jet_truth_matched_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_truth_matched_phi = new TH1F("track_jet_truth_matched_phi","track_jet_truth_matched_phi", phibins, phimin, phimax);
h_track_jet_truth_matched_E   = new TH1F("track_jet_truth_matched_E","track_jet_truth_matched_E",100, Emin, Emax);

h_track_jet_reco_fake_pt  = new TH1F("track_jet_reco_fake_pt","track_jet_reco_fake_pt", ptbins, ptmin, ptmax);
h_track_jet_reco_fake_eta = new TH1F("track_jet_reco_fake_eta","track_jet_reco_fake_eta", etabins_central, etamin_central, etamax_central);
h_track_jet_reco_fake_phi = new TH1F("track_jet_reco_fake_phi","track_jet_reco_fake_phi", phibins, phimin, phimax);
h_track_jet_reco_fake_nb  = new TH1F("track_jet_reco_fake_nb","track_jet_reco_fake_nb",19, -1, 18);
        
h_track_jet_reco_fake_vs_eta = new TH2F("track_jet_reco_fake_vs_eta","Faked Track Jets vs eta", 65, -3.2, 3.2, 2, -0.5, 1.5);
h_track_jet_reco_fake_vs_pt  = new TH2F("track_jet_reco_fake_vs_pt", "Faked Track Jets vs pt", ptbins, ptmin, ptmax, 2, -0.5, 1.5);     

h_track_jet_truth_notmatched_pt  = new TH1F("track_jet_truth_notmatched_pt","track_jet_truth_notmatched_pt",  ptbins, ptmin, ptmax);
h_track_jet_truth_notmatched_eta = new TH1F("track_jet_truth_notmatched_eta","track_jet_truth_notmatched_eta",  etabins_central, etamin_central, etamax_central);
h_track_jet_truth_notmatched_phi = new TH1F("track_jet_truth_notmatched_phi","track_jet_truth_notmatched_phi",  phibins, phimin, phimax);
h_track_jet_truth_notmatched_nb  = new TH1F("track_jet_truth_notmatched_nb","track_jet_truth_notmatched_nb",19, -1, 18);
     
h_track_jet_matched_dpt  = new TH1F("track_jet_matched_dpt","track_jet_matched_dpt", 200, -d_ptmax, d_ptmax);
h_track_jet_matched_deta = new TH1F("track_jet_matched_deta","track_jet_matched_deta", 200, -d_etamax, d_etamax);
h_track_jet_matched_dphi = new TH1F("track_jet_matched_dphi","track_jet_matched_dphi", 200, -d_phimax, d_phimax);
h_track_jet_matched_dr   = new TH1F("track_jet_matched_dr","track_jet_matched_dr", 200, 0.0, 0.1);
h_track_jet_matched_dE   = new TH1F("track_jet_matched_dE","track_jet_matched_dE", 200, -d_Emax, d_Emax);

h_track_jet_matched_resolution_pt  = new TH2F("track_jet_matched_resolution_pt","track_jet_matched_resolution_pt",
                                             20, ptmin, ptmax, 100, -1.0, 1.0);
h_track_jet_matched_resolution_eta = new TH2F("track_jet_matched_resolution_eta","track_jet_matched_resolution_eta",
                                              20, etamin_central, etamax_central, 100, -1.0, 1.0);
h_track_jet_matched_resolution_phi = new TH2F("track_jet_matched_resolution_phi","track_jet_matched_resolution_phi",
                                              20, phimin, phimax, 100, -1.0, 1.0);

//=========
   
h_track_jet_matched_Eresolution = new TH1F("track_jet_matched_Eresolution", 
                                           "E RESOLUTION MATCHED TRACK_JETS (GeV)", 
                                           80, -200, 200);
   
h_track_jet_matched_Eresolution_vs_pt = new TH2F("track_jet_matched_Eresolution_vs_pt", 
                                                 "E RES MATCHED TRACK_JETS (GeV) VS PT (GeV)", 
                                                 ptbins, ptmin, ptmax, 80, -200, 200);
   
h_track_jet_matched_Eresolution_vs_eta = new TH2F("track_jet_matched_Eresolution_vs_eta", 
                                                  "E RES MATCHED TRACK_JETS (GeV) VS ETA", 
                                                  20, etamin_central, etamax_central, 80, -200, 200);
   
h_track_jet_matched_Eresolution_vs_phi = new TH2F("track_jet_matched_Eresolution_vs_phi", 
                                                  "E RES MATCHED TRACK_JETS (GeV) VS PHI", 
                                                  14, phimin, phimax, 80, -200, 200);
   
h_track_jet_matched_Eresolution_vs_E = new TH2F("track_jet_matched_Eresolution_vs_E", 
                                                "E RES MATCHED TRACK_JETS (GeV) VS E (GeV)", 
                                                Ebins, Emin, Emax, 80, -200, 200);
//=========
   
h_track_jet_matched_relEresolution = new TH1F("track_jet_matched_relEresolution", 
                                              "Rel E RESOLUTION MATCHED TRACK_JETS (GeV)", 
                                              80, -1, 1);
   
h_track_jet_matched_relEresolution_vs_pt = new TH2F("track_jet_matched_relEresolution_vs_pt", 
                                                    "Rel E RES MATCHED TRACK_JETS VS PT (GeV)", 
                                                    ptbins, ptmin, ptmax, 80, -1, 1);
   
h_track_jet_matched_relEresolution_vs_eta = new TH2F("track_jet_matched_relEresolution_vs_eta", 
                                                     "Rel E RES MATCHED TRACK_JETS VS ETA", 
                                                     20, etamin_central, etamax_central, 80, -1, 1);
   
h_track_jet_matched_relEresolution_vs_phi = new TH2F("track_jet_matched_relEresolution_vs_phi", 
                                                     "Rel E RES MATCHED TRACK_JETS VS PHI", 
                                                     14, phimin, phimax, 80, -1, 1);

h_track_jet_matched_relEresolution_vs_E = new TH2F("track_jet_matched_relEresolution_vs_E", 
                                                   "Rel E RES MATCHED TRACK_JETS VS E (GeV)", 
                                                   Ebins, Emin, Emax, 80, -1, 1);

//=========
h_track_jet_matched_invptresolution = new TH1F("track_jet_matched_invptresolution", 
                                               "track_jet_matched_invptresolution",
                                               80, -0.02, 0.02);
h_track_jet_matched_relinvptresolution = new TH1F("track_jet_matched_relinvptresolution", 
                                                  "track_jet_matched_relinvptresolution",
                                                  80, -1, 1);
h_track_jet_matched_invEresolution = new TH1F("track_jet_matched_invEresolution", 
                                              "track_jet_matched_invEresolution",
                                              80, -0.02, 0.02);
h_track_jet_matched_relinvEresolution = new TH1F("track_jet_matched_relinvEresolution", 
                                                 "track_jet_matched_relinvEresolution",
                                                 80, -1, 1);
h_track_jet_matched_invptresolution_vs_pt = new TH2F("track_jet_matched_invptresolution_vs_pt", 
                                                     "track_jet_matched_invptresolution_vs_pt", 
                                                     ptbins, ptmin, ptmax, 80, -0.02, 0.02);
h_track_jet_matched_relinvptresolution_vs_pt = new TH2F("track_jet_matched_relinvptresolution_vs_pt", 
                                                        "track_jet_matched_relinvptresolution_vs_pt", 
                                                        ptbins, ptmin, ptmax, 80, -1, 1);
h_track_jet_matched_invEresolution_vs_E = new TH2F("track_jet_matched_invEresolution_vs_E", 
                                                   "track_jet_matched_invEresolution_vs_E", 
                                                   Ebins, Emin, Emax, 80, -0.02, 0.02);
h_track_jet_matched_relinvEresolution_vs_E = new TH2F("track_jet_matched_relinvEresolution_vs_E", 
                                                      "track_jet_matched_relinvEresolution_vs_E", 
                                                      Ebins, Emin, Emax, 80, -1, 1);

//=========


//track jets with standard jets

h_track_std_jet_matched_pt  = new TH1F("track_std_jet_matched_pt","track_std_jet_matched_pt", ptbins, ptmin, ptmax);
h_track_std_jet_matched_eta = new TH1F("track_std_jet_matched_eta","track_std_jet_matched_eta", etabins_central, etamin_central, etamax_central);
h_track_std_jet_matched_phi = new TH1F("track_std_jet_matched_phi","track_std_jet_matched_phi", phibins, phimin, phimax);
h_track_std_jet_matched_E   = new TH1F("track_std_jet_matched_E","track_std_jet_matched_E",100, Emin, Emax);

h_track_std_jet_notmatched_pt  = new TH1F("track_std_jet_notmatched_pt","track_std_jet_notmatched_pt",  ptbins, ptmin, ptmax);
h_track_std_jet_notmatched_eta = new TH1F("track_std_jet_notmatched_eta","track_std_jet_notmatched_eta",  etabins_central, etamin_central, etamax_central);
h_track_std_jet_notmatched_phi = new TH1F("track_std_jet_notmatched_phi","track_std_jet_notmatched_phi",  phibins, phimin, phimax);
h_track_std_jet_notmatched_E  = new TH1F("track_std_jet_notmatched_E","track_std_jet_notmatched_E",100, Emin, Emax);

h_track_std_jet_matched_dpt  = new TH1F("track_std_jet_matched_dpt","track_std_jet_matched_dpt", 200, -d_ptmax, d_ptmax);
h_track_std_jet_matched_deta = new TH1F("track_std_jet_matched_deta","track_std_jet_matched_deta", 200, -d_etamax, d_etamax);
h_track_std_jet_matched_dphi = new TH1F("track_std_jet_matched_dphi","track_std_jet_matched_dphi", 200, -d_phimax, d_phimax);
h_track_std_jet_matched_dr   = new TH1F("track_std_jet_matched_dr","track_std_jet_matched_dr", 200, 0.0, 0.1);
h_track_std_jet_matched_dE   = new TH1F("track_std_jet_matched_dE","track_std_jet_matched_dE", 200, -d_Emax, d_Emax);

h_track_std_jet_matched_resolution_pt  = new TH2F("track_std_jet_matched_resolution_pt","track_std_jet_matched_resolution_pt",
                                             20, ptmin, ptmax, 100, -1.0, 1.0);
h_track_std_jet_matched_resolution_eta = new TH2F("track_std_jet_matched_resolution_eta","track_std_jet_matched_resolution_eta",
                                              20, etamin_central, etamax_central, 100, -1.0, 1.0);
h_track_std_jet_matched_resolution_phi = new TH2F("track_std_jet_matched_resolution_phi","track_std_jet_matched_resolution_phi",
                                              20, phimin, phimax, 100, -1.0, 1.0);
        
h_track_std_jet_matched_vs_eta = new TH2F("track_std_jet_matched_vs_eta","track_std_jet_matched_vs_eta", 65, -3.2, 3.2, 2, -0.5, 1.5);
h_track_std_jet_matched_vs_pt  = new TH2F("track_std_jet_matched_vs_pt", "track_std_jet_matched_vs_pt", ptbins, ptmin, ptmax, 2, -0.5, 1.5);    

h_track_std_jet_matched_jetptratio = new TH1F("track_std_jet_matched_jetptratio", 
                                              "track_std_jet_matched_jetptratio", 
                                              60, 0, 2);
h_track_std_jet_matched_jetptratio_corr = new TH1F("track_std_jet_matched_jetptratio_corr", 
                                                   "track_std_jet_matched_jetptratio_corr", 
                                                   60, 0, 2);
//h_track_std_jet_matched_2Dpt = new TH2F("track_std_jet_matched_pt2D", "track_std_jet_matched_2Dpt", 
//                                  ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);        

//h_track_std_jet_matched_pt_ratio_vs_pttrk = new TH2F("track_std_jet_matched_pt_ratio_vs_pttrk", "track_std_jet_matched_pt_ratio_vs_pttrk", 
//                                                ptbins, ptmin, ptmax, 60, 0, 2);

//h_track_std_jet_matched_pt_ratio_vs_ptstd = new TH2F("track_std_jet_matched_pt_ratio_vs_ptstd", "track_std_jet_matched_pt_ratio_vs_ptstd", 
//                                                ptbins, ptmin, ptmax, 60, 0, 2);


h_track_std_jet_matched_pt_vs_jetptratio= new TH2F("track_std_jet_matched_pt_vs_jetptratio", 
                                                   "track_std_jet_matched_pt_vs_jetptratio",
                                                   ptbins, ptmin, ptmax, 60, 0, 2);
h_track_std_jet_matched_pt_vs_jetptratio_corr= new TH2F("track_std_jet_matched_pt_vs_jetptratio_corr", 
                                                        "track_std_jet_matched_pt_vs_jetptratio_corr",
                                                        ptbins, ptmin, ptmax, 60, 0, 2);
h_track_std_jet_matched_pt_vs_stdjetpt= new TH2F("track_std_jet_matched_pt_vs_stdjetpt", 
                                                 "track_std_jet_matched_pt_vs_stdjetpt",
                                                 ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);
h_track_std_jet_matched_pt_vs_stdjetpt_corr= new TH2F("track_std_jet_matched_pt_vs_stdjetpt_corr", 
                                                      "track_std_jet_matched_pt_vs_stdjetpt_corr",
                                                      ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);


h_std_track_jet_matched_pt  = new TH1F("std_track_jet_matched_pt","std_track_jet_matched_pt", ptbins, ptmin, ptmax);
h_std_track_jet_matched_eta = new TH1F("std_track_jet_matched_eta","std_track_jet_matched_eta", etabins_central, etamin_central, etamax_central);
h_std_track_jet_matched_phi = new TH1F("std_track_jet_matched_phi","std_track_jet_matched_phi", phibins, phimin, phimax);
h_std_track_jet_matched_E   = new TH1F("std_track_jet_matched_E","std_track_jet_matched_E",100, Emin, Emax);

h_std_track_jet_notmatched_pt    = new TH1F("std_track_jet_notmatched_pt","std_track_jet_notmatched_pt",  ptbins, ptmin, ptmax);
h_std_track_jet_notmatched_eta   = new TH1F("std_track_jet_notmatched_eta","std_track_jet_notmatched_eta",  etabins_central, etamin_central, etamax_central);
h_std_track_jet_notmatched_phi   = new TH1F("std_track_jet_notmatched_phi","std_track_jet_notmatched_phi",  phibins, phimin, phimax);
h_std_track_jet_notmatched_E     = new TH1F("std_track_jet_notmatched_E","std_track_jet_notmatched_E",100, Emin, Emax);
h_std_track_jet_notmatched_ntrks = new TH1F("std_track_jet_notmatched_ntrks","std_track_jet_notmatched_trks", 20, 0, 20);

h_std_track_jet_matched_dpt  = new TH1F("std_track_jet_matched_dpt","std_track_jet_matched_dpt", 200, -d_ptmax, d_ptmax);
h_std_track_jet_matched_deta = new TH1F("std_track_jet_matched_deta","std_track_jet_matched_deta", 200, -d_etamax, d_etamax);
h_std_track_jet_matched_dphi = new TH1F("std_track_jet_matched_dphi","std_track_jet_matched_dphi", 200, -d_phimax, d_phimax);
h_std_track_jet_matched_dr   = new TH1F("std_track_jet_matched_dr","std_track_jet_matched_dr", 200, 0.0, 0.1);
h_std_track_jet_matched_dE   = new TH1F("std_track_jet_matched_dE","std_track_jet_matched_dE", 200, -d_Emax, d_Emax);

h_std_track_jet_matched_resolution_pt  = new TH2F("std_track_jet_matched_resolution_pt","std_track_jet_matched_resolution_pt",
                                             20, ptmin, ptmax, 100, -1.0, 1.0);
h_std_track_jet_matched_resolution_eta = new TH2F("std_track_jet_matched_resolution_eta","std_track_jet_matched_resolution_eta",
                                              20, etamin_central, etamax_central, 100, -1.0, 1.0);
h_std_track_jet_matched_resolution_phi = new TH2F("std_track_jet_matched_resolution_phi","std_track_jet_matched_resolution_phi",
                                              20, phimin, phimax, 100, -1.0, 1.0);
        
h_std_track_jet_matched_vs_eta = new TH2F("std_track_jet_matched_vs_eta","std_track_jet_matched_vs_eta", 65, -3.2, 3.2, 2, -0.5, 1.5);
h_std_track_jet_matched_vs_pt  = new TH2F("std_track_jet_matched_vs_pt", "std_track_jet_matched_vs_pt", ptbins, ptmin, ptmax, 2, -0.5, 1.5);    



h_std_track_jet_matched_pt_corr = new TH1F("std_track_jet_matched_pt_corr",
                                           "std_track_jet_matched_pt_corr", 
                                           ptbins, ptmin, ptmax);
h_std_track_jet_matched_trkptratio= new TH1F("std_track_jet_matched_trkptratio",
                                             "std_track_jet_matched_trkptratio",
                                             100, 0,2);
h_std_track_jet_matched_cmstrkptratio= new TH1F("std_track_jet_matched_cmstrkptratio",
                                                "std_track_jet_matched_cmstrkptratio",
                                                100, 0,2);
h_std_track_jet_matched_pt_vs_jetptratio= new TH2F("std_track_jet_matched_pt_vs_jetptratio", 
                                                   "std_track_jet_matched_pt_vs_jetptratio",
                                                   ptbins, ptmin, ptmax, 60, 0, 2);
h_std_track_jet_matched_pt_vs_jetptratio_corr= new TH2F("std_track_jet_matched_pt_vs_jetptratio_corr", 
                                                        "std_track_jet_matched_pt_vs_jetptratio_corr",
                                                        ptbins, ptmin, ptmax, 60, 0, 2);
h_std_track_jet_matched_pt_vs_trackjetpt= new TH2F("std_track_jet_matched_pt_vs_trackjetpt", 
                                                   "std_track_jet_matched_pt_vs_trackjetpt",
                                                   ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);
h_std_track_jet_matched_pt_vs_trackjetpt_corr= new TH2F("std_track_jet_matched_pt_vs_trackjetpt_corr", 
                                                        "std_track_jet_matched_pt_vs_trackjetpt_corr",
                                                        ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);



//truth jet matching

h_truth_track_std_jet_matched_pt  = new TH1F("truth_track_std_jet_matched_pt","truth_track_std_jet_matched_pt", ptbins, ptmin, ptmax);
h_truth_track_std_jet_matched_eta = new TH1F("truth_track_std_jet_matched_eta","truth_track_std_jet_matched_eta", etabins_central, etamin_central, etamax_central);
h_truth_track_std_jet_matched_phi = new TH1F("truth_track_std_jet_matched_phi","truth_track_std_jet_matched_phi", phibins, phimin, phimax);
h_truth_track_std_jet_matched_E   = new TH1F("truth_track_std_jet_matched_E","truth_track_std_jet_matched_E",100, Emin, Emax);

h_truth_track_std_jet_notmatched_pt  = new TH1F("truth_track_std_jet_notmatched_pt","truth_track_std_jet_notmatched_pt",  ptbins, ptmin, ptmax);
h_truth_track_std_jet_notmatched_eta = new TH1F("truth_track_std_jet_notmatched_eta","truth_track_std_jet_notmatched_eta",  etabins_central, etamin_central, etamax_central);
h_truth_track_std_jet_notmatched_phi = new TH1F("truth_track_std_jet_notmatched_phi","truth_track_std_jet_notmatched_phi",  phibins, phimin, phimax);
h_truth_track_std_jet_notmatched_E  = new TH1F("truth_track_std_jet_notmatched_E","truth_track_std_jet_notmatched_E",100, Emin, Emax);


h_truth_track_std_jet_matched_vs_eta = new TH2F("truth_track_std_jet_matched_vs_eta","truth_track_std_jet_matched_vs_eta", 65, -3.2, 3.2, 2, -0.5, 1.5);
h_truth_track_std_jet_matched_vs_pt  = new TH2F("truth_track_std_jet_matched_vs_pt", "truth_track_std_jet_matched_vs_pt", ptbins, ptmin, ptmax, 2, -0.5, 1.5);  

h_truth_track_std_jet_matched_jetptratio = new TH1F("truth_track_std_jet_matched_jetptratio", 
                                                    "truth_track_std_jet_matched_jetptratio", 
                                                    60, 0, 2);
//h_truth_track_std_jet_matched_2Dpt = new TH2F("truth_track_std_jet_matched_pt2D", "h_truth_track_std_jet_matched_2Dpt", 
//                                            ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);        
//h_truth_track_std_jet_matched_pt_ratio_vs_pttrk = new TH2F("truth_track_std_jet_matched_pt_ratio_vs_pttrk", "truth_track_std_jet_matched_pt_ratio_vs_pttrk", 
//                                                      ptbins, ptmin, ptmax, 60, 0, 2);
//h_truth_track_std_jet_matched_pt_ratio_vs_ptstd = new TH2F("truth_track_std_jet_matched_pt_ratio_vs_ptstd", "truth_track_std_jet_matched_pt_ratio_vs_ptstd", 
//                                                      ptbins, ptmin, ptmax, 60, 0, 2);

h_truth_track_std_jet_matched_pt_vs_stdjetpt = new TH2F("truth_track_std_jet_matched_pt_vs_stdjetpt", 
                                                        "truth_track_std_jet_matched_pt_vs_stdjetpt", 
                                                        ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);        
h_truth_track_std_jet_matched_pt_vs_jetptratio = new TH2F("truth_track_std_jet_matched_pt_vs_jetptratio",
                                                          "truth_track_std_jet_matched_pt_vs_jetptratio", 
                                                          ptbins, ptmin, ptmax, 60, 0, 2);
//h_truth_track_std_jet_matched_pt_ratio_vs_ptstd = new TH2F("truth_track_std_jet_matched_pt_ratio_vs_ptstd", "truth_track_std_jet_matched_pt_ratio_vs_ptstd", 
//                                                         ptbins, ptmin, ptmax, 60, 0, 2);


h_truth_std_track_jet_matched_pt  = new TH1F("truth_std_track_jet_matched_pt","truth_std_track_jet_matched_pt", ptbins, ptmin, ptmax);
h_truth_std_track_jet_matched_eta = new TH1F("truth_std_track_jet_matched_eta","truth_std_track_jet_matched_eta", etabins_central, etamin_central, etamax_central);
h_truth_std_track_jet_matched_phi = new TH1F("truth_std_track_jet_matched_phi","truth_std_track_jet_matched_phi", phibins, phimin, phimax);
h_truth_std_track_jet_matched_E   = new TH1F("truth_std_track_jet_matched_E","truth_std_track_jet_matched_E",100, Emin, Emax);

h_truth_std_track_jet_notmatched_pt  = new TH1F("truth_std_track_jet_notmatched_pt","truth_std_track_jet_notmatched_pt",  ptbins, ptmin, ptmax);
h_truth_std_track_jet_notmatched_eta = new TH1F("truth_std_track_jet_notmatched_eta","truth_std_track_jet_notmatched_eta",  etabins_central, etamin_central, etamax_central);
h_truth_std_track_jet_notmatched_phi = new TH1F("truth_std_track_jet_notmatched_phi","truth_std_track_jet_notmatched_phi",  phibins, phimin, phimax);
h_truth_std_track_jet_notmatched_E  = new TH1F("truth_std_track_jet_notmatched_E","truth_std_track_jet_notmatched_E",100, Emin, Emax);

h_truth_std_track_jet_matched_vs_eta = new TH2F("truth_std_track_jet_matched_vs_eta","truth_std_track_jet_matched_vs_eta", 65, -3.2, 3.2, 2, -0.5, 1.5);
h_truth_std_track_jet_matched_vs_pt  = new TH2F("truth_std_track_jet_matched_vs_pt", "truth_std_track_jet_matched_vs_pt", ptbins, ptmin, ptmax, 2, -0.5, 1.5);  


h_truth_std_track_jet_matched_jetptratio = new TH1F("truth_std_track_jet_matched_pt_ratio",
                                                    "truth_std_track_jet_matched_pt_ratio", 
                                                    60, 0, 2);
h_truth_std_track_jet_matched_pt_vs_trackjetpt = new TH2F("truth_std_track_jet_matched_pt_vs_trackjetpt", 
                                                          "truth_std_track_jet_matched_pt_vs_trackjetpt", 
                                                          ptbins, ptmin, ptmax, ptbins, ptmin, ptmax);        
h_truth_std_track_jet_matched_pt_vs_jetptratio = new TH2F("truth_std_track_jet_matched_pt_vs_jetptratio",
                                                          "truth_std_track_jet_matched_pt_vs_jetptratio", 
                                                          ptbins, ptmin, ptmax, 60, 0, 2);


 //  h_jet1_pt = new TH1F("jet1_pt", 
//                       "PT OF LEADING JET", 
//                       pt_histo_bins, -0.5,  pt_histo_hilim);

//     h_jet2_pt = new TH1F("jet2_pt", 
//                       "PT OF SUBLEADING JET", 
//                       pt_histo_bins, -0.5,  pt_histo_hilim);

//     h_jet3_pt = new TH1F("jet3_pt", 
//                       "PT OF THIRD JET", 
//                       pt_histo_bins, -0.5,  pt_histo_hilim);
    
//     h_tagjets_deltaeta = new TH1F("tagjets_deltaeta", 
//                                "DELTA ETA OF TAG JETS", 
//                                20, -0, 10);
    
//     h_tagjets_eta3star = new TH1F("tagjets_eta3star", 
//                                "ZEPPENFELD PLOT OF TAG JETS", 
//                                20, -5, 5);
    
//     //=========

//     h_jet_EtavsPt = new TH2F("jet_EtavsPt", 
//                           "ETA vs PT (GeV) OF JETS",
//                           50, -5, 5, 50, -0.5, 250.5);
   
//     h_jet_EtavsE  = new TH2F("jet_EtavsE",
//                           "ETA vs E (GeV) OF JETS",
//                           50, -5, 5, 50, -0.5, 250.5);



h_remove_double_dpt = new TH1F("remove_double_dpt","remove_double_dpt", 200, -d_ptmax, d_ptmax);
h_remove_double_dr = new TH1F("remove_double_dr","remove_double_dr", 200, 0.0, 0.5);
h_remove_double_nb_before = new TH1F("remove_double_nb_before","remove_double_nb_before",19, -1, 18);
h_remove_double_nb_after = new TH1F("remove_double_nb_after","remove_double_nb_after",19, -1, 18);
h_remove_double_nb_doubles = new TH1F("remove_double_nb_doubles","remove_double_nb_doubles",19, -1, 18);


// histos for the analysis

   char char_help_array[50];    
   
   for (int j=0; j<nb_cuts_max+2; j++){
        string nb;
        sprintf(char_help_array,"%02d",j);
        nb = char_help_array;
        string s1  = "";
        
        s1 = "etmiss_ana_"+nb;
        v_h_etmiss_ana.push_back(new TH1F(s1.c_str(),s1.c_str(), ptbins, ptmin, ptmax));
        
        s1 = "etmiss_x_ana_"+nb;
        v_h_etmiss_x_ana.push_back(new TH1F(s1.c_str(),s1.c_str(), ptbins, ptmin, ptmax));

        s1 = "etmiss_y_ana_"+nb;
        v_h_etmiss_y_ana.push_back(new TH1F(s1.c_str(),s1.c_str(), ptbins, ptmin, ptmax));

        s1 = "etmiss_ana_truth_"+nb;
        v_h_etmiss_ana_truth.push_back(new TH1F(s1.c_str(),s1.c_str(), ptbins, ptmin, ptmax));
        
        s1 = "etmiss_ana_reco_truth_dpt_"+nb;
        v_h_etmiss_ana_reco_truth_dpt.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.2*ptmax, 0.2*ptmax));
        
        s1 = "etmiss_ana_reco_truth_dx_"+nb;
        v_h_etmiss_ana_reco_truth_dx.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.2*ptmax, 0.2*ptmax));
        
        s1 = "etmiss_ana_reco_truth_dy_"+nb;
        v_h_etmiss_ana_reco_truth_dy.push_back(new TH1F(s1.c_str(),s1.c_str(),100, -0.2*ptmax, 0.2*ptmax));
        
        s1 = "e_ana_nb_"+nb;
        v_h_e_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
          
        s1 = "e_ana_pt_"+nb;
        v_h_e_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
          
        s1 = "e_ana_eta_"+nb;
        v_h_e_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
          
        s1 = "e_ana_phi_"+nb;
        v_h_e_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
          
        s1 = "mu_ana_nb_"+nb;
        v_h_mu_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
          
        s1 = "mu_ana_pt_"+nb;
        v_h_mu_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
          
        s1 = "mu_ana_eta_"+nb;
        v_h_mu_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
          
        s1 = "mu_ana_phi_"+nb;
        v_h_mu_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
          
        s1 = "mu_ana_fitChi2_"+nb;
        v_h_mu_ana_fitChi2.push_back(new TH1F(s1.c_str(),s1.c_str() , 100, -1, 50));
       
        s1 = "mu_ana_fitChi2OverDoF_"+nb;
        v_h_mu_ana_fitChi2OverDoF.push_back(new TH1F(s1.c_str(),s1.c_str() , 100, -1, 9));
        
        s1 = "mu_ana_matchChi2_"+nb;
        v_h_mu_ana_matchChi2.push_back(new TH1F(s1.c_str(),s1.c_str() , 100, -1, 50));
        
        s1 = "mu_ana_matchChi2OverDoF_"+nb;
        v_h_mu_ana_matchChi2OverDoF.push_back(new TH1F(s1.c_str(),s1.c_str() , 100, -1, 9));
       
        s1 = "lepton_ana_nb_"+nb;
        v_h_lepton_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
          
        s1 = "lepton_ana_pt_"+nb;
        v_h_lepton_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
          
        s1 = "lepton_ana_eta_"+nb;
        v_h_lepton_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
          
        s1 = "lepton_ana_phi_"+nb;
        v_h_lepton_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
          
        s1 = "lpair_ana_nb_"+nb;
        v_h_2l_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str(),19, -1, 18));
        
        s1 = "lpair_ana_eta_"+nb;
        v_h_2l_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
        
        s1 = "lpair_ana_phi_"+nb;
        v_h_2l_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
        
        s1 = "lpair_ana_pt_"+nb;
        v_h_2l_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
        
        s1 = "lpair_ana_costheta_"+nb;
        v_h_2l_ana_costheta.push_back(new TH1F(s1.c_str(),s1.c_str() , 80, -1., 1.));

        s1 = "lpair_ana_deta_"+nb;
        v_h_2l_ana_deta.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 10.));
        
        s1 = "lpair_ana_dphi_"+nb;
        v_h_2l_ana_dphi.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 4.));
        
        s1 = "lpair_ana_dpt_"+nb;
        v_h_2l_ana_dpt.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0, 100.));
        
        s1 = "lpair_ana_dr_"+nb;
        v_h_2l_ana_dr.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 8.));
        
        s1 = "lpair_ana_m_"+nb;
        v_h_2l_ana_m.push_back(new TH1F(s1.c_str(),s1.c_str() ,100, 0., 1000.));
        
        s1 = "jet_ana_nb_"+nb;
        v_h_jet_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
          
        s1 = "jet_ana_pt_"+nb;
        v_h_jet_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
          
        s1 = "jet_ana_eta_"+nb;
        v_h_jet_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
          
        s1 = "jet_ana_phi_"+nb;
        v_h_jet_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
        
        s1 = "jet_ana_bweight_"+nb;
        v_h_jet_ana_bweight.push_back(new TH1F(s1.c_str(),s1.c_str() , 30, -15, 15));

        s1 = "jpair_ana_nb_"+nb;
        v_h_2j_ana_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
        
        s1 = "jpair_ana_eta_"+nb;
        v_h_2j_ana_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
        
        s1 = "jpair_ana_phi_"+nb;
        v_h_2j_ana_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));
        
        s1 = "jpair_ana_pt_"+nb;
        v_h_2j_ana_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
        
        s1 = "jpair_ana_deta_"+nb;
        v_h_2j_ana_deta.push_back(new TH1F(s1.c_str(),s1.c_str() , 60, 0., 10.));
        
        s1 = "jpair_ana_dphi_"+nb;
        v_h_2j_ana_dphi.push_back(new TH1F(s1.c_str(),s1.c_str() , 60, 0., 4.));
        
        s1 = "jpair_ana_dpt_"+nb;
        v_h_2j_ana_dpt.push_back(new TH1F(s1.c_str(),s1.c_str() , 80, 0., 100.));
        
        s1 = "jpair_ana_dr_"+nb;
        v_h_2j_ana_dr.push_back(new TH1F(s1.c_str(),s1.c_str() ,60, 0.0, 8.));
        
        s1 = "jpair_ana_m_"+nb;
        v_h_2j_ana_m.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 1500.));

        s1 = "ana_pttot_"+nb;
        v_h_ana_pttot.push_back(new TH1F(s1.c_str(),s1.c_str(), ptbins, ptmin, 2*ptmax));

 //        s1 = "ana_pt_vis_"+nb;
//      v_h_ana_pt_vis.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 2000.));

//         s1 = "ana_pt_miss_"+nb;
//      v_h_ana_pt_miss.push_back(new TH1F(s1.c_str(),s1.c_str() ,80, 0., 2000.));

        s1 = "candidates_nb_"+nb;
        v_h_candidates_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,30, 0., 30.));


        s1 = "lpair_x_tau1_"+nb;
        v_h_2l_x_tau1.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
        
        s1 = "lpair_x_tau2_"+nb;
        v_h_2l_x_tau2.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
        
        s1 = "lpair_x_tau12_"+nb;
        v_h_2l_x_tau12.push_back(new TH1F(s1.c_str(),s1.c_str() ,100,-3.,3.));
        
        s1 = "lpair_x_tau1vs2_"+nb;
        v_h_2l_x_tau1vs2.push_back(new TH2F(s1.c_str(),s1.c_str() ,100,-3.,3., 100,-3.,3.));

        s1 = "lpair_tautau_Z_dm_"+nb;
        v_h_2l_tautau_Z_dm.push_back(new TH1F(s1.c_str(),s1.c_str() ,160,0.,100.));

        s1 = "ana_mt_llnu_"+nb;
        v_h_ana_mt_llnu.push_back(new TH1F(s1.c_str(),s1.c_str(), 40, ptmin, ptmax));

        s1 = "ana_mt_"+nb;
        v_h_ana_mt.push_back(new TH1F(s1.c_str(),s1.c_str(), 40, ptmin, ptmax));
        
        //all central jet histos

        s1 = "centraljet_all_nb_"+nb;
        v_h_centraljet_all_nb.push_back(new TH1F(s1.c_str(),s1.c_str() ,19, -1, 18));
          
        s1 = "centraljet_all_pt_"+nb;
        v_h_centraljet_all_pt.push_back(new TH1F(s1.c_str(),s1.c_str() , ptbins, ptmin, ptmax));
          
        s1 = "centraljet_all_eta_"+nb;
        v_h_centraljet_all_eta.push_back(new TH1F(s1.c_str(),s1.c_str() , etabins, etamin, etamax));
          
        s1 = "centraljet_all_phi_"+nb;
        v_h_centraljet_all_phi.push_back(new TH1F(s1.c_str(),s1.c_str() , phibins, phimin, phimax));


//      s1 = ""+nb;
//      v_h_.push_back(new TH1F(s1.c_str(),s1.c_str() ,));
   }



        h_e_eff    = new TH1F("e_eff","electron efficiency", 2, -0.5, 1.5);
        h_e_fake   = new TH1F("e_fake","electron fakerate", 2, -0.5, 1.5);
        h_mu_eff   = new TH1F("mu_eff","muon efficiency", 2, -0.5, 1.5);
        h_mu_fake  = new TH1F("mu_fake","muon fakerate", 2, -0.5, 1.5);
        h_jet_eff  = new TH1F("jet_eff","jet efficiency", 2, -0.5, 1.5);
        h_jet_fake = new TH1F("jet_fake","jet fakerate", 2, -0.5, 1.5);

        h_track_jet_eff  = new TH1F("track_jet_eff","track_jet efficiency", 2, -0.5, 1.5);
        h_track_jet_fake = new TH1F("track_jet_fake","track_jet fakerate", 2, -0.5, 1.5);

        h_track_std_jet_match  = new TH1F("track_std_jet_match","track_std_jet_match", 2, -0.5, 1.5);
        h_std_track_jet_match  = new TH1F("std_track_jet_match","std_track_jet_match", 2, -0.5, 1.5);
        
        h_truth_track_std_jet_match  = new TH1F("truth_track_std_jet_match","truth_track_std_jet_match", 2, -0.5, 1.5);
        h_truth_std_track_jet_match  = new TH1F("truth_std_track_jet_match","truth_std_track_jet_match", 2, -0.5, 1.5);

        h_b_eff = new TH1F("b_eff","btag efficiency", 2, -0.5, 1.5);
        h_b_fake  = new TH1F("b_fake","btag fakerate", 2, -0.5, 1.5);
        h_bjet_ratio = new TH1F("bjet_ratio","bjet ratio", 2, -0.5, 1.5);

        h_b_truth_pt = new TH1F("b_truth_pt","b_truth_pt", ptbins, ptmin, ptmax);
//h_e_truth_presel_eta = new TH1F("e_truth_presel_eta","e_truth_presel_eta", etabins, etamin, etamax);
//h_e_truth_presel_phi = new TH1F("e_truth_presel_phi","e_truth_presel_phi", phibins, phimin, phimax);
        h_b_truth_nb = new TH1F("b_truth_nb","b_truth_nb",19, -1, 18);
        h_b_truth_mother_pdgID = new TH1F("b_truth_mother_pdgID","b_truth_mother_pdgID",3001, -1, 3000);
        h_b_truth_status = new TH1F("b_truth_status","b_truth_status",1001, -1, 1001);

        h_fwj_cj_dr = new TH1F("fwj_cj_dr","Delta R between central jet and forward jets", 200, 0.0, 6.);
        h_fwj_cj_dr_nocut = new TH1F("fwj_cj_dr_nocut","Delta R between central jet and forward jets (nocut)", 200, 0.0, 6.);

       

// central jets
        h_centraljet_ana_nb = new TH1F("centraljet_ana_nb","centraljet_ana_nb",19, -1, 18);
        h_centraljet_ana_pt = new TH1F("centraljet_ana_pt","centraljet_ana_pt",  ptbins, ptmin, ptmax);
        h_centraljet_ana_eta = new TH1F("centraljet_ana_eta","centraljet_ana_eta",  etabins, etamin, etamax);
        h_centraljet_ana_phi = new TH1F("centraljet_ana_phi","centraljet_ana_phi",  phibins, phimin, phimax);
        h_centraljet_ana_zeppvar = new TH1F("centraljet_ana_zeppvar","centraljet_ana_zeppvar",  etabins, etamin, etamax);
        h_centraljet_ana_ntrks   = new TH1F("centraljet_ana_ntrks","number of tracks per central jet", 20, 0.0, 20.);

        h_centraljet_cand_nb = new TH1F("centraljet_cand_nb","centraljet_cand_nb",19, -1, 18);
        h_centraljet_cand_pt = new TH1F("centraljet_cand_pt","centraljet_cand_pt",  ptbins, ptmin, ptmax);
        h_centraljet_cand_eta = new TH1F("centraljet_cand_eta","centraljet_cand_eta",  etabins, etamin, etamax);
        h_centraljet_cand_phi = new TH1F("centraljet_cand_phi","centraljet_cand_phi",  phibins, phimin, phimax);
 
        h_centraljet_std_cand_nb = new TH1F("centraljet_std_cand_nb","centraljet_std_cand_nb",19, -1, 18);
        h_centraljet_std_cand_pt = new TH1F("centraljet_std_cand_pt","centraljet_std_cand_pt",  ptbins, ptmin, ptmax);
        h_centraljet_std_cand_eta = new TH1F("centraljet_std_cand_eta","centraljet_std_cand_eta",  etabins, etamin, etamax);
        h_centraljet_std_cand_phi = new TH1F("centraljet_std_cand_phi","centraljet_std_cand_phi",  phibins, phimin, phimax);
        h_centraljet_std_cand_ntrks   = new TH1F("centraljet_std_cand_ntrks","number of tracks per central jet", 20, 0.0, 20.);

        h_centraljet_trk_cand_nb = new TH1F("centraljet_trk_cand_nb","centraljet_trk_cand_nb",19, -1, 18);
        h_centraljet_trk_cand_pt = new TH1F("centraljet_trk_cand_pt","centraljet_trk_cand_pt",  ptbins, ptmin, ptmax);
        h_centraljet_trk_cand_eta = new TH1F("centraljet_trk_cand_eta","centraljet_trk_cand_eta",  etabins, etamin, etamax);
        h_centraljet_trk_cand_phi = new TH1F("centraljet_trk_cand_phi","centraljet_trk_cand_phi",  phibins, phimin, phimax);
        
      
        h_jl_mineta = new TH1F("jl_mineta","jl_mineta",50,-5.,5.);
        h_jl_maxeta = new TH1F("jl_maxeta","jl_maxeta",50,-5.,5.);

// h_higgs_mass = new TH1F("higgs_mass","higgs_mass",40,0.,300.);
// h_lmissinget_ana_mt = new TH1F("lmissinget_ana_mt","lmissinget_ana_mt", ptbins, ptmin, ptmax);
// h_lmissinget_ana_mt2 = new TH1F("lmissinget_ana_mt2","lmissinget_ana_mt2", ptbins, ptmin, ptmax);

        h_trigger_online = new TH1F("trigger_online","trigger_online",19, -1, 18);
        h_trigger_offline = new TH1F("trigger_offline","trigger_offline",19, -1, 18);
        h_trigger_combined = new TH1F("trigger_combined","trigger_combined",19, -1, 18);

---