G4ecpssrFormFactorKxsModel Class Reference

#include <G4ecpssrFormFactorKxsModel.hh>

Inheritance diagram for G4ecpssrFormFactorKxsModel:

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

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(), python.hepunit::MeV, 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