gdetectorConstruction.h

Go to the documentation of this file.

#pragma once
 
// geant4
#include "G4VUserDetectorConstruction.hh"
 
// gemc
#include "goptions.h"
#include "g4world.h"
#include "gmagneto.h"
#include "gsd.h"
#include "gbase.h"
 
// c++
#include <vector>
 
class GDetectorConstruction : public GBase<GDetectorConstruction>, public G4VUserDetectorConstruction {
public:
    explicit GDetectorConstruction(std::shared_ptr<GOptions> gopts);
 
public:
    G4VPhysicalVolume* Construct() override;
 
    void ConstructSDandField() override;
 
    [[nodiscard]] bool is_empty() const { return g4world->is_empty(); }
 
    void reload_geometry(SystemList sl);
 
    std::shared_ptr<GDynamicDigitization> get_digitization_routines_for_sdname(const std::string &sd_name) const {
        return digitization_routines_map->at(sd_name);
    }
 
    std::shared_ptr<gdynamicdigitization::dRoutinesMap> get_digitization_routines_map() const {
        return digitization_routines_map;
    }
 
    [[nodiscard]] std::unordered_map<std::string, G4Volume*> get_g4volumes_map() const {
        if (g4world == nullptr) {
            log->error(ERR_G4VOLUMEBUILDFAILED, "No g4world ");
        }
        return g4world->get_g4volumes_map();
    }
 
private:
    std::shared_ptr<GOptions> gopt;    // need options inside Constructs() methods
 
    std::shared_ptr<GWorld>   gworld;
 
    std::shared_ptr<G4World>  g4world;
 
    std::shared_ptr<gdynamicdigitization::dRoutinesMap> digitization_routines_map;
 
    static G4ThreadLocal GMagneto* gmagneto;
 
    void loadDigitizationPlugins();
 
    SystemList gsystems;
};