G4NeutronCaptureAtRest Class Reference

#include <G4NeutronCaptureAtRest.hh>

Inheritance diagram for G4NeutronCaptureAtRest:

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

---

Constructor & Destructor Documentation

G4NeutronCaptureAtRest::G4NeutronCaptureAtRest	(	const G4String &	processName = "NeutronCaptureAtRest",
		G4ProcessType	aType = fHadronic
	)

Definition at line 46 of file G4NeutronCaptureAtRest.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   massProton(G4Proton::Proton()->GetPDGMass()/GeV),
00050   massNeutron(G4Neutron::Neutron()->GetPDGMass()/GeV),
00051   massElectron(G4Electron::Electron()->GetPDGMass()/GeV),
00052   massDeuteron(G4Deuteron::Deuteron()->GetPDGMass()/GeV),
00053   massAlpha(G4Alpha::Alpha()->GetPDGMass()/GeV),
00054   pdefGamma(G4Gamma::Gamma()),
00055   pdefNeutron(G4Neutron::Neutron())
00056 {
00057   G4HadronicDeprecate("G4NeutronCaptureAtRest");
00058   if (verboseLevel>0) {
00059     G4cout << GetProcessName() << " is created "<< G4endl;
00060   }
00061   SetProcessSubType(fHadronAtRest);
00062   pv   = new G4GHEKinematicsVector [MAX_SECONDARIES+1];
00063   eve  = new G4GHEKinematicsVector [MAX_SECONDARIES];
00064   gkin = new G4GHEKinematicsVector [MAX_SECONDARIES];
00065 
00066   G4HadronicProcessStore::Instance()->RegisterExtraProcess(this);
00067 }
 

G4NeutronCaptureAtRest::~G4NeutronCaptureAtRest

(

)

Definition at line 71 of file G4NeutronCaptureAtRest.cc.

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

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

---

Member Function Documentation

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

Reimplemented from G4VRestProcess.

Definition at line 139 of file G4NeutronCaptureAtRest.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 << "G4NeutronCaptureAtRest::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 Neutron

//   clear InteractionLengthLeft

  ResetNumberOfInteractionLengthLeft();

  return &aParticleChange;

}

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

Reimplemented from G4VRestProcess.

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

  return theNumberOfInteractionLengthLeft * currentInteractionLength;

}

void G4NeutronCaptureAtRest::BuildPhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 84 of file G4NeutronCaptureAtRest.cc.

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

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

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

Implements G4VRestProcess.

Definition at line 73 of file G4NeutronCaptureAtRest.hh.

Referenced by AtRestGetPhysicalInteractionLength().

                                                  {return 0.0;}

G4int G4NeutronCaptureAtRest::GetNumberOfSecondaries

(

)

Definition at line 100 of file G4NeutronCaptureAtRest.cc.

{
  return ( ngkine );

}

G4GHEKinematicsVector * G4NeutronCaptureAtRest::GetSecondaryKinematics

(

)

Definition at line 107 of file G4NeutronCaptureAtRest.cc.

{
  return ( &gkin[0] );

}

G4bool G4NeutronCaptureAtRest::IsApplicable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 91 of file G4NeutronCaptureAtRest.cc.

{
   return ( &particle == pdefNeutron );

}

void G4NeutronCaptureAtRest::PreparePhysicsTable

(

const G4ParticleDefinition &

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 79 of file G4NeutronCaptureAtRest.cc.

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

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

---

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

• G4NeutronCaptureAtRest.hh
• G4NeutronCaptureAtRest.cc

---