#define ww_lepton_id_cxx
#include "ww_lepton_id.h"
#include "TStyle.h"
#include "TCanvas.h" 
#include "TH1.h"
#include "TH2.h"
#include "TF1.h"
#include <TText.h>
#include <TFile.h>
#include <TSystem.h>
#include <TPaveLabel.h>
//
void ww_lepton_id::NumberOfEvents()
{
  if (fChain == 0) return;
  cout << "\n Number of events = " << Int_t(fChain->GetEntries()) << " \n" << endl;
}

 
//


void ww_lepton_id::Loop(int     NumEvt,
                       int     flag,
                       double  cut1,
                       double  cut2,
                       char   *file_name,
                       char   *directory) 
{
//   In a ROOT session, you can do:
//      Root > .L ww_lepton_id.C
//      Root > ww_lepton_id t
//      Root > t.GetEntry(12); // Fill t data members with entry number 12
//      Root > t.Show();       // Show values of entry 12
//      Root > t.Show(16);     // Read and show values of entry 16
//      Root > t.Loop();       // Loop on all entries
//

//     This is the loop skeleton
//       To read only selected branches, Insert statements like:
// METHOD1:
//    fChain->SetBranchStatus("*",0);  // disable all branches
//    fChain->SetBranchStatus("branchname",1);  // activate branchname
// METHOD2: replace line
//    fChain->GetEntry(i);  // read all branches
//by  b_branchname->GetEntry(i); //read only this branch
   if (fChain == 0) return;

  

   histoLeptonId histoww;
   hww=histoww.GetLeptonId();
   histoww.BookLeptonId();
  

   Int_t nentries = Int_t(fChain->GetEntriesFast());


   //Int_t nentries = Int_t(fChain->GetEntries());

   Int_t  IdRecEl[MAXLEPTON];
   Int_t  IdRecMu[MAXLEPTON];
   Int_t  NRecEl;
   Int_t  NRecMu;
   Int_t  Irec;
   Int_t  Irecflag;

   Int_t Ievprint=0;
   
   IevAna=0;
   IFileAna=0;

   Nttee7b=0;
   Nttee8b=0;
   Nttee9b=0;
   //
    for(Int_t i=0;i<MAXttbar;i++){
      Nttee[i]=0;
      //
      Nttemu[i]=0;
      Nttemu20[i]=0;
      Nttemu21[i]=0;
      Nttemu22[i]=0;
      Nttemu23[i]=0;
      //
      Nttmumu[i]=0;
      Nttmumu10[i]=0;
      Nttmumu11[i]=0;
      Nttmumu12[i]=0;
      Nttmumu13[i]=0;
      Nttmumu14[i]=0;
      Nttmumu15[i]=0;
      Nttmumu16[i]=0;
      Nttmumu17[i]=0;
      Nttmumu18[i]=0;
      Nttmumu19[i]=0;
      //
      Nwwee[i]=0;
      //
      Nwwemu[i]=0;
      Nwwemu20[i]=0;
      Nwwemu21[i]=0;
      Nwwemu22[i]=0;
      Nwwemu23[i]=0;
      //
      Nwwmumu[i]=0;
      Nwwmumu10[i]=0;
      Nwwmumu11[i]=0;
      Nwwmumu12[i]=0;
      Nwwmumu13[i]=0;
      Nwwmumu14[i]=0;
      Nwwmumu15[i]=0;
      Nwwmumu16[i]=0;
      Nwwmumu17[i]=0;
      Nwwmumu18[i]=0;
      Nwwmumu19[i]=0;
      //
      N1jee[i]=0;
      //
      N1jemu[i]=0;
      N1jemu20[i]=0;
      N1jemu21[i]=0;
      N1jemu22[i]=0;
      N1jemu23[i]=0;
      //
      N1jmumu[i]=0;
      N1jmumu10[i]=0;
      N1jmumu11[i]=0;
      N1jmumu12[i]=0;
      N1jmumu13[i]=0;
      N1jmumu14[i]=0;
      N1jmumu15[i]=0;
      N1jmumu16[i]=0;
      N1jmumu17[i]=0;
      N1jmumu18[i]=0;
      N1jmumu19[i]=0;
      //
    }
   
   //
   Nttee_ht_100=0;
   Nttee_ht_100_os=0;
   Nttee_ht_150=0;
   Nttee_ht_150_os=0;
   Nttee_ht_200=0;
   Nttee_ht_200_os=0;
   Nttee_ht_250=0;
   Nttee_ht_250_os=0;    
   //
   Nttmumu_ht_100=0;
   Nttmumu_ht_100_os=0;
   Nttmumu_ht_150=0;
   Nttmumu_ht_150_os=0;
   Nttmumu_ht_200=0;
   Nttmumu_ht_200_os=0;
   Nttmumu_ht_250=0;
   Nttmumu_ht_250_os=0;    
   //
   Nttmumu10_ht_100=0;
   Nttmumu10_ht_100_os=0;
   Nttmumu10_ht_150=0;
   Nttmumu10_ht_150_os=0;
   Nttmumu10_ht_200=0;
   Nttmumu10_ht_200_os=0;
   Nttmumu10_ht_250=0;
   Nttmumu10_ht_250_os=0;    
   //
   Nttmumu11_ht_100=0;
   Nttmumu11_ht_100_os=0;
   Nttmumu11_ht_150=0;
   Nttmumu11_ht_150_os=0;
   Nttmumu11_ht_200=0;
   Nttmumu11_ht_200_os=0;
   Nttmumu11_ht_250=0;
   Nttmumu11_ht_250_os=0;    
   //
   Nttmumu12_ht_100=0;
   Nttmumu12_ht_100_os=0;
   Nttmumu12_ht_150=0;
   Nttmumu12_ht_150_os=0;
   Nttmumu12_ht_200=0;
   Nttmumu12_ht_200_os=0;
   Nttmumu12_ht_250=0;
   Nttmumu12_ht_250_os=0;    
   //
   Nttmumu16_ht_100=0;
   Nttmumu16_ht_100_os=0;
   Nttmumu16_ht_150=0;
   Nttmumu16_ht_150_os=0;
   Nttmumu16_ht_200=0;
   Nttmumu16_ht_200_os=0;
   Nttmumu16_ht_250=0;
   Nttmumu16_ht_250_os=0;    
   //
   Nttmumu17_ht_100=0;
   Nttmumu17_ht_100_os=0;
   Nttmumu17_ht_150=0;
   Nttmumu17_ht_150_os=0;
   Nttmumu17_ht_200=0;
   Nttmumu17_ht_200_os=0;
   Nttmumu17_ht_250=0;
   Nttmumu17_ht_250_os=0;    
   //
   Nttmumu18_ht_100=0;
   Nttmumu18_ht_100_os=0;
   Nttmumu18_ht_150=0;
   Nttmumu18_ht_150_os=0;
   Nttmumu18_ht_200=0;
   Nttmumu18_ht_200_os=0;
   Nttmumu18_ht_250=0;
   Nttmumu18_ht_250_os=0;    
   //
   Nttmumu19_ht_100=0;
   Nttmumu19_ht_100_os=0;
   Nttmumu19_ht_150=0;
   Nttmumu19_ht_150_os=0;
   Nttmumu19_ht_200=0;
   Nttmumu19_ht_200_os=0;
   Nttmumu19_ht_250=0;
   Nttmumu19_ht_250_os=0;    
   //
   //
   Nttemu_ht_100=0;
   Nttemu_ht_100_os=0;
   Nttemu_ht_150=0;
   Nttemu_ht_150_os=0;
   Nttemu_ht_200=0;
   Nttemu_ht_200_os=0;
   Nttemu_ht_250=0;
   Nttemu_ht_250_os=0;    
   //
   //
   Nttemu20_ht_100=0;
   Nttemu20_ht_100_os=0;
   Nttemu20_ht_150=0;
   Nttemu20_ht_150_os=0;
   Nttemu20_ht_200=0;
   Nttemu20_ht_200_os=0;
   Nttemu20_ht_250=0;
   Nttemu20_ht_250_os=0;    
   //
   //
   Nttemu21_ht_100=0;
   Nttemu21_ht_100_os=0;
   Nttemu21_ht_150=0;
   Nttemu21_ht_150_os=0;
   Nttemu21_ht_200=0;
   Nttemu21_ht_200_os=0;
   Nttemu21_ht_250=0;
   Nttemu21_ht_250_os=0;    
   //
   //
   Nttemu22_ht_100=0;
   Nttemu22_ht_100_os=0;
   Nttemu22_ht_150=0;
   Nttemu22_ht_150_os=0;
   Nttemu22_ht_200=0;
   Nttemu22_ht_200_os=0;
   Nttemu22_ht_250=0;
   Nttemu22_ht_250_os=0;    
   //
   //
   Nttemu23_ht_100=0;
   Nttemu23_ht_100_os=0;
   Nttemu23_ht_150=0;
   Nttemu23_ht_150_os=0;
   Nttemu23_ht_200=0;
   Nttemu23_ht_200_os=0;
   Nttemu23_ht_250=0;
   Nttemu23_ht_250_os=0;
   //
   Nwwee_ht_150_os=0;
   Nwwee_ht_200_os=0;
   Nwwee_ht_250_os=0;
   Nwwee_ht_300_os=0;
   //
   Nwwmumu_ht_150_os=0;
   Nwwmumu_ht_200_os=0;
   Nwwmumu_ht_250_os=0;
   Nwwmumu_ht_300_os=0;
   //
   Nwwemu_ht_150_os=0;
   Nwwemu_ht_200_os=0;
   Nwwemu_ht_250_os=0;
   Nwwemu_ht_300_os=0;
   //
   Nwwee_htow_150_os=0;
   Nwwee_htow_200_os=0;
   Nwwee_htow_250_os=0;
   Nwwee_htow_300_os=0;
   //
   Nwwmumu_htow_150_os=0;
   Nwwmumu_htow_200_os=0;
   Nwwmumu_htow_250_os=0;
   Nwwmumu_htow_300_os=0;
   //
   Nwwemu_htow_150_os=0;
   Nwwemu_htow_200_os=0;
   Nwwemu_htow_250_os=0;
   Nwwemu_htow_300_os=0;
   //
   Nwwee_ht1j_150_os=0;
   Nwwee_ht1j_200_os=0;
   Nwwee_ht1j_250_os=0;
   Nwwee_ht1j_300_os=0;
   //
   Nwwmumu_ht1j_150_os=0;
   Nwwmumu_ht1j_200_os=0;
   Nwwmumu_ht1j_250_os=0;
   Nwwmumu_ht1j_300_os=0;
   //
   Nwwemu_ht1j_150_os=0;
   Nwwemu_ht1j_200_os=0;
   Nwwemu_ht1j_250_os=0;
   Nwwemu_ht1j_300_os=0;
   //
   //
   ssee=0;
   ssmumu=0;
   ssemu=0; 
   //
   GeomTowEta();
   PrintGeomTowEta();
   //
   Int_t nbytes = 0, nb = 0;
   for (Int_t jentry=0; jentry<nentries;jentry++) {
     //if(jentry==6450||jentry==9232){jentry++;}
      Int_t ientry = LoadTree(jentry); //in case of a TChain, ientry is the entry number in the current file
      if (ientry < 0) break;
      nb = fChain->GetEntry(jentry);   nbytes += nb;
      //
      //
      //WWDileptonFilter();
      WWDileptonFilterMadGraph();
      //      
      //
      //
      // DileptonCategoryAlgorithm(IdRecEl,IdRecMu,NRecEl,NRecMu);
      //
      //      PrintDileptonCategoryAlgorithm(IdRecEl,IdRecMu,NRecEl,NRecMu);
      //
      // Irec=ClassifyDileptonCategory(NRecEl,NRecMu);
      //
      // Irecflag=LabelDileptonCategory(IdRecEl,IdRecMu,NRecEl,NRecMu,Irec);
      //
      // Rec_Source(IdRecEl,IdRecMu,NRecEl,NRecMu,Irecflag,Irec);
      //
      // DebugBaseLinePlugElectron();
      // DebugBaseLinePlugElectronCalDriven();
      //
      if(jentry % 1000==0) { 
	cout << jentry << " -> events analyzed, file= " << fCurrent << endl; 
      }
      if(jentry % 5000==0) { 
        cout << "--------------------------------" << endl;  
      }
      //
      Ievprint++;
      IevAna++;
   }
   //
   //
   //
   cout << "--------------------------------" << endl; 
   cout << Ievprint << " -> events analyzed " << endl; 
   cout << "--------------------------------" << endl; 
   acceptancettbar(Ievprint);    
   //
   histoww.WriteHistoFile();
   
   
}