StructsDUNE.h File Reference

#include <iostream>
#include "TVector3.h"
#include "TLorentzVector.h"

Include dependency graph for StructsDUNE.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	dunemc_base

Enumerations
enum	MaCh3_Mode {   kMaCh3_CCQE = 0 , kMaCh3_CC_Single_Kaon = 1 , kMaCh3_CC_DIS = 2 , kMaCh3_CC_RES = 3 ,   kMaCh3_CC_COH = 4 , kMaCh3_CC_Diffractive = 5 , kMaCh3_CC_Nue_EL = 6 , kMaCh3_CC_IMD = 7 ,   kMaCh3_CC_AMnuGamma = 8 , kMaCh3_CC_MEC = 9 , kMaCh3_CC_COHEL = 10 , kMaCh3_CC_IBD = 11 ,   kMaCh3_CC_GlashowRES = 12 , kMaCh3_CC_IMDAnnihalation = 13 , kMaCh3_NCQE = 14 , kMaCh3_NC_DIS = 15 ,   kMaCh3_NC_RES = 16 , kMaCh3_NC_COH = 17 , kMaCh3_NC_Diffractive = 18 , kMaCh3_NC_Nue_EL = 19 ,   kMaCh3_NC_IMD = 20 , kMaCh3_NC_AMnuGamma = 21 , kMaCh3_NC_MEC = 22 , kMaCh3_NC_COHEL = 23 ,   kMaCh3_NC_IBD = 24 , kMaCh3_NC_GlashowRES = 25 , kMaCh3_NC_IMDAnnihalation = 26 , kMaCh3_nModes }
enum	SIMB_Mode {   kUnknownInteraction = -1 , kQE = 0 , kRes = 1 , kDIS = 2 ,   kCoh = 3 , kCohElastic = 4 , kElectronScattering = 5 , kIMDAnnihilation = 6 ,   kInverseBetaDecay = 7 , kGlashowResonance = 8 , kAMNuGamma = 9 , kMEC = 10 ,   kDiffractive = 11 , kEM = 12 , kWeakMix = 13 , kSIMB_nModes }
enum	GENIE_Mode {   gUnknownInteraction = -1 , gQE = 1 , gSingleKaon = 2 , gDIS = 3 ,   gRes = 4 , gCoh = 5 , gDiffractive = 6 , gElectronScattering = 7 ,   gIMD = 8 , gAMNuGamma = 9 , gMEC = 10 , gCohElastic = 11 ,   gInverseBetaDecay = 12 , gGlashowResonance = 13 , gIMDAnnihilation = 14 , gGENIE_nModes }
enum	MaCh3_Spline_Modes {   kMaCh3_Spline_CCQE = 0 , kMaCh3_Spline_CC_Single_Kaon = 1 , kMaCh3_Spline_CC_DIS = 2 , kMaCh3_Spline_CC_RES = 3 ,   kMaCh3_Spline_CC_COH = 4 , kMaCh3_Spline_CC_Diffractive = 5 , kMaCh3_Spline_CC_Nue_EL = 6 , kMaCh3_Spline_CC_IMD = 7 ,   kMaCh3_Spline_CC_AMnuGamma = 8 , kMaCh3_Spline_CC_MEC = 9 , kMaCh3_Spline_CC_COHEL = 10 , kMaCh3_Spline_CC_IBD = 11 ,   kMaCh3_Spline_CC_GlashowRES = 12 , kMaCh3_Spline_CC_IMDAnnihalation = 13 , kMaCh3_Spline_NCQE = 14 , kMaCh3_Spline_NC_DIS = 15 ,   kMaCh3_Spline_NC_RES = 16 , kMaCh3_Spline_NC_COH = 17 , kMaCh3_Spline_NC_Diffractive = 18 , kMaCh3_Spline_NC_Nue_EL = 19 ,   kMaCh3_Spline_NC_IMD = 20 , kMaCh3_Spline_NC_AMnuGamma = 21 , kMaCh3_Spline_NC_MEC = 22 , kMaCh3_Spline_NC_COHEL = 23 ,   kMaCh3_Spline_NC_IBD = 24 , kMaCh3_Spline_NC_GlashowRES = 25 , kMaCh3_Spline_NC_IMDAnnihalation = 26 , kMaCh3_Spline_nModes }

