GenieGenerator.h

Go to the documentation of this file.

 
#ifndef SIMCORE_GENIE_GENERATOR_H
#define SIMCORE_GENIE_GENERATOR_H
 
//----------//
//   ROOT   //
//----------//
#include "TRandom.h"
#include "TRandomGen.h"
 
//------------//
//   GENIE   //
//------------//
#include <string>
#include <vector>
 
#include "GENIE/Framework/EventGen/GEVGDriver.h"
#include "GENIE/Framework/EventGen/HepMC3Converter.h"
 
//------------//
//   LDMX     //
//------------//
#include "Framework/Logger.h"
#include "SimCore/G4User/UserEventInformation.h"
#include "SimCore/Generators/PrimaryGenerator.h"
 
// Forward declarations
class G4Event;
 
namespace simcore {
namespace generators {
 
class GenieGenerator : public simcore::PrimaryGenerator {
 public:
  GenieGenerator(const std::string& name,
                 const framework::config::Parameters& parameters);
 
  ~GenieGenerator();
 
  void GeneratePrimaryVertex(G4Event* event) final override;
 
  void RecordConfig(const std::string& id, ldmx::RunHeader& rh) final override;
 
 private:
  std::vector<genie::GEVGDriver> evg_drivers_;
  genie::HepMC3Converter hep_mc3_converter_;
 
  double energy_;
  std::vector<int> targets_;
  std::vector<double> abundances_;
  std::vector<double> position_;
  std::vector<double> beam_size_;
  float target_thickness_;
  float time_;
  std::vector<double> direction_;
 
  std::string tune_;
  std::string spline_file_;
 
  std::string message_threshold_file_;
 
  std::vector<float> ev_weighting_integral_;
  size_t n_events_generated_;
  std::vector<size_t> n_events_by_target_;
  std::vector<float> xsec_by_target_;
 
  float xsec_total_;
 
  void fillConfig(
      const framework::config::Parameters&);  
  bool validateConfig();  
 
  void initializeGENIE();   
  void calculateTotalXS();  
 
  enableLogging("GenieGenerator")
 
};  // GenieGenerator
 
}  // namespace generators
}  // namespace simcore
 
#endif  // SIMCORE_GENIE_GENERATOR_H