G4IonYangFluctuationModel Class Reference

#include <G4IonYangFluctuationModel.hh>

Inheritance diagram for G4IonYangFluctuationModel:

Collaboration diagram for G4IonYangFluctuationModel:

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

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

)

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 G4Material *  material  ) const

virtual

Implements G4VLowEnergyModel.

Definition at line 107 of file G4IonYangFluctuationModel.cc.

{
  return 1.0*TeV ;
}

HighEnergyLimit() [2/2]

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

virtual

Implements G4VLowEnergyModel.

Definition at line 123 of file G4IonYangFluctuationModel.cc.

{
  return 1.0*TeV ;
}

IsInCharge() [1/2]

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

virtual

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

virtual

Implements G4VLowEnergyModel.

Definition at line 146 of file G4IonYangFluctuationModel.cc.

{
  return true ;
}

LowEnergyLimit() [1/2]

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

virtual

Implements G4VLowEnergyModel.

Definition at line 115 of file G4IonYangFluctuationModel.cc.

{
  return 0.0 ;
}

LowEnergyLimit() [2/2]

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

virtual

Implements G4VLowEnergyModel.

Definition at line 130 of file G4IonYangFluctuationModel.cc.

{
  return 0.0 ;
}

TheValue() [1/2]

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

virtual

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

Here is the call graph for this function:

TheValue() [2/2]

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

virtual

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

Here is the call graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4IonYangFluctuationModel.hh
• G4IonYangFluctuationModel.cc