Functions
void	TotalEScaleND (const double *par, double *erec, double erecHad, double erecLep, bool NotCCnumu)
void	TotalEScaleSqrtND (const double *par, double *erec, double erecHad, double erecLep, double sqrtErecHad, double sqrtErecLep, bool NotCCnumu)
void	TotalEScaleInvSqrtND (const double *par, double *erec, double erecHad, double erecLep, double invSqrtErecHad, double invSqrtErecLep, bool NotCCnumu)
void	HadEScaleND (const double *par, double *erec, double sumEhad)
void	HadEScaleSqrtND (const double *par, double *erec, double sumEhad, double sqrtSumEhad)
void	HadEScaleInvSqrtND (const double *par, double *erec, double sumEhad, double invSqrtSumEhad)
void	MuEScaleND (const double *par, double *erec, double erecLep, bool CCnumu)
void	MuEScaleSqrtND (const double *par, double *erec, double erecLep, double sqrtErecLep, bool CCnumu)
void	MuEScaleInvSqrtND (const double *par, double *erec, double erecLep, double invSqrtErecLep, bool CCnumu)
void	NEScaleND (const double *par, double *erec, double eRecoN)
void	NEScaleSqrtND (const double *par, double *erec, double eRecoN, double sqrteRecoN)
void	NEScaleInvSqrtND (const double *par, double *erec, double eRecoN, double invSqrteRecoN)
void	EMEScaleND (const double *par, double *erec, double eRecoPi0, double erecLep, bool CCnue)
void	EMEScaleSqrtND (const double *par, double *erec, double eRecoPi0, double erecLep, double sqrtErecLep, double sqrteRecoPi0, bool CCnue)
void	EMEScaleInvSqrtND (const double *par, double *erec, double eRecoPi0, double erecLep, double invSqrtErecLep, double invSqrteRecoPi0, bool CCnue)
void	HadResND (const double *par, double *erec, double eRecoP, double eRecoPip, double eRecoPim, double eP, double ePip, double ePim)
void	MuResND (const double *par, double *erec, double erecLep, double LepE, bool CCnumu)
void	NResND (const double *par, double *erec, double eRecoN, double eN)
void	EMResND (const double *par, double *erec, double eRecoPi0, double ePi0, double erecLep, double LepE, bool CCnue)
void	TotalEScaleFD (const double *par, double *erec, double erecHad, double erecLep, bool NotCCnumu)
void	TotalEScaleSqrtFD (const double *par, double *erec, double erecHad, double erecLep, double sqrtErecHad, double sqrtErecLep, bool NotCCnumu)
void	TotalEScaleInvSqrtFD (const double *par, double *erec, double erecHad, double erecLep, double invSqrtErecHad, double invSqrtErecLep, bool NotCCnumu)
void	HadEScaleFD (const double *par, double *erec, double sumEhad)
void	HadEScaleSqrtFD (const double *par, double *erec, double sumEhad, double sqrtSumEhad)
void	HadEScaleInvSqrtFD (const double *par, double *erec, double sumEhad, double invSqrtSumEhad)
void	MuEScaleFD (const double *par, double *erec, double erecLep, bool CCnumu)
void	MuEScaleSqrtFD (const double *par, double *erec, double erecLep, double sqrtErecLep, bool CCnumu)
void	MuEScaleInvSqrtFD (const double *par, double *erec, double erecLep, double invSqrtErecLep, bool CCnumu)
void	NEScaleFD (const double *par, double *erec, double eRecoN)
void	NEScaleSqrtFD (const double *par, double *erec, double eRecoN, double sqrteRecoN)
void	NEScaleInvSqrtFD (const double *par, double *erec, double eRecoN, double invSqrteRecoN)
void	EMEScaleFD (const double *par, double *erec, double eRecoPi0, double erecLep, bool CCnue)
void	EMEScaleSqrtFD (const double *par, double *erec, double eRecoPi0, double erecLep, double sqrtErecLep, double sqrteRecoPi0, bool CCnue)
void	EMEScaleInvSqrtFD (const double *par, double *erec, double eRecoPi0, double erecLep, double invSqrtErecLep, double invSqrteRecoPi0, bool CCnue)
void	HadResFD (const double *par, double *erec, double eRecoP, double eRecoPip, double eRecoPim, double eP, double ePip, double ePim)
void	MuResFD (const double *par, double *erec, double erecLep, double LepE, bool CCnumu)
void	NResFD (const double *par, double *erec, double eRecoN, double eN)
void	EMResFD (const double *par, double *erec, double eRecoPi0, double ePi0, double erecLep, double LepE, bool CCnue)
void	CVNNumuFD (const double *par, double *cvnnumu)
void	CVNNueFD (const double *par, double *cvnnue)
int	SIMBMode_ToMaCh3Mode (int SIMB_mode, int isCC)
int	GENIEMode_ToMaCh3Mode (int GENIE_mode, int isCC)
std::string	MaCh3mode_ToDUNEString (MaCh3_Mode i)
int	MaCh3Mode_to_SplineMode (int iMode)

Enumeration Type Documentation

GENIE_Mode

enum GENIE_Mode

Enumerator
gUnknownInteraction
gQE
gSingleKaon
gDIS
gRes
gCoh
gDiffractive
gElectronScattering
gIMD
gAMNuGamma
gMEC
gCohElastic
gInverseBetaDecay
gGlashowResonance
gIMDAnnihilation
gGENIE_nModes

Definition at line 724 of file StructsDUNE.h.

                {
  // ***************************
 
  // Unknown
  gUnknownInteraction        =   -1,
  // QE
  gQE                        =    1,
  // Single Kaon
  gSingleKaon                =    2,
  // DIS
  gDIS                       =    3,
  // RES
  gRes                       =    4,
  // Coherent
  gCoh                       =    5,
  // Diffractive
  gDiffractive                =    6,
  // Nu-e Elastic
  gElectronScattering        =    7, 
  // Inverse Muon Decay
  gIMD                       =    8,
  // Atmospheric Muon Nu Gamma
  gAMNuGamma                 =    9,
  // MEC aka 2p2h
  gMEC                       =   10,
  // Coherent Elastic
  gCohElastic                =   11,
  // Inverse Beta Decay
  gInverseBetaDecay          =   12, 
  // Glasgow Resonance
  gGlashowResonance          =   13,
  // Inverse Muon Decay Annihliation
  gIMDAnnihilation           =   14,
  // Just keep a counter of the number of modes
  gGENIE_nModes
};

MaCh3_Mode

enum MaCh3_Mode

Enumerator
kMaCh3_CCQE
kMaCh3_CC_Single_Kaon
kMaCh3_CC_DIS
kMaCh3_CC_RES
kMaCh3_CC_COH
kMaCh3_CC_Diffractive
kMaCh3_CC_Nue_EL
kMaCh3_CC_IMD
kMaCh3_CC_AMnuGamma
kMaCh3_CC_MEC
kMaCh3_CC_COHEL
kMaCh3_CC_IBD
kMaCh3_CC_GlashowRES
kMaCh3_CC_IMDAnnihalation
kMaCh3_NCQE
kMaCh3_NC_DIS
kMaCh3_NC_RES
kMaCh3_NC_COH
kMaCh3_NC_Diffractive
kMaCh3_NC_Nue_EL
kMaCh3_NC_IMD
kMaCh3_NC_AMnuGamma
kMaCh3_NC_MEC
kMaCh3_NC_COHEL
kMaCh3_NC_IBD
kMaCh3_NC_GlashowRES
kMaCh3_NC_IMDAnnihalation
kMaCh3_nModes

