13 Double_t waferThickness = 200 *
eic::um;
16 Double_t waferWidth = 250 *
eic::mm;
17 Double_t waferSpacing = 100 *
eic::mm;
20 Double_t containerVolumeLength = (waferNum-1) * waferSpacing + waferThickness + (10 *
eic::mm);
21 Double_t containerVolumeWidth = waferWidth + (10 *
eic::mm);
24 TGeoBBox *container =
new TGeoBBox(
"ContainerVolume",
25 containerVolumeWidth/2,
26 containerVolumeWidth/2,
27 containerVolumeLength/2);
29 TGeoVolume *vcontainer =
new TGeoVolume(
"ContainerVolume", container, fst->
GetMedium(
"air"));
32 TGeoBBox *wafer =
new TGeoBBox(
"SiliconWafer", waferWidth/2, waferWidth/2, waferThickness/2);
33 TGeoVolume *vwafer =
new TGeoVolume(
"SiliconWafer", wafer, fst->
GetMedium(
"silicon"));
48 for(
unsigned wf=0; wf<waferNum; wf++) {
49 double offset = (wf - (waferNum-1)/2.)*waferSpacing;
51 UInt_t geant[1] = {wf}, group = 0, logical[3] = {0, 0, wf};
54 cout <<
"Failed to set mapping table entry!" << endl;
59 vcontainer->AddNode(vwafer, wf,
new TGeoCombiTrans(0.0, 0.0, offset,
new TGeoRotation()));
64 fst->
GetTopVolume()->AddNode(vcontainer, 0,
new TGeoCombiTrans(0.0, 0.0, 0.0,
new TGeoRotation()));