Belle II Software development
Public Member Functions | Private Member Functions | Private Attributes | List of all members
GeoECLCreator Class Reference

The GeoECLCreator class. More...

#include <GeoECLCreator.h>

Inheritance diagram for GeoECLCreator:
CreatorBase

Public Member Functions

 GeoECLCreator ()
 Constructor of the GeoECLCreator class.
 
 ~GeoECLCreator ()
 The destructor of the GeoECLCreator class.
 
virtual void create (const GearDir &content, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
 Function to actually create the geometry, has to be overridden by derived classes.
 
virtual void createFromDB (const std::string &name, G4LogicalVolume &topVolume, geometry::GeometryTypes type) override
 Function to create the geometry from the Database.
 
virtual void createPayloads (const GearDir &content, const IntervalOfValidity &iov) override
 Function to create the geometry database.
 
 BELLE2_DEFINE_EXCEPTION (DBNotImplemented, "Cannot create geometry from Database.")
 Exception that will be thrown in createFromDB if member is not yet implemented by creator.
 

Private Member Functions

void barrel (G4LogicalVolume &)
 Make the ECL barrel and then place elements inside it.
 
void backward (G4LogicalVolume &)
 Place elements inside the backward endcap.
 
void forward (G4LogicalVolume &)
 Place elements inside the forward endcap.
 
G4LogicalVolume * wrapped_crystal (const shape_t *s, const std::string &endcap, double wrapthickness)
 Wrapped crystal.
 
void defineVisAttributes ()
 Define visual attributes.
 
const G4VisAttributes * att (const std::string &n) const
 Define visual attributes.
 
G4LogicalVolume * get_preamp () const
 Get Logical volume of preamplifier.
 
double get_pa_box_height () const
 Getter for preamplifier box height (hard-coded to be 2)
 

Private Attributes

const ECLCrystalsShapeAndPosition * m_sap
 pointer to a storage with crystal shapes and positions
 
Simulation::SensitiveDetectorBase * m_sensitive
 Sensitive detector.
 
Simulation::SensitiveDetectorBase * m_sensediode
 Sensitive diode.
 
std::map< std::string, G4VisAttributes * > m_atts
 Vector of background-Sensitive detectors.
 
int m_overlap
 overlap
 

Detailed Description

The GeoECLCreator class.

The creator for the ECL geometry of the Belle II detector.

Definition at line 34 of file GeoECLCreator.h.

Constructor & Destructor Documentation

◆ GeoECLCreator()

GeoECLCreator ( )

Constructor of the GeoECLCreator class.

Definition at line 43 of file GeoECLCreator.cc.

43 : m_sap(0), m_overlap(0)
44{
45 m_sensediode = nullptr;
46 m_sensitive = new SensitiveDetector("ECLSensitiveDetector", (2 * 24)*CLHEP::eV, 10 * CLHEP::MeV);
47 G4SDManager::GetSDMpointer()->AddNewDetector(m_sensitive);
48 defineVisAttributes();
49}
Belle2::ECL::GeoECLCreator::defineVisAttributes
void defineVisAttributes()
Define visual attributes.
Definition GeoECLCreator.cc:119
Belle2::ECL::GeoECLCreator::m_sensediode
Simulation::SensitiveDetectorBase * m_sensediode
Sensitive diode.
Definition GeoECLCreator.h:96
Belle2::ECL::GeoECLCreator::m_sap
const ECLCrystalsShapeAndPosition * m_sap
pointer to a storage with crystal shapes and positions
Definition GeoECLCreator.h:91
Belle2::ECL::GeoECLCreator::m_sensitive
Simulation::SensitiveDetectorBase * m_sensitive
Sensitive detector.
Definition GeoECLCreator.h:94
Belle2::PXD::SensitiveDetector
VXD::SensitiveDetector< PXDSimHit, PXDTrueHit > SensitiveDetector
The PXD Sensitive Detector class.
Definition SensitiveDetector.h:23

◆ ~GeoECLCreator()

~GeoECLCreator ( )

The destructor of the GeoECLCreator class.

Definition at line 52 of file GeoECLCreator.cc.

53{
54 for (auto a : m_atts) delete a.second;
55}
Belle2::ECL::GeoECLCreator::m_atts
std::map< std::string, G4VisAttributes * > m_atts
Vector of background-Sensitive detectors.
Definition GeoECLCreator.h:98

Member Function Documentation

◆ att()

const G4VisAttributes * att ( const std::string & n) const
private

Define visual attributes.

Parameters
nAttribute name

Definition at line 141 of file GeoECLCreator.cc.

142{
143 auto p = m_atts.find(n);
144 assert(p != m_atts.end());
145 return p->second;
146}

◆ backward()

void backward ( G4LogicalVolume & _top)
private

Place elements inside the backward endcap.

Definition at line 42 of file backward.cc.

43{
44 G4LogicalVolume* top = &_top;
45
46 const bool sec = 0;
47 const double phi0 = 0;
48 const double dphi = sec ? M_PI / 16 : 2 * M_PI;
49
50 const bool b_inner_support_ring = true;
51 const bool b_outer_support_ring = true;
52 const bool b_support_wall = true;
53 const bool b_ribs = true;
54 const bool b_septum_wall = true;
55 const bool b_crystals = true;
56 const bool b_preamplifier = true;
57 const bool b_support_leg = true;
58 const int overlap = m_overlap;
59
60 int npoints = 1000 * 1000;
61
62 // vector<cplacement_t> bp = load_placements("/ecl/data/crystal_placement_backward.dat");
63 vector<cplacement_t> bp = load_placements(m_sap, ECLParts::backward);
64 vector<cplacement_t>::iterator fp = find_if(bp.begin(), bp.end(), [](const cplacement_t& p) {
65 const int ECL_backward_part = 1000; // lookup part
66 return p.nshape == ECL_backward_part;
67 });
68 // global transformation before placing the whole backward part in the top logical volume
69 G4Transform3D gTrans = (fp == bp.end()) ? G4Translate3D(0, 0, -1020) : get_transform(*fp);
70 if (fp != bp.end()) bp.erase(fp); // now not needed
71
72 if (b_inner_support_ring) {
73 zr_t vc1[] = {{0., 452.3 + 3}, {0., 452.3}, {3., 474.9 - 20 / cosd(27.81)}, {434., 702.27 - 20 / cosd(27.81)}, {434., 702.27}, {3., 474.9}, {3., 452.3 + 3}};
74 std::vector<zr_t> contour1(vc1, vc1 + sizeof(vc1) / sizeof(zr_t));
75 G4VSolid* part1solid = new BelleLathe("part1solid", phi0, dphi, contour1);
76 G4LogicalVolume* part1logical = new G4LogicalVolume(part1solid, Materials::get("SUS304"), "part1logical", 0, 0, 0);
77 part1logical->SetVisAttributes(att("iron"));
78 auto pv = new G4PVPlacement(gTrans * G4RotateY3D(M_PI), part1logical, "ECL_BackwardSupport_part1physical", top, false, 0, 0);
79 if (overlap) pv->CheckOverlaps(npoints);
80 }
81
82 if (b_support_wall) {
83 // solving equation to get L : 3+L*cosd(52.90)+1.6*cosd(52.90+90) = 435 - 202.67
84 double L = (435 - 202.67 - 3 - 1.6 * cosd(52.90 + 90)) / cosd(52.90);
85 zr_t vc23[] = {{0, 452.3 + 3}, {3, 452.3 + 3}, {3, 1190.2}, {3 + L * cosd(52.90), 1190.2 + L * sind(52.90)},
86 {3 + L * cosd(52.90) + 1.6 * cosd(52.90 + 90), 1190.2 + L * sind(52.90) + 1.6 * sind(52.90 + 90)}, {3 + 1.6 * cosd(52.90 + 90), 1190.2 + 1.6 * sind(52.90 + 90)}, {0, 1190.2}
87 };
88 std::vector<zr_t> contour23(vc23, vc23 + sizeof(vc23) / sizeof(zr_t));
89 G4VSolid* part23solid = new BelleLathe("part23solid", phi0, dphi, contour23);
90 G4LogicalVolume* part23logical = new G4LogicalVolume(part23solid, Materials::get("A5052"), "part23logical", 0, 0, 0);
91 part23logical->SetVisAttributes(att("alum"));
92 auto pv = new G4PVPlacement(gTrans * G4RotateY3D(M_PI), part23logical, "ECL_BackwardSupport_part23physical", top, false, 0, 0);
93 if (overlap) pv->CheckOverlaps(npoints);
94 }
95
96 if (b_outer_support_ring) {
97 zr_t vc4[] = {{434 - 214.8, 1496 - 20}, {434, 1496 - 20}, {434, 1496 - 5}, {434 + 5, 1496 - 5}, {434 + 5, 1496}, {434 - 199.66, 1496}};
98 std::vector<zr_t> contour4(vc4, vc4 + sizeof(vc4) / sizeof(zr_t));
99 G4VSolid* part4solid = new BelleLathe("part4solid", phi0, dphi, contour4);
100 G4LogicalVolume* part4logical = new G4LogicalVolume(part4solid, Materials::get("SUS304"), "part4logical", 0, 0, 0);
101 part4logical->SetVisAttributes(att("iron"));
102 auto pv = new G4PVPlacement(gTrans * G4RotateY3D(M_PI), part4logical, "ECL_BackwardSupport_part4physical", top, false, 0, 0);
103 if (overlap) pv->CheckOverlaps(npoints);
104 }
105
106 zr_t cont_array_in[] = {{3., 474.9}, {434., 702.27}, {434, 1496 - 20}, {434 - 214.8, 1496 - 20}, {3, 1190.2}};
107 std::vector<zr_t> contour_in(cont_array_in, cont_array_in + sizeof(cont_array_in) / sizeof(zr_t));
108 G4VSolid* innervolume_solid = new BelleLathe("innervolume_solid", 0, 2 * M_PI, contour_in);
109 G4LogicalVolume* innervolume_logical = new G4LogicalVolume(innervolume_solid, Materials::get("G4_AIR"),
110 "innervolume_logical", 0, 0, 0);
111 innervolume_logical->SetVisAttributes(att("air"));
112
113 // Set up region for production cuts
114 G4Region* aRegion = new G4Region("ECLBackwardEnvelope");
115 innervolume_logical->SetRegion(aRegion);
116 aRegion->AddRootLogicalVolume(innervolume_logical);
117
118 auto gpvbp = new G4PVPlacement(gTrans * G4RotateY3D(M_PI),
119 innervolume_logical, "ECLBackwardPhysical",
120 top, false, 0, 0);
121 if (overlap) gpvbp->CheckOverlaps(npoints);
122
123 G4VSolid* innervolumesector_solid = new BelleLathe("innervolumesector_solid", -M_PI / 8, M_PI / 4, contour_in);
124 G4LogicalVolume* innervolumesector_logical = new G4LogicalVolume(innervolumesector_solid, Materials::get("G4_AIR"),
125 "innervolumesector_logical", 0, 0, 0);
126 innervolumesector_logical->SetVisAttributes(att("air"));
127 new G4PVReplica("ECLBackwardSectorPhysical", innervolumesector_logical, innervolume_logical, kPhi, 8, M_PI / 4, 0);
128
129 if (b_ribs) {
130 double H = 60, W = 20;
131 double X0 = 702.27 + 0.001, X1 = 1496 - 20;
132 G4TwoVector r0o(X1, 0), r1o(X1 * sqrt(1 - pow(W / X1, 2)), W);
133 double beta = asin(W / X0);
134 G4TwoVector r0i(X0 / cos(beta / 2), 0), r1i(X0 * cos(beta / 2) - tan(beta / 2) * (W - X0 * sin(beta / 2)), W);
135 double dxymzp = (r0o - r0i).x(), dxypzp = (r1o - r1i).x();
136 double theta = atan(tand(27.81) / 2);
137 double dxymzm = dxymzp + tand(27.81) * H, dxypzm = dxypzp + tand(27.81) * H;
138
139 G4TwoVector m0 = (r0i + r0o) * 0.5, m1 = (r1i + r1o) * 0.5, dm = m1 - m0;
140 double alpha = atan(dm.x() / dm.y());
141
142 G4VSolid* solid6_p1 = new G4Trap("solid6_p1", H / 2, theta, 0, W / 2, dxymzm / 2, dxypzm / 2, alpha, W / 2, dxymzp / 2, dxypzp / 2,
143 alpha);
144 G4LogicalVolume* lsolid6_p1 = new G4LogicalVolume(solid6_p1, Materials::get("SUS304"), "lsolid6", 0, 0, 0);
145 lsolid6_p1->SetVisAttributes(att("asolid6"));
146 G4Transform3D tsolid6_p1(G4Translate3D(X0 * cos(beta / 2) + (dxymzp / 2 + dxypzp / 2) / 2 - tan(theta)*H / 2, W / 2, 434 - H / 2));
147 auto pv61 = new G4PVPlacement(G4RotateZ3D(-M_PI / 8)*tsolid6_p1, lsolid6_p1, "psolid6_p1", innervolumesector_logical, false, 0, 0);
148 if (overlap) pv61->CheckOverlaps(npoints);
149 auto pv62 = new G4PVPlacement(G4RotateZ3D(0)*tsolid6_p1, lsolid6_p1, "psolid6_p2", innervolumesector_logical, false, 0, 0);
150 if (overlap) pv62->CheckOverlaps(npoints);
151
152 H = 40;
153 dxymzm = dxymzp + tand(27.81) * H, dxypzm = dxypzp + tand(27.81) * H;
154 G4VSolid* solid6_p2 = new G4Trap("solid6_p2", H / 2, theta, 0, W / 2, dxypzm / 2, dxymzm / 2, -alpha, W / 2, dxypzp / 2, dxymzp / 2,
155 -alpha);
156 G4LogicalVolume* lsolid6_p2 = new G4LogicalVolume(solid6_p2, Materials::get("SUS304"), "lsolid6", 0, 0, 0);
157 lsolid6_p2->SetVisAttributes(att("asolid6"));
158 G4Transform3D tsolid6_p2(G4Translate3D(X0 * cos(beta / 2) + (dxymzp / 2 + dxypzp / 2) / 2 - tan(theta)*H / 2, -W / 2, 434 - H / 2));
159 auto pv63 = new G4PVPlacement(G4RotateZ3D(0)*tsolid6_p2, lsolid6_p2, "psolid6_p3", innervolumesector_logical, false, 0, 0);
160 if (overlap) pv63->CheckOverlaps(npoints);
161 auto pv64 = new G4PVPlacement(G4RotateZ3D(M_PI / 8)*tsolid6_p2, lsolid6_p2, "psolid6_p4", innervolumesector_logical, false, 0, 0);
162 if (overlap) pv64->CheckOverlaps(npoints);
163
164 G4VSolid* solid7_p8 = new G4Box("solid7_p8", 171. / 2, (140. - 40) / 2 / 2, 40. / 2);
165 G4LogicalVolume* lsolid7 = new G4LogicalVolume(solid7_p8, Materials::get("SUS304"), "lsolid7", 0, 0, 0);
166 lsolid7->SetVisAttributes(att("asolid7"));
167 double dx = sqrt(X1 * X1 - 70 * 70) - 171. / 2;
168 G4Transform3D tsolid7_p1(G4Translate3D(dx, -20 - 25, 434 - 40. / 2));
169 auto pv71 = new G4PVPlacement(tsolid7_p1, lsolid7, "psolid7_p1", innervolumesector_logical, false, 0, 0);
170 if (overlap) pv71->CheckOverlaps(npoints);
171 G4Transform3D tsolid7_p2(G4Translate3D(dx, 20 + 25, 434 - 40. / 2));
172 auto pv72 = new G4PVPlacement(tsolid7_p2, lsolid7, "psolid7_p2", innervolumesector_logical, false, 0, 0);
173 if (overlap) pv72->CheckOverlaps(npoints);
174
175 double L = X1 - (X0 - tand(27.81) * 40) - 10;
176 G4VSolid* solid13 = new G4Box("solid13", L / 2, 5. / 2, 18. / 2);
177 G4LogicalVolume* lsolid13 = new G4LogicalVolume(solid13, Materials::get("SUS304"), "lsolid13", 0, 0, 0);
178 lsolid13->SetVisAttributes(att("asolid13"));
179 G4Transform3D tsolid13(G4TranslateZ3D(434 - 60 + 18. / 2)*G4TranslateY3D(-5. / 2 - 0.5 / 2)*G4TranslateX3D(X0 - tand(
180 27.81) * 40 + L / 2 + 5));
181 auto pv131 = new G4PVPlacement(tsolid13, lsolid13, "psolid13_p1", innervolumesector_logical, false, 0, 0);
182 if (overlap) pv131->CheckOverlaps(npoints);
183 auto pv132 = new G4PVPlacement(G4RotateZ3D(M_PI / 8)*tsolid13, lsolid13, "psolid13_p2", innervolumesector_logical, false, 0, 0);
184 if (overlap) pv132->CheckOverlaps(npoints);
185 }
186
187
188 double zsep = 135;
189 if (b_septum_wall) {
190 double d = 5;
191 Point_t vin[] = {{434. - zsep, 702.27 - tand(27.81)* zsep}, {434. - 60, 702.27 - tand(27.81) * 60}, {434. - 60, 1496 - 20 - d}, {434. - zsep, 1496 - 20 - d}};
192 const int n = sizeof(vin) / sizeof(Point_t);
193 Point_t c = centerofgravity(vin, vin + n);
194 G4ThreeVector contour_swall[n * 2];
195 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2);
196 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2);
197
198 G4VSolid* septumwall_solid = new BelleCrystal("septumwall_solid", n, contour_swall);
199
200 G4LogicalVolume* septumwall_logical = new G4LogicalVolume(septumwall_solid, Materials::get("A5052"),
201 "septumwall_logical", 0, 0, 0);
202 septumwall_logical->SetVisAttributes(att("alum2"));
203 auto pv = new G4PVPlacement(G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, 0), septumwall_logical,
204 "septumwall_physical", innervolumesector_logical, false, 0, 0);
205 if (overlap) pv->CheckOverlaps(npoints);
206
207 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2 / 2);
208 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2 / 2);
209
210 G4VSolid* septumwall2_solid = new BelleCrystal("septumwall2_solid", n, contour_swall);
211
212 G4LogicalVolume* septumwall2_logical = new G4LogicalVolume(septumwall2_solid, Materials::get("A5052"),
213 "septumwall2_logical", 0, 0, 0);
214 septumwall2_logical->SetVisAttributes(att("alum2"));
215 auto pv0 = new G4PVPlacement(G4RotateZ3D(-M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y,
216 0.5 / 2 / 2),
217 septumwall2_logical, "septumwall2_physical", innervolumesector_logical, false, 0, 0);
218 if (overlap) pv0->CheckOverlaps(npoints);
219 auto pv1 = new G4PVPlacement(G4RotateZ3D(M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y,
220 -0.5 / 2 / 2),
221 septumwall2_logical, "septumwall2_physical", innervolumesector_logical, false, 1, 0);
222 if (overlap) pv1->CheckOverlaps(npoints);
223 }
224
225 zr_t vcr[] = {{3., 474.9}, {434. - zsep, 702.27 - tand(27.81)* zsep}, {434 - zsep, 1496 - 20}, {434 - 214.8, 1496 - 20}, {3, 1190.2}};
226 std::vector<zr_t> ccr(vcr, vcr + sizeof(vcr) / sizeof(zr_t));
227 G4VSolid* crystalvolume_solid = new BelleLathe("crystalvolume_solid", 0, M_PI / 8, ccr);
228 G4LogicalVolume* crystalvolume_logical = new G4LogicalVolume(crystalvolume_solid, Materials::get("G4_AIR"),
229 "crystalvolume_logical", 0, 0, 0);
230 crystalvolume_logical->SetVisAttributes(att("air"));
231 auto gpv0 = new G4PVPlacement(G4RotateZ3D(-M_PI / 8), crystalvolume_logical, "ECLBackwardCrystalSectorPhysical_0",
232 innervolumesector_logical,
233 false, 0, 0);
234 if (overlap) gpv0->CheckOverlaps(npoints);
235 auto gpv1 = new G4PVPlacement(G4RotateZ3D(0), crystalvolume_logical, "ECLBackwardCrystalSectorPhysical_1",
236 innervolumesector_logical, false, 1,
237 0);
238 if (overlap) gpv1->CheckOverlaps(npoints);
239
240 if (b_septum_wall) {
241 double d = 5, dr = 0.001;
242 Point_t vin[] = {{3., 474.9}, {434. - zsep, 702.27 - tand(27.81)* zsep}, {434 - zsep, 1496 - 20 - d - dr}, {434 - 214.8 - d / tand(52.90), 1496 - 20 - d - dr}, {3, 1190.2 - dr}};
243 const int n = sizeof(vin) / sizeof(Point_t);
244 Point_t c = centerofgravity(vin, vin + n);
245 G4ThreeVector contour_swall[n * 2];
246
247 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2 / 2);
248 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2 / 2);
249
250 G4VSolid* septumwall3_solid = new BelleCrystal("septumwall3_solid", n, contour_swall);
251
252 G4LogicalVolume* septumwall3_logical = new G4LogicalVolume(septumwall3_solid, Materials::get("A5052"),
253 "septumwall3_logical", 0, 0, 0);
254 septumwall3_logical->SetVisAttributes(att("alum2"));
255 auto pv0 = new G4PVPlacement(G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, 0.5 / 2 / 2),
256 septumwall3_logical,
257 "septumwall3_physical_0", crystalvolume_logical, false, 0, overlap);
258 if (overlap) pv0->CheckOverlaps(npoints);
259 auto pv1 = new G4PVPlacement(G4RotateZ3D(M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y,
260 -0.5 / 2 / 2),
261 septumwall3_logical, "septumwall3_physical_1", crystalvolume_logical, false, 1, overlap);
262 if (overlap) pv1->CheckOverlaps(npoints);
263 }
264
265 if (b_crystals) {
266 // vector<shape_t*> cryst = load_shapes("/ecl/data/crystal_shape_backward.dat");
267 vector<shape_t*> cryst = load_shapes(m_sap, ECLParts::backward);
268 vector<G4LogicalVolume*> wrapped_crystals;
269 for (auto it = cryst.begin(); it != cryst.end(); it++) {
270 shape_t* s = *it;
271 wrapped_crystals.push_back(wrapped_crystal(s, "backward", 0.20 - 0.02));
272 }
273 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
274 const cplacement_t& t = *it;
275 auto s = find_if(cryst.begin(), cryst.end(), [&t](const shape_t* shape) {return shape->nshape == t.nshape;});
276 if (s == cryst.end()) continue;
277
278 G4Transform3D twc = G4Translate3D(0, 0, 3) * get_transform(t);
279 int indx = it - bp.begin();
280 auto pv = new G4PVPlacement(twc, wrapped_crystals[s - cryst.begin()], suf("ECLBackwardWrappedCrystal_Physical", indx),
281 crystalvolume_logical,
282 false, (ECLElementNumbers::c_NCrystalsForwardBarrel) / 16 + indx, 0);
283 if (overlap)pv->CheckOverlaps(npoints);
284 }
285
286 for (shape_t* shape : cryst) {
287 delete shape;
288 }
289 }
290
291 if (b_preamplifier) {
292 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
293 G4Transform3D twc = G4Translate3D(0, 0, 3) * get_transform(*it);
294 int indx = it - bp.begin();
295 auto pv = new G4PVPlacement(twc * G4TranslateZ3D(300 / 2 + 0.20 + get_pa_box_height() / 2)*G4RotateZ3D(-M_PI / 2), get_preamp(),
296 suf("phys_backward_preamplifier", indx), crystalvolume_logical, false, (ECLElementNumbers::c_NCrystalsForwardBarrel) / 16 + indx,
297 0);
298 if (overlap)pv->CheckOverlaps(npoints);
299 }
300 }
301
302 if (b_support_leg) {
303 const G4VisAttributes* batt = att("iron");
304
305 G4VSolid* s1 = new G4Box("leg_p1", 130. / 2, 185. / 2, (40. - 5) / 2);
306 G4LogicalVolume* l1 = new G4LogicalVolume(s1, Materials::get("SUS304"), "l1", 0, 0, 0);
307 G4Transform3D t1 = G4Translate3D(0, 185. / 2, (40. - 5) / 2);
308 l1->SetVisAttributes(batt);
309
310 Point_t v3[] = {{ -212. / 2, -135. / 2}, {212. / 2 - 30, -135. / 2}, {212. / 2, -135. / 2 + 30}, {212. / 2, 135. / 2}, { -212. / 2, 135. / 2}};
311 const int n3 = sizeof(v3) / sizeof(Point_t);
312 G4ThreeVector c3[n3 * 2];
313
314 for (int i = 0; i < n3; i++) c3[i + 0] = G4ThreeVector(v3[i].x, v3[i].y, -60. / 2);
315 for (int i = 0; i < n3; i++) c3[i + n3] = G4ThreeVector(v3[i].x, v3[i].y, 60. / 2);
316
317 G4VSolid* s3 = new BelleCrystal("leg_p3", n3, c3);
318 G4LogicalVolume* l3 = new G4LogicalVolume(s3, Materials::get("SUS304"), "l3", 0, 0, 0);
319 G4Transform3D t3 = G4Translate3D(0, 135. / 2 + 35, 40. - 5. + 212. / 2) * G4RotateY3D(-M_PI / 2);
320 l3->SetVisAttributes(batt);
321
322 G4VSolid* s6 = new G4Box("leg_p6", 140. / 2, 189. / 2, 160. / 2);
323 G4LogicalVolume* l6 = new G4LogicalVolume(s6, Materials::get("G4_AIR"), "l6", 0, 0, 0);
324 G4Transform3D t6 = G4Translate3D(0, 170. + 189. / 2, 57. + 35. + 160. / 2);
325 l6->SetVisAttributes(att("air"));
326
327 G4VSolid* s6a = new G4Box("leg_p6a", 140. / 2, (189. - 45.) / 2, 160. / 2);
328 G4LogicalVolume* l6a = new G4LogicalVolume(s6a, Materials::get("SUS304"), "l6a", 0, 0, 0);
329 l6a->SetVisAttributes(batt);
330 new G4PVPlacement(G4TranslateY3D(-45. / 2), l6a, "l6a_physical", l6, false, 0, overlap);
331
332 G4VSolid* s6b = new G4Box("leg_p6b", 60. / 2, 45. / 2, 160. / 2);
333 G4LogicalVolume* l6b = new G4LogicalVolume(s6b, Materials::get("SUS304"), "l6b", 0, 0, 0);
334 l6b->SetVisAttributes(batt);
335 double dy = 189. / 2 - 45 + 45. / 2;
336 new G4PVPlacement(G4TranslateY3D(dy), l6b, "l6b_physical", l6, false, 0, overlap);
337
338 G4VSolid* s6c = new G4Box("leg_p6c", 40. / 2, 45. / 2, 22.5 / 2);
339 G4LogicalVolume* l6c = new G4LogicalVolume(s6c, Materials::get("SUS304"), "l6c", 0, 0, 0);
340 l6c->SetVisAttributes(batt);
341 new G4PVPlacement(G4Translate3D(30 + 20, dy, 20 + 22.5 / 2), l6c, "l6c_physical", l6, false, 0, overlap);
342 new G4PVPlacement(G4Translate3D(30 + 20, dy, -20 - 22.5 / 2), l6c, "l6c_physical", l6, false, 1, overlap);
343 new G4PVPlacement(G4Translate3D(-30 - 20, dy, 20 + 22.5 / 2), l6c, "l6c_physical", l6, false, 2, overlap);
344 new G4PVPlacement(G4Translate3D(-30 - 20, dy, -20 - 22.5 / 2), l6c, "l6c_physical", l6, false, 3, overlap);
345
346 G4AssemblyVolume* support_leg = new G4AssemblyVolume();
347
348 support_leg->AddPlacedVolume(l1, t1);
349 // support_leg->AddPlacedVolume(l2,t2);
350 support_leg->AddPlacedVolume(l3, t3);
351 // support_leg->AddPlacedVolume(l4,t4);
352 // support_leg->AddPlacedVolume(l5,t5);
353 support_leg->AddPlacedVolume(l6, t6);
354
355 G4VSolid* s_all = new G4Box("leg_all", 140. / 2, 359. / 2, (257. - 5.) / 2);
356 G4LogicalVolume* l_all = new G4LogicalVolume(s_all, Materials::get("G4_AIR"), "l_all", 0, 0, 0);
357 l_all->SetVisAttributes(att("air"));
358 G4Transform3D tp = G4Translate3D(0, -359. / 2, -(257. - 5.) / 2);
359 support_leg->MakeImprint(l_all, tp, 0, overlap);
360
361
362 for (int i = 0; i < 8; i++)
363 new G4PVPlacement(gTrans * G4RotateX3D(M_PI)*G4RotateZ3D(-M_PI / 2 + M_PI / 8 + i * M_PI / 4)*G4Translate3D(0,
364 1496 - 185 + 359. / 2,
365 434 + 5 + (257. - 5.) / 2), l_all, "ECL_BackwardSupport_lall_physical", top, false, i, overlap);
366
367
368 G4VSolid* s1a = new G4Box("leg_p1a", 130. / 2, 178. / 2, 5. / 2);
369 G4LogicalVolume* l1a = new G4LogicalVolume(s1a, Materials::get("SUS304"), "l1a", 0, 0, 0);
370 l1a->SetVisAttributes(batt);
371 for (int i = 0; i < 8; i++)
372 new G4PVPlacement(gTrans * G4RotateX3D(M_PI)*G4RotateZ3D(-M_PI / 2 + M_PI / 8 + i * M_PI / 4)*G4Translate3D(0,
373 1496 - 185 + 178. / 2,
374 434 + 5 - 5. / 2), l1a, "l1a_physical", top, false, i, overlap);
375
376 delete support_leg;
377 }
378
379}
Belle2::geometry::Materials::get
static G4Material * get(const std::string &name)
Find given material.
Definition Materials.h:63
Belle2::tan
double tan(double a)
tan for double
Definition beamHelpers.h:31
Belle2::atan
double atan(double a)
atan for double
Definition beamHelpers.h:34
Belle2::sqrt
double sqrt(double a)
sqrt for double
Definition beamHelpers.h:28
Belle2::ECLElementNumbers::c_NCrystalsForwardBarrel
const int c_NCrystalsForwardBarrel
Number of crystals in the forward and barrel ECL.
Definition ECLElementNumbers.h:35

