GGSPrimaryDisAction Class Reference

Action which finds the point of disappearance points for each primary particle. More...

#include <GGSPrimaryDisAction.h>

Inheritance diagram for GGSPrimaryDisAction:

Collaboration diagram for GGSPrimaryDisAction:

Public Member Functions
	GGSPrimaryDisAction ()
	Constructor.
void	PreUserTrackingAction (const G4Track *track)
	Initialization of primary track. More...
void	PostUserTrackingAction (const G4Track *track)
	Check the track end . More...
void	BeginOfEventAction (const G4Event *event)
	Clears the interaction arrays. More...
void	EndOfEventAction (const G4Event *event)
	Fills the output tree. More...
void	BeginOfRunAction (const G4Run *run)
	Opens the output file for the current run and prepares the output tree. More...
void	EndOfRunAction (const G4Run *run)
	Closes the output file for the current run. More...
Public Member Functions inherited from GGSUserAction
	GGSUserAction ()
	Constructor. More...
virtual	~GGSUserAction ()
	Destructor.
G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *)
	Override of the ClassifyNewTrack method. More...

Detailed Description

Action which finds the point of disappearance points for each primary particle.

This class follows the propagation of primary particles until they disappear from the simulation, either because of an interaction or because they exit from the world volume.

Definition at line 39 of file GGSPrimaryDisAction.h.

Member Function Documentation

void GGSPrimaryDisAction::BeginOfEventAction

(

const G4Event *

event

)

Clears the interaction arrays.

Parameters

event

The current event.

Definition at line 97 of file GGSPrimaryDisAction.cpp.

{ _disInfo.Clear("C"); }

void GGSPrimaryDisAction::BeginOfRunAction

(

const G4Run *

run

)

Opens the output file for the current run and prepares the output tree.

Parameters

run	The current run.

Definition at line 108 of file GGSPrimaryDisAction.cpp.

                                                           {

  _outRootFile = GGSRootFileService::GetInstance().GetFileForThisRun(_outBase, run);
  if (_outTreeName == "") {
    _outTree = GGSRootFileService::GetInstance().GetDefaultTree(_outRootFile);
    _outTree->SetTitle(TString(_outTree->GetTitle()) + "DisInfo ");
  } else
    _outTree = new TTree(_outTreeName.data(), "Disappearance of primary particles");

  // Create inelastic interactions branch
  _outTree->Branch("primaryDisInfo", "TClonesArray", &_disInfo);
}

void GGSPrimaryDisAction::EndOfEventAction

(

const G4Event *

event

)

Fills the output tree.

Parameters

event

The current event.

Definition at line 99 of file GGSPrimaryDisAction.cpp.

                                                          {

  if (_outTreeName != "") {
    // Call fill only if not using the default tree
    _outRootFile->cd();
    _outTree->Fill();
  }
}

void GGSPrimaryDisAction::EndOfRunAction

(

const G4Run *

run

)

Closes the output file for the current run.

Parameters

run	The current run.

Definition at line 121 of file GGSPrimaryDisAction.cpp.

                                                      {

  if (_outTreeName != "") {
    // Write the tree if we are not using the default one
    _outRootFile->cd();
    _outTree->Write();
  }
  GGSRootFileService::GetInstance().CloseFileForThisRun(_outBase);
  _outRootFile = nullptr;
  _outTree = nullptr;
}

void GGSPrimaryDisAction::PostUserTrackingAction

(

const G4Track *

track

)

Check the track end .

For each followed particle, this routine will store information about the termination of the track.

Parameters

track

The current track.

Definition at line 44 of file GGSPrimaryDisAction.cpp.

                                                                     {

  if (_tracks.size() > 0) {
    // See if the track is in the tracking list
    std::map<const G4Track *, int>::iterator trackIter = _tracks.find(track);
    if (trackIter != _tracks.end()) {
      GGSTPrimaryDisInfo *newElement = (GGSTPrimaryDisInfo *)(_disInfo.ConstructedAt(_disInfo.GetEntries()));
      newElement->Clear();

      // Fill the new element
      newElement->primary.PDGCode = track->GetDefinition()->GetPDGEncoding();
      const G4Step *step = track->GetStep();
      newElement->primary.trackID = step->GetTrack()->GetTrackID();
      for (int i = 0; i < 3; i++) {
        newElement->primary.pos[i] = track->GetPosition()[i] / cm;
        newElement->primary.mom[i] = track->GetMomentum()[i] / GeV;
      }
      newElement->primary.time = track->GetGlobalTime() / ns;
      newElement->processName = step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName().data();
      newElement->volumeName = step->GetPostStepPoint()->GetPhysicalVolume()->GetName();
      if (_outputProducts) {
        if (newElement->products) {
          newElement->products->Clear("C");
        } else {
          newElement->products = new TClonesArray("GGSTParticle");
        }
        auto *secondaries = step->GetSecondaryInCurrentStep();
        unsigned int nSecondaries = secondaries->size();
        for (unsigned int is = 0; is < nSecondaries; is++) {
          GGSTParticle *sec = new ((*(newElement->products))[is]) GGSTParticle;
          auto *currentSecondaryTrack = (*secondaries)[is];
          sec->PDGCode = currentSecondaryTrack->GetDefinition()->GetPDGEncoding();
          sec->trackID = currentSecondaryTrack->GetTrackID();
          for (int i = 0; i < 3; i++) {
            sec->pos[i] = currentSecondaryTrack->GetPosition()[i] / cm;
            sec->mom[i] = currentSecondaryTrack->GetMomentum()[i] / GeV;
          }
          sec->time = currentSecondaryTrack->GetGlobalTime() / ns;
        }

      } else {
        delete newElement->products;
        newElement->products = nullptr;
      } // if(_outputProducts)

      // Delete the currently followed particle from list
      _tracks.erase(trackIter);

    } // if (trackIter != _tracks.end())

  } // if (_tracks.size() > 0)
}

void GGSPrimaryDisAction::PreUserTrackingAction

(

const G4Track *

track

)

Initialization of primary track.

This routine checks if the current track belongs to a primary particle, and in affirmative case it adds it to the followed particles vector

Parameters

track

The current track.

Definition at line 35 of file GGSPrimaryDisAction.cpp.

                                                                    {

  // Check whether it is a primary particle
  G4PrimaryParticle *primary = track->GetDynamicParticle()->GetPrimaryParticle();
  if (primary != nullptr) {
    _tracks[track] = track->GetTrackID();
  }
}

---

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

• GGSPrimaryDisAction.h
• GGSPrimaryDisAction.cpp