TGRooTrackerVtx.hxx

Go to the documentation of this file.

#ifndef GRooTrackerVtx_hh_seen
#define GRooTrackerVtx_hh_seen

#include <iostream>

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

#include "TJNuBeamFlux.hxx"

using std::ostream;

///
/// This is a simple event class which is essentially an objectified version
/// of the GENIE gRooTracker output format (based on StdHep format).
/// This can be used with the GRooTracker utils class in AnalysisTools as well
/// as GENIE reweighting tools.
///

namespace ND {

const int kGStdHepNPmax = 350;
const int kGStdHepIdxPx = 0;
const int kGStdHepIdxPy = 1;
const int kGStdHepIdxPz = 2;
const int kGStdHepIdxE = 3;
const int kGStdHepIdxX = 0;
const int kGStdHepIdxY = 1;
const int kGStdHepIdxZ = 2;
const int kGStdHepIdxT = 3;

class GRooTrackerVtx : public JNuBeamFlux {
 public:
  GRooTrackerVtx();
  ~GRooTrackerVtx();

  void Reset(void);
  void Init(void);

  // Using methods from TObject to remove 'hidden' compiler warnings
  using TObject::Copy;
  void Copy(const GRooTrackerVtx* event);
  void Print(const Option_t* option = "") const;

  // Define the output rootracker tree branches
  TObjString* EvtCode;  ///< generator-specific string with 'event code'
  int EvtNum;           ///< event num.
  double EvtXSec;       ///< cross section for selected event (1E-38 cm2)
  double EvtDXSec;   ///< cross section for selected event kinematics (1E-38 cm2{K^n})
  double EvtWght;    ///< weight for that event
  double EvtProb;    ///< probability for that event (given cross section, pathlengths, etc)
  double EvtVtx[4];  ///< event vertex position in detector coord syst (SI)
  int StdHepN;       ///< number of particles in particle array  stdhep-like particle array
  int StdHepPdgTemp[kGStdHepNPmax];  //! pdg codes (& generator specific codes
                                     //! for pseudoparticles)
  int* StdHepPdg;  //[StdHepN] /< pdg codes (& generator specific codes for
                   //pseudoparticles)

  int StdHepRescatTemp[kGStdHepNPmax];  //! generator-specific status code
  int* StdHepRescat;  //[StdHepN] /< generator-specific status code

  int StdHepStatusTemp[kGStdHepNPmax];  //! generator-specific status code
  int* StdHepStatus;  //[StdHepN] /< generator-specific status code

  double
      StdHepX4[kGStdHepNPmax]
              [4];  ///< 4-x (x, y, z, t) of particle in hit nucleus frame (fm)
  double StdHepP4[kGStdHepNPmax]
                 [4];  ///< 4-p (px,py,pz,E) of particle in LAB frame (GeV)
  double StdHepPolz[kGStdHepNPmax][3];  ///< polarization vector

  int StdHepFdTemp[kGStdHepNPmax];  //! first daughter. First daughter of -1
                                    //! indicates final state particle (after
                                    //! FSI)
  int* StdHepFd;  //[StdHepN] /< first daughter. First daughter of -1 indicates
                  //final state particle (after FSI)

  int StdHepLdTemp[kGStdHepNPmax];  //! last  daughter
  int* StdHepLd;                    //[StdHepN] /< last  daughter

  int StdHepFmTemp[kGStdHepNPmax];  //! first mother. First mother of -1
                                    //! indicates initial particle (e.g.
                                    //! incoming neutrino, target nucleus)
  int* StdHepFm;  //[StdHepN] /< first mother. First mother of -1 indicates
                  //initial particle (e.g. incoming neutrino, target nucleus)

  int StdHepLmTemp[kGStdHepNPmax];  //! last  mother
  int* StdHepLm;                    //[StdHepN] /< last  mother

  //
  // neutrino parent info (passed-through from the beam-line MC / quantities in
  // 'jnubeam' units)
  // is defined in TJNuBeamFlux from which this class inherits.
  //

  //
  // etc
  //
  int G2NeutEvtCode;     ///< NEUT-like reaction code for the GENIE event
  TObjString* GeomPath;  ///< Geometry path of where the vertex is

  // Some pass through info
  TObjString* GeneratorName;  ///< genie, also gives some of the running options(e.g. filling events with a fixed or meannumber of vertices)
  TObjString* OrigFileName;   ///< Filename of the gnmc file
  TObjString* OrigTreeName;   ///< Tree name within the gnmc file
  int OrigEvtNum;             ///< Event number within the gnmc file
  int OrigTreeEntries;        ///< Number of entries in the gnmc file
  double OrigTreePOT;         ///< Number of POT in the gnmc file
  double TimeInSpill;         ///< Time of the vertex, relative to spill time

  int TruthVertexID;  ///< A link to the TTruthVertex tree

  ClassDef(ND::GRooTrackerVtx, 2);
};

}  // nd280 namespace
#endif