ATLAS Offline Software
GeoPixelLadderServices.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 //
6 // This class builds the Ladder services (Omega + Al tube + glue + part of pigtail + connector)
7 //
8 //
9 #include "GeoPixelLadderServices.h"
10 #include "GeoPixelFluid.h"
11 #include "GeoPixelCable.h"
12 #include "GeoModelKernel/GeoTubs.h"
13 #include "GeoModelKernel/GeoBox.h"
14 #include "GeoModelKernel/GeoPara.h"
15 #include "GeoModelKernel/GeoLogVol.h"
16 #include "GeoModelKernel/GeoPhysVol.h"
17 #include "GeoModelKernel/GeoMaterial.h"
18 #include "GeoModelKernel/GeoTransform.h"
19 #include "GeoModelKernel/GeoShapeShift.h"
20 #include "GeoModelKernel/GeoShapeUnion.h"
21 #include "GaudiKernel/PhysicalConstants.h"
22 
23 #include <cmath>
24 #include <utility>
25 #include <vector>
26 
27 GeoPixelLadderServices::GeoPixelLadderServices(InDetDD::PixelDetectorManager* m_DDmgr,
28  PixelGeometryManager* mgr,
29  GeoModelIO::ReadGeoModel* sqliteReader,
30  std::shared_ptr<std::map<std::string, GeoFullPhysVol*>> mapFPV,
31  std::shared_ptr<std::map<std::string, GeoAlignableTransform*>> mapAX,
32  int ladderType)
33  : GeoVPixelFactory(m_DDmgr, mgr, sqliteReader, std::move(mapFPV), std::move(mapAX)),
34  m_ladderType(ladderType)
35 {
36  //std::cout << "Building GeoPixelLadderServices with ladder type : " << ladderType << std::endl;
37  //
38  // Define the log volume in the constructor, so I do it only once.
39  //
40  // Length and thickness of the Omega are in the db
41  //
42  double halflength =0.5* m_gmt_mgr->PixelLadderLength()+m_epsilon;
43  //
44  const GeoMaterial* air = m_mat_mgr->getMaterial("std::Air");
45  //
46  //const GeoBox* ladderSvcBox = new GeoBox(thickness/2.,width/2.,halflength);
47  // Quick fix - we hardwire the numbers. Need to work out a way to extract this from the database numbers.
48  double safety = 0.01*Gaudi::Units::mm;
49  double xBase = 0;
50 
51  // ConnA: Part to fit Connector
52  // ConnC: Part to fit Cable
53 
54  double xOffsetConnA = xBase + m_gmt_mgr->PixelLadderCableOffsetX() - m_gmt_mgr->PixelLadderServicesX();
55  double xOffsetConnC = xOffsetConnA;
56 
57 
58  //double xMaxConnA = 1.250 * Gaudi::Units::cm + xOffsetConnA + safety;
59  // double xMaxConnC = 1.6575 * Gaudi::Units::cm + xOffsetConnC + safety;
60  // max offset is 12.5mm + 1/2 thickness of cables
61  // double xMaxConnA = 1.5075 * Gaudi::Units::cm + 0.5* 0.15*Gaudi::Units::cm + xOffsetConnA + safety;
62  double xMaxConnA = m_gmt_mgr->PixelConnectorPosX(1) + 0.5*m_gmt_mgr->PixelConnectorWidthX(1) + xOffsetConnA + safety;
63  double xMaxConnC = 1.25 * Gaudi::Units::cm + 0.5* 0.0125*Gaudi::Units::cm + xOffsetConnC + safety;
64  double xMaxOmegaBase = 0.055 * Gaudi::Units::cm + xBase + 1*Gaudi::Units::mm; // The 1 mm is just extra safety.
65  double yWidthConnA = 1.0 * Gaudi::Units::cm;
66  double yWidthConnC = 0.2 * Gaudi::Units::cm;
67  double yPosConnA = m_gmt_mgr->PixelLadderCableOffsetY() - m_gmt_mgr->PixelLadderServicesY();
68  double yPosConnC = yPosConnA;
69  double xCenter = 0;
70  double xWidthOmegaBase = xMaxOmegaBase - xBase;
71  double xWidthConnA = xMaxConnA - xBase;
72  double xWidthConnC = xMaxConnC - xBase;
73  double yWidthOmega = 1.2*Gaudi::Units::cm + m_epsilon;
74 
75  const GeoBox* omegaBaseEnv = new GeoBox(0.5*xWidthOmegaBase, 0.5*yWidthOmega, halflength);
76  const GeoBox* connAEnv = new GeoBox(0.5*xWidthConnA, 0.5*yWidthConnA + safety, halflength);
77  const GeoBox* connCEnv = new GeoBox(0.5*xWidthConnC, 0.5*yWidthConnC + safety, halflength);
78  GeoTrf::Transform3D omegaBaseShift = GeoTrf::Translate3D(xBase + 0.5*xWidthOmegaBase - xCenter, 0., 0.);
79  GeoTrf::Transform3D connAShift = GeoTrf::Translate3D(xBase + 0.5*xWidthConnA - xCenter, yPosConnA, 0.);
80  GeoTrf::Transform3D connCShift = GeoTrf::Translate3D(xBase + 0.5*xWidthConnC - xCenter, yPosConnC, 0.);
81 
82  const GeoShape & ladderSvcShape = ((*omegaBaseEnv) << omegaBaseShift)
83  .add(*connAEnv << connAShift)
84  .add(*connCEnv << connCShift);
85 
86  m_ladderServicesLV = new GeoLogVol("LadderSvc",&ladderSvcShape,air);
87  // store coordinates of corners of envelope.
88  m_envelopeCornerA1 = GeoTrf::Vector3D(xMaxConnA - xCenter, yPosConnA + 0.5* yWidthConnA + safety, 0.);
89  m_envelopeCornerA2 = GeoTrf::Vector3D(xMaxConnA - xCenter, yPosConnA - 0.5* yWidthConnA + safety, 0.);
90  m_envelopeCornerC1 = GeoTrf::Vector3D(xMaxConnC - xCenter, yPosConnC + 0.5* yWidthConnC + safety, 0.);
91  m_envelopeCornerC2 = GeoTrf::Vector3D(xMaxConnC - xCenter, yPosConnC - 0.5* yWidthConnC + safety, 0.);
92 
93 
94  m_xOffset = xBase;
95  m_yOffset = 0;
96 
97 }
98 
99 
100 
101 
102 GeoVPhysVol* GeoPixelLadderServices::Build() {
103  GeoPhysVol* ladderSvcPhys = new GeoPhysVol(m_ladderServicesLV);
104  //double thickness = m_gmt_mgr->PixelLadderThickness()+m_gmt_mgr->PixelCableThickness();
105  //double thickness = m_gmt_mgr->PixelLadderThickness() + m_gmt_mgr->PixelCableThickness() + 0.25*Gaudi::Units::cm; // plus 0.25 cm New DC3 ???
106  //double thickness = m_gmt_mgr->PixelLadderThickness()+ 6.5; // m_gmt_mgr->PixelCableThickness() was 0.4 cm, plus 0.25 cm New DC3
107  //
108  // The Glue
109  BuildGlue(ladderSvcPhys);
110 
111  // The aluminium tube
112  // Positioned at 0,0,0 relative to parent.
113  ladderSvcPhys->add(this->BuildAlTube());
114 
115  // The carbon-carbon Omega
116  // Positioned at 0,0,0 relative to parent.
117  ladderSvcPhys->add(this->BuildOmega());
118 
119  //
120  // Put the fluid/gaz mixture
121  //
122  GeoPixelFluid fluid(m_DDmgr, m_gmt_mgr, m_sqliteReader, m_mapFPV, m_mapAX, m_ladderType);
123  double xpos = fluid.posX() + m_xOffset;
124  double ypos = fluid.posY() + m_yOffset;
125  double zpos = fluid.posZ();
126 
127  GeoVPhysVol* fluidPhys = fluid.Build();
128  GeoTransform* xform = new GeoTransform(GeoTrf::Translate3D(xpos,ypos,zpos));
129  ladderSvcPhys->add(xform);
130  ladderSvcPhys->add(fluidPhys);
131 
132  //
133  // Place Pigtail + connector
134  //
135  BuildPigtailAndConnector(ladderSvcPhys);
136 
137 
138  //
139  // Cables
140  //
141  GeoPixelCable cable (m_DDmgr, m_gmt_mgr, m_sqliteReader, m_mapFPV, m_mapAX);
142  for (int ii = 0; ii < cable.numElements(); ii++) {
143  cable.setElement(ii);
144  GeoVPhysVol* phys = cable.Build();
145  double xpos = m_xOffset + cable.getStackOffset() + m_gmt_mgr->PixelLadderCableOffsetX() - m_gmt_mgr->PixelLadderServicesX();
146  double ypos = m_yOffset + m_gmt_mgr->PixelLadderCableOffsetY() - m_gmt_mgr->PixelLadderServicesY();
147  double zpos = cable.zpos();
148  GeoTransform* xform = new GeoTransform(GeoTrf::Translate3D(xpos,ypos,zpos));
149  ladderSvcPhys->add(xform);
150  ladderSvcPhys->add(phys);
151  }
152 
153  return ladderSvcPhys;
154 }
155 
156 GeoVPhysVol* GeoPixelLadderServices::BuildOmega() {
157 
158  double xOffset = m_xOffset;
159  double yOffset = m_yOffset;
160  /*
161  double xUpperBend = xOffset + 2.9*Gaudi::Units::mm;
162  double yUpperBend = yOffset + 0;
163  double radUpperBend = 2.3*Gaudi::Units::mm;
164  double xLowerBend = xOffset + 0.9;
165  double yLowerBend = yOffset + 3.35*Gaudi::Units::mm;
166  double radLowerBend = 0.8*Gaudi::Units::mm;
167  double yStart= yOffset + (4.675+0.5*2.65)*Gaudi::Units::mm;
168  double yEnd = yOffset -yStart;
169  double thick = 0.3*Gaudi::Units::mm;
170  double length = 816*Gaudi::Units::mm;
171  double zOffset = 0;
172  */
173  double xUpperBend = xOffset + m_gmt_mgr->PixelOmegaUpperBendX();
174  double yUpperBend = yOffset + m_gmt_mgr->PixelOmegaUpperBendY();
175  double radUpperBend = m_gmt_mgr->PixelOmegaUpperBendRadius();
176  double xLowerBend = xOffset + m_gmt_mgr->PixelOmegaLowerBendX();
177  double yLowerBend = yOffset + m_gmt_mgr->PixelOmegaLowerBendY();
178  double radLowerBend = m_gmt_mgr->PixelOmegaLowerBendRadius();
179  double yStart= yOffset + m_gmt_mgr->PixelOmegaStartY();
180  double yEnd = yOffset + m_gmt_mgr->PixelOmegaEndY();
181  double thick = m_gmt_mgr->PixelOmegaWallThickness();
182  double length = m_gmt_mgr->PixelOmegaLength();
183  double zOffset = 0;
184 
185  double sepX = (xUpperBend - xLowerBend);
186  double sepY = (yUpperBend - yLowerBend);
187  double sep2 = sepX*sepX+sepY*sepY;
188  double sep = sqrt(sep2);
189  double radDist = radUpperBend+radLowerBend-thick;
190  double alpha = atan(-sepY/sepX) - asin(radDist/sep);
191  double upperStraightLength = sqrt(sep2 - radDist*radDist);
192  double xLowerStraight = xLowerBend - (radLowerBend-0.5*thick)*sin(alpha);
193  double yLowerStraight = yLowerBend - (radLowerBend-0.5*thick)*cos(alpha);
194  double xUpperStraight = xUpperBend + (radUpperBend-0.5*thick)*sin(alpha);
195  double yUpperStraight = yUpperBend + (radUpperBend-0.5*thick)*cos(alpha);
196 
197 
198  // Tube sector for upper bend
199  GeoTubs * upperBendShape = new GeoTubs(radUpperBend - thick, radUpperBend, 0.5* length, alpha-0.5*Gaudi::Units::pi, Gaudi::Units::pi - 2*alpha);
200 
201  // Tube sector for lower bend (+y)
202  GeoTubs * lowerBendShapeP = new GeoTubs(radLowerBend - thick, radLowerBend, 0.5* length, Gaudi::Units::pi, 0.5*Gaudi::Units::pi-alpha);
203 
204  // Tube sector for lower bend (-y)
205  GeoTubs * lowerBendShapeM = new GeoTubs(radLowerBend - thick, radLowerBend, 0.5* length, 0.5*Gaudi::Units::pi + alpha, 0.5*Gaudi::Units::pi-alpha);
206 
207  // Lower Straight section (+y)
208  GeoBox * lowerStraightBoxP = new GeoBox(0.5*thick, 0.5*(yStart - yLowerBend), 0.5*length);
209 
210  // Lower Straight section (-y)
211  GeoBox * lowerStraightBoxM = new GeoBox(0.5*thick, 0.5*(-yLowerBend - yEnd), 0.5*length);
212 
213  // Upper Straight section (+-y)
214  GeoBox * upperStraightBox = new GeoBox(0.5*thick, 0.5*upperStraightLength, 0.5*length);
215 
216  // Put it all together
217  const GeoShape & omegaShape =
218  (*lowerStraightBoxP << GeoTrf::Translate3D(xLowerBend-radLowerBend+0.5*thick,0.5*(yLowerBend+yStart),zOffset) )
219  .add(*lowerBendShapeP << GeoTrf::Translate3D(xLowerBend,yLowerBend,zOffset) )
220  .add(*upperStraightBox << GeoTrf::Translate3D(0.5*(xLowerStraight+xUpperStraight),0.5*(yLowerStraight+yUpperStraight),zOffset)*GeoTrf::RotateZ3D(0.5*Gaudi::Units::pi-alpha) )
221  .add(*upperBendShape << GeoTrf::Translate3D(xUpperBend,yUpperBend,zOffset) )
222  .add(*upperStraightBox << GeoTrf::Translate3D(0.5*(xLowerStraight+xUpperStraight),-0.5*(yLowerStraight+yUpperStraight),zOffset)*GeoTrf::RotateZ3D(0.5*Gaudi::Units::pi+alpha) )
223  .add(*lowerBendShapeM << GeoTrf::Translate3D(xLowerBend,-yLowerBend,zOffset) )
224  .add(*lowerStraightBoxM << GeoTrf::Translate3D(xLowerBend-radLowerBend+0.5*thick,0.5*(-yLowerBend+yEnd),zOffset) );
225 
226  double totVolume =
227  lowerStraightBoxP->volume()
228  + lowerBendShapeP->volume()
229  + 2*upperStraightBox->volume()
230  + upperBendShape->volume()
231  + lowerBendShapeM->volume()
232  + lowerStraightBoxM->volume();
233 
234  std::string matName = m_gmt_mgr->getMaterialName("Omega", m_gmt_mgr->GetLD());
235  const GeoMaterial* omegaMat = m_mat_mgr->getMaterialForVolume(matName,totVolume);
236  GeoLogVol* omegaLV = new GeoLogVol("Omega",&omegaShape,omegaMat);
237  return new GeoPhysVol(omegaLV);
238 
239 }
240 
241 GeoVPhysVol* GeoPixelLadderServices::BuildAlTube() {
242 
243  double xOffset = m_xOffset;
244  double yOffset = m_yOffset;
245 
246  /*
247  double xUpperBend = xOffset + 2.7*Gaudi::Units::mm;
248  double yUpperBend = yOffset;
249  double radUpperBend = 2.0*Gaudi::Units::mm;
250  double xLowerBend = xOffset + 0.55;
251  double yLowerBend = yOffset+1.925*Gaudi::Units::mm;
252  double radLowerBend = 0.5*Gaudi::Units::mm;
253  double thick = 0.2*Gaudi::Units::mm;
254  double length = 838*Gaudi::Units::mm;
255  double zOffset = 0;
256  */
257 
258  double xUpperBend = xOffset + m_gmt_mgr->PixelAlTubeUpperBendX();
259  double yUpperBend = yOffset + m_gmt_mgr->PixelAlTubeUpperBendY();
260  double radUpperBend = m_gmt_mgr->PixelAlTubeUpperBendRadius();
261  double xLowerBend = xOffset + m_gmt_mgr->PixelAlTubeLowerBendX();
262  double yLowerBend = yOffset + m_gmt_mgr->PixelAlTubeLowerBendY();
263  double radLowerBend = m_gmt_mgr->PixelAlTubeLowerBendRadius();
264  double thick = m_gmt_mgr->PixelAlTubeWallThickness();
265  double length = m_gmt_mgr->PixelAlTubeLength();
266  double zOffset = 0;
267 
268  double sepX = (xUpperBend - xLowerBend);
269  double sepY = (yUpperBend - yLowerBend);
270  double sep2 = sepX*sepX+sepY*sepY;
271  double sep = sqrt(sep2);
272  double radDist = radUpperBend-radLowerBend;
273  double alpha = atan(-sepY/sepX) - asin(radDist/sep);
274  double upperStraightLength = sqrt(sep2 - radDist*radDist);
275  double xLowerStraight = xLowerBend + (radLowerBend-0.5*thick)*sin(alpha);
276  double yLowerStraight = yLowerBend + (radLowerBend-0.5*thick)*cos(alpha);
277  double xUpperStraight = xUpperBend + (radUpperBend-0.5*thick)*sin(alpha);
278  double yUpperStraight = yUpperBend + (radUpperBend-0.5*thick)*cos(alpha);
279 
280 
281  // Tube sector for upper bend
282  GeoTubs * upperBendShape = new GeoTubs(radUpperBend - thick, radUpperBend, 0.5* length, alpha-0.5*Gaudi::Units::pi, Gaudi::Units::pi - 2*alpha);
283 
284  // Tube sector for lower bend (+y)
285  GeoTubs * lowerBendShapeP = new GeoTubs(radLowerBend - thick, radLowerBend, 0.5* length, 0.5*Gaudi::Units::pi-alpha, 0.5*Gaudi::Units::pi+alpha);
286 
287  // Tube sector for lower bend (-y)
288  GeoTubs * lowerBendShapeM = new GeoTubs(radLowerBend - thick, radLowerBend, 0.5* length, Gaudi::Units::pi, 0.5*Gaudi::Units::pi+alpha);
289 
290  // Lower Straight section
291  GeoBox * lowerStraightBox = new GeoBox(0.5*thick, yLowerBend, 0.5*length);
292 
293  // Upper Straight section (+-y)
294  GeoBox * upperStraightBox = new GeoBox(0.5*thick, 0.5*upperStraightLength, 0.5*length);
295 
296  // Put it all together
297  const GeoShape & alTubeShape =
298  (*lowerStraightBox << GeoTrf::Translate3D(xLowerBend-radLowerBend+0.5*thick,0,zOffset) )
299  .add(*lowerBendShapeP << GeoTrf::Translate3D(xLowerBend,yLowerBend,zOffset) )
300  .add(*upperStraightBox << GeoTrf::Translate3D(0.5*(xLowerStraight+xUpperStraight),0.5*(yLowerStraight+yUpperStraight),zOffset)*GeoTrf::RotateZ3D(0.5*Gaudi::Units::pi-alpha) )
301  .add(*upperBendShape << GeoTrf::Translate3D(xUpperBend,yUpperBend,zOffset) )
302  .add(*upperStraightBox << GeoTrf::Translate3D(0.5*(xLowerStraight+xUpperStraight),-0.5*(yLowerStraight+yUpperStraight),zOffset)*GeoTrf::RotateZ3D(0.5*Gaudi::Units::pi+alpha) )
303  .add(*lowerBendShapeM << GeoTrf::Translate3D(xLowerBend,-yLowerBend,zOffset) );
304 
305  double totVolume =
306  lowerStraightBox->volume()
307  + lowerBendShapeP->volume()
308  + 2*upperStraightBox->volume()
309  + upperBendShape->volume()
310  + lowerBendShapeM->volume();
311 
312 
313  std::string matName = m_gmt_mgr->getMaterialName("AlTube", m_gmt_mgr->GetLD());
314  const GeoMaterial* alMat = m_mat_mgr->getMaterialForVolume(matName,totVolume);
315  GeoLogVol* tubeLV = new GeoLogVol("AlTube",&alTubeShape,alMat);
316  return new GeoPhysVol(tubeLV);
317 
318 
319 }
320 
321 void GeoPixelLadderServices::BuildGlue(GeoPhysVol * parent) {
322 
323 
324  int nGlueElements = m_gmt_mgr->PixelNumOmegaGlueElements();
325  for (int i = 0; i < nGlueElements; ++i) {
326  double x1 = m_gmt_mgr->PixelOmegaGlueStartX(i);
327  double thickness = m_gmt_mgr->PixelOmegaGlueThickness(i);
328  double y1 = m_gmt_mgr->PixelOmegaGlueStartY(i);
329  double y2 = m_gmt_mgr->PixelOmegaGlueEndY(i);
330  double length = m_gmt_mgr->PixelOmegaGlueLength(i);
331  double zOffset = m_gmt_mgr->PixelOmegaGluePosZ(i);
332  int typeNum = m_gmt_mgr->PixelOmegaGlueTypeNum(i);
333 
334  GeoBox * glueShape = new GeoBox(0.5*thickness, 0.5*std::abs(y1-y2), 0.5*length);
335  std::string matName = m_gmt_mgr->getMaterialName("GlueOmegaStave", m_gmt_mgr->GetLD(), typeNum);
336  const GeoMaterial* glueMat = m_mat_mgr->getMaterialForVolume(matName, glueShape->volume());
337  GeoLogVol* glueLV = new GeoLogVol("Glue",glueShape,glueMat);
338  GeoPhysVol* gluePV = new GeoPhysVol(glueLV);
339 
340  GeoTransform * xform = new GeoTransform(GeoTrf::Translate3D(m_xOffset+x1+0.5*thickness, m_yOffset-0.5*(y1+y2), zOffset));
341 
342  parent->add(xform);
343  parent->add(gluePV);
344  }
345 }
346 
347 
348 void GeoPixelLadderServices::BuildPigtailAndConnector(GeoPhysVol * parent) {
349 
350  double xOffset = m_xOffset + m_gmt_mgr->PixelLadderCableOffsetX() - m_gmt_mgr->PixelLadderServicesX();
351  double yOffset = m_yOffset + m_gmt_mgr->PixelLadderCableOffsetY() - m_gmt_mgr->PixelLadderServicesY();
352 
353  // Pigtail flat section
354  double xPosPigtail = xOffset + m_gmt_mgr->PixelPigtailPosX();
355  double yPosPigtail = yOffset + 0.5*(m_gmt_mgr->PixelPigtailStartY() + m_gmt_mgr->PixelPigtailEndY());
356  double zPosPigtail = m_gmt_mgr->PixelPigtailPosZ();
357  double xWidthPigtail = m_gmt_mgr->PixelPigtailThickness();
358  double yWidthPigtail = std::abs(m_gmt_mgr->PixelPigtailStartY() - m_gmt_mgr->PixelPigtailEndY());
359  double zWidthPigtail = m_gmt_mgr->PixelPigtailFlatWidthZ();
360  std::string matNamePigtail = m_gmt_mgr->getMaterialName("PigtailFlat", m_gmt_mgr->GetLD());
361  GeoBox * pigtailShape = new GeoBox(0.5*xWidthPigtail, 0.5*yWidthPigtail, 0.5*zWidthPigtail);
362  const GeoMaterial* pigtailMat = m_mat_mgr->getMaterialForVolume(matNamePigtail,pigtailShape->volume());
363  GeoLogVol * pigtailLV = new GeoLogVol("PigtailFlat",pigtailShape,pigtailMat);
364  GeoPhysVol * pigtailPhys = new GeoPhysVol(pigtailLV);
365 
366  // Connector
367  int numConn = m_gmt_mgr->PixelNumConnectorElements();
368  std::vector<GeoPhysVol *> connectorPhysVols(numConn);
369  std::vector<GeoTrf::Transform3D> connectorTransforms(numConn);
370  for (int iConn = 0; iConn < numConn; iConn++) {
371  double xPosConnector = xOffset + m_gmt_mgr->PixelConnectorPosX(iConn);
372  double yPosConnector = yOffset + m_gmt_mgr->PixelConnectorPosY(iConn);
373  double zPosConnector = m_gmt_mgr->PixelConnectorPosZ(iConn);
374  double xWidthConnector = m_gmt_mgr->PixelConnectorWidthX(iConn);
375  double yWidthConnector = m_gmt_mgr->PixelConnectorWidthY(iConn);
376  double zWidthConnector = m_gmt_mgr->PixelConnectorWidthZ(iConn);
377  std::string matNameConnector = m_gmt_mgr->getMaterialName("Connector", m_gmt_mgr->GetLD(), iConn);
378  GeoBox * connectorShape = new GeoBox(0.5*xWidthConnector, 0.5*yWidthConnector, 0.5*zWidthConnector);
379  const GeoMaterial* connectorMat = m_mat_mgr->getMaterialForVolume(matNameConnector,connectorShape->volume());
380  GeoLogVol * connectorLV = new GeoLogVol("Connector",connectorShape,connectorMat);
381  connectorPhysVols[iConn] = new GeoPhysVol(connectorLV);
382  connectorTransforms[iConn] = GeoTrf::Translate3D(xPosConnector, yPosConnector, zPosConnector);
383  }
384 
385  for (int iModule = 0; iModule<m_gmt_mgr->PixelNModule(); iModule++) {
386  int moduleEta = m_gmt_mgr->PixelModuleEtaFromIndex(iModule);
387  double zShift = m_gmt_mgr->PixelModuleZPosition(moduleEta);
388 
389  // Place pigtail
390  GeoTransform * xformPigtail = new GeoTransform(GeoTrf::Translate3D(xPosPigtail, yPosPigtail, zShift+zPosPigtail));
391  parent->add(xformPigtail);
392  parent->add(pigtailPhys);
393 
394  // Place connector
395  for (int iConn = 0; iConn < numConn; iConn++) {
396  GeoTransform * xformConnector = new GeoTransform(GeoTrf::TranslateZ3D(zShift)*connectorTransforms[iConn]);
397  parent->add(xformConnector);
398  parent->add(connectorPhysVols[iConn]);
399  }
400  }
401 }
PixelGeometryManager::PixelPigtailPosX
virtual double PixelPigtailPosX()=0
PixelGeometryManager::PixelLadderCableOffsetY
virtual double PixelLadderCableOffsetY()=0
GeoPixelCable::zpos
double zpos()
Definition: GeoPixelCable.cxx:91
plotBeamSpotCompare.x1
x1
Definition: plotBeamSpotCompare.py:216
GeoPixelLadderServices::m_envelopeCornerA1
GeoTrf::Vector3D m_envelopeCornerA1
Definition: GeoPixelLadderServices.h:42
PixelGeometryManager::PixelOmegaStartY
virtual double PixelOmegaStartY()=0
PixelGeometryManager::PixelOmegaGlueThickness
virtual double PixelOmegaGlueThickness(int index)=0
GeoPixelFluid
Definition: GeoPixelFluid.h:13
PixelGeometryManager::PixelNumConnectorElements
virtual int PixelNumConnectorElements()=0
GeoPixelCable::Build
virtual GeoVPhysVol * Build() override
Definition: GeoPixelCable.cxx:50
PixelGeometryManager::PixelNumOmegaGlueElements
virtual int PixelNumOmegaGlueElements()=0
GeoPixelCable
Definition: GeoPixelCable.h:11
PixelGeometryManager
Definition: PixelGeometryManager.h:28
GeoPixelLadderServices::m_ladderType
int m_ladderType
Definition: GeoPixelLadderServices.h:40
PixelGeometryManager::PixelConnectorPosZ
virtual double PixelConnectorPosZ(int index)=0
GeoPixelLadderServices.h
PixelGeometryManager::PixelPigtailThickness
virtual double PixelPigtailThickness()=0
GeoPixelLadderServices::GeoPixelLadderServices
GeoPixelLadderServices(InDetDD::PixelDetectorManager *ddmgr, PixelGeometryManager *mgr, GeoModelIO::ReadGeoModel *sqliteReader, std::shared_ptr< std::map< std::string, GeoFullPhysVol * >> mapFPV, std::shared_ptr< std::map< std::string, GeoAlignableTransform * >> mapAX, int ladderType)
Definition: GeoPixelLadderServices.cxx:27
PixelGeometryManager::PixelLadderServicesY
virtual double PixelLadderServicesY()=0
PixelGeometryManager::PixelOmegaGlueEndY
virtual double PixelOmegaGlueEndY(int index)=0
PixelGeometryManager::PixelLadderServicesX
virtual double PixelLadderServicesX()=0
PixelGeometryManager::PixelOmegaLength
virtual double PixelOmegaLength()=0
PixelGeometryManager::PixelOmegaGlueLength
virtual double PixelOmegaGlueLength(int index)=0
GeoPixelLadderServices::BuildPigtailAndConnector
void BuildPigtailAndConnector(GeoPhysVol *parent)
Definition: GeoPixelLadderServices.cxx:348
PixelGeometryManager::PixelOmegaGlueTypeNum
virtual int PixelOmegaGlueTypeNum(int index)=0
PixelGeometryManager::PixelPigtailStartY
virtual double PixelPigtailStartY()=0
GeoPixelFluid::Build
virtual GeoVPhysVol * Build() override
Definition: GeoPixelFluid.cxx:35
GeoPixelLadderServices::m_envelopeCornerC1
GeoTrf::Vector3D m_envelopeCornerC1
Definition: GeoPixelLadderServices.h:44
GeoPixelLadderServices::BuildAlTube
GeoVPhysVol * BuildAlTube()
Definition: GeoPixelLadderServices.cxx:241
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
GeoPixelCable::setElement
void setElement(int i)
Definition: GeoPixelCable.h:20
GeoVPixelFactory::m_mapAX
std::shared_ptr< std::map< std::string, GeoAlignableTransform * > > m_mapAX
Definition: GeoVPixelFactory.h:48
GeoVPixelFactory::m_sqliteReader
GeoModelIO::ReadGeoModel * m_sqliteReader
Definition: GeoVPixelFactory.h:46
PixelGeometryManager::PixelPigtailFlatWidthZ
virtual double PixelPigtailFlatWidthZ()=0
GeoVPixelFactory::m_DDmgr
InDetDD::PixelDetectorManager * m_DDmgr
Definition: GeoVPixelFactory.h:45
makeTRTBarrelCans.y1
tuple y1
Definition: makeTRTBarrelCans.py:15
pi
#define pi
Definition: TileMuonFitter.cxx:65
drawFromPickle.atan
atan
Definition: drawFromPickle.py:36
GeoVPixelFactory::m_gmt_mgr
PixelGeometryManager * m_gmt_mgr
Definition: GeoVPixelFactory.h:43
PixelGeometryManager::PixelAlTubeLowerBendY
virtual double PixelAlTubeLowerBendY()=0
GeoVPixelFactory::m_mapFPV
std::shared_ptr< std::map< std::string, GeoFullPhysVol * > > m_mapFPV
Definition: GeoVPixelFactory.h:47
PixelGeometryManager::PixelModuleEtaFromIndex
virtual int PixelModuleEtaFromIndex(int index)=0
PixelGeometryManager::PixelOmegaEndY
virtual double PixelOmegaEndY()=0
BchCleanup.mgr
mgr
Definition: BchCleanup.py:294
cm
const double cm
Definition: Simulation/ISF/ISF_FastCaloSim/ISF_FastCaloSimParametrization/tools/FCAL_ChannelMap.cxx:25
GeoPixelLadderServices::BuildOmega
GeoVPhysVol * BuildOmega()
Definition: GeoPixelLadderServices.cxx:156
PixelGeometryManager::PixelAlTubeUpperBendRadius
virtual double PixelAlTubeUpperBendRadius()=0
GeoPixelLadderServices::m_envelopeCornerC2
GeoTrf::Vector3D m_envelopeCornerC2
Definition: GeoPixelLadderServices.h:45
PixelGeometryManager::PixelOmegaLowerBendRadius
virtual double PixelOmegaLowerBendRadius()=0
lumiFormat.i
int i
Definition: lumiFormat.py:85
GeoPixelLadderServices::Build
virtual GeoVPhysVol * Build() override
Definition: GeoPixelLadderServices.cxx:102
PixelGeometryManager::PixelAlTubeUpperBendX
virtual double PixelAlTubeUpperBendX()=0
GeoPixelLadderServices::m_xOffset
double m_xOffset
Definition: GeoPixelLadderServices.h:47
makeTRTBarrelCans.y2
tuple y2
Definition: makeTRTBarrelCans.py:18
Amg::Transform3D
Eigen::Affine3d Transform3D
Definition: GeoPrimitives.h:46
PixelGeometryManager::PixelConnectorWidthZ
virtual double PixelConnectorWidthZ(int index)=0
InDetMaterialManager::getMaterialForVolume
const GeoMaterial * getMaterialForVolume(const std::string &materialName, double volume, const std::string &newName="")
Create and get material with a density calculated to give weight in predefined weight table.
Definition: InDetMaterialManager.cxx:460
test_pyathena.parent
parent
Definition: test_pyathena.py:15
add
bool add(const std::string &hname, TKey *tobj)
Definition: fastadd.cxx:55
PixelGeometryManager::PixelNModule
virtual int PixelNModule()=0
PixelGeometryManager::PixelAlTubeLowerBendRadius
virtual double PixelAlTubeLowerBendRadius()=0
PixelGeometryManager::PixelConnectorPosY
virtual double PixelConnectorPosY(int index)=0
GeoPixelFluid::posY
double posY()
Definition: GeoPixelFluid.h:23
PixelGeometryManager::PixelConnectorPosX
virtual double PixelConnectorPosX(int index)=0
PixelGeometryManager::PixelAlTubeLength
virtual double PixelAlTubeLength()=0
PixelGeometryManager::PixelConnectorWidthX
virtual double PixelConnectorWidthX(int index)=0
GeoPixelFluid::posZ
double posZ()
Definition: GeoPixelFluid.h:24
PixelGeometryManager::PixelAlTubeUpperBendY
virtual double PixelAlTubeUpperBendY()=0
GeoPixelLadderServices::m_ladderServicesLV
GeoIntrusivePtr< const GeoLogVol > m_ladderServicesLV
Definition: GeoPixelLadderServices.h:39
grepfile.sep
sep
Definition: grepfile.py:38
RCU::Check::typeNum
const int typeNum
Definition: Assert.h:75
InDetDD::PixelDetectorManager
Definition: PixelDetectorManager.h:47
PixelGeometryManager::PixelOmegaLowerBendX
virtual double PixelOmegaLowerBendX()=0
PixelGeometryManager::PixelOmegaGluePosZ
virtual double PixelOmegaGluePosZ(int index)=0
GeoVPixelFactory::m_epsilon
const double m_epsilon
Definition: GeoVPixelFactory.h:49
GeoPixelLadderServices::m_envelopeCornerA2
GeoTrf::Vector3D m_envelopeCornerA2
Definition: GeoPixelLadderServices.h:43
PixelGeometryManager::PixelLadderLength
virtual double PixelLadderLength()=0
GeoPixelCable::numElements
int numElements() const
Definition: GeoPixelCable.h:19
PixelGeometryManager::PixelOmegaGlueStartY
virtual double PixelOmegaGlueStartY(int index)=0
GeoVPixelFactory::m_mat_mgr
InDetMaterialManager * m_mat_mgr
Definition: GeoVPixelFactory.h:44
PixelGeometryManager::PixelOmegaWallThickness
virtual double PixelOmegaWallThickness()=0
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
GeoPixelCable.h
python.SystemOfUnits.mm
int mm
Definition: SystemOfUnits.py:83
PixelGeometryManager::PixelOmegaUpperBendX
virtual double PixelOmegaUpperBendX()=0
PixelGeometryManager::PixelPigtailEndY
virtual double PixelPigtailEndY()=0
GeoPixelFluid.h
PixelGeometryManager::PixelConnectorWidthY
virtual double PixelConnectorWidthY(int index)=0
PixelGeometryManager::PixelPigtailPosZ
virtual double PixelPigtailPosZ()=0
PixelGeometryManager::PixelOmegaUpperBendRadius
virtual double PixelOmegaUpperBendRadius()=0
GeoPixelLadderServices::m_yOffset
double m_yOffset
Definition: GeoPixelLadderServices.h:48
GeoPixelLadderServices::BuildGlue
void BuildGlue(GeoPhysVol *parent)
Definition: GeoPixelLadderServices.cxx:321
GeoVPixelFactory
This is the base class for all the pieces of the Pixel detector.
Definition: GeoVPixelFactory.h:31
PixelGeometryManager::PixelAlTubeWallThickness
virtual double PixelAlTubeWallThickness()=0
PixelGeometryManager::PixelOmegaLowerBendY
virtual double PixelOmegaLowerBendY()=0
PixelGeometryManager::getMaterialName
virtual std::string getMaterialName(const std::string &volumeName, int layerdisk=0, int typenum=0)=0
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
PixelGeometryManager::PixelOmegaGlueStartX
virtual double PixelOmegaGlueStartX(int index)=0
InDetMaterialManager::getMaterial
const GeoMaterial * getMaterial(const std::string &materialName)
Get material. First looks for locally defined material and if not found looks in GeoModel material ma...
Definition: InDetMaterialManager.cxx:96
PixelGeometryManager::PixelAlTubeLowerBendX
virtual double PixelAlTubeLowerBendX()=0
GeoPixelCable::getStackOffset
double getStackOffset()
Definition: GeoPixelCable.cxx:99
GeoPixelFluid::posX
double posX()
Definition: GeoPixelFluid.h:22
length
double length(const pvec &v)
Definition: FPGATrackSimLLPDoubletHoughTransformTool.cxx:26
PixelGeometryManager::PixelOmegaUpperBendY
virtual double PixelOmegaUpperBendY()=0
PixelGeometryManager::PixelLadderCableOffsetX
virtual double PixelLadderCableOffsetX()=0
PixelGeometryManager::PixelModuleZPosition
virtual double PixelModuleZPosition(int)=0
PixelGeometryManager::GetLD
virtual int GetLD()=0
Generated on Sun Dec 22 2024 21:10:50 for ATLAS Offline Software by   doxygen 1.8.18