47 G4LogicalVolume* logicalVolume =
48 TG4GeometryServices::Instance()->FindLogicalVolume(
"LAYE");
50 G4cerr <<
"Logical volume LAYE not found" << G4endl;
54 G4Region* testRegion =
new G4Region(
"Test_region");
55 testRegion->AddRootLogicalVolume(logicalVolume);
56 testRegion->SetProductionCuts(
new G4ProductionCuts());
58 std::vector<double> cuts;
59 cuts.push_back(1.0 * cm);
60 cuts.push_back(1.0 * cm);
61 cuts.push_back(1.0 * cm);
62 testRegion->GetProductionCuts()->SetProductionCuts(cuts);