Hello,
I have been trying to subtract two solid boxes, one is larger and one is smaller, I want to subtract the smaller box from the larger one. I tried to G4SubtractionSolid class to do this task but even after subtraction, I am seeing both the boxes! Not sure where did I make a mistake. Could anyone point me to where did I make a mistake by looking at the code provided below? Below is a copy-paste of my detectorconstruction.cc file. Thank you for your help in advance.
#include "DetectorConstruction.hh"
DetectorConstruction::DetectorConstruction()
{;}
DetectorConstruction::~DetectorConstruction()
{;}
G4VPhysicalVolume* DetectorConstruction::Construct()
{
// Material Information imported from NIST database
G4NistManager* nist = G4NistManager::Instance();
G4Material* Water = nist->FindOrBuildMaterial("G4_WATER");
//G4Material* Air = nist->FindOrBuildMaterial("G4_AIR");
//Polystyrene (C8H8)n (This is a material used in Cells Wells plate)
G4Material* Pstyrene = new G4Material("Polystyrene", 1.03*g/cm3, 2);
Pstyrene->AddElement(nist->FindOrBuildElement("C"), 8);
Pstyrene->AddElement(nist->FindOrBuildElement("H"), 8);
// Define a World Volume
G4VSolid* solidWorld = new G4Box("solidbox", 120*mm,120*mm,120*mm); // half length of sides
G4LogicalVolume* logicWorld = new G4LogicalVolume(solidWorld, Water, "logicWorld");
G4VPhysicalVolume* physWorld = new G4PVPlacement(0, // no rotation of daughter frame wrt mother frame
G4ThreeVector(), // position in mother volume at (0,0,0)
logicWorld, // logical volume
"physWorld", // name
0, // mother volume
false, // no boolean operation, set it to false
0, // copy number, unique arbitary number
true); // optional overlap check
// Define Cell well outer box
G4VSolid* outerbox = new G4Box("phantombox", 63.88*mm, 7.3*mm, 42.74*mm); //half length of cell well box
G4LogicalVolume* logicCell_well = new G4LogicalVolume(outerbox, Pstyrene, "logicCell_well"); // logical volume - assign properties to objexct material
G4VPhysicalVolume* fPhantomPhysOuter = new G4PVPlacement(0,
G4ThreeVector(),
logicCell_well, "fPhantomPhysOuter", // physical voloume - place the created object in the coordinate system of mother volume
logicWorld, false, 0);
// Visualization attributes -- outer box
//logicWorld->SetVisAttributes(G4VisAttributes::GetInvisible()); // world volume box will dissappear
// G4VisAttributes* outerBoxVisAtt= new G4VisAttributes(G4Colour(1, 1, 1)); // white volume
// outerBoxVisAtt-> SetVisibility(true);
// outerBoxVisAtt-> SetForceSolid(true);
// logicCell_well-> SetVisAttributes(outerBoxVisAtt);
// Define inner box
G4VSolid* innerbox = new G4Box("innerBox", 48.88*mm, 5.09*mm, 32.74*mm); // 97.76/2, 10.18/2, 65.47/2
G4LogicalVolume* logicCell_well_inner = new G4LogicalVolume(innerbox, Pstyrene, "logicCell_well_inner");
G4VPhysicalVolume* fPhantomPhysInner = new G4PVPlacement(0,
G4ThreeVector(0., 3.24*mm, 0.), //14.6 minus 10.18/2
logicCell_well_inner,
"fPhantomPhysInner",
logicWorld,
false, 0);
// Visualization attributes -- inner box
// G4VisAttributes* innerBoxVisAtt= new G4VisAttributes(G4Colour(1, 1, 0));
// innerBoxVisAtt-> SetVisibility(true);
// innerBoxVisAtt-> SetForceSolid(true);
// logicCell_well_inner-> SetVisAttributes(innerBoxVisAtt);
// Subtract inner box from outer box
G4SubtractionSolid* solidShape = new G4SubtractionSolid("outer-inner",
outerbox,
innerbox,
0, // no rotation
G4ThreeVector(0., 0., 0.)); // no translation
G4LogicalVolume* logicCellWell = new G4LogicalVolume(solidShape, Pstyrene, "logicCellWell", 0, 0, 0);
G4VPhysicalVolume* physCellWell = new G4PVPlacement(0,
G4ThreeVector(),
logicCellWell,
"physCellWell",
logicWorld, // its mother volume
false, // Boolean operation
0, true); // check overlaps
// Visualization attributes -- subtracted volume
G4VisAttributes* subtractedBoxVisAtt= new G4VisAttributes(G4Colour(0.4, 0.4, 0.4));
subtractedBoxVisAtt-> SetVisibility(true);
subtractedBoxVisAtt-> SetForceSolid(true);
logicCellWell-> SetVisAttributes(subtractedBoxVisAtt);
return physWorld;
}