IBDXSecMap.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
 
 Corey Reed <cjreed \at nikhef.nl>
 Nikhef 
*/
//____________________________________________________________________________
 
#include "IBDXSecMap.h"
 
#include "Framework/Algorithm/AlgConfigPool.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
 
#include <sstream>
 
using namespace genie;
 
ClassImp(IBDXSecMap)
 
//____________________________________________________________________________
IBDXSecMap::IBDXSecMap() :
   XSecAlgorithmI("genie::IBDXSecMap"),
   fDefaultModel(0)
{
 
}
//____________________________________________________________________________
IBDXSecMap::IBDXSecMap(string config) :
   XSecAlgorithmI("genie::IBDXSecMap", config),
   fDefaultModel(0)
{
 
}
//____________________________________________________________________________
IBDXSecMap::~IBDXSecMap()
{
 
}
 
//____________________________________________________________________________
void IBDXSecMap::Configure(const Registry & config)
{
   Algorithm::Configure(config);
   this->LoadConfig();
}
//____________________________________________________________________________
void IBDXSecMap::Configure(string config)
{
   Algorithm::Configure(config);
   this->LoadConfig();
}
//____________________________________________________________________________
void IBDXSecMap::LoadConfig(void)
{
   // build the default xsec model according to the options contained
   // in IBDXSecMap.xml and/or UserPhysicsOptions.xml
 
   // load default global model (should work for all nuclei)
   RgAlg dgmodel ;
   GetParam( "IBDNucXSecModel", dgmodel ) ;
   LOG("IBD", pINFO)
      << "Default IBD cross section model: " << dgmodel;
 
   fDefaultModel =
      dynamic_cast<const XSecAlgorithmI*>( this -> SubAlg("IBDNucXSecModel") );
   assert(fDefaultModel!=0);
 
   // check whether to map according to specific isotopes
   GetParam("IsotopesUseSameModel", fIsotopesUseSameModel ) ;
 
   // load refined models for specific nuclei
   for(int Z=1; Z<140; Z++) {
     for(int A=Z; A<3*Z; A++) {
         std::ostringstream key;
         const int nucpdg = pdg::IonPdgCode(A,Z);
         key << "IBDNucXSecModel@Pdg=" << nucpdg;
         RgKey rgkey = key.str();
         if ( GetConfig().Exists(rgkey) )  {
                 RgAlg rgmodel = GetConfig().GetAlg(rgkey);
                 LOG("IBD", pNOTICE)
                 << "Nucleus =" << nucpdg
                         << " -> refined nuclear model: " << rgmodel;
                 const XSecAlgorithmI* model =
                                 dynamic_cast<const XSecAlgorithmI*>(this->SubAlg(rgkey));
                 assert(model);
                 const int mapkeyval = (fIsotopesUseSameModel) ? Z : nucpdg;
                 fRefinedModels.
                         insert(map<int,const XSecAlgorithmI*>::value_type(mapkeyval,
                                         model));
         }
     }
   }
}
//____________________________________________________________________________
const XSecAlgorithmI* IBDXSecMap::SelectModel(const Target & t) const
{
   // search the map for the PDG code of the target
   // if a refined model is found, return it
   // otherwise return the default xsec model
 
   map<int,const XSecAlgorithmI*>::const_iterator it = fRefinedModels.find(t.Pdg());
 
   if(it != fRefinedModels.end()) return it->second;
   else return fDefaultModel;
}
//____________________________________________________________________________
double IBDXSecMap::XSec(const Interaction * i, KinePhaseSpace_t k) const
{
   const XSecAlgorithmI* xs = this->SelectModel(i->InitState().Tgt());
   if (xs!=0) {
      return xs->XSec(i, k);
   } else {
      LOG("IBD", pERROR) << "No IBD XSec model found for target "
                         << i->InitState().TgtPdg();
      return 0;
   }
}
//____________________________________________________________________________
double IBDXSecMap::Integral(const Interaction * i) const
{
   const XSecAlgorithmI* xs = this->SelectModel(i->InitState().Tgt());
   if (xs!=0) {
      return xs->Integral(i);
   } else {
      LOG("IBD", pERROR) << "No IBD XSec model found for target "
                         << i->InitState().TgtPdg();
      return 0;
   }
}
//____________________________________________________________________________
bool   IBDXSecMap::ValidProcess(const Interaction * i) const
{
   const XSecAlgorithmI* xs = this->SelectModel(i->InitState().Tgt());
   if (xs!=0) {
      return xs->ValidProcess(i);
   } else {
      LOG("IBD", pERROR) << "No IBD XSec model found for target "
                         << i->InitState().TgtPdg();
      return 0;
   }
}
//____________________________________________________________________________
bool   IBDXSecMap::ValidKinematics(const Interaction * i) const
{
   const XSecAlgorithmI* xs = this->SelectModel(i->InitState().Tgt());
   if (xs!=0) {
      return xs->ValidKinematics(i);
   } else {
      LOG("IBD", pERROR) << "No IBD XSec model found for target "
                         << i->InitState().TgtPdg();
      return 0;
   }
}