DetectorConstruction.cxx

/*~~~~~~~~~~~~~~~*/
/*   Framework   */
/*~~~~~~~~~~~~~~~*/
#include "Framework/Exception/Exception.h"
 
/*~~~~~~~~~~~~~*/
/*   SimCore   */
/*~~~~~~~~~~~~~*/
#include "SimCore/DetectorConstruction.h"
#include "SimCore/G4User/VolumeChecks.h"
#include "SimCore/SensitiveDetector.h"
#include "SimCore/XsecBiasingOperator.h"
 
namespace simcore {
namespace logical_volume_tests {
 
using Test = bool (*)(G4LogicalVolume*, const std::string&);
 
}  // namespace logical_volume_tests
 
DetectorConstruction::DetectorConstruction(
    simcore::geo::Parser* parser, framework::config::Parameters& parameters,
    ConditionsInterface& ci)
    : parser_(parser), parameters_{parameters}, conditions_interface_{ci} {}
DetectorConstruction::DetectorConstruction( {…}
 
G4VPhysicalVolume* DetectorConstruction::Construct() {
  return parser_->GetWorldVolume();
}
G4VPhysicalVolume* DetectorConstruction::Construct() {…}
 
void DetectorConstruction::ConstructSDandField() {
  auto sens_dets{
      parameters_.getParameter<std::vector<framework::config::Parameters>>(
          "sensitive_detectors", {})};
  for (auto& det : sens_dets) {
    // create
    auto sd = SensitiveDetector::Factory::get().make(
        det.getParameter<std::string>("class_name"),
        det.getParameter<std::string>("instance_name"), conditions_interface_,
        det);
    // attach to volumes
    for (G4LogicalVolume* volume : *G4LogicalVolumeStore::GetInstance()) {
      if (sd->isSensDet(volume)) {
        ldmx_log(debug) << "Attaching " << sd->GetName() << " to "
                        << volume->GetName();
        volume->SetSensitiveDetector(sd);
      }
    }
  }
 
  // Biasing operators were created in RunManager::setupPhysics
  //  which is called before G4RunManager::Initialize
  //  which is where this method ends up being called.
  simcore::XsecBiasingOperator::Factory::get().apply([&](auto bop) {
    logical_volume_tests::Test includeVolumeTest{nullptr};
    if (bop->getVolumeToBias().compare("ecal") == 0) {
      includeVolumeTest = &simcore::g4user::volumechecks::isInEcal;
    } else if (bop->getVolumeToBias().compare("old_ecal") == 0) {
      includeVolumeTest = &simcore::g4user::volumechecks::isInEcalOld;
    } else if (bop->getVolumeToBias().compare("target") == 0) {
      includeVolumeTest = &simcore::g4user::volumechecks::isInTargetOnly;
    } else if (bop->getVolumeToBias().compare("target_region") == 0) {
      includeVolumeTest = &simcore::g4user::volumechecks::isInTargetRegion;
    } else if (bop->getVolumeToBias().compare("hcal") == 0) {
      includeVolumeTest = &simcore::g4user::volumechecks::isInHcal;
    } else {
      std::cerr << "[ DetectorConstruction ] : "
                << "WARN - Requested volume to bias '" << bop->getVolumeToBias()
                << "' is not recognized. Will attach volumes based on if their"
                << " name contains the volume to bias." << std::endl;
      includeVolumeTest = &simcore::g4user::volumechecks::nameContains;
    }
 
    for (G4LogicalVolume* volume : *G4LogicalVolumeStore::GetInstance()) {
      auto volume_name = volume->GetName();
      if (includeVolumeTest(volume, bop->getVolumeToBias())) {
        bop->AttachTo(volume);
        ldmx_log(debug) << "Attaching biasing operator " << bop->GetName()
                        << " to volume " << volume->GetName();
      }  // BOP attached to target or ecal
    }    // loop over volumes
  });    // loop over biasing operators
}
}  // namespace simcore