PhotoNuclear.h

 
#ifndef SIMCORE_BIASOPERATORS_PHOTONUCLEAR_H_
#define SIMCORE_BIASOPERATORS_PHOTONUCLEAR_H_
 
#include "SimCore/BiasOperators/XsecBiasingOperator.h"
#include "SimCore/G4User/PtrRetrieval.h"
 
namespace simcore {
namespace biasoperators {
 
class PhotoNuclear : public XsecBiasingOperator {
 public:
  PhotoNuclear(std::string name, const framework::config::Parameters& p);
 
  void StartRun() override;
 
  G4VBiasingOperation* ProposeOccurenceBiasingOperation(
      const G4Track* track,
      const G4BiasingProcessInterface* callingProcess) override;
 
  std::string getProcessToBias() const override { return "photonNuclear"; }
 
  std::string getParticleToBias() const override { return "gamma"; }
 
  std::string getVolumeToBias() const override { return volume_; }
 
  void RecordConfig(ldmx::RunHeader& h) const override {
    h.setStringParameter("BiasOperators::PhotoNuclear::Volume", volume_);
    h.setFloatParameter("BiasOperators::PhotoNuclear::Threshold", threshold_);
    h.setFloatParameter("BiasOperators::PhotoNuclear::Factor", factor_);
    h.setIntParameter("BiasOperators::PhotoNuclear::Bias Conv Down",
                      down_bias_conv_);
    h.setIntParameter("BiasOperators::PhotoNuclear::Only Children Of Primary",
                      only_children_of_primary_);
  }
 
 private:
  static const std::string CONVERSION_PROCESS;
 
  G4BOptnChangeCrossSection* em_xsec_operation_{nullptr};
 
  double pn_xsec_unbiased_{0};
 
  double pn_xsec_biased_{0};
 
  std::string volume_;
 
  double threshold_;
 
  double factor_;
 
  bool down_bias_conv_;
 
  bool only_children_of_primary_;
 
};  // PhotoNuclear
}  // namespace biasoperators
}  // namespace simcore
 
#endif  // SIMCORE_BIASOPERATORS_PHOTONUCLEAR_H_