G4PaulKxsModel Class Reference

#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

(

)

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 ( )

virtual

Definition at line 77 of file G4PaulKxsModel.cc.

{ 

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

}

Member Function Documentation

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

Definition at line 86 of file G4PaulKxsModel.cc.

References G4Alpha::Alpha(), G4ParticleDefinition::GetPDGMass(), python.hepunit::MeV, and G4Proton::Proton().

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

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

---

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

• G4PaulKxsModel.hh
• G4PaulKxsModel.cc