44{
45
46
47 G4LogicalVolume* logicalVolume =
48 TG4GeometryServices::Instance()->FindLogicalVolume("LAYE");
49 if (!logicalVolume) {
50 G4cerr << "Logical volume LAYE not found" << G4endl;
51 return;
52 }
53
54 G4Region* testRegion = new G4Region("Test_region");
55 testRegion->AddRootLogicalVolume(logicalVolume);
56 testRegion->SetProductionCuts(new G4ProductionCuts());
57
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);
63}