IBDXSecMap.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

 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;
   }
}