RE01StackingAction Class Reference

#include <RE01StackingAction.hh>

Inheritance diagram for RE01StackingAction:

Collaboration 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)

Private Member Functions
G4VHitsCollection *	GetCalCollection ()

Private Attributes
G4int	fStage
G4int	fCalorimeterHitsColID

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::RE01StackingAction

(

)

Definition at line 48 of file RE01StackingAction.cc.

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

~RE01StackingAction()

RE01StackingAction::~RE01StackingAction ( )

virtual

Definition at line 54 of file RE01StackingAction.cc.

{;}

Here is the call graph for this function:

Member Function Documentation

ClassifyNewTrack()

G4ClassificationOfNewTrack RE01StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 59 of file RE01StackingAction.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetCalCollection()

G4VHitsCollection * RE01StackingAction::GetCalCollection ( )

private

Definition at line 89 of file RE01StackingAction.cc.

{
  G4SDManager* SDMan = G4SDManager::GetSDMpointer();
  G4RunManager* runMan = G4RunManager::GetRunManager();
  if(fCalorimeterHitsColID<0)
  { fCalorimeterHitsColID = SDMan->GetCollectionID("calCollection"); }
  if(fCalorimeterHitsColID>=0)
  {
    const G4Event* currentEvent = runMan->GetCurrentEvent();
    G4HCofThisEvent* HCE = currentEvent->GetHCofThisEvent();
    return HCE->GetHC(fCalorimeterHitsColID);
  }
  return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

NewStage()

void RE01StackingAction::NewStage ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 105 of file RE01StackingAction.cc.

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

Here is the call graph for this function:

PrepareNewEvent()

void RE01StackingAction::PrepareNewEvent ( )

virtual

Reimplemented from G4UserStackingAction.

Definition at line 172 of file RE01StackingAction.cc.

{ 
  fStage = 0; 
}

Member Data Documentation

fCalorimeterHitsColID

G4int RE01StackingAction::fCalorimeterHitsColID

private

Definition at line 56 of file RE01StackingAction.hh.

fStage

G4int RE01StackingAction::fStage

private

Definition at line 55 of file RE01StackingAction.hh.

---

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

• RE01StackingAction.hh
• RE01StackingAction.cc