RE01StackingAction Class Reference

#include <RE01StackingAction.hh>

Inheritance diagram for RE01StackingAction:

Public Member Functions
	RE01StackingAction ()
virtual	~RE01StackingAction ()
virtual G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *aTrack)
virtual void	NewStage ()
virtual void	PrepareNewEvent ()
Public Member Functions inherited from G4UserStackingAction
	G4UserStackingAction ()
virtual	~G4UserStackingAction ()
void	SetStackManager (G4StackManager *value)

Additional Inherited Members
Protected Attributes inherited from G4UserStackingAction
G4StackManager *	stackManager

Detailed Description

Definition at line 41 of file RE01StackingAction.hh.

Constructor & Destructor Documentation

RE01StackingAction::RE01StackingAction

(

)

Definition at line 48 of file RE01StackingAction.cc.

  :G4UserStackingAction(),
   fStage(0),fCalorimeterHitsColID(-1)
{;}

RE01StackingAction::~RE01StackingAction ( )

virtual

Definition at line 54 of file RE01StackingAction.cc.

{;}

Member Function Documentation

G4ClassificationOfNewTrack RE01StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 59 of file RE01StackingAction.cc.

References fSuspend, fUrgent, fWaiting, G4Track::GetParentID(), G4Track::GetTrackStatus(), G4Track::GetUserInformation(), RE01TrackInformation::SetSourceTrackInformation(), RE01TrackInformation::SetTrackingStatus(), and G4Track::SetUserInformation().

{
  G4ClassificationOfNewTrack classification = fUrgent;

  if(fStage==0)
  {
    RE01TrackInformation* trackInfo;
    if(aTrack->GetTrackStatus()==fSuspend) // Track reached to calorimeter
    {
      trackInfo = (RE01TrackInformation*)(aTrack->GetUserInformation());
      trackInfo->SetTrackingStatus(0);
      trackInfo->SetSourceTrackInformation(aTrack);
//      G4cout << "Track " << aTrack->GetTrackID() 
//             << " (parentID " << aTrack->GetParentID() 
//             << ") has reached to calorimeter and has been suspended at " 
//             << aTrack->GetPosition() << G4endl;
      classification = fWaiting;
    }
    else if(aTrack->GetParentID()==0) // Primary particle
    {
      trackInfo = new RE01TrackInformation(aTrack);
      trackInfo->SetTrackingStatus(1);
      G4Track* theTrack = (G4Track*)aTrack;
      theTrack->SetUserInformation(trackInfo);
    }
  }
  return classification;
}

void RE01StackingAction::NewStage ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 105 of file RE01StackingAction.cc.

References G4THitsCollection< T >::entries(), fUrgent, fWaiting, G4cout, G4endl, G4StackManager::GetNUrgentTrack(), python.hepunit::GeV, G4UserStackingAction::stackManager, G4StackManager::TransferOneStackedTrack(), and G4StackManager::TransferStackedTracks().

{
  G4cout << "+++++++++++ Stage " << fStage << G4endl;
  if(fStage==0)
  {
    // display trajetory information in the tracking region
    G4cout << G4endl;
    G4cout << "Tracks in tracking region have been processed. -- Stage 0 over." 
           << G4endl;
    G4cout << G4endl;
  }
  else
  {
    // display calorimeter information caused by a "source" track 
    // in the tracker region
    //    G4cout << G4endl;
    //    G4cout << "Processing one shower originated by a source track "
    //           << "in tracker region is over." << G4endl;
    //    G4cout << G4endl;
    RE01CalorimeterHitsCollection* CHC = 
      (RE01CalorimeterHitsCollection*)GetCalCollection();
    if(CHC)
    { 
      int n_hit = CHC->entries();
      G4double totE = 0;
      G4int n_hitByATrack = 0;
      for(int i=0;i<n_hit;i++)
      {
        G4double edepByATrack = (*CHC)[i]->GetEdepByATrack();
        if(edepByATrack>0.)
        {
          totE += edepByATrack;
          if(n_hitByATrack==0)
          { (*CHC)[i]->GetTrackInformation()->Print(); }
          n_hitByATrack++;
          G4cout << "Cell[" << (*CHC)[i]->GetZ() << "," 
                 << (*CHC)[i]->GetPhi() << "]    " 
                 << edepByATrack/GeV << " [GeV]" << G4endl;
          (*CHC)[i]->ClearEdepByATrack();
        }
      }
      if(n_hitByATrack>0)
      {
        G4cout << 
          "###  Total energy deposition in calorimeter by a source track in "
               << n_hitByATrack << " cells : " << totE / GeV << " (GeV)" 
               << G4endl;
        G4cout << G4endl;
      }
    }
  }

  if(stackManager->GetNUrgentTrack())
  {
    // Transfer all tracks in Urgent stack to Waiting stack, since all tracks
    // in Waiting stack have already been transfered to Urgent stack before
    // invokation of this method.
    stackManager->TransferStackedTracks(fUrgent,fWaiting);

    // Then, transfer only one track to Urgent stack.
    stackManager->TransferOneStackedTrack(fWaiting,fUrgent);

    fStage++;
  }
}

void RE01StackingAction::PrepareNewEvent ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 172 of file RE01StackingAction.cc.

{ 
  fStage = 0; 
}

---

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

• RE01StackingAction.hh
• RE01StackingAction.cc