How to overlap volumes?

Geometry, Fields and Transportation
1

I have a partially coaxial HPGe detector. Where I have a cylindrical germanium crystal with a cylindrical hole at the center, like this:

Screenshot_2021-04-03-17-36-00-67
Screenshot_2021-04-03-17-36-00-67869×587 17.6 KB

When I try to draw a cylinder of HPGe (inner radius = 0) and then draw another cylinder of vaccum (height is less than HPGe cylinder height), it shows me the warning like this for every:

In file included from /home/subhrod/G4_project/Examples_G4/zzHPGe_design/src/hpgeDetectorConstruction.cc:1:
/home/subhrod/G4_project/Examples_G4/zzHPGe_design/include/hpgeDetectorConstruction.hh:33:24: note: shadowed declaration is here
   33 |     G4VPhysicalVolume* physiWorld;
      |                        ^~~~~~~~~~

Will the vaccum be placed at there, even if it’s giving me the warning? And how to remove this warning?

I have set checkOverlap to true and false both but it is still showing.

2

The error you are showing only says that you are shadowing the physiWorld definition.

In my opinion, the best solution is to subtract two volumes from each other. This is what the G4SubtractionSolid class is for.

1 Like
3

Hi @MKasztelan,
It was on my mind. I have tried to do the cutting of the material for ease of coding using G4SubtractionSolid. I have written the code like:

G4double hpgeRi = 0. * cm;
    G4double hpgeRo = 3.05 * cm;
    G4double hpgeHz = 1.6 * cm;
    G4double hpgeAo = 0. * deg;
    G4double hpgeAs = 360. * deg;

    G4VSolid *hpgeTube = new G4Tubs("hpgeTube",
                                    hpgeRi,
                                    hpgeRo,
                                    hpgeHz,
                                    hpgeAo,
                                    hpgeAs);

    ///////
    //Bore Hole========
    G4double boreRi = 0. * mm;
    G4double boreRo = 7.072 * mm;
    G4double boreHz = 6.337 * mm;
    G4double boreAo = 0. * deg;
    G4double boreAs = 360. * deg;

    G4VSolid *boreTube = new G4Tubs("boreTube",
                                    boreRi,
                                    boreRo,
                                    boreHz,
                                    boreAo,
                                    boreAs);
    //position of the hole wrt hpge....
    G4double boreposx = 0. * cm;
    G4double boreposy = 0. * cm;
    G4double boreposz = -(hpgeHz - boreHz);
    G4VSolid* hpgeWithBore = new G4SubtractionSolid("hpgeWithBore",                              //name of the solid
                                                    hpgeTube,                                    //parent volume
                                                    boreTube,                                    //substractor volume
                                                    0,                                           // spin of subtractor w.r.t parent
                                                    G4ThreeVector(boreposx, boreposy, boreposz)); //position w.r.t. parent

But now I am getting the error:

$ make
Scanning dependencies of target hpge
[ 12%] Building CXX object CMakeFiles/hpge.dir/hpge_main.cc.o
[ 25%] Building CXX object CMakeFiles/hpge.dir/src/hpgeActionInitialization.cc.o
[ 37%] Building CXX object CMakeFiles/hpge.dir/src/hpgeDetectorConstruction.cc.o
/home/subhrod/G4_project/zzHPGe_design/src/hpgeDetectorConstruction.cc: In member function ‘virtual G4VPhysicalVolume* hpgeDetectorConstruction::Construct()’:
/home/subhrod/G4_project/zzHPGe_design/src/hpgeDetectorConstruction.cc:217:30: error: ‘G4SubtractionSolid’ was not declared in this scope
  217 |     G4VSolid* hpgeWithBore = G4SubtractionSolid("hpgeWithBore",                              //name of the solid
      |                              ^~~~~~~~~~~~~~~~~~
make[2]: *** [CMakeFiles/hpge.dir/build.make:89: CMakeFiles/hpge.dir/src/hpgeDetectorConstruction.cc.o] Error 1
make[1]: *** [CMakeFiles/Makefile2:104: CMakeFiles/hpge.dir/all] Error 2
make: *** [Makefile:130: all] Error 2

Can you please advise me on how to solve this issue?
Thank you.

4

#include "G4SubtractionSolid.hh". C++ is (or was :frowning: ) a strongly typed language. You need explicitly declare a class before attempting to use it.

2 Likes
5

Thank you @mkelsey and @MKasztelan
I thought G4SubtractionSolid is already included in G4Solid.hh header file.

1 Like
6

Each class has its own .hh file. G4VSolid.hh is the abstract base class for all of the concrete solids.

2 Likes
7

Yes sir,
Now I understand it fully.
Thank you so much for your support.

1 Like