Definition at line 617 of file StructsDUNE.h.

                {
  // ***************************
 
  // CCQE
  kMaCh3_CCQE              = 0,
  // 1 Kaon
  kMaCh3_CC_Single_Kaon     = 1,
  // CC DIS
  kMaCh3_CC_DIS            = 2,
  // CC RES
  kMaCh3_CC_RES             = 3,
  // CC COH
  kMaCh3_CC_COH             = 4,
  // CC Diffractive
  kMaCh3_CC_Diffractive     = 5,
  // CC Electron EL
  kMaCh3_CC_Nue_EL          = 6,
  // CC Inverse Muon Decay
  kMaCh3_CC_IMD             = 7,
  // CC Atmospheric Neutrino Gamma
  kMaCh3_CC_AMnuGamma       = 8,
  // CC MEC (aka 2p2h)
  kMaCh3_CC_MEC             = 9,
  // CC Coherent Elastic
  kMaCh3_CC_COHEL           = 10,
  // CC Inverse Beta Decay
  kMaCh3_CC_IBD             = 11,
  // CC GlashowRES 
  kMaCh3_CC_GlashowRES      = 12,
  // CC IMD Annihilation
  kMaCh3_CC_IMDAnnihalation = 13,
  
  // NCQE
  kMaCh3_NCQE              = 14,
  // NC DIS
  kMaCh3_NC_DIS            = 15,
  // NC RES
  kMaCh3_NC_RES             = 16,
  // NC COH
  kMaCh3_NC_COH             = 17,
  // CC Diffractive
  kMaCh3_NC_Diffractive     = 18,
  // NC Electron EL
  kMaCh3_NC_Nue_EL          = 19,
  // NC Inverse Muon Decay
  kMaCh3_NC_IMD             = 20,
  // NC Atmospheric Neutrino Gamma
  kMaCh3_NC_AMnuGamma       = 21,
  // NC MEC (aka 2p2h)
  kMaCh3_NC_MEC             = 22,
  // NC Coherent Elastic
  kMaCh3_NC_COHEL           = 23,
  // NC Inverse Beta Decay
  kMaCh3_NC_IBD             = 24,
  // NC GlashowRES 
  kMaCh3_NC_GlashowRES      = 25,
  // NC IMD Annihilation
  kMaCh3_NC_IMDAnnihalation = 26,
 
  // Keep a counter of the number of MaCh3 modes we have
  kMaCh3_nModes
};

MaCh3_Spline_Modes

enum MaCh3_Spline_Modes

Enumerator
kMaCh3_Spline_CCQE
kMaCh3_Spline_CC_Single_Kaon
kMaCh3_Spline_CC_DIS
kMaCh3_Spline_CC_RES
kMaCh3_Spline_CC_COH
kMaCh3_Spline_CC_Diffractive
kMaCh3_Spline_CC_Nue_EL
kMaCh3_Spline_CC_IMD
kMaCh3_Spline_CC_AMnuGamma
kMaCh3_Spline_CC_MEC
kMaCh3_Spline_CC_COHEL
kMaCh3_Spline_CC_IBD
kMaCh3_Spline_CC_GlashowRES
kMaCh3_Spline_CC_IMDAnnihalation
kMaCh3_Spline_NCQE
kMaCh3_Spline_NC_DIS
kMaCh3_Spline_NC_RES
kMaCh3_Spline_NC_COH
kMaCh3_Spline_NC_Diffractive
kMaCh3_Spline_NC_Nue_EL
kMaCh3_Spline_NC_IMD
kMaCh3_Spline_NC_AMnuGamma
kMaCh3_Spline_NC_MEC
kMaCh3_Spline_NC_COHEL
kMaCh3_Spline_NC_IBD
kMaCh3_Spline_NC_GlashowRES
kMaCh3_Spline_NC_IMDAnnihalation
kMaCh3_Spline_nModes

Definition at line 1138 of file StructsDUNE.h.

                        {
 
  // ***************************
  // CCQE
  kMaCh3_Spline_CCQE              = 0,
  // 1 Kaon
  kMaCh3_Spline_CC_Single_Kaon     = 1,
  // CC DIS
  kMaCh3_Spline_CC_DIS            = 2,
  // CC RES
  kMaCh3_Spline_CC_RES             = 3,
  // CC COH
  kMaCh3_Spline_CC_COH             = 4,
  // CC Diffractive
  kMaCh3_Spline_CC_Diffractive     = 5,
  // CC Electron EL
  kMaCh3_Spline_CC_Nue_EL          = 6,
  // CC Inverse Muon Decay
  kMaCh3_Spline_CC_IMD             = 7,
  // CC Atmospheric Neutrino Gamma
  kMaCh3_Spline_CC_AMnuGamma       = 8,
  // CC MEC (aka 2p2h)
  kMaCh3_Spline_CC_MEC             = 9,
  // CC Coherent Elastic
  kMaCh3_Spline_CC_COHEL           = 10,
  // CC Inverse Beta Decay
  kMaCh3_Spline_CC_IBD             = 11,
  // CC GlashowRES 
  kMaCh3_Spline_CC_GlashowRES      = 12,
  // CC IMD Annihilation
  kMaCh3_Spline_CC_IMDAnnihalation = 13,
  
  // NCQE
  kMaCh3_Spline_NCQE              = 14,
  // NC DIS
  kMaCh3_Spline_NC_DIS            = 15,
  // NC RES
  kMaCh3_Spline_NC_RES             = 16,
  // NC COH
  kMaCh3_Spline_NC_COH             = 17,
  // CC Diffractive
  kMaCh3_Spline_NC_Diffractive     = 18,
  // NC Electron EL
  kMaCh3_Spline_NC_Nue_EL          = 19,
  // NC Inverse Muon Decay
  kMaCh3_Spline_NC_IMD             = 20,
  // NC Atmospheric Neutrino Gamma
  kMaCh3_Spline_NC_AMnuGamma       = 21,
  // NC MEC (aka 2p2h)
  kMaCh3_Spline_NC_MEC             = 22,
  // NC Coherent Elastic
  kMaCh3_Spline_NC_COHEL           = 23,
  // NC Inverse Beta Decay
  kMaCh3_Spline_NC_IBD             = 24,
  // NC GlashowRES 
  kMaCh3_Spline_NC_GlashowRES      = 25,
  // NC IMD Annihilation
  kMaCh3_Spline_NC_IMDAnnihalation = 26,
 
  // Keep a counter of the number of MaCh3 modes we have
  kMaCh3_Spline_nModes
 
};

