TargetBremFilter.h

#ifndef BIASING_TARGETBREMFILTER_H
#define BIASING_TARGETBREMFILTER_H
 
//----------------//
//   C++ StdLib   //
//----------------//
#include <algorithm>
 
/*~~~~~~~~~~~~~*/
/*   SimCore   */
/*~~~~~~~~~~~~~*/
#include "SimCore/UserAction.h"
 
/*~~~~~~~~~~~~~~~*/
/*   Framework   */
/*~~~~~~~~~~~~~~~*/
#include "Framework/Configure/Parameters.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 recoilMaxPThreshold_{1500};  // MeV
 
  double bremEnergyThreshold_{2500};
 
  bool killRecoil_{false};
 
};  // TargetBremFilter
}  // namespace biasing
 
#endif  // BIASING_TARGETBREMFILTER_H