SRNeutrino.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \author  gsdavies@iastate.edu
// \date    February 2013
////////////////////////////////////////////////////////////////////////
#ifndef SRNEUTRINO_H
#define SRNEUTRINO_H

#include "StandardRecord/SRBeam.h"
#include "StandardRecord/SRTrueMichelE.h"
#include "StandardRecord/SRTrueParticle.h"
#include "StandardRecord/SRNuGenParticle.h"
#include "StandardRecord/SRVector3D.h"
#include "StandardRecord/SRLorentzVector.h"
#include "StandardRecord/SRMCReweight.h"
#include "StandardRecord/SREnums.h"
#include "StandardRecord/SRGeant4Weights.h"

namespace caf
{
  /// The SRNeutrino is a representation of neutrino interaction information
  class SRNeutrino
    {
    public:
      SRNeutrino();
      ~SRNeutrino() {  };

      short          pdg;           ///< pdg code
      float          E;             ///< True energy [GeV]
      float          visE;          ///< Sum of FLS hits that made CellHits from this neutrino [GeV]
      float          visEinslc;     ///< Sum of FLS hits that made CellHits  from this neutrino in this subevent [GeV]
      float          visENeutron;          ///< Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons
      float          visEinslcNeutron;     ///< Sum of FLS hits that made CellHits  from this neutrino in this subevent [GeV] that were daughters of neutrons
      float          visEBirks;          ///< Sum of FLS hits that made CellHits from this neutrino [GeV] with birks suppression
      float          visEinslcBirks;     ///< Sum of FLS hits that made CellHits  from this neutrino in this subevent [GeV] with birks suppression
      float          visENeutronBirks;          ///< Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons with birks suppression
      float          visEinslcNeutronBirks;     ///< Sum of FLS hits that made CellHits  from this neutrino in this subevent [GeV] that were daughters of neutrons with birks suppression
      float          L;             ///< True distance from hadron/muon decay to neutrino interaction [km]
      bool           isvtxcont;     ///< Checks if neutrino true vertex is within detector

      float          eff;           ///< Slicer efficiency for this truth interaction
      float          pur;           ///< Slicer purity for this truth interaction
      unsigned int   nhitslc;       ///< Number of hits recorded in this slice by this truth interaction
      unsigned int   nhittot;       ///< Total number of hits recorded for this truth interaction

      float          time;          ///< interaction time.
      SRLorentzVector p;             ///< True momentum [GeV]
      SRVector3D     vtx;           ///< Vertex position in detector coordinates [cm]
      std::vector<SRTrueMichelE> michel;  ///< Vector of true Michel electrons

      short          pdgorig;       ///< Unoscillated (unswapped) pdg code

      /// PDG code of struck nucleon (or, in the case of MEC, struck nucleon-nucleon pair).
      /// For MEC, the codes are: 2000000200 --> nn,  2000000201 --> np,  2000000202 --> "pp"
      int            hitnuc;
      SRLorentzVector hitnucp;      ///< Initial state 4-momentum of the struck nucleon 

      float          woscdumb;      ///< Simplest possible oscillation weight

      int            mode;          ///< interaction mode from enum mode_type::[QE, RES, COH, ...]
      bool           iscc;          ///< true if charged-current interaction, false if not.
      int            inttype;       ///< Interaction type enum int_type::[...]

      unsigned int   npiplus;       ///< Number of \f$\pi^+\f$'s after neutrino reaction, before FSI
      unsigned int   npiminus;      ///< Number of \f$\pi^1\f$'s after neutrino reaction, before FSI
      unsigned int   npizero;       ///< Number of \f$\pi^0\f$'s after neutrino reaction, before FSI
      unsigned int   nproton;       ///< Number of protons after neutrino reaction, before FSI
      unsigned int   nneutron;      ///< Number of neutrons after neutrino reaction, before FSI

      bool           ischarm;       ///< Is a charmed quark in interaction
      bool           isseaquark;    ///< Did neutrino scatter off a sea quark
      int            resnum;        ///< Straight from GENIE, resonance number

      float          xsec;          ///< xsec for thrown interaction, in 1/GeV^2, as stored by the GENIE spline
      int            tgtZ;         ///< Z of target nucleus
      int            tgtA;         ///< A of target nucleus

      float          q2;            ///< Squared momentum transfer [GeV^2]
      float          x;             ///< Bjorken x = (k-k')^2 / (2*p.q), [Dimensionless]
      float          y;             ///< Bjorken y =  (p.q) / (k.p), fractional energy loss of incoming particle [Dimensionless]
      float          W2;            ///< Invariant mass of final state squared. [GeV^2]

      float          genweight;     ///< Weight, if any, assigned by the generator

      std::vector<SRTrueParticle> prim;  ///< Primary daughters, lepton comes first in vector.
      std::vector<SRNuGenParticle> prefsi;  ///< Pre-FSI hadrons

      generator_                 generator;          ///< The generator that created this neutrino interaction
      std::vector<unsigned int>  genVersion;         ///< Version of the generator that created this neutrino interaction
      std::string                genConfigString;    ///< String associated with generator configuration.  (For GENIE 3, this is the "Comprehensive Model Configuration", sometimes a.k.a. "tune"

      SRMCReweight rwgt;

      SRBeam beam; ///< Information about neutrino production

      unsigned int   run;       ///< run number
      unsigned int   subrun;    ///< subrun number
      unsigned int   evt;       ///< ART event number, indexes trigger windows.
      int            cycle;     ///< MC repition index on a given run,subrun

      // Needed for flux systematics
      bool   isFHC;
      bool   is0HC;
      bool   isRHC;
      Det_t  det;

    };

} // end namespace

#endif // SRNEUTRINO_H
//////////////////////////////////////////////////////////////////////////////