#include <G4ANSTOecpssrKxsModel.hh>

Inheritance diagram for G4ANSTOecpssrKxsModel:

Public Member Functions
	G4ANSTOecpssrKxsModel ()

virtual	~G4ANSTOecpssrKxsModel ()

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

Public Member Functions inherited from G4VecpssrKModel
	G4VecpssrKModel ()

virtual	~G4VecpssrKModel ()

	G4VecpssrKModel (const G4VecpssrKModel &)=delete

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

Detailed Description

Definition at line 48 of file G4ANSTOecpssrKxsModel.hh.

Constructor & Destructor Documentation

◆ G4ANSTOecpssrKxsModel()

G4ANSTOecpssrKxsModel::G4ANSTOecpssrKxsModel ( )

Definition at line 54 of file G4ANSTOecpssrKxsModel.cc.

{
 
  G4cout << "Using ANSTO K Cross Sections! "<< G4endl;
 
  interpolation = new G4LogLogInterpolation();
  for (G4int i=6; i<93; i++)
  {
      protonDataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonDataSetMap[i]->LoadData("pixe_ANSTO/proton/k-");
  }
 
  for (G4int i=6; i<93; i++)
  {
      alphaDataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaDataSetMap[i]->LoadData("pixe_ANSTO/alpha/k-");
  }
 
}

◆ ~G4ANSTOecpssrKxsModel()

G4ANSTOecpssrKxsModel::~G4ANSTOecpssrKxsModel ( )

virtual

Definition at line 76 of file G4ANSTOecpssrKxsModel.cc.

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

Member Function Documentation

◆ CalculateCrossSection()

G4double G4ANSTOecpssrKxsModel::CalculateCrossSection	(	G4int	zTarget,
		G4double	massIncident,
		G4double	energyIncident )

virtual

Implements G4VecpssrKModel.

Definition at line 85 of file G4ANSTOecpssrKxsModel.cc.

{
  G4Proton* aProton = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();
  G4double sigma = 0;
 
  if (massIncident == aProton->GetPDGMass())
    {
      if (energyIncident > 0.2*MeV && energyIncident < 5.*MeV && zTarget < 93 && zTarget > 5) {
 
            sigma = protonDataSetMap[zTarget]->FindValue(energyIncident/MeV);
            if (sigma !=0 && energyIncident > protonDataSetMap[zTarget]->GetEnergies(0).back()*MeV) return 0.;
        }
    }
    
  else if (massIncident == aAlpha->GetPDGMass())
        {
          if (energyIncident > 0.2*MeV && energyIncident < 40.*MeV && zTarget < 93 && zTarget > 5) {
 
            sigma = alphaDataSetMap[zTarget]->FindValue(energyIncident/MeV);
            if (sigma !=0 && energyIncident > alphaDataSetMap[zTarget]->GetEnergies(0).back()*MeV) return 0.;
          }
         }
          
   else
          {
              sigma = 0.;
          }
        
  // sigma is in internal units: it has been converted from
  // the input file in barns bt the EmDataset
  return sigma;
}