AthenaStackingAction.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef G4ATLASALG_ATHENASTACKINGACTION_H
#define G4ATLASALG_ATHENASTACKINGACTION_H
 
#include "G4UserStackingAction.hh"
 
class PrimaryParticleInformation;
 
namespace G4UA
{

  class AthenaStackingAction : public G4UserStackingAction
  {
 
    public:

      struct Config
      {
        bool killAllNeutrinos;
        double photonEnergyCut;
        bool applyNRR;
        double russianRouletteNeutronThreshold;
        double russianRouletteNeutronWeight;
        bool applyPRR;
        double russianRoulettePhotonThreshold;
        double russianRoulettePhotonWeight;
        bool isISFJob;
      };
      
      AthenaStackingAction(const Config& config);

      virtual G4ClassificationOfNewTrack
      ClassifyNewTrack(const G4Track* track) override;
 
    protected:

      Config m_config;

      bool isNeutrino(const G4Track*) const;

      bool isGamma(const G4Track*) const;

      bool isNeutron(const G4Track*) const;

      PrimaryParticleInformation* getPrimaryParticleInformation(const G4Track *track) const;
 
      // one over m_config.russianRouletteNeutronWeight
      double m_oneOverWeightNeutron;
      
      // one over m_config.russianRoulettePhotonWeight
      double m_oneOverWeightPhoton;
 
  }; // class AthenaStackingAction
 
} // namespace G4UA
 
#endif