#ifndef histoLeptonId_h
#define histoLeptonId_h

#include "TH1.h"
#include "TH2.h"
#include "TF1.h"

#include <TPostScript.h>
#include <TCanvas.h>

#include <TStyle.h>
#include <TText.h>
#include <TFile.h>
#include <TSystem.h>
#include <TPaveLabel.h>


#include <stdlib.h>
#include <iostream.h>
#include <string.h>


class histoLeptonId{
 public:

  histoLeptonId(){};
  virtual ~histoLeptonId(){};

  struct Hist_t{
    
    TH1F *hele_delR;
    TH1F *hmu_delR;
    TH1F *hele0_delR;
    TH1F *hmu0_delR;
    TH1F *he_n;
    TH1F *he_corrtotalet;
    TH1F *he_rawtotalet;
    TH1F *he_rawtotale;
    TH1F *hnelec;
    TH1F *hnmuon;

    TH1F *hbe_et;
    TH1F *hbe_pt;
    TH1F *hbe_eop;
    TH1F *hbe_z0;
    TH1F *hbe_hadoem;
    TH1F *hbe_lshr;
    TH1F *hbe_csdeltax;
    TH1F *hbe_deltaz;
    TH1F *hbe_chi2;
    TH1F *hbe_fcaliso;
    TH1F *hbe_trkphi;
    TH1F *hbe_trketa;

    TH1F *hbec_et;
    TH1F *hbec_pt;
    TH1F *hbec_eop;
    TH1F *hbec_z0;
    TH1F *hbec_hadoem;
    TH1F *hbec_lshr;
    TH1F *hbec_csdeltax;
    TH1F *hbec_deltaz;
    TH1F *hbec_chi2;
    TH1F *hbec_fcaliso;
    TH1F *hbec_trkphi;
    TH1F *hbec_trketa;

    TH1F *hbm_pt;
    TH1F *hbm_z0;
    TH1F *hbm_em;
    TH1F *hbm_had;
    TH1F *hbm_d0ws;
    TH1F *hbm_d0wos;
    TH1F *hbm_fciso;
    TH1F *hbm_deltax_cmu;
    TH1F *hbm_deltax_cmp;
    TH1F *hbm_deltax_cmx;
    TH1F *hbm_trketa_cmup;
    TH1F *hbm_trketa_cmu;
    TH1F *hbm_trketa_cmp;
    TH1F *hbm_trketa_cmx;
    TH1F *hbm_trkphi_cmup;
    TH1F *hbm_trkphi_cmu;
    TH1F *hbm_trkphi_cmp;
    TH1F *hbm_trkphi_cmx;

    TH1F *hbmc_pt;
    TH1F *hbmc_z0;
    TH1F *hbmc_em;
    TH1F *hbmc_had;
    TH1F *hbmc_d0ws;
    TH1F *hbmc_d0wos;
    TH1F *hbmc_fciso;
    TH1F *hbmc_deltax_cmu;
    TH1F *hbmc_deltax_cmp;
    TH1F *hbmc_deltax_cmx;
    TH1F *hbmc_trketa_cmup;
    TH1F *hbmc_trketa_cmu;
    TH1F *hbmc_trketa_cmp;
    TH1F *hbmc_trketa_cmx;
    TH1F *hbmc_trkphi_cmup;
    TH1F *hbmc_trkphi_cmu;
    TH1F *hbmc_trkphi_cmp;
    TH1F *hbmc_trkphi_cmx;

    TH1F *hdilepflag;

    TH1F *hjets_njets;
    TH1F *hjets_et1;
    TH1F *hjets_eta1;
    TH1F *hjets_phi1;
    TH1F *hjets_et2;
    TH1F *hjets_eta2;
    TH1F *hjets_phi2;
    TH1F *hjets_et3;
    TH1F *hjets_eta3;
    TH1F *hjets_phi3;

    TH1F *hmet_j0_ee;
    TH1F *hmet_j1_ee;
    TH1F *hmet_j2_ee;

