PhotoNuclearTopologyFilters.h

#ifndef PHOTONUCLEARTOPOLOGYFILTERS_H
#define PHOTONUCLEARTOPOLOGYFILTERS_H
 
/*~~~~~~~~~~~~~*/
/*   SimCore   */
/*~~~~~~~~~~~~~*/
#include "SimCore/UserAction.h"
#include "SimCore/UserTrackInformation.h"
/*~~~~~~~~~~~~~~~*/
/*   Framework   */
/*~~~~~~~~~~~~~~~*/
#include <G4RunManager.hh>
#include <G4Step.hh>
 
#include "Framework/Configure/Parameters.h"
// Forward declaration
 
namespace biasing {
 
class PhotoNuclearTopologyFilter : public simcore::UserAction {
 public:
  PhotoNuclearTopologyFilter(const std::string& name,
                             framework::config::Parameters& parameters);
 
  ~PhotoNuclearTopologyFilter() = default;
 
  void stepping(const G4Step* step) override;
 
  virtual bool rejectEvent(const std::vector<G4Track*>& secondaries) const = 0;
 
  std::vector<simcore::TYPE> getTypes() override {
    return {simcore::TYPE::STEPPING};
  }
 
 protected:
  constexpr bool isLightIon(const int pdgCode) const {
    if (pdgCode > 1000000000) {
      // Check if the atomic number is less than or equal to 4
      return ((pdgCode / 10) % 1000) <= 4;
    }
    return false;
  }
 
  constexpr bool skipCountingParticle(const int pdgcode) const {
    return !(pdgcode < 10000 || (count_light_ions_ && isLightIon(pdgcode)));
  }
 
  constexpr bool isNeutron(const int pdgID) const { return pdgID == 2112; }
 
  bool count_light_ions_;
  double hard_particle_threshold_;
};  // PhotoNuclearTopologyFilter
 
class SingleNeutronFilter : public PhotoNuclearTopologyFilter {
 public:
  SingleNeutronFilter(const std::string& name,
                      framework::config::Parameters& parameters)
      : PhotoNuclearTopologyFilter{name, parameters} {}
 
  bool rejectEvent(const std::vector<G4Track*>& secondaries) const override;
};
 
class NothingHardFilter : public PhotoNuclearTopologyFilter {
 public:
  NothingHardFilter(const std::string& name,
                    framework::config::Parameters& parameters)
      : PhotoNuclearTopologyFilter{name, parameters} {}
 
  bool rejectEvent(const std::vector<G4Track*>& secondaries) const override;
};
}  // namespace biasing
 
#endif /* PHOTONUCLEARTOPOLOGYFILTERS_H */