SIMB_Mode

enum SIMB_Mode

Enumerator
kUnknownInteraction
kQE
kRes
kDIS
kCoh
kCohElastic
kElectronScattering
kIMDAnnihilation
kInverseBetaDecay
kGlashowResonance
kAMNuGamma
kMEC
kDiffractive
kEM
kWeakMix
kSIMB_nModes

Definition at line 684 of file StructsDUNE.h.

               {
  // ***************************
 
  // Unknown
  kUnknownInteraction        =   -1,
  // QE
  kQE                        =    0,
  // Resonant
  kRes                       =    1,
  // DIS
  kDIS                       =    2,
  // Coherent
  kCoh                       =    3,
  // Coherent Elastic
  kCohElastic                =    4,
  // Electron Scattering
  kElectronScattering        =    5,
  // Inverse Muon Decay Annihliation
  kIMDAnnihilation           =    6,
  // Inverse Beta Decay
  kInverseBetaDecay          =    7, 
  // Glasgow Resonance
  kGlashowResonance          =    8,
  // Atmospheric Muon Nu Gamma
  kAMNuGamma                 =    9,
  // MEC aka 2p2h
  kMEC                       =   10,
  // Diffractive
  kDiffractive               =   11,
  //
  kEM                        =   12,
  //
  kWeakMix                   =   13,
  // Just keep a counter of the number of modes
  kSIMB_nModes
};

Function Documentation

CVNNueFD()

void CVNNueFD ( const double * par, double * cvnnue )

inline

Definition at line 608 of file StructsDUNE.h.

                                                          {
 
  (*cvnnue) += (*par);
 
}

CVNNumuFD()

void CVNNumuFD ( const double * par, double * cvnnumu )

inline

Definition at line 601 of file StructsDUNE.h.

                                                            {
 
  (*cvnnumu) += (*par);
 
}

EMEScaleFD()

void EMEScaleFD ( const double * par, double * erec, double eRecoPi0, double erecLep, bool CCnue )

inline

Definition at line 500 of file StructsDUNE.h.

                                                                                                       {
 
  (*erec) += (*par) * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * erecLep;
  }
 
}

EMEScaleInvSqrtFD()

void EMEScaleInvSqrtFD ( const double * par, double * erec, double eRecoPi0, double erecLep, double invSqrtErecLep, double invSqrteRecoPi0, bool CCnue )

inline

Definition at line 526 of file StructsDUNE.h.

                                                                                                                                                             {
 
  (*erec) += (*par) * invSqrteRecoPi0 * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
 
}

EMEScaleInvSqrtND()

void EMEScaleInvSqrtND ( const double * par, double * erec, double eRecoPi0, double erecLep, double invSqrtErecLep, double invSqrteRecoPi0, bool CCnue )

inline

Definition at line 276 of file StructsDUNE.h.

                                                                                                                                                             {
 
  (*erec) += (*par) * invSqrteRecoPi0 * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
 
}

EMEScaleND()

void EMEScaleND ( const double * par, double * erec, double eRecoPi0, double erecLep, bool CCnue )

inline

Definition at line 250 of file StructsDUNE.h.

                                                                                                       {
 
  (*erec) += (*par) * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * erecLep;
  }
 
}

EMEScaleSqrtFD()

void EMEScaleSqrtFD ( const double * par, double * erec, double eRecoPi0, double erecLep, double sqrtErecLep, double sqrteRecoPi0, bool CCnue )

inline

Definition at line 513 of file StructsDUNE.h.

                                                                                                                                                    {
 
  (*erec) += (*par) * sqrteRecoPi0 * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
 
}

EMEScaleSqrtND()

void EMEScaleSqrtND ( const double * par, double * erec, double eRecoPi0, double erecLep, double sqrtErecLep, double sqrteRecoPi0, bool CCnue )

inline

Definition at line 263 of file StructsDUNE.h.

                                                                                                                                                    {
 
  (*erec) += (*par) * sqrteRecoPi0 * eRecoPi0;
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
 
}

EMResFD()

void EMResFD ( const double * par, double * erec, double eRecoPi0, double ePi0, double erecLep, double LepE, bool CCnue )

inline

Definition at line 584 of file StructsDUNE.h.

                                                                                                                              {
 
  (*erec) += (*par) * (ePi0 - eRecoPi0);
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * (LepE - erecLep);
  }
 
}

EMResND()

void EMResND ( const double * par, double * erec, double eRecoPi0, double ePi0, double erecLep, double LepE, bool CCnue )

inline