    TH1F *hmet_j0_emu;
    TH1F *hmet_j1_emu;
    TH1F *hmet_j2_emu;

    TH1F *hmet_j0_mumu;
    TH1F *hmet_j1_mumu;
    TH1F *hmet_j2_mumu;

    TH1F *hmee;
    TH1F *hmmumu;
    TH1F *hmee0;
    TH1F *hmmumu0;

    TH2F *heeDphimljMetj;
    TH2F *heeDphimljMetj0;
    TH2F *heeDphimljMetj1;
    TH2F *heeDphimljMetj2;

    TH2F *hemuDphimljMetj;
    TH2F *hemuDphimljMetj0;
    TH2F *hemuDphimljMetj1;
    TH2F *hemuDphimljMetj2;

    TH2F *hmumuDphimljMetj;
    TH2F *hmumuDphimljMetj0;
    TH2F *hmumuDphimljMetj1;
    TH2F *hmumuDphimljMetj2;


    TH2F *hhteeDphimljMetj;
    TH2F *hhteeDphimljMetj0;
    TH2F *hhteeDphimljMetj1;
    TH2F *hhteeDphimljMetj2;

    TH2F *hhtemuDphimljMetj;
    TH2F *hhtemuDphimljMetj0;
    TH2F *hhtemuDphimljMetj1;
    TH2F *hhtemuDphimljMetj2;

    TH2F *hhtmumuDphimljMetj;
    TH2F *hhtmumuDphimljMetj0;
    TH2F *hhtmumuDphimljMetj1;
    TH2F *hhtmumuDphimljMetj2;

    TH2F *hqeeDphimljMetj;
    TH2F *hqeeDphimljMetj0;
    TH2F *hqeeDphimljMetj1;
    TH2F *hqeeDphimljMetj2;
    TH2F *hqeeDphillj2;
    TH2F *hqeeDphillj0;

    TH2F *hqemuDphimljMetj;
    TH2F *hqemuDphimljMetj0;
    TH2F *hqemuDphimljMetj1;
    TH2F *hqemuDphimljMetj2;
    TH2F *hqemuDphillj2;
    TH2F *hqemuDphillj0;

    TH2F *hqmumuDphimljMetj;
    TH2F *hqmumuDphimljMetj0;
    TH2F *hqmumuDphimljMetj1;
    TH2F *hqmumuDphimljMetj2;
    TH2F *hqmumuDphillj2;
    TH2F *hqmumuDphillj0;

    TH2F *nothhteeDphimljMetj;
    TH2F *nothhteeDphimljMetj0;
    TH2F *nothhteeDphimljMetj1;
    TH2F *nothhteeDphimljMetj2;

    TH2F *nothhtemuDphimljMetj;
    TH2F *nothhtemuDphimljMetj0;
    TH2F *nothhtemuDphimljMetj1;
    TH2F *nothhtemuDphimljMetj2;

    TH2F *nothhtmumuDphimljMetj;
    TH2F *nothhtmumuDphimljMetj0;
    TH2F *nothhtmumuDphimljMetj1;
    TH2F *nothhtmumuDphimljMetj2;

    TH2F *nothqeeDphimljMetj;
    TH2F *nothqeeDphimljMetj0;
    TH2F *nothqeeDphimljMetj1;
    TH2F *nothqeeDphimljMetj2;
    TH2F *nothqeeDphillj2;
    TH2F *nothqeeDphillj0;

    TH2F *nothqemuDphimljMetj;
    TH2F *nothqemuDphimljMetj0;
    TH2F *nothqemuDphimljMetj1;
    TH2F *nothqemuDphimljMetj2;
    TH2F *nothqemuDphillj2;
    TH2F *nothqemuDphillj0;

    TH2F *nothqmumuDphimljMetj;
    TH2F *nothqmumuDphimljMetj0;
    TH2F *nothqmumuDphimljMetj1;
    TH2F *nothqmumuDphimljMetj2;
    TH2F *nothqmumuDphillj2;
    TH2F *nothqmumuDphillj0;

