G4KaonMinusAbsorption Class Reference

#include <G4KaonMinusAbsorption.hh>

Inheritance diagram for G4KaonMinusAbsorption:

Public Member Functions
	G4KaonMinusAbsorption (const G4String &processName="KaonMinusAbsorption", G4ProcessType aType=fHadronic)
	~G4KaonMinusAbsorption ()
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 47 of file G4KaonMinusAbsorption.hh.

---

Constructor & Destructor Documentation

G4KaonMinusAbsorption::G4KaonMinusAbsorption	(	const G4String &	processName = "KaonMinusAbsorption",
		G4ProcessType	aType = fHadronic
	)

Definition at line 46 of file G4KaonMinusAbsorption.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   massKaonMinus(G4KaonMinus::KaonMinus()->GetPDGMass()/GeV),
00050   massGamma(G4Gamma::Gamma()->GetPDGMass()/GeV),
00051   massPionZero(G4PionZero::PionZero()->GetPDGMass()/GeV),
00052   massProton(G4Proton::Proton()->GetPDGMass()/GeV),
00053   massLambda(G4Lambda::Lambda()->GetPDGMass()/GeV),
00054   pdefKaonMinus(G4KaonMinus::KaonMinus()),
00055   pdefGamma(G4Gamma::Gamma()),
00056   pdefPionZero(G4PionZero::PionZero()),
00057   pdefProton(G4Proton::Proton()),
00058   pdefNeutron(G4Neutron::Neutron()),
00059   pdefLambda(G4Lambda::Lambda()),
00060   pdefDeuteron(G4Deuteron::Deuteron()),
00061   pdefTriton(G4Triton::Triton()),
00062   pdefAlpha(G4Alpha::Alpha())
00063 {
00064   G4HadronicDeprecate("G4KaonMinusAbsorption");
00065   if (verboseLevel>0) {
00066     G4cout << GetProcessName() << " is created "<< G4endl;
00067   }
00068   SetProcessSubType(fHadronAtRest);
00069   pv   = new G4GHEKinematicsVector [MAX_SECONDARIES+1];
00070   eve  = new G4GHEKinematicsVector [MAX_SECONDARIES];
00071   gkin = new G4GHEKinematicsVector [MAX_SECONDARIES];
00072 
00073   G4HadronicProcessStore::Instance()->RegisterExtraProcess(this);
00074 }
 

G4KaonMinusAbsorption::~G4KaonMinusAbsorption

(

)

Definition at line 78 of file G4KaonMinusAbsorption.cc.

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

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

---

Member Function Documentation

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

Reimplemented from G4VRestProcess.

Definition at line 146 of file G4KaonMinusAbsorption.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 << "G4KaonMinusAbsorption::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 KaonMinus

//   clear InteractionLengthLeft

  ResetNumberOfInteractionLengthLeft();

  return &aParticleChange;

}

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

Reimplemented from G4VRestProcess.

Definition at line 120 of file G4KaonMinusAbsorption.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 << "G4KaonMinusAbsorptionProcess::AtRestGetPhysicalInteractionLength ";
    G4cout << "[ " << GetProcessName() << "]" <<G4endl;
    track.GetDynamicParticle()->DumpInfo();
    G4cout << " in Material  " << track.GetMaterial()->GetName() <<G4endl;
    G4cout << "MeanLifeTime = " << currentInteractionLength/ns << "[ns]" <<G4endl;
  }

  return theNumberOfInteractionLengthLeft * currentInteractionLength;

}

void G4KaonMinusAbsorption::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 91 of file G4KaonMinusAbsorption.cc.

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

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

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

Implements G4VRestProcess.

Definition at line 72 of file G4KaonMinusAbsorption.hh.

Referenced by AtRestGetPhysicalInteractionLength().

                                                  {return 0.0;}

G4int G4KaonMinusAbsorption::GetNumberOfSecondaries

(

)

Definition at line 107 of file G4KaonMinusAbsorption.cc.

{
  return ( ngkine );

}

G4GHEKinematicsVector * G4KaonMinusAbsorption::GetSecondaryKinematics

(

)

Definition at line 114 of file G4KaonMinusAbsorption.cc.

{
  return ( &gkin[0] );

}

G4bool G4KaonMinusAbsorption::IsApplicable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 98 of file G4KaonMinusAbsorption.cc.

{
   return ( &particle == pdefKaonMinus );

}

void G4KaonMinusAbsorption::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 86 of file G4KaonMinusAbsorption.cc.

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

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

---

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

• G4KaonMinusAbsorption.hh
• G4KaonMinusAbsorption.cc

---