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

 Original code contributed by J. Schwehr, D. Cherdack, R. Gran
 Substantial code refactorizations by the core GENIE group.
*/
//_________________________________________________________________________

#include "Framework/Algorithm/AlgConfigPool.h"
#include "Framework/Messenger/Messenger.h"
#include "Physics/Common/XSecLinearCombinations.h"

using namespace genie;

//_________________________________________________________________________
XSecLinearCombinations::XSecLinearCombinations() :
  XSecAlgorithmI("genie::XSecLinearCombinations")
{

}
//_________________________________________________________________________
XSecLinearCombinations::XSecLinearCombinations(string config) :
  XSecAlgorithmI("genie::XSecLinearCombinations", config)
{

}
//_________________________________________________________________________
XSecLinearCombinations::~XSecLinearCombinations()
{

}

//_________________________________________________________________________

double XSecLinearCombinations::XSec(
  const Interaction * interaction, KinePhaseSpace_t kps) const
{
  double xsec = 0 ; 
  // Store xsec prediction for each entry :
  for( unsigned int i = 0 ; i < fXSections.size() ; ++ i ) {
    if( fXSections[i] -> ValidKinematics( interaction ) ) { 
      xsec += ( fXSections[i] -> XSec( interaction, kps ) ) * fLinearCoefficients[i] ;
    }
  }

  return xsec;
}

//_________________________________________________________________________

double XSecLinearCombinations::Integral( const Interaction * interaction ) const
{

  double xsec = 0 ; 
  // Store xsec prediction for each entry :
  for( unsigned int i = 0 ; i < fXSections.size() ; ++ i ) {
    xsec += ( fXSections[i] -> Integral( interaction ) ) * fLinearCoefficients[i] ; 
  }

  return xsec;
}

//_________________________________________________________________________

bool XSecLinearCombinations::ValidProcess( const Interaction * interaction) const
{
  bool is_valid = true ;

  // CHECK THAT THE XSEC USED CONSIDER A VALID PROCESS
  for( unsigned int i = 0 ; i < fXSections.size() ; ++ i ) {
    if( ! fXSections[i] -> ValidProcess( interaction ) ) {
      is_valid = false ;
      break ; 
    }
  }
  
  return is_valid ;
}

//_________________________________________________________________________

void XSecLinearCombinations::Configure(const Registry & config)
{
    Algorithm::Configure(config);
    this->LoadConfig();
}

//____________________________________________________________________________

void XSecLinearCombinations::Configure(string config)
{
    Algorithm::Configure(config);
    this->LoadConfig();
}

//_________________________________________________________________________

void XSecLinearCombinations::LoadConfig(void)
{
  bool good_config = true ;
  GetParamVect( "LinearCoefficients", fLinearCoefficients ) ; 

  bool Normalise = false ;
  GetParam( "Normalise", Normalise ) ;

  if( Normalise ) {
    double sum = 0 ; 
    for( unsigned int i = 0 ; i < fLinearCoefficients.size() ; ++i ) {
      sum += fLinearCoefficients[i] ;
    }

    if( sum == 0 ) {
      good_config = false ;
      LOG("XSecLinearCombinations", pERROR) << "Sum of linear coefficients is zero. Cannot normalize " ;
    }

    for( unsigned int j = 0 ; j < fLinearCoefficients.size() ; ++j ) {
      fLinearCoefficients[j] /= sum ;
    }
  }

  std::vector<RgKey> XSectionsKeys ;
  GetParamVectKeys( "CrossSection", XSectionsKeys ) ;

  if( XSectionsKeys.size() != fLinearCoefficients.size() ) {
    good_config = false ; 
    LOG("XSecLinearCombinations", pERROR) << "Entries don't match" ;
    LOG("XSecLinearCombinations", pERROR) << "Cross Sections size: " << XSectionsKeys.size() ;
    LOG("XSecLinearCombinations", pERROR) << "Linear Coefficients size: " << fLinearCoefficients.size() ;
  }

  for( unsigned int i = 0 ; i < XSectionsKeys.size() ; ++i ) {
    fXSections.push_back( dynamic_cast<const XSecAlgorithmI *> ( this->SubAlg( XSectionsKeys[i] ) ) );
    if ( ! fXSections[i] ) {
      good_config = false ;
      LOG("XSecLinearCombinations", pERROR) << "The SubAlg is not a XSecAlgorithmI";
    }
  }

  if( ! good_config ) {
    LOG("XSecLinearCombinations", pERROR) << "XSecLinearCombilation Configuration has failed.";
    exit(78) ;
  }

}