#define inclphoton_cxx
#include "inclphoton.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <TLorentzVector.h>
#include <TRandom1.h>
#include <TRandom2.h>
#include <TRandom3.h>
#include <TMath.h>
//#include "FourMomUtils/xAODP4Helpers.h"
void inclphoton::Loop()
{
TRandom *r1=new TRandom();
TFile *inclphoton_Plots= new TFile("inclphoton_Plots.root","RECREATE");

//**********************************************Declare histograms****************************************************

TH1F *PholeadpT  = new TH1F("PhotonpT", "Photon_{p_{T}} (GeV)", 50, 0.0, 500.0);
  std::vector<float>* phleadingpt;

//******************************************************************************************************************
   if (fChain == 0) return;
   Long64_t nentries = fChain->GetEntriesFast();
   Long64_t nbytes = 0, nb = 0;
   for (Long64_t ipp=0; ipp<nentries;ipp++) {
      Long64_t ientry = LoadTree(ipp);
      if (ientry < 0) break;
      nb = fChain->GetEntry(ipp);   nbytes += nb;
      std::cout<< "photon Pt:"<<ttdata_presel_photonPt->at(ipp)<<std::endl;

          PholeadpT->Fill(ttdata_presel_photonPt->at(ipp));
   }

    inclphoton_Plots->Write();
}