G4ecpssrFormFactorLixsModel Class Reference

#include <G4ecpssrFormFactorLixsModel.hh>

Inheritance diagram for G4ecpssrFormFactorLixsModel:

Collaboration diagram for G4ecpssrFormFactorLixsModel:

Public Member Functions
	G4ecpssrFormFactorLixsModel ()
virtual	~G4ecpssrFormFactorLixsModel ()
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 ()

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

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

Detailed Description

Definition at line 47 of file G4ecpssrFormFactorLixsModel.hh.

Constructor & Destructor Documentation

G4ecpssrFormFactorLixsModel() [1/2]

G4ecpssrFormFactorLixsModel::G4ecpssrFormFactorLixsModel

(

)

Definition at line 53 of file G4ecpssrFormFactorLixsModel.cc.

{ 
  interpolation = new G4LinInterpolation();

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

      protonL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL2DataSetMap[i]->LoadData("pixe/ecpssr/proton/l2-");

      protonL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      protonL3DataSetMap[i]->LoadData("pixe/ecpssr/proton/l3-");
  }

  for (G4int i=6; i<93; i++) 
  {
      alphaL1DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL1DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l1-");
         
      alphaL2DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL2DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l2-");
         
      alphaL3DataSetMap[i] = new G4EMDataSet(i,interpolation);
      alphaL3DataSetMap[i]->LoadData("pixe/ecpssr/alpha/l3-");
  }

}

~G4ecpssrFormFactorLixsModel()

G4ecpssrFormFactorLixsModel::~G4ecpssrFormFactorLixsModel ( )

virtual

Definition at line 85 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 &  )

private

Member Function Documentation

CalculateL1CrossSection()

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

virtual

Implements G4VecpssrLiModel.

Definition at line 101 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 > 5) {

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

Here is the call graph for this function:

CalculateL2CrossSection()

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

virtual

Implements G4VecpssrLiModel.

Definition at line 132 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 > 5) {

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

Here is the call graph for this function:

CalculateL3CrossSection()

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

virtual

Implements G4VecpssrLiModel.

Definition at line 163 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 > 5) {

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

Here is the call graph for this function:

operator=()

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

private

Member Data Documentation

alphaL1DataSetMap

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

private

Definition at line 69 of file G4ecpssrFormFactorLixsModel.hh.

alphaL2DataSetMap

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

private

Definition at line 70 of file G4ecpssrFormFactorLixsModel.hh.

alphaL3DataSetMap

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

private

Definition at line 71 of file G4ecpssrFormFactorLixsModel.hh.

interpolation

G4VDataSetAlgorithm* G4ecpssrFormFactorLixsModel::interpolation

private

Definition at line 64 of file G4ecpssrFormFactorLixsModel.hh.

protonL1DataSetMap

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

private

Definition at line 66 of file G4ecpssrFormFactorLixsModel.hh.

protonL2DataSetMap

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

private

Definition at line 67 of file G4ecpssrFormFactorLixsModel.hh.

protonL3DataSetMap

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

private

Definition at line 68 of file G4ecpssrFormFactorLixsModel.hh.

---

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

• G4ecpssrFormFactorLixsModel.hh
• G4ecpssrFormFactorLixsModel.cc