#include <G4PaulKxsModel.hh>

Public Member Functions
	G4PaulKxsModel ()

virtual	~G4PaulKxsModel ()

G4double	CalculateKCrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)

	G4PaulKxsModel (const G4PaulKxsModel &)=delete

G4PaulKxsModel &	operator= (const G4PaulKxsModel &right)=delete

Detailed Description

Definition at line 50 of file G4PaulKxsModel.hh.

Constructor & Destructor Documentation

◆ G4PaulKxsModel() [1/2]

G4PaulKxsModel::G4PaulKxsModel ( )

explicit

Definition at line 57 of file G4PaulKxsModel.cc.

{   
  interpolation = new G4LogLogInterpolation();
 
  for (G4int i=4; i<93; i++) {
      protonDataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonDataSetMap[i]->LoadData("pixe/kpcsPaul/kcs-");
    }
    for (G4int i=6; i<93; i++) {
      alphaDataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaDataSetMap[i]->LoadData("pixe/kacsPaul/kacs-");
    }
}

◆ ~G4PaulKxsModel()

G4PaulKxsModel::~G4PaulKxsModel ( )

virtual

Definition at line 71 of file G4PaulKxsModel.cc.

{ 
  protonDataSetMap.clear();
  alphaDataSetMap.clear();
  delete interpolation;
}

◆ G4PaulKxsModel() [2/2]

G4PaulKxsModel::G4PaulKxsModel ( const G4PaulKxsModel & )

delete

Member Function Documentation

◆ CalculateKCrossSection()

G4double G4PaulKxsModel::CalculateKCrossSection	(	G4int	zTarget,
		G4double	massIncident,
		G4double	energyIncident
	)

Definition at line 79 of file G4PaulKxsModel.cc.

{
  
  G4Proton* aProtone = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();
  G4double sigma = 0;
 
  if (massIncident == aProtone->GetPDGMass() && zTarget < 93 && zTarget > 3)
    {
      if (energyIncident > protonDataSetMap[zTarget]->GetEnergies(0).back() ||
      energyIncident < protonDataSetMap[zTarget]->GetEnergies(0).front() )
    {sigma = 0;}
      else {     
    sigma = protonDataSetMap[zTarget]->FindValue(energyIncident/MeV); 
      }
    }
  else
    {
      if (massIncident == aAlpha->GetPDGMass() && zTarget < 93 && zTarget > 5)
    {
      if (energyIncident > alphaDataSetMap[zTarget]->GetEnergies(0).back() ||
          energyIncident < alphaDataSetMap[zTarget]->GetEnergies(0).front() )
        {sigma = 0;}
      else {
        sigma = alphaDataSetMap[zTarget]->FindValue(energyIncident/MeV); 
      }
    }
      else
    { 
      sigma = 0.;
    }
    }
  
  // sigma is in internal units (mm^2)
  return sigma;
}

Referenced by G4empCrossSection::CrossSection(), and G4empCrossSection::GetCrossSection().

◆ operator=()

G4PaulKxsModel & G4PaulKxsModel::operator= ( const G4PaulKxsModel & right )