    TH1F *hdphijele;
    TH1F *hRjele;

    TH1F *hallHtj0;
    TH1F *hallHtj2;
    TH1F *hallHtj; 
    TH1F *hallHtj1inc; 

    TH1F *heeHt;
    TH1F *heeHtTow;
    TH1F *heeHtj0;
    TH1F *heeHtj1;
    TH1F *heeHtj1inc;
    TH1F *heeHtj2;

    TH1F *hmumuHt;
    TH1F *hmumuHtTow;
    TH1F *hmumuHtj0;
    TH1F *hmumuHtj1;
    TH1F *hmumuHtj1inc;
    TH1F *hmumuHtj2;

    TH1F *hemuHt;
    TH1F *hemuHtTow;
    TH1F *hemuHtj0;
    TH1F *hemuHtj1;
    TH1F *hemuHtj1inc;
    TH1F *hemuHtj2;


    //Histograms questions pre-blessing: 

    TH1F *h_Ete_j0met20;
    TH1F *h_Ht_j0met20;
    TH1F *h_Htbis_j0met20;
    TH1F *h_dRlj;

    TH2F *h_Ht_Met_j0;
    TH2F *h_Etee_Met_j0;
    TH2F *h1_Etee_PhiMet_j0;
    TH2F *h2_Etee_PhiMet_j0;
    TH2F *h_Ht_EteeMet_j0;

    TH2F *h_Ete1_Ete2;

    //ww candidates after all cuts
    TH1F *h_dphimetdee_wwall;
    TH1F *h_dphillmet_wwall;
    TH1F *h_met_wwall;
    TH1F *h_zmass_wwall;
    TH1F *h_met_wwallbht;
    TH1F *h_zmass_wwallbht;
    TH1F *h_ptllPair_wwall;
    TH1F *h_ptl_wwall;
 
    TH1F *h_dphimetdee_wwee;
    TH1F *h_dphieemet_wwee;
    TH1F *h_met_wwee;
    TH1F *h_zmass_wwee;
    TH1F *h_met_wweebht;
    TH1F *h_zmass_wweebht;
    TH1F *h_pteePair_wwee;
    TH1F *h_pte_wwee;
 
    TH1F *h_dphimetdee_wwmumu;
    TH1F *h_dphimumumet_wwmumu;
    TH1F *h_met_wwmumu;
    TH1F *h_zmass_wwmumu;
    TH1F *h_met_wwmumubht;
    TH1F *h_zmass_wwmumubht;
    TH1F *h_ptmumuPair_wwmumu;
    TH1F *h_ptmu_wwmumu;
 
    TH1F *h_zmass_wwemu;
    TH1F *h_zmass_wwemubht;
    TH1F *h_dphimetdee_wwemu;
    TH1F *h_dphiemumet_wwemu;
    TH1F *h_met_wwemu;
    TH1F *h_met_wwemubht;
    TH1F *h_ptemuPair_wwemu;
    TH1F *h_ptl_wwemu;

    //tt candidates after all cuts
    TH1F *h_dphimetdee_ttall;
    TH1F *h_dphillmet_ttall;
    TH1F *h_met_ttall;
    TH1F *h_zmass_ttall;
    TH1F *h_met_ttallbht;
    TH1F *h_zmass_ttallbht;
    TH1F *h_ptllPair_ttall;
    TH1F *h_ptl_ttall;
 
    TH1F *h_dphimetdee_ttee;
    TH1F *h_dphieemet_ttee;
    TH1F *h_met_ttee;
    TH1F *h_zmass_ttee;
    TH1F *h_met_tteebht;
    TH1F *h_zmass_tteebht;
    TH1F *h_pteePair_ttee;
    TH1F *h_pte_ttee;
 
