GenieElectroNuclearProcess.h

Go to the documentation of this file.

 
#ifndef SIMCORE_GENIEELECTRONUCLEARPROCESS_H
#define SIMCORE_GENIEELECTRONUCLEARPROCESS_H
 
#include <map>
#include <memory>
#include <set>
#include <string>
#include <vector>
 
#include "Framework/Configure/Parameters.h"
#include "Framework/Logger.h"
#include "G4VDiscreteProcess.hh"
#include "GENIE/Framework/EventGen/GEVGDriver.h"
#include "GENIE/Framework/EventGen/HepMC3Converter.h"
 
class G4Element;
 
namespace simcore {
 
class GenieElectroNuclearProcess : public G4VDiscreteProcess {
 public:
  static const std::string PROCESS_NAME;
 
  GenieElectroNuclearProcess(const framework::config::Parameters& params);
 
  ~GenieElectroNuclearProcess() override;
 
  G4bool IsApplicable(const G4ParticleDefinition& p) override;
 
  G4VParticleChange* PostStepDoIt(const G4Track& track,
                                  const G4Step& step) override;
 
  void setOnlyOnePerEvent(bool val) { only_one_per_event_ = val; }
 
 protected:
  G4double GetMeanFreePath(const G4Track& track, G4double prevStepSize,
                           G4ForceCondition* condition) override;
 
 private:
  void initializeGENIE();
 
  void loadAvailableTargets();
 
  bool splineAvailable(int target_code) const;
 
  void setupDrivers();
 
  void discoverIsotopesForElement(const G4Element* element);
 
  void discoverFromVolume();
 
  std::vector<std::unique_ptr<genie::GEVGDriver>> evg_drivers_;
 
  genie::HepMC3Converter hep_mc3_converter_;
 
  std::vector<int> targets_;
 
  std::vector<double> abundances_;
 
  std::string tune_;
 
  std::string spline_file_;
 
  std::string message_threshold_file_;
 
  bool only_one_per_event_{true};
 
  bool genie_initialized_{false};
 
  bool auto_discover_{false};
 
  std::string discover_volume_;
 
  std::set<int> discovered_z_;
 
  std::set<int> available_targets_;
 
  std::map<int, std::vector<std::pair<int, double>>> z_to_targets_;
 
  std::vector<double> partial_sum_sigma_;
 
  enableLogging("GenieElectroNuclearProcess")
};
 
}  // namespace simcore
 
#endif  // SIMCORE_GENIEELECTRONUCLEARPROCESS_H