NeutronInelastic.h

#ifndef SIMCORE_BIASOPERATORS_NEUTRONINELASTIC_H_
#define SIMCORE_BIASOPERATORS_NEUTRONINELASTIC_H_
 
#include "SimCore/XsecBiasingOperator.h"
 
namespace simcore {
namespace biasoperators {
 
class NeutronInelastic : public XsecBiasingOperator {
 public:
  NeutronInelastic(std::string name, const framework::config::Parameters& p);
 
  virtual ~NeutronInelastic() = default;
 
  G4VBiasingOperation* ProposeOccurenceBiasingOperation(
      const G4Track* track,
      const G4BiasingProcessInterface* callingProcess) override;
 
  std::string getProcessToBias() const override { return "neutronInelastic"; }
 
  std::string getParticleToBias() const override { return "neutron"; }
 
  std::string getVolumeToBias() const override { return volume_; }
 
  void RecordConfig(ldmx::RunHeader& header) const override {
    header.setStringParameter("BiasOperator::NeutronInelastic::Volume",
                              volume_);
    header.setFloatParameter("BiasOperator::NeutronInelastic::Factor", factor_);
    header.setFloatParameter("BiasOperator::NeutronInelastic::Threshold",
                             threshold_);
  }
 
 private:
  std::string volume_;
 
  double factor_;
 
  double threshold_;
};
 
}  // namespace biasoperators
}  // namespace simcore
 
#endif  // SIMCORE_BIASOPERATORS_NEUTRONINELASTIC_H_