◆ barrel()

void barrel ( G4LogicalVolume & _top)
private

Make the ECL barrel and then place elements inside it.

Definition at line 63 of file barrel.cc.

64{
65 G4LogicalVolume* top = &_top;
66
67 // vector<cplacement_t> bp = load_placements("/ecl/data/crystal_placement_barrel.dat");
68 vector<cplacement_t> bp = load_placements(m_sap, ECLParts::barrel);
69
70 vector<cplacement_t>::iterator fp = find_if(bp.begin(), bp.end(), [](const cplacement_t& p) {
71 const int ECL_barrel_part = 1001;
72 return p.nshape == ECL_barrel_part;
73 });
74 // global transformation before placing the whole forward part in the top logical volume
75 G4Transform3D gT = (fp == bp.end()) ? G4Translate3D(0, 0, 0) : get_transform(*fp);
76 if (fp != bp.end()) bp.erase(fp); // now not needed
77
78 const bool b_crystals = true;
79 const bool b_forward_support_legs = true;
80 const bool b_backward_support_legs = true;
81 const bool b_forward_support_ring = true;
82 const bool b_backward_support_ring = true;
83 const bool b_forward_support_wall = true;
84 const bool b_backward_support_wall = true;
85 const bool b_crystal_support = true;
86 const bool b_preamplifier = true;
87 const bool b_septum_walls = true;
88 const bool b_inner_support_wall = true;
89 const bool b_support_ribs = true;
90 const bool b_outer_plates = true;
91 const bool b_forward_part5 = true;
92 const bool b_crystal_holder = true;
93 const bool b_pipe = true;
94
95 const double phi0 = M_PI / 2 - M_PI / 72;
96 const double dphi = 2 * M_PI; //2*2*M_PI/72;
97
98 const int nseg = 72;
99
100 int overlap = m_overlap;
101
102 G4LogicalVolume* sectorlogical; // tilted sector
103 {
104 double R = 1251.6;
105 double zp = 4.9, zm = 2.5;
106 // zr_t bint[] = {{-1000.8-zm, 0}, {1982.6 + zp,0}, {2288.,198.2 - zp*tand(32.982)}, {2288.,404 + 12}, {-1219.,404 + 12}, {-1219., 404 -112.6-zm*tand(52.897)-3.2}};
107 zr_t bint[] = {{ -1000.8 - zm, 0}, {1982.6 + zp, 0}, {2288., 198.2 - zp * tand(32.982)}, {2288., 1582 - R}, {2380, 1582 - R}, {2380, 404 + 12}, { -1330., 404 + 12}, { -1330., 1543 + 84 - R}, { -1219., 1543 + 84 - R}, { -1219., 404 - 112.6 - zm * tand(52.897) - 3.2}};
108 for (unsigned int i = 0; i < sizeof(bint) / sizeof(zr_t); i++) bint[i].r += R;
109 std::vector<zr_t> contourb(bint, bint + sizeof(bint) / sizeof(zr_t));
110 G4VSolid* sect0 = new BelleLathe("sect0", phi0, 2 * M_PI, contourb);
111
112 G4Vector3D n0(cos(phi0), sin(phi0), 0);
113 // shift r0 by half of septum wall thickness so from one side of the sector the wall will be completely inside the sector volume
114 G4Point3D r0 = R * n0 + 0.5 / 2 * (G4RotateZ3D(M_PI / 2) * n0), rt0 = r0 + G4RotateZ3D(M_PI / 144 - M_PI / 72) * (G4Vector3D(0, 417,
115 0)), r1 = G4RotateZ3D(2 * M_PI / 72) * r0, rt1 = G4RotateZ3D(2 * M_PI / 72) * rt0;
116 G4Transform3D t0 = G4RotateZ3D(-M_PI / 144) * G4Translate3D(-r0);
117 G4Point3D tr0 = t0 * r0, trt0 = t0 * rt0, tr1 = t0 * r1, trt1 = t0 * rt1;
118 G4Vector3D tn1 = trt1 - tr1;
119 double v0 = -(tr1.y() - tr0.y()) / tn1.y();
120 double v1 = -(tr1.y() - trt0.y()) / tn1.y();
121
122 G4Point3D u0 = tr1 + v0 * tn1, u1 = tr1 + v1 * tn1;
123 G4VSolid* sect1 = new G4Trd("sect1", -u0.x() / 2, (trt0 - u1).x() / 2, 2000, 2000, trt0.y() / 2);
124 G4IntersectionSolid* sector = new G4IntersectionSolid("sector", sect0, sect1,
125 G4Translate3D(r0)*G4RotateZ3D(M_PI / 144)*G4Translate3D(u0.x() / 2, trt0.y() / 2, (2288 - 1219.) / 2)*G4RotateX3D(-M_PI / 2));
126 sectorlogical = new G4LogicalVolume(sector, Materials::get("G4_AIR"), "ECLBarrelSectorLogical", 0, 0, 0);
127 sectorlogical->SetVisAttributes(att("air"));
128
129 // Set up region for production cuts
130 G4Region* aRegion = new G4Region("ECLBarrelSector");
131 sectorlogical->SetRegion(aRegion);
132 aRegion->AddRootLogicalVolume(sectorlogical);
133
134 for (int i = 0; i < nseg; i++) {
135 double phi = i * M_PI / 36 - M_PI / 2;
136 new G4PVPlacement(gT * G4RotateZ3D(phi), sectorlogical, suf("ECLBarrelSectorPhysical", i), top, false, i, overlap);
137 }
138 }
139
140 // phi septum wall
141 if (b_septum_walls) {
142 double zh = (0.5 - 0.001) / 2;
143 G4ThreeVector psw[] = {G4ThreeVector(-1000.8, 0, -zh), G4ThreeVector(1982.6, 0, -zh), G4ThreeVector(2288, 198.2, -zh),
144 G4ThreeVector(2288, 404, -zh), G4ThreeVector(-1219., 404, -zh), G4ThreeVector(-1219., 404 - 112.6, -zh),
145 G4ThreeVector(-1000.8, 0, zh), G4ThreeVector(1982.6, 0, zh), G4ThreeVector(2288, 198.2, zh),
146 G4ThreeVector(2288, 404, zh), G4ThreeVector(-1219., 404, zh), G4ThreeVector(-1219., 404 - 112.6, zh)
147 };
148 G4VSolid* pswsolid = new BelleCrystal("pswsolid", 6, psw);
149 G4LogicalVolume* pswlogical = new G4LogicalVolume(pswsolid, Materials::get("A5052"), "pswlogical", 0, 0, 0);
150 pswlogical->SetVisAttributes(att("alum"));
151 new G4PVPlacement(G4RotateZ3D(2.5 / 180 * M_PI)*G4Translate3D(0, 1251.6 + 0.01,
152 0)*G4RotateZ3D(1.25 / 180 * M_PI)*G4RotateY3D(-M_PI / 2), pswlogical, "phiwallphysical", sectorlogical, false, 0, overlap);
153 }
154
155 if (b_forward_support_ring) {
156 zr_t vc4[] = {{0, 1449.2}, {0, 1582}, {90, 1582}, {90, 1582 - 75.7}};
157 std::vector<zr_t> contour4(vc4, vc4 + 4);
158 G4VSolid* part4solid = new BelleLathe("part4solid", phi0, dphi, contour4);
159 G4LogicalVolume* part4logical = new G4LogicalVolume(part4solid, Materials::get("SUS304"), "part4logical", 0, 0, 0);
160 part4logical->SetVisAttributes(att("iron"));
161 new G4PVPlacement(gT * G4Translate3D(0, 0, 2290), part4logical, "ECL_BarrelSupport_part4physical", top, false, 0, overlap);
162 }
163
164 if (b_backward_support_ring) {
165 zr_t vc6[] = {{0, 1543}, { -105, 1543}, { -105, 1543 + 84}, {0, 1543 + 84}};
166 std::vector<zr_t> contour6(vc6, vc6 + 4);
167 G4VSolid* part6solid = new BelleLathe("part6solid", phi0, dphi, contour6);
168 G4LogicalVolume* part6logical = new G4LogicalVolume(part6solid, Materials::get("SUS304"), "part6logical", 0, 0, 0);
169 part6logical->SetVisAttributes(att("iron"));
170 new G4PVPlacement(gT * G4Translate3D(0, 0, -1225), part6logical, "part6physical", top, false, 0, overlap);
171 }
172
173 if (b_forward_support_wall) {
174 zr_t vc2[] = {{0, 1180.1}, {0, 1249.5}, {54.3, 1249.5}, {445.8, 1503.6}, {445.8 + 4 * sind(32.982), 1503.6 - 4 * cosd(32.982)}, {54.3 + 4 * sind(32.982 / 2), 1249.5 - 4}, {4, 1249.5 - 4}, {4, 1180.1}};
175 std::vector<zr_t> contour2(vc2, vc2 + 8);
176 G4VSolid* part2solid = new BelleLathe("part2solid", phi0, dphi, contour2);
177 G4LogicalVolume* part2logical = new G4LogicalVolume(part2solid, Materials::get("A5083"), "part2logical", 0, 0, 0);
178 part2logical->SetVisAttributes(att("alum"));
179 new G4PVPlacement(gT * G4Translate3D(0, 0, 1930), part2logical, "ECL_BarrelSupport_part2physical", top, false, 0, overlap);
180 }
181
182 if (b_backward_support_wall) {
183 zr_t vc3[] = {{0, 1180.1}, { -4, 1180.1}, { -4, 1249.5 - 4}, { -61.9 - 4 * sind(52.897 / 2), 1249.5 - 4}, { -285, 1539.1 - 4 * sind(52.897 / 2)}, { -285, 1618.8}, { -285 + 4, 1618.8}, { -285 + 4, 1539.1}, { -61.9, 1249.5}, {0, 1249.5}};
184 std::vector<zr_t> contour3(vc3, vc3 + 10);
185 G4VSolid* part3solid = new BelleLathe("part3solid", phi0, dphi, contour3);
186 G4LogicalVolume* part3logical = new G4LogicalVolume(part3solid, Materials::get("A5083"), "part3logical", 0, 0, 0);
187 part3logical->SetVisAttributes(att("alum"));
188 new G4PVPlacement(gT * G4Translate3D(0, 0, -940), part3logical, "ECL_BarrelSupport_part3physical", top, false, 0, overlap);
189 }
190
191 if (b_inner_support_wall) {
192 G4VSolid* part1solid = new G4Tubs("part1solid", 1250.1, 1250.1 + 1.5, (988.5 + 1972.5) / 2, phi0, dphi);
193 G4LogicalVolume* part1logical = new G4LogicalVolume(part1solid, Materials::get("A5083"), "part1logical", 0, 0, 0);
194 part1logical->SetVisAttributes(att("alum2"));
195 new G4PVPlacement(gT * G4Translate3D(0, 0, (-988.5 + 1972.5) / 2), part1logical, "ECL_BarrelSupport_part1physical", top, false, 0,
196 overlap);
197 }
198
199 // cppcheck-suppress knownConditionTrueFalse
200 if (b_support_ribs || b_outer_plates || b_forward_part5) {
201 // double zh = (1219.+2288.)/2, dz = (-1219.+2288.)/2;
202 double zh = (1330. + 2380.) / 2, dz = (-1330. + 2380.) / 2;
203
204 if (b_support_ribs) {
205 G4ThreeVector vc9[] = {G4ThreeVector(0, 0, -zh), G4ThreeVector(25, 0, -zh), G4ThreeVector(25, 24.5 - 25.0 * tand(2.5), -zh), G4ThreeVector(0, 24.5, -zh),
206 G4ThreeVector(0, 0, zh), G4ThreeVector(25, 0, zh), G4ThreeVector(25, 24.5 - 25.0 * tand(2.5), zh), G4ThreeVector(0, 24.5, zh)
207 };
208 G4VSolid* part9solid = new BelleCrystal("part9solid", 4, vc9);
209 G4LogicalVolume* part9logical = new G4LogicalVolume(part9solid, Materials::get("SUS304"), "part9logical", 0, 0, 0);
210 part9logical->SetVisAttributes(att("silv"));
211 G4Transform3D tpart9a = G4RotateZ3D(2.5 / 180 * M_PI) * G4Translate3D(0, 1251.6,
212 0) * G4RotateZ3D(1.25 / 180 * M_PI) * G4Translate3D(0.5 / 2, 404 - 25 - 3.7 + 0.2, dz);
213 G4Transform3D tpart9b = G4RotateZ3D(-2.5 / 180 * M_PI) * G4Translate3D(0, 1251.6,
214 0) * G4RotateZ3D(1.25 / 180 * M_PI) * G4Translate3D(-0.5 / 2, 404 - 25 - 3.7 + 0.5, dz) * G4RotateY3D(M_PI);
215 new G4PVPlacement(tpart9a, part9logical, "part9aphysical", sectorlogical, false, 0, overlap);
216 new G4PVPlacement(tpart9b, part9logical, "part9bphysical", sectorlogical, false, 0, overlap);
217 }
218
219 if (b_outer_plates) {
220 G4VSolid* part10solid = new G4Box("part10solid", 143. / 2, 8. / 2, zh);
221 G4LogicalVolume* part10logical = new G4LogicalVolume(part10solid, Materials::get("SUS304"), "part10logical", 0, 0, 0);
222 part10logical->SetVisAttributes(att("iron2"));
223 G4Transform3D tpart10 = G4Translate3D(0, 1251.6, 0) * G4Translate3D(-404 * tand(1.25), 404 - 1.1,
224 dz) * G4RotateZ3D(0.3 / 180 * M_PI);
225 new G4PVPlacement(tpart10, part10logical, "part10physical", sectorlogical, false, 0, overlap);
226 }
227
228 // G4VSolid *part5solid = new G4Box("part5solid", 60./2, 45./2, 90./2);
229 // G4LogicalVolume* part5logical = new G4LogicalVolume(part5solid, Materials::get("SUS304"), "part5logical", 0, 0, 0);
230 // part5logical->SetVisAttributes(att("iron2"));
231 // new G4PVPlacement(G4RotateZ3D(2.796/180*M_PI)*G4Translate3D(0,1582+45./2,2380.-90./2), part5logical, "part5physical", sectorlogical, false, 1, overlap);
232
233
234 if (b_forward_part5) { // slice the element to fit into the sector
235 double R = 1251.6, phi1 = phi0 + 1.25 * M_PI / 180, phi2 = (2.796 + 90 - 5) / 180 * M_PI;
236 G4Vector3D n0(cos(phi0), sin(phi0), 0), n1(cos(phi1), sin(phi1), 0), n2(cos(phi2), sin(phi2), 0);
237 // shift r0 by half of septum wall thickness so from one side of the sector the wall will be completely inside the sector volume
238 G4Point3D r0 = R * n0 + 0.5 / 2 * (G4RotateZ3D(M_PI / 2) * n0);
239
240 // n2*(r0 + t*n1 - n2*1582) = 0;
241 double t0 = -n2 * (r0 - n2 * 1582) / (n2 * n1);
242 double t1 = -n2 * (r0 - n2 * (1582 + 45)) / (n2 * n1);
243
244 G4Point3D x0 = r0 + t0 * n1;
245 G4Point3D x1 = (t1 - t0) * n1;
246
247 Point_t c5a[] = {{0, 0}, {30 * cosd(5 - 2.796), -30 * sind(5 - 2.796)}, {30 * cosd(5 - 2.796) + 45 * cos(phi2), -30 * sind(5 - 2.796) + 45 * sin(phi2)}, {x1.x(), x1.y()}};
248 G4ThreeVector vc5a[8];
249 for (int i = 0; i < 4; i++) vc5a[i] = G4ThreeVector(c5a[i].x, c5a[i].y, -90. / 2);
250 for (int i = 0; i < 4; i++) vc5a[i + 4] = G4ThreeVector(c5a[i].x, c5a[i].y, 90. / 2);
251 G4VSolid* part5asolid = new BelleCrystal("part5asolid", 4, vc5a);
252 G4LogicalVolume* part5alogical = new G4LogicalVolume(part5asolid, Materials::get("SUS304"), "part5alogical", 0, 0, 0);
253 part5alogical->SetVisAttributes(att("iron2"));
254 new G4PVPlacement(G4RotateZ3D(2 * M_PI / 72)*G4Translate3D(x0.x(), x0.y(), 2380. - 90. / 2), part5alogical, "part5aphysical",
255 sectorlogical, false, 0, overlap);
256
257 Point_t c5b[] = {{ -30 * cosd(5 - 2.796), 30 * sind(5 - 2.796)}, {0, 0}, {x1.x(), x1.y()}, { -30 * cosd(5 - 2.796) + 44.8 * cos(phi2), 30 * sind(5 - 2.796) + 44.8 * sin(phi2)}};
258 G4ThreeVector vc5b[8];
259 for (int i = 0; i < 4; i++) vc5b[i] = G4ThreeVector(c5b[i].x, c5b[i].y, -90. / 2);
260 for (int i = 0; i < 4; i++) vc5b[i + 4] = G4ThreeVector(c5b[i].x, c5b[i].y, 90. / 2);
261 G4VSolid* part5bsolid = new BelleCrystal("part5bsolid", 4, vc5b);
262 G4LogicalVolume* part5blogical = new G4LogicalVolume(part5bsolid, Materials::get("SUS304"), "part5blogical", 0, 0, 0);
263 part5blogical->SetVisAttributes(att("iron2"));
264 new G4PVPlacement(G4Translate3D(x0.x(), x0.y(), 2380. - 90. / 2), part5blogical, "part5bphysical", sectorlogical, false, 0,
265 overlap);
266 }
267
268 }
269
270 if (b_forward_support_legs) {
271 G4ThreeVector vc11a[] = {G4ThreeVector(0, 0, -63. / 2), G4ThreeVector(35, 24.5, -63. / 2), G4ThreeVector(35, 158 - 35, -63. / 2), G4ThreeVector(0, 158 - 35, -63. / 2),
272 G4ThreeVector(0, 0, 63. / 2), G4ThreeVector(35, 24.5, 63. / 2), G4ThreeVector(35, 158 - 35, 63. / 2), G4ThreeVector(0, 158 - 35, 63. / 2)
273 };
274 G4VSolid* part11asolid = new BelleCrystal("part11asolid", 4, vc11a);
275 G4VSolid* part11bsolid = new G4Box("part11bsolid", 63. / 2, 35. / 2, 100. / 2);
276 G4LogicalVolume* part11alogical = new G4LogicalVolume(part11asolid, Materials::get("SUS304"), "part11alogical", 0, 0, 0);
277 part11alogical->SetVisAttributes(att("iron2"));
278 G4LogicalVolume* part11blogical = new G4LogicalVolume(part11bsolid, Materials::get("SUS304"), "part11blogical", 0, 0, 0);
279 part11blogical->SetVisAttributes(att("iron"));
280
281 double Ro = 158.4;
282 zr_t vsf[] = {{0, 0}, {35, 24.5}, {35, 158 - 35}, {100, 158 - 35}, {100, Ro}, {0, Ro}};
283 for (zr_t* it = vsf; it != vsf + 6; it++) {it->r += 1668 - 158;}
284 std::vector<zr_t> csf(vsf, vsf + 6);
285 G4VSolid* sf = new BelleLathe("sf", 0, 2 * M_PI, csf);
286 G4LogicalVolume* sfl = new G4LogicalVolume(sf, Materials::get("G4_AIR"), "supportfwd", 0, 0, 0);
287 sfl->SetVisAttributes(att("air"));
288 new G4PVPlacement(gT * G4Translate3D(0, 0, 2380), sfl, "supportfwdphysical", top, false, 0, overlap);
289
290 G4VSolid* sfs = new BelleLathe("sfs", -M_PI / 72, 2 * M_PI / 72, csf);
291 G4LogicalVolume* sfsl = new G4LogicalVolume(sfs, Materials::get("G4_AIR"), "supportfwdsector", 0, 0, 0);
292 sfsl->SetVisAttributes(att("air"));
293
294 new G4PVPlacement(G4Translate3D(1668 - 158, 0, 0)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2), part11alogical, "part11aphysical",
295 sfsl, false, 0, overlap);
296 new G4PVPlacement(G4Translate3D(1668 - 158, 0, 0)*G4RotateZ3D(-M_PI / 2)*G4Translate3D(0, (158. - 35.) + 35. / 2, 100. / 2),
297 part11blogical, "part11bphysical", sfsl, false, 0, overlap);
298 new G4PVReplica("supportfwdrep", sfsl, sfl, kPhi, nseg, 2 * M_PI / 72, (2.796 - 2.5)*M_PI / 180);
299 }
300
301 if (b_backward_support_legs) {
302 G4VSolid* part12asolid = new G4Box("part12asolid", 90. / 2, (118. - 35.) / 2, 35. / 2);
303 G4VSolid* part12bsolid = new G4Box("part12bsolid", 90. / 2, 35. / 2, 120. / 2);
304 G4LogicalVolume* part12alogical = new G4LogicalVolume(part12asolid, Materials::get("SUS304"), "part12alogical", 0, 0, 0);
305 part12alogical->SetVisAttributes(att("iron2"));
306 G4LogicalVolume* part12blogical = new G4LogicalVolume(part12bsolid, Materials::get("SUS304"), "part12blogical", 0, 0, 0);
307 part12blogical->SetVisAttributes(att("iron"));
308
309 double Rb = 118.7;
310 zr_t vsb[] = {{0, 0}, {0, Rb}, { -120, Rb}, { -120, 118 - 35}, { -35, 118 - 35}, { -35, 0}};
311 for (zr_t* it = vsb; it != vsb + 6; it++) {it->r += 1668 - 118;}
312 std::vector<zr_t> csb(vsb, vsb + 6);
313 G4VSolid* sb = new BelleLathe("sb", 0, 2 * M_PI, csb);
314 G4LogicalVolume* sbl = new G4LogicalVolume(sb, Materials::get("G4_AIR"), "supportbkw", 0, 0, 0);
315 sbl->SetVisAttributes(att("air"));
316 new G4PVPlacement(gT * G4Translate3D(0, 0, -1330), sbl, "supportbkwphysical", top, false, 0, overlap);
317
318 G4VSolid* sbs = new BelleLathe("sbs", -M_PI / 72, 2 * M_PI / 72, csb);
319 G4LogicalVolume* sbsl = new G4LogicalVolume(sbs, Materials::get("G4_AIR"), "supportbkwsector", 0, 0, 0);
320 sbsl->SetVisAttributes(att("air"));
321
322 new G4PVPlacement(G4Translate3D(1668 - 118, 0, 0)*G4RotateZ3D(-M_PI / 2)*G4Translate3D(0, (118. - 35.) / 2, -35. / 2),
323 part12alogical, "part12aphysical", sbsl, false, 0, overlap);
324 new G4PVPlacement(G4Translate3D(1668 - 118, 0, 0)*G4RotateZ3D(-M_PI / 2)*G4Translate3D(0, (118. - 35.) + 35. / 2, -120. / 2),
325 part12blogical, "part12bphysical", sbsl, false, 0, overlap);
326
327 new G4PVReplica("supportbkwrep", sbsl, sbl, kPhi, nseg, 2 * M_PI / 72, (2.796 - 2.5)*M_PI / 180);
328 }
329
330 if (b_septum_walls) {
331 struct zwall_t { double zlow, zup;};
332 zwall_t zs[] = {{ -682.1, -883.9}, { -445.9, -571.7}, { -220.8, -274.2}, {0, 0}, {220.8, 274.2}, {445.9, 571.7}, {682.1, 883.9}, {936.6, 1220.5}, {1217.1, 1591.2}, {1532., 2007.4}};
333
334 const double dz = 0.5 / 2; // half of theta fin thickness
335
336 G4Vector3D n0(sind(90 + 2.5 + 1.25), -cosd(90 + 2.5 + 1.25), 0); // normal vector of the left sector wall pointing inside the sector
337 G4Point3D r0 = G4Point3D(-1251.6 * sind(2.5), 0, 0) + dz * n0; // point in the left sector wall shifted by phi fin
338 G4Vector3D n1(-sind(90 - 2.5 + 1.25), cosd(90 - 2.5 + 1.25),
339 0); // normal vector of the right sector wall pointing inside the sector
340 G4Point3D r1 = G4Point3D(1251.6 * sind(2.5), 0, 0) + dz * n1; // point in the right sector wall shifted by phi fin
341
342 G4Vector3D n2(0, 1, 0); G4Point3D r2(0, 0, 0), r3(0, 310, 0);
343
344 for (unsigned int i = 0; i < sizeof(zs) / sizeof(zwall_t); i++) {
345 double th = atan2(404, -(zs[i].zlow - zs[i].zup));
346
347 G4Vector3D nf(0, -cos(th), sin(th));
348 // take into account finite thickness of the theta fin
349 G4Point3D rf(0, 0, copysign(dz / sin(th), cos(th)));
350
351 // calculate three plane intersection point
352 auto inter = [](const G4Vector3D & l_n0, const G4Point3D & l_r0,
353 const G4Vector3D & l_n1, const G4Point3D & l_r1,
354 const G4Vector3D & l_n2, const G4Point3D & l_r2) -> G4Point3D {
355 CLHEP::HepMatrix A(3, 3);
356 CLHEP::HepVector B(3);
357 A[0][0] = l_n0.x(), A[0][1] = l_n0.y(), A[0][2] = l_n0.z();
358 A[1][0] = l_n1.x(), A[1][1] = l_n1.y(), A[1][2] = l_n1.z();
359 A[2][0] = l_n2.x(), A[2][1] = l_n2.y(), A[2][2] = l_n2.z();
360
361 B[0] = l_r0 * l_n0;
362 B[1] = l_r1 * l_n1;
363 B[2] = l_r2 * l_n2;
364
365 CLHEP::HepVector r = A.inverse() * B;
366 G4Point3D res(r[0], r[1], r[2]);
367 return res;
368 };
369
370 G4Transform3D tfin(G4RotateX3D(-th + M_PI / 2)), tfini(G4RotateX3D(th - M_PI / 2));
371 G4ThreeVector t0 = tfini * inter(n0, r0, n2, r2, nf, rf);
372 G4ThreeVector t1 = tfini * inter(n1, r1, n2, r2, nf, rf);
373 G4ThreeVector t2 = tfini * inter(n0, r0, n2, r3, nf, rf);
374 G4ThreeVector t3 = tfini * inter(n1, r1, n2, r3, nf, rf);
375 // G4cout<<t0<<" "<<t1<<" "<<t2<<" "<<t3<<" "<<G4RotateX3D(th-M_PI/2)*G4Point3D(0,404,zs[i].zup-zs[i].zlow)<<G4endl;
376
377 G4ThreeVector thfin[8];
378 thfin[0] = G4ThreeVector(t0.x(), t0.y(), -dz);
379 thfin[1] = G4ThreeVector(t1.x(), t1.y(), -dz);
380 thfin[2] = G4ThreeVector(t3.x(), t3.y(), -dz);
381 thfin[3] = G4ThreeVector(t2.x(), t2.y(), -dz);
382 thfin[4] = G4ThreeVector(t0.x(), t0.y(), +dz);
383 thfin[5] = G4ThreeVector(t1.x(), t1.y(), +dz);
384 thfin[6] = G4ThreeVector(t3.x(), t3.y(), +dz);
385 thfin[7] = G4ThreeVector(t2.x(), t2.y(), +dz);
386
387 G4VSolid* thfinsolid = new BelleCrystal("thfinsolid", 4, thfin);
388 G4LogicalVolume* thfinlogical = new G4LogicalVolume(thfinsolid, Materials::get("A5052"), "thfinlogical", 0, 0, 0);
389 thfinlogical->SetVisAttributes(att("iron2"));
390 new G4PVPlacement(G4Translate3D(0, 1251.6, zs[i].zlow)*tfin, thfinlogical, "thfinphysical", sectorlogical, false, 0, overlap);
391 }
392 }
393
394 if (b_pipe) {
395 double zh = (1219. + 2288.) / 2, dz = (-1219. + 2288.) / 2;
396 G4VSolid* tubec = new G4Tubs("tubec", 0, 12. / 2, zh, 0, 2 * M_PI);
397 G4VSolid* tubew = new G4Tubs("tubew", 0, 10. / 2, zh, 0, 2 * M_PI);
398 G4LogicalVolume* tubeclogical = new G4LogicalVolume(tubec, Materials::get("C1220"), "tubec", 0, 0, 0);
399 G4LogicalVolume* tubewlogical = new G4LogicalVolume(tubew, Materials::get("G4_WATER"), "tubew", 0, 0, 0);
400 new G4PVPlacement(G4Translate3D(0, 0, 0), tubewlogical, "tubewphysical", tubeclogical, false, 1, overlap);
401 tubeclogical->SetVisAttributes(att("iron2"));
402 new G4PVPlacement(G4Translate3D(25, 1640, dz), tubeclogical, "tubecphysical", sectorlogical, false, 0, overlap);
403 new G4PVPlacement(G4Translate3D(-35, 1640, dz), tubeclogical, "tubecphysical", sectorlogical, false, 1, overlap);
404 }
405
406 if (b_preamplifier) {
407 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
408 G4Transform3D twc = get_transform(*it);
409 int indx = it - bp.begin();
410 auto pv = new G4PVPlacement(twc * G4TranslateZ3D(300 / 2 + 0.250 + get_pa_box_height() / 2), get_preamp(),
411 suf("phys_barrel_preamplifier", indx),
412 sectorlogical, false, 72 + 9 * (indx / 2) + (indx % 2), 0);
413 if (overlap)pv->CheckOverlaps(1000);
414 }
415 }
416
417 if (b_crystal_holder) {
418 G4VSolid* holder0 = new G4Box("holder0", 19. / 2, 38. / 2, 6. / 2);
419 G4VSolid* holder1 = new G4Box("holder1", 7. / 2, 26. / 2, 7. / 2);
420 G4VSolid* holder = new G4SubtractionSolid("holder", holder0, holder1);
421 G4LogicalVolume* holderlogical = new G4LogicalVolume(holder, Materials::get("A5052"), "holderlogical", 0, 0, 0);
422 holderlogical->SetVisAttributes(att("holder"));
423 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
424 G4Transform3D twc = get_transform(*it);
425 int indx = it - bp.begin();
426 auto pv = new G4PVPlacement(twc * G4TranslateZ3D(300 / 2 + 0.250 + get_pa_box_height() + 38. / 2)*G4RotateZ3D(M_PI / 2)*G4RotateX3D(
427 M_PI / 2), holderlogical, "holderphysical", sectorlogical, false, indx, 0);
428 if (overlap)pv->CheckOverlaps(1000);
429 }
430 }
431
432 if (b_crystals) {
433 // vector<shape_t*> cryst = load_shapes("/ecl/data/crystal_shape_barrel.dat");
434 vector<shape_t*> cryst = load_shapes(m_sap, ECLParts::barrel);
435 vector<G4LogicalVolume*> wrapped_crystals;
436 for (auto it = cryst.begin(); it != cryst.end(); it++) {
437 shape_t* s = *it;
438 wrapped_crystals.push_back(wrapped_crystal(s, "barrel", 0.17 - 0.0325));
439 }
440
441 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
442 const cplacement_t& t = *it;
443 auto s = find_if(cryst.begin(), cryst.end(), [&t](const shape_t* shape) {return shape->nshape == t.nshape;});
444 if (s == cryst.end()) continue;
445
446 G4Transform3D twc = get_transform(t);
447 int indx = it - bp.begin();
448 new G4PVPlacement(twc, wrapped_crystals[s - cryst.begin()], suf("ECLBarrelWrappedCrystal_Physical", indx), sectorlogical, false,
449 72 + 9 * (indx / 2) + (indx % 2), overlap);
450 }
451
452 for (shape_t* shape : cryst) {
453 delete shape;
454 }
455 }
456
457 if (b_crystal_support) {
458 G4LogicalVolume* tl[8];
459 for (int i = 0; i < 8; i++) {
460 G4VSolid* t = get_crystal_support(i);
461 tl[i] = new G4LogicalVolume(t, Materials::get("SUS304"), "crystal_support_logical", 0, 0, 0);
462 tl[i]->SetVisAttributes(att("iron2"));
463 }
464
465 double offset[] = { -945, -607, -305, -45, 45, 305, 607, 909, 1240, 1628, 2040};
466 const int ns[] = {8, 5, 6, 7, 7, 6, 5, 4, 3, 2, 1};
467 const double offsetb[] = {263.4 - 210 - 30, 27.4, 10.6, 16.1, 26.2, 21.7, 35, 22 + 240};
468 bool flip[] = {1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1};
469 for (int i = 0; i < 11; i++) {
470 G4Transform3D tflip = flip[i] ? G4RotateY3D(-M_PI / 2) : G4RotateY3D(M_PI / 2);
471 G4Transform3D tr0(G4RotateZ3D(2.5 / 180 * M_PI)*G4Translate3D(0, 1251.6,
472 0)*G4RotateZ3D(1.25 / 180 * M_PI)*G4Translate3D(25 + 3. / 2 + 0.5 / 2, 404 - 25 - 3.7 + 24.5 - 25.0 * tand(2.5) - 0.5,
473 offset[i])*tflip * G4Translate3D(-offsetb[ns[i] - 1], 0, 0)*G4RotateX3D(M_PI));
474 new G4PVPlacement(tr0, tl[ns[i] - 1], "crystal_support_physical", sectorlogical, false, 0, overlap);
475 G4Transform3D tr1(G4RotateZ3D(-2.5 / 180 * M_PI)*G4Translate3D(0, 1251.6,
476 0)*G4RotateZ3D(1.25 / 180 * M_PI)*G4Translate3D(-25 - 3. / 2 - 0.5 / 2, 404 - 25 - 3.7 + 24.5 - 25.0 * tand(2.5) - 0.5,
477 offset[i])*tflip * G4Translate3D(-offsetb[ns[i] - 1], 0, 0)*G4RotateX3D(M_PI));
478 new G4PVPlacement(tr1, tl[ns[i] - 1], "crystal_support_physical", sectorlogical, false, 1, overlap);
479 }
480 }
481
482}
R
double R
typedef autogenerated by FFTW
Definition DiscreteCosineTransform_31points.cc:35
Belle2::ECL::GeoECLCreator::get_pa_box_height
double get_pa_box_height() const
Getter for preamplifier box height (hard-coded to be 2)
Definition GeoECLCreator.h:88
Belle2::ECL::GeoECLCreator::att
const G4VisAttributes * att(const std::string &n) const
Define visual attributes.
Definition GeoECLCreator.cc:141
Belle2::ECL::GeoECLCreator::get_preamp
G4LogicalVolume * get_preamp() const
Get Logical volume of preamplifier.
Definition GeoECLCreator.cc:148
Belle2::ECL::GeoECLCreator::wrapped_crystal
G4LogicalVolume * wrapped_crystal(const shape_t *s, const std::string &endcap, double wrapthickness)
Wrapped crystal.
Definition GeoECLCreator.cc:90