    TH1F *h_dphimetdee_ttmumu;
    TH1F *h_dphimumumet_ttmumu;
    TH1F *h_met_ttmumu;
    TH1F *h_zmass_ttmumu;
    TH1F *h_met_ttmumubht;
    TH1F *h_zmass_ttmumubht;
    TH1F *h_ptmumuPair_ttmumu;
    TH1F *h_ptmu_ttmumu;
 
    TH1F *h_zmass_ttemu;
    TH1F *h_zmass_ttemubht;
    TH1F *h_dphimetdee_ttemu;
    TH1F *h_dphiemumet_ttemu;
    TH1F *h_met_ttemu;
    TH1F *h_met_ttemubht;
    TH1F *h_ptemuPair_ttemu;
    TH1F *h_ptl_ttemu;

    //1j candidates after all cuts
    TH1F *h_dphimetdee_1jall;
    TH1F *h_dphillmet_1jall;
    TH1F *h_met_1jall;
    TH1F *h_zmass_1jall;
    TH1F *h_ptllPair_1jall;
    TH1F *h_ptl_1jall;
 
    TH1F *h_dphimetdee_1jee;
    TH1F *h_dphieemet_1jee;
    TH1F *h_met_1jee;
    TH1F *h_zmass_1jee;
    TH1F *h_pteePair_1jee;
    TH1F *h_pte_1jee;
 
    TH1F *h_dphimetdee_1jmumu;
    TH1F *h_dphimumumet_1jmumu;
    TH1F *h_met_1jmumu;
    TH1F *h_zmass_1jmumu;
    TH1F *h_ptmumuPair_1jmumu;
    TH1F *h_ptmu_1jmumu;
 
    TH1F *h_zmass_1jemu;
    TH1F *h_dphimetdee_1jemu;
    TH1F *h_dphiemumet_1jemu;
    TH1F *h_met_1jemu;
    TH1F *h_ptemuPair_1jemu;
    TH1F *h_ptl_1jemu;

    TH2F *h_pte1pte2_wwee;
    TH2F *h_ptmu1ptmu2_wwmumu; 
    TH2F *h_pteptmu_wwemu;
    TH2F *h_pte1pte2_1jee;
    TH2F *h_ptmu1ptmu2_1jmumu; 
    TH2F *h_pteptmu_1jemu;
    TH2F *h_pte1pte2_ttee;
    TH2F *h_ptmu1ptmu2_ttmumu; 
    TH2F *h_pteptmu_ttemu;
    TH2F *h_pte1pte2_tteebht;
    TH2F *h_ptmu1ptmu2_ttmumubht; 
    TH2F *h_pteptmu_ttemubht;

    TH1F *hepg_im34;
    TH1F *hepg_im56;
    TH1F *hepg_pt34;
    TH1F *hepg_pt56;
    TH1F *hepg_im3456;
    TH1F *hepg_pt3456;
    TH1F *hepg_imw;
    TH1F *hepg_ptw;
    TH1F *hepg_E3;
    TH1F *hepg_E4;
    TH1F *hepg_E5;
    TH1F *hepg_E6;
    TH1F *hepg_Pt3;
    TH1F *hepg_Pt4;
    TH1F *hepg_Pt5;
    TH1F *hepg_Pt6;
    TH1F *hepg_Ptl;
    TH1F *hepg_Ptmin;
    TH1F *hepg_Ptmax;
    TH1F *hepg_Ptnunu;
    TH1F *hepg_dR35;
    TH1F *hepg_dphi35;
    TH1F *hepg_dphi565;
    TH1F *hepg_eta56arest;
	// Dean adds:
    TH1F *hepg_im35;
    TH1F *hepg_pt35;

  };

  struct Title{
    char plotT[80];
  };

  Title titleList[10];

  TCanvas *c1;
  TPostScript *ps;

  Hist_t* GetLeptonId(){return &fHist;}

  void BookLeptonId(void);
  void WriteHistoFile(void);

 private:

  Hist_t fHist;

  void LegendFit(const TH1F* histo, Double_t x, Double_t y);
};
#endif