G4PaulKxsModel Class Reference

#include <G4PaulKxsModel.hh>

Collaboration diagram for G4PaulKxsModel:

Public Member Functions
	G4PaulKxsModel ()
virtual	~G4PaulKxsModel ()
G4double	CalculateKCrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)

Private Member Functions
	G4PaulKxsModel (const G4PaulKxsModel &)
G4PaulKxsModel &	operator= (const G4PaulKxsModel &right)

Private Attributes
G4VDataSetAlgorithm *	interpolation
std::map< G4int, G4VEMDataSet *>	protonDataSetMap
std::map< G4int, G4VEMDataSet *>	alphaDataSetMap

Detailed Description

Definition at line 51 of file G4PaulKxsModel.hh.

Constructor & Destructor Documentation

G4PaulKxsModel() [1/2]

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;

}

G4PaulKxsModel() [2/2]

G4PaulKxsModel::G4PaulKxsModel ( const G4PaulKxsModel &  )

private

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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

Member Data Documentation

alphaDataSetMap

std::map< G4int , G4VEMDataSet* > G4PaulKxsModel::alphaDataSetMap

private

Definition at line 73 of file G4PaulKxsModel.hh.

interpolation

G4VDataSetAlgorithm* G4PaulKxsModel::interpolation

private

Definition at line 69 of file G4PaulKxsModel.hh.

protonDataSetMap

std::map< G4int , G4VEMDataSet* > G4PaulKxsModel::protonDataSetMap

private

Definition at line 71 of file G4PaulKxsModel.hh.

---

The documentation for this class was generated from the following files:

• G4PaulKxsModel.hh
• G4PaulKxsModel.cc