G4PiMinusAbsorptionAtRest Class Reference

#include <G4PiMinusAbsorptionAtRest.hh>

Inheritance diagram for G4PiMinusAbsorptionAtRest:

Public Member Functions
	G4PiMinusAbsorptionAtRest (const G4String &processName="PiMinusAbsorptionAtRest", G4ProcessType aType=fHadronic)
	~G4PiMinusAbsorptionAtRest ()
G4bool	IsApplicable (const G4ParticleDefinition &particle)
void	PreparePhysicsTable (const G4ParticleDefinition &)
void	BuildPhysicsTable (const G4ParticleDefinition &)
G4VParticleChange *	AtRestDoIt (const G4Track &aTrack, const G4Step &aStep)
void	SetDeexcitationAlgorithm (G4int index)
Protected Member Functions
G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *)

---

Detailed Description

Definition at line 54 of file G4PiMinusAbsorptionAtRest.hh.

---

Constructor & Destructor Documentation

G4PiMinusAbsorptionAtRest::G4PiMinusAbsorptionAtRest	(	const G4String &	processName = "PiMinusAbsorptionAtRest",
		G4ProcessType	aType = fHadronic
	)

Definition at line 65 of file G4PiMinusAbsorptionAtRest.cc.

References fHadronAtRest, G4cout, G4endl, G4HadronicDeprecate, G4VProcess::GetProcessName(), G4HadronicProcessStore::Instance(), G4HadronicProcessStore::RegisterExtraProcess(), G4VProcess::SetProcessSubType(), and G4VProcess::verboseLevel.

                                                                          :
  G4VRestProcess (processName, aType)
{
  G4HadronicDeprecate("G4PiMinusAbsorptionAtRest");

  SetProcessSubType(fHadronAtRest);

  _indexDeexcitation = 0;

  if (verboseLevel>0) 
    { G4cout << GetProcessName() << " is created "<< G4endl; }
  G4HadronicProcessStore::Instance()->RegisterExtraProcess(this);
}

G4PiMinusAbsorptionAtRest::~G4PiMinusAbsorptionAtRest

(

)

Definition at line 83 of file G4PiMinusAbsorptionAtRest.cc.

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

{
  G4HadronicProcessStore::Instance()->DeRegisterExtraProcess(this);
}

---

Member Function Documentation

G4VParticleChange * G4PiMinusAbsorptionAtRest::AtRestDoIt	(	const G4Track &	aTrack,
		const G4Step &	aStep
	)			[virtual]

Reimplemented from G4VRestProcess.

Definition at line 98 of file G4PiMinusAbsorptionAtRest.cc.

References G4ParticleChange::AddSecondary(), G4VProcess::aParticleChange, G4PiMinusStopAbsorption::DoAbsorption(), G4StopDeexcitation::DoBreakUp(), G4PiMinusStopAbsorption::Energy(), FatalException, fStopAndKill, G4cerr, G4cout, G4endl, G4Exception(), G4UniformRand, G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4Material::GetElementVector(), G4Material::GetFractionVector(), G4ParticleChange::GetGlobalTime(), G4Track::GetMaterial(), G4Material::GetNumberOfElements(), G4DynamicParticle::GetTotalEnergy(), G4ParticleChange::Initialize(), IsApplicable(), G4PiMinusStopAbsorption::NNeutrons(), G4PiMinusStopAbsorption::NProtons(), G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeTrackStatus(), G4PiMinusStopAbsorption::RecoilMomentum(), G4VParticleChange::SetNumberOfSecondaries(), G4PiMinusStopAbsorption::SetVerboseLevel(), and G4VProcess::verboseLevel.

