Placing wires in the designated wireGroups for easy retrieval later, first and last are placed separately.
Stacking strip, wireGroup and pad layers Redefining the gasgap center for L3 (diamond) sector: Unlike the run 3 (legacy) code, we are NOT offsetting the gasgap center in the local y coordinate to match the base of the cutout. As a result, localToGlobal transformations and the local y coordinate of the pads, strips and wireGroups will not match the run 3. Additionally, the global y coordinate of the wireGroup is shifted by ~25 mm due to this center redefinition. The redefinition will not affect reconstruction and digitization, and will simplify the sTgc geometry code.
65 {
67 <<std::endl<<std::endl<<
m_geoUtilTool->dumpVolume(define.physVol));
68 const GeoShape* shape =
m_geoUtilTool->extractShape(define.physVol);
69 if (!shape) {
71 return StatusCode::FAILURE;
72 }
76
77 define.sHalfChamberLength = modPars.shortWidth;
78 define.lHalfChamberLength = modPars.longWidth;
79 define.halfChamberHeight = modPars.halfHeight;
80 define.halfChamberTck = modPars.thickness;
81 define.yCutout = modPars.yCutOut;
82
83 ATH_MSG_DEBUG(
"chamber length (L/S) is: " << 2*define.lHalfChamberLength <<
"/"
84 << 2*define.sHalfChamberLength << " chamber height is: "
85 << 2*define.halfChamberHeight << " chamber thickness is: " << 2*define.halfChamberTck << " yCutout value is: " << define.yCutout);
86
88 std::vector<physVolWithTrans> allGasGaps =
m_geoUtilTool->findAllLeafNodesByName(define.physVol,
"sTgcGas");
89 if (allGasGaps.empty()) {
92 return StatusCode::FAILURE;
93 }
94
96 define.numLayers = allGasGaps.size();
98 FactoryCache::ParamBookTable::const_iterator parBookItr = factoryCache.parameterBook.find(define.chambDesign);
99 if (parBookItr == factoryCache.parameterBook.end()) {
100 ATH_MSG_FATAL(
"The chamber "<<define.chambDesign<<
" is not part of the WSTGC table");
101 return StatusCode::FAILURE;
102 }
103
104 const wSTGCTable& paramBook{parBookItr->second};
105
106 define.
gasTck = paramBook.gasTck;
111 WireDesignPtr wireGroupDesign = std::make_unique<WireGroupDesign>();
113
115
117 double firstStripPos = -gapPars.halfHeight + paramBook.firstStripPitch[
gasGap] - 0.5 * paramBook.stripPitch;
118 define.firstStripPitch = paramBook.firstStripPitch;
119 ATH_MSG_DEBUG(
"FirstStripPos for stripPitch: " << paramBook.firstStripPitch[gasGap] <<
" is: " <<
Amg::toString(firstStripPos, 2) <<
" and the half height is: " << gapPars.halfHeight);
120
122 unsigned int numWireGroups = paramBook.numWireGroups[
gasGap];
124 double firstWirePos = paramBook.firstWirePos[
gasGap];
125
126 if(gapPars.yCutOut) {
128 stripDesign->defineDiamond(gapPars.shortWidth, gapPars.longWidth, gapPars.halfHeight + 0.001, paramBook.yCutoutCathode);
129 ATH_MSG_VERBOSE(
"The yCutout of the active area is: " << gapPars.yCutOut);
130 stripDesign->defineStripLayout(
Amg::Vector2D{firstStripPos, 0.},
131 paramBook.stripPitch, paramBook.stripWidth, paramBook.numStrips);
134 wireGroupDesign->defineDiamond(0.5*paramBook.sPadLength, 0.5*paramBook.lPadLength, gapPars.halfHeight, paramBook.yCutoutCathode);
135 wireGroupDesign->flipTrapezoid();
136 wireGroupDesign->defineStripLayout(
Amg::Vector2D{firstWirePos, 0.},
137 paramBook.wirePitch,
138 paramBook.wireWidth,
139 numWireGroups);
140 ATH_MSG_VERBOSE(
"Created new diamond wireGroup design "<<(*wireGroupDesign));
141
143 padDesign->defineDiamond(0.5*paramBook.sPadLength, 0.5*paramBook.lPadLength, gapPars.halfHeight, paramBook.yCutoutCathode);
144 padDesign->flipTrapezoid();
146 } else {
148 stripDesign->defineTrapezoid(gapPars.shortWidth, gapPars.longWidth, gapPars.halfHeight + 0.01);
149 stripDesign->defineStripLayout(
Amg::Vector2D{firstStripPos, 0.},
150 paramBook.stripPitch, paramBook.stripWidth, paramBook.numStrips);
152
154 wireGroupDesign->defineTrapezoid(0.5*paramBook.sPadLength, 0.5*paramBook.lPadLength, gapPars.halfHeight);
155 wireGroupDesign->flipTrapezoid();
156 wireGroupDesign->defineStripLayout(
Amg::Vector2D{firstWirePos, 0.},
157 paramBook.wirePitch,
158 paramBook.wireWidth,
159 numWireGroups);
161
163 padDesign->defineTrapezoid(0.5*paramBook.sPadLength, 0.5*paramBook.lPadLength, gapPars.halfHeight);
164 padDesign->flipTrapezoid();
166 }
167
169 wireGroupDesign->declareGroup(paramBook.firstWireGroupWidth[gasGap]);
170 for (
uint wireGr=2; wireGr<numWireGroups; wireGr++){
171 wireGroupDesign->declareGroup(paramBook.wireGroupWidth);
172 }
173 unsigned int lastWireGroup = (paramBook.numWires[
gasGap] - wireGroupDesign->nAllWires());
174 wireGroupDesign->declareGroup(lastWireGroup);
175 wireGroupDesign->defineWireCutout(paramBook.wireCutout[gasGap]);
176
178 double beamlineRadius = (define.physVol->getAbsoluteTransform() * gapVol.transform).translation().perp();
179 padDesign->defineBeamlineRadius(beamlineRadius);
181 padDesign->definePadRow(paramBook.firstPadPhiDivision[gasGap],
182 paramBook.numPadPhi[gasGap],
183 paramBook.anglePadPhi,
184 paramBook.PadPhiShift[gasGap]);
185 padDesign->definePadColumn(paramBook.firstPadHeight[gasGap],
186 paramBook.numPadEta[gasGap],
187 paramBook.padHeight[gasGap]);
188
199 stripDesign = (*factoryCache.stripDesigns.emplace(stripDesign).first);
200 wireGroupDesign = (*factoryCache.wireGroupDesigns.emplace(wireGroupDesign).first);
201 StripLayer stripLayer{factoryCache.trfNodeMaker.makeTransform(gapVol.transform
204 stripDesign, wireGroupDesign,
206 stripLayer.flipPhiRotation();
208 define.stripLayers.push_back(std::move(stripLayer));
209
211 padDesign = (*factoryCache.padDesigns.emplace(padDesign).first);
212 StripLayer padLayer(factoryCache.trfNodeMaker.makeTransform(gapVol.transform
216 padLayer.flipPhiRotation();
218 define.padLayers.push_back(std::move(padLayer));
219 }
220 return StatusCode::SUCCESS;
221}
static IdentifierHash createHash(const unsigned gasGap, const unsigned channelType, const unsigned channel, const unsigned wireInGrp=0)
Create a measurement hash from the Identifier fields.
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Amg::Transform3D getRotateZ3D(double angle)
get a rotation transformation around Z-axis
Eigen::Matrix< double, 2, 1 > Vector2D
Amg::Transform3D getRotateY3D(double angle)
get a rotation transformation around Y-axis
GeoModel::TransientConstSharedPtr< WireGroupDesign > WireDesignPtr
IMuonGeoUtilityTool::physVolWithTrans physVolWithTrans
GeoModel::TransientConstSharedPtr< StripDesign > StripDesignPtr
GeoModel::TransientConstSharedPtr< PadDesign > PadDesignPtr
1.14.0