G4SmartTrackStack Class Reference

#include <G4SmartTrackStack.hh>

Public Member Functions
	G4SmartTrackStack ()
	~G4SmartTrackStack ()
void	PushToStack (const G4StackedTrack &aStackedTrack)
G4StackedTrack	PopFromStack ()
void	clear ()
void	clearAndDestroy ()
void	TransferTo (G4TrackStack *aStack)
G4double	getEnergyOfStack (G4TrackStack *aTrackStack)
void	dumpStatistics ()
G4int	GetNTrack () const
G4int	GetMaxNTrack () const

Detailed Description

Definition at line 46 of file G4SmartTrackStack.hh.

Constructor & Destructor Documentation

G4SmartTrackStack::G4SmartTrackStack

(

)

Definition at line 45 of file G4SmartTrackStack.cc.

    :fTurn(0), nTurn(5), maxNTracks(0), nTracks(0)
{
    for(int i=0;i<nTurn;i++)
    {
      stacks[i] = new G4TrackStack(5000);
      energies[i] = 0.;
    }
}

G4SmartTrackStack::~G4SmartTrackStack

(

)

Definition at line 55 of file G4SmartTrackStack.cc.

{
    for (int i  = 0; i < nTurn; i++) {
        delete stacks[i];
    }
}

Member Function Documentation

void G4SmartTrackStack::clear

(

void

)

Definition at line 142 of file G4SmartTrackStack.cc.

{
  for (int i = 0; i < nTurn; i++) {
    stacks[i]->clear();
    energies[i] = 0.0;
    fTurn = 0;
  }
  nTracks = 0;
}

void G4SmartTrackStack::clearAndDestroy

(

)

Definition at line 152 of file G4SmartTrackStack.cc.

References G4TrackStack::clearAndDestroy().

{
  for (int i = 0; i < nTurn; i++) {
    stacks[i]->clearAndDestroy();
    energies[i] = 0.0;
    fTurn = 0;
  }
  nTracks = 0;
}

void G4SmartTrackStack::dumpStatistics

(

)

Definition at line 34 of file G4SmartTrackStack.cc.

References G4cerr, G4endl, G4TrackStack::GetNTrack(), and G4TrackStack::getTotalEnergy().

{
    // Print to stderr so that we can split stats output from normal
    // output of Geant4 which is typically being printed to stdout
    for (int i = 0; i < nTurn; i++) {
        G4cerr << stacks[i]->GetNTrack() << " ";
        G4cerr << stacks[i]->getTotalEnergy() << " ";
    }
    G4cerr << G4endl;
}

G4double G4SmartTrackStack::getEnergyOfStack

(

G4TrackStack *

aTrackStack

)

G4int G4SmartTrackStack::GetMaxNTrack ( ) const

inline

Definition at line 82 of file G4SmartTrackStack.hh.

{ return maxNTracks; }

G4int G4SmartTrackStack::GetNTrack ( ) const

inline

Definition at line 81 of file G4SmartTrackStack.hh.

Referenced by PopFromStack(), and PushToStack().

{ return nTracks; }

G4StackedTrack G4SmartTrackStack::PopFromStack

(

)

Definition at line 83 of file G4SmartTrackStack.cc.

References G4Track::GetDynamicParticle(), GetNTrack(), G4DynamicParticle::GetTotalEnergy(), G4StackedTrack::GetTrack(), and G4TrackStack::PopFromStack().

{
    G4StackedTrack aStackedTrack;

    if (nTracks) {
        while (true) {
            if (stacks[fTurn]->GetNTrack()) {
                aStackedTrack = stacks[fTurn]->PopFromStack();
                energies[fTurn] -= aStackedTrack.GetTrack()->GetDynamicParticle()->GetTotalEnergy();
        nTracks--;
        break;
            } else {
                fTurn = (fTurn+1) % nTurn;
            }
        }
    }

    // dumpStatistics();
    return aStackedTrack;
}

void G4SmartTrackStack::PushToStack

(

const G4StackedTrack &

aStackedTrack

)

Definition at line 108 of file G4SmartTrackStack.cc.

References electronCode, gammaCode, G4Track::GetDynamicParticle(), G4TrackStack::GetNTrack(), GetNTrack(), G4Track::GetParentID(), G4DynamicParticle::GetPDGcode(), G4TrackStack::GetSafetyValve1(), G4TrackStack::GetSafetyValve2(), G4DynamicParticle::GetTotalEnergy(), G4StackedTrack::GetTrack(), neutronCode, positronCode, and G4TrackStack::PushToStack().

Referenced by G4TrackStack::TransferTo().

{

  G4int iDest = 0;
  if (aStackedTrack.GetTrack()->GetParentID()) {
    G4int code = aStackedTrack.GetTrack()->GetDynamicParticle()->GetPDGcode();
    if (code == electronCode)
      iDest = 2;
    else if (code == gammaCode)
      iDest = 3;
    else if (code == positronCode)
      iDest = 4;
    else if (code == neutronCode)
      iDest = 1;
  } else {
    // We have a primary track, which should go first.
    fTurn = 0; // reseting the turn
  }
  stacks[iDest]->PushToStack(aStackedTrack);
  energies[iDest] += aStackedTrack.GetTrack()->GetDynamicParticle()->GetTotalEnergy();
  nTracks++;
  
  G4int dy1 = stacks[iDest]->GetNTrack() - stacks[iDest]->GetSafetyValve1();
  G4int dy2 = stacks[fTurn]->GetNTrack() - stacks[fTurn]->GetSafetyValve2();
  
  if (dy1 > 0 || dy1 > dy2 ||
      (iDest == 2 &&
       stacks[iDest]->GetNTrack() < 50 && energies[iDest] < energies[fTurn])) {
        fTurn = iDest;
      }
  
  if (nTracks > maxNTracks) maxNTracks = nTracks;
}

void G4SmartTrackStack::TransferTo

(

G4TrackStack *

aStack

)

Definition at line 75 of file G4SmartTrackStack.cc.

References G4TrackStack::TransferTo().

{
    for (int i = 0; i < nTurn; i++) {
        stacks[i]->TransferTo(aStack);
    }
  nTracks = 0;
}

---

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

• G4SmartTrackStack.hh
• G4SmartTrackStack.cc