G4PionMinusAbsorptionAtRest Class Reference

#include <G4PionMinusAbsorptionAtRest.hh>

Inheritance diagram for G4PionMinusAbsorptionAtRest:

Public Member Functions
	G4PionMinusAbsorptionAtRest (const G4String &processName="PionMinusAbsorptionAtRest", G4ProcessType aType=fHadronic)
	~G4PionMinusAbsorptionAtRest ()
G4bool	IsApplicable (const G4ParticleDefinition &)
void	PreparePhysicsTable (const G4ParticleDefinition &)
void	BuildPhysicsTable (const G4ParticleDefinition &)
G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
G4double	GetMeanLifeTime (const G4Track &, G4ForceCondition *)
G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
G4int	GetNumberOfSecondaries ()
G4GHEKinematicsVector *	GetSecondaryKinematics ()

---

Detailed Description

Definition at line 46 of file G4PionMinusAbsorptionAtRest.hh.

---

Constructor & Destructor Documentation

G4PionMinusAbsorptionAtRest::G4PionMinusAbsorptionAtRest	(	const G4String &	processName = "PionMinusAbsorptionAtRest",
		G4ProcessType	aType = fHadronic
	)

Definition at line 46 of file G4PionMinusAbsorptionAtRest.cc.

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

00047                                                               :
00048   G4VRestProcess (processName, aType),       // initialization
00049   massPionMinus(G4PionMinus::PionMinus()->GetPDGMass()/GeV),
00050   pdefGamma(G4Gamma::Gamma()),
00051   pdefPionZero(G4PionZero::PionZero()),
00052   pdefPionMinus(G4PionMinus::PionMinus()),
00053   pdefProton(G4Proton::Proton()),
00054   pdefNeutron(G4Neutron::Neutron()),
00055   pdefDeuteron(G4Deuteron::Deuteron()),
00056   pdefTriton(G4Triton::Triton()),
00057   pdefAlpha(G4Alpha::Alpha())
00058 {
00059   G4HadronicDeprecate("G4PiMinusAbsorptionAtRest");
00060 
00061   if (verboseLevel>0) {
00062     G4cout << GetProcessName() << " is created "<< G4endl;
00063   }
00064   SetProcessSubType(fHadronAtRest);
00065   pv   = new G4GHEKinematicsVector [MAX_SECONDARIES+1];
00066   eve  = new G4GHEKinematicsVector [MAX_SECONDARIES];
00067   gkin = new G4GHEKinematicsVector [MAX_SECONDARIES];
00068 
00069   G4HadronicProcessStore::Instance()->RegisterExtraProcess(this);
00070 }
 

G4PionMinusAbsorptionAtRest::~G4PionMinusAbsorptionAtRest

(

)

Definition at line 74 of file G4PionMinusAbsorptionAtRest.cc.

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

{
  G4HadronicProcessStore::Instance()->DeRegisterExtraProcess(this);
  delete [] pv;
  delete [] eve;
  delete [] gkin;
}

---

Member Function Documentation

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

Reimplemented from G4VRestProcess.

Definition at line 142 of file G4PionMinusAbsorptionAtRest.cc.

References G4ParticleChange::AddSecondary(), G4VProcess::aParticleChange, fStopAndKill, G4cout, G4endl, G4UniformRand, G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), G4Track::GetGlobalTime(), G4Track::GetMaterial(), G4Material::GetNumberOfElements(), G4Track::GetPosition(), G4GHEKinematicsVector::GetTOF(), G4Track::GetTouchableHandle(), G4ParticleChange::Initialize(), position, G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeTrackStatus(), G4VProcess::ResetNumberOfInteractionLengthLeft(), G4VParticleChange::SetNumberOfSecondaries(), G4Track::SetTouchableHandle(), and G4VProcess::verboseLevel.

