G4ecpssrFormFactorLixsModel Class Reference

#include <G4ecpssrFormFactorLixsModel.hh>

Inheritance diagram for G4ecpssrFormFactorLixsModel:

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

Detailed Description

Definition at line 47 of file G4ecpssrFormFactorLixsModel.hh.

Constructor & Destructor Documentation

G4ecpssrFormFactorLixsModel() [1/2]

G4ecpssrFormFactorLixsModel::G4ecpssrFormFactorLixsModel ( )

explicit

Definition at line 51 of file G4ecpssrFormFactorLixsModel.cc.

{ 
  interpolation = new G4LinInterpolation();
 
  for (G4int i=11; i<93; i++) 
  {
      protonL1DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL1DataSetMap[i]->LoadData("pixe/ecpssr/proton/l1-i01m001c01-");
 
      protonL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL2DataSetMap[i]->LoadData("pixe/ecpssr/proton/l2-i01m001c01-");
 
      protonL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL3DataSetMap[i]->LoadData("pixe/ecpssr/proton/l3-i01m001c01-");
  }
 
  for (G4int i=11; i<93; i++) 
  {
      alphaL1DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL1DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l1-i02m004c02-");
         
      alphaL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL2DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l2-i02m004c02-");
         
      alphaL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL3DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l3-i02m004c02-");
  }
 
}

Referenced by G4ecpssrFormFactorLixsModel(), and operator=().

~G4ecpssrFormFactorLixsModel()

G4ecpssrFormFactorLixsModel::~G4ecpssrFormFactorLixsModel ( )

virtual

Definition at line 83 of file G4ecpssrFormFactorLixsModel.cc.

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

G4ecpssrFormFactorLixsModel() [2/2]

G4ecpssrFormFactorLixsModel::G4ecpssrFormFactorLixsModel ( const G4ecpssrFormFactorLixsModel & )

delete

Member Function Documentation

CalculateL1CrossSection()

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

overridevirtual

Implements G4VecpssrLiModel.

Definition at line 99 of file G4ecpssrFormFactorLixsModel.cc.

{
  G4Proton* aProton = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();  
  G4double sigma = 0;
 
  if (energyIncident > 0.1*MeV && energyIncident < 100.*MeV && zTarget < 93 && zTarget > 10) {
 
    if (massIncident == aProton->GetPDGMass())
      {      
    sigma = protonL1DataSetMap[zTarget]->FindValue(energyIncident/MeV);  
        if (sigma !=0 && energyIncident > protonL1DataSetMap[zTarget]->GetEnergies(0).back()*MeV) return 0.;
      }
    else if (massIncident == aAlpha->GetPDGMass())
      {
        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 G4ecpssrFormFactorLixsModel::CalculateL2CrossSection ( G4int zTarget, G4double massIncident, G4double energyIncident )

overridevirtual

Implements G4VecpssrLiModel.

Definition at line 130 of file G4ecpssrFormFactorLixsModel.cc.

{
  G4Proton* aProton = G4Proton::Proton();
  G4Alpha* aAlpha = G4Alpha::Alpha();  
  G4double sigma = 0;
 
  if (energyIncident > 0.1*MeV && energyIncident < 100.*MeV && zTarget < 93 && zTarget > 10) {
 
    if (massIncident == aProton->GetPDGMass())
      {      
    sigma = protonL2DataSetMap[zTarget]->FindValue(energyIncident/MeV);  
        if (sigma !=0 && energyIncident > protonL2DataSetMap[zTarget]->GetEnergies(0).back()*MeV) return 0.;
      }
    else if (massIncident == aAlpha->GetPDGMass())
      {
        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 G4ecpssrFormFactorLixsModel::CalculateL3CrossSection ( G4int zTarget, G4double massIncident, G4double energyIncident )

overridevirtual

Implements G4VecpssrLiModel.

Definition at line 161 of file G4ecpssrFormFactorLixsModel.cc.

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

operator=()

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

delete

---

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

• G4ecpssrFormFactorLixsModel.hh
• G4ecpssrFormFactorLixsModel.cc