doxygen/Interaction_8cxx_source.html

 //____________________________________________________________________________
 /*
  Copyright (c) 2003-2020, The GENIE Collaboration
  For the full text of the license visit http://copyright.genie-mc.org

  Costas Andreopoulos <constantinos.andreopoulos \at cern.ch>
  University of Liverpool & STFC Rutherford Appleton Laboratory

  Changes required to implement the GENIE Boosted Dark Matter module
  were installed by Josh Berger (Univ. of Wisconsin)

  Changes required to implement the GENIE Dark Neutrino module
  were installed by Iker de Icaza (Univ. of Sussex)
 */
 //____________________________________________________________________________

 #include <sstream>

 #include <TRootIOCtor.h>

 #include "Framework/Conventions/Constants.h"
 #include "Framework/Interaction/Interaction.h"
 #include "Framework/Messenger/Messenger.h"
 #include "Framework/ParticleData/PDGLibrary.h"
 #include "Framework/ParticleData/PDGCodes.h"
 #include "Framework/ParticleData/PDGUtils.h"

 using namespace genie;
 using namespace genie::constants;

 using std::endl;
 using std::ostringstream;

 ClassImp(Interaction)

 //____________________________________________________________________________
 namespace genie {
  ostream & operator<< (ostream& stream, const Interaction & interaction)
  {
    interaction.Print(stream);
    return stream;
  }
 }
 //___________________________________________________________________________
 Interaction::Interaction() :
 TObject()
 {
   this->Init();
 }
 //___________________________________________________________________________
 Interaction::Interaction(const InitialState & ist, const ProcessInfo & prc) :
 TObject()
 {
   this->Init();

   fInitialState -> Copy (ist);
   fProcInfo     -> Copy (prc);
 }
 //___________________________________________________________________________
 Interaction::Interaction(const Interaction & interaction) :
 TObject()
 {
   this->Init();
   this->Copy(interaction);
 }
 //___________________________________________________________________________
 Interaction::Interaction(TRootIOCtor*) :
 TObject(),
 fInitialState(0),
 fProcInfo(0),
 fKinematics(0),
 fExclusiveTag(0),
 fKinePhSp(0)
 {

 }
 //___________________________________________________________________________
 Interaction::~Interaction()
 {
   this->CleanUp();
 }
 //___________________________________________________________________________
 void Interaction::Reset(void)
 {
   this->CleanUp();
   this->Init();
 }
 //___________________________________________________________________________
 void Interaction::Init(void)
 {
   fInitialState = new InitialState ();
   fProcInfo     = new ProcessInfo  ();
   fKinematics   = new Kinematics   ();
   fExclusiveTag = new XclsTag      ();
   fKinePhSp     = new KPhaseSpace  (this);
 }
 //___________________________________________________________________________
 void Interaction::CleanUp(void)
 {
   if ( fInitialState ) delete fInitialState;
   if ( fProcInfo     ) delete fProcInfo;
   if ( fKinematics   ) delete fKinematics;
   if ( fExclusiveTag ) delete fExclusiveTag;
   if ( fKinePhSp     ) delete fKinePhSp;

   fInitialState = 0;
   fProcInfo     = 0;
   fKinematics   = 0;
   fExclusiveTag = 0;
   fKinePhSp     = 0;
 }
 //___________________________________________________________________________
 void Interaction::Copy(const Interaction & interaction)
 {
   const InitialState & init = *interaction.fInitialState;
   const ProcessInfo &  proc = *interaction.fProcInfo;
   const Kinematics &   kine = *interaction.fKinematics;
   const XclsTag &      xcls = *interaction.fExclusiveTag;

   fInitialState -> Copy (init);
   fProcInfo     -> Copy (proc);
   fKinematics   -> Copy (kine);
   fExclusiveTag -> Copy (xcls);
 }
 //___________________________________________________________________________
 TParticlePDG * Interaction::FSPrimLepton(void) const
 {
   int pdgc = this->FSPrimLeptonPdg();

   if(pdgc) return PDGLibrary::Instance()->Find(pdgc);
   else     return 0;
 }
 //___________________________________________________________________________
 int Interaction::FSPrimLeptonPdg(void) const
 {
   const ProcessInfo &  proc_info  = this -> ProcInfo();
   const InitialState & init_state = this -> InitState();
   const XclsTag &      xclstag    = this -> ExclTag();

   int pdgc = init_state.ProbePdg();

   LOG("Interaction", pDEBUG) << "Probe PDG code: " << pdgc;

   if (proc_info.IsNuElectronElastic())
     return kPdgElectron;

   // vN (Weak-NC) or eN (EM)
   if (proc_info.IsWeakNC() || proc_info.IsEM() || proc_info.IsWeakMix() || proc_info.IsDarkMatter()) return pdgc;  // EDIT: DM does not change in FS

   // vN (Weak-CC)
   else if (proc_info.IsWeakCC()) {
     int clpdgc;
     if (proc_info.IsIMDAnnihilation())
       clpdgc = kPdgMuon;
     else if (proc_info.IsGlashowResonance()) {
       if      ( pdg::IsMuon(xclstag.FinalLeptonPdg())     ) clpdgc = kPdgMuon;
       else if ( pdg::IsTau(xclstag.FinalLeptonPdg())      ) clpdgc = kPdgTau;
       else if ( pdg::IsElectron(xclstag.FinalLeptonPdg()) ) clpdgc = kPdgElectron;
       else if ( pdg::IsPion(xclstag.FinalLeptonPdg())     ) clpdgc = kPdgPiP;
     }
     else
       clpdgc = pdg::Neutrino2ChargedLepton(pdgc);
     return clpdgc;
   }
   else if (proc_info.IsDarkNeutralCurrent()){
     return kPdgDarkNeutrino;
   }
   LOG("Interaction", pWARN)
         << "Could not figure out the final state primary lepton pdg code!!";

   return 0;
 }
 //___________________________________________________________________________
 TParticlePDG * Interaction::RecoilNucleon(void) const
 {
   int rnuc = this->RecoilNucleonPdg();

   if(rnuc) return PDGLibrary::Instance()->Find(rnuc);
   else     return 0;
 }
 //___________________________________________________________________________
 int Interaction::RecoilNucleonPdg(void) const
 {
 // Determine the recoil nucleon PDG code

   const Target & target = fInitialState->Tgt();

   int recoil_nuc = 0;
   int struck_nuc = target.HitNucPdg();

   if (fProcInfo->IsQuasiElastic() || fProcInfo->IsInverseBetaDecay() || fProcInfo->IsDarkMatterElastic()) {
     bool struck_is_nuc = pdg::IsNucleon(struck_nuc);
     bool is_weak = fProcInfo->IsWeak();
     bool is_em   = fProcInfo->IsEM();
     bool is_dm   = fProcInfo->IsDarkMatter();
     assert(struck_is_nuc && (is_weak || is_em || is_dm));
     if(fProcInfo->IsWeakCC()) {
        recoil_nuc = pdg::SwitchProtonNeutron(struck_nuc); // CC
     } else {
        recoil_nuc = struck_nuc; // NC, EM
     }
   }

   if (fProcInfo->IsMEC()) {
     bool struck_is_2nuc_cluster = pdg::Is2NucleonCluster(struck_nuc);
     bool is_weak = fProcInfo->IsWeak();
     bool is_em   = fProcInfo->IsEM();
     assert(struck_is_2nuc_cluster && (is_weak || is_em));
     if(fProcInfo->IsWeakCC()) {
        bool isnu = pdg::IsNeutrino(fInitialState->ProbePdg());
        // nucleon cluster charge should be incremented by +1 for
        // neutrino CC and by -1 for antineutrino CC
        int dQ = (isnu) ? +1 : -1;
        recoil_nuc = pdg::ModifyNucleonCluster(struck_nuc,dQ); // CC
     }
     else {
        recoil_nuc = struck_nuc; // NC, EM
     }
   }

   LOG("Interaction", pDEBUG) << "Recoil nucleon PDG = " << recoil_nuc;
   return recoil_nuc;
 }
 //___________________________________________________________________________
 void Interaction::SetInitState(const InitialState & init_state)
 {
   if (!fInitialState) fInitialState = new InitialState();
   fInitialState->Copy(init_state);
 }
 //___________________________________________________________________________
 void Interaction::SetProcInfo(const ProcessInfo & proc_info)
 {
   if (!fProcInfo) fProcInfo = new ProcessInfo();
   fProcInfo->Copy(proc_info);
 }
 //___________________________________________________________________________
 void Interaction::SetKine(const Kinematics & kinematics)
 {
   if (!fKinematics) fKinematics = new Kinematics();
   fKinematics->Copy(kinematics);
 }
 //___________________________________________________________________________
 void Interaction::SetExclTag(const XclsTag & xcls_tag)
 {
   if (!fExclusiveTag) fExclusiveTag = new XclsTag();
   fExclusiveTag->Copy(xcls_tag);
 }
 //___________________________________________________________________________
 string Interaction::AsString(void) const
 {
 // Code-ify the interaction in a string to be used as (part of a) cache
 // branch key.
 // Template:
 // nu:x;tgt:x;N:x;q:x(s/v);proc:x;xclv_tag

   const Target & tgt = fInitialState->Tgt();

   ostringstream interaction;

   // If the probe has non-zero mass, then it is DM
   if (fInitialState->Probe()->PdgCode() == kPdgDarkMatter) {
     interaction << "dm;";
   }
   else if (fInitialState->Probe()->PdgCode() == kPdgAntiDarkMatter) {
     interaction << "dmb;";
   }
   else {
     interaction << "nu:"  << fInitialState->ProbePdg() << ";";
   }
   interaction << "tgt:" << tgt.Pdg() << ";";

   if(tgt.HitNucIsSet()) {
     interaction << "N:" << tgt.HitNucPdg() << ";";
   }
   if(tgt.HitQrkIsSet()) {
     interaction << "q:" << tgt.HitQrkPdg()
                 << (tgt.HitSeaQrk() ? "(s)" : "(v)") << ";";
   }

   interaction << "proc:" << fProcInfo->InteractionTypeAsString()
               << "," << fProcInfo->ScatteringTypeAsString()  << ";";

   string xcls = fExclusiveTag->AsString();
   interaction << xcls;
   if(xcls.size()>0) interaction << ";";

   return interaction.str();
 }
 //___________________________________________________________________________
 void Interaction::Print(ostream & stream) const
 {
   const string line(110, '-');

   stream << endl;
   stream << line << endl;

   stream << "GENIE Interaction Summary" << endl;
   stream << line << endl;

   stream << *fInitialState << endl; // print initial state
   stream << *fProcInfo;             // print process info
   stream << *fKinematics;           // print scattering parameters
   stream << *fExclusiveTag;         // print exclusive process tag

   stream << line << endl;
 }
 //___________________________________________________________________________
 Interaction & Interaction::operator = (const Interaction & interaction)
 {
   this->Copy(interaction);
   return (*this);
 }
 //___________________________________________________________________________
 //
 //       **** Methods using the "named constructor" C++ idiom ****
 //
 //___________________________________________________________________________
 Interaction * Interaction::Create(
             int target, int probe, ScatteringType_t st, InteractionType_t it)
 {
   InitialState init_state (target, probe);
   ProcessInfo  proc_info  (st, it);

   Interaction * interaction = new Interaction(init_state, proc_info);
   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISCC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
              Interaction::Create(target,probe,kScDeepInelastic, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISCC(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe, double E)
 {
   Interaction* interaction = Interaction::DISCC(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISCC(
    int target, int hitnuc, int hitqrk, bool fromsea, int fqrk, int probe, double E)
 {
   Interaction* interaction = Interaction::DISCC(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   XclsTag * xclstag = interaction->ExclTagPtr();
   xclstag->SetFinalQuark(fqrk);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISCC(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScDeepInelastic, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISCC(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe,
    const TLorentzVector & p4probe)
 {
   Interaction* interaction = Interaction::DISCC(target,hitnuc,probe,p4probe);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISNC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
              Interaction::Create(target,probe,kScDeepInelastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISNC(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe, double E)
 {
   Interaction* interaction = Interaction::DISNC(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISNC(
    int target, int hitnuc, int hitqrk, bool fromsea, int fqrk, int probe, double E)
 {
   Interaction* interaction = Interaction::DISNC(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   XclsTag * xclstag = interaction->ExclTagPtr();
   xclstag->SetFinalQuark(fqrk);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISNC(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScDeepInelastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISNC(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe,
    const TLorentzVector & p4probe)
 {
   Interaction * interaction = Interaction::DISNC(target,hitnuc,probe,p4probe);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISEM(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
              Interaction::Create(target,probe,kScDeepInelastic, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISEM(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe, double E)
 {
   Interaction* interaction = Interaction::DISEM(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISEM(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScDeepInelastic, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DISEM(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe,
    const TLorentzVector & p4probe)
 {
   Interaction * interaction = Interaction::DISEM(target,hitnuc,probe,p4probe);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELCC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELCC(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELNC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELNC(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELEM(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::QELEM(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScQuasiElastic, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::IBD(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScInverseBetaDecay,kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::IBD(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
       Interaction::Create(target,probe,kScInverseBetaDecay,kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESCC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESCC(
            int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESNC(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESNC(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESEM(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::RESEM(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScResonant, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DFRCC(int tgt,int hitnuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScDiffractive, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DFRCC(
    int tgt, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScDiffractive, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::COHCC(int tgt, int probe, unsigned int prod_pdg, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentProduction, kIntWeakCC);

   XclsTag * xcl = interaction -> ExclTagPtr() ;
   if ( pdg::IsPion( prod_pdg ) ) {
     if ( pdg::IsNeutrino( probe ) ) xcl -> SetNPions( 1,0,0 ) ;
     else if ( pdg::IsAntiNeutrino( probe ) ) xcl -> SetNPions( 0,0,1 ) ;
   }

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::COHCC(
     int tgt, int probe, unsigned int prod_pdg, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentProduction, kIntWeakCC);

   XclsTag * xcl = interaction -> ExclTagPtr() ;
   if ( pdg::IsPion( prod_pdg ) ) {
     if ( pdg::IsNeutrino( probe ) ) xcl -> SetNPions( 1,0,0 ) ;
     else if ( pdg::IsAntiNeutrino( probe ) ) xcl -> SetNPions( 0,0,1 ) ;
   }

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::COHNC(int tgt, int probe, unsigned int prod_pdg, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentProduction, kIntWeakNC);

   XclsTag * xcl = interaction -> ExclTagPtr() ;
   if ( pdg::IsPion( prod_pdg ) ) xcl -> SetNPions( 0,1,0 ) ;
   else if ( prod_pdg == kPdgGamma )  xcl -> SetNSingleGammas(1) ;

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::COHNC(
    int tgt, int probe, unsigned int prod_pdg, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentProduction, kIntWeakNC);

   XclsTag * xcl = interaction -> ExclTagPtr() ;
   if ( pdg::IsPion( prod_pdg ) ) xcl -> SetNPions( 0,1,0 ) ;
   else if ( prod_pdg == kPdgGamma )  xcl -> SetNSingleGammas(1) ;

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::CEvNS(int tgt, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentElastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::CEvNS(
    int tgt, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScCoherentElastic, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::IMD(int target, double E)
 {
   Interaction * interaction =
      Interaction::Create(target,kPdgNuMu,kScInverseMuDecay, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::IMD(int target, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,kPdgNuMu,kScInverseMuDecay, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::AMNuGamma(int tgt, int nuc, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScAMNuGamma, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(nuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::AMNuGamma(
    int tgt, int nuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScAMNuGamma, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(nuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECCC(int tgt, int ncluster, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECCC(
    int tgt, int ncluster, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECCC(int tgt, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECCC(
    int tgt, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }

 //___________________________________________________________________________
 Interaction * Interaction::MECNC(int tgt, int ncluster, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECNC(
    int tgt, int ncluster, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntWeakNC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECEM(int tgt, int probe, double E)
 {

   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECEM(int tgt, int ncluster, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::MECEM(
    int tgt, int ncluster, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, probe, kScMEC, kIntEM);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(ncluster);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::GLR(int tgt, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt, kPdgAntiNuE, kScGlashowResonance, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(0);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::GLR(int tgt, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt, kPdgAntiNuE, kScGlashowResonance, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(0);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::NDecay(int tgt, int decay_mode, int decayed_nucleon)
 {
   Interaction * interaction =
      Interaction::Create(tgt, 0, kScNull, kIntNDecay);
   interaction->ExclTagPtr()->SetDecayMode(decay_mode);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->TgtPtr()->SetHitNucPdg(decayed_nucleon);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::NOsc(int tgt, int annihilation_mode)
 {
   Interaction * interaction =
     Interaction::Create(tgt, 0, kScNull, kIntNOsc);
   interaction->ExclTagPtr()->SetDecayMode(annihilation_mode);
   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::NHL(double E, int decayed_mode)
 {
   Interaction * interaction =
     Interaction::Create(0, 0, kScNull, kIntNHL);
   interaction->ExclTagPtr()->SetDecayMode(decayed_mode);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::ASK(int tgt, int probe, double E)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScSingleKaon, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::ASK(
     int tgt, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(tgt,probe,kScSingleKaon, kIntWeakCC);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DME(int target, int hitnuc, int probe, double E)
 {
   // EDIT: need to be able to create dark matter elastic
   Interaction * interaction =
      Interaction::Create(target,probe,kScDarkMatterElastic, kIntDarkMatter);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DME(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   // EDIT: need to be able to create dark matter elastic
   Interaction * interaction =
      Interaction::Create(target,probe,kScDarkMatterElastic, kIntDarkMatter);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DMDI(int target, int hitnuc, int probe, double E)
 {
   Interaction * interaction =
              Interaction::Create(target,probe,kScDarkMatterDeepInelastic, kIntDarkMatter);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeE(E);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DMDI(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe, double E)
 {
   Interaction* interaction = Interaction::DMDI(target,hitnuc,probe,E);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DMDI(
    int target, int hitnuc, int probe, const TLorentzVector & p4probe)
 {
   Interaction * interaction =
      Interaction::Create(target,probe,kScDarkMatterDeepInelastic, kIntDarkMatter);

   InitialState * init_state = interaction->InitStatePtr();
   init_state->SetProbeP4(p4probe);
   init_state->TgtPtr()->SetHitNucPdg(hitnuc);

   return interaction;
 }
 //___________________________________________________________________________
 Interaction * Interaction::DMDI(
    int target, int hitnuc, int hitqrk, bool fromsea, int probe,
    const TLorentzVector & p4probe)
 {
   Interaction * interaction = Interaction::DMDI(target,hitnuc,probe,p4probe);

   Target * tgt = interaction->InitStatePtr()->TgtPtr();
   tgt -> SetHitQrkPdg (hitqrk);
   tgt -> SetHitSeaQrk (fromsea);

   return interaction;
 }
genie::pdg::IsPion
bool IsPion(int pdgc)
Definition: PDGUtils.cxx:323

genie::constants
Basic constants.

genie::ProcessInfo::IsWeak
bool IsWeak(void) const
Definition: ProcessInfo.cxx:178

genie::Interaction::IMD
static Interaction * IMD(int tgt, double E=0)
Definition: Interaction.cxx:797

genie::ProcessInfo::IsWeakMix
bool IsWeakMix(void) const
Definition: ProcessInfo.cxx:198

genie::Target::HitSeaQrk
bool HitSeaQrk(void) const
Definition: Target.cxx:299

genie::ProcessInfo::IsWeakCC
bool IsWeakCC(void) const
Definition: ProcessInfo.cxx:183

genie::kIntWeakNC
Definition: InteractionType.h:38

genie::XclsTag::FinalLeptonPdg
int FinalLeptonPdg(void) const
Definition: XclsTag.h:74

genie::pdg::IsNeutrino
bool IsNeutrino(int pdgc)
Definition: PDGUtils.cxx:107

genie::InitialState::SetProbeP4
void SetProbeP4(const TLorentzVector &P4)
Definition: InitialState.cxx:183

genie
THE MAIN GENIE PROJECT NAMESPACE
Definition: AlgCmp.h:25

genie::Interaction::DME
static Interaction * DME(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:1034

genie::ProcessInfo::ScatteringTypeAsString
string ScatteringTypeAsString(void) const
Definition: ProcessInfo.cxx:226

genie::Interaction::GLR
static Interaction * GLR(int tgt, double E=0)
Definition: Interaction.cxx:955

genie::Target::HitNucPdg
int HitNucPdg(void) const
Definition: Target.cxx:304

genie::Interaction::RecoilNucleon
TParticlePDG * RecoilNucleon(void) const
recoil nucleon
Definition: Interaction.cxx:174

genie::ProcessInfo::IsDarkNeutralCurrent
bool IsDarkNeutralCurrent(void) const
Definition: ProcessInfo.cxx:136

genie::pdg::IsNucleon
bool IsNucleon(int pdgc)
Definition: PDGUtils.cxx:343

genie::Interaction::RecoilNucleonPdg
int RecoilNucleonPdg(void) const
recoil nucleon pdg
Definition: Interaction.cxx:182

genie::Target::HitQrkPdg
int HitQrkPdg(void) const
Definition: Target.cxx:242

genie::Interaction::fInitialState
InitialState * fInitialState
Initial State info.
Definition: Interaction.h:179

genie::ProcessInfo::IsQuasiElastic
bool IsQuasiElastic(void) const
Definition: ProcessInfo.cxx:69

genie::kPdgAntiNuE
const int kPdgAntiNuE
Definition: PDGCodes.h:29

genie::kPdgDarkMatter
const int kPdgDarkMatter
Definition: PDGCodes.h:218

genie::Interaction::Reset
void Reset(void)
Definition: Interaction.cxx:83

genie::Target::Pdg
int Pdg(void) const
Definition: Target.h:71

genie::kPdgNuMu
const int kPdgNuMu
Definition: PDGCodes.h:30

genie::Interaction::CEvNS
static Interaction * CEvNS(int tgt, int probe, double E=0)
Definition: Interaction.cxx:774

genie::kScInverseBetaDecay
Definition: ScatteringType.h:50

genie::kScDarkMatterElastic
Definition: ScatteringType.h:53

genie::Kinematics
Generated/set kinematical variables for an event.
Definition: Kinematics.h:39

genie::pdg::SwitchProtonNeutron
int SwitchProtonNeutron(int pdgc)
Definition: PDGUtils.cxx:353

genie::kScGlashowResonance
Definition: ScatteringType.h:51

genie::XclsTag::Copy
void Copy(const XclsTag &xcls)
copy input XclsTag object
Definition: XclsTag.cxx:173

genie::ProcessInfo::IsInverseBetaDecay
bool IsInverseBetaDecay(void) const
Definition: ProcessInfo.cxx:141

genie::InitialState::Probe
TParticlePDG * Probe(void) const
Definition: InitialState.cxx:147

genie::Interaction::IBD
static Interaction * IBD(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:583

genie::Interaction::NDecay
static Interaction * NDecay(int tgt, int decay_mode=-1, int decayed_nucleon=0)
Definition: Interaction.cxx:979

genie::InitialState::SetProbeE
void SetProbeE(double E)
Definition: InitialState.cxx:175

genie::ProcessInfo::IsIMDAnnihilation
bool IsIMDAnnihilation(void) const
Definition: ProcessInfo.cxx:126

Constants.h

genie::Interaction::SetKine
void SetKine(const Kinematics &kine)
Definition: Interaction.cxx:237

genie::kPdgElectron
const int kPdgElectron
Definition: PDGCodes.h:35

train.init
init
Definition: train.py:42

genie::Interaction::fExclusiveTag
XclsTag * fExclusiveTag
Additional info for exclusive channels.
Definition: Interaction.h:182

genie::Interaction::QELCC
static Interaction * QELCC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:508

genie::Interaction::MECNC
static Interaction * MECNC(int tgt, int nuccluster, int probe, double E=0)
Definition: Interaction.cxx:893

generate_CCQE_events.target
target
Definition: generate_CCQE_events.py:44

genie::Interaction::AsString
string AsString(void) const
Definition: Interaction.cxx:249

genie::Interaction::RESNC
static Interaction * RESNC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:633

genie::Interaction::COHCC
static Interaction * COHCC(int tgt, int probe, unsigned int prod_pdg, double E=0)
Definition: Interaction.cxx:708

genie::XclsTag
Contains minimal information for tagging exclusive processes.
Definition: XclsTag.h:39

genie::Interaction::FSPrimLeptonPdg
int FSPrimLeptonPdg(void) const
final state primary lepton pdg
Definition: Interaction.cxx:134

genie::kIntNHL
Definition: InteractionType.h:43

genie::kScCoherentElastic
Definition: ScatteringType.h:49

genie::pdg::ModifyNucleonCluster
int ModifyNucleonCluster(int pdgc, int dQ)
Definition: PDGUtils.cxx:361

genie::kScMEC
Definition: ScatteringType.h:48

genie::Interaction
Summary information for an interaction.
Definition: Interaction.h:56

genie::Interaction::DMDI
static Interaction * DMDI(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:1061

genie::kIntNOsc
Definition: InteractionType.h:42

Interaction.h

genie::Interaction::RESEM
static Interaction * RESEM(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:658

genie::ProcessInfo::IsWeakNC
bool IsWeakNC(void) const
Definition: ProcessInfo.cxx:188

LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96

genie::kPdgTau
const int kPdgTau
Definition: PDGCodes.h:39

genie::kIntDarkMatter
Definition: InteractionType.h:40

genie::ProcessInfo::IsNuElectronElastic
bool IsNuElectronElastic(void) const
Definition: ProcessInfo.cxx:116

genie::InitialState::Copy
void Copy(const InitialState &init_state)
Definition: InitialState.cxx:131

cvn::interaction
Definition: InteractionType.h:60

genie::kPdgAntiDarkMatter
const int kPdgAntiDarkMatter
Definition: PDGCodes.h:219

genie::XclsTag::AsString
string AsString(void) const
pack into a string code
Definition: XclsTag.cxx:212

genie::Interaction::MECEM
static Interaction * MECEM(int tgt, int nuccluster, int probe, double E=0)
Definition: Interaction.cxx:930

genie::Interaction::QELNC
static Interaction * QELNC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:533

genie::kScDeepInelastic
Definition: ScatteringType.h:41

genie::ProcessInfo
A class encapsulating an enumeration of interaction types (EM, Weak-CC, Weak-NC) and scattering types...
Definition: ProcessInfo.h:46

genie::Kinematics::Copy
void Copy(const Kinematics &kine)
Definition: Kinematics.cxx:83

genie::kPdgGamma
const int kPdgGamma
Definition: PDGCodes.h:189

genie::pdg::IsAntiNeutrino
bool IsAntiNeutrino(int pdgc)
Definition: PDGUtils.cxx:115

genie::Target
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
Definition: Target.h:40

genie::InitialState::ProbePdg
int ProbePdg(void) const
Definition: InitialState.h:64

genie::KPhaseSpace
Kinematical phase space.
Definition: KPhaseSpace.h:33

genie::Interaction::CleanUp
void CleanUp(void)
Definition: Interaction.cxx:98

genie::pdg::IsTau
bool IsTau(int pdgc)
Definition: PDGUtils.cxx:205

genie::kScCoherentProduction
Definition: ScatteringType.h:43

genie::ProcessInfo::InteractionTypeAsString
string InteractionTypeAsString(void) const
Definition: ProcessInfo.cxx:232

genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:158

generate_datataset.stream
stream
Definition: generate_datataset.py:30

genie::ProcessInfo::IsDarkMatterElastic
bool IsDarkMatterElastic(void) const
Definition: ProcessInfo.cxx:74

genie::Interaction::Copy
void Copy(const Interaction &i)
Definition: Interaction.cxx:113

genie::pdg::IsMuon
bool IsMuon(int pdgc)
Definition: PDGUtils.cxx:195

genie::kIntNDecay
Definition: InteractionType.h:41

Messenger.h

genie::XclsTag::SetDecayMode
void SetDecayMode(int decay_mode)
Definition: XclsTag.cxx:133

pWARN
#define pWARN
Definition: Messenger.h:60

genie::Interaction::FSPrimLepton
TParticlePDG * FSPrimLepton(void) const
final state primary lepton
Definition: Interaction.cxx:126

genie::ProcessInfo::IsMEC
bool IsMEC(void) const
Definition: ProcessInfo.cxx:156

genie::ProcessInfo::IsEM
bool IsEM(void) const
Definition: ProcessInfo.cxx:173

PDGUtils.h

genie::Interaction::ExclTagPtr
XclsTag * ExclTagPtr(void) const
Definition: Interaction.h:77

E
Definition: pointersEqual_t.cc:12

genie::pdg::Is2NucleonCluster
bool Is2NucleonCluster(int pdgc)
Definition: PDGUtils.cxx:399

genie::Interaction::SetExclTag
void SetExclTag(const XclsTag &xcls)
Definition: Interaction.cxx:243

genie::kIntEM
Definition: InteractionType.h:36

genie::kScDarkMatterDeepInelastic
Definition: ScatteringType.h:54

genie::Interaction::fKinematics
Kinematics * fKinematics
kinematical variables
Definition: Interaction.h:181

genie::PDGLibrary::Instance
static PDGLibrary * Instance(void)
Definition: PDGLibrary.cxx:57

genie::Interaction::ASK
static Interaction * ASK(int tgt, int probe, double E=0)
Definition: Interaction.cxx:1011

PDGLibrary.h

genie::Target::HitNucIsSet
bool HitNucIsSet(void) const
Definition: Target.cxx:283

genie::Target::HitQrkIsSet
bool HitQrkIsSet(void) const
Definition: Target.cxx:292

genie::Interaction::Print
void Print(ostream &stream) const
Definition: Interaction.cxx:290

genie::Interaction::QELEM
static Interaction * QELEM(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:558

genie::ScatteringType_t
enum genie::EScatteringType ScatteringType_t

genie::operator<<
ostream & operator<<(ostream &stream, const AlgConfigPool &config_pool)
Definition: AlgConfigPool.cxx:41

genie::ProcessInfo::IsDarkMatter
bool IsDarkMatter(void) const
Definition: ProcessInfo.cxx:193

genie::Target::SetHitNucPdg
void SetHitNucPdg(int pdgc)
Definition: Target.cxx:171

genie::Interaction::ExclTag
const XclsTag & ExclTag(void) const
Definition: Interaction.h:72

pduneana::line
void line(double t, double *p, double &x, double &y, double &z)
Definition: PDSPAnalyzer_module.cc:4741

genie::kScInverseMuDecay
Definition: ScatteringType.h:46

genie::Interaction::DISCC
static Interaction * DISCC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:328

genie::InitialState::TgtPtr
Target * TgtPtr(void) const
Definition: InitialState.h:67

genie::Interaction::Init
void Init(void)
Definition: Interaction.cxx:89

genie::Interaction::DFRCC
static Interaction * DFRCC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:683

genie::Interaction::~Interaction
~Interaction()
Definition: Interaction.cxx:78

genie::pdg::Neutrino2ChargedLepton
int Neutrino2ChargedLepton(int pdgc)
Definition: PDGUtils.cxx:215

genie::Interaction::InitStatePtr
InitialState * InitStatePtr(void) const
Definition: Interaction.h:74

genie::kScDiffractive
Definition: ScatteringType.h:44

genie::kIntWeakCC
Definition: InteractionType.h:37

genie::Interaction::InitState
const InitialState & InitState(void) const
Definition: Interaction.h:69

genie::Interaction::ProcInfo
const ProcessInfo & ProcInfo(void) const
Definition: Interaction.h:70

ClassImp
ClassImp(Interaction) namespace genie
Definition: Interaction.cxx:34

genie::Interaction::operator=
Interaction & operator=(const Interaction &i)
copy
Definition: Interaction.cxx:308

genie::kScAMNuGamma
Definition: ScatteringType.h:47

genie::PDGLibrary::Find
TParticlePDG * Find(int pdgc, bool must_exist=true)
Definition: PDGLibrary.cxx:75

genie::Interaction::fProcInfo
ProcessInfo * fProcInfo
Process info (scattering, weak current,...)
Definition: Interaction.h:180

genie::InitialState::Tgt
const Target & Tgt(void) const
Definition: InitialState.h:66

genie::kScNull
Definition: ScatteringType.h:38

genie::kPdgMuon
const int kPdgMuon
Definition: PDGCodes.h:37

genie::Interaction::NHL
static Interaction * NHL(double E=0, int decayed_mode=-1)
Definition: Interaction.cxx:999

genie::Interaction::RESCC
static Interaction * RESCC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:608

genie::kScResonant
Definition: ScatteringType.h:42

genie::ProcessInfo::IsGlashowResonance
bool IsGlashowResonance(void) const
Definition: ProcessInfo.cxx:146

genie::Interaction::MECCC
static Interaction * MECCC(int tgt, int nuccluster, int probe, double E=0)
Definition: Interaction.cxx:844

genie::Interaction::COHNC
static Interaction * COHNC(int tgt, int probe, unsigned int prod_pdg, double E=0)
Definition: Interaction.cxx:743

genie::XclsTag::SetFinalQuark
void SetFinalQuark(int finalquark_pdgc=0)
Definition: XclsTag.cxx:138

genie::Interaction::DISEM
static Interaction * DISEM(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:458

genie::Interaction::Create
static Interaction * Create(int tgt, int probe, ScatteringType_t st, InteractionType_t it)
Definition: Interaction.cxx:318

genie::kScQuasiElastic
Definition: ScatteringType.h:39

genie::kPdgDarkNeutrino
const int kPdgDarkNeutrino
Definition: PDGCodes.h:222

PDGCodes.h
Most commonly used PDG codes. A set of utility functions to handle PDG codes is provided in PDGUtils...

genie::InteractionType_t
enum genie::EInteractionType InteractionType_t

genie::Interaction::SetInitState
void SetInitState(const InitialState &init)
Definition: Interaction.cxx:225

genie::Interaction::AMNuGamma
static Interaction * AMNuGamma(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:819

genie::Interaction::SetProcInfo
void SetProcInfo(const ProcessInfo &proc)
Definition: Interaction.cxx:231

genie::ProcessInfo::Copy
void Copy(const ProcessInfo &proc)
Definition: ProcessInfo.cxx:252

genie::pdg::IsElectron
bool IsElectron(int pdgc)
Definition: PDGUtils.cxx:185

genie::kScSingleKaon
Definition: ScatteringType.h:40

endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030

genie::Interaction::fKinePhSp
KPhaseSpace * fKinePhSp
Kinematic phase space.
Definition: Interaction.h:183

genie::Interaction::NOsc
static Interaction * NOsc(int tgt, int annihilation_mode=-1)
Definition: Interaction.cxx:991

genie::Interaction::DISNC
static Interaction * DISNC(int tgt, int nuc, int probe, double E=0)
Definition: Interaction.cxx:393

genie::InitialState
Initial State information.
Definition: InitialState.h:48

pDEBUG
#define pDEBUG
Definition: Messenger.h:63

genie::Interaction::Interaction
Interaction()
Definition: Interaction.cxx:45