MCNeutrino.h

Go to the documentation of this file.

#ifndef SIMB_MCNEUTRINO_H
#define SIMB_MCNEUTRINO_H

#include "nusimdata/SimulationBase/MCParticle.h"

namespace simb {

  //......................................................................

  class MCNeutrino {
  public:

    MCNeutrino();

  private:

    simb::MCParticle fNu;              
    simb::MCParticle fLepton;          
    int              fMode;            
    int              fInteractionType; 
    int              fCCNC;            
    int              fTarget;          
    int              fHitNuc;          
    int              fHitQuark;        
    double           fW;               
    double           fX;               
    double           fY;               
    double           fQSqr;            

  public:

    MCNeutrino(simb::MCParticle &nu,
               simb::MCParticle &lep,
               int CCNC,
               int mode,
               int interactionType,
               int target,
               int nucleon,
               int quark,
               double w,
               double x,
               double y,
               double qsqr);

    const  simb::MCParticle& Nu()              const;
    const  simb::MCParticle& Lepton()          const;
    int                      CCNC()            const;
    int                      Mode()            const;
    int                      InteractionType() const;
    int                      Target()          const;
    int                      HitNuc()          const;
    int                      HitQuark()        const;
    double                   W()               const;
    double                   X()               const;
    double                   Y()               const;
    double                   QSqr()            const;
    double                   Pt()              const; 
    double                   Theta()           const; 
    friend std::ostream&  operator<< (std::ostream& output, const simb::MCNeutrino &mcnu);
  };
}

namespace simb{
  enum curr_type_{
    kCC,
    kNC
  };

  enum int_type_{
    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,
    kNuanceOffset              = 1000,                
    kCCQE                      = kNuanceOffset +  1,  
    kNCQE                      = kNuanceOffset +  2,  
    kResCCNuProtonPiPlus       = kNuanceOffset +  3,  
    kResCCNuNeutronPi0         = kNuanceOffset +  4,  
    kResCCNuNeutronPiPlus      = kNuanceOffset +  5,  
    kResNCNuProtonPi0          = kNuanceOffset +  6,  
    kResNCNuProtonPiPlus       = kNuanceOffset +  7,  
    kResNCNuNeutronPi0         = kNuanceOffset +  8,  
    kResNCNuNeutronPiMinus     = kNuanceOffset +  9,  
    kResCCNuBarNeutronPiMinus  = kNuanceOffset + 10,  
    kResCCNuBarProtonPi0       = kNuanceOffset + 11,  
    kResCCNuBarProtonPiMinus   = kNuanceOffset + 12,  
    kResNCNuBarProtonPi0       = kNuanceOffset + 13,  
    kResNCNuBarProtonPiPlus    = kNuanceOffset + 14,  
    kResNCNuBarNeutronPi0      = kNuanceOffset + 15,  
    kResNCNuBarNeutronPiMinus  = kNuanceOffset + 16,  
    kResCCNuDeltaPlusPiPlus    = kNuanceOffset + 17,
    kResCCNuDelta2PlusPiMinus  = kNuanceOffset + 21,
    kResCCNuBarDelta0PiMinus   = kNuanceOffset + 28,
    kResCCNuBarDeltaMinusPiPlus= kNuanceOffset + 32,
    kResCCNuProtonRhoPlus      = kNuanceOffset + 39,
    kResCCNuNeutronRhoPlus     = kNuanceOffset + 41,
    kResCCNuBarNeutronRhoMinus = kNuanceOffset + 46,
    kResCCNuBarNeutronRho0     = kNuanceOffset + 48,
    kResCCNuSigmaPlusKaonPlus  = kNuanceOffset + 53,
    kResCCNuSigmaPlusKaon0     = kNuanceOffset + 55,
    kResCCNuBarSigmaMinusKaon0 = kNuanceOffset + 60,
    kResCCNuBarSigma0Kaon0     = kNuanceOffset + 62,
    kResCCNuProtonEta          = kNuanceOffset + 67,
    kResCCNuBarNeutronEta      = kNuanceOffset + 70,
    kResCCNuKaonPlusLambda0    = kNuanceOffset + 73,
    kResCCNuBarKaon0Lambda0    = kNuanceOffset + 76,
    kResCCNuProtonPiPlusPiMinus= kNuanceOffset + 79,
    kResCCNuProtonPi0Pi0       = kNuanceOffset + 80,
    kResCCNuBarNeutronPiPlusPiMinus = kNuanceOffset + 85,
    kResCCNuBarNeutronPi0Pi0   = kNuanceOffset + 86,
    kResCCNuBarProtonPi0Pi0    = kNuanceOffset + 90,
    kCCDIS                     = kNuanceOffset + 91,  
    kNCDIS                     = kNuanceOffset + 92,  
    kUnUsed1                   = kNuanceOffset + 93,
    kUnUsed2                   = kNuanceOffset + 94,
    kCCQEHyperon               = kNuanceOffset + 95,
    kNCCOH                     = kNuanceOffset + 96,
    kCCCOH                     = kNuanceOffset + 97,  
    kNuElectronElastic         = kNuanceOffset + 98,  
    kInverseMuDecay            = kNuanceOffset + 99,  
    kMEC2p2h                   = kNuanceOffset + 100  
  };
}

inline const  simb::MCParticle& simb::MCNeutrino::Nu()              const { return fNu;              }
inline const  simb::MCParticle& simb::MCNeutrino::Lepton()          const { return fLepton;          }
inline        int               simb::MCNeutrino::CCNC()            const { return fCCNC;            }
inline        int               simb::MCNeutrino::Mode()            const { return fMode;            }
inline        int               simb::MCNeutrino::InteractionType() const { return fInteractionType; }
inline        int               simb::MCNeutrino::Target()          const { return fTarget;          }
inline        int               simb::MCNeutrino::HitNuc()          const { return fHitNuc;          }
inline        int               simb::MCNeutrino::HitQuark()        const { return fHitQuark;        }
inline        double            simb::MCNeutrino::W()               const { return fW;               }
inline        double            simb::MCNeutrino::X()               const { return fX;               }
inline        double            simb::MCNeutrino::Y()               const { return fY;               }
inline        double            simb::MCNeutrino::QSqr()            const { return fQSqr;            }

#endif //SIMB_MCNEUTRINO_H