G4IonYangFluctuationModel Class Reference

#include <G4IonYangFluctuationModel.hh>

Inheritance diagram for G4IonYangFluctuationModel:

Public Member Functions
	G4IonYangFluctuationModel (const G4String &name)
G4double	HighEnergyLimit (const G4ParticleDefinition *aParticle) const override
G4double	HighEnergyLimit (const G4ParticleDefinition *aParticle, const G4Material *material) const override
G4bool	IsInCharge (const G4DynamicParticle *particle, const G4Material *material) const override
G4bool	IsInCharge (const G4ParticleDefinition *aParticle, const G4Material *material) const override
G4double	LowEnergyLimit (const G4ParticleDefinition *aParticle) const override
G4double	LowEnergyLimit (const G4ParticleDefinition *aParticle, const G4Material *material) const override
G4double	TheValue (const G4DynamicParticle *particle, const G4Material *material) override
G4double	TheValue (const G4ParticleDefinition *aParticle, const G4Material *material, G4double kineticEnergy) override
	~G4IonYangFluctuationModel ()

Private Member Functions
G4double	YangFluctuationModel (const G4Material *material, G4double kineticEnergy, G4double particleMass, G4double charge) const

Detailed Description

Definition at line 56 of file G4IonYangFluctuationModel.hh.

Constructor & Destructor Documentation

G4IonYangFluctuationModel()

G4IonYangFluctuationModel::G4IonYangFluctuationModel ( const G4String &  name )

explicit

Definition at line 68 of file G4IonYangFluctuationModel.cc.

  : G4VLowEnergyModel(name)
{;}

~G4IonYangFluctuationModel()

G4IonYangFluctuationModel::~G4IonYangFluctuationModel

(

)

Definition at line 74 of file G4IonYangFluctuationModel.cc.

{;}

Member Function Documentation

HighEnergyLimit() [1/2]

G4double G4IonYangFluctuationModel::HighEnergyLimit ( const G4ParticleDefinition *  aParticle ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 123 of file G4IonYangFluctuationModel.cc.

{
  return 1.0*TeV ;
}

References TeV.

HighEnergyLimit() [2/2]

G4double G4IonYangFluctuationModel::HighEnergyLimit ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 107 of file G4IonYangFluctuationModel.cc.

{
  return 1.0*TeV ;
}

References TeV.

IsInCharge() [1/2]

G4bool G4IonYangFluctuationModel::IsInCharge ( const G4DynamicParticle *  particle, const G4Material *  material  ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 138 of file G4IonYangFluctuationModel.cc.

{
  return true ;
}

IsInCharge() [2/2]

G4bool G4IonYangFluctuationModel::IsInCharge ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 146 of file G4IonYangFluctuationModel.cc.

{
  return true ;
}

LowEnergyLimit() [1/2]

G4double G4IonYangFluctuationModel::LowEnergyLimit ( const G4ParticleDefinition *  aParticle ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 130 of file G4IonYangFluctuationModel.cc.

{
  return 0.0 ;
}

LowEnergyLimit() [2/2]

G4double G4IonYangFluctuationModel::LowEnergyLimit ( const G4ParticleDefinition *  aParticle, const G4Material *  material  ) const

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 115 of file G4IonYangFluctuationModel.cc.

{
  return 0.0 ;
}

TheValue() [1/2]

G4double G4IonYangFluctuationModel::TheValue ( const G4DynamicParticle *  particle, const G4Material *  material  )

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 79 of file G4IonYangFluctuationModel.cc.

{
  G4double energy = particle->GetKineticEnergy() ;
  G4double mass = particle->GetMass() ;
  G4double charge = (particle->GetCharge())/eplus ;
 
  G4double q = YangFluctuationModel(material,energy,mass,charge) ;
 
  return q ;
}

References G4INCL::KinematicsUtils::energy(), eplus, G4DynamicParticle::GetCharge(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMass(), eplot::material, and YangFluctuationModel().

TheValue() [2/2]

G4double G4IonYangFluctuationModel::TheValue ( const G4ParticleDefinition *  aParticle, const G4Material *  material, G4double  kineticEnergy  )

overridevirtual

Implements G4VLowEnergyModel.

Definition at line 93 of file G4IonYangFluctuationModel.cc.

{
  G4double mass = aParticle->GetPDGMass() ;
  G4double charge = (aParticle->GetPDGCharge())/eplus ;
 
  G4double q = YangFluctuationModel(material,kineticEnergy,mass,charge);
 
  return q ;
}

References eplus, G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), eplot::material, and YangFluctuationModel().

YangFluctuationModel()

G4double G4IonYangFluctuationModel::YangFluctuationModel ( const G4Material *  material, G4double  kineticEnergy, G4double  particleMass, G4double  charge  ) const

private

Definition at line 154 of file G4IonYangFluctuationModel.cc.

{
  // The aproximation of energy loss fluctuations
  // Q.Yang et al., NIM B61(1991)149-155.
 
  // Reduced energy in MeV/AMU
  G4double energy = kineticEnergy *amu_c2/(particleMass*MeV) ;
 
  G4int i = 0 ;
  G4double factor = 1.0 ;
 
  // The index of set of parameters i = 0 for protons(hadrons) in gases
  //                                    1 for protons(hadrons) in solids
  //                                    2 for ions in atomic gases
  //                                    3 for ions in molecular gases
  //                                    4 for ions in solids
  static const G4double b[5][4] = {
  {0.1014,  0.3700,  0.9642,  3.987},
  {0.1955,  0.6941,  2.522,   1.040},
  {0.05058, 0.08975, 0.1419, 10.80},
  {0.05009, 0.08660, 0.2751,  3.787},
  {0.01273, 0.03458, 0.3951,  3.812}
  } ;
 
  // protons (hadrons)
  if(1.5 > charge) {
    if( kStateGas != material->GetState() ) i = 1 ;
 
  // ions
  } else {
    G4double zeff = (material->GetElectronDensity())/
                    (material->GetTotNbOfAtomsPerVolume()) ;
    factor = charge * std::pow(charge/zeff, 0.3333) ;
 
    if( kStateGas == material->GetState() ) {
      energy /= (charge * std::sqrt(charge)) ;
 
      if(1 == (material->GetNumberOfElements())) {
        i = 2 ;
      } else {
        i = 3 ;
      }
 
    } else {
      energy /= (charge * std::sqrt(charge*zeff)) ;
      i = 4 ;
    }
  }
 
  G4double x = b[i][2] * (1.0 - G4Exp( - energy * b[i][3] )) ;
 
  G4double q = factor * x * b[i][0] /
             ((energy - b[i][1])*(energy - b[i][1]) + x*x) ;
 
  return q ;
}

References source.hepunit::amu_c2, G4INCL::KinematicsUtils::energy(), G4Exp(), kStateGas, eplot::material, and MeV.

Referenced by TheValue().

---

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

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