SHiP Geometry
SHiP experiment geometry implementation using GeoModel.
Loading...
Searching...
No Matches
TargetFactory.h
Go to the documentation of this file.
1// SPDX-License-Identifier: LGPL-3.0-or-later
2// Copyright (C) CERN for the benefit of the SHiP Collaboration
3
4#pragma once
5
6#include <GeoModelKernel/Units.h>
7
8#include <array>
9
10class GeoPhysVol;
11class GeoLogVol;
12class GeoShape;
13
14namespace SHiPGeometry {
15
16class SHiPMaterials;
17
43class TargetFactory {
44 public:
45 explicit TargetFactory(SHiPMaterials& materials);
46 ~TargetFactory() = default;
47
52 [[nodiscard]] GeoPhysVol* build();
53
54 private:
55 SHiPMaterials& m_materials;
56
57 // Helper methods
58 GeoPhysVol* createProximityShielding();
59 GeoPhysVol* createTopShielding();
60 GeoPhysVol* createBottomShielding();
61 GeoPhysVol* createShieldingPedestal();
62 GeoPhysVol* createHeVolume();
63 const GeoShape* createSteelCoreShape();
64
65 // Unit conversion helpers
66 static constexpr double cm = GeoModelKernelUnits::cm;
67 static constexpr double mm = GeoModelKernelUnits::mm;
68 static constexpr double deg = GeoModelKernelUnits::degree;
69
70 // Target vacuum box dimensions (half-sizes)
71 static constexpr double s_vacuumBoxHalfX = 80.0 * cm;
72 static constexpr double s_vacuumBoxHalfY = 113.55 * cm;
73 static constexpr double s_vacuumBoxHalfZ = 150.0 * cm;
74
75 // Proximity shielding
76 static constexpr double s_proxEnvHalfX = 80.0 * cm;
77 static constexpr double s_proxEnvHalfY = 56.3 * cm;
78 static constexpr double s_proxEnvHalfZ = 150.0 * cm;
79 static constexpr double s_proxInnerHalfX = 55.0 * cm;
80 static constexpr double s_proxInnerHalfY = 56.3 * cm;
81 static constexpr double s_proxInnerHalfZ = 110.0 * cm;
82 static constexpr double s_proxInnerOffsetZ = 15.0 * cm;
83 static constexpr double s_proxHoleRadius = 10.0 * cm;
84 static constexpr double s_proxHoleHalfZ = 27.5 * cm;
85 static constexpr double s_proxHoleOffsetY = 17.2 * cm;
86 static constexpr double s_proxHoleOffsetZ = -122.5 * cm;
87 static constexpr double s_proxPosY = -2.75 * cm;
88
89 // Top shielding
90 static constexpr double s_topShieldHalfX = 80.0 * cm;
91 static constexpr double s_topShieldHalfY = 30.0 * cm;
92 static constexpr double s_topShieldHalfZ = 150.0 * cm;
93 static constexpr double s_topShieldPosY = 83.55 * cm;
94
95 // Bottom shielding
96 static constexpr double s_bottomShieldHalfX = 80.0 * cm;
97 static constexpr double s_bottomShieldHalfY = 27.25 * cm;
98 static constexpr double s_bottomShieldHalfZ = 150.0 * cm;
99 static constexpr double s_bottomShieldPosY = -86.3 * cm;
100
101 // Shielding pedestal
102 static constexpr double s_pedestalHalfX = 53.5 * cm;
103 static constexpr double s_pedestalHalfY = 7.5 * cm;
104 static constexpr double s_pedestalHalfZ = 108.5 * cm;
105 static constexpr double s_pedestalPosY = -51.55 * cm;
106 static constexpr double s_pedestalPosZ = 15.0 * cm;
107
108 // TargetArea position within vacuum box: the target frame (z = 0 at the
109 // front face of the first disk) sits at this offset in the vacuum box
110 static constexpr double s_targetAreaPosY = 14.45 * cm;
111 static constexpr double s_targetAreaPosZ = -43.25 * cm;
112
113 // ---- 2026 BDF target (CATIA ST1A07710_01_AB.02) ----
114 // All z values below are in the target frame (z = 0 at disk-1 front face).
115
116 // HeVolume: container for disks, steel core, jacket, flanges, beam
117 // window, cover-plate bore rings and endcap
118 static constexpr double s_heRadius = 237.0 * mm;
119 static constexpr double s_heZMin = -37.8 * mm; // upstream face of the cover plate
120 static constexpr double s_heZMax = 1509.7 * mm; // downstream end of rear endcap
121 static constexpr double s_heCentreZ = 0.5 * (s_heZMin + s_heZMax);
122
123 // Tungsten disks: [z start, z end] per disk; 4-5 mm He slits in between
124 static constexpr int s_numDisks = 33;
125 static constexpr double s_diskRadius = 125.0 * mm;
126 static constexpr double s_lastDiskRadius = 157.0 * mm; // rear block
127 static constexpr std::array<std::array<double, 2>, s_numDisks> s_diskZ = {{
128 {0.0 * mm, 30.0 * mm}, {34.0 * mm, 44.0 * mm}, {48.0 * mm, 56.0 * mm},
129 {60.0 * mm, 68.0 * mm}, {72.0 * mm, 80.0 * mm}, {84.0 * mm, 92.0 * mm},
130 {96.0 * mm, 104.0 * mm}, {108.0 * mm, 116.0 * mm}, {120.0 * mm, 128.0 * mm},
131 {132.0 * mm, 140.0 * mm}, {144.0 * mm, 152.0 * mm}, {156.0 * mm, 164.0 * mm},
132 {168.0 * mm, 176.0 * mm}, {180.0 * mm, 188.0 * mm}, {192.0 * mm, 200.0 * mm},
133 {204.0 * mm, 212.0 * mm}, {216.0 * mm, 224.0 * mm}, {228.0 * mm, 237.0 * mm},
134 {241.0 * mm, 252.0 * mm}, {256.0 * mm, 268.0 * mm}, {272.0 * mm, 285.0 * mm},
135 {289.0 * mm, 302.0 * mm}, {306.0 * mm, 322.0 * mm}, {326.0 * mm, 346.0 * mm},
136 {350.0 * mm, 370.0 * mm}, {375.0 * mm, 396.0 * mm}, {401.0 * mm, 431.0 * mm},
137 {435.0 * mm, 475.0 * mm}, {479.0 * mm, 534.0 * mm}, {538.0 * mm, 611.0 * mm},
138 {615.0 * mm, 721.0 * mm}, {725.0 * mm, 1001.0 * mm}, {1005.0 * mm, 1460.0 * mm},
139 }};
140
141 // Steel core (two clamp halves modelled as one axisymmetric piece)
142 static constexpr double s_coreZMin = 7.2 * mm;
143 static constexpr double s_coreZMax = 1443.7 * mm;
144 static constexpr double s_coreBoreR1 = 125.0 * mm; // matches disks 1..32
145 static constexpr double s_coreBoreR2 = 157.0 * mm; // matches the rear block
146 static constexpr double s_coreBoreStepZ = 1005.0 * mm; // bore radius change
147 static constexpr double s_coreOuterR = 207.0 * mm;
148 static constexpr double s_coreFrontOuterR = 195.0 * mm; // front, inside the flange
149 static constexpr double s_coreFrontZMax = 55.2 * mm; // end of front step
150 static constexpr double s_coreRearOuterR = 190.0 * mm; // rear, inside rear flange
151 static constexpr double s_coreRearZMin = 1233.7 * mm; // start of rear step
152
153 // Serpentine He cooling grooves: arcs in the bore, centred on the vertical
154 // axis, staggered between the upper and lower half; He flows between them
155 // through the inter-disk slits
156 static constexpr double s_grooveRmin = 120.0 * mm; // below bore: clean subtraction
157 static constexpr double s_grooveRmax = 153.0 * mm;
158 static constexpr double s_groovePhiWidth = 61.0 * deg;
159 static constexpr std::array<std::array<double, 2>, 5> s_groovesTop = {{
160 {12.2 * mm, 57.0 * mm},
161 {95.0 * mm, 153.0 * mm},
162 {191.0 * mm, 257.0 * mm},
163 {300.0 * mm, 464.0 * mm},
164 {559.0 * mm, 933.2 * mm},
165 }};
166 static constexpr std::array<std::array<double, 2>, 4> s_groovesBottom = {{
167 {45.0 * mm, 104.7 * mm},
168 {142.7 * mm, 200.7 * mm},
169 {238.7 * mm, 334.7 * mm},
170 {385.7 * mm, 1005.0 * mm},
171 }};
172 // Groove around the rear block (upper half only), ends 2 mm before the
173 // core rear face
174 static constexpr double s_rearGrooveRmin = 152.0 * mm; // below bore R2
175 static constexpr double s_rearGrooveRmax = 182.0 * mm;
176 static constexpr double s_rearGroovePhiWidth = 49.0 * deg;
177 static constexpr double s_rearGrooveZMin = 1005.0 * mm;
178 static constexpr double s_rearGrooveZMax = 1441.7 * mm;
179
180 // Jacket tube and flanges; He annulus between core (207) and jacket (217).
181 // The flange bore matches the core front step, so the flange/jacket
182 // boundary is placed at the end of the step. The flange proper starts at
183 // the cover-plate step; upstream of it the outer radius is the r 195
184 // nose that slides into the cover-plate bore.
185 static constexpr double s_jacketRmin = 217.0 * mm;
186 static constexpr double s_jacketRmax = 225.0 * mm;
187 static constexpr double s_jacketZMin = s_coreFrontZMax;
188 static constexpr double s_jacketZMax = 1198.7 * mm;
189 static constexpr double s_flangeFrontRmin = 195.0 * mm;
190 static constexpr double s_flangeFrontZMin = -23.8 * mm; // cover-plate step
191 static constexpr double s_noseZMin = -35.8 * mm; // upstream face of the vessel
192 static constexpr double s_flangeRearRmax = 237.0 * mm;
193 static constexpr double s_flangeRearZMax = 1263.7 * mm;
194
195 // Beam window: dished membrane closing the bore, simplified to a flat
196 // disc preserving the 8 mm of steel on the beam axis; the r 141-195 nose
197 // ring connects the window rim to the flange
198 static constexpr double s_windowRmax = 141.0 * mm;
199 static constexpr double s_windowZMin = -33.5 * mm;
200 static constexpr double s_windowZMax = -25.5 * mm;
201
202 // Upstream cover plate (600 x 650 x 20 mm, stepped bore r 196/226 mm,
203 // asymmetric about the beam axis: +250 mm above, -400 mm below): the
204 // part within the He container is modelled as two bore rings, the
205 // rectangular remainder as a box-minus-tube outside it
206 static constexpr double s_coverRing1Rmin = 196.0 * mm;
207 static constexpr double s_coverRing2Rmin = 226.0 * mm;
208 static constexpr double s_coverZMax = -17.8 * mm;
209 static constexpr double s_coverPlateHalfX = 300.0 * mm;
210 static constexpr double s_coverPlateHalfY = 325.0 * mm;
211 static constexpr double s_coverPlateHalfZ = 10.0 * mm;
212 static constexpr double s_coverPlateOffsetY = -75.0 * mm; // plate centre below the axis
213
214 // Rear endcap: 8 mm shell, cylindrical section closed by a domed head;
215 // polycone planes (z, rmin, rmax) approximating the dome with three cone
216 // segments and an apex disc
217 static constexpr std::array<std::array<double, 3>, 6> s_endcapPlanes = {{
218 {1263.7 * mm, 229.0 * mm, 237.0 * mm},
219 {1413.7 * mm, 229.0 * mm, 237.0 * mm},
220 {1450.0 * mm, 212.0 * mm, 223.0 * mm},
221 {1480.0 * mm, 142.0 * mm, 166.5 * mm},
222 {1505.0 * mm, 0.0 * mm, 67.0 * mm},
223 {1509.7 * mm, 0.0 * mm, 67.0 * mm},
224 }};
225};
226
227} // namespace SHiPGeometry
SHiPGeometry::SHiPMaterials
Central material manager for the SHiP detector.
Definition SHiPMaterials.h:20
SHiPGeometry::TargetFactory
Factory for the Target (proton target and shielding) geometry.
Definition TargetFactory.h:43
SHiPGeometry::TargetFactory::build
GeoPhysVol * build()
Build the Target geometry.
Definition TargetFactory.cpp:27
SHiPGeometry
Definition ConfigPath.h:9
Generated by doxygen 1.9.8