IBDInteractionListGenerator.cxx

Go to the documentation of this file.

//____________________________________________________________________________
/*
 Copyright (c) 2003-2025, The GENIE Collaboration
 For the full text of the license visit http://copyright.genie-mc.org
 
 
 Author: Corey Reed <cjreed \at nikhef.nl>
         Nikhef
 
 For the class documentation see the corresponding header file.
*/
//____________________________________________________________________________
 
#include "Framework/EventGen/InteractionList.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Physics/InverseBetaDecay/EventGen/IBDInteractionListGenerator.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
 
using namespace genie;
 
//___________________________________________________________________________
IBDInteractionListGenerator::IBDInteractionListGenerator() :
   InteractionListGeneratorI(
      "genie::IBDInteractionListGenerator")
{
 
}
//___________________________________________________________________________
IBDInteractionListGenerator::IBDInteractionListGenerator(
   string config):
   InteractionListGeneratorI(
      "genie::IBDInteractionListGenerator", config)
{
 
}
//___________________________________________________________________________
IBDInteractionListGenerator::~IBDInteractionListGenerator()
{
 
}
//___________________________________________________________________________
InteractionList * IBDInteractionListGenerator::CreateInteractionList(
   const InitialState & init_state) const
{
   LOG("IBD", pINFO) << "InitialState = " << init_state.AsString();
 
   const int  nupdg   = init_state.ProbePdg();
   const bool isnu    = pdg::IsNeutrino     (nupdg);
   const bool isnubar = pdg::IsAntiNeutrino (nupdg);
   
   const Target& target = init_state.Tgt();
   const int     tgtpdg = target.Pdg();
   const bool    hasP   = (target.Z() > 0);
   const bool    hasN   = (target.N() > 0);
   
   InteractionList * intlist = 0;
 
   if (isnu && hasN) {
      intlist = new InteractionList;
      intlist->push_back( Interaction::IBD(tgtpdg,kPdgNeutron, nupdg,0) );
   } else if (isnubar && hasP) {
      intlist = new InteractionList;
      intlist->push_back( Interaction::IBD(tgtpdg,kPdgProton,nupdg,0) );
   } else {
      LOG("IBD", pWARN)
         << "Returning NULL InteractionList for init-state: "
         << init_state.AsString();
   }
 
   return intlist;
}
//____________________________________________________________________________
void IBDInteractionListGenerator::Configure(const Registry & config)
{
   Algorithm::Configure(config);
   this->LoadConfigData();
}
//____________________________________________________________________________
void IBDInteractionListGenerator::Configure(string config)
{
   Algorithm::Configure(config);
   this->LoadConfigData();
}
//____________________________________________________________________________
void IBDInteractionListGenerator::LoadConfigData(void)
{
   // (no data to be loaded)
}
//____________________________________________________________________________