TargetBremFilter.h

#ifndef BIASING_TARGETBREMFILTER_H
#define BIASING_TARGETBREMFILTER_H
 
/*~~~~~~~~~~~~~*/
/*   SimCore   */
/*~~~~~~~~~~~~~*/
#include "SimCore/G4User/UserAction.h"
 
/*~~~~~~~~~~~~~~~*/
/*   Framework   */
/*~~~~~~~~~~~~~~~*/
#include "Framework/Configure/Parameters.h"
#include "Framework/Logger.h"
 
namespace biasing {
 
class TargetBremFilter : public simcore::UserAction {
 public:
  TargetBremFilter(const std::string& name,
                   framework::config::Parameters& parameters);
 
  ~TargetBremFilter();
 
  void stepping(const G4Step* step) override;
 
  void endOfEventAction(const G4Event*) override;
 
  G4ClassificationOfNewTrack classifyNewTrack(
      const G4Track* aTrack,
      const G4ClassificationOfNewTrack& currentTrackClass) override;
 
  std::vector<simcore::TYPE> getTypes() override {
    return {simcore::TYPE::EVENT, simcore::TYPE::STACKING,
            simcore::TYPE::STEPPING};
  }
 
 private:
  double recoil_max_p_threshold_{1500};  // MeV
 
  double brem_energy_threshold_{2500};
 
  double brem_theta_min_;
  double brem_theta_max_;
 
  double dral_min_;
  double dral_max_;
 
  bool kill_recoil_{false};
 
};  // TargetBremFilter
}  // namespace biasing
 
#endif  // BIASING_TARGETBREMFILTER_H