GLRESInteractionListGenerator.cxx

Go to the documentation of this file.

//____________________________________________________________________________
/*
 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
*/
//____________________________________________________________________________

#include "Framework/EventGen/InteractionList.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Physics/GlashowResonance/EventGen/GLRESInteractionListGenerator.h"

using namespace genie;

//___________________________________________________________________________
GLRESInteractionListGenerator::GLRESInteractionListGenerator() :
InteractionListGeneratorI("genie::GLRESInteractionListGenerator")
{

}
//___________________________________________________________________________
GLRESInteractionListGenerator::GLRESInteractionListGenerator(string config) :
InteractionListGeneratorI("genie::GLRESInteractionListGenerator", config)
{

}
//___________________________________________________________________________
GLRESInteractionListGenerator::~GLRESInteractionListGenerator()
{

}
//___________________________________________________________________________
InteractionList *
   GLRESInteractionListGenerator::CreateInteractionList(
                                       const InitialState & init_state) const
{
// channels:
// nuebar + e- -> W- -> nuebar + e-
// nuebar + e- -> W- -> nuebar + mu-
// nuebar + e- -> W- -> nuebar + tau-
// nuebar + e- -> W- -> hadrons

  if(init_state.ProbePdg() != kPdgAntiNuE) {
     LOG("IntLst", pDEBUG)
          << "Return *null* interaction list";
     return 0;
  }

  InitialState init(init_state);
  init_state.TgtPtr()->SetHitNucPdg(0);  

  ProcessInfo   proc_info(kScGlashowResonance, kIntWeakCC);

  InteractionList * intlist = new InteractionList;

  if (fIsMu) {
    Interaction * interaction = new Interaction(init_state, proc_info);
    XclsTag exclusive_tag;
    exclusive_tag.SetFinalLepton(kPdgMuon);
    interaction->SetExclTag(exclusive_tag);
    intlist->push_back(interaction);  
  }
  else if (fIsTau) {
    Interaction * interaction = new Interaction(init_state, proc_info);
    XclsTag exclusive_tag;
    exclusive_tag.SetFinalLepton(kPdgTau);
    interaction->SetExclTag(exclusive_tag);
    intlist->push_back(interaction);  
  }
  else if (fIsEle) {
    Interaction * interaction = new Interaction(init_state, proc_info);
    XclsTag exclusive_tag;
    exclusive_tag.SetFinalLepton(kPdgElectron);
    interaction->SetExclTag(exclusive_tag);
    intlist->push_back(interaction);  
  }
  else if (fIsHad) {
    Interaction * interaction = new Interaction(init_state, proc_info);
    XclsTag exclusive_tag;
    exclusive_tag.SetFinalLepton(kPdgPiP);
    interaction->SetExclTag(exclusive_tag);
    intlist->push_back(interaction);  
  }

  if(intlist->size() == 0) {
     LOG("IntLst", pERROR)
         << "Returning NULL InteractionList for init-state: "
                                                  << init_state.AsString();
     delete intlist;
     return 0;
  }

  return intlist;
}
//___________________________________________________________________________
void GLRESInteractionListGenerator::Configure(const Registry & config)
{
  Algorithm::Configure(config);
  this->LoadConfigData();
}
//____________________________________________________________________________
void GLRESInteractionListGenerator::Configure(string config)
{
  Algorithm::Configure(config);
  this->LoadConfigData();
}
//____________________________________________________________________________
void GLRESInteractionListGenerator::LoadConfigData(void)
{

  GetParamDef("is-Mu",  fIsMu,  false ) ;
  GetParamDef("is-Tau", fIsTau, false ) ;
  GetParamDef("is-Ele", fIsEle, false ) ;
  GetParamDef("is-Had", fIsHad, false ) ;

}