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 ()

Detailed Description

Definition at line 47 of file G4ecpssrFormFactorLixsModel.hh.

Constructor & Destructor Documentation

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 ( )

virtual

Definition at line 85 of file G4ecpssrFormFactorLixsModel.cc.

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

Member Function Documentation

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:

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:

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:

---

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

• source/geant4.10.03.p02/source/processes/electromagnetic/lowenergy/include/G4ecpssrFormFactorLixsModel.hh
• source/geant4.10.03.p02/source/processes/electromagnetic/lowenergy/src/G4ecpssrFormFactorLixsModel.cc