G4NIELCalculator Class Reference

#include <G4NIELCalculator.hh>

Public Member Functions
	G4NIELCalculator (G4VEmModel *, G4int verb)
	~G4NIELCalculator ()
void	Initialise ()
G4double	ComputeNIEL (const G4Step *)
G4double	RecoilEnergy (const G4Step *)
void	AddEmModel (G4VEmModel *)

Detailed Description

Definition at line 60 of file G4NIELCalculator.hh.

Constructor & Destructor Documentation

G4NIELCalculator()

G4NIELCalculator::G4NIELCalculator ( G4VEmModel *  mod, G4int  verb  )

explicit

Definition at line 60 of file G4NIELCalculator.cc.

  : fModel(mod), fVerbose(verb)
{
  G4LossTableManager::Instance()->SetNIELCalculator(this);
  if(fVerbose > 0) {
    G4cout << "G4NIELCalculator: is created with the model <" 
           << fModel->GetName() << ">" << G4endl;
  }
}

~G4NIELCalculator()

G4NIELCalculator::~G4NIELCalculator

(

)

Definition at line 72 of file G4NIELCalculator.cc.

{}

Member Function Documentation

AddEmModel()

void G4NIELCalculator::AddEmModel

(

G4VEmModel *

mod

)

Definition at line 77 of file G4NIELCalculator.cc.

{
  if(mod && mod != fModel) { 
    fModel = mod; 
    if(fVerbose > 0) {
      G4cout << "G4NIELCalculator: new model <" << fModel->GetName()
             << "> is added" << G4endl;
    }
  }
}

ComputeNIEL()

G4double G4NIELCalculator::ComputeNIEL

(

const G4Step *

step

)

Definition at line 95 of file G4NIELCalculator.cc.

{
  G4double niel = 0.0;
  G4double T2 = step->GetPostStepPoint()->GetKineticEnergy();
  if(fModel && T2 > 0.) {
    const G4Track* track = step->GetTrack();
    const G4ParticleDefinition* part = track->GetParticleDefinition();
    G4double length = step->GetStepLength(); 
 
    if(length > 0.0 && part->GetPDGMass() > 100*CLHEP::MeV) {
 
      // primary
      G4double T1= step->GetPreStepPoint()->GetKineticEnergy();
      G4double T = 0.5*(T1 + T2);
      const G4MaterialCutsCouple* couple = 
    step->GetPreStepPoint()->GetMaterialCutsCouple(); 
      niel = length*fModel->ComputeDEDXPerVolume(couple->GetMaterial(),part,T);
      niel = std::min(niel, T1);
    }
  }
  return niel;
}

Initialise()

void G4NIELCalculator::Initialise

(

)

Definition at line 90 of file G4NIELCalculator.cc.

{}

Referenced by G4LossTableManager::BuildPhysicsTable().

RecoilEnergy()

G4double G4NIELCalculator::RecoilEnergy

(

const G4Step *

step

)

Definition at line 120 of file G4NIELCalculator.cc.

{
  G4double erec = 0.0;
  const std::vector<const G4Track*>* sec = step->GetSecondaryInCurrentStep(); 
 
  if(sec) {
    for(auto track : *sec) {
      const G4ParticleDefinition* part = track->GetParticleDefinition();
      if(part->IsGeneralIon()) {
        erec += track->GetKineticEnergy();
      }
    }
  }
  return erec;
}

---

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

• G4NIELCalculator.hh
• G4NIELCalculator.cc