Definition at line 334 of file StructsDUNE.h.

                                                                                                                              {
 
  (*erec) += (*par) * (ePi0 - eRecoPi0);
 
  //if true CC nue AND reco nue event
  if (CCnue)
  {
    (*erec) += (*par) * (LepE - erecLep);
  }
 
}

GENIEMode_ToMaCh3Mode()

int GENIEMode_ToMaCh3Mode ( int GENIE_mode, int isCC )

inline

Definition at line 896 of file StructsDUNE.h.

                                                           {
  // **********************
 
  int ReturnMode = kMaCh3_nModes;
  
  if (isCC == 1) {
    switch (GENIE_mode) {
        //Unknown
      case gUnknownInteraction:
        ReturnMode = kMaCh3_nModes;
            break;
        //QE
      case gQE:
        ReturnMode = kMaCh3_CCQE; // CC QE in MaCh3 DUNE
        break;
        // DIS
      case gDIS:
        ReturnMode = kMaCh3_CC_DIS; //CC DIS in MaCh3 DUNE
        break;
        // RES
      case gRes:
        ReturnMode = kMaCh3_CC_RES; // CC RES in MaCh3 DUNE
        break;
        // Coherent
      case gCoh:
        ReturnMode = kMaCh3_CC_COH; // CC Coh in MaCh3 DUNE
        break;
        // Diffractive
      case gDiffractive:
        ReturnMode = kMaCh3_CC_Diffractive; // CC multi-pi in MaCh3
        break;
        // Nue Elastic
      case gElectronScattering:
        ReturnMode = kMaCh3_CC_Nue_EL; // CC Nue scattering in MaCh3 DUNE
        break;
        // Atmospheric Mu Gamma
      case gAMNuGamma:
        ReturnMode = kMaCh3_CC_AMnuGamma; // CC Am Nu Mu in MaCh3 DUNE
        break;
        // MEC
      case gMEC:
        ReturnMode = kMaCh3_CC_MEC; // CC MEC in MaCh3 DUNE
        break;
        // Coherent Elastic
      case gCohElastic:
        ReturnMode = kMaCh3_CC_COHEL; // CC Coherent Elastic in MaCh3 DUNE
        break;
        // Inverse Beta Decay
      case gInverseBetaDecay:
        ReturnMode = kMaCh3_CC_IBD; // CC Inverse Beta Decay in MaCh3 DUNE
        break;
        // Glashow Resonance
      case gGlashowResonance:
        ReturnMode = kMaCh3_CC_GlashowRES; // NC 1 gamma
        break;
        // IMD Annihilation
      case gIMDAnnihilation:
        ReturnMode = kMaCh3_CC_IMDAnnihalation; // NC other in MaCh3
        break;
      case gIMD:
        ReturnMode = kMaCh3_CC_IMDAnnihalation; // Stick Inverse Muon Decay into above
        break;
      case gSingleKaon:
        ReturnMode = kMaCh3_CC_IMDAnnihalation; // Stick Single Kaon into above
        break;
      default:
        ReturnMode = kMaCh3_nModes; // Something else in MaCh3 (sand?)
        break;
    }
 
  }
 
 
  if (isCC == 0) {
    switch (GENIE_mode) {
        //Unknown
      case kUnknownInteraction:
        ReturnMode = kMaCh3_nModes;
            break;
        //QE
      case gQE:
        ReturnMode = kMaCh3_NCQE; // NC QE in MaCh3 DUNE
        break;
        // DIS
      case gDIS:
        ReturnMode = kMaCh3_NC_DIS; // NC DIS in MaCh3 DUNE
        break;
        // RES
      case gRes:
        ReturnMode = kMaCh3_NC_RES; // NC RES in MaCh3 DUNE
        break;
        // Coherent
      case gCoh:
        ReturnMode = kMaCh3_NC_COH; // NC Coh in MaCh3 DUNE
        break;
        // Diffractive
      case gDiffractive:
        ReturnMode = kMaCh3_NC_Diffractive; // CC multi-pi in MaCh3
        break;
        // Nue Elastic
      case gElectronScattering:
        ReturnMode = kMaCh3_NC_Nue_EL; // NC Nue scattering in MaCh3 DUNE
        break;
        // Atmospheric Mu Gamma
      case gAMNuGamma:
        ReturnMode = kMaCh3_NC_AMnuGamma; // NC Am Nu Mu in MaCh3 DUNE
        break;
        // MEC
      case gMEC:
        ReturnMode = kMaCh3_NC_MEC; // NC MEC in MaCh3 DUNE
        break;
        // Coherent Elastic
      case gCohElastic:
        ReturnMode = kMaCh3_NC_COHEL; // NC Coherent Elastic in MaCh3 DUNE
        break;
        // Inverse Beta Decay
      case gInverseBetaDecay:
        ReturnMode = kMaCh3_NC_IBD; // Inverse Beta Decay in MaCh3 DUNE
        break;
        // Glashow Resonance
      case gGlashowResonance:
        ReturnMode = kMaCh3_NC_GlashowRES;
        break;
        // IMD Annihilation
      case gIMDAnnihilation:
        ReturnMode = kMaCh3_NC_IMDAnnihalation;
        break;
      case gIMD:
        ReturnMode = kMaCh3_NC_IMDAnnihalation; 
        break;
      case gSingleKaon:
        ReturnMode = kMaCh3_NC_IMDAnnihalation; 
        break;
      default:
        ReturnMode = kMaCh3_nModes; // Something else in MaCh3 (sand?)
        break;
    }
 
  }
 
  return ReturnMode;
};

HadEScaleFD()

void HadEScaleFD ( const double * par, double * erec, double sumEhad )

inline

Definition at line 406 of file StructsDUNE.h.

                                                                           {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * sumEhad;
  
}

HadEScaleInvSqrtFD()

