G4ANSTOecpssrLixsModel Class Reference

#include <G4ANSTOecpssrLixsModel.hh>

Inheritance diagram for G4ANSTOecpssrLixsModel:

Public Member Functions
G4double	CalculateL1CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
G4double	CalculateL2CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
G4double	CalculateL3CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
	G4ANSTOecpssrLixsModel ()
virtual	~G4ANSTOecpssrLixsModel ()

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

Private Attributes
std::map< G4int, G4VEMDataSet * >	alphaL1DataSetMap
std::map< G4int, G4VEMDataSet * >	alphaL2DataSetMap
std::map< G4int, G4VEMDataSet * >	alphaL3DataSetMap
std::map< G4int, G4VEMDataSet * >	carbonL1DataSetMap
std::map< G4int, G4VEMDataSet * >	carbonL2DataSetMap
std::map< G4int, G4VEMDataSet * >	carbonL3DataSetMap
G4VDataSetAlgorithm *	interpolation
std::map< G4int, G4VEMDataSet * >	protonL1DataSetMap
std::map< G4int, G4VEMDataSet * >	protonL2DataSetMap
std::map< G4int, G4VEMDataSet * >	protonL3DataSetMap

Detailed Description

Definition at line 48 of file G4ANSTOecpssrLixsModel.hh.

Constructor & Destructor Documentation

G4ANSTOecpssrLixsModel() [1/2]

G4ANSTOecpssrLixsModel::G4ANSTOecpssrLixsModel

(

)

Definition at line 54 of file G4ANSTOecpssrLixsModel.cc.

{
  G4cout << "Using ANSTO L Cross Sections! "<< G4endl;
 
  interpolation = new G4LinInterpolation();
 
  for (G4int i=26; i<93; i++)
  {
      protonL1DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL1DataSetMap[i]->LoadData("pixe_ANSTO/proton/l1-");
 
      protonL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL2DataSetMap[i]->LoadData("pixe_ANSTO/proton/l2-");
 
      protonL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL3DataSetMap[i]->LoadData("pixe_ANSTO/proton/l3-");
  }
 
  for (G4int i=26; i<93; i++)
  {
      alphaL1DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL1DataSetMap[i]->LoadData("pixe_ANSTO/alpha/l1-");
 
      alphaL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL2DataSetMap[i]->LoadData("pixe_ANSTO/alpha/l2-");
 
      alphaL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL3DataSetMap[i]->LoadData("pixe_ANSTO/alpha/l3-");
  }
 
}

References alphaL1DataSetMap, alphaL2DataSetMap, alphaL3DataSetMap, G4cout, G4endl, interpolation, protonL1DataSetMap, protonL2DataSetMap, and protonL3DataSetMap.

~G4ANSTOecpssrLixsModel()

G4ANSTOecpssrLixsModel::~G4ANSTOecpssrLixsModel ( )

virtual

Definition at line 88 of file G4ANSTOecpssrLixsModel.cc.

{
  protonL1DataSetMap.clear();
  alphaL1DataSetMap.clear();
 
  protonL2DataSetMap.clear();
  alphaL2DataSetMap.clear();
 
  protonL3DataSetMap.clear();
  alphaL3DataSetMap.clear();
 
  delete interpolation;
}

References alphaL1DataSetMap, alphaL2DataSetMap, alphaL3DataSetMap, interpolation, protonL1DataSetMap, protonL2DataSetMap, and protonL3DataSetMap.

G4ANSTOecpssrLixsModel() [2/2]

G4ANSTOecpssrLixsModel::G4ANSTOecpssrLixsModel ( const G4ANSTOecpssrLixsModel &  )

private

Member Function Documentation

CalculateL1CrossSection()

G4double G4ANSTOecpssrLixsModel::CalculateL1CrossSection ( G4int  zTarget, G4double  massIncident, G4double  energyIncident  )

virtual

Implements G4VecpssrLiModel.

Definition at line 104 of file G4ANSTOecpssrLixsModel.cc.

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

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

CalculateL2CrossSection()

G4double G4ANSTOecpssrLixsModel::CalculateL2CrossSection ( G4int  zTarget, G4double  massIncident, G4double  energyIncident  )

virtual

Implements G4VecpssrLiModel.

Definition at line 140 of file G4ANSTOecpssrLixsModel.cc.

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

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

CalculateL3CrossSection()

G4double G4ANSTOecpssrLixsModel::CalculateL3CrossSection ( G4int  zTarget, G4double  massIncident, G4double  energyIncident  )

virtual

Implements G4VecpssrLiModel.

Definition at line 174 of file G4ANSTOecpssrLixsModel.cc.

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

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

operator=()

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

private

Field Documentation

alphaL1DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::alphaL1DataSetMap

private

Definition at line 70 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL1CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

alphaL2DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::alphaL2DataSetMap

private

Definition at line 71 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL2CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

alphaL3DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::alphaL3DataSetMap

private

Definition at line 72 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL3CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

carbonL1DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::carbonL1DataSetMap

private

Definition at line 73 of file G4ANSTOecpssrLixsModel.hh.

carbonL2DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::carbonL2DataSetMap

private

Definition at line 74 of file G4ANSTOecpssrLixsModel.hh.

carbonL3DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::carbonL3DataSetMap

private

Definition at line 75 of file G4ANSTOecpssrLixsModel.hh.

interpolation

G4VDataSetAlgorithm* G4ANSTOecpssrLixsModel::interpolation

private

Definition at line 65 of file G4ANSTOecpssrLixsModel.hh.

Referenced by G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

protonL1DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::protonL1DataSetMap

private

Definition at line 67 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL1CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

protonL2DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::protonL2DataSetMap

private

Definition at line 68 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL2CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

protonL3DataSetMap

std::map< G4int , G4VEMDataSet* > G4ANSTOecpssrLixsModel::protonL3DataSetMap

private

Definition at line 69 of file G4ANSTOecpssrLixsModel.hh.

Referenced by CalculateL3CrossSection(), G4ANSTOecpssrLixsModel(), and ~G4ANSTOecpssrLixsModel().

---

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

• source/processes/electromagnetic/lowenergy/include/G4ANSTOecpssrLixsModel.hh
• source/processes/electromagnetic/lowenergy/src/G4ANSTOecpssrLixsModel.cc