SRNeutrino.h
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 // \author gsdavies@iastate.edu
3 // \date February 2013
4 ////////////////////////////////////////////////////////////////////////
5 #ifndef SRNEUTRINO_H
6 #define SRNEUTRINO_H
7 
15 #include "StandardRecord/SREnums.h"
17 
18 namespace caf
19 {
20  /// The SRNeutrino is a representation of neutrino interaction information
21  class SRNeutrino
22  {
23  public:
24  SRNeutrino();
25  ~SRNeutrino() { };
26 
27  short pdg; ///< pdg code
28  float E; ///< True energy [GeV]
29  float visE; ///< Sum of FLS hits that made CellHits from this neutrino [GeV]
30  float visEinslc; ///< Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV]
31  float visENeutron; ///< Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons
32  float visEinslcNeutron; ///< Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] that were daughters of neutrons
33  float visEBirks; ///< Sum of FLS hits that made CellHits from this neutrino [GeV] with birks suppression
34  float visEinslcBirks; ///< Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] with birks suppression
35  float visENeutronBirks; ///< Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons with birks suppression
36  float visEinslcNeutronBirks; ///< Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] that were daughters of neutrons with birks suppression
37  float L; ///< True distance from hadron/muon decay to neutrino interaction [km]
38  bool isvtxcont; ///< Checks if neutrino true vertex is within detector
39 
40  float eff; ///< Slicer efficiency for this truth interaction
41  float pur; ///< Slicer purity for this truth interaction
42  unsigned int nhitslc; ///< Number of hits recorded in this slice by this truth interaction
43  unsigned int nhittot; ///< Total number of hits recorded for this truth interaction
44 
45  float time; ///< interaction time.
46  SRLorentzVector p; ///< True momentum [GeV]
47  SRVector3D vtx; ///< Vertex position in detector coordinates [cm]
48  std::vector<SRTrueMichelE> michel; ///< Vector of true Michel electrons
49 
50  short pdgorig; ///< Unoscillated (unswapped) pdg code
51 
52  /// PDG code of struck nucleon (or, in the case of MEC, struck nucleon-nucleon pair).
53  /// For MEC, the codes are: 2000000200 --> nn, 2000000201 --> np, 2000000202 --> "pp"
54  int hitnuc;
55  SRLorentzVector hitnucp; ///< Initial state 4-momentum of the struck nucleon
56 
57  float woscdumb; ///< Simplest possible oscillation weight
58 
59  int mode; ///< interaction mode from enum mode_type::[QE, RES, COH, ...]
60  bool iscc; ///< true if charged-current interaction, false if not.
61  int inttype; ///< Interaction type enum int_type::[...]
62 
63  unsigned int npiplus; ///< Number of \f$\pi^+\f$'s after neutrino reaction, before FSI
64  unsigned int npiminus; ///< Number of \f$\pi^1\f$'s after neutrino reaction, before FSI
65  unsigned int npizero; ///< Number of \f$\pi^0\f$'s after neutrino reaction, before FSI
66  unsigned int nproton; ///< Number of protons after neutrino reaction, before FSI
67  unsigned int nneutron; ///< Number of neutrons after neutrino reaction, before FSI
68 
69  bool ischarm; ///< Is a charmed quark in interaction
70  bool isseaquark; ///< Did neutrino scatter off a sea quark
71  int resnum; ///< Straight from GENIE, resonance number
72 
73  float xsec; ///< xsec for thrown interaction, in 1/GeV^2, as stored by the GENIE spline
74  int tgtZ; ///< Z of target nucleus
75  int tgtA; ///< A of target nucleus
76 
77  float q2; ///< Squared momentum transfer [GeV^2]
78  float x; ///< Bjorken x = (k-k')^2 / (2*p.q), [Dimensionless]
79  float y; ///< Bjorken y = (p.q) / (k.p), fractional energy loss of incoming particle [Dimensionless]
80  float W2; ///< Invariant mass of final state squared. [GeV^2]
81 
82  float genweight; ///< Weight, if any, assigned by the generator
83 
84  std::vector<SRTrueParticle> prim; ///< Primary daughters, lepton comes first in vector.
85  std::vector<SRNuGenParticle> prefsi; ///< Pre-FSI hadrons
86 
87  generator_ generator; ///< The generator that created this neutrino interaction
88  std::vector<unsigned int> genVersion; ///< Version of the generator that created this neutrino interaction
89  std::string genConfigString; ///< String associated with generator configuration. (For GENIE 3, this is the "Comprehensive Model Configuration", sometimes a.k.a. "tune"
90 
92 
93  SRBeam beam; ///< Information about neutrino production
94 
95  unsigned int run; ///< run number
96  unsigned int subrun; ///< subrun number
97  unsigned int evt; ///< ART event number, indexes trigger windows.
98  int cycle; ///< MC repition index on a given run,subrun
99 
100  // Needed for flux systematics
101  bool isFHC;
102  bool is0HC;
103  bool isRHC;
105 
106  };
107 
108 } // end namespace
109 
110 #endif // SRNEUTRINO_H
111 //////////////////////////////////////////////////////////////////////////////
float x
Bjorken x = (k-k&#39;)^2 / (2*p.q), [Dimensionless].
Definition: SRNeutrino.h:78
Information about the neutrino production. Docs from http://www.hep.utexas.edu/~zarko/wwwgnumi/v19/v1...
Definition: SRBeam.h:14
float visEinslcNeutron
Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] that were daughters of n...
Definition: SRNeutrino.h:32
float visEinslcNeutronBirks
Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] that were daughters of n...
Definition: SRNeutrino.h:36
Det_t
Which NOvA detector?
Definition: SREnums.h:7
SRVector3D vtx
Vertex position in detector coordinates [cm].
Definition: SRNeutrino.h:47
float visE
Sum of FLS hits that made CellHits from this neutrino [GeV].
Definition: SRNeutrino.h:29
int cycle
MC repition index on a given run,subrun.
Definition: SRNeutrino.h:98
float visENeutronBirks
Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons with birk...
Definition: SRNeutrino.h:35
A 3-vector with more efficient storage than TVector3.
Definition: SRVector3D.h:14
float visEinslc
Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV].
Definition: SRNeutrino.h:30
unsigned int npizero
Number of &#39;s after neutrino reaction, before FSI.
Definition: SRNeutrino.h:65
SRBeam beam
Information about neutrino production.
Definition: SRNeutrino.h:93
float q2
Squared momentum transfer [GeV^2].
Definition: SRNeutrino.h:77
int inttype
Interaction type enum int_type::[...].
Definition: SRNeutrino.h:61
float woscdumb
Simplest possible oscillation weight.
Definition: SRNeutrino.h:57
std::vector< SRTrueParticle > prim
Primary daughters, lepton comes first in vector.
Definition: SRNeutrino.h:84
unsigned int npiplus
Number of &#39;s after neutrino reaction, before FSI.
Definition: SRNeutrino.h:63
SRMCReweight rwgt
Definition: SRNeutrino.h:91
std::vector< unsigned int > genVersion
Version of the generator that created this neutrino interaction.
Definition: SRNeutrino.h:88
float visENeutron
Sum of FLS hits that made CellHits from this neutrino [GeV] that were daughters of neutrons...
Definition: SRNeutrino.h:31
int tgtA
A of target nucleus.
Definition: SRNeutrino.h:75
unsigned int evt
ART event number, indexes trigger windows.
Definition: SRNeutrino.h:97
unsigned int nproton
Number of protons after neutrino reaction, before FSI.
Definition: SRNeutrino.h:66
std::string genConfigString
String associated with generator configuration. (For GENIE 3, this is the "Comprehensive Model Config...
Definition: SRNeutrino.h:89
generator_
Known generators of neutrino interactions (extend as other generators are used)
Definition: SREnums.h:28
float genweight
Weight, if any, assigned by the generator.
Definition: SRNeutrino.h:82
unsigned int subrun
subrun number
Definition: SRNeutrino.h:96
bool isseaquark
Did neutrino scatter off a sea quark.
Definition: SRNeutrino.h:70
The SRNeutrino is a representation of neutrino interaction information.
Definition: SRNeutrino.h:21
float pur
Slicer purity for this truth interaction.
Definition: SRNeutrino.h:41
SRLorentzVector p
True momentum [GeV].
Definition: SRNeutrino.h:46
unsigned int npiminus
Number of &#39;s after neutrino reaction, before FSI.
Definition: SRNeutrino.h:64
int resnum
Straight from GENIE, resonance number.
Definition: SRNeutrino.h:71
bool iscc
true if charged-current interaction, false if not.
Definition: SRNeutrino.h:60
float L
True distance from hadron/muon decay to neutrino interaction [km].
Definition: SRNeutrino.h:37
float W2
Invariant mass of final state squared. [GeV^2].
Definition: SRNeutrino.h:80
float eff
Slicer efficiency for this truth interaction.
Definition: SRNeutrino.h:40
unsigned int run
run number
Definition: SRNeutrino.h:95
unsigned int nhitslc
Number of hits recorded in this slice by this truth interaction.
Definition: SRNeutrino.h:42
short pdgorig
Unoscillated (unswapped) pdg code.
Definition: SRNeutrino.h:50
int mode
interaction mode from enum mode_type::[QE, RES, COH, ...]
Definition: SRNeutrino.h:59
4-vector with more efficient storage than TLorentzVector
int tgtZ
Z of target nucleus.
Definition: SRNeutrino.h:74
short pdg
pdg code
Definition: SRNeutrino.h:25
bool isvtxcont
Checks if neutrino true vertex is within detector.
Definition: SRNeutrino.h:38
float visEBirks
Sum of FLS hits that made CellHits from this neutrino [GeV] with birks suppression.
Definition: SRNeutrino.h:33
float time
interaction time.
Definition: SRNeutrino.h:45
float xsec
xsec for thrown interaction, in 1/GeV^2, as stored by the GENIE spline
Definition: SRNeutrino.h:73
std::vector< SRNuGenParticle > prefsi
Pre-FSI hadrons.
Definition: SRNeutrino.h:85
float visEinslcBirks
Sum of FLS hits that made CellHits from this neutrino in this subevent [GeV] with birks suppression...
Definition: SRNeutrino.h:34
float E
True energy [GeV].
Definition: SRNeutrino.h:28
This module creates Common Analysis Files.
Definition: FileReducer.h:10
std::vector< SRTrueMichelE > michel
Vector of true Michel electrons.
Definition: SRNeutrino.h:48
float y
Bjorken y = (p.q) / (k.p), fractional energy loss of incoming particle [Dimensionless].
Definition: SRNeutrino.h:79
bool ischarm
Is a charmed quark in interaction.
Definition: SRNeutrino.h:69
unsigned int nneutron
Number of neutrons after neutrino reaction, before FSI.
Definition: SRNeutrino.h:67
generator_ generator
The generator that created this neutrino interaction.
Definition: SRNeutrino.h:87
Various weights for systematic reweights of MC.
Definition: SRMCReweight.h:17
unsigned int nhittot
Total number of hits recorded for this truth interaction.
Definition: SRNeutrino.h:43
enum BeamMode string
SRLorentzVector hitnucp
Initial state 4-momentum of the struck nucleon.
Definition: SRNeutrino.h:55