◆ create()

void create ( const GearDir & content,
G4LogicalVolume & topVolume,
geometry::GeometryTypes type )
overridevirtual

Function to actually create the geometry, has to be overridden by derived classes.

Parameters
contentGearDir pointing to the parameters which should be used for construction
topVolumeTop volume in which the geometry has to be placed
typeType of geometry to be build

Implements CreatorBase.

Definition at line 77 of file GeoECLCreator.cc.

78{
79 m_sensediode = new SensitiveDiode("ECLSensitiveDiode");
80 G4SDManager::GetSDMpointer()->AddNewDetector(m_sensediode);
81
82 ECLCrystalsShapeAndPosition crystals = loadCrystalsShapeAndPosition();
83 m_sap = &crystals;
84
85 forward(topVolume);
86 barrel(topVolume);
87 backward(topVolume);
88}
Belle2::ECL::GeoECLCreator::barrel
void barrel(G4LogicalVolume &)
Make the ECL barrel and then place elements inside it.
Definition barrel.cc:63
Belle2::ECL::GeoECLCreator::backward
void backward(G4LogicalVolume &)
Place elements inside the backward endcap.
Definition backward.cc:42
Belle2::ECL::GeoECLCreator::forward
void forward(G4LogicalVolume &)
Place elements inside the forward endcap.
Definition forward.cc:45

