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::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::~G4SmartTrackStack

(

)

Definition at line 55 of file G4SmartTrackStack.cc.

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

Member Function Documentation

clear()

void G4SmartTrackStack::clear

(

)

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;
}

clearAndDestroy()

void G4SmartTrackStack::clearAndDestroy

(

)

Definition at line 152 of file G4SmartTrackStack.cc.

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

dumpStatistics()

void G4SmartTrackStack::dumpStatistics

(

)

Definition at line 34 of file G4SmartTrackStack.cc.

{
    // 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;
}

getEnergyOfStack()

G4double G4SmartTrackStack::getEnergyOfStack

(

G4TrackStack *

aTrackStack

)

GetMaxNTrack()

G4int G4SmartTrackStack::GetMaxNTrack ( ) const

inline

Definition at line 82 of file G4SmartTrackStack.hh.

{ return maxNTracks; }

GetNTrack()

G4int G4SmartTrackStack::GetNTrack ( ) const

inline

Definition at line 81 of file G4SmartTrackStack.hh.

{ return nTracks; }

Referenced by PopFromStack(), and PushToStack().

PopFromStack()

G4StackedTrack G4SmartTrackStack::PopFromStack

(

)

Definition at line 83 of file G4SmartTrackStack.cc.

{
    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;
}

PushToStack()

void G4SmartTrackStack::PushToStack

(

const G4StackedTrack &

aStackedTrack

)

Definition at line 108 of file G4SmartTrackStack.cc.

{
 
  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;
}

Referenced by G4TrackStack::TransferTo().

TransferTo()

void G4SmartTrackStack::TransferTo

(

G4TrackStack *

aStack

)

Definition at line 75 of file G4SmartTrackStack.cc.

{
    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