{
  const G4DynamicParticle* stoppedHadron = track.GetDynamicParticle();
      
  // Check applicability
  if (! IsApplicable(*(stoppedHadron->GetDefinition())))
    {
      G4cerr  
      << "G4PiMinusAbsorptionAtRest: ERROR, particle must be a pion minus!" 
      << G4endl;
      return NULL;
    }

  // Get the current material
  const G4Material* material = track.GetMaterial();

  G4double A=-1;
  G4double Z=-1;
  G4double random = G4UniformRand();
  const G4ElementVector* theElementVector = material->GetElementVector();
  unsigned int i;
  G4double sum = 0;
  G4double totalsum=0;
  for(i=0; i<material->GetNumberOfElements(); ++i)
    {
      if((*theElementVector)[i]->GetZ()!=1) totalsum+=material->GetFractionVector()[i];
    }
  for (i = 0; i<material->GetNumberOfElements(); ++i)
    {
      if((*theElementVector)[i]->GetZ()!=1) sum += material->GetFractionVector()[i];
      if ( sum/totalsum > random )  
        { 
          A = (*theElementVector)[i]->GetA()*mole/g;
          Z = (*theElementVector)[i]->GetZ();
          break;
        }
    }

  // Do the interaction with the nucleon cluster
  
  G4PiMinusStopMaterial* algorithm = LoadAlgorithm(static_cast<G4int>(Z));
  G4PiMinusStopAbsorption stopAbsorption(algorithm,Z,A);
  stopAbsorption.SetVerboseLevel(verboseLevel);

  G4DynamicParticleVector* absorptionProducts = stopAbsorption.DoAbsorption();

  // Deal with the leftover nucleus

  G4double pionEnergy = stoppedHadron->GetTotalEnergy();
  G4double excitation = pionEnergy - stopAbsorption.Energy();
  if (excitation < 0.) 
  {
    G4Exception("G4PiMinusAbsorptionAtRest::AtRestDoIt()", "HAD_STOP_0000",
                FatalException, "Excitation energy < 0");
  }
  if (verboseLevel>0) { G4cout << " excitation " << excitation << G4endl; }

  G4StopDeexcitationAlgorithm* nucleusAlgorithm = LoadNucleusAlgorithm();
  G4StopDeexcitation stopDeexcitation(nucleusAlgorithm);

  G4double newZ = Z - stopAbsorption.NProtons();
  G4double newN = A - Z - stopAbsorption.NNeutrons();
  G4double newA = newZ + newN;
  G4ReactionProductVector* fragmentationProducts = stopDeexcitation.DoBreakUp(newA,newZ,excitation,stopAbsorption.RecoilMomentum());

  unsigned int nAbsorptionProducts = 0;
  if (absorptionProducts != 0)     
    { nAbsorptionProducts  =  absorptionProducts->size(); }

  unsigned int nFragmentationProducts = 0;
  if (fragmentationProducts != 0) 
    { nFragmentationProducts = fragmentationProducts->size(); }
  
  if (verboseLevel>0) 
    {
      G4cout << "nAbsorptionProducts = " << nAbsorptionProducts
             << "  nFragmentationProducts = " << nFragmentationProducts
             << G4endl;
    }

  // Deal with ParticleChange final state

  aParticleChange.Initialize(track);
  aParticleChange.SetNumberOfSecondaries(G4int(nAbsorptionProducts + nFragmentationProducts)); 
     
  for (i = 0; i<nAbsorptionProducts; i++)
    { aParticleChange.AddSecondary((*absorptionProducts)[i]); }

//  for (i = 0; i<nFragmentationProducts; i++)
//    { aParticleChange.AddSecondary(fragmentationProducts->at(i)); }
  for(i=0; i<nFragmentationProducts; i++)
  {
    G4DynamicParticle * aNew = 
       new G4DynamicParticle((*fragmentationProducts)[i]->GetDefinition(),
                             (*fragmentationProducts)[i]->GetMomentum());
    G4double newTime = aParticleChange.GetGlobalTime((*fragmentationProducts)[i]->GetFormationTime());
    aParticleChange.AddSecondary(aNew, newTime);
    delete (*fragmentationProducts)[i];
  }

  if (fragmentationProducts != 0)   delete fragmentationProducts;   

  if (_indexDeexcitation == 1) aParticleChange.ProposeLocalEnergyDeposit(excitation);

  // Kill the absorbed pion
  aParticleChange.ProposeTrackStatus(fStopAndKill); 

  return &aParticleChange;

}

void G4PiMinusAbsorptionAtRest::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 93 of file G4PiMinusAbsorptionAtRest.cc.

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

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

G4double G4PiMinusAbsorptionAtRest::GetMeanLifeTime	(	const G4Track &	aTrack,
		G4ForceCondition *
	)			[inline, protected, virtual]

Implements G4VRestProcess.

Definition at line 88 of file G4PiMinusAbsorptionAtRest.hh.

References DBL_MAX, G4Track::GetMaterial(), G4Material::GetNumberOfElements(), and G4Material::GetZ().

  {
     G4double result = 0;
     if(aTrack.GetMaterial()->GetNumberOfElements() == 1)
        if(aTrack.GetMaterial()->GetZ()<1.5) result = DBL_MAX;
     return result;
   }

G4bool G4PiMinusAbsorptionAtRest::IsApplicable

(

const G4ParticleDefinition &

particle

)

[inline, virtual]

Reimplemented from G4VProcess.

Definition at line 74 of file G4PiMinusAbsorptionAtRest.hh.

References G4PionMinus::PionMinus().

Referenced by AtRestDoIt().

    { return ( particle == *(G4PionMinus::PionMinus()) ); }

void G4PiMinusAbsorptionAtRest::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 88 of file G4PiMinusAbsorptionAtRest.cc.

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

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

void G4PiMinusAbsorptionAtRest::SetDeexcitationAlgorithm

(

G4int

index

)

Definition at line 287 of file G4PiMinusAbsorptionAtRest.cc.

{  
  _indexDeexcitation = index;
}

---

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

• G4PiMinusAbsorptionAtRest.hh
• G4PiMinusAbsorptionAtRest.cc

---