Fix segmentation fault due to G4 bug by setting G4MaterialPropertiesTable

source/geometries/GenericPhotosensor.cc

@@ -83,6 +83,7 @@ void GenericPhotosensor::DefineMaterials()
 {
   // Case /////
   case_mat_ = materials::FR4();
+  case_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Window /////
   // The optical properties of a given material, is common for the whole geometry so,
@@ -119,6 +120,7 @@ void GenericPhotosensor::DefineMaterials()
 
   // Sensitive /////
   sensitive_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Si");
+  sensitive_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // WLS coating /////
   wls_mat_ = materials::TPB();

source/geometries/Next1EL.cc

@@ -255,19 +255,25 @@ void Next1EL::DefineMaterials()
 {
   // AIR
   air_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR");
+  air_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   // GASEOUS XENON
   gxe_ = materials::GXe(pressure_, 303);
   gxe_->SetMaterialPropertiesTable(opticalprops::GXe(pressure_, 303, sc_yield_, e_lifetime_));
   // PTFE (TEFLON)
   teflon_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+  teflon_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   // STAINLESS STEEL
   steel_ = materials::Steel();
+  steel_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   // ALUMINUM
   aluminum_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Al");
+  aluminum_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   //LEAD
   lead_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Pb");
+  lead_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   //PLASTIC
   plastic_ = materials::PS();
+  plastic_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   //TPB
   tpb_ = materials::TPB();
   tpb_->SetMaterialPropertiesTable(opticalprops::TPB());
@@ -352,6 +358,7 @@ void Next1EL::BuildExtScintillator()
     new G4Tubs("SOURCE", 0., source_diam/2., source_thick/2., 0., twopi);
   G4Material* sodium22_mat =
     G4NistManager::Instance()->FindOrBuildMaterial("G4_Na");
+  sodium22_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   G4LogicalVolume* source_logic =
     new G4LogicalVolume(source_solid, sodium22_mat, "SOURCE");
 
@@ -400,6 +407,7 @@ void Next1EL::BuildExtScintillator()
   G4Tubs* sc_solid = new G4Tubs("NaI", 0., radius, length/2., 0., twopi);
   G4Material* mat =
     G4NistManager::Instance()->FindOrBuildMaterial("G4_SODIUM_IODIDE");
+  mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
   G4LogicalVolume* sc_logic = new G4LogicalVolume(sc_solid, mat, "NaI");
   sc_logic->SetUserLimits(new G4UserLimits(1.*mm));
 
@@ -955,8 +963,11 @@ void Next1EL::BuildFieldCage()
     new G4UnionSolid("SUPPORT_BAR", bar_base, addon, 0,
 		     G4ThreeVector(bar_thickn_, 0., (bar_length - bar_addon_length_)/2.));
 
+  G4Material* peek_mat = materials::PEEK();
+  peek_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
   G4LogicalVolume* bar_logic =
-    new G4LogicalVolume(bar_solid, materials::PEEK(), "SUPPORT_BAR");
+    new G4LogicalVolume(bar_solid, peek_mat, "SUPPORT_BAR");
 
 
   G4double pos_rad = ring_diam_/2. + ring_thickn_ + bar_thickn_/2.;

source/geometries/Next1ELDBO.cc

@@ -65,6 +65,7 @@ namespace nexus {
 
     G4Material* teflon =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+    teflon->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4LogicalVolume* board_logic =
       new G4LogicalVolume(board_solid, teflon, "DBO");

source/geometries/NextDemo.cc

@@ -75,9 +75,14 @@ void NextDemo::ConstructLab()
   G4Box* lab_solid_vol =
     new G4Box(lab_name, lab_size_/2., lab_size_/2., lab_size_/2.);
 
+  G4Material* air_mat = 
+    G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR");
+
+  air_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
   G4LogicalVolume* lab_logic_vol =
     new G4LogicalVolume(lab_solid_vol,
-                        G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR"),
+                        air_mat,
                         lab_name);
 
   // Make the volume invisible in the visualization

source/geometries/NextDemoEnergyPlane.cc

@@ -120,9 +120,12 @@ namespace nexus {
                                pmt_positions_[i]);
     }
 
+    G4Material* aluminum = G4NistManager::Instance()->FindOrBuildMaterial("G4_Al");
+    aluminum->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* carrier_plate_logic =
       new G4LogicalVolume(carrier_plate_solid,
-                          G4NistManager::Instance()->FindOrBuildMaterial("G4_Al"),
+                          aluminum,
                           "EP_PLATE");
     G4double carrier_plate_posz =
       GetELzCoord() + gate_support_surface_dist_ + carrier_plate_thickness_/2.;
@@ -139,6 +142,7 @@ namespace nexus {
     sapphire->SetMaterialPropertiesTable(opticalprops::Sapphire());
 
     G4Material* copper = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+    copper->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4Material* vacuum =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_Galactic");
@@ -247,9 +251,13 @@ namespace nexus {
     G4Tubs* pmt_base_solid =
       new G4Tubs("PMT_BASE", 0., pmt_base_diam_/2., pmt_base_thickness_,
                  0.,twopi);
+
+    G4Material* kapton = G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON");
+    kapton->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* pmt_base_logic =
       new G4LogicalVolume(pmt_base_solid,
-                          G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON"),
+                          kapton,
 			  "PMT_BASE");
     G4double pmt_base_pos = pmt_hole_length_/2. - pmt_base_thickness_;
     new G4PVPlacement(0, G4ThreeVector(0.,0., pmt_base_pos),

source/geometries/NextDemoFieldCage.cc

@@ -225,10 +225,12 @@ namespace nexus {
     temperature_ = gas_->GetTemperature();
 
     aluminum_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Al");
+    aluminum_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     /// Teflon for the light tube
     teflon_ =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+    teflon_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     /// Quartz
     quartz_ =  materials::FusedSilica();

source/geometries/NextDemoSiPMBoard.cc

@@ -179,9 +179,12 @@ void NextDemoSiPMBoard::Construct()
   G4Box* kapton_solid = new G4Box(kapton_name, board_size_x/2.,
                                   board_size_y/2., kapton_thickn_/2.);
 
+  G4Material* kapton_mat = G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON");
+  kapton_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
   G4LogicalVolume* kapton_logic =
     new G4LogicalVolume(kapton_solid,
-                        G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON"),
+                        kapton_mat,
                         kapton_name);
 
   new G4PVPlacement(nullptr, G4ThreeVector(0., 0., kapton_posz), kapton_logic,
@@ -196,8 +199,11 @@ void NextDemoSiPMBoard::Construct()
   G4Box* mask_solid = new G4Box(mask_name, board_size_x/2.,
                                 board_size_y/2., mask_thickn_/2.);
 
+  G4Material* teflon = G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+  teflon->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
   G4LogicalVolume* mask_logic =
-    new G4LogicalVolume(mask_solid, G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON"),
+    new G4LogicalVolume(mask_solid, teflon,
                         mask_name);
 
   // Adding the optical surface

source/geometries/NextElDB.cc

@@ -67,6 +67,8 @@ namespace nexus {
     G4Material* teflon =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
 
+    teflon->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* board_logic =
       new G4LogicalVolume(board_solid, teflon, "DICE_BOARD");
 

source/geometries/NextFlex.cc

@@ -142,9 +142,11 @@ void NextFlex::DefineMaterials()
 {
   // Copper
   copper_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+  copper_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Air
   air_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR");
+  air_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Defining the gas xenon
   if (gas_name_ == "naturalXe")

source/geometries/NextFlexEnergyPlane.cc

@@ -143,9 +143,11 @@ void NextFlexEnergyPlane::DefineMaterials()
 
   // Copper
   copper_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+  copper_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Teflon
   teflon_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+  teflon_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Sapphire
   sapphire_mat_ = materials::Sapphire();

source/geometries/NextFlexFieldCage.cc

@@ -292,6 +292,7 @@ void NextFlexFieldCage::DefineMaterials()
 
   // Teflon
   teflon_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+  teflon_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // UV shifting material
   if (wls_mat_name_ == "NONE") {

source/geometries/NextFlexTrackingPlane.cc

@@ -200,9 +200,11 @@ void NextFlexTrackingPlane::DefineMaterials()
 
   // Copper
   copper_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+  copper_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // Teflon
   teflon_mat_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+  teflon_mat_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
   // UV shifting material
   if (wls_mat_name_ == "NONE") {

source/geometries/NextNew.cc

@@ -179,13 +179,17 @@ namespace nexus {
 				   hallA_length, hallA_length);
       G4Material *vacuum =
 	G4NistManager::Instance()->FindOrBuildMaterial("G4_Galactic");
+      vacuum->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
       lab_logic_ = new G4LogicalVolume(lab_solid, vacuum, "LAB");
       this->SetSpan(2 * hallA_length);
     } else {
       G4Box* lab_solid =
 	new G4Box("LAB", lab_size_/2., lab_size_/2., lab_size_/2.);
 
-      lab_logic_ = new G4LogicalVolume(lab_solid, G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR"), "LAB");
+      G4Material* air_mat = G4NistManager::Instance()->FindOrBuildMaterial("G4_AIR");
+      air_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+      lab_logic_ = new G4LogicalVolume(lab_solid, air_mat, "LAB");
     }
     lab_logic_->SetVisAttributes(G4VisAttributes::GetInvisible());
 

source/geometries/NextNewEnergyPlane.cc

@@ -138,10 +138,12 @@ namespace nexus {
       new G4SubtractionSolid("CARRIER_PLATE", carrier_plate_solid, axial_port_hole_solid,
 			     0, G4ThreeVector(0., 0., - carrier_plate_front_buffer_thickness_ - axial_port_thickn_));
 
+    G4Material* copper_mat = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+    copper_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4LogicalVolume* carrier_plate_logic =
       new G4LogicalVolume(carrier_plate_solid,
-			  G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu"),
+			  copper_mat,
 			  "CARRIER_PLATE");
 
 

source/geometries/NextNewFieldCage.cc

@@ -281,9 +281,11 @@ namespace nexus {
 
     // Copper for field rings
     copper_ = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+    copper_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
     // Teflon for the light tube
     teflon_ =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+    teflon_->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
     // TPB coating
     tpb_ = materials::TPB();
     tpb_->SetMaterialPropertiesTable(opticalprops::TPB());

source/geometries/NextNewIcs.cc

@@ -125,10 +125,11 @@ namespace nexus {
     ics_solid = new G4SubtractionSolid("ICS", ics_solid, up_big_nozzle_hole_solid, rotx,
 				                               G4ThreeVector( 0., up_nozzle_y_pos_, center_nozzle_z_pos_ +up_nozzle_z_pos_));
 
+    G4Material* copper_mat = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+    copper_mat->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4LogicalVolume* ics_logic =
-      new G4LogicalVolume(ics_solid,
-     						          G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu"), "ICS");
+      new G4LogicalVolume(ics_solid, copper_mat, "ICS");
 
     new G4PVPlacement(0, G4ThreeVector(0.,0.,-body_zpos_), ics_logic, "ICS", mother_logic_, false, 0, false);
 

source/geometries/NextNewKDB.cc

@@ -111,6 +111,7 @@ namespace nexus {
 
     G4Material* kapton =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON");
+    kapton->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4LogicalVolume* board_logic =
       new G4LogicalVolume(board_solid, kapton, "DICE_BOARD");
@@ -125,6 +126,7 @@ namespace nexus {
 
     G4Material* teflon =
       G4NistManager::Instance()->FindOrBuildMaterial("G4_TEFLON");
+    teflon->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
 
     G4LogicalVolume* mask_logic =
       new G4LogicalVolume(mask_solid, teflon, "DICE_MASK");

source/geometries/NextNewMiniCastle.cc

@@ -116,9 +116,11 @@ namespace nexus {
       new G4UnionSolid("MINI_CASTLE", castle_solid, z_wall_solid_sipms, 0,
      		       G4ThreeVector(0., -y_/2.+thickness_/2., z_/2.-thickness_/2.));
 
+    G4Material* lead = G4NistManager::Instance()->FindOrBuildMaterial("G4_Pb");
+    lead->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* castle_logic =
-      new G4LogicalVolume (castle_solid,
-			   G4NistManager::Instance()->FindOrBuildMaterial("G4_Pb"),"CASTLE");
+      new G4LogicalVolume (castle_solid, lead, "CASTLE");
 
     G4double y_pos = pedestal_surf_y_ + y_ -thickness_/2.;
     new G4PVPlacement(0, G4ThreeVector(0., y_pos, 0.), castle_logic,
@@ -158,8 +160,11 @@ namespace nexus {
       new G4UnionSolid("MINI_CASTLE_STEEL", steel_solid, z_steel_solid_sipms, 0,
       		       G4ThreeVector(0., -steel_y/2.+steel_thickn_/2., steel_z/2.-steel_thickn_/2.));
 
+    G4Material* steel = materials::Steel();
+    steel->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* steel_logic =
-      new G4LogicalVolume (steel_solid, materials::Steel(), "MINI_CASTLE_STEEL");
+      new G4LogicalVolume (steel_solid, steel, "MINI_CASTLE_STEEL");
 
     y_pos = pedestal_surf_y_ + y_ -thickness_ - steel_thickn_/2.;
     new G4PVPlacement(0, G4ThreeVector(0., y_pos, 0.), steel_logic,

source/geometries/NextNewPedestal.cc

@@ -10,6 +10,7 @@
 #include "MaterialsList.h"
 #include "Visibilities.h"
 
+#include <G4Material.hh>
 #include <G4GenericMessenger.hh>
 #include <G4LogicalVolume.hh>
 #include <G4PVPlacement.hh>
@@ -53,8 +54,11 @@ namespace nexus {
     G4Box* table_solid =
       new G4Box("TABLE", table_x_/2., table_y_/2., table_z_/2.);
 
+    G4Material* steel_316_Ti = materials::Steel316Ti();
+    steel_316_Ti->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* table_logic = new G4LogicalVolume(table_solid,
-						       materials::Steel316Ti(),
+						       steel_316_Ti,
 						       "TABLE");
     new G4PVPlacement(0, G4ThreeVector(0.,y_pos_,0.), table_logic, "TABLE", mother_logic_, false, 0,false);
 

source/geometries/NextNewPmtEnclosure.cc

@@ -104,9 +104,12 @@ namespace nexus{
     G4UnionSolid* enclosure_solid =
       new G4UnionSolid("ENCLOSURE_SOLID", enclosure_body, enclosure_endcap, nullptr, transl);
 
+    G4Material* copper = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu");
+    copper->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* enclosure_logic =
       new G4LogicalVolume(enclosure_solid,
-                          G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu"),
+                          copper,
                           "ENCLOSURE");
     GeometryBase::SetLogicalVolume(enclosure_logic);
 
@@ -198,9 +201,13 @@ namespace nexus{
     // Adding the PMT base
     G4Tubs* pmt_base_solid =
       new G4Tubs("PMT_BASE", 0., pmt_base_diam_/2., pmt_base_thickness_, 0.,twopi);
+
+    G4Material* kapton = G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON");
+    kapton->SetMaterialPropertiesTable(new G4MaterialPropertiesTable());
+
     G4LogicalVolume* pmt_base_logic =
       new G4LogicalVolume(pmt_base_solid,
-                          G4NistManager::Instance()->FindOrBuildMaterial("G4_KAPTON"),
+                          kapton,
                           "PMT_BASE");
     new G4PVPlacement(nullptr, G4ThreeVector(0.,0., -pmt_base_z_),
                       pmt_base_logic, "PMT_BASE", enclosure_gas_logic, false, 0, false);