void HadEScaleInvSqrtFD ( const double * par, double * erec, double sumEhad, double invSqrtSumEhad )

inline

Definition at line 422 of file StructsDUNE.h.

                                                                                                         {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * invSqrtSumEhad * sumEhad;
  
}

HadEScaleInvSqrtND()

void HadEScaleInvSqrtND ( const double * par, double * erec, double sumEhad, double invSqrtSumEhad )

inline

Definition at line 172 of file StructsDUNE.h.

                                                                                                         {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * invSqrtSumEhad * sumEhad;
  
}

HadEScaleND()

void HadEScaleND ( const double * par, double * erec, double sumEhad )

inline

Definition at line 156 of file StructsDUNE.h.

                                                                           {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * sumEhad;
  
}

HadEScaleSqrtFD()

void HadEScaleSqrtFD ( const double * par, double * erec, double sumEhad, double sqrtSumEhad )

inline

Definition at line 414 of file StructsDUNE.h.

                                                                                                   {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * sqrtSumEhad * sumEhad;
  
}

HadEScaleSqrtND()

void HadEScaleSqrtND ( const double * par, double * erec, double sumEhad, double sqrtSumEhad )

inline

Definition at line 164 of file StructsDUNE.h.

                                                                                                   {
 
  // Protons + Positive Pions + Negative Pions
  (*erec) += (*par) * sqrtSumEhad * sumEhad;
  
}

HadResFD()

void HadResFD ( const double * par, double * erec, double eRecoP, double eRecoPip, double eRecoPim, double eP, double ePip, double ePim )

inline

Definition at line 545 of file StructsDUNE.h.

                                                                                                                                              {
 
  // Reco Sum: Protons + Positive Pions + Negative Pions
  double recoSum = eRecoP + eRecoPip + eRecoPim;
 
  // True Sum: Protons + Positive Pions + Negative Pions
  double trueSum = eP + ePip + ePim;
 
  (*erec) += (*par) * (trueSum - recoSum);
 
}

HadResND()

void HadResND ( const double * par, double * erec, double eRecoP, double eRecoPip, double eRecoPim, double eP, double ePip, double ePim )

inline

Definition at line 295 of file StructsDUNE.h.

                                                                                                                                              {
 
  // Reco Sum: Protons + Positive Pions + Negative Pions
  double recoSum = eRecoP + eRecoPip + eRecoPim;
 
  // True Sum: Protons + Positive Pions + Negative Pions
  double trueSum = eP + ePip + ePim;
 
  (*erec) += (*par) * (trueSum - recoSum);
 
}

MaCh3Mode_to_SplineMode()

int MaCh3Mode_to_SplineMode ( int iMode )

inline

Definition at line 1202 of file StructsDUNE.h.

                                             {
  //No grouping of modes in MaCh3
  return iMode;
}

MaCh3mode_ToDUNEString()

std::string MaCh3mode_ToDUNEString ( MaCh3_Mode i )

inline

Definition at line 1041 of file StructsDUNE.h.

                                                      {
  // *******************************
  std::string name;
 
  switch(i) {
  case kMaCh3_CCQE:
    name = "ccqe";
    break;
  case kMaCh3_CC_Single_Kaon:
    name = "unknown";
    break;
  case kMaCh3_CC_DIS:
    name = "ccdis";
    break;
  case kMaCh3_CC_RES:
    name = "ccres";
    break;
  case kMaCh3_CC_COH:
    name = "cccoh";
    break;
  case kMaCh3_CC_COHEL:
    name = "cccohel";
    break;
  case kMaCh3_CC_Diffractive:
    name = "ccdiff";
    break;
  case kMaCh3_CC_Nue_EL:
    name = "ccnueel";
    break;
  case kMaCh3_CC_MEC:
    name = "ccmec";
    break;
  case kMaCh3_CC_IMD:
    name = "ccIMD";
    break;
  case kMaCh3_CC_GlashowRES:
    name = "ccglasres";
    break;
  case kMaCh3_CC_IMDAnnihalation:
    name = "ccimdannihilation";
    break;
  case kMaCh3_CC_AMnuGamma:
    name = "ccamnugamma";
    break;
  case kMaCh3_CC_IBD:
    name = "ccibd";
    break;
  case kMaCh3_NCQE:
    name = "ncqe";
    break;
  case kMaCh3_NC_DIS:
    name = "ncdis";
    break;
  case kMaCh3_NC_RES:
    name = "ncres";
    break;
  case kMaCh3_NC_COH:
    name = "nccoh";
    break;
  case kMaCh3_NC_COHEL:
    name = "nccohel";
    break;
  case kMaCh3_NC_Diffractive:
    name = "ncdiff";
    break;
  case kMaCh3_NC_Nue_EL:
    name = "ncnueel";
    break;
  case kMaCh3_NC_MEC:
    name = "ncmec";
    break;
  case kMaCh3_NC_IMD:
    name = "ncIMD";
    break;
  case kMaCh3_NC_GlashowRES:
    name = "ncglasres";
    break;
  case kMaCh3_NC_IMDAnnihalation:
    name = "ncimdannihilation";
    break;
  case kMaCh3_NC_AMnuGamma:
    name = "ncamnugamma";
    break;
  case kMaCh3_NC_IBD:
    name = "ncibd";
    break;
  case kMaCh3_nModes:
    name = "unknown";
    break;
  default:
    std::cerr << "Did not find the MaCh3 mode you specified " << i << std::endl;
    name = "UNKNOWN_BAD";
  }
  
  return name;
}

MuEScaleFD()

void MuEScaleFD ( const double * par, double * erec, double erecLep, bool CCnumu )

inline

Definition at line 436 of file StructsDUNE.h.

                                                                                       {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * erecLep;
  }
  
}

