30 #include "B3DetectorConstruction.hh"
91 G4double cosdPhi = std::cos(half_dPhi);
92 G4double tandPhi = std::tan(half_dPhi);
94 G4double ring_R1 = 0.5*cryst_dY/tandPhi;
95 G4double ring_R2 = (ring_R1+cryst_dZ)/cosdPhi;
97 G4double detector_dZ = nb_rings*cryst_dX;
106 G4double world_sizeXY = 2.4*ring_R2;
107 G4double world_sizeZ = 1.2*detector_dZ;
111 0.5*world_sizeXY, 0.5*world_sizeXY, 0.5*world_sizeZ);
132 new G4Tubs(
"Ring", ring_R1, ring_R2, 0.5*cryst_dX, 0.,
twopi);
144 G4Box* solidCryst =
new G4Box(
"crystal", dX/2, dY/2, cryst_dZ/2);
153 for (
G4int icrys = 0; icrys < nb_cryst ; icrys++) {
175 new G4Tubs(
"Detector", ring_R1, ring_R2, 0.5*detector_dZ, 0.,
twopi);
185 G4double OG = -0.5*(detector_dZ + cryst_dX);
186 for (
G4int iring = 0; iring < nb_rings ; iring++) {
218 new G4Tubs(
"Patient", 0., patient_radius, 0.5*patient_dZ, 0.,
twopi);