G4ecpssrFormFactorKxsModel Class Reference

#include <G4ecpssrFormFactorKxsModel.hh>

Inheritance diagram for G4ecpssrFormFactorKxsModel:

Public Member Functions
	G4ecpssrFormFactorKxsModel ()
virtual	~G4ecpssrFormFactorKxsModel ()
G4double	CalculateCrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)

---

Detailed Description

Definition at line 47 of file G4ecpssrFormFactorKxsModel.hh.

---

Constructor & Destructor Documentation

G4ecpssrFormFactorKxsModel::G4ecpssrFormFactorKxsModel

(

)

Definition at line 53 of file G4ecpssrFormFactorKxsModel.cc.

{ 
  interpolation = new G4LogLogInterpolation();

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

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

}

G4ecpssrFormFactorKxsModel::~G4ecpssrFormFactorKxsModel

(

)

[virtual]

Definition at line 73 of file G4ecpssrFormFactorKxsModel.cc.

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

---

Member Function Documentation

G4double G4ecpssrFormFactorKxsModel::CalculateCrossSection	(	G4int	zTarget,
		G4double	massIncident,
		G4double	energyIncident
	)			[virtual]

Implements G4VecpssrKModel.

Definition at line 82 of file G4ecpssrFormFactorKxsModel.cc.

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

{
  G4Proton* aProton = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();  
  G4double sigma = 0;

  if (energyIncident > 0.1*MeV && energyIncident < 100.*MeV && zTarget < 93 && zTarget > 5) {

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

---

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

• G4ecpssrFormFactorKxsModel.hh
• G4ecpssrFormFactorKxsModel.cc

---