#ifndef _NEUTVECT_H_
#define _NEUTVECT_H_

#include <TArrayI.h>

#include <TObject.h>
#include <TLorentzVector.h>

#include "neutpart.h"
#include "neutfsipart.h"
#include "neutfsivert.h"

class NeutVect : public TObject {

 public:
  NeutVect(Int_t np = 0, Int_t nfsip = 0, Int_t nfsiv = 0);

  ~NeutVect();

  void ClearVars();

  void Dump();

  /* Generated vector information */

  Int_t          EventNo;  // Event number

  /* Target parameters */
  Int_t          TargetA;   // Atomic number of the target nucleus
  Int_t          TargetZ;   // Mass number of the target nucleus
  Int_t          TargetH;   // Number of free protons
  Int_t          Ibound;    // Interaction on nucleus (1) or free proton (0)

  /* Nuclear potentials */
  Double_t       VNuclIni;  // Nuclear potential for the initial state (MeV)
                            //                                     
  Double_t       VNuclFin;  // Nuclear potential for the final state (MeV)
                            //                                     

  /* Fermi momentum related */
  Double_t       PFSurf;    // Fermi surface momentum (MeV/c)
  Double_t       PFMax;     // Maximum value of the Fermi momentum (MeV/c)
                            //                                 

  /* Flux ID including detector code?*/
  Int_t          FluxID;    // Flux ID code
                            //  0: Randomly generated

  // Interaction model
  Int_t          QEModel;
  Int_t          SPIModel;
  Int_t          COHModel;
  Int_t          DISModel;

  // 1 pi & the other meson productions form factor type
  Int_t          SPIForm;
  Int_t          RESForm;

  // 1 pi & the other meson productions NR type ( form factor detail )
  Int_t          SPINRType;
  Int_t          RESNRType;

  Int_t          QEVForm;
  Float_t        QEMA;
  Float_t        SPIMA;
  Float_t        RESMA;
  Float_t        QEMV;
  Float_t        SPIMV;
  Float_t        RESMV;
  Float_t        KAPPA;
  Float_t        COHMA;
  Float_t        COHR0;
  

  // 1 pi CA5I and BG scale
  Float_t        SPICA5I;
  Float_t        SPIBGScale;
	

  // NEUT version
  Int_t          COREVer;
  Int_t          NUCEVer;
  Int_t          NUCCVer;
  // NEUT card configuration
  Int_t          FrmFlg;
  Int_t          PauFlg;
  Int_t          NefO16;
  Int_t          ModFlg;
  Int_t          SelMod;
  Int_t          FormLen;
  Int_t          IPilessDcy;
  Float_t        RPilessDcy;
  Int_t          NucScat;
  Float_t        NucFac;

  Float_t        NuceffFactorPIQE;
  Float_t        NuceffFactorPIInel;
  Float_t        NuceffFactorPIAbs;
  Float_t        NuceffFactorPIQEH;
  Float_t        NuceffFactorPICX;
  Float_t        NuceffFactorPICXH;
  Float_t        NuceffFactorPICoh;

  Float_t        NuceffFactorPIQEHKin;
  Float_t        NuceffFactorPIQELKin;
  Float_t        NuceffFactorPICXKin;

  /////////////////////////////////////////////////////////////////////////
  /* Ineraction mode */
  Int_t          Mode;

  /* Cross sections */
  Float_t        Totcrs;

  Float_t        CrsEnergy;
  Float_t        DifCrsNE[8];
  Float_t        Crsx;
  Float_t        Crsy;
  Float_t        Crsz;
  Float_t        Crsphi;
  Float_t        Crsq2;

  /* Number of particles */
  Int_t          Npart(void)        const {return fNpart; }
  void           SetNpart(Int_t npart);

  /* Number of primary particles */
  Int_t          Nprimary(void)     const {return fNprimary; }
  void           SetNprimary(Int_t nprimary);

  /* Vertex ID */
  Int_t          VertexID(int idx)   const { return fVertexID[idx]; };
  void           SetVertexID(Int_t idx, Int_t VtxID);
  void           SetVertexID(Int_t np, Int_t *VtxID_array);

  /* Particle information (PID, 4 momentum etc.) */
  NeutPart       *PartInfo(Int_t idx)
	{ return (NeutPart *)((fNpart>idx) ? (fPartInfo->At(idx)) : NULL );};
  void SetPartInfo(Int_t idx      , NeutPart PInfo);
  void SetPartInfo(Int_t npart   , NeutPart *PInfo_array);

  /* Parent particle's index */
  Int_t          ParentIdx(int idx)  const { return fParentIdx[idx]; };
  void           SetParentIdx(Int_t idx, Int_t ParentIdx);
  void           SetParentIdx(Int_t np, Int_t *ParentIdx_array);


  /* Number of FSI particles */
  Int_t          NfsiPart(void)        const {return fNfsiPart; }
  void           SetNfsiPart(Int_t npart);

  /* FSI Particle Information */
  NeutFsiPart       *FsiPartInfo(Int_t idx)
	{ return (NeutFsiPart *)((fNfsiPart>idx) ? (fFsiPartInfo->At(idx)) : NULL );};
  void SetFsiPartInfo(Int_t idx      , NeutFsiPart PInfo);
  void SetFsiPartInfo(Int_t npart   , NeutFsiPart *PInfo_array);


  /* Number of FSI Vertices */
  Int_t          NfsiVert(void)        const {return fNfsiVert; }
  void           SetNfsiVert(Int_t nvert);

  /* FSI Vertex Information */
  NeutFsiVert       *FsiVertInfo(Int_t idx)
	{ return (NeutFsiVert *)((fNfsiVert>idx) ? (fFsiVertInfo->At(idx)) : NULL );};
  void SetFsiVertInfo(Int_t idx      , NeutFsiVert VInfo);
  void SetFsiVertInfo(Int_t nvert   , NeutFsiVert *VInfo_array);

  Float_t Fsiprob;

 private:

  Int_t          fNpart;     // Number of particles 

  Int_t          fNprimary;  // Number of primary particles

  TObjArray      *fPartInfo; // ->

  TArrayI        fVertexID;  //         Vertex Index ID 
                             //         ( To refer NeutVertex )

  TArrayI        fParentIdx;  //        Idx of parent particle

  Int_t          fRandomSeed[16]; // Random seeds

  Int_t          fNfsiPart;  // Number of FSI particles
  Int_t          fNfsiVert;  // Number of FSI vertices
  TObjArray      *fFsiPartInfo; // ->
  TObjArray      *fFsiVertInfo; // ->

  TLorentzVector fZeroVect;  // !
  
  ClassDef(NeutVect, 2)
};

#endif