#include <INukeOset.h>

Inheritance diagram for INukeOset:

Public Member Functions
	INukeOset ()
	contructor More...

virtual void	setupOset (const double &density, const double &pionTk, const int &pionPDG, const double &protonFraction)=0
	use to set up Oset class (assign pion Tk, nuclear density etc) More...

double	getTotalCrossSection () const
	return total = (qel+cex+abs) cross section More...

double	getCexCrossSection () const
	return cex cross section More...

double	getAbsorptionCrossSection () const
	return absorption cross section More...

double	getCexFraction () const
	return fraction of cex events More...

double	getAbsorptionFraction () const
	return fraction of absorption events More...

Protected Member Functions
virtual void	setCrossSections ()=0
	calculalte cross sections for each channel More...

void	setCrossSections (const int &pionPDG, const double &protonFraction)
	calculate avg cross sections according to proton / neutron fraction More...

Protected Attributes
double	fNuclearDensity
	nuclear density in fm-3 More...

double	fPionKineticEnergy
	pion kinetic energy in MeV More...

double	fTotalCrossSection
	el+cex+abs cross section (averaged over proton / neutron fraction) More...

double	fCexCrossSection
	cex cross section (averaged over proton / neutron fraction) More...

double	fAbsorptionCrossSection
	absorption cross section (averaged over proton / neutron fraction) More...

double	fQelCrossSections [fNChannels]
	total qel (el+cex) cross section for each channel More...

double	fCexCrossSections [fNChannels]
	cex cross section for each channel More...

Static Protected Attributes
static const unsigned int	fNChannels = 3
	number of possible channels: pi+n, pi+p, pi0 More...

Detailed Description

Definition at line 24 of file INukeOset.h.

Constructor & Destructor Documentation

INukeOset::INukeOset ( )

contructor

set up initial density and enegry values; set up pointer to current instance

Definition at line 8 of file INukeOset.cxx.

                         : fNuclearDensity (-1.0), fPionKineticEnergy (-1.0)
 {
   osetUtils::currentInstance = this;
 }

Member Function Documentation

double INukeOset::getAbsorptionCrossSection ( ) const

inline

return absorption cross section

Definition at line 47 of file INukeOset.h.

   {
     return fAbsorptionCrossSection;
   }

double INukeOset::getAbsorptionFraction ( ) const

inline

return fraction of absorption events

Definition at line 59 of file INukeOset.h.

Referenced by genie::HNIntranuke2018::HadronFateOset().

   {
     return fAbsorptionCrossSection / fTotalCrossSection;
   }

double INukeOset::getCexCrossSection ( ) const

inline

return cex cross section

Definition at line 41 of file INukeOset.h.

   {
     return fCexCrossSection;
   }

double INukeOset::getCexFraction ( ) const

inline

return fraction of cex events

Definition at line 53 of file INukeOset.h.

Referenced by genie::HNIntranuke2018::HadronFateOset().

   {
     return fCexCrossSection / fTotalCrossSection;
   }

double INukeOset::getTotalCrossSection ( ) const

inline

return total = (qel+cex+abs) cross section

Definition at line 35 of file INukeOset.h.

Referenced by genie::utils::intranuke2018::sigmaTotalOset().

   {
     return fTotalCrossSection;
   }

virtual void INukeOset::setCrossSections ( )

protectedpure virtual

calculalte cross sections for each channel

Implemented in INukeOsetFormula, and INukeOsetTable.

Referenced by INukeOsetFormula::setupOset(), and INukeOsetTable::setupOset().

void INukeOset::setCrossSections	(	const int &	pionPDG,
		const double &	protonFraction
	)

protected

calculate avg cross sections according to proton / neutron fraction

Definition at line 13 of file INukeOset.cxx.

References fAbsorptionCrossSection, fCexCrossSection, fCexCrossSections, fQelCrossSections, fTotalCrossSection, genie::kPdgPi0, genie::kPdgPiM, and genie::kPdgPiP.

 {  
   if (pionPDG == kPdgPi0)
   {
       fCexCrossSection = fCexCrossSections[2];
     fTotalCrossSection = fQelCrossSections[2] + fAbsorptionCrossSection;
   }
   else
   {
     // set channel for pi on proton and pi on neutron
     const int channelIndexOnProton  = (pionPDG == kPdgPiP); // 0 = pi-, 1 = pi+
     const int channelIndexOnNeutron = (pionPDG == kPdgPiM); // 0 = pi+, 1 = pi-
 
     // total xsec = (Z * xsec_proton + (A-Z) * xsec_neutron) / A
     fCexCrossSection = protonFraction * fCexCrossSections[channelIndexOnProton] +
                        (1.0 - protonFraction) * fCexCrossSections[channelIndexOnNeutron];
 
     fTotalCrossSection = protonFraction * fQelCrossSections[channelIndexOnProton] +
                         (1.0 - protonFraction) * fQelCrossSections[channelIndexOnNeutron] + fAbsorptionCrossSection;
   }
 }  