◆ createFromDB()

void createFromDB ( const std::string & name,
G4LogicalVolume & topVolume,
geometry::GeometryTypes type )
overridevirtual

Function to create the geometry from the Database.

Parameters
namename of the component in the database, could be used to disambiguate multiple components created with the same creator
topVolumeTop volume in which the geometry has to be placed
typeType of geometry to be build

Reimplemented from CreatorBase.

Definition at line 57 of file GeoECLCreator.cc.

58{
59 m_sensediode = new SensitiveDiode("ECLSensitiveDiode");
60 G4SDManager::GetSDMpointer()->AddNewDetector(m_sensediode);
61
62 DBObjPtr<ECLCrystalsShapeAndPosition> crystals;
63 if (!crystals.isValid()) B2FATAL("No crystal's data in the database.");
64 m_sap = &(*crystals);
65
66 forward(topVolume);
67 barrel(topVolume);
68 backward(topVolume);
69}
Belle2::DBAccessorBase::isValid
bool isValid() const
Check whether a valid object was obtained from the database.
Definition DBAccessorBase.h:72

◆ createPayloads()

void createPayloads ( const GearDir & content,
const IntervalOfValidity & iov )
overridevirtual

Function to create the geometry database.

This function should be implemented to convert Gearbox parameters to one or more database payloads

Parameters
contentGearDir pointing to the parameters which should be used for construction
iovinterval of validity to use when generating payloads

Reimplemented from CreatorBase.

Definition at line 71 of file GeoECLCreator.cc.

72{
73 ECLCrystalsShapeAndPosition crystals = loadCrystalsShapeAndPosition();
74 Database::Instance().storeData<ECLCrystalsShapeAndPosition>(&crystals, iov);
75}
Belle2::Database::Instance
static Database & Instance()
Instance of a singleton Database.
Definition Database.cc:41
Belle2::Database::storeData
bool storeData(const std::string &name, TObject *object, const IntervalOfValidity &iov)
Store an object in the database.
Definition Database.cc:140

◆ defineVisAttributes()

void defineVisAttributes ( )
private

Define visual attributes.

Definition at line 119 of file GeoECLCreator.cc.

120{
121 m_atts["wrap"] = new G4VisAttributes(G4Colour(0.5, 0.5, 1.0));
122 m_atts["cryst"] = new G4VisAttributes(G4Colour(0.7, 0.7, 1.0));
123
124 m_atts["iron"] = new G4VisAttributes(G4Colour(1., 0.1, 0.1));
125 m_atts["iron2"] = new G4VisAttributes(G4Colour(1., 0.5, 0.5));
126 m_atts["alum"] = new G4VisAttributes(G4Colour(0.25, 0.25, 1.0, 0.5));
127 m_atts["alum2"] = new G4VisAttributes(G4Colour(0.5, 0.5, 1.0));
128 m_atts["silv"] = new G4VisAttributes(G4Colour(0.9, 0., 0.9));
129 m_atts["air"] = new G4VisAttributes(G4Colour(1., 1., 1.));
130 m_atts["air"]->SetVisibility(false);
131 m_atts["preamp"] = new G4VisAttributes(G4Colour(0.1, 0.1, 0.8));
132 m_atts["plate"] = new G4VisAttributes(G4Colour(0.2, 0.8, 0.2));
133 m_atts["connector"] = new G4VisAttributes(G4Colour(0.1, 0.1, 0.1));
134 m_atts["capacitor"] = new G4VisAttributes(G4Colour(0.1, 0.1, 0.8));
135 m_atts["holder"] = new G4VisAttributes(G4Colour(0.4, 0.8, 0.8));
136 m_atts["asolid6"] = new G4VisAttributes(G4Colour(1., 0.3, 0.2));
137 m_atts["asolid7"] = new G4VisAttributes(G4Colour(1., 0.3, 0.2));
138 m_atts["asolid13"] = new G4VisAttributes(G4Colour(1., 0.5, 0.5));
139}

◆ forward()

void forward ( G4LogicalVolume & _top)
private

Place elements inside the forward endcap.

Definition at line 45 of file forward.cc.

