G4SmartTrackStack Class Reference

#include <G4SmartTrackStack.hh>

Public Member Functions
	G4SmartTrackStack ()
	~G4SmartTrackStack ()
G4SmartTrackStack &	operator= (const G4SmartTrackStack &)=delete
G4bool	operator== (const G4SmartTrackStack &) const =delete
G4bool	operator!= (const G4SmartTrackStack &) const =delete
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 44 of file G4SmartTrackStack.hh.

Constructor & Destructor Documentation

G4SmartTrackStack()

G4SmartTrackStack::G4SmartTrackStack

(

)

Definition at line 47 of file G4SmartTrackStack.cc.

{
  for(G4int i=0; i<nTurn; ++i)
  {
    stacks[i] = new G4TrackStack(5000);
    energies[i] = 0.;
  }
}

~G4SmartTrackStack()

G4SmartTrackStack::~G4SmartTrackStack

(

)

Definition at line 56 of file G4SmartTrackStack.cc.

{
  for (auto* sp : stacks)
  {
    delete sp;
  }
}

Member Function Documentation

clear()

void G4SmartTrackStack::clear

(

)

Definition at line 139 of file G4SmartTrackStack.cc.

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

clearAndDestroy()

void G4SmartTrackStack::clearAndDestroy

(

)

Definition at line 150 of file G4SmartTrackStack.cc.

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

dumpStatistics()

void G4SmartTrackStack::dumpStatistics

(

)

Definition at line 35 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 (G4int 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 64 of file G4SmartTrackStack.hh.

{ return maxNTracks; }

GetNTrack()

G4int G4SmartTrackStack::GetNTrack ( ) const

inline

Definition at line 63 of file G4SmartTrackStack.hh.

{ return nTracks; }

Referenced by PopFromStack(), and PushToStack().

operator!=()

G4bool G4SmartTrackStack::operator!= ( const G4SmartTrackStack &  ) const

delete

operator=()

G4SmartTrackStack & G4SmartTrackStack::operator= ( const G4SmartTrackStack &  )

delete

operator==()

G4bool G4SmartTrackStack::operator== ( const G4SmartTrackStack &  ) const

delete

PopFromStack()

G4StackedTrack G4SmartTrackStack::PopFromStack

(

)

Definition at line 73 of file G4SmartTrackStack.cc.

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

PushToStack()

void G4SmartTrackStack::PushToStack

(

const G4StackedTrack &

aStackedTrack

)

Definition at line 101 of file G4SmartTrackStack.cc.

{
 
  G4int iDest = 0;
  if (aStackedTrack.GetTrack()->GetParentID() != 0)
  {
    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;
  
  G4long dy1 = stacks[iDest]->GetNTrack() - stacks[iDest]->GetSafetyValue1();
  G4long dy2 = stacks[fTurn]->GetNTrack() - stacks[fTurn]->GetSafetyValue2();
  
  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 64 of file G4SmartTrackStack.cc.

{
  for (auto* sp : stacks)
  {
    sp->TransferTo(aStack);
  }
  nTracks = 0;
}

---

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

• G4SmartTrackStack.hh
• G4SmartTrackStack.cc