#include <G4PaulKxsModel.hh>

Public Member Functions
	G4PaulKxsModel ()

virtual	~G4PaulKxsModel ()

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

Detailed Description

Definition at line 51 of file G4PaulKxsModel.hh.

Constructor & Destructor Documentation

◆ G4PaulKxsModel()

G4PaulKxsModel::G4PaulKxsModel ( )

Definition at line 56 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 77 of file G4PaulKxsModel.cc.

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

Member Function Documentation

◆ CalculateKCrossSection()

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

Definition at line 86 of file G4PaulKxsModel.cc.

{
  
  G4Proton* aProtone = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();
  
  G4double sigma = 0;
 
  if (massIncident == aProtone->GetPDGMass() && zTarget < 93 && zTarget > 3)
    {
 
      //      G4EMDataSet* currentDataset =  protonDataSetMap[zTarget];
      //      currentDataset->GetEnergies
 
  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
    { 
      
// G4Exception("G4PaulKxsModel::CalculateKCrossSection()","de0004",JustWarning, "Energy deposited locally");
      sigma = 0.;
 
    }
    }
  
  
  // sigma is in internal units (mm^2)
  return sigma;
}

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