{

//   Initialize ParticleChange
//     all members of G4VParticleChange are set to equal to 
//     corresponding member in G4Track

  aParticleChange.Initialize(track);

//   Store some global quantities that depend on current material and particle

  globalTime = track.GetGlobalTime()/s;
  G4Material * aMaterial = track.GetMaterial();
  const G4int numberOfElements = aMaterial->GetNumberOfElements();
  const G4ElementVector* theElementVector = aMaterial->GetElementVector();

  const G4double* theAtomicNumberDensity = aMaterial->GetAtomicNumDensityVector();
  G4double normalization = 0;
  for ( G4int i1=0; i1 < numberOfElements; i1++ )
  {
    normalization += theAtomicNumberDensity[i1] ; // change when nucleon specific
                                                  // probabilities are included.
  }
  G4double runningSum= 0.;
  G4double random = G4UniformRand()*normalization;
  for ( G4int i2=0; i2 < numberOfElements; i2++ )
  {
    runningSum += theAtomicNumberDensity[i2]; // change when nucleon specific
                                              // probabilities are included.
    if (random<=runningSum)
    {
      targetCharge = G4double((*theElementVector)[i2]->GetZ());
      targetAtomicMass = (*theElementVector)[i2]->GetN();
    }
  }
  if (random>runningSum)
  {
    targetCharge = G4double((*theElementVector)[numberOfElements-1]->GetZ());
    targetAtomicMass = (*theElementVector)[numberOfElements-1]->GetN();

  }

  if (verboseLevel>1) {
    G4cout << "G4PionMinusAbsorptionAtRest::AtRestDoIt is invoked " <<G4endl;
    }

  G4ParticleMomentum momentum;
  G4float localtime;

  G4ThreeVector   position = track.GetPosition();

  GenerateSecondaries(); // Generate secondaries

  aParticleChange.SetNumberOfSecondaries( ngkine ); 

  for ( G4int isec = 0; isec < ngkine; isec++ ) {
    G4DynamicParticle* aNewParticle = new G4DynamicParticle;
    aNewParticle->SetDefinition( gkin[isec].GetParticleDef() );
    aNewParticle->SetMomentum( gkin[isec].GetMomentum() * GeV );

    localtime = globalTime + gkin[isec].GetTOF();

    G4Track* aNewTrack = new G4Track( aNewParticle, localtime*s, position );
                aNewTrack->SetTouchableHandle(track.GetTouchableHandle());
    aParticleChange.AddSecondary( aNewTrack );

  }

  aParticleChange.ProposeLocalEnergyDeposit( 0.0*GeV );

  aParticleChange.ProposeTrackStatus(fStopAndKill); // Kill the incident PionMinus

//   clear InteractionLengthLeft

  ResetNumberOfInteractionLengthLeft();

  return &aParticleChange;

}

G4double G4PionMinusAbsorptionAtRest::AtRestGetPhysicalInteractionLength	(	const G4Track &	,
		G4ForceCondition *
	)			[virtual]

Reimplemented from G4VRestProcess.

Definition at line 116 of file G4PionMinusAbsorptionAtRest.cc.

References G4VProcess::currentInteractionLength, G4DynamicParticle::DumpInfo(), G4cout, G4endl, G4Track::GetDynamicParticle(), G4Track::GetMaterial(), GetMeanLifeTime(), G4Material::GetName(), G4VProcess::GetProcessName(), NotForced, ns, G4VProcess::ResetNumberOfInteractionLengthLeft(), G4VProcess::theNumberOfInteractionLengthLeft, and G4VProcess::verboseLevel.

{
  // beggining of tracking 
  ResetNumberOfInteractionLengthLeft();

  // condition is set to "Not Forced"
  *condition = NotForced;

  // get mean life time
  currentInteractionLength = GetMeanLifeTime(track, condition);

  if ((currentInteractionLength <0.0) || (verboseLevel>2)){
    G4cout << "G4PionMinusAbsorptionAtRestProcess::AtRestGetPhysicalInteractionLength ";
    G4cout << "[ " << GetProcessName() << "]" <<G4endl;
    track.GetDynamicParticle()->DumpInfo();
    G4cout << " in Material  " << track.GetMaterial()->GetName() <<G4endl;
    G4cout << "MeanLifeTime = " << currentInteractionLength/ns << "[ns]" <<G4endl;
  }

  return theNumberOfInteractionLengthLeft * currentInteractionLength;

}

void G4PionMinusAbsorptionAtRest::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 87 of file G4PionMinusAbsorptionAtRest.cc.

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

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

G4double G4PionMinusAbsorptionAtRest::GetMeanLifeTime	(	const G4Track &	,
		G4ForceCondition *
	)			[inline, virtual]

Implements G4VRestProcess.

Definition at line 72 of file G4PionMinusAbsorptionAtRest.hh.

Referenced by AtRestGetPhysicalInteractionLength().

                                                  {return 0.0;}

G4int G4PionMinusAbsorptionAtRest::GetNumberOfSecondaries

(

)

Definition at line 103 of file G4PionMinusAbsorptionAtRest.cc.

{
  return ( ngkine );

}

G4GHEKinematicsVector * G4PionMinusAbsorptionAtRest::GetSecondaryKinematics

(

)

Definition at line 110 of file G4PionMinusAbsorptionAtRest.cc.

{
  return ( &gkin[0] );

}

G4bool G4PionMinusAbsorptionAtRest::IsApplicable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 94 of file G4PionMinusAbsorptionAtRest.cc.

{
   return ( &particle == pdefPionMinus );

}

void G4PionMinusAbsorptionAtRest::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 82 of file G4PionMinusAbsorptionAtRest.cc.

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

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

---

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

• G4PionMinusAbsorptionAtRest.hh
• G4PionMinusAbsorptionAtRest.cc

---