MuEScaleInvSqrtFD()

void MuEScaleInvSqrtFD ( const double * par, double * erec, double erecLep, double invSqrtErecLep, bool CCnumu )

inline

Definition at line 458 of file StructsDUNE.h.

                                                                                                                     {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
  
}

MuEScaleInvSqrtND()

void MuEScaleInvSqrtND ( const double * par, double * erec, double erecLep, double invSqrtErecLep, bool CCnumu )

inline

Definition at line 208 of file StructsDUNE.h.

                                                                                                                     {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
  
}

MuEScaleND()

void MuEScaleND ( const double * par, double * erec, double erecLep, bool CCnumu )

inline

Definition at line 186 of file StructsDUNE.h.

                                                                                       {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * erecLep;
  }
  
}

MuEScaleSqrtFD()

void MuEScaleSqrtFD ( const double * par, double * erec, double erecLep, double sqrtErecLep, bool CCnumu )

inline

Definition at line 447 of file StructsDUNE.h.

                                                                                                               {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
  
}

MuEScaleSqrtND()

void MuEScaleSqrtND ( const double * par, double * erec, double erecLep, double sqrtErecLep, bool CCnumu )

inline

Definition at line 197 of file StructsDUNE.h.

                                                                                                               {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
  
}

MuResFD()

void MuResFD ( const double * par, double * erec, double erecLep, double LepE, bool CCnumu )

inline

Definition at line 560 of file StructsDUNE.h.

                                                                                                 {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * (LepE - erecLep);
  }
 
}

MuResND()

void MuResND ( const double * par, double * erec, double erecLep, double LepE, bool CCnumu )

inline

Definition at line 310 of file StructsDUNE.h.

                                                                                                 {
 
  //if true CC numu AND reco numu event
  if (CCnumu)
  {
    (*erec) += (*par) * (LepE - erecLep);
  }
 
}

NEScaleFD()

void NEScaleFD ( const double * par, double * erec, double eRecoN )

inline

Definition at line 474 of file StructsDUNE.h.

                                                                        {
 
  (*erec) += (*par) * eRecoN;
  
}

NEScaleInvSqrtFD()

void NEScaleInvSqrtFD ( const double * par, double * erec, double eRecoN, double invSqrteRecoN )

inline

Definition at line 488 of file StructsDUNE.h.

                                                                                                     {
 
  (*erec) += (*par) * invSqrteRecoN * eRecoN;
  
}

NEScaleInvSqrtND()

void NEScaleInvSqrtND ( const double * par, double * erec, double eRecoN, double invSqrteRecoN )

inline

Definition at line 238 of file StructsDUNE.h.

                                                                                                     {
 
  (*erec) += (*par) * invSqrteRecoN * eRecoN;
  
}

NEScaleND()

void NEScaleND ( const double * par, double * erec, double eRecoN )

inline

Definition at line 224 of file StructsDUNE.h.

                                                                        {
 
  (*erec) += (*par) * eRecoN;
  
}

NEScaleSqrtFD()

void NEScaleSqrtFD ( const double * par, double * erec, double eRecoN, double sqrteRecoN )

inline

Definition at line 481 of file StructsDUNE.h.

                                                                                               {
 
  (*erec) += (*par) * sqrteRecoN * eRecoN;
  
}

NEScaleSqrtND()

void NEScaleSqrtND ( const double * par, double * erec, double eRecoN, double sqrteRecoN )

inline

Definition at line 231 of file StructsDUNE.h.

                                                                                               {
 
  (*erec) += (*par) * sqrteRecoN * eRecoN;
  
}

NResFD()

void NResFD ( const double * par, double * erec, double eRecoN, double eN )

inline

Definition at line 574 of file StructsDUNE.h.

                                                                                {
 
  (*erec) += (*par) * (eN - eRecoN);
  
}

NResND()

void NResND ( const double * par, double * erec, double eRecoN, double eN )

inline

Definition at line 324 of file StructsDUNE.h.

                                                                                {
 
  (*erec) += (*par) * (eN - eRecoN);
  
}

SIMBMode_ToMaCh3Mode()

int SIMBMode_ToMaCh3Mode ( int SIMB_mode, int isCC )

inline

