#include <G4ecpssrFormFactorKxsModel.hh>

Inheritance diagram for G4ecpssrFormFactorKxsModel:

Public Member Functions
	G4ecpssrFormFactorKxsModel ()

virtual	~G4ecpssrFormFactorKxsModel ()

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

	G4ecpssrFormFactorKxsModel (const G4ecpssrFormFactorKxsModel &)=delete

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

Public Member Functions inherited from G4VecpssrKModel
	G4VecpssrKModel ()

virtual	~G4VecpssrKModel ()

	G4VecpssrKModel (const G4VecpssrKModel &)=delete

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

Detailed Description

Definition at line 47 of file G4ecpssrFormFactorKxsModel.hh.

Constructor & Destructor Documentation

◆ G4ecpssrFormFactorKxsModel() [1/2]

G4ecpssrFormFactorKxsModel::G4ecpssrFormFactorKxsModel ( )

explicit

Definition at line 53 of file G4ecpssrFormFactorKxsModel.cc.

{ 
  interpolation = new G4LogLogInterpolation();
 
  for (G4int i=3; i<93; i++) 
  {
      protonDataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonDataSetMap[i]->LoadData("pixe/ecpssr/proton/k-i01m001c01-");
  }
 
  for (G4int i=3; i<93; i++) 
  {
      alphaDataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaDataSetMap[i]->LoadData("pixe/ecpssr/alpha/k-i02m004c02-");
  }
}

◆ ~G4ecpssrFormFactorKxsModel()

G4ecpssrFormFactorKxsModel::~G4ecpssrFormFactorKxsModel ( )

virtual

Definition at line 72 of file G4ecpssrFormFactorKxsModel.cc.

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

◆ G4ecpssrFormFactorKxsModel() [2/2]

G4ecpssrFormFactorKxsModel::G4ecpssrFormFactorKxsModel ( const G4ecpssrFormFactorKxsModel & )

delete

Member Function Documentation

◆ CalculateCrossSection()

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

overridevirtual

Implements G4VecpssrKModel.

Definition at line 81 of file G4ecpssrFormFactorKxsModel.cc.

{
  G4Proton* aProton = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();  
  G4double sigma = 0;
 
  if (energyIncident > 0.1*MeV && energyIncident < 100.*MeV && zTarget < 93 && zTarget > 2) {
 
    if (massIncident == aProton->GetPDGMass())
      {         
    sigma = protonDataSetMap[zTarget]->FindValue(energyIncident/MeV);  
    if (sigma !=0 && energyIncident > protonDataSetMap[zTarget]->GetEnergies(0).back()*MeV) return 0.;
      }
    else if (massIncident == aAlpha->GetPDGMass())
      {
    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;
}

◆ operator=()

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