68 G4VisAttributes* FoamVisAtt =
new G4VisAttributes(G4Colour::Blue());
69 G4VisAttributes* KaptonVisAtt =
new G4VisAttributes(G4Colour::Red());
70 G4VisAttributes* GlueVisAtt =
new G4VisAttributes(G4Colour::White());
71 G4VisAttributes* MetalVisAtt =
new G4VisAttributes(G4Colour::Gray());
72 G4VisAttributes* SubstrateVisAtt =
new G4VisAttributes(G4Colour::Green());
73 G4VisAttributes* MetalizedVisAtt =
new G4VisAttributes(G4Colour::Cyan());
74 G4VisAttributes* EpitaxialVisAtt =
new G4VisAttributes(G4Colour(1.0, 1.0, 0.0));
77 double stepSize = content.getLength(
"stepSize", 5 * CLHEP::um);
81 GearDir activeParams(content,
"Ladder");
82 G4double SensorLengthX = (activeParams.getInt(
"nPixelsX") * activeParams.getLength(
"pitchX")) * CLHEP::cm;
83 G4double SensorLengthY = (activeParams.getInt(
"nPixelsY") * activeParams.getLength(
"pitchY")) * CLHEP::cm;
85 G4double AirGap = activeParams.getLength(
"AirGap") * CLHEP::cm;
87 G4double SubstrateThickness = activeParams.getLength(
"SubstrateThickness") * CLHEP::cm;
88 G4double EpitaxialThickness = activeParams.getLength(
"EpitaxialThickness") * CLHEP::cm;
89 G4double MetalThickness = activeParams.getLength(
"MetalThickness") * CLHEP::cm;
90 G4double GlueThickness = activeParams.getLength(
"GlueThickness") * CLHEP::cm;
91 G4double KaptonThickness = activeParams.getLength(
"KaptonThickness") * CLHEP::cm;
92 G4double AluminiumThickness = activeParams.getLength(
"AluminiumThickness") * CLHEP::cm;
93 G4double foamlen = activeParams.getLength(
"Foam/length") * CLHEP::cm;
94 G4double foamwid = activeParams.getLength(
"Foam/width") * CLHEP::cm;
95 G4double foamthick = activeParams.getLength(
"Foam/thickness") * CLHEP::cm;
96 G4double foamZ = activeParams.getLength(
"Foam/zPosition") * CLHEP::cm;
97 G4double foamY = activeParams.getLength(
"Foam/yShift") * CLHEP::cm;
98 G4double flexlen = activeParams.getLength(
"Flex/length") * CLHEP::cm;
99 G4double flexwid = activeParams.getLength(
"Flex/width") * CLHEP::cm;
100 G4double flexZ = activeParams.getLength(
"Flex/zPosition") * CLHEP::cm;
101 G4double allen1 = activeParams.getLength(
"AlBat/length1") * CLHEP::cm;
102 G4double allen2 = activeParams.getLength(
"AlBat/length2") * CLHEP::cm;
103 G4double allen3 = activeParams.getLength(
"AlBat/length3") * CLHEP::cm;
104 G4double alwid = activeParams.getLength(
"AlBat/width") * CLHEP::cm;
105 G4double althick1 = activeParams.getLength(
"AlBat/thickness1") * CLHEP::cm;
106 G4double althick2 = activeParams.getLength(
"AlBat/thickness2") * CLHEP::cm;
107 G4double althick3 = activeParams.getLength(
"AlBat/thickness3") * CLHEP::cm;
129 G4double dz_ladder = foamthick / 2. + KaptonThickness + GlueThickness +
130 AluminiumThickness + 3. * AirGap;
131 G4double dz_sensor = (SubstrateThickness + EpitaxialThickness +
132 MetalThickness + 2. * AirGap) / 2.;
134 G4double fullLength = allen1 + allen2 + allen3 - (foamZ - foamlen / 2.);
135 double zshift = fullLength / 2. - (allen1 + allen2 + allen3);
138 G4AssemblyVolume* support =
new G4AssemblyVolume();
139 G4Transform3D transl;;
142 G4Box* s_foam =
new G4Box(
"s_foam", foamlen / 2., foamwid / 2., foamthick / 2.);
144 transl = G4Translate3D(foamZ + zshift, foamY, 0);
145 support->AddPlacedVolume(l_foam, transl);
148 G4Box* s_glue =
new G4Box(
"s_glue", flexlen / 2., foamwid / 2., GlueThickness / 2.);
150 G4double r_glue = foamthick / 2. + AirGap + GlueThickness / 2.;
151 transl = G4Translate3D(flexZ + zshift, foamY, -r_glue);
152 support->AddPlacedVolume(l_glue, transl);
153 transl = G4Translate3D(flexZ + zshift, foamY, r_glue);
154 support->AddPlacedVolume(l_glue, transl);
157 G4Box* s_kapton =
new G4Box(
"s_kapton", flexlen / 2., flexwid / 2., KaptonThickness / 2.);
159 G4double r_Kapton = r_glue + AirGap + (GlueThickness + KaptonThickness) / 2.;
160 transl = G4Translate3D(flexZ + zshift, 0, -r_Kapton);
161 support->AddPlacedVolume(l_kapton, transl);
162 transl = G4Translate3D(flexZ + zshift, 0, r_Kapton);
163 support->AddPlacedVolume(l_kapton, transl);
166 G4Box* s_metal =
new G4Box(
"s_metal", flexlen / 2., flexwid / 2., AluminiumThickness / 2.);
168 G4double r_metal = r_Kapton + AirGap + (KaptonThickness + AluminiumThickness) / 2.;
169 transl = G4Translate3D(flexZ + zshift, 0, -r_metal);
170 support->AddPlacedVolume(l_metal, transl);
171 transl = G4Translate3D(flexZ + zshift, 0, r_metal);
172 support->AddPlacedVolume(l_metal, transl);
175 G4Box* s_bat1p =
new G4Box(
"s_bat1p", allen1 / 2., alwid / 2., althick1 / 2.);
176 double incut = foamZ + foamlen / 2.;
177 G4Box* s_bat1t =
new G4Box(
"s_bat1t", incut, foamwid / 2., althick1 + 0.2 / 2.);
179 G4Transform3D tt = G4Translate3D(-allen1 / 2., foamY, 0);
180 G4SubtractionSolid* s_bat1 =
new G4SubtractionSolid(
"s_bat1", s_bat1p, s_bat1t, tt);
181 G4Box* s_bat2 =
new G4Box(
"s_bat2", allen2 / 2., alwid / 2., althick2 / 2.);
182 G4Box* s_bat3 =
new G4Box(
"s_bat3", allen3 / 2., alwid / 2., althick3 / 2.);
186 transl = G4Translate3D(allen1 / 2. + zshift, 0, 0);
187 support->AddPlacedVolume(l_bat1, transl);
188 transl = G4Translate3D(allen1 + allen2 / 2. + zshift, 0, 0);
189 support->AddPlacedVolume(l_bat2, transl);
190 transl = G4Translate3D(allen1 + allen2 + allen3 / 2. + zshift, 0, 0);
191 support->AddPlacedVolume(l_bat3, transl);
193 l_foam->SetVisAttributes(FoamVisAtt);
194 l_kapton->SetVisAttributes(KaptonVisAtt);
195 l_glue->SetVisAttributes(GlueVisAtt);
196 l_metal->SetVisAttributes(MetalVisAtt);
199 G4String symbol, name;
203 a = 1.01 * CLHEP::g / CLHEP::mole;
204 G4Element* elH =
new G4Element(name =
"Hydrogen", symbol =
"H", zz = 1., a);
205 a = 12.01 * CLHEP::g / CLHEP::mole;
206 G4Element* elC =
new G4Element(name =
"Carbon", symbol =
"C", zz = 6., a);
207 a = 16.00 * CLHEP::g / CLHEP::mole;
208 G4Element* elO =
new G4Element(name =
"Oxygen", symbol =
"O", zz = 8., a);
209 density = 1.31 * CLHEP::g / CLHEP::cm3;
210 G4Material* peekMat =
new G4Material(name =
"plumePeek", density, ncomponents = 3);
211 peekMat->AddElementByMassFraction(elC, 0.76);
212 peekMat->AddElementByMassFraction(elH, 0.08);
213 peekMat->AddElementByMassFraction(elO, 0.16);
214 density = 1.5 * CLHEP::g / CLHEP::cm3;
215 G4Material* carbMat =
new G4Material(name =
"plumeCarb", density, ncomponents = 1);
216 carbMat->AddElementByMassFraction(elC, 1.0);
218 G4Box* s_sensor =
new G4Box(
"s_sensor", SensorLengthX / 2., SensorLengthY / 2., dz_sensor);
222 G4Box* s_substrate =
new G4Box(
"s_substrate", SensorLengthX / 2., SensorLengthY / 2., SubstrateThickness / 2.);
223 G4LogicalVolume* l_substrate =
new G4LogicalVolume(s_substrate,
geometry::Materials::get(
"Silicon"),
"PLUME.l_substrate");
224 new G4PVPlacement(0, G4ThreeVector(0, 0, -dz_sensor + SubstrateThickness / 2.), l_substrate,
"p_substrate", l_sensor,
false, 1);
227 G4Box* s_epitaxial =
new G4Box(
"s_epitaxial", SensorLengthX / 2., SensorLengthY / 2., EpitaxialThickness / 2.);
228 G4LogicalVolume* l_epitaxial =
new G4LogicalVolume(s_epitaxial,
geometry::Materials::get(
"Silicon"),
"PLUME.l_epitaxial", 0,
230 new G4PVPlacement(0, G4ThreeVector(0, 0, -dz_sensor + AirGap + SubstrateThickness + EpitaxialThickness / 2.), l_epitaxial,
231 "p_epitaxial", l_sensor,
false, 1);
233 G4Box* s_metalized =
new G4Box(
"s_metalized", SensorLengthX / 2., SensorLengthY / 2., MetalThickness / 2.);
234 G4LogicalVolume* l_metalized =
new G4LogicalVolume(s_metalized,
geometry::Materials::get(
"SiO2Al"),
"PLUME.l_metalized");
235 new G4PVPlacement(0, G4ThreeVector(0, 0, -dz_sensor + 2 * AirGap + SubstrateThickness + EpitaxialThickness + MetalThickness / 2.),
236 l_metalized,
"p_metalized", l_sensor,
false, 1);
248 l_epitaxial->SetUserLimits(
new G4UserLimits(stepSize));
249 l_substrate->SetVisAttributes(SubstrateVisAtt);
250 l_epitaxial->SetVisAttributes(EpitaxialVisAtt);
251 l_metalized->SetVisAttributes(MetalizedVisAtt);
254 double zSens = dz_ladder + dz_sensor + AirGap;
256 G4ThreeVector Ta(0, 0, 0);
257 G4Transform3D tra(ra, Ta);
260 GearDir pipeParams(content,
"CoolingPipes");
261 double pipeLen = pipeParams.getLength(
"length") * CLHEP::cm;
262 double pipeInR = pipeParams.getLength(
"rIn") * CLHEP::cm;
263 double pipeOutR = pipeParams.getLength(
"rOut") * CLHEP::cm;
264 G4Tubs* pipe_s =
new G4Tubs(
"pipe_s", pipeInR, pipeOutR, pipeLen / 2., 0, 2 * M_PI);
266 G4LogicalVolume* pipe_l =
new G4LogicalVolume(pipe_s, carbMat,
"PLUME.l_pipe");
267 double pipeX = pipeParams.getLength(
"x") * CLHEP::cm;
268 double pipeY = pipeParams.getLength(
"y") * CLHEP::cm;
269 double pipeZ = pipeParams.getLength(
"z") * CLHEP::cm;
270 G4double alphaPipe = pipeParams.getAngle(
"alpha");
273 G4Box* s_sensors =
new G4Box(
"s_sensors", flexlen / 2., flexwid / 2., dz_sensor + AirGap);
275 for (
auto ladder : content.getNodes(
"Placements/Ladder")) {
277 std::string
id = ladder.getString(
"@id");
279 G4AssemblyVolume* assemblyLadder =
new G4AssemblyVolume();
280 assemblyLadder->AddPlacedAssembly(support, tra);
282 for (
auto pars : content.getNodes(
"SensorAlignment/ladder")) {
283 if (pars.getString(
"@id") ==
id) {
284 GearDir alignPars(pars);
286 for (
auto sidePars : alignPars.getNodes(
"side")) {
289 G4Transform3D transformSens;
290 if (sidePars.getString(
"@id") ==
"mirror") mirror = 1;
294 for (
auto sensorPars : sidePars.getNodes(
"sensor")) {
295 double x = sensorPars.getLength(
"x") * CLHEP::cm + zshift;
296 double y = sensorPars.getLength(
"y") * CLHEP::cm;
297 int ids = sensorPars.getInt(
"id");
300 double alpha1 = sensorPars.getAngle(
"alpha") * CLHEP::rad;
301 if (mirror) transformSens = G4Translate3D(x, y,
302 0) * G4RotateZ3D(alpha1) * G4RotateX3D(M_PI);
303 else transformSens = G4Translate3D(x, y, 0) * G4RotateZ3D(-alpha1);
305 new G4PVPlacement(transformSens, l_sensor,
"p_sensor", l_sensors,
true,
306 (ids + mirror * 6 + LadderID * 12));
308 transformSens = G4Translate3D(0, 0, mirror ? -zSens : zSens);
309 assemblyLadder->AddPlacedVolume(l_sensors, transformSens);
316 G4double thetaZ = ladder.getAngle(
"ThetaZ");
317 G4double r = ladder.getLength(
"r_plume") * CLHEP::cm - zshift * sin(thetaZ);
318 G4double z = ladder.getLength(
"z_plume") * CLHEP::cm - zshift * cos(thetaZ);
319 G4double phi = ladder.getAngle(
"Phi");
320 G4double alpha = ladder.getAngle(
"Alpha");
324 G4Transform3D transformPipe = G4Translate3D(pipeZ + zshift, pipeY, pipeX) * G4RotateZ3D(alphaPipe) * G4RotateY3D(+ M_PI / 2.);
325 assemblyLadder->AddPlacedVolume(pipe_l, transformPipe);
326 transformPipe = G4Translate3D(pipeZ + zshift
327 , pipeY, -pipeX) * G4RotateZ3D(alphaPipe) * G4RotateY3D(+ M_PI / 2.);
328 assemblyLadder->AddPlacedVolume(pipe_l, transformPipe);
331 G4Transform3D transform1;
335 transform1 = G4Translate3D(-fullLength / 2. - 18.7, 0., 1.2) * G4RotateY3D(thetaZ) * G4RotateX3D(M_PI / 2.) * G4RotateY3D(M_PI);
336 assemblyLadder->AddPlacedAssembly(sup1, transform1);
337 transform1 = G4Translate3D(fullLength / 2. + 19.5, 0, -2.6) * G4RotateY3D(thetaZ) * G4RotateX3D(M_PI / 2.);
338 assemblyLadder->AddPlacedAssembly(sup2, transform1);
344 transform1 = G4Translate3D(fullLength / 2. + 5, 0,
345 7.8) * G4RotateY3D(-thetaZ) * G4RotateZ3D(M_PI / 2.) * G4RotateX3D(M_PI / 2.) * G4RotateY3D(M_PI / 2.);
346 assemblyLadder->AddPlacedAssembly(sup1, transform1);
347 transform1 = G4Translate3D(-fullLength / 2. - 31, 0.,
348 -21.7) * G4RotateY3D(-thetaZ) * G4RotateZ3D(M_PI / 2.) * G4RotateX3D(-M_PI / 2.) * G4RotateY3D(M_PI);
349 assemblyLadder->AddPlacedAssembly(sup2, transform1);
352 G4Transform3D transform = G4RotateZ3D(phi) * G4Translate3D(r, 0,
353 z) * G4RotateY3D(thetaZ) * G4RotateZ3D(alpha + M_PI) * G4RotateY3D(- M_PI / 2.0);
355 assemblyLadder->MakeImprint(&topVolume, transform);
static G4Material * get(const std::string &name)
Find given material.
G4AssemblyVolume * buildSupport1()
Volume of support1.
G4AssemblyVolume * buildSupport2()
Volume of support2.
G4AssemblyVolume * buildSupport4()
Volume of support4.
G4AssemblyVolume * buildSupport3()
Volume of support3.
SensitiveDetector * m_sensitive
SensitiveDetector PLUME.
VXD::SensitiveDetector< PXDSimHit, PXDTrueHit > SensitiveDetector
The PXD Sensitive Detector class.