G4UA::FixG4CreatorProcess Class Reference

#include <FixG4CreatorProcess.h>

Inheritance diagram for G4UA::FixG4CreatorProcess:

Collaboration diagram for G4UA::FixG4CreatorProcess:

Public Member Functions
	FixG4CreatorProcess ()
virtual void	UserSteppingAction (const G4Step *) override final

Detailed Description

Definition at line 14 of file FixG4CreatorProcess.h.

Constructor & Destructor Documentation

FixG4CreatorProcess()

G4UA::FixG4CreatorProcess::FixG4CreatorProcess

(

)

Definition at line 18 of file FixG4CreatorProcess.cxx.

{}

Member Function Documentation

UserSteppingAction()

void G4UA::FixG4CreatorProcess::UserSteppingAction ( const G4Step *  aStep )

finaloverridevirtual

Definition at line 22 of file FixG4CreatorProcess.cxx.

  {
    //
    //  M. Novak, M. Bandieramonte: 10 November 2021
    //  
    //  Make sure that the creator process of the secondary tracks of the current step
    //  is the same as the one that limited the step. These can be different in case of
    //  using `wrapper` processes e.g. G4GammaGeneralProcess or G4HepEm, etc.
    if (G4EmParameters::Instance()->GeneralProcessActive()) {
        const G4ParticleDefinition* thePrimaryParticle = aStep->GetTrack()->GetParticleDefinition();
        const std::size_t           numThisSecondaries = aStep->GetNumberOfSecondariesInCurrentStep();
        if (thePrimaryParticle==G4Gamma::GammaDefinition() && numThisSecondaries>0) {
                // Get the pointer to the process that limited the step: i.e. the one that
                // created the secondaries of the current step
            const G4VProcess* theLimiterProcess = aStep->GetPostStepPoint()->GetProcessDefinedStep();
            // note: this is a vector of secondaries containing all secondaries created
            // along the tracking of the current `primary` track (i.e. not only
            // secondaries created in this step)
            const G4TrackVector* theSecTrackVector = aStep->GetSecondary();
            const std::size_t    numAllSecondaries = theSecTrackVector->size();
            for (std::size_t it = numAllSecondaries-numThisSecondaries; it<numAllSecondaries; ++it) {
 
                G4Track* secTrack = (*theSecTrackVector)[it];
                if (secTrack->GetCreatorProcess()!=theLimiterProcess) {
 
                    secTrack->SetCreatorProcess(theLimiterProcess);
                }
                }
            }
 
    }   
 }

---

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

• FixG4CreatorProcess.h
• FixG4CreatorProcess.cxx