#ifndef JNuBeamFlux_hh_seen
#define JNuBeamFlux_hh_seen

#include <iostream>

#include "TObject.h"
#include "TObjString.h"

///
/// JNuBeam flux pass-through info. Both GRooTracker and NRooTracker inherit
/// from this so as to avoid duplication of code. Currently up to date with 
/// 11a version of pass-through info.
///
namespace COMET {

  const int kNgmax = 12; 

  // This is the base class for all RooTracker type objects. For now
  // it is  empty and just so that an pass round pointers to
  // all derived classes (will be usefull for reweighting applications etc...)
  class IRooTrackerVtxBase : public TObject {
    public:
      IRooTrackerVtxBase();
      virtual ~IRooTrackerVtxBase(){}
      ClassDef(COMET::IRooTrackerVtxBase, 1);
  };

  class IJNuBeamFlux : public IRooTrackerVtxBase {
    
  public:
    IJNuBeamFlux();
   virtual ~IJNuBeamFlux();
    
    void Reset   (void);
    // Using methods from TObject to remove 'hidden' compiler warnings
    using TObject::Copy;
    void Copy    (const IJNuBeamFlux * flux);
    
    //
    // Meta information
    //
    long        NuFluxEntry;        ///< corresponding entry in orig flux file
    TObjString* NuFileName;         ///< orig flux file name
    //
    // neutrino parent info
    //
    int         NuParentPdg;        ///< parent hadron pdg code
    int         NuParentDecMode;    ///< parent hadron decay mode
    double      NuParentDecP4 [4];  ///< parent hadron 4-momentum at decay
    double      NuParentDecX4 [4];  ///< parent hadron 4-position at decay 
    float       NuCospibm;          ///< parent hadron cos(theta) relative to beam at decay
    float       NuNorm;             ///< Normalization weight 
    double      NuParentProP4 [4];  ///< parent hadron 4-momentum at production 
    double      NuParentProX4 [4];  ///< parent hadron 4-position at production  
    float       NuCospi0bm;         ///< parent hadron cos(theta) relative to beam at production
    float       NuRnu;              ///< neutrino r position at ND5/6 plane 
    float       NuXnu [2];          ///< neutrino (x,y) position at ND5/6 plane 
    int         NuIdfd;             ///< Detector ID
    //
    // primary particle information
    // 
    int         NuGipart;           ///< primary particle ID  
    float       NuGpos0[3];         ///< primary particle starting point 
    float       NuGvec0[3];         ///< primary particle direction at starting point 
    float       NuGamom0;           ///< momentum of the primary particle at the starting point 
    //
    // Interation History Information
    //
    int         NuNg;               ///< Number of interaction steps 
    float       NuGp[kNgmax][3];    ///< Momentum of ancestor particle 
    float       NuGcosbm[kNgmax];   ///< Ancestor particle cos(theta) relative to beam
    float       NuGv[kNgmax][3];    ///< Position of ancestor particle 
    int         NuGpid[kNgmax];     ///< PDG ancestor particle ID 
    int         NuGmec[kNgmax];     ///< particle production mechanism of ancestor particle
     //
    // Normalization and Transfer Matrix Information
    //
    float       NuEnusk;            ///< Neutrino Energy at SK 
    float       NuNormsk;           ///< Normalization weight at SK 
    float       NuAnorm;            ///< ND Acceptance Weight 
    //
    // Out-of-target Secondary Interactions
    //
    int         NuGmat[kNgmax];     ///< material in which the particle originates
    float       NuGdistc[kNgmax];   ///< distance traveled through carbon
    float       NuGdistal[kNgmax];  ///< distance traveled through aluminum
    float       NuGdistti[kNgmax];  ///< distance traveled through titanium
    float       NuGdistfe[kNgmax];  ///< distance traveled through iron
    //                                  
    // Beam parameter information 
    //                                     
    float       NuVersion;          ///< jnubeam flux version 
    int         NuTuneid;           ///< beam tune ID # 
    int         NuNtrig;            ///< Number of triggers
    int         NuPint;             ///< Interaction model ID 
    float       NuBpos[2];          ///< Beam center position 
    float       NuBtilt[2];         ///< Beam angle  
    float       NuBrms[2];          ///< Beam RMS width 
    float       NuEmit[2];          ///< Beam Emittance 
    float       NuAlpha[2];         ///< Beam Alpha 
    float       NuHcur[3];          ///< Horn currents 
    int         NuRand;             ///< Random seed
//    int         NuRseed[2];         // Random seed

    ClassDef(COMET::IJNuBeamFlux, 1);
  };
  
} // comet namespace
#endif