46{
47 G4LogicalVolume* top = &_top;
48
49 const bool sec = 0;
50 const double phi0 = 0;
51 const double dphi = sec ? M_PI / 16 : 2 * M_PI;
52
53 const bool b_inner_support_ring = true;
54 const bool b_outer_support_ring = true;
55 const bool b_support_wall = true;
56 const bool b_ribs = true;
57 const bool b_septum_wall = true;
58 const bool b_crystals = true;
59 const bool b_preamplifier = true;
60 const bool b_support_leg = true;
61 const bool b_support_structure_13 = true;
62 const bool b_support_structure_15 = true;
63 bool b_connectors = true;
64 bool b_boards = true;
65 const bool b_cover = true;
66
67 b_connectors &= b_support_structure_15;
68 b_boards &= b_support_structure_15;
69
70 int overlap = m_overlap;
71
72 // vector<cplacement_t> bp = load_placements("/ecl/data/crystal_placement_forward.dat");
73 vector<cplacement_t> bp = load_placements(m_sap, ECLParts::forward);
74 vector<cplacement_t>::iterator fp = find_if(bp.begin(), bp.end(), [](const cplacement_t& p) {
75 const int ECL_forward_part = 1000;
76 return p.nshape == ECL_forward_part;
77 });
78 // global transformation before placing the whole forward part in the top logical volume
79 G4Transform3D gTrans = (fp == bp.end()) ? G4Translate3D(0, 0, 1960) : get_transform(*fp);
80 // since the calorimeter supporting structures attach to the yoke at
81 // the same place we need to modify supporting legs by milling
82 // necessary thickness of steel when moving forward part outward in
83 // z-direction
84 double milled_thickness = gTrans.dz() - 1960; // [mm]
85 if (fp != bp.end()) bp.erase(fp); // now not needed
86
87 const double th0 = 13.12, th1 = 32.98;
88 const double ZT = 437, ZI = 434, RI = 431, RIp = 532.2, RC = 1200.4, RT = 1415;
89 const double thinnerPart_translation = 1.95; // translation of the thin inner ring of the support ring
90
91 if (b_inner_support_ring) {
92 // This object has a (newly) translated part to eliminate overlaps with ARICH - 2021-07.
93 zr_t vc1[] = {{ZI - 487, 410 - thinnerPart_translation}, {ZT - (RIp - 410 - 20 / cosd(th0)) / tand(th0), 410 - thinnerPart_translation}, {ZT - (RIp - 410 - 20 / cosd(th0)) / tand(th0), 410}, {ZT, RIp - 20 / cosd(th0)}, {ZT, RIp}, {3., RI}, {3., 418 - thinnerPart_translation}, {ZI - 487, 418 - thinnerPart_translation}};
94 std::vector<zr_t> contour1(vc1, vc1 + sizeof(vc1) / sizeof(zr_t));
95 G4VSolid* part1solid = new BelleLathe("fwd_part1solid", phi0, dphi, contour1);
96 G4LogicalVolume* part1logical = new G4LogicalVolume(part1solid, Materials::get("SUS304"), "part1logical", 0, 0, 0);
97 part1logical->SetVisAttributes(att("iron"));
98 new G4PVPlacement(gTrans, part1logical, "ECL_ForwardSupport_part1physical", top, false, 0, overlap);
99 }
100
101 if (b_support_wall) {
102 // solving equation to get L : 3+L*cosd(th1)+1.6*cosd(th1+90) = 434 + 3 - 107.24
103 double L = (ZT - 107.24 - 3 - 1.6 * cosd(th1 + 90)) / cosd(th1);
104 double R0 = 418, R1 = RC;
105 zr_t vc23[] = {{0, R0}, {3, R0}, {3, R1}, {3 + L * cosd(th1), R1 + L * sind(th1)},
106 {3 + L * cosd(th1) + 1.6 * cosd(th1 + 90), R1 + L * sind(th1) + 1.6 * sind(th1 + 90)}, {3 + 1.6 * cosd(th1 + 90), R1 + 1.6 * sind(th1 + 90)}, {0, R1}
107 };
108 std::vector<zr_t> contour23(vc23, vc23 + sizeof(vc23) / sizeof(zr_t));
109 G4VSolid* part23solid = new BelleLathe("fwd_part23solid", phi0, dphi, contour23);
110 G4LogicalVolume* part23logical = new G4LogicalVolume(part23solid, Materials::get("A5052"), "part23logical", 0, 0, 0);
111 part23logical->SetVisAttributes(att("alum"));
112 new G4PVPlacement(gTrans, part23logical, "ECL_ForwardSupport_part23physical", top, false, 0, overlap);
113 }
114
115 if (b_outer_support_ring) {
116 zr_t vc4[] = {{3 + (RT - 20 - RC) / tand(th1), RT - 20}, {ZT, RT - 20}, {ZT, RT}, {3 + (RT - RC) / tand(th1), RT}};
117 std::vector<zr_t> contour4(vc4, vc4 + sizeof(vc4) / sizeof(zr_t));
118 G4VSolid* part4solid = new BelleLathe("fwd_part4solid", phi0, dphi, contour4);
119 G4LogicalVolume* part4logical = new G4LogicalVolume(part4solid, Materials::get("SUS304"), "part4logical", 0, 0, 0);
120 part4logical->SetVisAttributes(att("iron"));
121 new G4PVPlacement(gTrans, part4logical, "ECL_ForwardSupport_part4physical", top, false, 0, overlap);
122 }
123
124 zr_t cont_array_in[] = {{3., RI}, {ZT, RIp}, {ZT, RT - 20}, {3 + (RT - 20 - RC) / tand(th1), RT - 20}, {3, RC}};
125 std::vector<zr_t> contour_in(cont_array_in, cont_array_in + sizeof(cont_array_in) / sizeof(zr_t));
126 G4VSolid* innervolume_solid = new BelleLathe("fwd_innervolume_solid", 0, 2 * M_PI, contour_in);
127 G4LogicalVolume* innervolume_logical = new G4LogicalVolume(innervolume_solid, Materials::get("G4_AIR"),
128 "innervolume_logical", 0, 0, 0);
129 innervolume_logical->SetVisAttributes(att("air"));
130
131 // Set up region for production cuts
132 G4Region* aRegion = new G4Region("ECLForwardEnvelope");
133 innervolume_logical->SetRegion(aRegion);
134 aRegion->AddRootLogicalVolume(innervolume_logical);
135
136 new G4PVPlacement(gTrans, innervolume_logical, "ECLForwardPhysical", top, false, 0, overlap);
137
138 G4VSolid* innervolumesector_solid = new BelleLathe("fwd_innervolumesector_solid", -M_PI / 8, M_PI / 4, contour_in);
139 G4LogicalVolume* innervolumesector_logical = new G4LogicalVolume(innervolumesector_solid, Materials::get("G4_AIR"),
140 "innervolumesector_logical", 0, 0, 0);
141 innervolumesector_logical->SetVisAttributes(att("air"));
142 new G4PVReplica("ECLForwardSectorPhysical", innervolumesector_logical, innervolume_logical, kPhi, 8, M_PI / 4, 0);
143
144 if (b_ribs) {
145 double H = 60, W = 20;
146 double X0 = RIp, X1 = RT - 20;
147 G4TwoVector r0o(X1, 0), r1o(X1 * sqrt(1 - pow(W / X1, 2)), W);
148 double beta = asin(W / X0);
149 G4TwoVector r0i(X0 / cos(beta / 2), 0), r1i(X0 * cos(beta / 2) - tan(beta / 2) * (W - X0 * sin(beta / 2)), W);
150 double dxymzp = (r0o - r0i).x(), dxypzp = (r1o - r1i).x();
151 double theta = atan(tand(th0) / 2);
152 double dxymzm = dxymzp + tand(th0) * H, dxypzm = dxypzp + tand(th0) * H;
153
154 G4TwoVector m0 = (r0i + r0o) * 0.5, m1 = (r1i + r1o) * 0.5, dm = m1 - m0;
155 double alpha = atan(dm.x() / dm.y());
156
157 G4VSolid* solid6_p1 = new G4Trap("fwd_solid6_p1", H / 2, theta, 0, W / 2, dxymzm / 2, dxypzm / 2, alpha, W / 2, dxymzp / 2,
158 dxypzp / 2,
159 alpha);
160 G4LogicalVolume* lsolid6_p1 = new G4LogicalVolume(solid6_p1, Materials::get("SUS304"), "lsolid6", 0, 0, 0);
161 lsolid6_p1->SetVisAttributes(att("iron"));
162 G4Transform3D tsolid6_p1(G4Translate3D(X0 * cos(beta / 2) + (dxymzp / 2 + dxypzp / 2) / 2 - tan(theta)*H / 2, W / 2, ZT - H / 2));
163 new G4PVPlacement(G4RotateZ3D(-M_PI / 8)*tsolid6_p1, lsolid6_p1, "psolid6_p1", innervolumesector_logical, false, 0, overlap);
164 new G4PVPlacement(G4RotateZ3D(0)*tsolid6_p1, lsolid6_p1, "psolid6_p2", innervolumesector_logical, false, 0, overlap);
165
166 H = 40;
167 dxymzm = dxymzp + tand(th0) * H, dxypzm = dxypzp + tand(th0) * H;
168
169 G4Transform3D tsolid6_p2(G4Translate3D(X0 * cos(beta / 2) + (dxymzp / 2 + dxypzp / 2) / 2 - tan(theta)*H / 2, -W / 2, ZT - H / 2));
170
171 // Added intersection to avoid protrusion from the mother volume
172 G4VSolid* solid6_p2_x = new G4Trap("fwd_solid6_p2_x", H / 2, theta, 0, W / 2, dxypzm / 2, dxymzm / 2, -alpha, W / 2, dxypzp / 2,
173 dxymzp / 2,
174 -alpha);
175 G4VSolid* solid6_p2 = new G4IntersectionSolid("fwd_solid6_p2", solid6_p2_x, innervolumesector_solid, tsolid6_p2.inverse());
176
177 G4LogicalVolume* lsolid6_p2 = new G4LogicalVolume(solid6_p2, Materials::get("SUS304"), "lsolid6", 0, 0, 0);
178 lsolid6_p2->SetVisAttributes(att("iron"));
179 new G4PVPlacement(G4RotateZ3D(0)*tsolid6_p2, lsolid6_p2, "psolid6_p3", innervolumesector_logical, false, 0, overlap);
180 new G4PVPlacement(G4RotateZ3D(M_PI / 8)*tsolid6_p2, lsolid6_p2, "psolid6_p4", innervolumesector_logical, false, 0, overlap);
181
182 double hpad = 148.4;
183 G4VSolid* solid7_p8 = new G4Box("fwd_solid7_p8", hpad / 2, (140. - 40) / 2 / 2, 40. / 2);
184 G4LogicalVolume* lsolid7 = new G4LogicalVolume(solid7_p8, Materials::get("SUS304"), "lsolid7", 0, 0, 0);
185 lsolid7->SetVisAttributes(att("iron"));
186 double dx = sqrt(X1 * X1 - 70 * 70) - hpad / 2;
187 G4Transform3D tsolid7_p1(G4Translate3D(dx, -20 - 25, ZT - 40. / 2));
188 new G4PVPlacement(tsolid7_p1, lsolid7, "psolid7_p1", innervolumesector_logical, false, 0, overlap);
189 G4Transform3D tsolid7_p2(G4Translate3D(dx, 20 + 25, ZT - 40. / 2));
190 new G4PVPlacement(tsolid7_p2, lsolid7, "psolid7_p2", innervolumesector_logical, false, 0, overlap);
191
192 double L = X1 - (X0 - tand(th0) * 40) - 10;
193 G4VSolid* solid13 = new G4Box("fwd_solid13", L / 2, 5. / 2, 18. / 2);
194 G4LogicalVolume* lsolid13 = new G4LogicalVolume(solid13, Materials::get("SUS304"), "lsolid13", 0, 0, 0);
195 lsolid13->SetVisAttributes(att("iron"));
196 G4Transform3D tsolid13(G4TranslateZ3D(ZT - 60 + 18. / 2)*G4TranslateY3D(-5. / 2 - 0.5 / 2)*G4TranslateX3D(X0 - tand(
197 th0) * 40 + L / 2 + 5));
198 new G4PVPlacement(tsolid13, lsolid13, "psolid13_p1", innervolumesector_logical, false, 0, overlap);
199 new G4PVPlacement(G4RotateZ3D(M_PI / 8)*tsolid13, lsolid13, "psolid13_p2", innervolumesector_logical, false, 0, overlap);
200 }
201
202
203 double zsep = 125;
204 if (b_septum_wall) {
205 double d = 5;
206 Point_t vin[] = {{ZT - zsep, RIp - tand(th0)* zsep}, {ZT - 60, RIp - tand(th0) * 60}, {ZT - 60, RT - 20 - d}, {ZT - zsep, RT - 20 - d}};
207 const int n = sizeof(vin) / sizeof(Point_t);
208 Point_t c = centerofgravity(vin, vin + n);
209 G4ThreeVector contour_swall[n * 2];
210 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2);
211 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2);
212
213 G4VSolid* septumwall_solid = new BelleCrystal("fwd_septumwall_solid", n, contour_swall);
214
215 G4LogicalVolume* septumwall_logical = new G4LogicalVolume(septumwall_solid, Materials::get("A5052"),
216 "septumwall_logical", 0, 0, 0);
217 septumwall_logical->SetVisAttributes(att("alum2"));
218 new G4PVPlacement(G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, 0), septumwall_logical,
219 "septumwall_physical", innervolumesector_logical, false, 0, overlap);
220
221 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2 / 2);
222 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2 / 2);
223
224 G4VSolid* septumwall2_solid = new BelleCrystal("fwd_septumwall2_solid", n, contour_swall);
225
226 G4LogicalVolume* septumwall2_logical = new G4LogicalVolume(septumwall2_solid, Materials::get("A5052"),
227 "septumwall2_logical", 0, 0, 0);
228 septumwall2_logical->SetVisAttributes(att("alum2"));
229 new G4PVPlacement(G4RotateZ3D(-M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, 0.5 / 2 / 2),
230 septumwall2_logical, "septumwall2_physical", innervolumesector_logical, false, 0, overlap);
231 new G4PVPlacement(G4RotateZ3D(M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, -0.5 / 2 / 2),
232 septumwall2_logical, "septumwall2_physical", innervolumesector_logical, false, 1, overlap);
233 }
234
235 zr_t vcr[] = {{3., RI}, {ZT - zsep, RIp - tand(th0)* zsep}, {ZT - zsep, RT - 20}, {3 + (RT - 20 - RC) / tand(th1), RT - 20}, {3, RC}};
236 std::vector<zr_t> ccr(vcr, vcr + sizeof(vcr) / sizeof(zr_t));
237 G4VSolid* crystalvolume_solid = new BelleLathe("fwd_crystalvolume_solid", 0, M_PI / 8, ccr);
238 G4LogicalVolume* crystalvolume_logical = new G4LogicalVolume(crystalvolume_solid, Materials::get("G4_AIR"),
239 "crystalvolume_logical", 0, 0, 0);
240 crystalvolume_logical->SetVisAttributes(att("air"));
241 new G4PVPlacement(G4RotateZ3D(-M_PI / 8), crystalvolume_logical, "ECLForwardCrystalSectorPhysical_0", innervolumesector_logical,
242 false, 0, overlap);
243 new G4PVPlacement(G4RotateZ3D(0), crystalvolume_logical, "ECLForwardCrystalSectorPhysical_1", innervolumesector_logical, false, 1,
244 overlap);
245
246 if (b_septum_wall) {
247 double d = 5, aRC = RC - 30e-6;
248 Point_t vin[] = {{3., RI}, {ZT - zsep, RIp - tand(th0)* zsep}, {ZT - zsep, RT - 20 - d}, {3 + (RT - 20 - d - aRC) / tand(th1), RT - 20 - d}, {3, aRC}};
249 const int n = sizeof(vin) / sizeof(Point_t);
250 Point_t c = centerofgravity(vin, vin + n);
251 G4ThreeVector contour_swall[n * 2];
252
253 for (int i = 0; i < n; i++) contour_swall[i + 0] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, -0.5 / 2 / 2);
254 for (int i = 0; i < n; i++) contour_swall[i + n] = G4ThreeVector(vin[i].x - c.x, vin[i].y - c.y, 0.5 / 2 / 2);
255
256 G4VSolid* septumwall3_solid = new BelleCrystal("fwd_septumwall3_solid", n, contour_swall);
257
258 G4LogicalVolume* septumwall3_logical = new G4LogicalVolume(septumwall3_solid, Materials::get("A5052"),
259 "septumwall3_logical", 0, 0, 0);
260 septumwall3_logical->SetVisAttributes(att("alum2"));
261 new G4PVPlacement(G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, 0.5 / 2 / 2), septumwall3_logical,
262 "septumwall3_physical_0", crystalvolume_logical, false, 0, overlap);
263 new G4PVPlacement(G4RotateZ3D(M_PI / 8)*G4RotateZ3D(-M_PI / 2)*G4RotateY3D(-M_PI / 2)*G4Translate3D(c.x, c.y, -0.5 / 2 / 2),
264 septumwall3_logical, "septumwall3_physical_1", crystalvolume_logical, false, 1, overlap);
265 }
266
267 if (b_crystals) {
268 // vector<shape_t*> cryst = load_shapes("/ecl/data/crystal_shape_forward.dat");
269 vector<shape_t*> cryst = load_shapes(m_sap, ECLParts::forward);
270 vector<G4LogicalVolume*> wrapped_crystals;
271 for (auto it = cryst.begin(); it != cryst.end(); it++) {
272 shape_t* s = *it;
273 wrapped_crystals.push_back(wrapped_crystal(s, "forward", 0.20 - 0.02));
274 }
275
276 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
277 const cplacement_t& t = *it;
278 auto s = find_if(cryst.begin(), cryst.end(), [&t](const shape_t* shape) {return shape->nshape == t.nshape;});
279 if (s == cryst.end()) continue;
280
281 G4Transform3D twc = G4Translate3D(0, 0, 3) * get_transform(t);
282 int indx = it - bp.begin();
283 new G4PVPlacement(twc, wrapped_crystals[s - cryst.begin()], suf("ECLForwardWrappedCrystal_Physical", indx), crystalvolume_logical,
284 false, indx, overlap);
285 }
286
287 for (shape_t* shape : cryst) {
288 delete shape;
289 }
290 }
291
292 if (b_preamplifier) {
293 for (vector<cplacement_t>::const_iterator it = bp.begin(); it != bp.end(); ++it) {
294 G4Transform3D twc = G4Translate3D(0, 0, 3) * get_transform(*it);
295 int indx = it - bp.begin();
296 auto pv = new G4PVPlacement(twc * G4TranslateZ3D(300 / 2 + 0.20 + get_pa_box_height() / 2)*G4RotateZ3D(-M_PI / 2), get_preamp(),
297 suf("phys_forward_preamplifier", indx), crystalvolume_logical, false, indx, 0);
298 if (overlap)pv->CheckOverlaps(1000);
299 }
300 }
301
302 if (b_support_leg) {
303 const G4VisAttributes* batt = att("iron");
304
305 G4VSolid* s1 = new G4Box("fwd_leg_p1", 130. / 2, 170. / 2, (40. - milled_thickness) / 2);
306 G4LogicalVolume* l1 = new G4LogicalVolume(s1, Materials::get("SUS304"), "l1", 0, 0, 0);
307 G4Transform3D t1 = G4Translate3D(0, 170. / 2, (40. - milled_thickness) / 2 + milled_thickness);
308 l1->SetVisAttributes(batt);
309
310 G4VSolid* s2 = new G4Box("fwd_leg_p2", 60. / 2, 130. / 2, 137. / 2);
311 G4LogicalVolume* l2 = new G4LogicalVolume(s2, Materials::get("SUS304"), "l2", 0, 0, 0);
312 G4Transform3D t2 = G4Translate3D(0, 130. / 2 + 35, 40. + 137. / 2);
313 l2->SetVisAttributes(batt);
314
315 Point_t v3[] = {{ -75. / 2, -265. / 2}, {75. / 2 - 30, -265. / 2}, {75. / 2, -265. / 2 + 30}, {75. / 2, 265. / 2}, { -75. / 2, 265. / 2}};
316 const int n3 = sizeof(v3) / sizeof(Point_t);
317 G4ThreeVector c3[n3 * 2];
318
319 for (int i = 0; i < n3; i++) c3[i + 0] = G4ThreeVector(v3[i].x, v3[i].y, -140. / 2);
320 for (int i = 0; i < n3; i++) c3[i + n3] = G4ThreeVector(v3[i].x, v3[i].y, 140. / 2);
321
322 G4VSolid* s3 = new BelleCrystal("fwd_leg_p3", n3, c3);
323 G4LogicalVolume* l3 = new G4LogicalVolume(s3, Materials::get("SUS304"), "l3", 0, 0, 0);
324 G4Transform3D t3 = G4Translate3D(0, 265. / 2 + 35, 40. + 137. + 75. / 2) * G4RotateY3D(-M_PI / 2);
325 l3->SetVisAttributes(batt);
326
327 G4VSolid* s4 = new G4Box("fwd_leg_p4", 130. / 2, 5. / 2, (5. + milled_thickness) / 2);
328 G4LogicalVolume* l4 = new G4LogicalVolume(s4, Materials::get("SUS304"), "l4", 0, 0, 0);
329 G4Transform3D t4 = G4Translate3D(0, 170. - 5. / 2, (milled_thickness - 5.) / 2);
330 l4->SetVisAttributes(batt);
331
332 G4VSolid* s5 = new G4Box("fwd_leg_p5", 140. / 2, 130. / 2, 80. / 2);
333 G4LogicalVolume* l5 = new G4LogicalVolume(s5, Materials::get("SUS304"), "l5", 0, 0, 0);
334 G4Transform3D t5 = G4Translate3D(0, 180. + 130. / 2, 97. + 80. / 2);
335 l5->SetVisAttributes(batt);
336
337 G4VSolid* s6 = new G4Box("fwd_leg_p6", 140. / 2, 110. / 2, 160. / 2);
338 G4LogicalVolume* l6 = new G4LogicalVolume(s6, Materials::get("G4_AIR"), "l6", 0, 0, 0);
339 G4Transform3D t6 = G4Translate3D(0, 310. + 110. / 2, 97. + 160. / 2);
340 l6->SetVisAttributes(att("air"));
341
342 G4VSolid* s6a = new G4Box("fwd_leg_p6a", 140. / 2, (110. - 45.) / 2, 160. / 2);
343 G4LogicalVolume* l6a = new G4LogicalVolume(s6a, Materials::get("SUS304"), "l6a", 0, 0, 0);
344 l6a->SetVisAttributes(batt);
345 new G4PVPlacement(G4TranslateY3D(-45. / 2), l6a, "l6a_physical", l6, false, 0, overlap);
346
347 G4VSolid* s6b = new G4Box("fwd_leg_p6b", 60. / 2, 45. / 2, 160. / 2);
348 G4LogicalVolume* l6b = new G4LogicalVolume(s6b, Materials::get("SUS304"), "l6b", 0, 0, 0);
349 l6b->SetVisAttributes(batt);
350 double dy = 110. / 2 - 45 + 45. / 2;
351 new G4PVPlacement(G4TranslateY3D(dy), l6b, "l6b_physical", l6, false, 0, overlap);
352
353 G4VSolid* s6c = new G4Box("fwd_leg_p6c", 40. / 2, 45. / 2, 22.5 / 2);
354 G4LogicalVolume* l6c = new G4LogicalVolume(s6c, Materials::get("SUS304"), "l6c", 0, 0, 0);
355 l6c->SetVisAttributes(batt);
356 new G4PVPlacement(G4Translate3D(30 + 20, dy, 20 + 22.5 / 2), l6c, "l6c_physical", l6, false, 0, overlap);
357 new G4PVPlacement(G4Translate3D(30 + 20, dy, -20 - 22.5 / 2), l6c, "l6c_physical", l6, false, 1, overlap);
358 new G4PVPlacement(G4Translate3D(-30 - 20, dy, 20 + 22.5 / 2), l6c, "l6c_physical", l6, false, 2, overlap);
359 new G4PVPlacement(G4Translate3D(-30 - 20, dy, -20 - 22.5 / 2), l6c, "l6c_physical", l6, false, 3, overlap);
360
361 G4AssemblyVolume* support_leg = new G4AssemblyVolume();
362
363 support_leg->AddPlacedVolume(l1, t1);
364 support_leg->AddPlacedVolume(l2, t2);
365 support_leg->AddPlacedVolume(l3, t3);
366 support_leg->AddPlacedVolume(l4, t4);
367 support_leg->AddPlacedVolume(l5, t5);
368 support_leg->AddPlacedVolume(l6, t6);
369
370 for (int i = 0; i < 8; i++) {
371 G4Transform3D tp = gTrans * G4RotateZ3D(-M_PI / 2 + M_PI / 8 + i * M_PI / 4) *
372 G4Translate3D(0, 1415 - 165 + 420. / 2, ZT + (97. + 160.) / 2) * G4Translate3D(0, -420. / 2, -(97. + 160.) / 2);
373 support_leg->MakeImprint(top, tp, 0, overlap);
374 }
375
376 // G4VSolid* s_all = new G4Box("leg_all", 140. / 2, 420. / 2, (97. + 160) / 2);
377 // G4LogicalVolume* l_all = new G4LogicalVolume(s_all, Materials::get("G4_AIR"), "l_all", 0, 0, 0);
378 // l_all->SetVisAttributes(att("silv"));
379 // G4Transform3D tp = G4Translate3D(0, -420. / 2, -(97. + 160.) / 2);
380 // support_leg->MakeImprint(l_all, tp, 0, overlap);
381
382 // for (int i = 0; i < 8; i++)
383 // new G4PVPlacement(G4RotateZ3D(-M_PI / 2 + M_PI / 8 + i * M_PI / 4)*G4Translate3D(0, 1415 - 165 + 420. / 2,
384 // 1960 + ZT + (97. + 160.) / 2), l_all, suf("support_leg_physical", i), top, false, i, overlap);
385
386 // for (int i = 0; i < 8; i++)
387 // new G4PVPlacement(G4RotateZ3D(-M_PI / 2 + M_PI / 8 + i * M_PI / 4)*G4Translate3D(0, 1415 - 165 + 420. / 2,
388 // 1960 + ZT + (97. + 160.) / 2)*tp * t4, l4, suf("support_leg_p4_physical", i), top, false, i, overlap);
389
390 delete support_leg;
391 }
392
393
394 if (b_support_structure_13) { // numbering scheme as in ECL-004K102.pdf page 13
395
396 // Define one layer as one assembly volume
397 G4AssemblyVolume* acs = new G4AssemblyVolume();
398
399 double Z0 = 434, Ro = 1415 - 20;
400
401 // double R2_0 = 667, R2_1 = 1245/cos(M_PI/16), dR2 = R2_1 - R2_0, dz = dR2*cos(M_PI/16);
402 // G4VSolid* sv_crystal_support1 = new G4Trd("sv_crystal_support1", 30/2, 30/2, R2_0*sin(M_PI/16)-40*cos(M_PI/16), R2_1*sin(M_PI/16)-40*cos(M_PI/16), dz/2);
403 // G4LogicalVolume* lv_crystal_support1 = new G4LogicalVolume(sv_crystal_support1, world_mat, "lv_crystal_support1", 0, 0, 0);
404 // // lv_crystal_support->SetVisAttributes(airvol);
405 // new G4PVPlacement(G4Translate3D(R2_0*cos(M_PI/16)+dz/2+3,0,Z0-95)*G4RotateY3D(M_PI/2), lv_crystal_support1, "phys_crystal_support1", crystalSectorLogical, false, 0, overlaps);
406
407 // double dx=425+33, dy1 = dx*tan(M_PI/16), dz = 630+20-67, dy2 = dy1+dz*tan(M_PI/16);
408 // G4VSolid* sv_crystal_support1 = new G4Trd("sv_crystal_support1", 30/2, 30/2, dy1, dy2, dz/2);
409 // G4LogicalVolume* lv_crystal_support1 = new G4LogicalVolume(sv_crystal_support1, world_mat, "lv_crystal_support1", 0, 0, 0);
410 // // lv_crystal_support->SetVisAttributes(airvol);
411 // new G4PVPlacement(G4Translate3D(40/sin(M_PI/16)+dx+dz/2,0,Z0-95)*G4RotateY3D(M_PI/2), lv_crystal_support1, "phys_crystal_support1", crystalSectorLogical, false, 0, overlaps);
412
413 G4VSolid* solid10_p1 = new G4Box("fwd_solid10_p1", 558. / 2 + 9.5, 20. / 2, 105. / 2);
414 G4VSolid* solid10_p2 = new G4Box("fwd_solid10_p2", 559. / 2 + 9.5, 11. / 2, 71. / 2);
415 G4VSolid* solid10_p3 = new G4SubtractionSolid("fwd_solid10_p3", solid10_p1, solid10_p2, G4Translate3D(0, -11. / 2, 71. / 2 - 17.5));
416
417 G4LogicalVolume* lsolid10 = new G4LogicalVolume(solid10_p3, Materials::get("A6063"), "lsolid10", 0, 0, 0);
418 lsolid10->SetVisAttributes(att("alum"));
419 G4Transform3D tsolid10_p1(G4RotateZ3D(M_PI / 16)*G4Translate3D(954.5 + 2.55 - 1, -30, Z0 - 105. / 2 - 5));
420 acs->AddPlacedVolume(lsolid10, tsolid10_p1);
421 // new G4PVPlacement(tsolid10_p1, lsolid10, "psolid10_p1", crystalSectorLogical, false, 0, overlaps);
422 G4Transform3D tsolid10_p2(G4RotateZ3D(-M_PI / 16)*G4Translate3D(954.5 + 2.55 - 1, 30, Z0 - 105. / 2 - 5)*G4RotateZ3D(M_PI));
423 acs->AddPlacedVolume(lsolid10, tsolid10_p2);
424 // new G4PVPlacement(tsolid10_p2, lsolid10, "psolid10_p2", crystalSectorLogical, false, 0, overlaps);
425
426 G4VSolid* solid1_p1 = new G4Box("fwd_solid1_p1", 100. / 2, 30. / 2, 30. / 2);
427 G4VSolid* solid1_p2 = new G4Box("fwd_solid1_p2", 80. / 2, 10. / 2, 31. / 2);
428 G4VSolid* solid1_p3 = new G4SubtractionSolid("fwd_solid1_p3", solid1_p1, solid1_p2, G4Transform3D::Identity);
429
430 G4LogicalVolume* lsolid1 = new G4LogicalVolume(solid1_p3, Materials::get("A5052"), "lsolid1", 0, 0, 0);
431 lsolid1->SetVisAttributes(att("alum"));
432
433 G4Transform3D tsolid1_p1(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 3 * 140, 0, Z0 - 95));
434 acs->AddPlacedVolume(lsolid1, tsolid1_p1);
435 // new G4PVPlacement(tsolid1_p1, lsolid1, "psolid1_p1", crystalSectorLogical, false, 0, overlaps);
436 G4Transform3D tsolid1_p2(G4RotateZ3D(M_PI / 16 * (1 - 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 3 * 140, 0, Z0 - 95));
437 acs->AddPlacedVolume(lsolid1, tsolid1_p2);
438 // new G4PVPlacement(tsolid1_p2, lsolid1, "psolid1_p2", crystalSectorLogical, false, 0, overlaps);
439
440 G4Transform3D tsolid1_p3(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 2 * 140, 0, Z0 - 95));
441 acs->AddPlacedVolume(lsolid1, tsolid1_p3);
442 // new G4PVPlacement(tsolid1_p3, lsolid1, "psolid1_p3", crystalSectorLogical, false, 0, overlaps);
443 G4Transform3D tsolid1_p4(G4RotateZ3D(M_PI / 16 * (1 - 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 2 * 140, 0, Z0 - 95));
444 acs->AddPlacedVolume(lsolid1, tsolid1_p4);
445 // new G4PVPlacement(tsolid1_p4, lsolid1, "psolid1_p4", crystalSectorLogical, false, 0, overlaps);
446
447 G4Transform3D tsolid1_p5(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 1 * 140, 0, Z0 - 95));
448 acs->AddPlacedVolume(lsolid1, tsolid1_p5);
449 // new G4PVPlacement(tsolid1_p5, lsolid1, "psolid1_p5", crystalSectorLogical, false, 0, overlaps);
450 G4Transform3D tsolid1_p6(G4RotateZ3D(M_PI / 16 * (1 - 2 / 3.))*G4Translate3D(Ro - 8 - 50 - 1 * 140, 0, Z0 - 95));
451 acs->AddPlacedVolume(lsolid1, tsolid1_p6);
452 // new G4PVPlacement(tsolid1_p6, lsolid1, "psolid1_p6", crystalSectorLogical, false, 0, overlaps);
453
454 G4Transform3D tsolid1_p8(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3.))*G4Translate3D(Ro - 8 - 50, 0, Z0 - 100));
455 acs->AddPlacedVolume(lsolid1, tsolid1_p8);
456 // new G4PVPlacement(tsolid1_p8, lsolid1, "psolid1_p8", crystalSectorLogical, false, 0, overlaps);
457 G4Transform3D tsolid1_p9(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3. + (1 + 1. / 3) / 3))*G4Translate3D(Ro - 8 - 50, 0, Z0 - 100));
458 acs->AddPlacedVolume(lsolid1, tsolid1_p9);
459 // new G4PVPlacement(tsolid1_p9, lsolid1, "psolid1_p9", crystalSectorLogical, false, 0, overlaps);
460 G4Transform3D tsolid1_p10(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3. + 2 * (1 + 1. / 3) / 3))*G4Translate3D(Ro - 8 - 50, 0, Z0 - 100));
461 acs->AddPlacedVolume(lsolid1, tsolid1_p10);
462 // new G4PVPlacement(tsolid1_p10, lsolid1, "psolid1_p10", crystalSectorLogical, false, 0, overlaps);
463
464 G4VSolid* solid1_p11 = new G4Box("fwd_solid1_p1", 100. / 2, 10. / 2, 30. / 2);
465 G4LogicalVolume* lsolid1_p2 = new G4LogicalVolume(solid1_p11, Materials::get("A5052"), "lsolid1_p2", 0, 0, 0);
466 lsolid1_p2->SetVisAttributes(att("alum"));
467 G4Transform3D tsolid1_p11(G4RotateZ3D(M_PI / 16 * (-1 + 2 / 3. - 1. / 3))*G4Translate3D(Ro - 8 - 50, 0, Z0 - 100));
468 acs->AddPlacedVolume(lsolid1_p2, tsolid1_p11);
469 // new G4PVPlacement(tsolid1_p11, lsolid1_p2, "psolid1_p11", crystalSectorLogical, false, 0, overlaps);
470
471 G4VSolid* solid1_p12 = new G4Box("fwd_solid1_p1", 86. / 2, 10. / 2, 30. / 2);
472 G4LogicalVolume* lsolid1_p3 = new G4LogicalVolume(solid1_p12, Materials::get("A5052"), "lsolid1_p3", 0, 0, 0);
473 lsolid1_p3->SetVisAttributes(att("alum"));
474 double alpha_p12 = M_PI / 16 * (-1 + 1. / 3);
475 G4Transform3D tsolid1_p12(G4Translate3D(532.2 + 43, 0, Z0 - 75));
476 acs->AddPlacedVolume(lsolid1_p3, tsolid1_p12);
477 // new G4PVPlacement(tsolid1_p12, lsolid1_p3, "psolid1_p12", crystalSectorLogical, false, 0, overlaps);
478 G4Transform3D tsolid1_p13(G4RotateZ3D(alpha_p12)*G4Translate3D(532.2 + 43, 0, Z0 - 75));
479 acs->AddPlacedVolume(lsolid1_p3, tsolid1_p13);
480 // new G4PVPlacement(tsolid1_p13, lsolid1_p3, "psolid1_p13", crystalSectorLogical, false, 0, overlaps);
481 G4Transform3D tsolid1_p14(G4RotateZ3D(-alpha_p12)*G4Translate3D(532.2 + 43, 0, Z0 - 75));
482 acs->AddPlacedVolume(lsolid1_p3, tsolid1_p14);
483 // new G4PVPlacement(tsolid1_p14, lsolid1_p3, "psolid1_p14", crystalSectorLogical, false, 0, overlaps);
484
485 G4VSolid* solid1_p4 = new G4Box("fwd_solid1_p4", 160. / 2, 30. / 2, 30. / 2);
486 G4VSolid* solid1_p5 = new G4Box("fwd_solid1_p5", 140. / 2, 10. / 2, 31. / 2);
487 G4VSolid* solid1_p7 = new G4SubtractionSolid("fwd_solid1_p7", solid1_p4, solid1_p5, G4Transform3D::Identity);
488 G4LogicalVolume* lsolid1_p7 = new G4LogicalVolume(solid1_p7, Materials::get("A5052"), "lsolid1_p7", 0, 0, 0);
489 lsolid1_p7->SetVisAttributes(att("alum"));
490 G4Transform3D tsolid1_p7(G4Translate3D(Ro - 8 - 80 - 4 * 140 + 4, 0, Z0 - 95));
491 acs->AddPlacedVolume(lsolid1_p7, tsolid1_p7);
492 // new G4PVPlacement(tsolid1_p7, lsolid1_p7, "psolid1_p7", crystalSectorLogical, false, 0, overlaps);
493
494 // int sol2count = 0;
495 auto get_bracket = [&](double L, double ang, G4Transform3D & lt) {
496 double thick = 3;
497 double dL = (ang > 0) ? 0 : thick * abs(tan(ang));
498
499 G4VSolid* solid2_p1 = new G4Box("fwd_solid2_p1", (L - 2 * dL) / 2, thick / 2, 30. / 2);
500 G4VSolid* solid2_p2 = new G4Box("fwd_solid2_p2", thick / 2, (15. - dL) / 2, 30. / 2);
501 double dx = thick / 2, y0 = (15. + dL) / 2;
502 G4Transform3D t1(G4Translate3D(L / 2, -dx, 0.)*G4RotateZ3D(-ang)*G4Translate3D(-dx, y0, 0));
503 G4Transform3D t2(G4Translate3D(-L / 2, -dx, 0.)*G4RotateZ3D(ang)*G4Translate3D(dx, y0, 0));
504 G4VSolid* solid2_p3 = new G4UnionSolid("fwd_solid2_p3", solid2_p1, solid2_p2, t1);
505 G4VSolid* solid2_p4 = new G4UnionSolid("fwd_solid2_p4", solid2_p3, solid2_p2, t2);
506
507 G4Transform3D u((ang > 0) ? G4Transform3D::Identity : G4RotateZ3D(M_PI));
508 lt = u * G4Translate3D(dx, 0, 0) * G4RotateZ3D(-M_PI / 2);
509
510 return solid2_p4;
511 };
512 double obj2_dz = Z0 - 95;
513 auto place_solid2 = [&](double dz, double L, double ang, double phi, double mx, double dy) {
514 G4Transform3D lt;
515 G4LogicalVolume* lsolid2 = new G4LogicalVolume(get_bracket(L, ang, lt), Materials::get("A5052"), "lsolid2", 0, 0, 0);
516 lsolid2->SetVisAttributes(att("alum"));
517
518 G4Transform3D tsolid2_p1(G4RotateZ3D(phi)*G4Translate3D(mx, dy, dz)*lt);
519 // string pname("psolid2_p"); pname += to_string(++sol2count);
520 acs->AddPlacedVolume(lsolid2, tsolid2_p1);
521 // new G4PVPlacement(tsolid2_p1, lsolid2, pname.c_str(), crystalSectorLogical, false, 0, overlaps);
522 };
523 auto place_solid3 = [&](double L, double ang, const G4Transform3D & t) {
524 G4Transform3D lt;
525 G4LogicalVolume* lsolid2 = new G4LogicalVolume(get_bracket(L, ang, lt), Materials::get("A5052"), "lsolid2", 0, 0, 0);
526 lsolid2->SetVisAttributes(att("alum"));
527
528 G4Transform3D tsolid2_p1(t * lt);
529 // string pname("psolid2_p"); pname += to_string(++sol2count);
530 acs->AddPlacedVolume(lsolid2, tsolid2_p1);
531 // new G4PVPlacement(tsolid2_p1, lsolid2, pname.c_str(), crystalSectorLogical, false, 0, overlaps);
532 };
533
534 G4Point3D aa(-50 + 15, 15, 0), bb(-50 + 15, -15, 0);
535
536 aa = tsolid1_p1 * G4Point3D(-50 + 15, 15, 0);
537 bb = tsolid1_p2 * G4Point3D(-50 + 15, -15, 0);
538 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
539
540 aa = tsolid1_p1 * G4Point3D(50 - 15, 15, 0);
541 bb = tsolid1_p2 * G4Point3D(50 - 15, -15, 0);
542 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
543
544 G4Point3D r0(0, 40 / cos(M_PI / 16), 0);
545 G4Vector3D np(sin(M_PI / 16), cos(M_PI / 16), 0), mid;
546 G4RotateZ3D rb(M_PI / 24);
547 double L, phi_uu;
548
549
550 phi_uu = (tsolid1_p1 * G4Vector3D(1, 0, 0)).angle(G4Vector3D(cos(M_PI / 16), -sin(M_PI / 16), 0));
551 G4Vector3D n = tsolid1_p1 * G4Vector3D(-sin(phi_uu / 2), -cos(phi_uu / 2), 0);
552 double xj = 50 - 15;
553 aa = tsolid1_p1 * G4Point3D(-xj, -15, 0);
554 L = -((aa - r0) * np) / (n * np);
555 place_solid3(L, -phi_uu / 2, tsolid1_p1 * G4Translate3D(-xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
556 place_solid3(L, -phi_uu / 2, tsolid1_p2 * G4Translate3D(-xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
557 aa = tsolid1_p1 * G4Point3D(xj, -15, 0);
558 L = -((aa - r0) * np) / (n * np);
559 place_solid3(L, phi_uu / 2, tsolid1_p1 * G4Translate3D(xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
560 place_solid3(L, phi_uu / 2, tsolid1_p2 * G4Translate3D(xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
561
562 aa = tsolid1_p3 * G4Point3D(-xj, -15, 0);
563 L = -((aa - r0) * np) / (n * np);
564 place_solid3(L, -phi_uu / 2, tsolid1_p3 * G4Translate3D(-xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
565 place_solid3(L, -phi_uu / 2, tsolid1_p4 * G4Translate3D(-xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
566 aa = tsolid1_p3 * G4Point3D(xj, -15, 0);
567 L = -((aa - r0) * np) / (n * np);
568 place_solid3(L, phi_uu / 2, tsolid1_p3 * G4Translate3D(xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
569 place_solid3(L, phi_uu / 2, tsolid1_p4 * G4Translate3D(xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
570
571 aa = tsolid1_p5 * G4Point3D(-xj, -15, 0);
572 L = -((aa - r0) * np) / (n * np);
573 place_solid3(L, -phi_uu / 2, tsolid1_p5 * G4Translate3D(-xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
574 place_solid3(L, -phi_uu / 2, tsolid1_p6 * G4Translate3D(-xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
575 aa = tsolid1_p5 * G4Point3D(xj, -15, 0);
576 L = -((aa - r0) * np) / (n * np);
577 place_solid3(L, phi_uu / 2, tsolid1_p5 * G4Translate3D(xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
578 place_solid3(L, phi_uu / 2, tsolid1_p6 * G4Translate3D(xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
579
580 phi_uu = M_PI / 16;
581 n = tsolid1_p7 * G4Vector3D(-sin(phi_uu / 2), -cos(phi_uu / 2), 0);
582 xj = 80 - 15;
583 aa = tsolid1_p7 * G4Point3D(-xj, -15, 0);
584 L = -((aa - r0) * np) / (n * np);
585 place_solid3(L, -phi_uu / 2, tsolid1_p7 * G4Translate3D(-xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
586 place_solid3(L, -phi_uu / 2, tsolid1_p7 * G4Translate3D(-xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
587 aa = tsolid1_p7 * G4Point3D(xj, -15, 0);
588 L = -((aa - r0) * np) / (n * np);
589 place_solid3(L, phi_uu / 2, tsolid1_p7 * G4Translate3D(xj, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
590 place_solid3(L, phi_uu / 2, tsolid1_p7 * G4Translate3D(xj, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
591
592
593 aa = tsolid1_p3 * G4Point3D(-50 + 15, 15, 0);
594 bb = tsolid1_p4 * G4Point3D(-50 + 15, -15, 0);
595 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
596
597 aa = tsolid1_p3 * G4Point3D(50 - 15, 15, 0);
598 bb = tsolid1_p4 * G4Point3D(50 - 15, -15, 0);
599 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
600
601
602 aa = tsolid1_p5 * G4Point3D(-50 + 15, 15, 0);
603 bb = tsolid1_p6 * G4Point3D(-50 + 15, -15, 0);
604 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
605
606 aa = tsolid1_p5 * G4Point3D(50 - 15, 15, 0);
607 bb = tsolid1_p6 * G4Point3D(50 - 15, -15, 0);
608 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
609
610
611 obj2_dz = Z0 - 100;
612 aa = tsolid1_p8 * G4Point3D(-50 + 15, 15, 0);
613 bb = tsolid1_p9 * G4Point3D(-50 + 15, -15, 0);
614 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 72, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
615
616 aa = tsolid1_p8 * G4Point3D(50 - 15, 15, 0);
617 bb = tsolid1_p9 * G4Point3D(50 - 15, -15, 0);
618 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 72, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
619
620 aa = tsolid1_p9 * G4Point3D(-50 + 15, 15, 0);
621 bb = tsolid1_p10 * G4Point3D(-50 + 15, -15, 0);
622 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 72, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
623
624 aa = tsolid1_p9 * G4Point3D(50 - 15, 15, 0);
625 bb = tsolid1_p10 * G4Point3D(50 - 15, -15, 0);
626 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 72, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
627
628 phi_uu = (tsolid1_p8 * G4Vector3D(1, 0, 0)).angle(tsolid1_p11 * G4Vector3D(1, 0, 0));
629 r0 = tsolid1_p11 * G4Point3D(-50 + 15, 5, 0);
630 np = tsolid1_p11 * G4Vector3D(0, 1, 0);
631
632 n = tsolid1_p8 * G4Vector3D(-sin(phi_uu / 2), -cos(phi_uu / 2), 0);
633 aa = tsolid1_p8 * G4Point3D(-50 + 15, -15, 0);
634 L = -((aa - r0) * np) / (n * np);
635 place_solid3(L, -phi_uu / 2, tsolid1_p8 * G4Translate3D(-50 + 15, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
636
637 aa = tsolid1_p8 * G4Point3D(50 - 15, -15, 0);
638 L = -((aa - r0) * np) / (n * np);
639 place_solid3(L, phi_uu / 2, tsolid1_p8 * G4Translate3D(50 - 15, -15, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
640
641 phi_uu = (tsolid1_p11 * G4Vector3D(1, 0, 0)).angle(G4RotateZ3D(-M_PI / 16) * G4Vector3D(1, 0, 0));
642 r0 = G4RotateZ3D(-M_PI / 16) * G4Point3D(0, 40, 0);
643 np = G4RotateZ3D(-M_PI / 16) * G4Vector3D(0, 1, 0);
644
645 n = tsolid1_p11 * G4Vector3D(-sin(phi_uu / 2), -cos(phi_uu / 2), 0);
646 aa = tsolid1_p11 * G4Point3D(-50 + 15, -5, 0);
647 L = -((aa - r0) * np) / (n * np);
648 place_solid3(L, -phi_uu / 2, tsolid1_p11 * G4Translate3D(-50 + 15, -5, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
649
650 aa = tsolid1_p11 * G4Point3D(50 - 15, -5, 0);
651 L = -((aa - r0) * np) / (n * np);
652 place_solid3(L, phi_uu / 2, tsolid1_p11 * G4Translate3D(50 - 15, -5, 0)*G4RotateZ3D(-phi_uu / 2)*G4Translate3D(0, -L / 2, 0));
653
654 G4Transform3D ttt(G4RotateZ3D(M_PI / 16)*G4Translate3D(Ro - 8 - 50, 0, Z0 - 100));
655 phi_uu = (tsolid1_p10 * G4Vector3D(1, 0, 0)).angle(ttt * G4Vector3D(1, 0, 0));
656 r0 = ttt * G4Point3D(0, -40, 0);
657 np = ttt * G4Vector3D(0, 1, 0);
658 n = tsolid1_p10 * G4Vector3D(-sin(phi_uu / 2), cos(phi_uu / 2), 0);
659 aa = tsolid1_p10 * G4Point3D(-50 + 15, 15, 0);
660 L = -((aa - r0) * np) / (n * np);
661 place_solid3(L, -phi_uu / 2, tsolid1_p10 * G4Translate3D(-50 + 15, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
662
663 aa = tsolid1_p10 * G4Point3D(50 - 15, 15, 0);
664 L = -((aa - r0) * np) / (n * np);
665 place_solid3(L, phi_uu / 2, tsolid1_p10 * G4Translate3D(50 - 15, 15, 0)*G4RotateZ3D(phi_uu / 2)*G4Translate3D(0, L / 2, 0));
666
667 obj2_dz = Z0 - 75;
668 aa = tsolid1_p12 * G4Point3D(-43 + 15, -5, 0);
669 bb = tsolid1_p13 * G4Point3D(-43 + 15, 5, 0);
670 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
671
672 aa = tsolid1_p12 * G4Point3D(43 - 15, -5, 0);
673 bb = tsolid1_p13 * G4Point3D(43 - 15, 5, 0);
674 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
675
676 aa = tsolid1_p12 * G4Point3D(-43 + 15, 5, 0);
677 bb = tsolid1_p14 * G4Point3D(-43 + 15, -5, 0);
678 place_solid2(obj2_dz, (aa - bb).mag(), -M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
679
680 aa = tsolid1_p12 * G4Point3D(43 - 15, 5, 0);
681 bb = tsolid1_p14 * G4Point3D(43 - 15, -5, 0);
682 place_solid2(obj2_dz, (aa - bb).mag(), M_PI / 48, (aa + bb).phi(), (aa + bb).rho() / 2, 0);
683
684 G4VSolid* solid11_p1 = new G4Box("fwd_solid11_p1", 80. / 2, 30. / 2, 40. / 2);
685 G4VSolid* solid11_p2 = new G4Box("fwd_solid11_p2", 81. / 2, 30. / 2, 40. / 2);
686 G4VSolid* solid11_p3 = new G4SubtractionSolid("fwd_solid11_p3", solid11_p1, solid11_p2, G4Translate3D(0, -3, 3));
687
688 G4LogicalVolume* lsolid11 = new G4LogicalVolume(solid11_p3, Materials::get("SUS304"), "lsolid11", 0, 0, 0);
689 lsolid11->SetVisAttributes(att("iron"));
690 G4Transform3D tsolid11_p1(G4RotateZ3D(M_PI / 16)*G4Translate3D(580, -35, Z0 - 40));
691 acs->AddPlacedVolume(lsolid11, tsolid11_p1);
692
693 G4Transform3D tsolid11_p2(G4RotateZ3D(-M_PI / 16)*G4Translate3D(580, 35, Z0 - 40)*G4RotateZ3D(M_PI));
694 acs->AddPlacedVolume(lsolid11, tsolid11_p2);
695
696 G4VSolid* solid12_p1 = new G4Box("fwd_solid12_p1", 120. / 2, 20. / 2, 115. / 2);
697 G4VSolid* solid12_p2 = new G4Box("fwd_solid12_p2", 121. / 2, 12. / 2, 86. / 2);
698 G4VSolid* solid12_p3 = new G4SubtractionSolid("fwd_solid12_p3", solid12_p1, solid12_p2, G4Translate3D(0, -11. / 2 + 1,
699 86. / 2 - 17.5 - 5));
700 G4LogicalVolume* lsolid12 = new G4LogicalVolume(solid12_p3, Materials::get("A6063"), "lsolid12", 0, 0, 0);
701 lsolid12->SetVisAttributes(att("alum"));
702
703 G4Transform3D tsolid12_p2(G4RotateZ3D(-M_PI / 16)*G4Translate3D(Ro - 8 - 60, 30, Z0 - 115. / 2)*G4RotateZ3D(M_PI));
704 acs->AddPlacedVolume(lsolid12, tsolid12_p2);
705
706 G4VSolid* solid12r_p1 = new G4Box("fwd_solid12r_p1", 100. / 2, 30. / 2, 75. / 2);
707 G4VSolid* solid12r_p2 = new G4Box("fwd_solid12r_p2", 101. / 2, 21. / 2, 41. / 2);
708 G4VSolid* solid12r_p3 = new G4SubtractionSolid("fwd_solid12r_p3", solid12r_p1, solid12r_p2, G4Translate3D(0, -21. / 2 + 5,
709 -41. / 2 - 75. / 2 + 40));
710 G4LogicalVolume* lsolid12r = new G4LogicalVolume(solid12r_p3, Materials::get("A6063"), "lsolid12r", 0, 0, 0);
711 lsolid12r->SetVisAttributes(att("alum"));
712 G4Transform3D tsolid12r_p1(G4RotateZ3D(M_PI / 16)*G4Translate3D(Ro - 8 - 50, -30 - 15, Z0 - 40 - 75. / 2));
713 acs->AddPlacedVolume(lsolid12r, tsolid12r_p1);
714
715 G4Transform3D tr = G4RotateZ3D(M_PI / 16) * G4ReflectY3D();
716 acs->MakeImprint(innervolumesector_logical, tr, 0, overlap);
717 tr = G4RotateZ3D(-M_PI / 16);
718 acs->MakeImprint(innervolumesector_logical, tr, 1, overlap);
719
720 delete acs;
721 }
722
723 if (b_support_structure_15) { // numbering scheme as in ECL-004K102.pdf page 15
724 G4AssemblyVolume* acs = new G4AssemblyVolume();
725
726 double Z0 = 434, Ro = 1415 - 20;
727
728 G4VSolid* solid1_p1 = new G4Box("fwd_solid1_p1", 10. / 2, 398. / 2, 5. / 2);
729 G4VSolid* solid1_p2 = new G4Box("fwd_solid1_p2", 4. / 2, 398. / 2 - 32, 6. / 2);
730 G4VSolid* solid1_p3 = new G4SubtractionSolid("fwd_solid1_p3", solid1_p1, solid1_p2, G4Transform3D(G4RotationMatrix(),
731 G4ThreeVector(-4,
732 0, 0)));
733
734 G4LogicalVolume* lsolid1 = new G4LogicalVolume(solid1_p3, Materials::get("SUS304"), "lsolid1", 0, 0, 0);
735 lsolid1->SetVisAttributes(att("iron"));
736 G4Transform3D tsolid1_p1(G4Translate3D(1350, -16, Z0 - 40 + 3 + 2.5));
737 acs->AddPlacedVolume(lsolid1, tsolid1_p1);
738
739 G4VSolid* solid1a_p1 = new G4Box("fwd_solid1a_p1", 10. / 2, 348.5 / 2, 5. / 2);
740 G4VSolid* solid1a_p2 = new G4Box("fwd_solid1a_p2", 4. / 2, 348.5 / 2 - 32, 6. / 2);
741 G4VSolid* solid1a_p3 = new G4SubtractionSolid("fwd_solid1a_p3", solid1a_p1, solid1a_p2, G4Transform3D(G4RotationMatrix(),
742 G4ThreeVector(4, 0, 0)));
743
744 G4LogicalVolume* lsolid1a = new G4LogicalVolume(solid1a_p3, Materials::get("SUS304"), "lsolid1a", 0, 0, 0);
745 lsolid1a->SetVisAttributes(att("iron"));
746 G4Transform3D tsolid1a_p1(G4Translate3D(1250, -16, Z0 - 40 + 3 + 2.5));
747 acs->AddPlacedVolume(lsolid1a, tsolid1a_p1);
748
749 G4VSolid* solid1b_p1 = new G4Box("fwd_solid1b_p1", 10. / 2, (210 + 210 + 16) / 2, 5. / 2);
750 G4VSolid* solid1b_p2 = new G4Box("fwd_solid1b_p2", 4. / 2, (210 + 210 + 16) / 2 - 16, 6. / 2);
751 G4VSolid* solid1b_p3 = new G4SubtractionSolid("fwd_solid1b_p3", solid1b_p1, solid1b_p2, G4Transform3D(G4RotationMatrix(),
752 G4ThreeVector(4, 0, 0)));
753
754 G4LogicalVolume* lsolid1b = new G4LogicalVolume(solid1b_p3, Materials::get("SUS304"), "lsolid1b", 0, 0, 0);
755 lsolid1b->SetVisAttributes(att("iron"));
756 G4Transform3D t1b(G4Translate3D(1204 - (210 + 210 + 16) / 2 + 16, 0, Z0 - 52 + 8 + 5. / 2));
757 G4Transform3D tsolid1b_p1(t1b * G4TranslateY3D(-10)*G4RotateZ3D(-M_PI / 2));
758 acs->AddPlacedVolume(lsolid1b, tsolid1b_p1);
759 G4Transform3D tsolid1b_p2(t1b * G4TranslateY3D(10)*G4RotateZ3D(M_PI / 2));
760 acs->AddPlacedVolume(lsolid1b, tsolid1b_p2);
761 G4Transform3D tsolid1b_p3(t1b * G4TranslateY3D(-110)*G4RotateZ3D(M_PI / 2));
762 acs->AddPlacedVolume(lsolid1b, tsolid1b_p3);
763 G4Transform3D tsolid1b_p4(t1b * G4TranslateY3D(110)*G4RotateZ3D(-M_PI / 2));
764 acs->AddPlacedVolume(lsolid1b, tsolid1b_p4);
765
766 G4VSolid* solid2_p1 = new G4Box("fwd_solid2_p1", 20. / 2, 210. / 2, 5. / 2);
767 G4VSolid* solid2_p2 = new G4Box("fwd_solid2_p2", 14. / 2, 210. / 2 - 16, 6. / 2);
768 G4VSolid* solid2_p3 = new G4SubtractionSolid("fwd_solid2_p3", solid2_p1, solid2_p2, G4Transform3D(G4RotationMatrix(),
769 G4ThreeVector(-4,
770 0, 0)));
771
772 G4LogicalVolume* lsolid2 = new G4LogicalVolume(solid2_p3, Materials::get("SUS304"), "lsolid2", 0, 0, 0);
773 lsolid2->SetVisAttributes(att("iron"));
774 G4Transform3D tsolid2_p1(G4Translate3D(1204 - (210 + 210) - 210 / 2 + 16, 58, Z0 - 52 + 8 + 5. / 2)*G4RotateZ3D(-M_PI / 2));
775 acs->AddPlacedVolume(lsolid2, tsolid2_p1);
776 // new G4PVPlacement(tsolid2_p1, lsolid2, "psolid2_p1", crystalSectorLogical, false, 0, overlap);
777 G4Transform3D tsolid2_p2(G4Translate3D(1204 - (210 + 210) - 210 / 2 + 16, -58, Z0 - 52 + 8 + 5. / 2)*G4RotateZ3D(M_PI / 2));
778 acs->AddPlacedVolume(lsolid2, tsolid2_p2);
779 // new G4PVPlacement(tsolid2_p2, lsolid2, "psolid2_p2", crystalSectorLogical, false, 0, overlap);
780
781 G4VSolid* solid3_p1 = new G4Box("fwd_solid3_p1", 32. / 2, 396. / 2, 8. / 2);
782 G4LogicalVolume* lsolid3 = new G4LogicalVolume(solid3_p1, Materials::get("SUS304"), "lsolid3", 0, 0, 0);
783 lsolid3->SetVisAttributes(att("iron"));
784 G4Transform3D tsolid3_p1(G4Translate3D(1204, 0, Z0 - 52 + 8. / 2));
785 acs->AddPlacedVolume(lsolid3, tsolid3_p1);
786 // new G4PVPlacement(tsolid3_p1, lsolid3, "psolid3_p1", crystalSectorLogical, false, 0, overlap);
787
788 G4VSolid* solid3n_p1 = new G4Box("fwd_solid3n_p1", 32. / 2, 230. / 2, 8. / 2);
789 G4LogicalVolume* lsolid3n = new G4LogicalVolume(solid3n_p1, Materials::get("SUS304"), "lsolid3n", 0, 0, 0);
790 lsolid3n->SetVisAttributes(att("iron"));
791 G4Transform3D tsolid3n_p1(G4Translate3D(1204 - 420, 0, Z0 - 52 + 8. / 2));
792 acs->AddPlacedVolume(lsolid3n, tsolid3n_p1);
793 // new G4PVPlacement(tsolid3n_p1, lsolid3n, "psolid3n_p1", crystalSectorLogical, false, 0, overlap);
794
795 G4VSolid* solid4_p1 = new G4Box("fwd_solid4_p1", 16. / 2, 160. / 2, 8. / 2);
796 G4LogicalVolume* lsolid4 = new G4LogicalVolume(solid4_p1, Materials::get("SUS304"), "lsolid4", 0, 0, 0);
797 lsolid4->SetVisAttributes(att("iron"));
798 G4Transform3D tsolid4_p1(G4Translate3D(598, 0, Z0 - 52 + 8. / 2));
799 acs->AddPlacedVolume(lsolid4, tsolid4_p1);
800 // new G4PVPlacement(tsolid4_p1, lsolid4, "psolid4_p1", crystalSectorLogical, false, 0, overlap);
801
802 G4VSolid* solid5_p1 = new G4Box("fwd_solid5_p1", 650. / 2, 30. / 2, 40. / 2);
803 G4VSolid* solid5_p2 = new G4Box("fwd_solid5_p2", 651. / 2, 30. / 2, 40. / 2);
804 G4VSolid* solid5_p3 = new G4SubtractionSolid("solid5_p3", solid5_p1, solid5_p2, G4Translate3D(0, -3, 3));
805
806 G4LogicalVolume* lsolid5 = new G4LogicalVolume(solid5_p3, Materials::get("SUS304"), "lsolid5", 0, 0, 0);
807 lsolid5->SetVisAttributes(att("iron"));
808 G4Transform3D tsolid5_p1(G4RotateZ3D(M_PI / 16)*G4Translate3D(910, -45, Z0 - 20 - 15));
809 acs->AddPlacedVolume(lsolid5, tsolid5_p1);
810 // new G4PVPlacement(tsolid5_p1, lsolid5, "psolid5_p1", crystalSectorLogical, false, 0, overlap);
811 G4Transform3D tsolid5_p3(G4RotateZ3D(-M_PI / 16)*G4Translate3D(910, 45, Z0 - 20 - 15)*G4RotateZ3D(M_PI));
812 acs->AddPlacedVolume(lsolid5, tsolid5_p3);
813 // new G4PVPlacement(tsolid5_p3, lsolid5, "psolid5_p3", crystalSectorLogical, false, 0, overlap);
814
815 G4VSolid* solid7_p1 = new G4Box("fwd_solid7_p1", 130. / 2, 30. / 2, 40. / 2);
816 G4VSolid* solid7_p2 = new G4Box("fwd_solid7_p2", 131. / 2, 30. / 2, 40. / 2);
817 G4VSolid* solid7_p3 = new G4SubtractionSolid("fwd_solid7_p3", solid7_p1, solid7_p2, G4Transform3D(G4RotationMatrix(),
818 G4ThreeVector(0,
819 -3, 3)));
820 G4LogicalVolume* lsolid7 = new G4LogicalVolume(solid7_p3, Materials::get("SUS304"), "lsolid7", 0, 0, 0);
821 lsolid7->SetVisAttributes(att("iron"));
822 G4Transform3D tsolid7_p1(G4RotateZ3D(-M_PI / 16)*G4Translate3D(Ro - 8 - 65, 54, Z0 - 40. / 2)*G4RotateZ3D(M_PI));
823 acs->AddPlacedVolume(lsolid7, tsolid7_p1);
824 // new G4PVPlacement(tsolid7_p1, lsolid7, "psolid7_p1", crystalSectorLogical, false, 0, overlap);
825 G4Transform3D tsolid7_p2(G4RotateZ3D(M_PI / 16)*G4Translate3D(Ro - 8 - 65 - 8, -85, Z0 - 40. / 2));
826 acs->AddPlacedVolume(lsolid7, tsolid7_p2);
827 // new G4PVPlacement(tsolid7_p2, lsolid7, "psolid7_p2", crystalSectorLogical, false, 0, overlap);
828
829 // G4Transform3D tsolid7_p3(G4RotateZ3D(M_PI/16)*G4TranslateZ3D(1960 + 438 - 40./2 - 1)*G4TranslateX3D(1315)*G4TranslateY3D(-85));
830 // new G4PVPlacement(tsolid7_p3, lsolid7, "psolid7_p3", crystalSectorLogical, false, 0, overlap);
831 // G4Transform3D tsolid7_p4(G4RotateZ3D(M_PI/16)*G4TranslateZ3D(1960 + 438 - 40./2 - 1)*G4TranslateX3D(1315)*G4TranslateY3D(85)*G4RotateZ3D(M_PI));
832 // new G4PVPlacement(tsolid7_p4, lsolid7, "psolid7_p4", crystalSectorLogical, false, 0, overlap);
833
834 G4VSolid* solid8_p1 = new G4Box("fwd_solid8_p1", 120. / 2, 10. / 2, 42. / 2);
835 G4LogicalVolume* lsolid8 = new G4LogicalVolume(solid8_p1, Materials::get("SUS304"), "lsolid8", 0, 0, 0);
836 lsolid8->SetVisAttributes(att("iron"));
837 G4Transform3D tsolid8_p1(G4RotateZ3D(-M_PI / 16)*G4Translate3D(Ro - 8 - 60, 34, Z0 - 42. / 2));
838 acs->AddPlacedVolume(lsolid8, tsolid8_p1);
839 // new G4PVPlacement(tsolid8_p1, lsolid8, "psolid8_p1", crystalSectorLogical, false, 0, overlap);
840 // G4Transform3D tsolid8_p2(G4RotateZ3D(-M_PI/16)*G4TranslateZ3D(1960 + 438 - 42./2 - 1)*G4TranslateX3D(1320)*G4TranslateY3D(34)*G4RotateZ3D(M_PI));
841 // new G4PVPlacement(tsolid8_p2, lsolid8, "psolid8_p2", crystalSectorLogical, false, 0, overlap);
842
843
844 G4VSolid* solid9_p1 = new G4Box("fwd_solid9_p1", 26. / 2, 12. / 2, 26. / 2);
845 G4LogicalVolume* lsolid9 = new G4LogicalVolume(solid9_p1, Materials::get("SUS304"), "lsolid9", 0, 0, 0);
846 lsolid9->SetVisAttributes(att("iron"));
847 G4Transform3D tsolid9_p1(G4RotateZ3D(-M_PI / 16)*G4Translate3D(Ro - 8 - 120 + 26. / 2 + 17, 20 + 9 + 10 + 9 + 3,
848 Z0 - 37 + 26. / 2));
849 acs->AddPlacedVolume(lsolid9, tsolid9_p1);
850 // new G4PVPlacement(tsolid9_p1, lsolid9, "psolid9", crystalSectorLogical, false, 0, overlap);
851
852
853 G4VSolid* solid10_p1 = new G4Box("fwd_solid10_p1", 26. / 2, 12. / 2, 42. / 2);
854 G4LogicalVolume* lsolid10 = new G4LogicalVolume(solid10_p1, Materials::get("SUS304"), "lsolid10", 0, 0, 0);
855 lsolid10->SetVisAttributes(att("iron"));
856 // G4Transform3D tsolid5_p3(G4RotateZ3D(-M_PI/16)*G4Translate3D(910, 45, Z0-20-15)*G4RotateZ3D(M_PI));
857 G4Transform3D tsolid10_p1(G4RotateZ3D(-M_PI / 16)*G4Translate3D(910 + 650. / 2 - 221, 20 + 10 + 6 + 6, Z0 - 52 + 42. / 2));
858 acs->AddPlacedVolume(lsolid10, tsolid10_p1);
859 // new G4PVPlacement(tsolid10_p1, lsolid10, "psolid10", crystalSectorLogical, false, 0, overlap);
860
861 if (b_connectors) {
862 double t = 2, h20 = 32;
863 G4VSolid* solid_connector = new G4Box("fwd_solid_connector", (110 + 2 * 20) / 2, (250 + 2 * 20) / 2, h20 / 2);
864 G4VSolid* solid_connector2 = new G4Box("fwd_solid_connector2", (20 + t) / 2, (20 + t) / 2, 1.01 * h20 / 2);
865 G4VSolid* solid_connector3 = new G4Box("fwd_solid_connector3", 20 / 2, (250 + 2 * 20 + 2) / 2, h20 / 2);
866 solid_connector = new G4SubtractionSolid("fwd_solid_connector", solid_connector, solid_connector2, G4Translate3D(55 + (20 + t) / 2,
867 -140, 0));
868 solid_connector = new G4SubtractionSolid("fwd_solid_connector", solid_connector, solid_connector3, G4Translate3D(-70, 0, -t - 1));
869 G4LogicalVolume* lsolid_connector = new G4LogicalVolume(solid_connector, Materials::get("G4_AIR"), "lsolid_connector", 0,
870 0, 0);
871 lsolid_connector->SetVisAttributes(att("air"));
872 G4Transform3D tsolid_connector(G4Translate3D(1360 - 60, 0, Z0 - h20 / 2));
873 acs->AddPlacedVolume(lsolid_connector, tsolid_connector);
874 // new G4PVPlacement(tsolid_connector, lsolid_connector, "psolid20", crystalSectorLogical, false, 0, overlap);
875
876 auto lvolume = [&](int part, double dx, double dy, double dz) {
877 ostringstream ost(""); ost << "solid20_p" << part;
878 G4VSolid* sv = new G4Box(ost.str().c_str(), dx / 2, dy / 2, dz / 2);
879 ost.str(""); ost << "lsolid20_p" << part;
880 return new G4LogicalVolume(sv, Materials::get("A5052"), ost.str().c_str(), 0, 0, 0);
881 };
882
883 auto place = [&](G4LogicalVolume * lv, const G4Translate3D & move, int n) {
884 lv->SetVisAttributes(att("alum"));
885 ostringstream ost(""); ost << "phys_" << lv->GetName();
886 new G4PVPlacement(move, lv, ost.str().c_str(), lsolid_connector, false, n, overlap);
887 };
888 G4LogicalVolume* lv1 = lvolume(1, 20, 250 + 2 * 20, t);
889 place(lv1, G4Translate3D(-55 - 10, 0, h20 / 2 - t / 2), 0);
890 G4LogicalVolume* lv1_2 = lvolume(1, 20, 250 + 20, t);
891 place(lv1_2, G4Translate3D(55 + 10, 10, h20 / 2 - t / 2), 1);
892
893 G4LogicalVolume* lv2 = lvolume(2, 110., 20., t);
894 place(lv2, G4Translate3D(0, 250 / 2 + 10, h20 / 2 - t / 2), 0);
895 place(lv2, G4Translate3D(0, -250 / 2 - 10, h20 / 2 - t / 2), 1);
896
897 G4LogicalVolume* lv3 = lvolume(3, 110., t, h20 - t);
898 place(lv3, G4Translate3D(0, 250 / 2 + t / 2, -t / 2), 0);
899 place(lv3, G4Translate3D(0, -250 / 2 - t / 2, -t / 2), 1);
900
901 G4LogicalVolume* lv4 = lvolume(4, t, 250 + 2 * t, h20 - t);
902 place(lv4, G4Translate3D(55 + t / 2, 0, -t / 2), 0);
903 place(lv4, G4Translate3D(-55 - t / 2, 0, -t / 2), 1);
904
905 G4LogicalVolume* lv5 = lvolume(5, 7, 250, t);
906 place(lv5, G4Translate3D(55 - 7. / 2, 0, -h20 / 2 + t / 2 + t), 0);
907 place(lv5, G4Translate3D(-55 + 7. / 2, 0, -h20 / 2 + t / 2 + t), 1);
908
909 G4LogicalVolume* lv6 = lvolume(6, 110 - 14, 7, t);
910 place(lv6, G4Translate3D(0, 250 / 2 - 7. / 2, -h20 / 2 + t / 2 + t), 0);
911 place(lv6, G4Translate3D(0, -250 / 2 + 7. / 2, -h20 / 2 + t / 2 + t), 1);
912
913 G4LogicalVolume* lv7 = lvolume(7, 110, 250, t);
914 place(lv7, G4Translate3D(0, 0, -h20 / 2 + t / 2), 0);
915
916 // G4LogicalVolume *lv8 = lvolume(8, 90, 10, 30);
917 G4VSolid* p8_1 = new G4Box("fwd_solid20_p8_1", 90. / 2, 10. / 2, 30. / 2);
918 G4VSolid* p8_2 = new G4Box("fwd_solid20_p8_2", 88. / 2, 8. / 2, 30. / 2);
919 G4VSolid* sp8 = new G4SubtractionSolid("fwd_solid20_p8", p8_1, p8_2, G4TranslateZ3D(-1));
920 G4LogicalVolume* lv8 = new G4LogicalVolume(sp8, Materials::get("A5052"), "lsolid20_p8", 0, 0, 0);
921 lv8->SetVisAttributes(att("alum2"));
922 for (int i = 0; i < 10; i++) place(lv8, G4Translate3D(0, 25 * (i - 4.5), -h20 / 2 + t + 30 / 2), i);
923 }
924
925 if (b_boards) {
926 double hbv = 30;
927 G4VSolid* solid_board = new G4Box("fwd_solid_board", (210) / 2, (110) / 2, hbv / 2);
928 G4LogicalVolume* lsolid_board = new G4LogicalVolume(solid_board, Materials::get("G4_AIR"), "lsolid_board", 0, 0, 0);
929 lsolid_board->SetVisAttributes(att("air"));
930 for (int i = 0; i < 1; i++) {
931 // G4Transform3D t0 = G4RotateZ3D(M_PI/16*(1-2*i))*G4Translate3D(598-8+210/2, 0, Z0-52+8+5+hbv/2);
932 G4Transform3D t0 = G4Translate3D(598 - 8 + 210 / 2, 0, Z0 - 52 + 8 + 5 + hbv / 2);
933 G4Transform3D t1 = t0 * G4Translate3D(210, 110 / 2 + 5, 0);
934 G4Transform3D t2 = t0 * G4Translate3D(210, -110 / 2 + 5, 0);
935 G4Transform3D t3 = t0 * G4Translate3D(2 * 210, 110 / 2 + 5, 0);
936 G4Transform3D t4 = t0 * G4Translate3D(2 * 210, -110 / 2 + 5, 0);
937 // new G4PVPlacement(t0, lsolid_board, "phys_solid_board0", crystalSectorLogical, false, 0, overlap);
938 // new G4PVPlacement(t1, lsolid_board, "phys_solid_board1", crystalSectorLogical, false, 1, overlap);
939 // new G4PVPlacement(t2, lsolid_board, "phys_solid_board2", crystalSectorLogical, false, 2, overlap);
940 // new G4PVPlacement(t3, lsolid_board, "phys_solid_board3", crystalSectorLogical, false, 3, overlap);
941 // new G4PVPlacement(t4, lsolid_board, "phys_solid_board4", crystalSectorLogical, false, 4, overlap);
942 acs->AddPlacedVolume(lsolid_board, t0);
943 acs->AddPlacedVolume(lsolid_board, t1);
944 acs->AddPlacedVolume(lsolid_board, t2);
945 acs->AddPlacedVolume(lsolid_board, t3);
946 acs->AddPlacedVolume(lsolid_board, t4);
947 }
948 G4Material* boxmaterial = Materials::get("G4_GLASS_PLATE");
949 auto lvolumeb = [&](int part, double dx, double dy, double dz) {
950 ostringstream ost(""); ost << "fwd_sboard_p" << part;
951 G4VSolid* sv = new G4Box(ost.str().c_str(), dx / 2, dy / 2, dz / 2);
952 ost.str(""); ost << "lboard_p" << part;
953 return new G4LogicalVolume(sv, boxmaterial, ost.str().c_str(), 0, 0, 0);
954 };
955
956 const G4VisAttributes* asolid20 = att("plate");
957
958 auto placeb = [&](G4LogicalVolume * lv, const G4Translate3D & move, int n) {
959 lv->SetVisAttributes(asolid20);
960 ostringstream ost(""); ost << "phys_" << lv->GetName();
961 new G4PVPlacement(move, lv, ost.str().c_str(), lsolid_board, false, n, overlap);
962 };
963
964 double hboard = 2;
965 G4LogicalVolume* lb1 = lvolumeb(1, 210, 110, hboard);
966 placeb(lb1, G4Translate3D(0, 0, -hbv / 2 + hboard / 2), 0);
967
968 double wcon = 20, hcon = 40, hc = hbv - hboard;
969 G4VSolid* sv_connector_bundle = new G4Box("fwd_sv_connector_bundle", 4 * wcon / 2, hcon / 2, hc / 2);
970 G4LogicalVolume* lv_connector_bundle = new G4LogicalVolume(sv_connector_bundle, Materials::get("G4_AIR"),
971 "lv_connector_bundle", 0, 0, 0);
972 lv_connector_bundle->SetVisAttributes(att("air"));
973 new G4PVPlacement(G4Translate3D(-210 / 2 + 10 + wcon * 2, -110 / 2 + hcon / 2, -hbv / 2 + hboard + hc / 2), lv_connector_bundle,
974 "pv_connector_bundle", lsolid_board, false, 0, overlap);
975 new G4PVPlacement(G4Translate3D(-210 / 2 + 10 + wcon * 2, 10 + hcon / 2, -hbv / 2 + hboard + hc / 2), lv_connector_bundle,
976 "pv_connector_bundle", lsolid_board, false, 0, overlap);
977 new G4PVPlacement(G4Translate3D(10 + wcon * 2, -110 / 2 + hcon / 2, -hbv / 2 + hboard + hc / 2), lv_connector_bundle,
978 "pv_connector_bundle", lsolid_board, false, 0, overlap);
979 new G4PVPlacement(G4Translate3D(10 + wcon * 2, 10 + hcon / 2, -hbv / 2 + hboard + hc / 2), lv_connector_bundle,
980 "pv_connector_bundle", lsolid_board, false, 0, overlap);
981
982 G4VSolid* sv_crystal_connector = new G4Box("fwd_sv_crystal_connector", wcon / 2, hcon / 2, hc / 2);
983 G4LogicalVolume* lv_crystal_connector = new G4LogicalVolume(sv_crystal_connector, Materials::get("G4_AIR"),
984 "lv_crystal_connector", 0, 0, 0);
985 lv_crystal_connector->SetVisAttributes(att("air"));
986
987 new G4PVPlacement(G4Translate3D(-1.5 * 20, 0, 0), lv_crystal_connector, "pv_crystal_connector", lv_connector_bundle, false, 0,
988 overlap);
989 new G4PVPlacement(G4Translate3D(-0.5 * 20, 0, 0), lv_crystal_connector, "pv_crystal_connector", lv_connector_bundle, false, 1,
990 overlap);
991 new G4PVPlacement(G4Translate3D(0.5 * 20, 0, 0), lv_crystal_connector, "pv_crystal_connector", lv_connector_bundle, false, 2,
992 overlap);
993 new G4PVPlacement(G4Translate3D(1.5 * 20, 0, 0), lv_crystal_connector, "pv_crystal_connector", lv_connector_bundle, false, 3,
994 overlap);
995
996 G4VSolid* sv_crystal_connector_p1 = new G4Box("fwd_sv_crystal_connector_p1", 8 / 2, 30 / 2, 20 / 2);
997 G4LogicalVolume* lv_crystal_connector_p1 = new G4LogicalVolume(sv_crystal_connector_p1, Materials::get("SUS304"),
998 "lv_crystal_connector_p1", 0, 0, 0);
999 lv_crystal_connector_p1->SetVisAttributes(att("connector"));
1000
1001 new G4PVPlacement(G4Translate3D(-5, 0, -hc / 2 + 20. / 2), lv_crystal_connector_p1, "pv_crystal_connector_p1", lv_crystal_connector,
1002 false, 0, overlap);
1003
1004 G4VSolid* sv_capacitor = new G4Tubs("fwd_sv_capacitor", 0, 5, 5. / 2, 0, 2 * M_PI);
1005 G4LogicalVolume* lv_capacitor = new G4LogicalVolume(sv_capacitor, Materials::get("SUS304"), "lv_capacitor", 0, 0, 0);
1006 lv_capacitor->SetVisAttributes(att("capacitor"));
1007
1008 new G4PVPlacement(G4Translate3D(5, -15, -hc / 2 + 5. / 2), lv_capacitor, "pv_capacitor", lv_crystal_connector, false, 0, overlap);
1009 new G4PVPlacement(G4Translate3D(5, -5, -hc / 2 + 5. / 2), lv_capacitor, "pv_capacitor", lv_crystal_connector, false, 1, overlap);
1010 new G4PVPlacement(G4Translate3D(5, 5, -hc / 2 + 5. / 2), lv_capacitor, "pv_capacitor", lv_crystal_connector, false, 2, overlap);
1011 new G4PVPlacement(G4Translate3D(5, 15, -hc / 2 + 5. / 2), lv_capacitor, "pv_capacitor", lv_crystal_connector, false, 3, overlap);
1012
1013 G4VSolid* sv_board_connector_p1 = new G4Box("fwd_sv_board_connector_p1", 80. / 2, 8. / 2, 20. / 2);
1014 G4LogicalVolume* lv_board_connector_p1 = new G4LogicalVolume(sv_board_connector_p1, Materials::get("SUS304"),
1015 "lv_board_connector_p1", 0, 0, 0);
1016 lv_board_connector_p1->SetVisAttributes(att("connector"));
1017
1018 new G4PVPlacement(G4Translate3D(-210 / 2 + 10 + 80 / 2, 0, -hbv / 2 + hboard + 20. / 2), lv_board_connector_p1,
1019 "pv_board_connector_p1", lsolid_board, false, 0, overlap);
1020 new G4PVPlacement(G4Translate3D(210 / 2 - 10 - 80 / 2, 0, -hbv / 2 + hboard + 20. / 2), lv_board_connector_p1,
1021 "pv_board_connector_p1", lsolid_board, false, 1, overlap);
1022 }
1023
1024 G4Transform3D tr = G4RotateZ3D(M_PI / 16) * G4ReflectY3D();
1025 acs->MakeImprint(innervolumesector_logical, tr, 0, overlap);
1026 tr = G4RotateZ3D(-M_PI / 16);
1027 acs->MakeImprint(innervolumesector_logical, tr, 1, overlap);
1028 delete acs;
1029 }// end of ECL-004K102.pdf page 15
1030
1031 if (b_cover) {
1032 G4VSolid* solid8_p1 = new G4Tubs("fwd_solid8_p1", RI + tand(13.12) * (434 + 1) - 20 / cosd(13.12), 1415, 1. / 2, -M_PI / 16,
1033 M_PI / 8);
1034 G4VSolid* solid8_p2 = new G4Box("fwd_solid8_p2", 130. / 2, 270. / 2, 2. / 2);
1035 G4VSolid* solid8_p3 = new G4Tubs("fwd_solid8_p3", 0, 16, 2, 0, 2 * M_PI);
1036 G4VSolid* solid8_p4 = new G4Box("fwd_solid8_p4", 130. / 2, (75. - 2 * 16.) / 2, 2. / 2);
1037 double width_p5 = 180;
1038 G4VSolid* solid8_p5 = new G4Box("fwd_solid8_p5", width_p5 / 2, 2.5 + 75. / 2, 2. / 2);
1039
1040 double xx0 = 1415 - 47.8715;
1041 G4VSolid* solid8 = new G4SubtractionSolid("fwd_solid8", solid8_p1, solid8_p2, G4TranslateX3D(xx0 - 130. / 2));
1042 double xx1 = xx0 + 1.7;
1043 solid8 = new G4SubtractionSolid("fwd_solid8", solid8, solid8_p3, G4Translate3D(xx1, 159.5, 0));
1044 solid8 = new G4SubtractionSolid("fwd_solid8", solid8, solid8_p3, G4Translate3D(xx1, 202.5, 0));
1045 solid8 = new G4SubtractionSolid("fwd_solid8", solid8, solid8_p4, G4Translate3D(xx1 - 16 + 130. / 2, (202.5 + 159.5) / 2, 0));
1046 solid8 = new G4SubtractionSolid("fwd_solid8", solid8, solid8_p5, G4Translate3D(xx1 + width_p5 / 2, (202.5 + 159.5) / 2, 0));
1047 solid8 = new G4SubtractionSolid("fwd_solid8", solid8, solid8_p5, G4Translate3D(xx1 - 130 + (1230.77 - 1230.88), -177.57,
1048 0)*G4RotateZ3D(-M_PI / 16)*G4Translate3D(width_p5 / 2, -75. / 2, 0));
1049
1050 // for(int i=0;i<100000;i++){
1051 // G4ThreeVector v = solid8->GetPointOnSurface();
1052 // G4cout<<v.x()<<" "<<v.y()<<" "<<v.z()<<"\n";
1053 // }
1054
1055 G4LogicalVolume* lsolid8 = new G4LogicalVolume(solid8, Materials::get("A5052"), "lsolid8", 0, 0, 0);
1056 lsolid8->SetVisAttributes(att("alum"));
1057 for (int i = 0; i < 8; i++) {
1058 G4Transform3D tc = gTrans * G4Translate3D(0, 0, 3 + 434 + 0.5) * G4RotateZ3D(M_PI / 8 + i * M_PI / 4);
1059 new G4PVPlacement(tc * G4RotateZ3D(M_PI / 16), lsolid8, suf("ECL_Forward_cover", 0 + 2 * i), top, false, 0 + 2 * i, overlap);
1060 G4ReflectionFactory::Instance()->Place(tc * G4RotateZ3D(-M_PI / 16)*G4ReflectY3D(), suf("ECL_Forward_cover", 0 + 2 * i), lsolid8,
1061 top, false,
1062 1 + 2 * i, overlap);
1063 }
1064 }
1065 // end of ECL-004K102.pdf page 11 - 12
1066
1067}

◆ get_pa_box_height()

double get_pa_box_height ( ) const
inlineprivate

Getter for preamplifier box height (hard-coded to be 2)

Definition at line 88 of file GeoECLCreator.h.

88{return 2;}

◆ get_preamp()

G4LogicalVolume * get_preamp ( ) const
private

Get Logical volume of preamplifier.

Definition at line 148 of file GeoECLCreator.cc.

149{
150 static G4LogicalVolume* lv_preamplifier = nullptr;
151 if (lv_preamplifier == nullptr) {
152 G4VSolid* sv_preamplifier = new G4Box("sv_preamplifier", 58. / 2, 51. / 2, get_pa_box_height() / 2);
153 lv_preamplifier = new G4LogicalVolume(sv_preamplifier, Materials::get("A5052"), "lv_preamplifier", 0, 0, 0);
154 G4VSolid* sv_diode = new G4Box("sv_diode", 20. / 2, 20. / 2, 0.3 / 2);
155 G4LogicalVolume* lv_diode = new G4LogicalVolume(sv_diode, Materials::get("G4_Si"), "lv_diode", 0, 0, 0);
156 lv_diode->SetUserLimits(new G4UserLimits(0.01));
157 lv_diode->SetSensitiveDetector(m_sensediode);
158 new G4PVPlacement(G4TranslateZ3D(-get_pa_box_height() / 2 + 0.3 / 2), lv_diode, "pv_diode", lv_preamplifier, false, 0, m_overlap);
159
160 // G4VSolid* sv_diode = new G4Box("sv_diode", 10. / 2, 20. / 2, 0.3 / 2);
161 // G4LogicalVolume* lv_diode = new G4LogicalVolume(sv_diode, Materials::get("G4_Si"), "lv_diode", 0, 0, 0);
162 // lv_diode->SetUserLimits(new G4UserLimits(0.01));
163 // new G4PVPlacement(G4Translate3D(-5, 0, -pa_box_height / 2 + 0.3 / 2), lv_diode, "pv_diode1", lv_preamplifier, false, 1, overlap);
164 // new G4PVPlacement(G4Translate3D(5, 0, -pa_box_height / 2 + 0.3 / 2), lv_diode, "pv_diode2", lv_preamplifier, false, 2, overlap);
165
166 lv_preamplifier->SetVisAttributes(att("preamp"));
167 }
168 return lv_preamplifier;
169}

◆ wrapped_crystal()

G4LogicalVolume * wrapped_crystal ( const shape_t * s,
const std::string & endcap,
double wrapthickness )
private

Wrapped crystal.

Definition at line 90 of file GeoECLCreator.cc.

91{
92 std::string prefix("sv_"); prefix += endcap; prefix += "_wrap";
93 G4Translate3D tw;
94 G4VSolid* wrapped_crystal = s->get_solid(prefix, wrapthickness, tw);
95 std::string name("lv_"); name += endcap + "_wrap_" + std::to_string(s->nshape);
96 G4Material* wrap = nullptr;
97 if (wrapthickness < 0.170)
98 wrap = Materials::get("WRAP170");
99 else if (wrapthickness < 0.200)
100 wrap = Materials::get("WRAP200");
101 else
102 wrap = Materials::get("WRAP250");
103 G4LogicalVolume* wrapped_logical = new G4LogicalVolume(wrapped_crystal, wrap, name.c_str(), 0, 0, 0);
104 wrapped_logical->SetVisAttributes(att("wrap"));
105
106 prefix = "sv_"; prefix += endcap; prefix += "_crystal";
107 G4Translate3D tc;
108 G4VSolid* crystal_solid = s->get_solid(prefix, 0, tc);
109 name = "lv_" + endcap + "_crystal_" + std::to_string(s->nshape);
110 G4LogicalVolume* crystal_logical = new G4LogicalVolume(crystal_solid, Materials::get("G4_CESIUM_IODIDE"), name.c_str(),
111 0, 0, 0);
112 crystal_logical->SetVisAttributes(att("cryst"));
113 crystal_logical->SetSensitiveDetector(m_sensitive);
114
115 new G4PVPlacement(nullptr, G4ThreeVector(), crystal_logical, name.c_str(), wrapped_logical, false, 0, 0);
116 return wrapped_logical;
117}

Member Data Documentation

◆ m_atts

std::map<std::string, G4VisAttributes*> m_atts
private

Vector of background-Sensitive detectors.

Definition at line 98 of file GeoECLCreator.h.

◆ m_overlap

int m_overlap
private

overlap

Definition at line 100 of file GeoECLCreator.h.

◆ m_sap

const ECLCrystalsShapeAndPosition* m_sap
private

pointer to a storage with crystal shapes and positions

Definition at line 91 of file GeoECLCreator.h.

◆ m_sensediode

Simulation::SensitiveDetectorBase* m_sensediode
private

Sensitive diode.

Definition at line 96 of file GeoECLCreator.h.

◆ m_sensitive

Simulation::SensitiveDetectorBase* m_sensitive
private

Sensitive detector.

Definition at line 94 of file GeoECLCreator.h.

The documentation for this class was generated from the following files: