#define MyClass_cxx
#include "MyClass.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>

void MyClass::Loop()
{
//   In a ROOT session, you can do:
//      Root > .L MyClass.C
//      Root > MyClass 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 where:
//    jentry is the global entry number in the chain
//    ientry is the entry number in the current Tree
//  Note that the argument to GetEntry must be:
//    jentry for TChain::GetEntry
//    ientry for TTree::GetEntry and TBranch::GetEntry
//
//       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(jentry);       //read all branches
//by  b_branchname->GetEntry(ientry); //read only this branch
   if (fChain == 0) return;

   Long64_t nentries = fChain->GetEntriesFast();

   TH1F *myHisto = new TH1F("myHisto", "NumberOfJets", 10, 0, 10);
   TH1F *myHisto2 = new TH1F("myHisto2", "Jet Energy", 100, 0, 3000);

   Long64_t nbytes = 0, nb = 0;
   for (Long64_t jentry=0; jentry<nentries;jentry++) {
      Long64_t ientry = LoadTree(jentry);
      if (ientry < 0) break;
      nb = fChain->GetEntry(jentry);   nbytes += nb;
      // if (Cut(ientry) < 0) continue;

      myHisto->Fill(njets); 
      
      // Loop Over Jets
      for (Int_t j = 0; j< jet_E->size(); j++){
	myHisto2->Fill(jet_E->at(j));
      }
   }

   // Where we are last time
   // myHisto2->Draw();

   //////////////////////////
   // Uncomment for slide 3//
   //////////////////////////

//   myHisto2->SetName("Jet Energy");
//   myHisto2->SetFillColor(2);
//   myHisto2->SetLineColor(1);
//   myHisto2->SetLineWidth(3);
//   myHisto2->GetXaxis()->SetTitle("Jet Energy");
//   myHisto2->GetXaxis()->SetTitleSize(0.04);
//   myHisto2->GetYaxis()->SetTitle("Number of Jet / Bin");
//   myHisto2->GetYaxis()->SetTitleSize(0.04);
//   myHisto2->GetYaxis()->SetRangeUser(0.01, 250000);
//   myHisto2->Draw();

   //////////////////////////
   // Uncomment for slide 4//
   //////////////////////////

//   std::cout<<"Histogram Mean "<< myHisto2->GetMean() <<std::endl;
//   std::cout<<"Histogram Standard Deviation "<< myHisto2->GetStdDev() <<std::endl;
//   std::cout<<"bin 1 has "<< myHisto2->GetBinContent(1)<<" jets" <<std::endl;
//   std::cout<<"bin 1 has an uncertainty of "<< myHisto2->GetBinError(1)<<" jets" <<std::endl;
//   
//   myHisto2->Divide(myHisto2);
//   myHisto2->Draw();


   //////////////////////////
   // Uncomment for slide 5//
   //////////////////////////

//   TCanvas *Canv_Jet = new TCanvas("Jet Energy", "Jet Envergy", 800, 600);
//
//   TH1F *H_LeadJetE = new TH1F("LeadJetE", "LeadJetE", 100, 0, 3000);
//   TH1F *H_SubLeadJetE  = new TH1F("SubLeadJetE", "SubLeadJetE", 100, 0, 3000);
//   H_LeadJetE->SetLineColor(1);
//   H_SubLeadJetE->SetLineColor(2);
//
//   Long64_t nbytes = 0, nb = 0;
//   for (Long64_t jentry=0; jentry<nentries;jentry++) {
//      Long64_t ientry = LoadTree(jentry);
//      if (ientry < 0) break;
//      nb = fChain->GetEntry(jentry);   nbytes += nb;
//      // if (Cut(ientry) < 0) continue;
//
//      // make sure the events have at least two jets
//      if (jet_E->size()>1){
//	// Histogram for leading jets
//	Int_t j = 0;
//	H_LeadJetE->Fill(jet_E->at(j));
//	j = 1;
//	H_SubLeadJetE->Fill(jet_E->at(j));
//      }
//   }

//   Canv_Jet->cd();
//   H_LeadJetE->Draw();
//   H_SubLeadJetE->Draw("same"); // first part
   
//   Canv_Jet->cd();
//   TPad *Pad1 = new TPad("Pad1", "Pad1", 0.0, 0.55, 1.0, 1.0, 0);
//   TPad *Pad2 = new TPad("Pad2", "Pad2", 0.0, 0.0, 1.0, 0.45, 0);
//   Pad1->Draw();
//   Pad2->Draw();
//   Pad1->cd();
//   H_LeadJetE->Draw();
//   Pad2->cd();
//   H_SubLeadJetE->Draw();

   //////////////////////////
   // Uncomment for slide 6//
   //////////////////////////

//   TH2F *H_LeadE_SubLeadE = new TH2F("LeadJetE_SubJetE", "LeadJetE_SubJetE", 10, 0, 3000, 10, 0, 3000);
//
//   Long64_t nbytes = 0, nb = 0;
//   for (Long64_t jentry=0; jentry<nentries;jentry++) {
//      Long64_t ientry = LoadTree(jentry);
//      if (ientry < 0) break;
//      nb = fChain->GetEntry(jentry);   nbytes += nb;
//      // if (Cut(ientry) < 0) continue;
//
//      // make sure the events have at least two jets
//      if (jet_E->size()>1){
//	// Histogram for leading jets
//	Int_t j = 0;
//	Int_t k = 1;
//	H_LeadE_SubLeadE->Fill(jet_E->at(j), jet_E->at(j));
//      }
//   }
//
//   H_LeadE_SubLeadE->Draw();


   //////////////////////////
   // Uncomment for slide 7//
   //////////////////////////

//   TH1F *H_DijetMass  = new TH1F("DijetMass", "DijetMass", 100, 0, 2000);
//   Long64_t nbytes = 0, nb = 0;
//   for (Long64_t jentry=0; jentry<nentries;jentry++) {
//      Long64_t ientry = LoadTree(jentry);
//      if (ientry < 0) break;
//      nb = fChain->GetEntry(jentry);   nbytes += nb;
//      // if (Cut(ientry) < 0) continue;
//
//      // make sure the events have at least two jets
//      if (jet_E->size()>1){
//	// Histogram for leading jets
//	TLorentzVector LeadJet;
//	TLorentzVector SubLeadJet;
//	TLorentzVector JetSum;
//	Int_t j = 0;
//	Int_t k = 1;
//	LeadJet->SetPtEtaPhiE(jet_pt->at(j), jet_eta->at(j), jet_phi->at(j), jet_E->at(j));
//	SubLeadJet->SetPtEtaPhiE(jet_pt->at(k), jet_eta->at(k), jet_phi->at(k), jet_E->at(k));
//	JetSum = LeadJet;
//	JetSum += SubLeadJet;
//	H_DijetMass->Fill(JetSum->M());
//      }
//   }
//   H_DijetMass->Draw();

}