ATLAS Offline Software
SCTGeometryManagerTool.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include "AthContainers/DataVector.h"
6 
7 #include "SCT_ReadoutGeometry/SCT_DetectorManager.h"
8 #include "InDetReadoutGeometry/SiDetectorElement.h"
9 
10 #include "InDetIdentifier/SCT_ID.h"
11 
12 #include "TrkAlignInterfaces/IAlignModuleTool.h"
13 
14 #include "TrkAlignEvent/AlignModule.h"
15 #include "TrkAlignEvent/AlignModuleList.h"
16 #include "TrkAlignEvent/AlignPar.h"
17 
18 
19 
20 #include "InDetAlignGeomTools/SCTGeometryManagerTool.h"
21 
22 
23 #include <iostream>
24 
25 using namespace InDetDD;
26 
27 namespace InDet {
28 
29  //________________________________________________________________________
30  SCTGeometryManagerTool::SCTGeometryManagerTool(const std::string& type,
31  const std::string& name,
32  const IInterface * parent)
33  : AthAlgTool(type,name,parent)
34  , m_detManager(nullptr)
35  , m_idHelper()
36  , m_alignModuleTool("Trk::AlignModuleTool/AlignModuleTool")
37  , m_idHashToAlignModuleMaps(Trk::AlignModule::NDetectorTypes,(Trk::AlignModuleList*)nullptr)
38  , m_alignParList(nullptr)
39  , m_fullAlignParList(nullptr)
40  {
41  declareInterface<IGeometryManagerTool>(this);
42  declareProperty("AlignModuleTool", m_alignModuleTool);
43 
44  // Pixel joboptions
45  declareProperty("AlignBarrel", m_alignBarrel = true);
46  declareProperty("AlignEndcaps", m_alignEndcaps = true);
47 
48  declareProperty("AlignBarrelX", m_alignBarrelX = true);
49  declareProperty("AlignBarrelY", m_alignBarrelY = true);
50  declareProperty("AlignBarrelZ", m_alignBarrelZ = true);
51  declareProperty("AlignBarrelRotX", m_alignBarrelRotX = true);
52  declareProperty("AlignBarrelRotY", m_alignBarrelRotY = true);
53  declareProperty("AlignBarrelRotZ", m_alignBarrelRotZ = true);
54 
55  declareProperty("AlignEndcapX", m_alignEndcapX = true);
56  declareProperty("AlignEndcapY", m_alignEndcapY = true);
57  declareProperty("AlignEndcapZ", m_alignEndcapZ = true);
58  declareProperty("AlignEndcapRotX", m_alignEndcapRotX = true);
59  declareProperty("AlignEndcapRotY", m_alignEndcapRotY = true);
60  declareProperty("AlignEndcapRotZ", m_alignEndcapRotZ = true);
61 
62  declareProperty("SetSigmaBarrelX", m_sigmaBarrelX = 1.);
63  declareProperty("SetSigmaBarrelY", m_sigmaBarrelY = 1.);
64  declareProperty("SetSigmaBarrelZ", m_sigmaBarrelZ = 1.);
65  declareProperty("SetSigmaBarrelRotX", m_sigmaBarrelRotX = 0.001);
66  declareProperty("SetSigmaBarrelRotY", m_sigmaBarrelRotY = 0.001);
67  declareProperty("SetSigmaBarrelRotZ", m_sigmaBarrelRotZ = 0.001);
68 
69  declareProperty("SetSigmaEndcapX", m_sigmaEndcapX = 1.);
70  declareProperty("SetSigmaEndcapY", m_sigmaEndcapY = 1.);
71  declareProperty("SetSigmaEndcapZ", m_sigmaEndcapZ = 1.);
72  declareProperty("SetSigmaEndcapRotX", m_sigmaEndcapRotX = 0.001);
73  declareProperty("SetSigmaEndcapRotY", m_sigmaEndcapRotY = 0.001);
74  declareProperty("SetSigmaEndcapRotZ", m_sigmaEndcapRotZ = 0.001);
75 
76  declareProperty("SetSoftCutBarrelX", m_softcutBarrelX = 1.);
77  declareProperty("SetSoftCutBarrelY", m_softcutBarrelY = 1.);
78  declareProperty("SetSoftCutBarrelZ", m_softcutBarrelZ = 1.);
79  declareProperty("SetSoftCutBarrelRotX", m_softcutBarrelRotX = 1.);
80  declareProperty("SetSoftCutBarrelRotY", m_softcutBarrelRotY = 1.);
81  declareProperty("SetSoftCutBarrelRotZ", m_softcutBarrelRotZ = 1.);
82 
83  declareProperty("SetSoftCutEndcapX", m_softcutEndcapX = 1.);
84  declareProperty("SetSoftCutEndcapY", m_softcutEndcapY = 1.);
85  declareProperty("SetSoftCutEndcapZ", m_softcutEndcapZ = 1.);
86  declareProperty("SetSoftCutEndcapRotX", m_softcutEndcapRotX = 1.);
87  declareProperty("SetSoftCutEndcapRotY", m_softcutEndcapRotY = 1.);
88  declareProperty("SetSoftCutEndcapRotZ", m_softcutEndcapRotZ = 1.);
89 
90  // defines alignment level
91  declareProperty("AlignmentLevel", m_alignLevel = 1);
92  declareProperty("AlignmentLevelBarrel", m_alignLevelBarrel = -1);
93  declareProperty("AlignmentLevelEndcaps", m_alignLevelEndcaps = -1);
94 
95  declareProperty("doModuleSelection", m_doModuleSelection = false);
96  declareProperty("ModuleSelection", m_moduleSelection);
97 
98  declareProperty("UseOldL2", m_useOldL2 = false);
99  declareProperty("DumpGeometry", m_dumpGeometry = true);
100 
101  m_hashCounter = 0;
102  m_logStream = nullptr;
103  }
104 
105  //________________________________________________________________________
106  SCTGeometryManagerTool::~SCTGeometryManagerTool()
107  {
108  ATH_MSG_DEBUG("deleting alignModuleList");
109  for (int i=0;i<(int)m_alignModuleList.size();i++)
110  delete m_alignModuleList[i];
111  m_alignModuleList.clear();
112 
113  ATH_MSG_DEBUG("deleting fullAlignParList");
114  delete m_fullAlignParList;
115  ATH_MSG_DEBUG("deleting alignParList");
116  delete m_alignParList;
117  }
118 
119  //________________________________________________________________________
120  StatusCode SCTGeometryManagerTool::initialize()
121  {
122  ATH_MSG_DEBUG("initialize() of SCTGeometryManagerTool");
123 
124  // retrieve AlignModuleTool
125  if ( m_alignModuleTool.retrieve().isFailure() ) {
126  msg(MSG::FATAL)<<"Could not get " << m_alignModuleTool << endmsg;
127  return StatusCode::FAILURE;
128  }
129  else
130  ATH_MSG_INFO("Retrieved " << m_alignModuleTool);
131 
132  // retrieve SCT helper
133  if ( detStore()->retrieve(m_idHelper).isFailure() ) {
134  msg(MSG::FATAL) << " Cannot retrieve SCT Helper " << endmsg;
135  return StatusCode::FAILURE;
136  }
137  else
138  ATH_MSG_INFO("retrieved Silicon SCT Helper");
139 
140  // retrieve SCT detector manager
141  if ( detStore()->retrieve(m_detManager,m_stripDetManagerName).isFailure() ) {
142  msg(MSG::FATAL) << " Cannot retrieve SCT Detector Manager " << endmsg;
143  return StatusCode::FAILURE;
144  }
145  else
146  ATH_MSG_INFO("retrieved SCT Detector Manager");
147 
148  // dump module selection
149  if(m_doModuleSelection && msgLvl(MSG::INFO)) {
150  msg(MSG::INFO)<<"Creating geometry for selected "<<m_moduleSelection.size()<<" modules:"<<endmsg;
151  for(unsigned int i=0;i<m_moduleSelection.size();i++)
152  msg(MSG::INFO)<<" "<<i<<". "<<m_moduleSelection.at(i)<<endmsg;
153  }
154 
155  // check the allowed geometry levels
156  if(!checkAlignLevel())
157  return StatusCode::FAILURE;
158 
159  return StatusCode::SUCCESS;
160  }
161 
162  //________________________________________________________________________
163  StatusCode SCTGeometryManagerTool::finalize()
164  {
165  ATH_MSG_DEBUG("finalize() of SCTGeometryManagerTool");
166 
167  return StatusCode::SUCCESS;
168  }
169 
170  //________________________________________________________________________
171  void SCTGeometryManagerTool::setAlignLevel(int level)
172  {
173  m_alignLevelBarrel = m_alignLevelEndcaps = -1;
174  m_alignLevel = level;
175  }
176 
177  //________________________________________________________________________
178  bool SCTGeometryManagerTool::checkAlignLevel()
179  {
180  if(m_alignLevel == 0) {
181  ATH_MSG_INFO("Alignment level for SCT is "<<m_alignLevel);
182  return true;
183  }
184 
185  if(m_alignLevel == 2 && m_useOldL2) {
186  ATH_MSG_INFO("Setting up OLD L2 alignment level for SCT");
187  m_alignLevelBarrel = m_alignLevel;
188  m_alignLevelEndcaps = m_alignLevel;
189  return true;
190  }
191 
192  if(m_alignLevelBarrel == -1)
193  m_alignLevelBarrel = m_alignLevel;
194 
195  if(m_alignLevelEndcaps == -1)
196  m_alignLevelEndcaps = m_alignLevel;
197 
198  bool ok = checkAlignLevelBarrel() && checkAlignLevelEndcaps();
199 
200  // if the alignment level for barrel and EC is set correctly
201  // reset the global alignment level
202  if(ok)
203  m_alignLevel = -1;
204 
205  return ok;
206  }
207 
208  //________________________________________________________________________
209  bool SCTGeometryManagerTool::checkAlignLevelBarrel()
210  {
211  if(!m_alignBarrel)
212  return true;
213 
214  // check whether geometry level is allowed
215  switch(m_alignLevelBarrel) {
216  case 1: case 2: case 27: case 3:
217  ATH_MSG_INFO("Alignment level for SCT Barrel is "<<m_alignLevelBarrel);
218  return true;
219  default:
220  msg(MSG::FATAL)<<"Alignment level "<<m_alignLevelBarrel<<" does not exist for SCT Barrel"<<endmsg;
221  return false;
222  }
223  }
224 
225  //________________________________________________________________________
226  bool SCTGeometryManagerTool::checkAlignLevelEndcaps()
227  {
228  if(!m_alignEndcaps)
229  return true;
230 
231  // check whether geometry level is allowed
232  switch(m_alignLevelEndcaps) {
233  case 1: case 2: case 25: case 3:
234  ATH_MSG_INFO("Alignment level for SCT Endcaps is "<<m_alignLevelBarrel);
235  return true;
236  default:
237  msg(MSG::FATAL)<<"Alignment level "<<m_alignLevelEndcaps<<" does not exist for SCT Endcaps"<<endmsg;
238  return false;
239  }
240  }
241 
242  //________________________________________________________________________
243  int SCTGeometryManagerTool::ReadGeometry(int solveLevel)
244  {
245  ATH_MSG_DEBUG("in ReadGeometry() solveLevel="<<solveLevel);
246 
247  // set pointers
248  m_idHashToAlignModuleMapsPtr = &m_idHashToAlignModuleMaps;
249  m_alignModuleListPtr = &m_alignModuleList;
250 
251  // build alignment geometry
252  buildGeometry();
253 
254  // now set the alignment parameters
255  // first prepare the parameter lists
256  m_alignParList = new DataVector< DataVector<Trk::AlignPar> >(SG::OWN_ELEMENTS);
257  m_fullAlignParList = new DataVector< DataVector<Trk::AlignPar> >(SG::OWN_ELEMENTS);
258  // loop over modules
259  ATH_MSG_DEBUG("Adding module parameters to modules");
260  std::vector<Trk::AlignModule *>::const_iterator imod = m_alignModuleList.begin();
261  std::vector<Trk::AlignModule *>::const_iterator imod_end = m_alignModuleList.end();
262  for( ; imod!=imod_end ; ++imod)
263  addModuleParameters(*imod,m_fullAlignParList,m_alignParList);
264 
265  // set alignModuleList and hash table in the alignModuleTool
266  m_alignModuleTool->setAlignModules(&m_alignModuleList, &m_idHashToAlignModuleMaps);
267  ATH_MSG_DEBUG(" geometry set in m_alignModuleTool");
268 
269  // set alignPar lists in the alignModuleTool
270  ATH_MSG_DEBUG(" alignParList = "<<m_alignParList);
271  ATH_MSG_DEBUG(" fullAlignParList = "<<m_fullAlignParList);
272  m_alignModuleTool->setAlignParLists(m_fullAlignParList, m_alignParList);
273  ATH_MSG_DEBUG(" AlignParLists set in m_alignModuleTool");
274 
275  // dump summary about the geometry setup
276  if (m_dumpGeometry)
277  dumpGeometry();
278 
279  int nDoF= m_alignModuleTool->nAlignParameters();
280  ATH_MSG_INFO("Total number of degrees of freedom: "<<nDoF);
281 
282  return nDoF;
283  }
284 
285  //_______________________________________________________________________
286  void SCTGeometryManagerTool::buildGeometry()
287  {
288  ATH_MSG_INFO("Preparing the SCT geometry");
289 
290  if(m_alignLevel == 0)
291  buildL0();
292  else if (m_alignLevel == 2 && m_useOldL2)
293  buildL2Old();
294  else {
295 
296  if(m_alignBarrel)
297  switch(m_alignLevelBarrel) {
298  case 1:
299  buildL1Barrel();
300  break;
301  case 2:
302  buildL2Barrel();
303  break;
304  case 27:
305  buildL27Barrel();
306  break;
307  case 3:
308  buildL3Barrel();
309  break;
310  default:
311  break;
312  }
313 
314  if(m_alignEndcaps)
315  switch(m_alignLevelEndcaps) {
316  case 1:
317  buildL1Endcaps();
318  break;
319  case 2:
320  buildL2Endcaps();
321  break;
322  case 25:
323  buildL25Endcaps();
324  break;
325  case 3:
326  buildL3Endcaps();
327  break;
328  default:
329  break;
330  }
331 
332  }
333  }
334 
335  //_______________________________________________________________________
336  void SCTGeometryManagerTool::buildL0()
337  {
338  ATH_MSG_INFO("Preparing the SCT geometry for L0: the whole SCT");
339  // ===================================================
340  // Level 0 is just one module containing the whole SCT
341 
342  Trk::AlignModule * sct = new Trk::AlignModule(this);
343  sct->setIdHash(getNextIDHash());
344  sct->setName("SCT");
345 
346  // use the identifier of the SCT barrel for the whole SCT at L0
347  sct->setIdentifier(m_idHelper->wafer_id(0,0,0,0,0));
348 
349  ATH_MSG_DEBUG("Created module "<<sct->name()<<" idHash: "<<sct->identifyHash()<<" identifier: "<<sct->identify());
350 
351  // for L0 alignment the alignment frame is equal to the global frame
352  // and since the SCT detector element positions are also stored
353  // in the global frame in DB, transform is identity
354  const Amg::Transform3D transform = Amg::Transform3D::Identity();
355 
356  // get maximum number of elements from the helper
357  unsigned int sctmaxHash = m_idHelper->wafer_hash_max();
358  ATH_MSG_DEBUG("maxHash for the SCT: "<<sctmaxHash);
359 
360  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
361  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(sctmaxHash),nullptr);
362  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
363 
364  // ================================================================
365  // loop over SCT elements and add them to respective alignModules
366  // ================================================================
367  for (unsigned int index = 0; index < sctmaxHash; index++) {
368  IdentifierHash idHash = index;
369  Identifier id = m_idHelper->wafer_id(idHash);
370 
371  ATH_MSG_DEBUG(" SCT DetectorElement idhash: "<<index);
372  ATH_MSG_DEBUG(" DetectorElement id: "<<id);
373 
374  // get the element via hash
375  const SiDetectorElement * element = m_detManager->getDetectorElement(id);
376  if (element) {
377  // add element to respective AlignModule
378  // add to the sct barrel structure
379  if(msgLvl(MSG::DEBUG)) {
380  if (m_idHelper->is_barrel(id))
381  msg(MSG::DEBUG)<<"... SCT barrel element"<<endmsg;
382  else
383  msg(MSG::DEBUG)<<"... SCT endcap element"<<endmsg;
384  }
385  sct->addDetElement(Trk::AlignModule::SCT,element,transform);
386 
387  // and fill the corresponding map
388  (*sctIdHashMap)[idHash] = sct;
389  }
390  else
391  ATH_MSG_DEBUG("No SCT detector with id: "<<id);
392  }
393 
394  // add created module to the geometry
395  m_alignModuleListPtr->push_back(sct);
396 
397  ATH_MSG_DEBUG("SCT L0 module successfully added to the geometry");
398  }
399 
400  //_______________________________________________________________________
401  void SCTGeometryManagerTool::buildL1Barrel()
402  {
403  // ===========================
404  // BARREL
405  // ===========================
406  ATH_MSG_INFO("Preparing the SCT Barrel geometry for L1: barrel + 2 endcaps");
407 
408  // ========================================
409  // Level 1 is whole SCT barrel
410 
411  // get maximum number of elements from the helper
412  unsigned int maxHash = m_idHelper->wafer_hash_max();
413  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
414 
415  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
416  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
417  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
418 
419  // use the 'real' identifier of the Pixel
420  // build barrel module
421  Trk::AlignModule * mod = new Trk::AlignModule(this);
422  mod->setIdHash(getNextIDHash());
423  mod->setName("SCT Barrel");
424  mod->setIdentifier(m_idHelper->wafer_id(0,0,0,0,0));
425 
426  // check if selected
427  if(!m_alignBarrel || !moduleSelected(mod)) {
428  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
429  delete mod;
430  return;
431  }
432 
433  // fill module with detector elements
434  ATH_MSG_DEBUG("Building module "<<mod->name());
435 
436  const Amg::Transform3D transform = Amg::Transform3D::Identity();
437 
438  for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
439  if (!m_detManager->numerology().useLayer(iLayer))
440  ATH_MSG_INFO(" Layer "<<iLayer<<" not present");
441 
442  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
443  ATH_MSG_DEBUG("iPhi "<<iPhi);
444  for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
445  ATH_MSG_DEBUG("iEta "<<iEta);
446 
447  if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
448  // iEta=0 is not defined for the SCT
449  continue;
450 
451  for (int is = 0; is < 2; is++) { // module side
452  const SiDetectorElement* element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
453 
454  if (element) {
455  // get element location for debugging
456  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
457  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
458  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
459 
460  // add element to the AlignModule
461  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
462  // and fill the corresponding map
463  (*sctIdHashMap)[element->identifyHash()] = mod;
464  }
465  }
466  }
467  }
468  }
469 
470  m_alignModuleListPtr->push_back(mod);
471 
472  ATH_MSG_DEBUG("SCT Barrel successfully added to the geometry");
473  }
474 
475  //_______________________________________________________________________
476  void SCTGeometryManagerTool::buildL1Endcaps()
477  {
478  // ===========================
479  // ENDCAPs
480  // ===========================
481  ATH_MSG_INFO("Preparing the SCT Endcap geometry for L1: 2 endcaps");
482 
483  // ========================================
484  // Level 1 is 2 SCT endcaps
485 
486  // get maximum number of elements from the helper
487  unsigned int maxHash = m_idHelper->wafer_hash_max();
488  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
489 
490  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
491  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
492  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
493 
494  for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
495  int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
496 
497  Trk::AlignModule * mod = new Trk::AlignModule(this);
498  mod->setIdHash(getNextIDHash());
499  mod->setIdentifier(m_idHelper->wafer_id(iSide,0,0,0,0));
500  if(iEndcapIndex == 0)
501  mod->setName("SCT EndCap A");
502  else
503  mod->setName("SCT EndCap C");
504 
505  if(!m_alignEndcaps || !moduleSelected(mod)) {
506  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
507  delete mod;
508  continue;
509  }
510 
511  ATH_MSG_DEBUG("Building module "<<mod->name());
512 
513  const Amg::Transform3D transform = Amg::Transform3D::Identity();
514 
515  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
516  if (!m_detManager->numerology().useDisk(iWheel))
517  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
518 
519  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
520  ATH_MSG_DEBUG("iEta "<<iEta);
521  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
522  ATH_MSG_DEBUG("iPhi "<<iPhi);
523  for(int is=0; is<2; is++) { // module side
524 
525  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
526  if (element) {
527  // get element location for debugging
528  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
529  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
530  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
531 
532  // add element to the AlignModule
533  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
534  // and fill the corresponding map
535  (*sctIdHashMap)[element->identifyHash()] = mod;
536  }
537  }
538  }
539  }
540  }
541 
542  // for endcap we move the CoG
543  Amg::Translation3D translation(mod->centerOfGravity());
544  ATH_MSG_DEBUG("Endcap:" << iSide << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
545  Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
546 
547  mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
548 
549  m_alignModuleListPtr->push_back(mod);
550  }
551 
552  ATH_MSG_DEBUG("SCT End-caps successfully added to the geometry");
553  }
554 
555  //_______________________________________________________________________
556  void SCTGeometryManagerTool::buildL2Barrel()
557  {
558  // ========================================
559  // BARREL
560  ATH_MSG_INFO("Preparing the SCT Barrel geometry for L2: 4 layers");
561 
562  // ========================================
563  // get all modules for Level 2 alignment of the Barrel
564  // Level 2 is 4 SCT barrel layers
565 
566  const Amg::Transform3D transform = Amg::Transform3D::Identity();
567 
568  unsigned int maxHash = m_idHelper->wafer_hash_max();
569  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
570 
571  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
572  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
573  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
574 
575  for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
576  if (!m_detManager->numerology().useLayer(iLayer))
577  ATH_MSG_INFO(" Layer "<<iLayer<<" not present");
578 
579  // create the AlignModule
580  Trk::AlignModule * mod = new Trk::AlignModule(this);
581  mod->setIdHash(getNextIDHash());
582  mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, 0, 0, 0));
583 
584  std::stringstream name;
585  name <<"SCT/Barrel/Layer_"<<iLayer;
586  mod->setName(name.str());
587 
588  if(!m_alignBarrel || !moduleSelected(mod)) {
589  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
590  delete mod;
591  continue;
592  }
593 
594  ATH_MSG_DEBUG("Building module "<<mod->name());
595 
596  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
597  ATH_MSG_DEBUG("iPhi "<<iPhi);
598  for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
599  ATH_MSG_DEBUG("iEta "<<iEta);
600  if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
601  // iEta=0 is not defined for the SCT
602  continue;
603  for (int is = 0; is < 2; is++) { // module side
604 
605  const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
606  if (element) {
607  // get element location for debugging
608  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
609  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
610  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
611 
612  // add element to the AlignModule
613  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
614  // and fill the corresponding map
615  (*sctIdHashMap)[element->identifyHash()] = mod;
616  }
617  }
618  }
619  }
620 
621  m_alignModuleListPtr->push_back(mod);
622  }
623  }
624 
625  //_______________________________________________________________________
626  void SCTGeometryManagerTool::buildL27Barrel()
627  {
628  // ========================================
629  // BARREL
630  ATH_MSG_INFO("Preparing the SCT Barrel geometry for L27: 176 staves (ladders)");
631 
632  if(!m_alignBarrel) {
633  ATH_MSG_DEBUG("Not aligning barrel");
634  return;
635  }
636 
637  // ========================================
638  // get all modules for Level 27 alignment of the Barrel
639  // 176 staves (ladders)
640 
641  const Amg::Transform3D transform = Amg::Transform3D::Identity();
642 
643  unsigned int maxHash = m_idHelper->wafer_hash_max();
644  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
645 
646  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
647  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
648  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
649 
650  for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
651  if (!m_detManager->numerology().useLayer(iLayer))
652  ATH_MSG_INFO(" Layer "<<iLayer<<" not present");
653  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
654  ATH_MSG_DEBUG("iPhi "<<iPhi);
655 
656  // create the AlignModule
657  Trk::AlignModule * mod = new Trk::AlignModule(this);
658  mod->setIdHash(getNextIDHash());
659 
660  // even though there is no iEta=0 module for SCT barrel, the
661  // Identifier is still valid so we use it for the stave
662  mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, iPhi, 0, 0));
663 
664  std::stringstream name;
665  name <<"SCT/Barrel/Layer_"<<iLayer<<"/PhiStave_"<<iPhi;
666  mod->setName(name.str());
667 
668  if(!moduleSelected(mod)) {
669  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
670  delete mod;
671  continue;
672  }
673 
674  ATH_MSG_DEBUG("Building module "<<mod->name());
675 
676  // rotation of the stave alignment frame
677  Amg::RotationMatrix3D rotation; rotation.setIdentity();
678 
679  for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
680  ATH_MSG_DEBUG("iEta "<<iEta);
681  if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
682  // iEta=0 is not defined for the SCT
683  continue;
684  for (int is = 0; is < 2; is++) { // module side
685 
686  const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
687  if (element) {
688  // get element location for debugging
689  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
690  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
691  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
692 
693  // add element to the AlignModule
694  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
695  // and fill the corresponding map
696  (*sctIdHashMap)[element->identifyHash()] = mod;
697  }
698 
699  // for the stave alignment frame rotation we use the one of the iEta=1
700  // non-stereo side (which is the module local frame)
701  if(iEta==1 && !element->isStereo())
702  rotation = element->moduleTransform().rotation();
703  }
704  }
705 
706  // we set the alignment frame to be the CoG of the stave
707  // with rotation being the one of the iEta=1 module set above
708  Amg::Translation3D translation(mod->centerOfGravity());
709  Amg::Transform3D localToGlobal = translation * rotation;
710 
711 
712  ATH_MSG_DEBUG("Prepared local to global transform :");
713  ATH_MSG_DEBUG(" - translation: "<<localToGlobal.translation().x()<<" "<<localToGlobal.translation().y()<<" "<<localToGlobal.translation().z());
714  ATH_MSG_DEBUG(" - rotation:");
715  ATH_MSG_DEBUG(" "<<localToGlobal.rotation()(0,0)<<" "<<localToGlobal.rotation()(0,1)<<" "<<localToGlobal.rotation()(0,2));
716  ATH_MSG_DEBUG(" "<<localToGlobal.rotation()(1,0)<<" "<<localToGlobal.rotation()(1,1)<<" "<<localToGlobal.rotation()(1,2));
717  ATH_MSG_DEBUG(" "<<localToGlobal.rotation()(2,0)<<" "<<localToGlobal.rotation()(2,1)<<" "<<localToGlobal.rotation()(2,2));
718 
719  mod->setGlobalFrameToAlignFrameTransform(localToGlobal.inverse());
720 
721  m_alignModuleListPtr->push_back(mod);
722  }
723  }
724  }
725 
726  //_______________________________________________________________________
727  void SCTGeometryManagerTool::buildL2Endcaps()
728  {
729  // ========================================
730  // END CAPS
731  ATH_MSG_INFO("Preparing the SCT Endcap geometry for L2: 2 x 9 disks");
732 
733  // ========================================
734  // get all modules for Level 2 alignment of the Endcaps
735  // Level 2 is 2 x 9 SCT Endcap disks
736 
737  const Amg::Transform3D transform = Amg::Transform3D::Identity();
738 
739  unsigned int maxHash = m_idHelper->wafer_hash_max();
740  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
741 
742  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
743  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
744  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
745 
746  for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
747  int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
748  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
749  if (!m_detManager->numerology().useDisk(iWheel))
750  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
751 
752  Trk::AlignModule * mod = new Trk::AlignModule(this);
753  mod->setIdHash(getNextIDHash());
754  mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
755 
756  std::stringstream name;
757  if(iEndcapIndex == 0)
758  name <<"SCT/EndcapA/Disk_" << iWheel;
759  else
760  name <<"SCT/EndcapC/Disk_" << iWheel;
761  mod->setName(name.str());
762 
763  if(!m_alignEndcaps || !moduleSelected(mod)) {
764  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
765  delete mod;
766  continue;
767  }
768 
769  ATH_MSG_DEBUG("Building module "<<mod->name());
770 
771  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
772  ATH_MSG_DEBUG("iEta "<<iEta);
773  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
774  ATH_MSG_DEBUG("iPhi "<<iPhi);
775  for(int is=0; is<2; is++) { // module side
776 
777  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
778  if (element) {
779  // get element location for debugging
780  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
781  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
782  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
783 
784  // add element to the AlignModule
785  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
786  // and fill the corresponding map
787  (*sctIdHashMap)[element->identifyHash()] = mod;
788  }
789  }
790  }
791  }
792 
793  // for endcap we move the CoG
794  Amg::Translation3D translation(mod->centerOfGravity());
795  ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
796  Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
797  mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
798 
799  m_alignModuleListPtr->push_back(mod);
800  }
801  }
802  }
803 
804  //_______________________________________________________________________
805  void SCTGeometryManagerTool::buildL25Endcaps()
806  {
807  // ========================================
808  // END CAPS
809  ATH_MSG_INFO("Preparing the SCT Endcap geometry for L25: 2 x 22 rings");
810 
811  if(!m_alignEndcaps) {
812  ATH_MSG_DEBUG("Not aligning endcaps");
813  return;
814  }
815 
816  // ========================================
817  // get all modules for Level 2 alignment of the Endcaps
818  // 2 x 22 SCT Endcap rings
819 
820  const Amg::Transform3D transform = Amg::Transform3D::Identity();
821 
822  unsigned int maxHash = m_idHelper->wafer_hash_max();
823  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
824 
825  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
826  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
827  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
828 
829  for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
830  int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
831  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
832  if (!m_detManager->numerology().useDisk(iWheel))
833  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
834  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
835  ATH_MSG_DEBUG("iEta "<<iEta);
836 
837  Trk::AlignModule * mod = new Trk::AlignModule(this);
838  mod->setIdHash(getNextIDHash());
839  mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,iEta,0,0));
840 
841  std::stringstream name;
842  if(iEndcapIndex == 0)
843  name <<"SCT/EndcapA/Disk_"<<iWheel<<"/Ring_"<<iEta;
844  else
845  name <<"SCT/EndcapC/Disk_"<<iWheel<<"/Ring_"<<iEta;
846  mod->setName(name.str());
847 
848  if(!moduleSelected(mod)) {
849  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
850  delete mod;
851  continue;
852  }
853 
854  ATH_MSG_DEBUG("Building module "<<mod->name());
855 
856  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
857  ATH_MSG_DEBUG("iPhi "<<iPhi);
858  for(int is=0; is<2; is++) { // module side
859 
860  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
861  if (element) {
862  // get element location for debugging
863  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
864  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
865  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
866 
867  // add element to the AlignModule
868  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
869  // and fill the corresponding map
870  (*sctIdHashMap)[element->identifyHash()] = mod;
871  }
872  }
873  }
874 
875  // for endcap we move the CoG
876 
877 
878 
879  Amg::Translation3D translation(mod->centerOfGravity());
880  ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
881  Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
882  mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
883 
884  m_alignModuleListPtr->push_back(mod);
885  }
886  }
887  }
888  }
889 
890  //_______________________________________________________________________
891  void SCTGeometryManagerTool::buildL2Old()
892  {
893  ATH_MSG_INFO("Preparing the old SCT Barrel geometry for L2: 4 barrel layers + 2 x 9 endcap disks");
894 
895  const Amg::Transform3D transform = Amg::Transform3D::Identity();
896 
897  unsigned int maxHash = m_idHelper->wafer_hash_max();
898  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
899 
900  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
901  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
902  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
903 
904  // first endcaap C
905  int iEndcapIndex = 1;
906  int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
907  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
908  if (!m_detManager->numerology().useDisk(iWheel))
909  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
910 
911  Trk::AlignModule * mod = new Trk::AlignModule(this);
912  mod->setIdHash(getNextIDHash());
913  mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
914 
915  std::stringstream name;
916  name <<"SCT/EndcapC/Disk_" << iWheel;
917  mod->setName(name.str());
918 
919  if(!m_alignEndcaps || !moduleSelected(mod)) {
920  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
921  delete mod;
922  continue;
923  }
924 
925  ATH_MSG_DEBUG("Building module "<<mod->name());
926 
927  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
928  ATH_MSG_DEBUG("iEta "<<iEta);
929  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
930  ATH_MSG_DEBUG("iPhi "<<iPhi);
931  for(int is=0; is<2; is++) { // module side
932  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
933  if (element) {
934  // get element location for debugging
935  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
936  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
937  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
938 
939  // add element to the AlignModule
940  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
941  // and fill the corresponding map
942  (*sctIdHashMap)[element->identifyHash()] = mod;
943  }
944  }
945  }
946  }
947 
948  // for endcap we move the CoG
949  Amg::Translation3D translation(mod->centerOfGravity());
950  ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
951  Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
952  mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
953 
954  m_alignModuleListPtr->push_back(mod);
955  }
956 
957  // now barrel
958  for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
959  if (!m_detManager->numerology().useLayer(iLayer))
960  ATH_MSG_INFO(" Layer "<<iLayer<<" not present");
961 
962  // create the AlignModule
963  Trk::AlignModule * mod = new Trk::AlignModule(this);
964  mod->setIdHash(getNextIDHash());
965  mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, 0, 0, 0));
966 
967  std::stringstream name;
968  name <<"SCT/Barrel/Layer_"<<iLayer;
969  mod->setName(name.str());
970 
971  if(!m_alignBarrel || !moduleSelected(mod)) {
972  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
973  delete mod;
974  continue;
975  }
976 
977  ATH_MSG_DEBUG("Building module "<<mod->name());
978 
979  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
980  ATH_MSG_DEBUG("iPhi "<<iPhi);
981  for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
982  ATH_MSG_DEBUG("iEta "<<iEta);
983  if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
984  // iEta=0 is not defined for the SCT
985  continue;
986  for (int is = 0; is < 2; is++) { // module side
987 
988  const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
989  if (element) {
990  // get element location for debugging
991  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
992  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
993  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
994 
995  // add element to the AlignModule
996  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
997  // and fill the corresponding map
998  (*sctIdHashMap)[element->identifyHash()] = mod;
999  }
1000  }
1001  }
1002  }
1003 
1004  m_alignModuleListPtr->push_back(mod);
1005  }
1006 
1007  // and endcaap A
1008  iEndcapIndex = 0;
1009  iSide = m_detManager->numerology().endcapId(iEndcapIndex);
1010  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
1011  if (!m_detManager->numerology().useDisk(iWheel))
1012  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
1013 
1014  Trk::AlignModule * mod = new Trk::AlignModule(this);
1015  mod->setIdHash(getNextIDHash());
1016  mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,0,0,0));
1017 
1018  std::stringstream name;
1019  name <<"SCT/EndcapA/Disk_" << iWheel;
1020  mod->setName(name.str());
1021 
1022  if(!m_alignEndcaps || !moduleSelected(mod)) {
1023  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
1024  delete mod;
1025  continue;
1026  }
1027 
1028  ATH_MSG_DEBUG("Building module "<<mod->name());
1029 
1030  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
1031  ATH_MSG_DEBUG("iEta "<<iEta);
1032  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
1033  ATH_MSG_DEBUG("iPhi "<<iPhi);
1034  for(int is=0; is<2; is++) { // module side
1035  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
1036  if (element) {
1037  // get element location for debugging
1038  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
1039  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
1040  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
1041 
1042  // add element to the AlignModule
1043  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
1044  // and fill the corresponding map
1045  (*sctIdHashMap)[element->identifyHash()] = mod;
1046  }
1047  }
1048  }
1049  }
1050 
1051  // for endcap we move the CoG
1052  Amg::Translation3D translation(mod->centerOfGravity());
1053  ATH_MSG_DEBUG("Endcap:" << iSide << " , disk: " << iWheel << " , CoG (" << translation.x() << " , " << translation.y() << " , " << translation.z() << " ) " );
1054 
1055  Amg::Transform3D localtoglobal = translation * Amg::RotationMatrix3D::Identity();
1056 
1057  mod->setGlobalFrameToAlignFrameTransform(localtoglobal.inverse());
1058 
1059  m_alignModuleListPtr->push_back(mod);
1060  }
1061  }
1062 
1063  //_______________________________________________________________________
1064  void SCTGeometryManagerTool::buildL3Barrel()
1065  {
1066  // ========================================
1067  // BARREL
1068  ATH_MSG_INFO("Preparing the SCT Barrel geometry for L3: module (bi-wafer) level");
1069 
1070  // =========================================
1071  // get all modules for SCT Barrel Level 3 alignment
1072  // 2112 modules
1073 
1074  unsigned int maxHash = m_idHelper->wafer_hash_max();
1075  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
1076 
1077  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
1078  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
1079  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
1080 
1081  for (int iLayer = 0; iLayer < m_detManager->numerology().numLayers(); iLayer++) {
1082  if (!m_detManager->numerology().useLayer(iLayer))
1083  ATH_MSG_INFO(" Layer "<<iLayer<<" not present");
1084  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForLayer(iLayer); iPhi++) {
1085  ATH_MSG_DEBUG("iPhi "<<iPhi);
1086  for (int iEta = m_detManager->numerology().beginEtaModuleForLayer(iLayer); iEta < m_detManager->numerology().endEtaModuleForLayer(iLayer); iEta++) {
1087  ATH_MSG_DEBUG("iEta "<<iEta);
1088  if (!iEta && m_detManager->numerology().skipEtaZeroForLayer(iLayer))
1089  // iEta=0 is not defined for the SCT
1090  continue;
1091 
1092  // create the AlignModule
1093  Trk::AlignModule * mod = new Trk::AlignModule(this);
1094  mod->setIdHash(getNextIDHash());
1095  mod->setIdentifier(m_idHelper->wafer_id(0, iLayer, iPhi, iEta, 0));
1096 
1097  std::stringstream name;
1098  name <<"SCT/Barrel/Layer_"<<iLayer<<"/Phi_"<<iPhi<<"/Eta_"<<iEta;
1099  mod->setName(name.str());
1100 
1101  if(!m_alignBarrel || !moduleSelected(mod)) {
1102  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
1103  delete mod;
1104  continue;
1105  }
1106 
1107  ATH_MSG_DEBUG("Building module "<<mod->name());
1108 
1109  for(int is=0;is<2;is++) { // module side
1110 
1111  const SiDetectorElement * element = m_detManager->getDetectorElement(0, iLayer, iPhi, iEta, is);
1112 
1113  if (element) {
1114  // get element location for debugging
1115  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
1116  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
1117  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
1118 
1119  // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1120  // should this be identity? if yes, it means that the
1121  // db frame constants are in the global frame. Or not?
1122  const Amg::Transform3D transform = Amg::Transform3D::Identity();
1123  // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1124 
1125  // add element to the AlignModule
1126  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
1127  // and fill the corresponding map
1128  (*sctIdHashMap)[element->identifyHash()] = mod;
1129 
1130  // the local frame for the SCT barrel modules is always that
1131  // of the Rphi element, NOT stereo !
1132  // in any case the correct frame can be obtained using moduleTransform() call
1133  // regardless of the side, but we only need it once so we retrieve it
1134  // for the non-stereo side
1135  if(!element->isStereo())
1136  mod->setGlobalFrameToAlignFrameTransform(element->moduleTransform().inverse());
1137  }
1138  }
1139 
1140  m_alignModuleListPtr->push_back(mod);
1141  }
1142  }
1143 
1144  }
1145  }
1146 
1147  //_______________________________________________________________________
1148  void SCTGeometryManagerTool::buildL3Endcaps()
1149  {
1150  // ========================================
1151  // ENDCAPs
1152  ATH_MSG_INFO("Preparing the SCT Endcap geometry for L3: module (bi-wafer) level");
1153 
1154  // =========================================
1155  // get all modules for SCT Endcap Level 3 alignment
1156  // 2 * 988
1157 
1158  unsigned int maxHash = m_idHelper->wafer_hash_max();
1159  ATH_MSG_DEBUG("maxHash for the SCT "<<maxHash);
1160 
1161  if(!m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT))
1162  m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT) = new Trk::AlignModuleList((size_t)(maxHash),nullptr);
1163  Trk::AlignModuleList * sctIdHashMap = m_idHashToAlignModuleMapsPtr->at(Trk::AlignModule::SCT);
1164 
1165  for (int iEndcapIndex = 0; iEndcapIndex < m_detManager->numerology().numEndcaps(); iEndcapIndex++) {
1166  int iSide = m_detManager->numerology().endcapId(iEndcapIndex);
1167  for (int iWheel = 0; iWheel < m_detManager->numerology().numDisks(); iWheel++) {
1168  if (!m_detManager->numerology().useDisk(iWheel))
1169  ATH_MSG_INFO(" Disk "<<iWheel<<" not present");
1170  for (int iEta = 0; iEta < m_detManager->numerology().numRingsForDisk(iWheel); iEta++) {
1171  ATH_MSG_DEBUG("iEta "<<iEta);
1172  for (int iPhi = 0; iPhi < m_detManager->numerology().numPhiModulesForDiskRing(iWheel,iEta); iPhi++) {
1173  ATH_MSG_DEBUG("iPhi "<<iPhi);
1174 
1175  Trk::AlignModule * mod = new Trk::AlignModule(this);
1176  mod->setIdHash(getNextIDHash());
1177  mod->setIdentifier(m_idHelper->wafer_id(iSide,iWheel,iPhi,iEta,0));
1178 
1179  std::stringstream name;
1180  if(iEndcapIndex == 0)
1181  name<<"SCT/EndcapA";
1182  else
1183  name<<"SCT/EndcapC";
1184  name<<"/Disk_"<<iWheel<<"/Phi_"<<iPhi<<"/Eta_"<<iEta;
1185  mod->setName(name.str());
1186 
1187  if(!m_alignEndcaps || !moduleSelected(mod)) {
1188  ATH_MSG_DEBUG("Module "<<mod->name()<<" NOT selected");
1189  delete mod;
1190  continue;
1191  }
1192 
1193  ATH_MSG_DEBUG("Building module "<<mod->name());
1194 
1195  for(int is=0; is<2; is++) { // module side
1196 
1197  const SiDetectorElement * element = m_detManager->getDetectorElement(iSide, iWheel, iPhi, iEta, is);
1198  if (element) {
1199  // get element location for debugging
1200  // HepGeom::Point3D<double> center = element->transform() * HepGeom::Point3D<double>();
1201  // ATH_MSG_DEBUG(" SCTDetectorElement id: " << id << " with center = " << center);
1202  // ATH_MSG_DEBUG(" Is Barrel: "<< m_idHelper->is_barrel(id));
1203 
1204  // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1205  // should this be identity? if yes, it means that the
1206  // db frame constants are in the global frame. Or not?
1207  const Amg::Transform3D transform = Amg::Transform3D::Identity();
1208  // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1209 
1210  // add element to the AlignModule
1211  mod->addDetElement(Trk::AlignModule::SCT,element,transform);
1212  // and fill the corresponding map
1213  (*sctIdHashMap)[element->identifyHash()] = mod;
1214 
1215  // the local frame for the SCT endcap modules is always that
1216  // of the Rphi element, NOT stereo !
1217  // in any case the correct frame can be obtained using moduleTransform() call
1218  // regardless of the side, but we only need it once so we retrieve it
1219  // for the non-stereo side
1220  if(!element->isStereo())
1221  mod->setGlobalFrameToAlignFrameTransform(element->moduleTransform().inverse());
1222  }
1223  }
1224 
1225  m_alignModuleListPtr->push_back(mod);
1226  }
1227  }
1228  }
1229  }
1230 
1231  }
1232 
1233  //________________________________________________________________________
1234  void SCTGeometryManagerTool::addModuleParameters(Trk::AlignModule * module, DataVector< DataVector<Trk::AlignPar> > * allFullModPars, DataVector< DataVector<Trk::AlignPar> > * allActiveModPars)
1235  {
1236  // prepare all parameters
1237  DataVector<Trk::AlignPar> * fullModPars = new DataVector<Trk::AlignPar>(SG::OWN_ELEMENTS);
1238  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransX));
1239  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransY));
1240  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::TransZ));
1241  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotX));
1242  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotY));
1243  fullModPars->push_back(new Trk::AlignPar(module,Trk::AlignModule::RotZ));
1244 
1245  // set sigmas
1246  setSigmas(module,fullModPars);
1247 
1248  // select active parameters based on jobOption properties
1249  DataVector<Trk::AlignPar> * activeModPars = new DataVector<Trk::AlignPar>(SG::VIEW_ELEMENTS);
1250  for(unsigned int ipar=0;ipar<fullModPars->size();++ipar) {
1251 
1252  Identifier AlimodID = module->identify();
1253  Identifier modID = ((module->detElementCollection(Trk::AlignModule::SCT))->at(0))->identify();
1254 
1255  if(m_idHelper->is_barrel(modID)) {
1256  ATH_MSG_DEBUG("SCT barrel module with id "<<AlimodID);
1257  if( (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransX && m_alignBarrelX)
1258  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransY && m_alignBarrelY)
1259  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransZ && m_alignBarrelZ)
1260  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotX && m_alignBarrelRotX)
1261  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotY && m_alignBarrelRotY)
1262  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotZ && m_alignBarrelRotZ) ) {
1263  ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is now active");
1264  activeModPars->push_back(fullModPars->at(ipar));
1265  }
1266  else
1267  ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is NOT active");
1268  }
1269  else {
1270  ATH_MSG_DEBUG("SCT endcap module with id "<<AlimodID);
1271  if( (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransX && m_alignEndcapX)
1272  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransY && m_alignEndcapY)
1273  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::TransZ && m_alignEndcapZ)
1274  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotX && m_alignEndcapRotX)
1275  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotY && m_alignEndcapRotY)
1276  || (fullModPars->at(ipar)->paramType() == Trk::AlignModule::RotZ && m_alignEndcapRotZ) ) {
1277  ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is now active");
1278  activeModPars->push_back(fullModPars->at(ipar));
1279  }
1280  else
1281  ATH_MSG_DEBUG("parameter type "<<fullModPars->at(ipar)->paramType()<<" is NOT active");
1282  }
1283  }
1284 
1285  // now add parameters to the list
1286  allFullModPars->push_back(fullModPars);
1287  allActiveModPars->push_back(activeModPars);
1288  }
1289 
1290  //________________________________________________________________________
1291  void SCTGeometryManagerTool::setSigmas(Trk::AlignModule * module, DataVector<Trk::AlignPar> * modPars)
1292  {
1293  for(unsigned int ipar=0;ipar<modPars->size();++ipar) {
1294  // barrel
1295  if(m_idHelper->is_barrel(module->identify()))
1296  switch(modPars->at(ipar)->paramType()) {
1297  case Trk::AlignModule::TransX:
1298  modPars->at(ipar)->setSigma(m_sigmaBarrelX);
1299  modPars->at(ipar)->setSoftCut(m_softcutBarrelX);
1300  break;
1301  case Trk::AlignModule::TransY:
1302  modPars->at(ipar)->setSigma(m_sigmaBarrelY);
1303  modPars->at(ipar)->setSoftCut(m_softcutBarrelY);
1304  break;
1305  case Trk::AlignModule::TransZ:
1306  modPars->at(ipar)->setSigma(m_sigmaBarrelZ);
1307  modPars->at(ipar)->setSoftCut(m_softcutBarrelZ);
1308  break;
1309  case Trk::AlignModule::RotX:
1310  modPars->at(ipar)->setSigma(m_sigmaBarrelRotX);
1311  modPars->at(ipar)->setSoftCut(m_softcutBarrelRotX);
1312  break;
1313  case Trk::AlignModule::RotY:
1314  modPars->at(ipar)->setSigma(m_sigmaBarrelRotY);
1315  modPars->at(ipar)->setSoftCut(m_softcutBarrelRotY);
1316  break;
1317  case Trk::AlignModule::RotZ:
1318  modPars->at(ipar)->setSigma(m_sigmaBarrelRotZ);
1319  modPars->at(ipar)->setSoftCut(m_softcutBarrelRotZ);
1320  break;
1321  default:
1322  break;
1323  }
1324  // end-caps
1325  else
1326  switch(modPars->at(ipar)->paramType()) {
1327  case Trk::AlignModule::TransX:
1328  modPars->at(ipar)->setSigma(m_sigmaEndcapX);
1329  modPars->at(ipar)->setSoftCut(m_softcutEndcapX);
1330  break;
1331  case Trk::AlignModule::TransY:
1332  modPars->at(ipar)->setSigma(m_sigmaEndcapY);
1333  modPars->at(ipar)->setSoftCut(m_softcutEndcapY);
1334  break;
1335  case Trk::AlignModule::TransZ:
1336  modPars->at(ipar)->setSigma(m_sigmaEndcapZ);
1337  modPars->at(ipar)->setSoftCut(m_softcutEndcapZ);
1338  break;
1339  case Trk::AlignModule::RotX:
1340  modPars->at(ipar)->setSigma(m_sigmaEndcapRotX);
1341  modPars->at(ipar)->setSoftCut(m_softcutEndcapRotX);
1342  break;
1343  case Trk::AlignModule::RotY:
1344  modPars->at(ipar)->setSigma(m_sigmaEndcapRotY);
1345  modPars->at(ipar)->setSoftCut(m_softcutEndcapRotY);
1346  break;
1347  case Trk::AlignModule::RotZ:
1348  modPars->at(ipar)->setSigma(m_sigmaEndcapRotZ);
1349  modPars->at(ipar)->setSoftCut(m_softcutEndcapRotZ);
1350  break;
1351  default:
1352  break;
1353  }
1354  }
1355  }
1356 
1357  //________________________________________________________________________
1358  bool SCTGeometryManagerTool::moduleSelected(unsigned long long id)
1359  {
1360  if(!m_doModuleSelection)
1361  return true;
1362 
1363  int nsel = m_moduleSelection.size();
1364  for(int i=0;i<nsel;++i)
1365  if(m_moduleSelection.at(i) == id)
1366  return true;
1367 
1368  return false;
1369  }
1370 
1371  //________________________________________________________________________
1372  bool SCTGeometryManagerTool::moduleSelected(Trk::AlignModule * mod)
1373  {
1374  return moduleSelected(mod->identify().get_compact());
1375  }
1376 
1377  //________________________________________________________________________
1378  void SCTGeometryManagerTool::dumpGeometry()
1379  {
1380  ATH_MSG_INFO("---------------------------------------------------");
1381  ATH_MSG_INFO("Summary of the alignment geometry");
1382  ATH_MSG_INFO("Number of alignable objects: "<< m_alignModuleList.size());
1383  for(unsigned int i=0;i<m_alignModuleList.size();i++) {
1384  const Trk::AlignModule* module = m_alignModuleList.at(i);
1385  ATH_MSG_INFO(i<<". "<< module->name());
1386  ATH_MSG_INFO(" - identifier: "<<module->identify());
1387  ATH_MSG_INFO(" - has "<<module->detElementCollection(Trk::AlignModule::SCT)->size()<<" SCT modules");
1388 
1389  Amg::Transform3D localtoglobal = (module->globalFrameToAlignFrame()).inverse();
1390  ATH_MSG_DEBUG(" - local to global : "<<std::setprecision(12)<<localtoglobal.translation()<<" "<<localtoglobal.rotation());
1391 
1392  DataVector<Trk::AlignPar> * pars = m_alignModuleTool->getAlignPars(module);
1393  int npars = pars->size();
1394  ATH_MSG_DEBUG(" - number of active transform parameters: "<<npars);
1395  for(int j=0;j<npars;j++)
1396  ATH_MSG_DEBUG(" * par "<<j<<" \'"<<(*pars)[j]->dumpType()<<"\' : sigma = "<<(*pars)[j]->sigma()<<" , softCut = "<<(*pars)[j]->softCut());
1397  }
1398  ATH_MSG_INFO("---------------------------------------------------");
1399  }
1400 
1401 } // end namespace
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
make_hlt_rep.pars
pars
Definition: make_hlt_rep.py:90
InDet::SCTGeometryManagerTool::buildL2Endcaps
void buildL2Endcaps()
creates L2 AlignModules for SCT endcaps
Definition: SCTGeometryManagerTool.cxx:727
InDet::SCTGeometryManagerTool::m_softcutEndcapY
double m_softcutEndcapY
Definition: SCTGeometryManagerTool.h:181
Trk::IGeometryManagerTool::m_alignLevelBarrel
int m_alignLevelBarrel
alignment level
Definition: IGeometryManagerTool.h:116
InDet::SCTGeometryManagerTool::buildL3Barrel
void buildL3Barrel()
creates L3 AlignModules for SCT barrel
Definition: SCTGeometryManagerTool.cxx:1064
xAOD::identify
const Identifier & identify(const UncalibratedMeasurement *meas)
Returns the associated identifier from the muon measurement.
Definition: MuonSpectrometer/MuonPhaseII/Event/xAOD/xAODMuonPrepData/Root/UtilFunctions.cxx:82
InDet::SCTGeometryManagerTool::m_idHelper
const SCT_ID * m_idHelper
pointer to SCT detector manager
Definition: SCTGeometryManagerTool.h:127
InDet::SCTGeometryManagerTool::moduleSelected
bool moduleSelected(Trk::AlignModule *mod)
check wheather module is selected for module pointer
Definition: SCTGeometryManagerTool.cxx:1372
SCT_ID.h
This is an Identifier helper class for the SCT subdetector. This class is a factory for creating comp...
InDet::SCTGeometryManagerTool::addModuleParameters
void addModuleParameters(Trk::AlignModule *module, DataVector< DataVector< Trk::AlignPar > > *allFullModPars, DataVector< DataVector< Trk::AlignPar > > *allActiveModPars)
adds alignment parameters for the module checks for active parameters and calls setSigmas()
Definition: SCTGeometryManagerTool.cxx:1234
pdg_comparison.sigma
sigma
Definition: pdg_comparison.py:324
InDetDD::SiDetectorManager::numerology
const SiNumerology & numerology() const
Access Numerology.
Definition: SiDetectorManager.h:128
InDet::SCTGeometryManagerTool::m_alignBarrel
bool m_alignBarrel
Definition: SCTGeometryManagerTool.h:142
Trk::AlignModule::TransZ
@ TransZ
Definition: AlignModule.h:54
python.Constants.FATAL
int FATAL
Definition: Control/AthenaCommon/python/Constants.py:18
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
InDet::SCTGeometryManagerTool::buildL25Endcaps
void buildL25Endcaps()
creates L25 AlignModules for SCT endcaps
Definition: SCTGeometryManagerTool.cxx:805
SCTGeometryManagerTool.h
Trk::AlignModuleList
std::vector< AlignModule * > AlignModuleList
Definition: AlignModuleList.h:37
SG::VIEW_ELEMENTS
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
Definition: OwnershipPolicy.h:18
Trk::IGeometryManagerTool::m_hashCounter
int m_hashCounter
variable for setting the idHash of the AlignModules
Definition: IGeometryManagerTool.h:111
index
Definition: index.py:1
InDet::SCTGeometryManagerTool::m_moduleSelection
std::vector< unsigned long long > m_moduleSelection
Definition: SCTGeometryManagerTool.h:140
InDet
Primary Vertex Finder.
Definition: VP1ErrorUtils.h:36
Trk::AlignModule::addDetElement
void addDetElement(AlignModule::DetectorType detType, const TrkDetElementBase *det, const Amg::Transform3D &transform, Identifier id=Identifier())
used to add a detector element to the align module with a align frame to detector element local frame...
Definition: AlignModule.cxx:126
InDet::SCTGeometryManagerTool::m_sigmaEndcapZ
double m_sigmaEndcapZ
Definition: SCTGeometryManagerTool.h:168
InDetDD::SiNumerology::useLayer
bool useLayer(int layer) const
Check if layer exists.
InDet::SCTGeometryManagerTool::m_sigmaEndcapY
double m_sigmaEndcapY
Definition: SCTGeometryManagerTool.h:167
Trk::AlignModule
Definition: AlignModule.h:45
InDet::SCTGeometryManagerTool::buildL27Barrel
void buildL27Barrel()
creates L27 AlignModules for SCT barrel
Definition: SCTGeometryManagerTool.cxx:626
AthCommonMsg< AlgTool >::msgLvl
bool msgLvl(const MSG::Level lvl) const
Definition: AthCommonMsg.h:30
Trk::IGeometryManagerTool::m_alignLevelEndcaps
int m_alignLevelEndcaps
alignment level
Definition: IGeometryManagerTool.h:117
InDet::SCTGeometryManagerTool::m_alignBarrelRotY
bool m_alignBarrelRotY
Definition: SCTGeometryManagerTool.h:149
InDet::SCTGeometryManagerTool::m_sigmaEndcapRotX
double m_sigmaEndcapRotX
Definition: SCTGeometryManagerTool.h:169
InDet::SCTGeometryManagerTool::m_alignEndcaps
bool m_alignEndcaps
Definition: SCTGeometryManagerTool.h:143
Trk::AlignPar::paramType
AlignModule::TransformParameters paramType() const
returns the type of parameter (i.e.
Definition: AlignPar.h:47
InDetDD::SCT_DetectorManager::getDetectorElement
virtual const SiDetectorElement * getDetectorElement(const Identifier &id) const override
access to individual elements via Identifier
Definition: SCT_DetectorManager.cxx:66
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
InDetDD::SolidStateDetectorElementBase::identifyHash
virtual IdentifierHash identifyHash() const override final
identifier hash (inline)
InDetDD::SiNumerology::skipEtaZeroForLayer
bool skipEtaZeroForLayer(int layer) const
Check if eta_module=0 exists.
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
InDet::SCTGeometryManagerTool::m_sigmaBarrelZ
double m_sigmaBarrelZ
Definition: SCTGeometryManagerTool.h:161
InDet::SCTGeometryManagerTool::m_softcutBarrelRotZ
double m_softcutBarrelRotZ
Definition: SCTGeometryManagerTool.h:178
Trk::IGeometryManagerTool::m_alignLevel
int m_alignLevel
alignment level
Definition: IGeometryManagerTool.h:115
InDet::SCTGeometryManagerTool::m_alignBarrelZ
bool m_alignBarrelZ
Definition: SCTGeometryManagerTool.h:147
InDet::SCTGeometryManagerTool::m_alignEndcapRotX
bool m_alignEndcapRotX
Definition: SCTGeometryManagerTool.h:155
InDet::SCTGeometryManagerTool::m_sigmaEndcapRotY
double m_sigmaEndcapRotY
Definition: SCTGeometryManagerTool.h:170
AlignModule
AlignModule is a grouping of TrkDetElementBase objects, grouped according to the type of alignment,...
python.PyAthena.module
module
Definition: PyAthena.py:131
InDet::SCTGeometryManagerTool::m_sigmaEndcapRotZ
double m_sigmaEndcapRotZ
Definition: SCTGeometryManagerTool.h:171
InDet::SCTGeometryManagerTool::buildL1Barrel
void buildL1Barrel()
creates L1 AlignModules for SCT barrel
Definition: SCTGeometryManagerTool.cxx:401
InDet::SCTGeometryManagerTool::buildL0
void buildL0()
creates L0 AlignModules for SCT
Definition: SCTGeometryManagerTool.cxx:336
InDet::SCTGeometryManagerTool::m_softcutBarrelX
double m_softcutBarrelX
Definition: SCTGeometryManagerTool.h:173
InDet::SCTGeometryManagerTool::~SCTGeometryManagerTool
virtual ~SCTGeometryManagerTool()
Definition: SCTGeometryManagerTool.cxx:106
Trk::AlignPar::setSigma
void setSigma(double sigma)
sets sigma
Definition: AlignPar.h:87
Trk::AlignModule::setName
void setName(const std::string &name)
Set and return name of align module (i.e.
Definition: AlignModule.h:88
Trk::IGeometryManagerTool::m_alignModuleListPtr
Trk::AlignModuleList * m_alignModuleListPtr
pointer to module list to which the modules are added
Definition: IGeometryManagerTool.h:106
maskDeadModules.mod
mod
Definition: maskDeadModules.py:36
Trk::AlignModule::setIdentifier
void setIdentifier(Identifier identifier)
Set and return identifier of module.
Definition: AlignModule.h:96
InDet::SCTGeometryManagerTool::ReadGeometry
int ReadGeometry(int solveLevel)
read the geometry Method is called from the main AlignAlg to build the geometry based on the requeste...
Definition: SCTGeometryManagerTool.cxx:243
SCT_ID::wafer_hash_max
size_type wafer_hash_max() const
Definition: SCT_ID.cxx:646
InDet::SCTGeometryManagerTool::buildL2Old
void buildL2Old()
creates L2 AlignModules for the full Pixel in the ordering equal to the old code
Definition: SCTGeometryManagerTool.cxx:891
InDet::SCTGeometryManagerTool::m_sigmaBarrelY
double m_sigmaBarrelY
Definition: SCTGeometryManagerTool.h:160
InDetDD::SiDetectorElement::moduleTransform
const Amg::Transform3D & moduleTransform() const
Module to global frame transform.
Definition: SiDetectorElement.cxx:173
lumiFormat.i
int i
Definition: lumiFormat.py:85
SG::OWN_ELEMENTS
@ OWN_ELEMENTS
this data object owns its elements
Definition: OwnershipPolicy.h:17
Trk::IGeometryManagerTool::getNextIDHash
virtual int getNextIDHash()
get next free IDHash usable for new AlignModule
Definition: IGeometryManagerTool.h:65
Trk::AlignModule::identifyHash
IdentifierHash identifyHash() const
Set and return index of module, used by alignment classes to keep track of order of align module.
Definition: AlignModule.h:92
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
InDet::SCTGeometryManagerTool::m_idHashToAlignModuleMaps
std::vector< Trk::AlignModuleList * > m_idHashToAlignModuleMaps
Definition: SCTGeometryManagerTool.h:134
InDet::SCTGeometryManagerTool::m_softcutEndcapRotX
double m_softcutEndcapRotX
Definition: SCTGeometryManagerTool.h:183
InDet::SCTGeometryManagerTool::checkAlignLevelBarrel
bool checkAlignLevelBarrel()
check whether the Barrel alignment level is correct
Definition: SCTGeometryManagerTool.cxx:209
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
xAOD::rotation
rotation
Definition: TrackSurface_v1.cxx:15
InDet::SCTGeometryManagerTool::m_sigmaBarrelRotZ
double m_sigmaBarrelRotZ
Definition: SCTGeometryManagerTool.h:164
InDet::SCTGeometryManagerTool::buildL1Endcaps
void buildL1Endcaps()
creates L1 AlignModules for SCT endcaps
Definition: SCTGeometryManagerTool.cxx:476
Amg::Transform3D
Eigen::Affine3d Transform3D
Definition: GeoPrimitives.h:46
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
Trk::AlignPar
Definition: AlignPar.h:25
test_pyathena.parent
parent
Definition: test_pyathena.py:15
InDetDD::SiNumerology::endcapId
int endcapId(int index) const
Barrel/endcap identifier for each endcap.
Trk::AlignModule::RotX
@ RotX
Definition: AlignModule.h:54
InDet::SCTGeometryManagerTool::m_alignEndcapRotY
bool m_alignEndcapRotY
Definition: SCTGeometryManagerTool.h:156
InDetDD::SiNumerology::useDisk
bool useDisk(int disk) const
Check if disk exists.
Trk::AlignModule::SCT
@ SCT
Definition: AlignModule.h:57
InDet::SCTGeometryManagerTool::m_softcutEndcapRotY
double m_softcutEndcapRotY
Definition: SCTGeometryManagerTool.h:184
InDet::SCTGeometryManagerTool::m_softcutEndcapRotZ
double m_softcutEndcapRotZ
Definition: SCTGeometryManagerTool.h:185
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
InDet::SCTGeometryManagerTool::m_dumpGeometry
bool m_dumpGeometry
Definition: SCTGeometryManagerTool.h:188
AlignModuleList.h
InDet::SCTGeometryManagerTool::checkAlignLevel
bool checkAlignLevel()
check whether the alignment level is correct
Definition: SCTGeometryManagerTool.cxx:178
Trk::AlignModule::name
const std::string & name() const
Definition: AlignModule.h:89
InDet::SCTGeometryManagerTool::m_alignEndcapX
bool m_alignEndcapX
Definition: SCTGeometryManagerTool.h:152
InDet::SCTGeometryManagerTool::m_stripDetManagerName
Gaudi::Property< std::string > m_stripDetManagerName
Definition: SCTGeometryManagerTool.h:190
InDet::SCTGeometryManagerTool::m_alignBarrelX
bool m_alignBarrelX
Definition: SCTGeometryManagerTool.h:145
Trk
Ensure that the ATLAS eigen extensions are properly loaded.
Definition: FakeTrackBuilder.h:9
InDet::SCTGeometryManagerTool::m_detManager
const InDetDD::SCT_DetectorManager * m_detManager
pointer to SCT detector manager
Definition: SCTGeometryManagerTool.h:125
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
AlignModule.h
InDet::SCTGeometryManagerTool::m_alignEndcapZ
bool m_alignEndcapZ
Definition: SCTGeometryManagerTool.h:154
InDetDD::SiDetectorElement
Definition: SiDetectorElement.h:109
InDet::SCTGeometryManagerTool::m_softcutBarrelRotX
double m_softcutBarrelRotX
Definition: SCTGeometryManagerTool.h:176
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
HI::TowerBins::numLayers
constexpr unsigned int numLayers()
Definition: HIEventDefs.h:23
InDet::SCTGeometryManagerTool::m_sigmaBarrelX
double m_sigmaBarrelX
Definition: SCTGeometryManagerTool.h:159
Trk::IGeometryManagerTool::m_logStream
std::ostream * m_logStream
logfile output stream
Definition: IGeometryManagerTool.h:113
InDet::SCTGeometryManagerTool::m_useOldL2
bool m_useOldL2
Definition: SCTGeometryManagerTool.h:187
SiDetectorElement.h
Trk::iPhi
@ iPhi
Definition: ParamDefs.h:47
InDet::SCTGeometryManagerTool::buildGeometry
void buildGeometry()
builds geometry for SCT alignment
Definition: SCTGeometryManagerTool.cxx:286
Trk::AlignModule::setIdHash
void setIdHash(IdentifierHash id)
Definition: AlignModule.h:93
Trk::AlignModule::TransX
@ TransX
Definition: AlignModule.h:54
InDet::SCTGeometryManagerTool::m_alignEndcapRotZ
bool m_alignEndcapRotZ
Definition: SCTGeometryManagerTool.h:157
DataVector.h
An STL vector of pointers that by default owns its pointed-to elements.
DeMoScan.index
string index
Definition: DeMoScan.py:362
InDet::SCTGeometryManagerTool::m_sigmaBarrelRotY
double m_sigmaBarrelRotY
Definition: SCTGeometryManagerTool.h:163
InDet::SCTGeometryManagerTool::m_softcutBarrelRotY
double m_softcutBarrelRotY
Definition: SCTGeometryManagerTool.h:177
InDet::SCTGeometryManagerTool::m_sigmaBarrelRotX
double m_sigmaBarrelRotX
Definition: SCTGeometryManagerTool.h:162
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
Trk::IGeometryManagerTool::m_idHashToAlignModuleMapsPtr
std::vector< Trk::AlignModuleList * > * m_idHashToAlignModuleMapsPtr
pointer to vector of hashMaps to which the elements are added
Definition: IGeometryManagerTool.h:109
InDet::SCTGeometryManagerTool::m_alignModuleList
Trk::AlignModuleList m_alignModuleList
Definition: SCTGeometryManagerTool.h:131
InDet::SCTGeometryManagerTool::m_softcutEndcapZ
double m_softcutEndcapZ
Definition: SCTGeometryManagerTool.h:182
python.Constants.INFO
int INFO
Definition: Control/AthenaCommon/python/Constants.py:15
InDetDD
Message Stream Member.
Definition: FakeTrackBuilder.h:8
Amg::RotationMatrix3D
Eigen::Matrix< double, 3, 3 > RotationMatrix3D
Definition: GeoPrimitives.h:49
Trk::AlignModule::RotY
@ RotY
Definition: AlignModule.h:54
InDet::SCTGeometryManagerTool::m_alignEndcapY
bool m_alignEndcapY
Definition: SCTGeometryManagerTool.h:153
InDet::SCTGeometryManagerTool::m_fullAlignParList
DataVector< DataVector< Trk::AlignPar > > * m_fullAlignParList
Definition: SCTGeometryManagerTool.h:137
InDetDD::SiNumerology::beginEtaModuleForLayer
int beginEtaModuleForLayer(int layer) const
First eta_module number for a layer.
DEBUG
#define DEBUG
Definition: page_access.h:11
Amg::Translation3D
Eigen::Translation< double, 3 > Translation3D
Definition: GeoPrimitives.h:44
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
IAlignModuleTool.h
InDet::SCTGeometryManagerTool::m_alignBarrelRotX
bool m_alignBarrelRotX
Definition: SCTGeometryManagerTool.h:148
AlignPar.h
InDet::SCTGeometryManagerTool::setAlignLevel
void setAlignLevel(int level)
set alignment level
Definition: SCTGeometryManagerTool.cxx:171
InDet::SCTGeometryManagerTool::buildL3Endcaps
void buildL3Endcaps()
creates L3 AlignModules for SCT endcaps
Definition: SCTGeometryManagerTool.cxx:1148
Trk::AlignPar::setSoftCut
void setSoftCut(double softCut)
sets soft-mode-cut
Definition: AlignPar.h:90
InDet::SCTGeometryManagerTool::m_alignBarrelRotZ
bool m_alignBarrelRotZ
Definition: SCTGeometryManagerTool.h:150
InDet::SCTGeometryManagerTool::m_sigmaEndcapX
double m_sigmaEndcapX
Definition: SCTGeometryManagerTool.h:166
InDet::SCTGeometryManagerTool::buildL2Barrel
void buildL2Barrel()
creates L2 AlignModules for SCT barrel
Definition: SCTGeometryManagerTool.cxx:556
InDet::SCTGeometryManagerTool::m_softcutBarrelY
double m_softcutBarrelY
Definition: SCTGeometryManagerTool.h:174
DataVector::at
const T * at(size_type n) const
Access an element, as an rvalue.
InDet::SCTGeometryManagerTool::m_softcutBarrelZ
double m_softcutBarrelZ
Definition: SCTGeometryManagerTool.h:175
InDet::SCTGeometryManagerTool::finalize
StatusCode finalize()
Definition: SCTGeometryManagerTool.cxx:163
InDet::SCTGeometryManagerTool::m_doModuleSelection
bool m_doModuleSelection
Definition: SCTGeometryManagerTool.h:139
InDet::SCTGeometryManagerTool::initialize
StatusCode initialize()
Definition: SCTGeometryManagerTool.cxx:120
AthAlgTool
Definition: AthAlgTool.h:26
SCT_ID::wafer_id
Identifier wafer_id(int barrel_ec, int layer_disk, int phi_module, int eta_module, int side) const
For a single side of module.
Definition: SCT_ID.h:464
SCT_DetectorManager.h
InDetDD::SiDetectorElement::isStereo
bool isStereo() const
Check if it is the stereo side (useful for SCT)
Definition: SiDetectorElement.cxx:300
Trk::AlignModule::RotZ
@ RotZ
Definition: AlignModule.h:54
xAOD::iEta
setScale setgFexType iEta
Definition: gFexJetRoI_v1.cxx:77
InDet::SCTGeometryManagerTool::setSigmas
void setSigmas(Trk::AlignModule *mod, DataVector< Trk::AlignPar > *modPars)
sets sigmas for modules
Definition: SCTGeometryManagerTool.cxx:1291
Trk::AlignModule::TransY
@ TransY
Definition: AlignModule.h:54
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
InDet::SCTGeometryManagerTool::m_softcutEndcapX
double m_softcutEndcapX
Definition: SCTGeometryManagerTool.h:180
Trk::AlignModule::identify
Identifier identify() const
Definition: AlignModule.h:97
InDet::SCTGeometryManagerTool::dumpGeometry
void dumpGeometry()
print basic geometry info to screen
Definition: SCTGeometryManagerTool.cxx:1378
InDet::SCTGeometryManagerTool::m_alignParList
DataVector< DataVector< Trk::AlignPar > > * m_alignParList
Definition: SCTGeometryManagerTool.h:136
InDet::SCTGeometryManagerTool::checkAlignLevelEndcaps
bool checkAlignLevelEndcaps()
check whether the Endcap alignment level is correct
Definition: SCTGeometryManagerTool.cxx:226
InDet::SCTGeometryManagerTool::m_alignBarrelY
bool m_alignBarrelY
Definition: SCTGeometryManagerTool.h:146
SCT_ID::is_barrel
bool is_barrel(const Identifier &id) const
Test for barrel - WARNING: id MUST be sct id, otherwise answer is not accurate. Use SiliconID for gen...
Definition: SCT_ID.h:721
InDet::SCTGeometryManagerTool::m_alignModuleTool
ToolHandle< Trk::IAlignModuleTool > m_alignModuleTool
pointer to AlignModuleTool
Definition: SCTGeometryManagerTool.h:129
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Mon Sep 29 2025 21:18:24 for ATLAS Offline Software by   doxygen 1.8.18