G4HadronStoppingProcess Class Reference

#include <G4HadronStoppingProcess.hh>

Inheritance diagram for G4HadronStoppingProcess:

Public Member Functions
	G4HadronStoppingProcess (const G4String &name="hadronCaptureAtRest")
virtual	~G4HadronStoppingProcess ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual void	PreparePhysicsTable (const G4ParticleDefinition &)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
virtual void	ProcessDescription (std::ostream &outFile) const
void	SetElementSelector (G4ElementSelector *ptr)
void	SetEmCascade (G4HadronicInteraction *ptr)
void	SetBoundDecay (G4HadronicInteraction *ptr)
Protected Member Functions
G4double	GetMeanLifeTime (const G4Track &, G4ForceCondition *)

---

Detailed Description

Definition at line 62 of file G4HadronStoppingProcess.hh.

---

Constructor & Destructor Documentation

G4HadronStoppingProcess::G4HadronStoppingProcess

(

const G4String &

name = "hadronCaptureAtRest"

)

[explicit]

Definition at line 63 of file G4HadronStoppingProcess.cc.

References G4VProcess::enableAtRestDoIt, G4VProcess::enablePostStepDoIt, G4HadronicProcessStore::Instance(), and G4HadronicProcessStore::RegisterExtraProcess().

  : G4HadronicProcess(name, fHadronAtRest)
{
  // Modify G4VProcess flags to emulate G4VRest instead of G4VDiscrete
  enableAtRestDoIt = true;
  enablePostStepDoIt = false;

  fElementSelector = new G4ElementSelector();
  fEmCascade = new G4EmCaptureCascade();        // Owned by InteractionRegistry
  fBoundDecay = 0;
  G4HadronicProcessStore::Instance()->RegisterExtraProcess(this);
}

G4HadronStoppingProcess::~G4HadronStoppingProcess

(

)

[virtual]

Definition at line 78 of file G4HadronStoppingProcess.cc.

References G4HadronicProcessStore::DeRegisterExtraProcess(), and G4HadronicProcessStore::Instance().

{
  G4HadronicProcessStore::Instance()->DeRegisterExtraProcess(this);
  delete fElementSelector;
  // NOTE: fEmCascade and fEmBoundDecay owned by registry, not locally
}

---

Member Function Documentation

