G4ANSTOecpssrLixsModel Class Reference

#include <G4ANSTOecpssrLixsModel.hh>

Inheritance diagram for G4ANSTOecpssrLixsModel:

Public Member Functions
	G4ANSTOecpssrLixsModel ()
virtual	~G4ANSTOecpssrLixsModel ()
G4double	CalculateL1CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
G4double	CalculateL2CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
G4double	CalculateL3CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)
Public Member Functions inherited from G4VecpssrLiModel
	G4VecpssrLiModel ()
virtual	~G4VecpssrLiModel ()
virtual G4double	CalculateL1CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)=0
virtual G4double	CalculateL2CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)=0
virtual G4double	CalculateL3CrossSection (G4int zTarget, G4double massIncident, G4double energyIncident)=0
	G4VecpssrLiModel (const G4VecpssrLiModel &)=delete
G4VecpssrLiModel &	operator= (const G4VecpssrLiModel &right)=delete

Detailed Description

Definition at line 48 of file G4ANSTOecpssrLixsModel.hh.

Constructor & Destructor Documentation

G4ANSTOecpssrLixsModel()

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

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

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

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

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

---

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

• G4ANSTOecpssrLixsModel.hh
• G4ANSTOecpssrLixsModel.cc