//STL 'n' stuff
#include <iostream>
//Root includes
#include <TROOT.h>
#include <TChain.h>

using namespace std;
int main(){

  const int kMaxm_genParticles = 4214;
  Int_t           m_genParticles_;
  Float_t         m_genParticles_m_px[kMaxm_genParticles];   //[m_genParticles_]
  Float_t         m_genParticles_m_py[kMaxm_genParticles];   //[m_genParticles_]
  Float_t         m_genParticles_m_pz[kMaxm_genParticles];   //[m_genParticles_]

  TBranch        *b_McEventCollection_p4_GEN_EVENT_m_genParticles_;   //!
  TBranch        *b_m_genParticles_m_px;   //!
  TBranch        *b_m_genParticles_m_py;   //!
  TBranch        *b_m_genParticles_m_pz;   //!

  TChain* chain = new TChain("CollectionTree"); 

  chain->Add("/Users/rasmack/root/ntuples/alpgen/mc08.107682.AlpgenJimmyWenuNp2_pt20.evgen.EVNT.e349/EVNT.024150._00001.pool.root.1");  
  chain->Add("/Users/rasmack/root/ntuples/alpgen/mc08.107682.AlpgenJimmyWenuNp2_pt20.evgen.EVNT.e349/EVNT.024150._00002.pool.root.1");

  chain->SetBranchAddress("m_genParticles", &m_genParticles_, &b_McEventCollection_p4_GEN_EVENT_m_genParticles_);
  chain->SetBranchAddress("m_genParticles.m_px", m_genParticles_m_px, &b_m_genParticles_m_px);
  chain->SetBranchAddress("m_genParticles.m_py", m_genParticles_m_py, &b_m_genParticles_m_py);
  chain->SetBranchAddress("m_genParticles.m_pz", m_genParticles_m_pz, &b_m_genParticles_m_pz);


  for (int i = 0;i!=chain->GetEntries();i++)
    {
      chain->GetEntry(i);
      cout<<m_genParticles_<<endl;
    }
  return 0;
}