G4AntiProtonAnnihilationAtRest Class Reference

#include <G4AntiProtonAnnihilationAtRest.hh>

Inheritance diagram for G4AntiProtonAnnihilationAtRest:

Public Member Functions
	G4AntiProtonAnnihilationAtRest (const G4String &processName="AntiProtonAnnihilationAtRest", G4ProcessType aType=fHadronic)
	~G4AntiProtonAnnihilationAtRest ()
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 G4AntiProtonAnnihilationAtRest.hh.

---

Constructor & Destructor Documentation

G4AntiProtonAnnihilationAtRest::G4AntiProtonAnnihilationAtRest	(	const G4String &	processName = "AntiProtonAnnihilationAtRest",
		G4ProcessType	aType = fHadronic
	)

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

G4AntiProtonAnnihilationAtRest::~G4AntiProtonAnnihilationAtRest

(

)

Definition at line 80 of file G4AntiProtonAnnihilationAtRest.cc.

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

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

---

Member Function Documentation

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

Reimplemented from G4VRestProcess.

Definition at line 148 of file G4AntiProtonAnnihilationAtRest.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 << "G4AntiProtonAnnihilationAtRest::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 AntiProton

//   clear InteractionLengthLeft

  ResetNumberOfInteractionLengthLeft();

  return &aParticleChange;

}

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

Reimplemented from G4VRestProcess.

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

  return theNumberOfInteractionLengthLeft * currentInteractionLength;

}

void G4AntiProtonAnnihilationAtRest::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 93 of file G4AntiProtonAnnihilationAtRest.cc.

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

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

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

Implements G4VRestProcess.

Definition at line 71 of file G4AntiProtonAnnihilationAtRest.hh.

Referenced by AtRestGetPhysicalInteractionLength().

                                                  {return 0.0;}

G4int G4AntiProtonAnnihilationAtRest::GetNumberOfSecondaries

(

)

Definition at line 109 of file G4AntiProtonAnnihilationAtRest.cc.

{
  return ( ngkine );

}

G4GHEKinematicsVector * G4AntiProtonAnnihilationAtRest::GetSecondaryKinematics

(

)

Definition at line 116 of file G4AntiProtonAnnihilationAtRest.cc.

{
  return ( &gkin[0] );

}

G4bool G4AntiProtonAnnihilationAtRest::IsApplicable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 100 of file G4AntiProtonAnnihilationAtRest.cc.

{
   return ( &particle == pdefAntiProton );

}

void G4AntiProtonAnnihilationAtRest::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 88 of file G4AntiProtonAnnihilationAtRest.cc.

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

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

---

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

• G4AntiProtonAnnihilationAtRest.hh
• G4AntiProtonAnnihilationAtRest.cc

---