G4VParticleChange * G4HadronStoppingProcess::AtRestDoIt	(	const G4Track &	,
		const G4Step &
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 126 of file G4HadronStoppingProcess.cc.

References G4HadFinalState::AddSecondaries(), G4ParticleChange::AddSecondary(), G4HadronicInteraction::ApplyYourself(), G4HadronicProcess::CheckEnergyMomentumConservation(), G4HadronicProcess::CheckResult(), G4HadronicProcess::ChooseHadronicInteraction(), G4HadFinalState::Clear(), G4HadronicProcess::DumpState(), G4HadronicProcess::epReportLevel, FatalException, fStopAndKill, G4endl, G4Exception(), G4Nucleus::GetA_asInt(), G4Track::GetGlobalTime(), G4HadFinalState::GetLocalEnergyDeposit(), G4Track::GetMaterial(), G4HadronicInteraction::GetModelName(), G4Element::GetName(), G4HadFinalState::GetNumberOfSecondaries(), G4Track::GetPosition(), G4HadFinalState::GetSecondary(), G4HadFinalState::GetStatusChange(), G4HadronicProcess::GetTargetNucleusPointer(), G4Track::GetTouchableHandle(), G4Track::GetWeight(), G4Nucleus::GetZ_asInt(), G4HadProjectile::Initialise(), G4ParticleChange::Initialize(), CLHEP::detail::n, G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeTrackStatus(), G4VParticleChange::ProposeWeight(), G4ElementSelector::SelectZandA(), G4HadProjectile::SetBoundEnergy(), G4VParticleChange::SetNumberOfSecondaries(), G4Track::SetTouchableHandle(), G4Track::SetWeight(), stopAndKill, G4HadronicProcess::thePro, and G4HadronicProcess::theTotalResult.

{
  // if primary is not Alive then do nothing
  theTotalResult->Initialize(track);

  G4Nucleus* nucleus = GetTargetNucleusPointer();
  G4Element* elm = fElementSelector->SelectZandA(track, nucleus);

  G4HadFinalState* result = 0;
  thePro.Initialise(track);
  G4double time0 = track.GetGlobalTime();
  G4bool nuclearCapture = true;

  // Do the electromagnetic cascade in the nuclear field.
  // EM cascade should keep G4HadFinalState object,
  // because it will not be deleted at the end of this method
  //
  result = fEmCascade->ApplyYourself(thePro, *nucleus);
  G4double ebound = result->GetLocalEnergyDeposit(); 
  G4double edep = 0.0; 
  G4int nSecondaries = result->GetNumberOfSecondaries();

  // Try decay from bound level 
  // For mu- the time of projectile should be changed.
  // Decay should keep G4HadFinalState object,
  // because it will not be deleted at the end of this method
  //
  thePro.SetBoundEnergy(ebound);
  if(fBoundDecay) {
    G4HadFinalState* resultDecay = 
      fBoundDecay->ApplyYourself(thePro, *nucleus);
    G4int n = resultDecay->GetNumberOfSecondaries();
    if(0 < n) {
      nSecondaries += n;
      result->AddSecondaries(resultDecay);
    } 
    if(resultDecay->GetStatusChange() == stopAndKill) {
      nuclearCapture = false; 
    }
    resultDecay->Clear();
  }

  if(nuclearCapture) {

    // select model
    G4HadronicInteraction* model = 0;
    try {
      model = ChooseHadronicInteraction(0.0, track.GetMaterial(), elm);
    }
    catch(G4HadronicException & aE) {
      G4ExceptionDescription ed;
      ed << "Target element "<<elm->GetName()<<"  Z= " 
         << nucleus->GetZ_asInt() << "  A= " 
         << nucleus->GetA_asInt() << G4endl;
      DumpState(track,"ChooseHadronicInteraction",ed);
      ed << " No HadronicInteraction found out" << G4endl;
      G4Exception("G4HadronStoppingProcess::AtRestDoIt", "had005", 
                  FatalException, ed);
    }

    G4HadFinalState* resultNuc = 0;
    G4int reentryCount = 0;
    do {
      // sample final state
      // nuclear interaction should keep G4HadFinalState object
      // model should define time of each secondary particle
      try {
        resultNuc = model->ApplyYourself(thePro, *nucleus);
        ++reentryCount;
      }
      catch(G4HadronicException aR) {
        G4ExceptionDescription ed;
        ed << "Call for " << model->GetModelName() << G4endl;
        ed << "Target element "<<elm->GetName()<<"  Z= " 
           << nucleus->GetZ_asInt() 
           << "  A= " << nucleus->GetA_asInt() << G4endl;
        DumpState(track,"ApplyYourself",ed);
        ed << " ApplyYourself failed" << G4endl;
        G4Exception("G4HadronStoppingProcess::AtRestDoIt", "had006", 
                    FatalException, ed);
      }

      // Check the result for catastrophic energy non-conservation
      resultNuc = CheckResult(thePro, *nucleus, resultNuc);

      if(reentryCount>100) {
        G4ExceptionDescription ed;
        ed << "Call for " << model->GetModelName() << G4endl;
        ed << "Target element "<<elm->GetName()<<"  Z= " 
           << nucleus->GetZ_asInt() 
           << "  A= " << nucleus->GetA_asInt() << G4endl;
        DumpState(track,"ApplyYourself",ed);
        ed << " ApplyYourself does not completed after 100 attempts" << G4endl;
        G4Exception("G4HadronStoppingProcess::AtRestDoIt", "had006", 
                    FatalException, ed);  
      }
    }
    while(!resultNuc);

    edep = resultNuc->GetLocalEnergyDeposit();
    nSecondaries += resultNuc->GetNumberOfSecondaries();
    result->AddSecondaries(resultNuc); 
    resultNuc->Clear();
  }

  // Fill results
  //
  theTotalResult->ProposeTrackStatus(fStopAndKill);
  theTotalResult->ProposeLocalEnergyDeposit(edep);  
  theTotalResult->SetNumberOfSecondaries(nSecondaries);
  G4double w  = track.GetWeight();
  theTotalResult->ProposeWeight(w);
  for(G4int i=0; i<nSecondaries; ++i) {
    G4HadSecondary* sec = result->GetSecondary(i);

    // add track global time to the reaction time
    G4double time = sec->GetTime();
    if(time < 0.0) { time = 0.0; }
    time += time0;

    // create secondary track
    G4Track* t = new G4Track(sec->GetParticle(),
                             time, 
                             track.GetPosition());
    t->SetWeight(w*sec->GetWeight());
    t->SetTouchableHandle(track.GetTouchableHandle());
    theTotalResult->AddSecondary(t);
  }
  result->Clear();

  if (epReportLevel != 0) { 
    CheckEnergyMomentumConservation(track, *nucleus);
  }
  return theTotalResult;
}

G4double G4HadronStoppingProcess::AtRestGetPhysicalInteractionLength	(	const G4Track &	track,
		G4ForceCondition *	condition
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 108 of file G4HadronStoppingProcess.cc.

References NotForced.

{
  *condition = NotForced;
  return 0.0;
}

void G4HadronStoppingProcess::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4HadronicProcess.

Definition at line 101 of file G4HadronStoppingProcess.cc.

References G4HadronicProcessStore::Instance(), and G4HadronicProcessStore::PrintInfo().

{
  G4HadronicProcessStore::Instance()->PrintInfo(&p);
}

G4double G4HadronStoppingProcess::GetMeanLifeTime	(	const G4Track &	,
		G4ForceCondition *
	)			[inline, protected]

Definition at line 98 of file G4HadronStoppingProcess.hh.

                                                            { return -1.0; }

G4bool G4HadronStoppingProcess::IsApplicable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Reimplemented in G4HadronicAbsorptionBertini, G4HadronicAbsorptionFritiof, and G4MuonMinusCapture.

Definition at line 87 of file G4HadronStoppingProcess.cc.

References G4ParticleDefinition::GetPDGCharge().

Referenced by G4HadronicAbsorptionBertini::IsApplicable().

{
  return (p.GetPDGCharge() < 0.0);
}

G4double G4HadronStoppingProcess::PostStepGetPhysicalInteractionLength	(	const G4Track &	track,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 117 of file G4HadronStoppingProcess.cc.

References DBL_MAX, and NotForced.

{
  *condition = NotForced;
  return DBL_MAX;
}

void G4HadronStoppingProcess::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4HadronicProcess.

Definition at line 94 of file G4HadronStoppingProcess.cc.

References G4HadronicProcessStore::Instance(), and G4HadronicProcessStore::RegisterParticleForExtraProcess().

{
  G4HadronicProcessStore::Instance()->RegisterParticleForExtraProcess(this, &p);
}

void G4HadronStoppingProcess::ProcessDescription

(

std::ostream &

outFile

)

const [virtual]

Reimplemented from G4HadronicProcess.

Reimplemented in G4HadronicAbsorptionBertini, G4HadronicAbsorptionFritiof, and G4MuonMinusCapture.

Definition at line 265 of file G4HadronStoppingProcess.cc.

{
  outFile << "Base process for negatively charged particle capture at rest.\n";
}

void G4HadronStoppingProcess::SetBoundDecay

(

G4HadronicInteraction *

ptr

)

[inline]

Definition at line 133 of file G4HadronStoppingProcess.hh.

Referenced by G4MuonMinusCapture::G4MuonMinusCapture().

{
  fBoundDecay = ptr;
}

void G4HadronStoppingProcess::SetElementSelector

(

G4ElementSelector *

ptr

)

[inline]

Definition at line 118 of file G4HadronStoppingProcess.hh.

{
  if(fElementSelector != ptr) {
    delete fElementSelector;
    fElementSelector = ptr;
  }
}

void G4HadronStoppingProcess::SetEmCascade

(

G4HadronicInteraction *

ptr

)

[inline]

Definition at line 127 of file G4HadronStoppingProcess.hh.

{
  fEmCascade = ptr;
}

---

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

• G4HadronStoppingProcess.hh
• G4HadronStoppingProcess.cc

---