MidShowerNuclearBkgdFilter.h

#ifndef BIASING_MIDSHOWERBKGDFILTER_H_
#define BIASING_MIDSHOWERBKGDFILTER_H_
 
/*~~~~~~~~~~~~*/
/*   SimCore  */
/*~~~~~~~~~~~~*/
#include "SimCore/UserAction.h"
 
class G4VProcess;
 
namespace biasing {
 
class MidShowerNuclearBkgdFilter : public simcore::UserAction {
 public:
  MidShowerNuclearBkgdFilter(const std::string& name,
                             framework::config::Parameters& parameters);
 
  ~MidShowerNuclearBkgdFilter() {}
 
  std::vector<simcore::TYPE> getTypes() override {
    return {simcore::TYPE::STACKING, simcore::TYPE::STEPPING,
            simcore::TYPE::EVENT};
  }
 
  void BeginOfEventAction(const G4Event* event) override;
 
  void stepping(const G4Step* step) override;
 
  void NewStage() override;
 
 private:
  bool isOutsideCalorimeterRegion(const G4Step* step) const;
 
  bool isNuclearProcess(const G4VProcess* proc) const;
 
  void save(const G4Track* track) const;
 
  void AbortEvent(const std::string& reason) const;
 
 private:
  double threshold_;
 
  std::vector<std::string> nuclear_processes_;
 
  double total_process_energy_{0.};
 
};  // MidShowerNuclearBkgdFilter
}  // namespace biasing
 
#endif  // BIASING_MIDSHOWERBKGDFILTER_H__