Definition at line 763 of file StructsDUNE.h.

                                                         {
  // **********************
 
  int ReturnMode = kMaCh3_nModes;
  
  if (isCC == 1) {
    switch (SIMB_mode) {
      //Unknown
    case kUnknownInteraction:
      std::cerr << "Unknown Interaction mode:" << SIMB_mode << std::endl;
      throw;
      //QE
    case kQE:
      ReturnMode = kMaCh3_CCQE; // CC QE in MaCh3 DUNE
      break;
      // DIS
    case kDIS:
      ReturnMode = kMaCh3_CC_DIS; //CC DIS in MaCh3 DUNE
      break;
      // RES
    case kRes:
      ReturnMode = kMaCh3_CC_RES; // CC RES in MaCh3 DUNE
      break;
      // Coherent
    case kCoh:
      ReturnMode = kMaCh3_CC_COH; // CC Coh in MaCh3 DUNE
      break;
      // Diffractive
    case kDiffractive:
      ReturnMode = kMaCh3_CC_Diffractive; // CC multi-pi in MaCh3
      break;
      // Nue Elastic
    case kElectronScattering:
      ReturnMode = kMaCh3_CC_Nue_EL; // CC Nue scattering in MaCh3 DUNE
      break;
      // Atmospheric Mu Gamma
    case kAMNuGamma:
      ReturnMode = kMaCh3_CC_AMnuGamma; // CC Am Nu Mu in MaCh3 DUNE
      break;
      // MEC
    case kMEC:
      ReturnMode = kMaCh3_CC_MEC; // CC MEC in MaCh3 DUNE
      break;
      // Coherent Elastic
    case kCohElastic:
      ReturnMode = kMaCh3_CC_COHEL; // CC Coherent Elastic in MaCh3 DUNE
      break;
      // Inverse Beta Decay
    case kInverseBetaDecay:
      ReturnMode = kMaCh3_CC_IBD; // CC Inverse Beta Decay in MaCh3 DUNE
      break;
      // Glashow Resonance
    case kGlashowResonance:
      ReturnMode = kMaCh3_CC_GlashowRES; // CC Glashow Reaonance in DUNE
      break;
      // IMD Annihilation
    case kIMDAnnihilation:
      ReturnMode = kMaCh3_CC_IMDAnnihalation; // CC Inverse Muon Decay Annihalation in DUNE
      break;
    default:
      std::cerr << "Invalid mode:" << SIMB_mode << std::endl;
      throw;
    }
    
  }
  
  
  if (isCC == 0) {
    switch (SIMB_mode) {
      //Unknown
    case kUnknownInteraction:
      std::cerr << "Unknown Interaction mode:" << SIMB_mode << std::endl;
      throw;
      //QE
    case kQE:
      ReturnMode = kMaCh3_NCQE; // NC QE in MaCh3 DUNE
      break;
      // DIS
    case kDIS:
      ReturnMode = kMaCh3_NC_DIS; // NC DIS in MaCh3 DUNE
      break;
      // RES
    case kRes:
      ReturnMode = kMaCh3_NC_RES; // NC RES in MaCh3 DUNE
      break;
      // Coherent
    case kCoh:
      ReturnMode = kMaCh3_NC_COH; // NC Coh in MaCh3 DUNE
      break;
      // Diffractive
    case kDiffractive:
      ReturnMode = kMaCh3_NC_Diffractive; // CC multi-pi in MaCh3
      break;
      // Nue Elastic
    case kElectronScattering:
      ReturnMode = kMaCh3_NC_Nue_EL; // NC Nue scattering in MaCh3 DUNE
      break;
      // Atmospheric Mu Gamma
    case kAMNuGamma:
      ReturnMode = kMaCh3_NC_AMnuGamma; // NC Am Nu Mu in MaCh3 DUNE
      break;
      // MEC
    case kMEC:
      ReturnMode = kMaCh3_NC_MEC; // NC MEC in MaCh3 DUNE
      break;
      // Coherent Elastic
    case kCohElastic:
      ReturnMode = kMaCh3_NC_COHEL; // NC Coherent Elastic in MaCh3 DUNE
      break;
      // Inverse Beta Decay
    case kInverseBetaDecay:
      ReturnMode = kMaCh3_NC_IBD; // Inverse Beta Decay in MaCh3 DUNE
      break;
      // Glashow Resonance
    case kGlashowResonance:
      ReturnMode = kMaCh3_NC_GlashowRES; 
      break;
      // IMD Annihilation
    case kIMDAnnihilation:
      ReturnMode = kMaCh3_NC_IMDAnnihalation; 
      break;
    default:
      std::cerr << "Invalid mode:" << SIMB_mode << std::endl;
      throw;
    }
    
  }
 
  return ReturnMode;
};

TotalEScaleFD()

void TotalEScaleFD ( const double * par, double * erec, double erecHad, double erecLep, bool NotCCnumu )

inline

Definition at line 357 of file StructsDUNE.h.

                                                                                                             {
 
  (*erec) += (*par) * erecHad;
 
  //if not true CC numu event AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * erecLep;
  }
 
}

TotalEScaleInvSqrtFD()

void TotalEScaleInvSqrtFD ( const double * par, double * erec, double erecHad, double erecLep, double invSqrtErecHad, double invSqrtErecLep, bool NotCCnumu )

inline

Definition at line 383 of file StructsDUNE.h.

                                                                                                                                                                  {
 
  (*erec) += (*par) * invSqrtErecHad * erecHad ;
 
  //if not true CC numu AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
 
}

TotalEScaleInvSqrtND()

void TotalEScaleInvSqrtND ( const double * par, double * erec, double erecHad, double erecLep, double invSqrtErecHad, double invSqrtErecLep, bool NotCCnumu )

inline

Definition at line 133 of file StructsDUNE.h.

                                                                                                                                                                  {
 
  (*erec) += (*par) * invSqrtErecHad * erecHad ;
 
  //if not true CC numu AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * invSqrtErecLep * erecLep;
  }
 
}

TotalEScaleND()

void TotalEScaleND ( const double * par, double * erec, double erecHad, double erecLep, bool NotCCnumu )

inline

Definition at line 107 of file StructsDUNE.h.

                                                                                                             {
 
  (*erec) += (*par) * erecHad;
 
  //if not true CC numu event AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * erecLep;
  }
 
}

TotalEScaleSqrtFD()

void TotalEScaleSqrtFD ( const double * par, double * erec, double erecHad, double erecLep, double sqrtErecHad, double sqrtErecLep, bool NotCCnumu )

inline

Definition at line 370 of file StructsDUNE.h.

                                                                                                                                                         {
 
  (*erec) += (*par) * sqrtErecHad * erecHad ;
 
  //if not true CC numu AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
 
}

TotalEScaleSqrtND()

void TotalEScaleSqrtND ( const double * par, double * erec, double erecHad, double erecLep, double sqrtErecHad, double sqrtErecLep, bool NotCCnumu )

inline

Definition at line 120 of file StructsDUNE.h.

                                                                                                                                                         {
 
  (*erec) += (*par) * sqrtErecHad * erecHad ;
 
  //if not true CC numu AND reco nue event
  if (NotCCnumu)
  {
    (*erec) += (*par) * sqrtErecLep * erecLep;
  }
 
}