#include <G4NIELCalculator.hh>

Public Member Functions
void	AddEmModel (G4VEmModel *)

G4double	ComputeNIEL (const G4Step *)

	G4NIELCalculator (const G4NIELCalculator &)=delete

	G4NIELCalculator (G4VEmModel *, G4int verb)

void	Initialise ()

G4NIELCalculator &	operator= (const G4NIELCalculator &right)=delete

G4double	RecoilEnergy (const G4Step *)

	~G4NIELCalculator ()

Private Attributes
G4VEmModel *	fModel

G4int	fVerbose

Detailed Description

Definition at line 60 of file G4NIELCalculator.hh.

Constructor & Destructor Documentation

◆ G4NIELCalculator() [1/2]

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

References fModel, fVerbose, G4cout, G4endl, G4VEmModel::GetName(), G4LossTableManager::Instance(), and G4LossTableManager::SetNIELCalculator().

◆ ~G4NIELCalculator()

G4NIELCalculator::~G4NIELCalculator ( )

Definition at line 72 of file G4NIELCalculator.cc.

73{}

◆ G4NIELCalculator() [2/2]

G4NIELCalculator::G4NIELCalculator ( const G4NIELCalculator & )

delete

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

References fModel, fVerbose, G4cout, G4endl, and G4VEmModel::GetName().

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

References G4VEmModel::ComputeDEDXPerVolume(), fModel, G4StepPoint::GetKineticEnergy(), G4MaterialCutsCouple::GetMaterial(), G4StepPoint::GetMaterialCutsCouple(), G4Track::GetParticleDefinition(), G4ParticleDefinition::GetPDGMass(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4Step::GetStepLength(), G4Step::GetTrack(), CLHEP::MeV, and G4INCL::Math::min().

◆ Initialise()

void G4NIELCalculator::Initialise ( )

Definition at line 90 of file G4NIELCalculator.cc.

91{}

Referenced by G4LossTableManager::BuildPhysicsTable().

◆ operator=()

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

delete

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

References G4Step::GetSecondaryInCurrentStep(), and G4ParticleDefinition::IsGeneralIon().

Field Documentation

◆ fModel

G4VEmModel* G4NIELCalculator::fModel

private

Definition at line 86 of file G4NIELCalculator.hh.

Referenced by AddEmModel(), ComputeNIEL(), and G4NIELCalculator().

◆ fVerbose

G4int G4NIELCalculator::fVerbose