CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
TrackerGeometryCompare Class Reference

#include <TrackerGeometryCompare.h>

Inheritance diagram for TrackerGeometryCompare:
edm::one::EDAnalyzer<> edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Types

typedef AlignTransform SurveyValue
 
typedef Alignments SurveyValues
 
- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 Read from DB and print survey info. More...
 
void endJob () override
 
 TrackerGeometryCompare (const edm::ParameterSet &)
 Do nothing. Required by framework. More...
 
- Public Member Functions inherited from edm::one::EDAnalyzer<>
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
const EDAnalyzer & operator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void addSurveyInfo (Alignable *ali)
 
void compareGeometries (Alignable *refAli, Alignable *curAli, const TrackerTopology *tTopo, const edm::EventSetup &iSetup)
 
void compareSurfaceDeformations (TTree *_inputTree11, TTree *_inputTree12)
 
void createROOTGeometry (const edm::EventSetup &iSetup)
 
void diffCommonTrackerSystem (Alignable *refAli, Alignable *curAli)
 
void fillIdentifiers (int subdetlevel, int rawid, const TrackerTopology *tTopo)
 
void fillTree (Alignable *refAli, const AlgebraicVector &diff, const TrackerTopology *tTopo, const edm::EventSetup &iSetup)
 
bool passIdCut (uint32_t)
 
void setCommonTrackerSystem ()
 
void surveyToTracker (AlignableTracker *ali, Alignments *alignVals, AlignmentErrorsExtended *alignErrors)
 

Private Attributes

TTree * alignTree_
 
float alphaVal_
 
int badModuleQuality_
 
float betaVal_
 
align::StructureType commonTrackerLevel_
 
const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord > cpvTokenDD4hep_
 
const edm::ESGetToken< DDCompactView, IdealGeometryRecord > cpvTokenDDD_
 
AlignableTracker * currentTracker
 
float dalphaVal_
 
float daVal_
 
float dbetaVal_
 
float dbVal_
 
int detDim_
 
bool detIdFlag_
 
std::string detIdFlagFile_
 
std::vector< uint32_t > detIdFlagVector_
 
float dgammaVal_
 
float dgVal_
 
float dphiVal_
 
float drVal_
 
AlignableTracker * dummyTracker
 
float duVal_
 
float dvVal_
 
float dwVal_
 
float dxVal_
 
float dyVal_
 
float dzVal_
 
float etaVal_
 
std::string filename_
 
bool firstEvent_
 
bool fromDD4hep_
 
float gammaVal_
 
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
 
int id_
 
uint32_t identifiers_ [6]
 
int inModuleList_
 
std::string inputFilename1_
 
std::string inputFilename2_
 
TFile * inputRootFile1_
 
TFile * inputRootFile2_
 
TTree * inputTree01_
 
TTree * inputTree02_
 
TTree * inputTree11_
 
TTree * inputTree12_
 
std::string inputTreenameAlign_
 
std::string inputTreenameDeform_
 
int level_
 
const std::vector< std::string > levelStrings_
 
std::map< std::string, TH1D * > m_h1_
 
int m_nBins_
 
double m_rangeHigh_
 
double m_rangeLow_
 
std::vector< align::StructureType > m_theLevels
 
std::vector< TrackerMap > m_vtkmap_
 
int mid_
 
int mlevel_
 
int moduleInList_
 
std::vector< int > moduleList_
 
std::ifstream moduleListFile_
 
std::string moduleListName_
 
SiPixelPI::phase phase_
 
float phiVal_
 
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
 
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
 
AlignableTracker * referenceTracker
 
float rVal_
 
std::string setCommonTrackerSystem_
 
const edm::ESGetToken< SiStripQuality, SiStripQualityRcd > stripQualityToken_
 
int sublevel_
 
double surfDeform_ [13]
 
std::string surfdir_
 
float surLength_
 
double surRot_ [9]
 
float surWidth_
 
TFile * theFile_
 
const SurveyErrors * theSurveyErrors
 
unsigned int theSurveyIndex
 
const Alignments * theSurveyValues
 
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
 
align::PositionType TrackerCommonCM_
 
align::GlobalVector TrackerCommonR_
 
align::GlobalVector TrackerCommonT_
 
int type_
 
int useDetId_
 
std::string weightBy_
 
bool weightById_
 
std::string weightByIdFile_
 
std::vector< unsigned int > weightByIdVector_
 
bool writeToDB_
 
float xVal_
 
float yVal_
 
float zVal_
 

Additional Inherited Members

- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Module that reads survey info from DB and prints them out.

Date
2012/12/02 22:13:12
Revision
1.14
Author
Nhan Tran

******** Including surface deformations in the geometry comparison ********

Definition at line 41 of file TrackerGeometryCompare.h.

Member Typedef Documentation

◆ SurveyValue

typedef AlignTransform TrackerGeometryCompare::SurveyValue

Definition at line 43 of file TrackerGeometryCompare.h.

◆ SurveyValues

typedef Alignments TrackerGeometryCompare::SurveyValues

Definition at line 44 of file TrackerGeometryCompare.h.

Constructor & Destructor Documentation

◆ TrackerGeometryCompare()

TrackerGeometryCompare::TrackerGeometryCompare ( const edm::ParameterSet &  cfg)

Do nothing. Required by framework.

Definition at line 56 of file TrackerGeometryCompare.cc.

References alignTree_, alphaVal_, badModuleQuality_, betaVal_, looper::cfg, dalphaVal_, daVal_, dbetaVal_, dbVal_, detDim_, detIdFlag_, detIdFlagFile_, detIdFlagVector_, dgammaVal_, dgVal_, dphiVal_, drVal_, duVal_, dvVal_, dwVal_, dxVal_, dyVal_, dzVal_, etaVal_, filename_, groupFilesInBlocks::fin, compareTotals::fs, gammaVal_, l1ctLayer2EG_cff::id, id_, identifiers_, cuy::ii, testHGCalDigi_cfg::inFile, inModuleList_, inputFilename1_, inputFilename2_, inputTreenameAlign_, inputTreenameDeform_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, level_, m_h1_, m_nBins_, m_rangeHigh_, m_rangeLow_, m_vtkmap_, TFileDirectory::make(), mid_, mlevel_, moduleListName_, phiVal_, rVal_, setCommonTrackerSystem_, AlCaHLTBitMon_QueryRunRegistry::string, sublevel_, surfDeform_, surfdir_, surLength_, surRot_, surWidth_, theFile_, type_, useDetId_, weightBy_, weightById_, weightByIdFile_, weightByIdVector_, xVal_, yVal_, and zVal_.

57  : cpvTokenDDD_(esConsumes()),
58  cpvTokenDD4hep_(esConsumes()),
59  topoToken_(esConsumes()),
60  geomDetToken_(esConsumes()),
61  ptpToken_(esConsumes()),
62  pixQualityToken_(esConsumes()),
63  stripQualityToken_(esConsumes()),
64  referenceTracker(nullptr),
65  dummyTracker(nullptr),
66  currentTracker(nullptr),
67  theSurveyIndex(0),
68  theSurveyValues(nullptr),
69  theSurveyErrors(nullptr),
70  levelStrings_(cfg.getUntrackedParameter<std::vector<std::string> >("levels")),
71  fromDD4hep_(cfg.getUntrackedParameter<bool>("fromDD4hep")),
72  writeToDB_(cfg.getUntrackedParameter<bool>("writeToDB")),
73  commonTrackerLevel_(align::invalid),
74  moduleListFile_(),
75  moduleList_(0),
76  inputRootFile1_(nullptr),
77  inputRootFile2_(nullptr),
78  inputTree01_(nullptr),
79  inputTree02_(nullptr),
80  inputTree11_(nullptr),
81  inputTree12_(nullptr),
82  m_nBins_(10000),
83  m_rangeLow_(-.1),
84  m_rangeHigh_(.1),
85  firstEvent_(true),
86  m_vtkmap_(13) {
87  moduleListName_ = cfg.getUntrackedParameter<std::string>("moduleList");
88 
89  //input is ROOT
90  inputFilename1_ = cfg.getUntrackedParameter<std::string>("inputROOTFile1");
91  inputFilename2_ = cfg.getUntrackedParameter<std::string>("inputROOTFile2");
92  inputTreenameAlign_ = cfg.getUntrackedParameter<std::string>("treeNameAlign");
93  inputTreenameDeform_ = cfg.getUntrackedParameter<std::string>("treeNameDeform");
94 
95  //output file
96  filename_ = cfg.getUntrackedParameter<std::string>("outputFile");
97  //output dir for surface deformations
98  surfdir_ = cfg.getUntrackedParameter<std::string>("surfDir");
99 
100  weightBy_ = cfg.getUntrackedParameter<std::string>("weightBy");
101  setCommonTrackerSystem_ = cfg.getUntrackedParameter<std::string>("setCommonTrackerSystem");
102  detIdFlag_ = cfg.getUntrackedParameter<bool>("detIdFlag");
103  detIdFlagFile_ = cfg.getUntrackedParameter<std::string>("detIdFlagFile");
104  weightById_ = cfg.getUntrackedParameter<bool>("weightById");
105  weightByIdFile_ = cfg.getUntrackedParameter<std::string>("weightByIdFile");
106 
107  // if want to use, make id cut list
108  if (detIdFlag_) {
109  std::ifstream fin;
110  fin.open(detIdFlagFile_.c_str());
111 
112  while (!fin.eof() && fin.good()) {
113  uint32_t id;
114  fin >> id;
115  detIdFlagVector_.push_back(id);
116  }
117  fin.close();
118  }
119 
120  // turn weightByIdFile into weightByIdVector
121  if (weightById_) {
122  std::ifstream inFile;
123  inFile.open(weightByIdFile_.c_str());
124  while (!inFile.eof()) {
125  unsigned int listId;
126  inFile >> listId;
127  inFile.ignore(256, '\n');
128 
129  weightByIdVector_.push_back(listId);
130  }
131  inFile.close();
132  }
133 
134  //root configuration
135  theFile_ = new TFile(filename_.c_str(), "RECREATE");
136  alignTree_ = new TTree("alignTree",
137  "alignTree"); //,"id:level:mid:mlevel:sublevel:x:y:z:r:phi:a:b:c:dx:dy:dz:dr:dphi:da:db:dc");
138  alignTree_->Branch("id", &id_, "id/I");
139  alignTree_->Branch("badModuleQuality", &badModuleQuality_, "badModuleQuality/I");
140  alignTree_->Branch("inModuleList", &inModuleList_, "inModuleList/I");
141  alignTree_->Branch("level", &level_, "level/I");
142  alignTree_->Branch("mid", &mid_, "mid/I");
143  alignTree_->Branch("mlevel", &mlevel_, "mlevel/I");
144  alignTree_->Branch("sublevel", &sublevel_, "sublevel/I");
145  alignTree_->Branch("x", &xVal_, "x/F");
146  alignTree_->Branch("y", &yVal_, "y/F");
147  alignTree_->Branch("z", &zVal_, "z/F");
148  alignTree_->Branch("r", &rVal_, "r/F");
149  alignTree_->Branch("phi", &phiVal_, "phi/F");
150  alignTree_->Branch("eta", &etaVal_, "eta/F");
151  alignTree_->Branch("alpha", &alphaVal_, "alpha/F");
152  alignTree_->Branch("beta", &betaVal_, "beta/F");
153  alignTree_->Branch("gamma", &gammaVal_, "gamma/F");
154  alignTree_->Branch("dx", &dxVal_, "dx/F");
155  alignTree_->Branch("dy", &dyVal_, "dy/F");
156  alignTree_->Branch("dz", &dzVal_, "dz/F");
157  alignTree_->Branch("dr", &drVal_, "dr/F");
158  alignTree_->Branch("dphi", &dphiVal_, "dphi/F");
159  alignTree_->Branch("dalpha", &dalphaVal_, "dalpha/F");
160  alignTree_->Branch("dbeta", &dbetaVal_, "dbeta/F");
161  alignTree_->Branch("dgamma", &dgammaVal_, "dgamma/F");
162  alignTree_->Branch("du", &duVal_, "du/F");
163  alignTree_->Branch("dv", &dvVal_, "dv/F");
164  alignTree_->Branch("dw", &dwVal_, "dw/F");
165  alignTree_->Branch("da", &daVal_, "da/F");
166  alignTree_->Branch("db", &dbVal_, "db/F");
167  alignTree_->Branch("dg", &dgVal_, "dg/F");
168  alignTree_->Branch("useDetId", &useDetId_, "useDetId/I");
169  alignTree_->Branch("detDim", &detDim_, "detDim/I");
170  alignTree_->Branch("surW", &surWidth_, "surW/F");
171  alignTree_->Branch("surL", &surLength_, "surL/F");
172  alignTree_->Branch("surRot", &surRot_, "surRot[9]/D");
173  alignTree_->Branch("identifiers", &identifiers_, "identifiers[6]/I");
174  alignTree_->Branch("type", &type_, "type/I");
175  alignTree_->Branch("surfDeform", &surfDeform_, "surfDeform[13]/D");
176 
177  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
178  it->setPalette(1);
179  it->addPixel(true);
180  }
181 
182  edm::Service<TFileService> fs;
183  TFileDirectory subDir_All = fs->mkdir("AllSubdetectors");
184  TFileDirectory subDir_PXB = fs->mkdir("PixelBarrel");
185  TFileDirectory subDir_PXF = fs->mkdir("PixelEndcap");
186  for (int ii = 0; ii < 13; ++ii) {
187  std::stringstream histname0;
188  histname0 << "SurfDeform_Par_" << ii;
189  m_h1_[histname0.str()] = subDir_All.make<TH1D>(
190  (histname0.str()).c_str(), (histname0.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
191 
192  std::stringstream histname1;
193  histname1 << "SurfDeform_PixelBarrel_Par_" << ii;
194  m_h1_[histname1.str()] = subDir_PXB.make<TH1D>(
195  (histname1.str()).c_str(), (histname1.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
196 
197  std::stringstream histname2;
198  histname2 << "SurfDeform_PixelEndcap_Par_" << ii;
199  m_h1_[histname2.str()] = subDir_PXF.make<TH1D>(
200  (histname2.str()).c_str(), (histname2.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
201  }
202 }
TrackerGeometryCompare::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: TrackerGeometryCompare.h:93
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:91
TrackerGeometryCompare::levelStrings_
const std::vector< std::string > levelStrings_
Definition: TrackerGeometryCompare.h:106
TrackerGeometryCompare::theSurveyValues
const Alignments * theSurveyValues
Definition: TrackerGeometryCompare.h:102
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:164
TrackerGeometryCompare::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: TrackerGeometryCompare.h:92
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85
TrackerGeometryCompare::dummyTracker
AlignableTracker * dummyTracker
Definition: TrackerGeometryCompare.h:98
TrackerGeometryCompare::m_h1_
std::map< std::string, TH1D * > m_h1_
Definition: TrackerGeometryCompare.h:166
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:130
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TrackerGeometryCompare::stripQualityToken_
const edm::ESGetToken< SiStripQuality, SiStripQualityRcd > stripQualityToken_
Definition: TrackerGeometryCompare.h:95
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
TrackerGeometryCompare::detIdFlagVector_
std::vector< uint32_t > detIdFlagVector_
Definition: TrackerGeometryCompare.h:123
TrackerGeometryCompare::cpvTokenDDD_
const edm::ESGetToken< DDCompactView, IdealGeometryRecord > cpvTokenDDD_
Definition: TrackerGeometryCompare.h:89
TrackerGeometryCompare::weightByIdVector_
std::vector< unsigned int > weightByIdVector_
Definition: TrackerGeometryCompare.h:120
testHGCalDigi_cfg.inFile
inFile
Definition: testHGCalDigi_cfg.py:73
cuy.ii
ii
Definition: cuy.py:589
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:99
TrackerGeometryCompare::cpvTokenDD4hep_
const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord > cpvTokenDD4hep_
Definition: TrackerGeometryCompare.h:90
TrackerGeometryCompare::commonTrackerLevel_
align::StructureType commonTrackerLevel_
Definition: TrackerGeometryCompare.h:124
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:97
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
TrackerGeometryCompare::theSurveyErrors
const SurveyErrors * theSurveyErrors
Definition: TrackerGeometryCompare.h:103
TrackerGeometryCompare::pixQualityToken_
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
Definition: TrackerGeometryCompare.h:94

Member Function Documentation

◆ addSurveyInfo()

void TrackerGeometryCompare::addSurveyInfo ( Alignable *  ali)
private

Definition at line 1008 of file TrackerGeometryCompare.cc.

References Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), relativeConstraints::error, Exception, Alignable::geomDetId(), mps_fire::i, AlignableSurface::length(), Alignments::m_align, SurveyErrors::m_surveyErrors, DetId::rawId(), makeMuonMisalignmentScenario::rot, AlignableSurface::setLength(), Alignable::setSurvey(), AlignableSurface::setWidth(), Alignable::surface(), theSurveyErrors, theSurveyIndex, theSurveyValues, and AlignableSurface::width().

1008  {
1009  const auto& comp = ali->components();
1010 
1011  unsigned int nComp = comp.size();
1012 
1013  for (unsigned int i = 0; i < nComp; ++i)
1014  addSurveyInfo(comp[i]);
1015 
1016  const SurveyError& error = theSurveyErrors->m_surveyErrors[theSurveyIndex];
1017 
1018  if (ali->geomDetId().rawId() != error.rawId() || ali->alignableObjectId() != error.structureType()) {
1019  throw cms::Exception("DatabaseError") << "Error reading survey info from DB. Mismatched id!";
1020  }
1021 
1022  const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[theSurveyIndex].translation();
1023  const CLHEP::HepRotation& rot = theSurveyValues->m_align[theSurveyIndex].rotation();
1024 
1025  AlignableSurface surf(
1026  align::PositionType(pos.x(), pos.y(), pos.z()),
1027  align::RotationType(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
1028 
1029  surf.setWidth(ali->surface().width());
1030  surf.setLength(ali->surface().length());
1031 
1032  ali->setSurvey(new SurveyDet(surf, error.matrix()));
1033 
1034  ++theSurveyIndex;
1035 }
mps_fire.i
i
Definition: mps_fire.py:429
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:132
TrackerGeometryCompare::theSurveyValues
const Alignments * theSurveyValues
Definition: TrackerGeometryCompare.h:102
TrackerGeometryCompare::addSurveyInfo
void addSurveyInfo(Alignable *ali)
Definition: TrackerGeometryCompare.cc:1008
AlignableSurface::width
align::Scalar width() const
Definition: AlignableSurface.h:30
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:249
Alignments::m_align
std::vector< AlignTransform > m_align
Definition: Alignments.h:19
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
Alignable::components
virtual const Alignables & components() const =0
Return vector of all direct components.
Alignable::geomDetId
const DetId & geomDetId() const
Definition: Alignable.h:177
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
Alignable::setSurvey
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:266
AlignableSurface::length
align::Scalar length() const
Definition: AlignableSurface.h:32
TrackerGeometryCompare::theSurveyErrors
const SurveyErrors * theSurveyErrors
Definition: TrackerGeometryCompare.h:103
SurveyErrors::m_surveyErrors
std::vector< SurveyError > m_surveyErrors
Definition: SurveyErrors.h:22

◆ analyze()

void TrackerGeometryCompare::analyze ( const edm::Event &  ,
const edm::EventSetup &  iSetup 
)
overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 246 of file TrackerGeometryCompare.cc.

References AlignableTracker::alignmentErrors(), AlignableTracker::alignments(), cond::service::PoolDBOutputService::beginOfTime(), compareGeometries(), compareSurfaceDeformations(), createROOTGeometry(), currentTracker, firstEvent_, edm::EventSetup::getData(), inputTree11_, inputTree12_, edm::Service< T >::isAvailable(), personalPlayback::level, levelStrings_, m_theLevels, AlignableTracker::objectIdProvider(), referenceTracker, setCommonTrackerSystem(), setCommonTrackerSystem_, AlignableObjectId::stringToId(), topoToken_, cond::service::PoolDBOutputService::writeOneIOV(), and writeToDB_.

246  {
247  if (firstEvent_) {
248  //Retrieve tracker topology from geometry
249  const TrackerTopology* const tTopo = &iSetup.getData(topoToken_);
250 
251  //upload the ROOT geometries
252  createROOTGeometry(iSetup);
253 
254  //setting the levels being used in the geometry comparator
255  edm::LogInfo("TrackerGeometryCompare") << "levels: " << levelStrings_.size();
256  for (const auto& level : levelStrings_) {
257  m_theLevels.push_back(currentTracker->objectIdProvider().stringToId(level));
258  edm::LogInfo("TrackerGeometryCompare") << "level: " << level;
259  edm::LogInfo("TrackerGeometryCompare")
260  << "structure type: " << currentTracker->objectIdProvider().stringToId(level);
261  }
262 
263  //set common tracker system first
264  // if setting the tracker common system
265  if (setCommonTrackerSystem_ != "NONE") {
266  setCommonTrackerSystem();
267  }
268 
269  //compare the goemetries
270  compareGeometries(referenceTracker, currentTracker, tTopo, iSetup);
271  compareSurfaceDeformations(inputTree11_, inputTree12_);
272 
273  //write out ntuple
274  //might be better to do within output module
275 
276  if (writeToDB_) {
277  Alignments* myAlignments = currentTracker->alignments();
278  AlignmentErrorsExtended* myAlignmentErrorsExtended = currentTracker->alignmentErrors();
279 
280  // 2. Store alignment[Error]s to DB
281  edm::Service<cond::service::PoolDBOutputService> poolDbService;
282  // Call service
283  if (!poolDbService.isAvailable()) // Die if not available
284  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
285 
286  poolDbService->writeOneIOV<Alignments>(*myAlignments, poolDbService->beginOfTime(), "TrackerAlignmentRcd");
287  poolDbService->writeOneIOV<AlignmentErrorsExtended>(
288  *myAlignmentErrorsExtended, poolDbService->beginOfTime(), "TrackerAlignmentErrorExtendedRcd");
289  }
290 
291  firstEvent_ = false;
292  }
293 }
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
AlignableTracker::alignments
Alignments * alignments() const override
Return alignments, sorted by DetId.
Definition: AlignableTracker.cc:42
TrackerGeometryCompare::compareGeometries
void compareGeometries(Alignable *refAli, Alignable *curAli, const TrackerTopology *tTopo, const edm::EventSetup &iSetup)
Definition: TrackerGeometryCompare.cc:659
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:91
TrackerGeometryCompare::m_theLevels
std::vector< align::StructureType > m_theLevels
Definition: TrackerGeometryCompare.h:60
TrackerGeometryCompare::levelStrings_
const std::vector< std::string > levelStrings_
Definition: TrackerGeometryCompare.h:106
cond::service::PoolDBOutputService::writeOneIOV
Hash writeOneIOV(const T &payload, Time_t time, const std::string &recordName)
Definition: PoolDBOutputService.h:64
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:99
AlignableTracker::alignmentErrors
AlignmentErrorsExtended * alignmentErrors() const override
Return alignment errors, sorted by DetId.
Definition: AlignableTracker.cc:59
AlignableObjectId::stringToId
align::StructureType stringToId(const char *) const
Definition: AlignableObjectId.cc:284
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:97
TrackerGeometryCompare::createROOTGeometry
void createROOTGeometry(const edm::EventSetup &iSetup)
Definition: TrackerGeometryCompare.cc:295
AlignableTracker::objectIdProvider
const AlignableObjectId & objectIdProvider() const
Return tracker alignable object ID provider derived from the tracker&#39;s geometry.
Definition: AlignableTracker.h:110
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:40
TrackerGeometryCompare::compareSurfaceDeformations
void compareSurfaceDeformations(TTree *_inputTree11, TTree *_inputTree12)
Definition: TrackerGeometryCompare.cc:483

◆ beginJob()

void TrackerGeometryCompare::beginJob ( )
overridevirtual

Read from DB and print survey info.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 226 of file TrackerGeometryCompare.cc.

References firstEvent_.

226 { firstEvent_ = true; }

◆ compareGeometries()

void TrackerGeometryCompare::compareGeometries ( Alignable *  refAli,
Alignable *  curAli,
const TrackerTopology *  tTopo,
const edm::EventSetup &  iSetup 
)
private

Definition at line 659 of file TrackerGeometryCompare.cc.

References Alignable::alignableObjectId(), RPCNoise_example::check, Alignable::components(), change_name::diff, align::diffAlignables(), HGC3DClusterGenMatchSelector_cfi::dR, MillePedeFileConverter_cfg::e, Exception, fillTree(), Alignable::geomDetId(), mps_fire::i, Alignable::id(), m_theLevels, PV3DBase< T, PVType, FrameType >::mag(), align::moveAlignable(), DetId::rawId(), makeMuonMisalignmentScenario::rot, tolerance, weightBy_, weightById_, and weightByIdVector_.

Referenced by analyze().

662  {
663  using namespace align;
664 
665  const auto& refComp = refAli->components();
666  const auto& curComp = curAli->components();
667 
668  unsigned int nComp = refComp.size();
669  //only perform for designate levels
670  bool useLevel = false;
671  for (unsigned int i = 0; i < m_theLevels.size(); ++i) {
672  if (refAli->alignableObjectId() == m_theLevels[i])
673  useLevel = true;
674  }
675 
676  //another added level for difference between det and detunit
677  //if ((refAli->alignableObjectId()==2)&&(nComp == 1)) useLevel = false;
678 
679  //coordinate matching, etc etc
680  if (useLevel) {
681  DetId detid(refAli->id());
682 
683  CLHEP::Hep3Vector Rtotal, Wtotal, lRtotal, lWtotal;
684  Rtotal.set(0., 0., 0.);
685  Wtotal.set(0., 0., 0.);
686  lRtotal.set(0., 0., 0.);
687  lWtotal.set(0., 0., 0.);
688 
689  bool converged = false;
690 
691  AlgebraicVector diff, check;
692 
693  for (int i = 0; i < 100; i++) {
694  // Get differences between alignments for rotations and translations
695  // both local and global
696  diff = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
697 
698  // 'diffAlignables' returns 'refAli - curAli' for translations and 'curAli - refAli' for rotations.
699  // The plan is to unify this at some point, but a simple change of the sign for one of them was postponed
700  // to do some further checks to understand the rotations better
701  //Updated July 2018: as requested the sign in the translations has been changed to match the one in rotations. A test was done to change the diffAlignables function and solve the issue there, but proved quite time consuming. To unify the sign convention in the least amount of time the choice was made to change the sign here.
702  CLHEP::Hep3Vector dR(-diff[0], -diff[1], -diff[2]);
703  CLHEP::Hep3Vector dW(diff[3], diff[4], diff[5]);
704  CLHEP::Hep3Vector dRLocal(-diff[6], -diff[7], -diff[8]);
705  CLHEP::Hep3Vector dWLocal(diff[9], diff[10], diff[11]);
706 
707  // Translations
708  Rtotal += dR;
709  lRtotal += dRLocal;
710 
711  //Rotations
712  CLHEP::HepRotation rot(Wtotal.unit(), Wtotal.mag());
713  CLHEP::HepRotation drot(dW.unit(), dW.mag());
714  rot *= drot;
715  Wtotal.set(rot.axis().x() * rot.delta(), rot.axis().y() * rot.delta(), rot.axis().z() * rot.delta());
716 
717  CLHEP::HepRotation rotLocal(lWtotal.unit(), lWtotal.mag());
718  CLHEP::HepRotation drotLocal(dWLocal.unit(), dWLocal.mag());
719  rotLocal *= drotLocal;
720  lWtotal.set(rotLocal.axis().x() * rotLocal.delta(),
721  rotLocal.axis().y() * rotLocal.delta(),
722  rotLocal.axis().z() * rotLocal.delta());
723 
724  // Move current alignable by shift and check if difference
725  // is smaller than tolerance value
726  // if true, break the loop
727  align::moveAlignable(curAli, diff);
728  float tolerance = 1e-7;
729  check = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
730  align::GlobalVector checkR(check[0], check[1], check[2]);
731  align::GlobalVector checkW(check[3], check[4], check[5]);
732  if ((checkR.mag() < tolerance) && (checkW.mag() < tolerance)) {
733  converged = true;
734  break;
735  }
736  }
737 
738  // give an exception if difference has not fallen below tolerance level
739  // i.e. method has not converged
740  if (!converged) {
741  edm::LogInfo("TrackerGeometryCompare")
742  << "Tolerance Exceeded!(alObjId: " << refAli->alignableObjectId()
743  << ", rawId: " << refAli->geomDetId().rawId() << ", subdetId: " << detid.subdetId() << "): " << diff << check;
744  throw cms::Exception("Tolerance in TrackerGeometryCompare exceeded");
745  }
746  AlgebraicVector TRtot(12);
747  // global
748  TRtot(1) = Rtotal.x();
749  TRtot(2) = Rtotal.y();
750  TRtot(3) = Rtotal.z();
751  TRtot(4) = Wtotal.x();
752  TRtot(5) = Wtotal.y();
753  TRtot(6) = Wtotal.z();
754  // local
755  TRtot(7) = lRtotal.x();
756  TRtot(8) = lRtotal.y();
757  TRtot(9) = lRtotal.z();
758  TRtot(10) = lWtotal.x();
759  TRtot(11) = lWtotal.y();
760  TRtot(12) = lWtotal.z();
761 
762  fillTree(refAli, TRtot, tTopo, iSetup);
763  }
764 
765  // another added level for difference between det and detunit
766  for (unsigned int i = 0; i < nComp; ++i)
767  compareGeometries(refComp[i], curComp[i], tTopo, iSetup);
768 }
mps_fire.i
i
Definition: mps_fire.py:429
RPCNoise_example.check
check
Definition: RPCNoise_example.py:71
TrackerGeometryCompare::compareGeometries
void compareGeometries(Alignable *refAli, Alignable *curAli, const TrackerTopology *tTopo, const edm::EventSetup &iSetup)
Definition: TrackerGeometryCompare.cc:659
TrackerGeometryCompare::m_theLevels
std::vector< align::StructureType > m_theLevels
Definition: TrackerGeometryCompare.h:60
tolerance
const double tolerance
Definition: HGCalGeomParameters.cc:29
align::diffAlignables
AlgebraicVector diffAlignables(Alignable *refAli, Alignable *curAli, const std::string &weightBy, bool weightById, const std::vector< unsigned int > &weightByIdVector)
Definition: AlignTools.cc:10
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
TrackerGeometryCompare::weightByIdVector_
std::vector< unsigned int > weightByIdVector_
Definition: TrackerGeometryCompare.h:120
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
change_name.diff
diff
Definition: change_name.py:13
Alignable::components
virtual const Alignables & components() const =0
Return vector of all direct components.
TrackerGeometryCompare::fillTree
void fillTree(Alignable *refAli, const AlgebraicVector &diff, const TrackerTopology *tTopo, const edm::EventSetup &iSetup)
Definition: TrackerGeometryCompare.cc:847
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
DetId
Definition: DetId.h:17
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:13
Alignable::id
align::ID id() const
Return the ID of Alignable, i.e. DetId of &#39;first&#39; component GeomDet(Unit).
Definition: Alignable.h:180
Alignable::geomDetId
const DetId & geomDetId() const
Definition: Alignable.h:177
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
HGC3DClusterGenMatchSelector_cfi.dR
dR
Definition: HGC3DClusterGenMatchSelector_cfi.py:7
align::moveAlignable
void moveAlignable(Alignable *ali, AlgebraicVector diff)
Moves the alignable by the AlgebraicVector.
Definition: AlignTools.cc:84

◆ compareSurfaceDeformations()

void TrackerGeometryCompare::compareSurfaceDeformations ( TTree *  _inputTree11,
TTree *  _inputTree12 
)
private

Definition at line 483 of file TrackerGeometryCompare.cc.

References cuy::ii, inputFilename1_, inputFilename2_, inputRootFile1_, inputRootFile2_, inputTreenameDeform_, createfilelist::int, m_h1_, m_nBins_, m_rangeHigh_, m_rangeLow_, m_vtkmap_, phase_, surfDeform_, SiPixelPI::two, and SiPixelPI::zero.

Referenced by analyze().

483  {
484  if (inputFilename1_ != "IDEAL" && inputFilename2_ != "IDEAL") {
485  int inputRawid1;
486  int inputRawid2;
487  int inputSubdetid1, inputSubdetid2;
488  int inputDtype1, inputDtype2;
489  std::vector<double> inputDpar1;
490  std::vector<double> inputDpar2;
491  std::vector<double>* p_inputDpar1 = &inputDpar1;
492  std::vector<double>* p_inputDpar2 = &inputDpar2;
493 
494  TTree* refTree = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
495  refTree->SetBranchAddress("irawid", &inputRawid1);
496  refTree->SetBranchAddress("subdetid", &inputSubdetid1);
497  refTree->SetBranchAddress("dtype", &inputDtype1);
498  refTree->SetBranchAddress("dpar", &p_inputDpar1);
499 
500  TTree* curTree = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
501  curTree->SetBranchAddress("irawid", &inputRawid2);
502  curTree->SetBranchAddress("subdetid", &inputSubdetid2);
503  curTree->SetBranchAddress("dtype", &inputDtype2);
504  curTree->SetBranchAddress("dpar", &p_inputDpar2);
505 
506  unsigned int nEntries11 = refTree->GetEntries();
507  unsigned int nEntries12 = curTree->GetEntries();
508 
509  if (nEntries11 != nEntries12) {
510  edm::LogError("TrackerGeometryCompare") << " Surface deformation parameters in two geometries differ!\n";
511  return;
512  }
513 
514  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
515  refTree->GetEntry(iEntry);
516  curTree->GetEntry(iEntry);
517  for (int ii = 0; ii < 13; ++ii) {
518  surfDeform_[ii] = -1.0;
519  }
520  for (int npar = 0; npar < int(inputDpar2.size()); ++npar) {
521  if (inputRawid1 == inputRawid2) {
522  surfDeform_[npar] = inputDpar2.at(npar) - inputDpar1.at(npar);
523  std::stringstream histname0;
524  histname0 << "SurfDeform_Par_" << npar;
525  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
526  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
527  if (inputSubdetid1 == 1 && inputSubdetid2 == 1) {
528  std::stringstream histname1;
529  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
530  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
531  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
532  }
533  if (inputSubdetid1 == 2 && inputSubdetid2 == 2) {
534  std::stringstream histname2;
535  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
536  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
537  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
538  }
539  if (phase_ < SiPixelPI::phase::two) {
540  if (phase_ == SiPixelPI::phase::zero) {
541  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
542  } else {
543  // old-style tracker map does not support Phase1 Pixel
544  if (!(inputSubdetid1 == 1 && inputSubdetid2 == 1) && !(inputSubdetid1 == 2 && inputSubdetid2 == 2)) {
545  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
546  }
547  }
548  }
549  }
550  }
551  }
552 
553  } else if (inputFilename1_ == "IDEAL" && inputFilename2_ != "IDEAL") {
554  int inputRawid2;
555  int inputSubdetid2;
556  int inputDtype2;
557  std::vector<double> inputDpar2;
558  std::vector<double>* p_inputDpar2 = &inputDpar2;
559 
560  TTree* curTree = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
561  curTree->SetBranchAddress("irawid", &inputRawid2);
562  curTree->SetBranchAddress("subdetid", &inputSubdetid2);
563  curTree->SetBranchAddress("dtype", &inputDtype2);
564  curTree->SetBranchAddress("dpar", &p_inputDpar2);
565 
566  unsigned int nEntries12 = curTree->GetEntries();
567  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
568  curTree->GetEntry(iEntry);
569  for (int ii = 0; ii < 12; ++ii) {
570  surfDeform_[ii] = -1.0;
571  }
572  for (int npar = 0; npar < int(inputDpar2.size()); ++npar) {
573  surfDeform_[npar] = inputDpar2.at(npar);
574  std::stringstream histname0;
575  histname0 << "SurfDeform_Par_" << npar;
576  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
577  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
578  if (inputSubdetid2 == 1) {
579  std::stringstream histname1;
580  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
581  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
582  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
583  }
584  if (inputSubdetid2 == 2) {
585  std::stringstream histname2;
586  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
587  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
588  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
589  }
590  if (phase_ < SiPixelPI::phase::two) {
591  if (phase_ == SiPixelPI::phase::zero) {
592  (m_vtkmap_.at(npar)).fill_current_val(inputRawid2, surfDeform_[npar]);
593  } else {
594  // old-style tracker map does not support Phase1 Pixel
595  if (inputSubdetid2 != 1 && inputSubdetid2 != 2) {
596  (m_vtkmap_.at(npar)).fill_current_val(inputRawid2, surfDeform_[npar]);
597  }
598  }
599  }
600  }
601  }
602 
603  } else if (inputFilename1_ != "IDEAL" && inputFilename2_ == "IDEAL") {
604  int inputRawid1;
605  int inputSubdetid1;
606  int inputDtype1;
607  std::vector<double> inputDpar1;
608  std::vector<double>* p_inputDpar1 = &inputDpar1;
609 
610  TTree* refTree = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
611  refTree->SetBranchAddress("irawid", &inputRawid1);
612  refTree->SetBranchAddress("subdetid", &inputSubdetid1);
613  refTree->SetBranchAddress("dtype", &inputDtype1);
614  refTree->SetBranchAddress("dpar", &p_inputDpar1);
615  unsigned int nEntries11 = refTree->GetEntries();
616 
617  for (unsigned int iEntry = 0; iEntry < nEntries11; ++iEntry) {
618  refTree->GetEntry(iEntry);
619  for (int ii = 0; ii < 12; ++ii) {
620  surfDeform_[ii] = -1.0;
621  }
622  for (int npar = 0; npar < int(inputDpar1.size()); ++npar) {
623  surfDeform_[npar] = -inputDpar1.at(npar);
624  std::stringstream histname0;
625  histname0 << "SurfDeform_Par_" << npar;
626  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
627  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
628  if (inputSubdetid1 == 1) {
629  std::stringstream histname1;
630  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
631  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
632  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
633  }
634  if (inputSubdetid1 == 2) {
635  std::stringstream histname2;
636  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
637  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
638  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
639  }
640  if (phase_ < SiPixelPI::phase::two) {
641  if (phase_ == SiPixelPI::phase::zero) {
642  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
643  } else {
644  // old-style tracker map does not support Phase1 Pixel
645  if (inputSubdetid1 != 1 && inputSubdetid1 != 2) {
646  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
647  }
648  }
649  }
650  }
651  }
652 
653  } else if (inputFilename1_ == "IDEAL" && inputFilename2_ == "IDEAL") {
654  edm::LogInfo("TrackerGeometryCompare") << ">>>> Comparing IDEAL with IDEAL: nothing to do! <<<<\n";
655  }
656  return;
657 }
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:164
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
TrackerGeometryCompare::m_h1_
std::map< std::string, TH1D * > m_h1_
Definition: TrackerGeometryCompare.h:166
cuy.ii
ii
Definition: cuy.py:589
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131

◆ createROOTGeometry()

void TrackerGeometryCompare::createROOTGeometry ( const edm::EventSetup &  iSetup)
private

Definition at line 295 of file TrackerGeometryCompare.cc.

References GeometryAligner::applyAlignments(), TrackerGeomBuilderFromGeometricDet::build(), cpvTokenDD4hep_, cpvTokenDDD_, SurfaceDeformationFactory::create(), currentTracker, Alignable::deepComponents(), fromDD4hep_, geomDetToken_, edm::EventSetup::getData(), edm::EventSetup::getTransientHandle(), mps_fire::i, inputFilename1_, inputFilename2_, inputRootFile1_, inputRootFile2_, inputTree01_, inputTree02_, inputTree11_, inputTree12_, inputTreenameAlign_, inputTreenameDeform_, TrackerGeometry::isThere(), mps_splice::line, Alignments::m_align, AlignmentErrorsExtended::m_alignError, moduleList_, moduleListFile_, moduleListName_, SiPixelPI::one, GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, phase_, ptpToken_, referenceTracker, jetUpdater_cfi::sort, AlCaHLTBitMon_QueryRunRegistry::string, topoToken_, SiPixelPI::two, and SiPixelPI::zero.

Referenced by analyze().

295  {
296  int inputRawId1, inputRawId2;
297  double inputX1, inputY1, inputZ1, inputX2, inputY2, inputZ2;
298  double inputAlpha1, inputBeta1, inputGamma1, inputAlpha2, inputBeta2, inputGamma2;
299 
300  //Retrieve tracker topology from geometry
301  const TrackerTopology* const tTopo = &iSetup.getData(topoToken_);
302 
303  // Fill module IDs from file into a list
304  moduleListFile_.open(moduleListName_);
305  if (moduleListFile_.is_open()) {
306  std::string line;
307  while (!moduleListFile_.eof()) {
308  std::getline(moduleListFile_, line);
309  moduleList_.push_back(std::atoi(line.c_str()));
310  }
311  } else {
312  edm::LogInfo("TrackerGeometryCompare") << "Error: Module list not found! Please verify that given list exists!";
313  }
314 
315  //declare alignments
316  Alignments* alignments1 = new Alignments();
317  AlignmentErrorsExtended* alignmentErrors1 = new AlignmentErrorsExtended();
318  if (inputFilename1_ != "IDEAL") {
319  inputRootFile1_ = new TFile(inputFilename1_.c_str());
320  TTree* inputTree01_ = (TTree*)inputRootFile1_->Get(inputTreenameAlign_.c_str());
321  inputTree01_->SetBranchAddress("rawid", &inputRawId1);
322  inputTree01_->SetBranchAddress("x", &inputX1);
323  inputTree01_->SetBranchAddress("y", &inputY1);
324  inputTree01_->SetBranchAddress("z", &inputZ1);
325  inputTree01_->SetBranchAddress("alpha", &inputAlpha1);
326  inputTree01_->SetBranchAddress("beta", &inputBeta1);
327  inputTree01_->SetBranchAddress("gamma", &inputGamma1);
328 
329  int nEntries1 = inputTree01_->GetEntries();
330  //fill alignments
331  for (int i = 0; i < nEntries1; ++i) {
332  inputTree01_->GetEntry(i);
333  CLHEP::Hep3Vector translation1(inputX1, inputY1, inputZ1);
334  CLHEP::HepEulerAngles eulerangles1(inputAlpha1, inputBeta1, inputGamma1);
335  uint32_t detid1 = inputRawId1;
336  AlignTransform transform1(translation1, eulerangles1, detid1);
337  alignments1->m_align.push_back(transform1);
338 
339  // dummy errors
340  // APE matrix was 3x3, now it's 6x6 (because of muons), see PR #6483
341  CLHEP::HepSymMatrix clhepSymMatrix(6, 0);
342  AlignTransformErrorExtended transformError(clhepSymMatrix, detid1);
343  alignmentErrors1->m_alignError.push_back(transformError);
344  }
345 
346  // to get the right order, sort by rawId
347  std::sort(alignments1->m_align.begin(), alignments1->m_align.end());
348  std::sort(alignmentErrors1->m_alignError.begin(), alignmentErrors1->m_alignError.end());
349  }
350  //------------------
351  Alignments* alignments2 = new Alignments();
352  AlignmentErrorsExtended* alignmentErrors2 = new AlignmentErrorsExtended();
353  if (inputFilename2_ != "IDEAL") {
354  inputRootFile2_ = new TFile(inputFilename2_.c_str());
355  TTree* inputTree02_ = (TTree*)inputRootFile2_->Get(inputTreenameAlign_.c_str());
356  inputTree02_->SetBranchAddress("rawid", &inputRawId2);
357  inputTree02_->SetBranchAddress("x", &inputX2);
358  inputTree02_->SetBranchAddress("y", &inputY2);
359  inputTree02_->SetBranchAddress("z", &inputZ2);
360  inputTree02_->SetBranchAddress("alpha", &inputAlpha2);
361  inputTree02_->SetBranchAddress("beta", &inputBeta2);
362  inputTree02_->SetBranchAddress("gamma", &inputGamma2);
363 
364  int nEntries2 = inputTree02_->GetEntries();
365  //fill alignments
366  for (int i = 0; i < nEntries2; ++i) {
367  inputTree02_->GetEntry(i);
368  CLHEP::Hep3Vector translation2(inputX2, inputY2, inputZ2);
369  CLHEP::HepEulerAngles eulerangles2(inputAlpha2, inputBeta2, inputGamma2);
370  uint32_t detid2 = inputRawId2;
371  AlignTransform transform2(translation2, eulerangles2, detid2);
372  alignments2->m_align.push_back(transform2);
373 
374  // dummy errors
375  // APE matrix was 3x3, now it's 6x6 (because of muons), see PR #6483
376  CLHEP::HepSymMatrix clhepSymMatrix(6, 0);
377  AlignTransformErrorExtended transformError(clhepSymMatrix, detid2);
378  alignmentErrors2->m_alignError.push_back(transformError);
379  }
380 
381  // to get the right order, sort by rawId
382  std::sort(alignments2->m_align.begin(), alignments2->m_align.end());
383  std::sort(alignmentErrors2->m_alignError.begin(), alignmentErrors2->m_alignError.end());
384  }
385 
386  //accessing the initial geometry
387  if (!fromDD4hep_) {
388  edm::ESTransientHandle<DDCompactView> cpv = iSetup.getTransientHandle(cpvTokenDDD_);
389  } else {
390  edm::ESTransientHandle<cms::DDCompactView> cpv = iSetup.getTransientHandle(cpvTokenDD4hep_);
391  }
392 
393  const GeometricDet* theGeometricDet = &iSetup.getData(geomDetToken_);
394  const PTrackerParameters* ptp = &iSetup.getData(ptpToken_);
395  TrackerGeomBuilderFromGeometricDet trackerBuilder;
396 
397  //reference tracker
398  TrackerGeometry* theRefTracker = trackerBuilder.build(theGeometricDet, *ptp, tTopo);
399  if (inputFilename1_ != "IDEAL") {
400  GeometryAligner aligner1;
401  aligner1.applyAlignments<TrackerGeometry>(
402  &(*theRefTracker), &(*alignments1), &(*alignmentErrors1), AlignTransform());
403  }
404  referenceTracker = new AlignableTracker(&(*theRefTracker), tTopo);
405  //referenceTracker->setSurfaceDeformation(surfDef1, true) ;
406 
407  //set tracker pixel Phase, assume current tracker is the same Phase
408  if (theRefTracker->isThere(GeomDetEnumerators::P2PXB) || theRefTracker->isThere(GeomDetEnumerators::P2PXEC)) {
409  phase_ = SiPixelPI::phase::two;
410  } else if (theRefTracker->isThere(GeomDetEnumerators::P1PXB) || theRefTracker->isThere(GeomDetEnumerators::P1PXEC)) {
411  phase_ = SiPixelPI::phase::one;
412  } else {
413  phase_ = SiPixelPI::phase::zero;
414  }
415 
416  int inputRawid1;
417  int inputRawid2;
418  int inputDtype1, inputDtype2;
419  std::vector<double> inputDpar1;
420  std::vector<double> inputDpar2;
421  std::vector<double>* p_inputDpar1 = &inputDpar1;
422  std::vector<double>* p_inputDpar2 = &inputDpar2;
423 
424  const auto& comp1 = referenceTracker->deepComponents();
425 
426  SurfaceDeformation* surfDef1;
427  if (inputFilename1_ != "IDEAL") {
428  TTree* inputTree11_ = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
429  inputTree11_->SetBranchAddress("irawid", &inputRawid1);
430  inputTree11_->SetBranchAddress("dtype", &inputDtype1);
431  inputTree11_->SetBranchAddress("dpar", &p_inputDpar1);
432 
433  unsigned int nEntries11 = inputTree11_->GetEntries();
434  edm::LogInfo("TrackerGeometryCompare") << " nentries11 = " << nEntries11;
435  for (unsigned int iEntry = 0; iEntry < nEntries11; ++iEntry) {
436  inputTree11_->GetEntry(iEntry);
437 
438  surfDef1 = SurfaceDeformationFactory::create(inputDtype1, inputDpar1);
439 
440  if (int(comp1[iEntry]->id()) == inputRawid1) {
441  comp1[iEntry]->setSurfaceDeformation(surfDef1, true);
442  }
443  }
444  }
445 
446  //currernt tracker
447  TrackerGeometry* theCurTracker = trackerBuilder.build(&*theGeometricDet, *ptp, tTopo);
448  if (inputFilename2_ != "IDEAL") {
449  GeometryAligner aligner2;
450  aligner2.applyAlignments<TrackerGeometry>(
451  &(*theCurTracker), &(*alignments2), &(*alignmentErrors2), AlignTransform());
452  }
453  currentTracker = new AlignableTracker(&(*theCurTracker), tTopo);
454 
455  const auto& comp2 = currentTracker->deepComponents();
456 
457  SurfaceDeformation* surfDef2;
458  if (inputFilename2_ != "IDEAL") {
459  TTree* inputTree12_ = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
460  inputTree12_->SetBranchAddress("irawid", &inputRawid2);
461  inputTree12_->SetBranchAddress("dtype", &inputDtype2);
462  inputTree12_->SetBranchAddress("dpar", &p_inputDpar2);
463 
464  unsigned int nEntries12 = inputTree12_->GetEntries();
465  edm::LogInfo("TrackerGeometryCompare") << " nentries12 = " << nEntries12;
466  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
467  inputTree12_->GetEntry(iEntry);
468 
469  surfDef2 = SurfaceDeformationFactory::create(inputDtype2, inputDpar2);
470 
471  if (int(comp2[iEntry]->id()) == inputRawid2) {
472  comp2[iEntry]->setSurfaceDeformation(surfDef2, true);
473  }
474  }
475  }
476 
477  delete alignments1;
478  delete alignmentErrors1;
479  delete alignments2;
480  delete alignmentErrors2;
481 }
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
TrackerGeometryCompare::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: TrackerGeometryCompare.h:93
mps_fire.i
i
Definition: mps_fire.py:429
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:91
GeometryAligner
Class to update a given geometry with a set of alignments.
Definition: GeometryAligner.h:33
TrackerGeometryCompare::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: TrackerGeometryCompare.h:92
TrackerGeomBuilderFromGeometricDet::build
TrackerGeometry * build(const GeometricDet *gd, const PTrackerParameters &ptp, const TrackerTopology *tTopo)
Definition: TrackerGeomBuilderFromGeometricDet.cc:44
Alignable::deepComponents
const Alignables & deepComponents() const
Definition: Alignable.h:72
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:130
Alignments::m_align
std::vector< AlignTransform > m_align
Definition: Alignments.h:19
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
GeometryAligner::applyAlignments
void applyAlignments(const C *geometry, const Alignments *alignments, const AlignmentErrorsExtended *alignmentErrors, const AlignTransform &globalCoordinates)
Definition: GeometryAligner.h:52
TrackerGeometry::isThere
bool isThere(GeomDetEnumerators::SubDetector subdet) const
Definition: TrackerGeometry.cc:219
TrackerGeometryCompare::cpvTokenDDD_
const edm::ESGetToken< DDCompactView, IdealGeometryRecord > cpvTokenDDD_
Definition: TrackerGeometryCompare.h:89
mps_splice.line
line
Definition: mps_splice.py:76
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:99
AlignmentErrorsExtended::m_alignError
std::vector< AlignTransformErrorExtended > m_alignError
Definition: AlignmentErrorsExtended.h:19
TrackerGeometryCompare::cpvTokenDD4hep_
const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord > cpvTokenDD4hep_
Definition: TrackerGeometryCompare.h:90
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:97
edm::EventSetup::getTransientHandle
ESTransientHandle< T > getTransientHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:141
SurfaceDeformationFactory::create
SurfaceDeformation * create(int type, const std::vector< double > &params)
Definition: SurfaceDeformationFactory.cc:40

◆ diffCommonTrackerSystem()

void TrackerGeometryCompare::diffCommonTrackerSystem ( Alignable *  refAli,
Alignable *  curAli 
)
private

Definition at line 813 of file TrackerGeometryCompare.cc.

References Alignable::alignableObjectId(), commonTrackerLevel_, Alignable::components(), change_name::diff, align::diffAlignables(), HGC3DClusterGenMatchSelector_cfi::dR, Alignable::globalPosition(), mps_fire::i, makeMuonMisalignmentScenario::rot, TrackerCommonCM_, TrackerCommonR_, TrackerCommonT_, weightBy_, weightById_, and weightByIdVector_.

Referenced by setCommonTrackerSystem().

813  {
814  const auto& refComp = refAli->components();
815  const auto& curComp = curAli->components();
816  unsigned int nComp = refComp.size();
817  //only perform for designate levels
818  bool useLevel = false;
819  if (refAli->alignableObjectId() == commonTrackerLevel_)
820  useLevel = true;
821 
822  //useLevel = false;
823  if (useLevel) {
824  CLHEP::Hep3Vector Rtotal, Wtotal;
825  Rtotal.set(0., 0., 0.);
826  Wtotal.set(0., 0., 0.);
827 
828  AlgebraicVector diff = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
829  CLHEP::Hep3Vector dR(diff[0], diff[1], diff[2]);
830  Rtotal += dR;
831  CLHEP::Hep3Vector dW(diff[3], diff[4], diff[5]);
832  CLHEP::HepRotation rot(Wtotal.unit(), Wtotal.mag());
833  CLHEP::HepRotation drot(dW.unit(), dW.mag());
834  rot *= drot;
835  Wtotal.set(rot.axis().x() * rot.delta(), rot.axis().y() * rot.delta(), rot.axis().z() * rot.delta());
836 
837  TrackerCommonT_ = align::GlobalVector(Rtotal.x(), Rtotal.y(), Rtotal.z());
838  TrackerCommonR_ = align::GlobalVector(Wtotal.x(), Wtotal.y(), Wtotal.z());
839  TrackerCommonCM_ = curAli->globalPosition();
840 
841  } else {
842  for (unsigned int i = 0; i < nComp; ++i)
843  diffCommonTrackerSystem(refComp[i], curComp[i]);
844  }
845 }
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometryCompare::TrackerCommonR_
align::GlobalVector TrackerCommonR_
Definition: TrackerGeometryCompare.h:126
align::GlobalVector
Vector3DBase< Scalar, GlobalTag > GlobalVector
Definition: Definitions.h:31
align::diffAlignables
AlgebraicVector diffAlignables(Alignable *refAli, Alignable *curAli, const std::string &weightBy, bool weightById, const std::vector< unsigned int > &weightByIdVector)
Definition: AlignTools.cc:10
TrackerGeometryCompare::diffCommonTrackerSystem
void diffCommonTrackerSystem(Alignable *refAli, Alignable *curAli)
Definition: TrackerGeometryCompare.cc:813
Alignable::globalPosition
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:135
TrackerGeometryCompare::weightByIdVector_
std::vector< unsigned int > weightByIdVector_
Definition: TrackerGeometryCompare.h:120
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
change_name.diff
diff
Definition: change_name.py:13
TrackerGeometryCompare::TrackerCommonT_
align::GlobalVector TrackerCommonT_
Definition: TrackerGeometryCompare.h:125
Alignable::components
virtual const Alignables & components() const =0
Return vector of all direct components.
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:13
HGC3DClusterGenMatchSelector_cfi.dR
dR
Definition: HGC3DClusterGenMatchSelector_cfi.py:7
TrackerGeometryCompare::commonTrackerLevel_
align::StructureType commonTrackerLevel_
Definition: TrackerGeometryCompare.h:124
TrackerGeometryCompare::TrackerCommonCM_
align::PositionType TrackerCommonCM_
Definition: TrackerGeometryCompare.h:127

◆ endJob()

void TrackerGeometryCompare::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 228 of file TrackerGeometryCompare.cc.

References alignTree_, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, m_vtkmap_, AlCaHLTBitMon_QueryRunRegistry::string, surfdir_, and theFile_.

Referenced by o2olib.O2ORunMgr::executeJob().

228  {
229  int iname(0);
230  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
231  std::stringstream mapname;
232  mapname << surfdir_ << "/TkMap_SurfDeform_0" << iname << ".png";
233  it->save(true, 0, 0, mapname.str());
234  mapname.str(std::string());
235  mapname.clear();
236  mapname << surfdir_ << "/TkMap_SurfDeform_1" << iname << ".pdf";
237  it->save(true, 0, 0, mapname.str());
238  ++iname;
239  }
240 
241  theFile_->cd();
242  alignTree_->Write();
243  theFile_->Close();
244 }
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:164
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48

◆ fillDescriptions()

void TrackerGeometryCompare::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 204 of file TrackerGeometryCompare.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

204  {
205  edm::ParameterSetDescription desc;
206  desc.setComment("Validates alignment payloads by comparing positions of tracker modules positiona and orientations");
207  desc.addUntracked<std::vector<std::string> >("levels", {});
208  desc.addUntracked<bool>("fromDD4hep", false);
209  desc.addUntracked<bool>("writeToDB", false);
210  desc.addUntracked<std::string>("moduleList", "moduleList.txt");
211  desc.addUntracked<std::string>("inputROOTFile1", "IDEAL");
212  desc.addUntracked<std::string>("inputROOTFile2", "idealtracker2.root");
213  desc.addUntracked<std::string>("treeNameAlign", "alignTree");
214  desc.addUntracked<std::string>("treeNameDeform", "alignTreeDeformations");
215  desc.addUntracked<std::string>("outputFile", "output.root");
216  desc.addUntracked<std::string>("surfDir", ".");
217  desc.addUntracked<std::string>("weightBy", "DetUnit");
218  desc.addUntracked<std::string>("setCommonTrackerSystem", "NONE");
219  desc.addUntracked<bool>("detIdFlag", false);
220  desc.addUntracked<std::string>("detIdFlagFile", "blah.txt");
221  desc.addUntracked<bool>("weightById", false);
222  desc.addUntracked<std::string>("weightByIdFile", "blah2.txt");
223  descriptions.addWithDefaultLabel(desc);
224 }
edm::ConfigurationDescriptions::addWithDefaultLabel
void addWithDefaultLabel(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:87
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256

◆ fillIdentifiers()

void TrackerGeometryCompare::fillIdentifiers ( int  subdetlevel,
int  rawid,
const TrackerTopology *  tTopo 
)
private

Definition at line 1049 of file TrackerGeometryCompare.cc.

References identifiers_, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfModule(), TrackerTopology::pxfPanel(), TrackerTopology::pxfSide(), TrackerTopology::tecModule(), TrackerTopology::tecPetalInfo(), TrackerTopology::tecRing(), TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tibModule(), TrackerTopology::tibStringInfo(), TrackerTopology::tidModuleInfo(), TrackerTopology::tidRing(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), TrackerTopology::tobLayer(), TrackerTopology::tobModule(), and TrackerTopology::tobRodInfo().

Referenced by fillTree().

1049  {
1050  switch (subdetlevel) {
1051  case 1: {
1052  identifiers_[0] = tTopo->pxbModule(rawid);
1053  identifiers_[1] = tTopo->pxbLadder(rawid);
1054  identifiers_[2] = tTopo->pxbLayer(rawid);
1055  identifiers_[3] = 999;
1056  identifiers_[4] = 999;
1057  identifiers_[5] = 999;
1058  break;
1059  }
1060  case 2: {
1061  identifiers_[0] = tTopo->pxfModule(rawid);
1062  identifiers_[1] = tTopo->pxfPanel(rawid);
1063  identifiers_[2] = tTopo->pxfBlade(rawid);
1064  identifiers_[3] = tTopo->pxfDisk(rawid);
1065  identifiers_[4] = tTopo->pxfSide(rawid);
1066  identifiers_[5] = 999;
1067  break;
1068  }
1069  case 3: {
1070  identifiers_[0] = tTopo->tibModule(rawid);
1071  identifiers_[1] = tTopo->tibStringInfo(rawid)[0];
1072  identifiers_[2] = tTopo->tibStringInfo(rawid)[1];
1073  identifiers_[3] = tTopo->tibStringInfo(rawid)[2];
1074  identifiers_[4] = tTopo->tibLayer(rawid);
1075  identifiers_[5] = 999;
1076  break;
1077  }
1078  case 4: {
1079  identifiers_[0] = tTopo->tidModuleInfo(rawid)[0];
1080  identifiers_[1] = tTopo->tidModuleInfo(rawid)[1];
1081  identifiers_[2] = tTopo->tidRing(rawid);
1082  identifiers_[3] = tTopo->tidWheel(rawid);
1083  identifiers_[4] = tTopo->tidSide(rawid);
1084  identifiers_[5] = 999;
1085  break;
1086  }
1087  case 5: {
1088  identifiers_[0] = tTopo->tobModule(rawid);
1089  identifiers_[1] = tTopo->tobRodInfo(rawid)[0];
1090  identifiers_[2] = tTopo->tobRodInfo(rawid)[1];
1091  identifiers_[3] = tTopo->tobLayer(rawid);
1092  identifiers_[4] = 999;
1093  identifiers_[5] = 999;
1094  break;
1095  }
1096  case 6: {
1097  identifiers_[0] = tTopo->tecModule(rawid);
1098  identifiers_[1] = tTopo->tecRing(rawid);
1099  identifiers_[2] = tTopo->tecPetalInfo(rawid)[0];
1100  identifiers_[3] = tTopo->tecPetalInfo(rawid)[1];
1101  identifiers_[4] = tTopo->tecWheel(rawid);
1102  identifiers_[5] = tTopo->tecSide(rawid);
1103  break;
1104  }
1105  default: {
1106  edm::LogInfo("TrackerGeometryCompare") << "Error: bad subdetid!!";
1107  break;
1108  }
1109  }
1110 }
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:149
TrackerTopology::pxbLayer
unsigned int pxbLayer(const DetId &id) const
Definition: TrackerTopology.h:146
TrackerTopology::tidModuleInfo
std::vector< unsigned int > tidModuleInfo(const DetId &id) const
Definition: TrackerTopology.h:470
TrackerTopology::pxfBlade
unsigned int pxfBlade(const DetId &id) const
Definition: TrackerTopology.h:483
TrackerTopology::tibModule
unsigned int tibModule(const DetId &id) const
Definition: TrackerTopology.h:174
TrackerTopology::tidSide
unsigned int tidSide(const DetId &id) const
Definition: TrackerTopology.h:192
TrackerTopology::pxfModule
unsigned int pxfModule(const DetId &id) const
Definition: TrackerTopology.h:165
TrackerTopology::tecPetalInfo
std::vector< unsigned int > tecPetalInfo(const DetId &id) const
Definition: TrackerTopology.h:443
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:203
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:200
TrackerTopology::pxbLadder
unsigned int pxbLadder(const DetId &id) const
Definition: TrackerTopology.h:157
TrackerTopology::tecRing
unsigned int tecRing(const DetId &id) const
ring id
Definition: TrackerTopology.h:219
TrackerTopology::tecModule
unsigned int tecModule(const DetId &id) const
Definition: TrackerTopology.h:171
TrackerTopology::tecSide
unsigned int tecSide(const DetId &id) const
Definition: TrackerTopology.h:186
TrackerTopology::pxfDisk
unsigned int pxfDisk(const DetId &id) const
Definition: TrackerTopology.h:482
TrackerTopology::pxfPanel
unsigned int pxfPanel(const DetId &id) const
Definition: TrackerTopology.h:486
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
TrackerTopology::pxfSide
unsigned int pxfSide(const DetId &id) const
Definition: TrackerTopology.h:194
TrackerTopology::tibStringInfo
std::vector< unsigned int > tibStringInfo(const DetId &id) const
Definition: TrackerTopology.h:457
TrackerTopology::tobRodInfo
std::vector< unsigned int > tobRodInfo(const DetId &id) const
Definition: TrackerTopology.h:228
TrackerTopology::tidRing
unsigned int tidRing(const DetId &id) const
Definition: TrackerTopology.h:220
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:152
TrackerTopology::tobModule
unsigned int tobModule(const DetId &id) const
Definition: TrackerTopology.h:168
TrackerTopology::pxbModule
unsigned int pxbModule(const DetId &id) const
Definition: TrackerTopology.h:162

◆ fillTree()

void TrackerGeometryCompare::fillTree ( Alignable *  refAli,
const AlgebraicVector &  diff,
const TrackerTopology *  tTopo,
const edm::EventSetup &  iSetup 
)
private

Definition at line 847 of file TrackerGeometryCompare.cc.

References align::AlignableDet, align::AlignableDetUnit, Alignable::alignableObjectId(), alignTree_, alphaVal_, badModuleQuality_, betaVal_, dalphaVal_, daVal_, dbetaVal_, dbVal_, detDim_, detIdFlag_, dgammaVal_, dgVal_, change_name::diff, dphiVal_, drVal_, duVal_, dvVal_, dwVal_, dxVal_, dyVal_, dzVal_, PV3DBase< T, PVType, FrameType >::eta(), etaVal_, fillIdentifiers(), g, gammaVal_, Alignable::geomDetId(), edm::EventSetup::getData(), Alignable::globalPosition(), Alignable::globalRotation(), mps_fire::i, Alignable::id(), id_, inModuleList_, SiStripQuality::IsModuleBad(), SiPixelQuality::IsModuleBad(), AlignableSurface::length(), level_, mid_, mlevel_, moduleList_, Alignable::mother(), passIdCut(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), phiVal_, pixQualityToken_, DetId::rawId(), makeMuonMisalignmentScenario::rot, hcal_runs::rt, rVal_, stripQualityToken_, DetId::subdetId(), sublevel_, Alignable::surface(), surLength_, surRot_, surWidth_, align::toAngles(), useDetId_, AlignableSurface::width(), PV3DBase< T, PVType, FrameType >::x(), xVal_, PV3DBase< T, PVType, FrameType >::y(), yVal_, PV3DBase< T, PVType, FrameType >::z(), and zVal_.

Referenced by compareGeometries(), and core.AutoFillTreeProducer.AutoFillTreeProducer::process().

850  {
851  //Get bad modules
852  const SiPixelQuality* SiPixelModules = &iSetup.getData(pixQualityToken_);
853  const SiStripQuality* SiStripModules = &iSetup.getData(stripQualityToken_);
854 
855  id_ = refAli->id();
856 
857  badModuleQuality_ = 0;
858  //check if module has a bad quality tag
859  if (SiPixelModules->IsModuleBad(id_)) {
860  badModuleQuality_ = 1;
861  }
862  if (SiStripModules->IsModuleBad(id_)) {
863  badModuleQuality_ = 1;
864  }
865 
866  //check if module is in a given list of bad/untouched etc. modules
867  inModuleList_ = 0;
868  for (unsigned int i = 0; i < moduleList_.size(); i++) {
869  if (moduleList_[i] == id_) {
870  inModuleList_ = 1;
871  break;
872  }
873  }
874 
875  level_ = refAli->alignableObjectId();
876  //need if ali has no mother
877  if (refAli->mother()) {
878  mid_ = refAli->mother()->geomDetId().rawId();
879  mlevel_ = refAli->mother()->alignableObjectId();
880  } else {
881  mid_ = -1;
882  mlevel_ = -1;
883  }
884  DetId detid(id_);
885  sublevel_ = detid.subdetId();
886  fillIdentifiers(sublevel_, id_, tTopo);
887  xVal_ = refAli->globalPosition().x();
888  yVal_ = refAli->globalPosition().y();
889  zVal_ = refAli->globalPosition().z();
890  align::GlobalVector vec(xVal_, yVal_, zVal_);
891  rVal_ = vec.perp();
892  phiVal_ = vec.phi();
893  etaVal_ = vec.eta();
894  align::RotationType rot = refAli->globalRotation();
895  align::EulerAngles eulerAngles = align::toAngles(rot);
896  alphaVal_ = eulerAngles[0];
897  betaVal_ = eulerAngles[1];
898  gammaVal_ = eulerAngles[2];
899  // global
900  dxVal_ = diff[0];
901  dyVal_ = diff[1];
902  dzVal_ = diff[2];
903  // local
904  duVal_ = diff[6];
905  dvVal_ = diff[7];
906  dwVal_ = diff[8];
907  //...TODO...
908  align::GlobalVector g(dxVal_, dyVal_, dzVal_);
909  //getting dR and dPhi
910  align::GlobalVector vRef(xVal_, yVal_, zVal_);
911  align::GlobalVector vCur(xVal_ + dxVal_, yVal_ + dyVal_, zVal_ + dzVal_);
912  drVal_ = vCur.perp() - vRef.perp();
913  dphiVal_ = vCur.phi() - vRef.phi();
914  // global
915  dalphaVal_ = diff[3];
916  dbetaVal_ = diff[4];
917  dgammaVal_ = diff[5];
918  // local
919  daVal_ = diff[9];
920  dbVal_ = diff[10];
921  dgVal_ = diff[11];
922  //detIdFlag
923  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
924  if (detIdFlag_) {
925  if ((passIdCut(refAli->id())) || (passIdCut(refAli->mother()->id()))) {
926  useDetId_ = 1;
927  } else {
928  useDetId_ = 0;
929  }
930  }
931  }
932  // det module dimension
933  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
934  if (refAli->mother()->alignableObjectId() != align::AlignableDet)
935  detDim_ = 1;
936  else if (refAli->mother()->alignableObjectId() == align::AlignableDet)
937  detDim_ = 2;
938  } else
939  detDim_ = 0;
940 
941  surWidth_ = refAli->surface().width();
942  surLength_ = refAli->surface().length();
943  align::RotationType rt = refAli->globalRotation();
944  surRot_[0] = rt.xx();
945  surRot_[1] = rt.xy();
946  surRot_[2] = rt.xz();
947  surRot_[3] = rt.yx();
948  surRot_[4] = rt.yy();
949  surRot_[5] = rt.yz();
950  surRot_[6] = rt.zx();
951  surRot_[7] = rt.zy();
952  surRot_[8] = rt.zz();
953 
954  //Fill
955  alignTree_->Fill();
956 }
Alignable::mother
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:91
mps_fire.i
i
Definition: mps_fire.py:429
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:132
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
SiStripQuality::IsModuleBad
bool IsModuleBad(uint32_t detid) const
Definition: SiStripQuality.cc:557
AlignableSurface::width
align::Scalar width() const
Definition: AlignableSurface.h:30
SiPixelQuality::IsModuleBad
bool IsModuleBad(const uint32_t &detid) const
Definition: SiPixelQuality.cc:93
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:130
g
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
TrackerGeometryCompare::stripQualityToken_
const edm::ESGetToken< SiStripQuality, SiStripQualityRcd > stripQualityToken_
Definition: TrackerGeometryCompare.h:95
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
Alignable::globalPosition
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:135
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
hcal_runs.rt
rt
Definition: hcal_runs.py:76
change_name.diff
diff
Definition: change_name.py:13
align::toAngles
EulerAngles toAngles(const RotationType &)
Convert rotation matrix to angles about x-, y-, z-axes (frame rotation).
Definition: Utilities.cc:8
DetId
Definition: DetId.h:17
Alignable::id
align::ID id() const
Return the ID of Alignable, i.e. DetId of &#39;first&#39; component GeomDet(Unit).
Definition: Alignable.h:180
align::EulerAngles
AlgebraicVector EulerAngles
Definition: Definitions.h:34
Alignable::geomDetId
const DetId & geomDetId() const
Definition: Alignable.h:177
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
AlignableSurface::length
align::Scalar length() const
Definition: AlignableSurface.h:32
TrackerGeometryCompare::fillIdentifiers
void fillIdentifiers(int subdetlevel, int rawid, const TrackerTopology *tTopo)
Definition: TrackerGeometryCompare.cc:1049
Alignable::globalRotation
const RotationType & globalRotation() const
Return the global orientation of the object.
Definition: Alignable.h:138
TrackerGeometryCompare::pixQualityToken_
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
Definition: TrackerGeometryCompare.h:94

◆ passIdCut()

bool TrackerGeometryCompare::passIdCut ( uint32_t  id)
private

Definition at line 1037 of file TrackerGeometryCompare.cc.

References detIdFlagVector_, and mps_fire::i.

Referenced by fillTree().

1037  {
1038  bool pass = false;
1039  int nEntries = detIdFlagVector_.size();
1040 
1041  for (int i = 0; i < nEntries; i++) {
1042  if (detIdFlagVector_[i] == id)
1043  pass = true;
1044  }
1045 
1046  return pass;
1047 }
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometryCompare::detIdFlagVector_
std::vector< uint32_t > detIdFlagVector_
Definition: TrackerGeometryCompare.h:123

◆ setCommonTrackerSystem()

void TrackerGeometryCompare::setCommonTrackerSystem ( )
private

Definition at line 770 of file TrackerGeometryCompare.cc.

References commonTrackerLevel_, currentTracker, diffCommonTrackerSystem(), Alignable::globalPosition(), align::moveAlignable(), AlignableTracker::objectIdProvider(), referenceTracker, makeMuonMisalignmentScenario::rot, setCommonTrackerSystem_, AlignableObjectId::stringToId(), align::toMatrix(), TrackerCommonCM_, TrackerCommonR_, TrackerCommonT_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

770  {
771  edm::LogInfo("TrackerGeometryCompare") << "Setting Common Tracker System....";
772 
773  // DM_534??AlignableObjectId dummy;
774  // DM_534??_commonTrackerLevel = dummy.nameToType(_setCommonTrackerSystem);
775  commonTrackerLevel_ = currentTracker->objectIdProvider().stringToId(setCommonTrackerSystem_);
776 
777  diffCommonTrackerSystem(referenceTracker, currentTracker);
778 
779  align::EulerAngles dOmega(3);
780  dOmega[0] = TrackerCommonR_.x();
781  dOmega[1] = TrackerCommonR_.y();
782  dOmega[2] = TrackerCommonR_.z();
783  align::RotationType rot = align::toMatrix(dOmega);
784  align::GlobalVector theR = TrackerCommonT_;
785 
786  edm::LogInfo("TrackerGeometryCompare") << "what we get from overlaying the pixels..." << theR << ", " << rot;
787 
788  //transform to the Tracker System
789  align::PositionType trackerCM = currentTracker->globalPosition();
790  align::GlobalVector cmDiff(
791  trackerCM.x() - TrackerCommonCM_.x(), trackerCM.y() - TrackerCommonCM_.y(), trackerCM.z() - TrackerCommonCM_.z());
792 
793  edm::LogInfo("TrackerGeometryCompare") << "Pixel CM: " << TrackerCommonCM_ << ", tracker CM: " << trackerCM;
794 
795  //adjust translational difference factoring in different rotational CM
796  //needed because rotateInGlobalFrame is about CM of alignable, not Tracker
797  const align::GlobalVector::BasicVectorType& lpvgf = cmDiff.basicVector();
798  align::GlobalVector moveV(rot.multiplyInverse(lpvgf) - lpvgf);
799  align::GlobalVector theRprime(theR + moveV);
800 
801  AlgebraicVector TrackerCommonTR(6);
802  TrackerCommonTR(1) = theRprime.x();
803  TrackerCommonTR(2) = theRprime.y();
804  TrackerCommonTR(3) = theRprime.z();
805  TrackerCommonTR(4) = TrackerCommonR_.x();
806  TrackerCommonTR(5) = TrackerCommonR_.y();
807  TrackerCommonTR(6) = TrackerCommonR_.z();
808 
809  edm::LogInfo("TrackerGeometryCompare") << "and after the transformation: " << TrackerCommonTR;
810  align::moveAlignable(currentTracker, TrackerCommonTR);
811 }
TrackerGeometryCompare::TrackerCommonR_
align::GlobalVector TrackerCommonR_
Definition: TrackerGeometryCompare.h:126
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
TrackerGeometryCompare::diffCommonTrackerSystem
void diffCommonTrackerSystem(Alignable *refAli, Alignable *curAli)
Definition: TrackerGeometryCompare.cc:813
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
Alignable::globalPosition
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:135
TrackerGeometryCompare::TrackerCommonT_
align::GlobalVector TrackerCommonT_
Definition: TrackerGeometryCompare.h:125
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:13
align::EulerAngles
AlgebraicVector EulerAngles
Definition: Definitions.h:34
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:99
TrackerGeometryCompare::commonTrackerLevel_
align::StructureType commonTrackerLevel_
Definition: TrackerGeometryCompare.h:124
AlignableObjectId::stringToId
align::StructureType stringToId(const char *) const
Definition: AlignableObjectId.cc:284
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:97
TrackerGeometryCompare::TrackerCommonCM_
align::PositionType TrackerCommonCM_
Definition: TrackerGeometryCompare.h:127
align::toMatrix
RotationType toMatrix(const EulerAngles &)
Convert rotation angles about x-, y-, z-axes to matrix.
Definition: Utilities.cc:34
AlignableTracker::objectIdProvider
const AlignableObjectId & objectIdProvider() const
Return tracker alignable object ID provider derived from the tracker&#39;s geometry.
Definition: AlignableTracker.h:110
align::moveAlignable
void moveAlignable(Alignable *ali, AlgebraicVector diff)
Moves the alignable by the AlgebraicVector.
Definition: AlignTools.cc:84

◆ surveyToTracker()

void TrackerGeometryCompare::surveyToTracker ( AlignableTracker *  ali,
Alignments *  alignVals,
AlignmentErrorsExtended *  alignErrors 
)
private

Definition at line 958 of file TrackerGeometryCompare.cc.

References AlCaHLTBitMon_QueryRunRegistry::comp, filterCSVwithJSON::copy, AlignableTracker::endcapGeomDets(), mps_fire::i, AlignableTracker::innerBarrelGeomDets(), dqmiolumiharvest::j, dqmdumpme::k, Alignments::m_align, AlignmentErrorsExtended::m_alignError, AlignableTracker::outerBarrelGeomDets(), AlignableTracker::pixelEndcapGeomDets(), AlignableTracker::pixelHalfBarrelGeomDets(), SurveyDet::position(), makeMuonMisalignmentScenario::rot, SurveyDet::rotation(), jetUpdater_cfi::sort, AlignableTracker::TIDGeomDets(), and HcalDetIdTransform::transform().

960  {
961  //getting the right alignables for the alignment record
962  auto detPB = ali->pixelHalfBarrelGeomDets();
963  auto detPEC = ali->pixelEndcapGeomDets();
964  auto detTIB = ali->innerBarrelGeomDets();
965  auto detTID = ali->TIDGeomDets();
966  auto detTOB = ali->outerBarrelGeomDets();
967  auto detTEC = ali->endcapGeomDets();
968 
969  align::Alignables allGeomDets;
970  std::copy(detPB.begin(), detPB.end(), std::back_inserter(allGeomDets));
971  std::copy(detPEC.begin(), detPEC.end(), std::back_inserter(allGeomDets));
972  std::copy(detTIB.begin(), detTIB.end(), std::back_inserter(allGeomDets));
973  std::copy(detTID.begin(), detTID.end(), std::back_inserter(allGeomDets));
974  std::copy(detTOB.begin(), detTOB.end(), std::back_inserter(allGeomDets));
975  std::copy(detTEC.begin(), detTEC.end(), std::back_inserter(allGeomDets));
976 
977  align::Alignables rcdAlis;
978  for (const auto& i : allGeomDets) {
979  if (i->components().size() == 1) {
980  rcdAlis.push_back(i);
981  } else if (i->components().size() > 1) {
982  rcdAlis.push_back(i);
983  const auto& comp = i->components();
984  for (const auto& j : comp)
985  rcdAlis.push_back(j);
986  }
987  }
988 
989  //turning them into alignments
990  for (const auto& k : rcdAlis) {
991  const SurveyDet* surveyInfo = k->survey();
992  const align::PositionType& pos(surveyInfo->position());
993  align::RotationType rot(surveyInfo->rotation());
994  CLHEP::Hep3Vector clhepVector(pos.x(), pos.y(), pos.z());
995  CLHEP::HepRotation clhepRotation(
996  CLHEP::HepRep3x3(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
997  AlignTransform transform(clhepVector, clhepRotation, k->id());
998  AlignTransformErrorExtended transformError(CLHEP::HepSymMatrix(3, 1), k->id());
999  alignVals->m_align.push_back(transform);
1000  alignErrors->m_alignError.push_back(transformError);
1001  }
1002 
1003  // to get the right order, sort by rawId
1004  std::sort(alignVals->m_align.begin(), alignVals->m_align.end());
1005  std::sort(alignErrors->m_alignError.begin(), alignErrors->m_alignError.end());
1006 }
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:30
AlignableTracker::pixelHalfBarrelGeomDets
Alignables & pixelHalfBarrelGeomDets()
Return pixel barrel GeomDets.
Definition: AlignableTracker.h:59
mps_fire.i
i
Definition: mps_fire.py:429
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:249
Alignments::m_align
std::vector< AlignTransform > m_align
Definition: Alignments.h:19
SurveyDet::rotation
const align::RotationType & rotation() const
Definition: SurveyDet.h:60
AlignableTracker::innerBarrelGeomDets
Alignables & innerBarrelGeomDets()
Return inner barrel GeomDets.
Definition: AlignableTracker.h:55
AlignableTracker::TIDGeomDets
Alignables & TIDGeomDets()
Return TID GeomDets.
Definition: AlignableTracker.h:63
SurveyDet::position
const align::PositionType & position() const
Definition: SurveyDet.h:58
AlignableTracker::endcapGeomDets
Alignables & endcapGeomDets()
Return endcap GeomDets.
Definition: AlignableTracker.h:61
AlignableTracker::pixelEndcapGeomDets
Alignables & pixelEndcapGeomDets()
Return pixel endcap GeomDets.
Definition: AlignableTracker.h:65
AlignableTracker::outerBarrelGeomDets
Alignables & outerBarrelGeomDets()
Return outer barrel GeomDets.
Definition: AlignableTracker.h:57
AlignmentErrorsExtended::m_alignError
std::vector< AlignTransformErrorExtended > m_alignError
Definition: AlignmentErrorsExtended.h:19
align::Alignables
std::vector< Alignable * > Alignables
Definition: Utilities.h:31
HcalDetIdTransform::transform
unsigned transform(const HcalDetId &id, unsigned transformCode)
Definition: HcalDetIdTransform.cc:7

Member Data Documentation

◆ alignTree_

TTree* TrackerGeometryCompare::alignTree_
private

Definition at line 137 of file TrackerGeometryCompare.h.

Referenced by endJob(), fillTree(), and TrackerGeometryCompare().

◆ alphaVal_

float TrackerGeometryCompare::alphaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ badModuleQuality_

int TrackerGeometryCompare::badModuleQuality_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ betaVal_

float TrackerGeometryCompare::betaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ commonTrackerLevel_

align::StructureType TrackerGeometryCompare::commonTrackerLevel_
private

Definition at line 124 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ cpvTokenDD4hep_

const edm::ESGetToken<cms::DDCompactView, IdealGeometryRecord> TrackerGeometryCompare::cpvTokenDD4hep_
private

Definition at line 90 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ cpvTokenDDD_

const edm::ESGetToken<DDCompactView, IdealGeometryRecord> TrackerGeometryCompare::cpvTokenDDD_
private

Definition at line 89 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ currentTracker

AlignableTracker* TrackerGeometryCompare::currentTracker
private

Definition at line 99 of file TrackerGeometryCompare.h.

Referenced by analyze(), createROOTGeometry(), and setCommonTrackerSystem().

◆ dalphaVal_

float TrackerGeometryCompare::dalphaVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ daVal_

float TrackerGeometryCompare::daVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dbetaVal_

float TrackerGeometryCompare::dbetaVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dbVal_

float TrackerGeometryCompare::dbVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detDim_

int TrackerGeometryCompare::detDim_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detIdFlag_

bool TrackerGeometryCompare::detIdFlag_
private

Definition at line 116 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detIdFlagFile_

std::string TrackerGeometryCompare::detIdFlagFile_
private

Definition at line 117 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ detIdFlagVector_

std::vector<uint32_t> TrackerGeometryCompare::detIdFlagVector_
private

Definition at line 123 of file TrackerGeometryCompare.h.

Referenced by passIdCut(), and TrackerGeometryCompare().

◆ dgammaVal_

float TrackerGeometryCompare::dgammaVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dgVal_

float TrackerGeometryCompare::dgVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dphiVal_

float TrackerGeometryCompare::dphiVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ drVal_

float TrackerGeometryCompare::drVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dummyTracker

AlignableTracker* TrackerGeometryCompare::dummyTracker
private

Definition at line 98 of file TrackerGeometryCompare.h.

◆ duVal_

float TrackerGeometryCompare::duVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dvVal_

float TrackerGeometryCompare::dvVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dwVal_

float TrackerGeometryCompare::dwVal_
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dxVal_

float TrackerGeometryCompare::dxVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dyVal_

float TrackerGeometryCompare::dyVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dzVal_

float TrackerGeometryCompare::dzVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ etaVal_

float TrackerGeometryCompare::etaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ filename_

std::string TrackerGeometryCompare::filename_
private

Definition at line 134 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ firstEvent_

bool TrackerGeometryCompare::firstEvent_
private

Definition at line 162 of file TrackerGeometryCompare.h.

Referenced by analyze(), and beginJob().

◆ fromDD4hep_

bool TrackerGeometryCompare::fromDD4hep_
private

Definition at line 112 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ gammaVal_

float TrackerGeometryCompare::gammaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ geomDetToken_

const edm::ESGetToken<GeometricDet, IdealGeometryRecord> TrackerGeometryCompare::geomDetToken_
private

Definition at line 92 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ id_

int TrackerGeometryCompare::id_
private

\ Tree variables

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ identifiers_

uint32_t TrackerGeometryCompare::identifiers_[6]
private

Definition at line 153 of file TrackerGeometryCompare.h.

Referenced by fillIdentifiers(), and TrackerGeometryCompare().

◆ inModuleList_

int TrackerGeometryCompare::inModuleList_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ inputFilename1_

std::string TrackerGeometryCompare::inputFilename1_
private

Definition at line 108 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), createROOTGeometry(), and TrackerGeometryCompare().

◆ inputFilename2_

std::string TrackerGeometryCompare::inputFilename2_
private

Definition at line 109 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), createROOTGeometry(), and TrackerGeometryCompare().

◆ inputRootFile1_

TFile* TrackerGeometryCompare::inputRootFile1_
private

Definition at line 138 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

◆ inputRootFile2_

TFile* TrackerGeometryCompare::inputRootFile2_
private

Definition at line 139 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

◆ inputTree01_

TTree* TrackerGeometryCompare::inputTree01_
private

Definition at line 140 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ inputTree02_

TTree* TrackerGeometryCompare::inputTree02_
private

Definition at line 141 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ inputTree11_

TTree* TrackerGeometryCompare::inputTree11_
private

Definition at line 142 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ inputTree12_

TTree* TrackerGeometryCompare::inputTree12_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ inputTreenameAlign_

std::string TrackerGeometryCompare::inputTreenameAlign_
private

Definition at line 110 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

◆ inputTreenameDeform_

std::string TrackerGeometryCompare::inputTreenameDeform_
private

Definition at line 111 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), createROOTGeometry(), and TrackerGeometryCompare().

◆ level_

int TrackerGeometryCompare::level_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ levelStrings_

const std::vector<std::string> TrackerGeometryCompare::levelStrings_
private

Definition at line 106 of file TrackerGeometryCompare.h.

Referenced by analyze().

◆ m_h1_

std::map<std::string, TH1D*> TrackerGeometryCompare::m_h1_
private

Definition at line 166 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_nBins_

int TrackerGeometryCompare::m_nBins_
private

Definition at line 158 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_rangeHigh_

double TrackerGeometryCompare::m_rangeHigh_
private

Definition at line 160 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_rangeLow_

double TrackerGeometryCompare::m_rangeLow_
private

Definition at line 159 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_theLevels

std::vector<align::StructureType> TrackerGeometryCompare::m_theLevels
private

Definition at line 60 of file TrackerGeometryCompare.h.

Referenced by analyze(), and compareGeometries().

◆ m_vtkmap_

std::vector<TrackerMap> TrackerGeometryCompare::m_vtkmap_
private

Definition at line 164 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), endJob(), and TrackerGeometryCompare().

◆ mid_

int TrackerGeometryCompare::mid_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ mlevel_

int TrackerGeometryCompare::mlevel_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ moduleInList_

int TrackerGeometryCompare::moduleInList_
private

Definition at line 131 of file TrackerGeometryCompare.h.

◆ moduleList_

std::vector<int> TrackerGeometryCompare::moduleList_
private

Definition at line 130 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and fillTree().

◆ moduleListFile_

std::ifstream TrackerGeometryCompare::moduleListFile_
private

Definition at line 129 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ moduleListName_

std::string TrackerGeometryCompare::moduleListName_
private

Definition at line 107 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

◆ phase_

SiPixelPI::phase TrackerGeometryCompare::phase_
private

Definition at line 121 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

◆ phiVal_

float TrackerGeometryCompare::phiVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ pixQualityToken_

const edm::ESGetToken<SiPixelQuality, SiPixelQualityRcd> TrackerGeometryCompare::pixQualityToken_
private

Definition at line 94 of file TrackerGeometryCompare.h.

Referenced by fillTree().

◆ ptpToken_

const edm::ESGetToken<PTrackerParameters, PTrackerParametersRcd> TrackerGeometryCompare::ptpToken_
private

Definition at line 93 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ referenceTracker

AlignableTracker* TrackerGeometryCompare::referenceTracker
private

Definition at line 97 of file TrackerGeometryCompare.h.

Referenced by analyze(), createROOTGeometry(), and setCommonTrackerSystem().

◆ rVal_

float TrackerGeometryCompare::rVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ setCommonTrackerSystem_

std::string TrackerGeometryCompare::setCommonTrackerSystem_
private

Definition at line 115 of file TrackerGeometryCompare.h.

Referenced by analyze(), setCommonTrackerSystem(), and TrackerGeometryCompare().

◆ stripQualityToken_

const edm::ESGetToken<SiStripQuality, SiStripQualityRcd> TrackerGeometryCompare::stripQualityToken_
private

Definition at line 95 of file TrackerGeometryCompare.h.

Referenced by fillTree().

◆ sublevel_

int TrackerGeometryCompare::sublevel_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surfDeform_

double TrackerGeometryCompare::surfDeform_[13]
private

Definition at line 156 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ surfdir_

std::string TrackerGeometryCompare::surfdir_
private

Definition at line 135 of file TrackerGeometryCompare.h.

Referenced by endJob(), and TrackerGeometryCompare().

◆ surLength_

float TrackerGeometryCompare::surLength_
private

Definition at line 152 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surRot_

double TrackerGeometryCompare::surRot_[9]
private

Definition at line 154 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surWidth_

float TrackerGeometryCompare::surWidth_
private

Definition at line 152 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ theFile_

TFile* TrackerGeometryCompare::theFile_
private

Definition at line 136 of file TrackerGeometryCompare.h.

Referenced by endJob(), and TrackerGeometryCompare().

◆ theSurveyErrors

const SurveyErrors* TrackerGeometryCompare::theSurveyErrors
private

Definition at line 103 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ theSurveyIndex

unsigned int TrackerGeometryCompare::theSurveyIndex
private

Definition at line 101 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ theSurveyValues

const Alignments* TrackerGeometryCompare::theSurveyValues
private

Definition at line 102 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TrackerGeometryCompare::topoToken_
private

Definition at line 91 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ TrackerCommonCM_

align::PositionType TrackerGeometryCompare::TrackerCommonCM_
private

Definition at line 127 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ TrackerCommonR_

align::GlobalVector TrackerGeometryCompare::TrackerCommonR_
private

Definition at line 126 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ TrackerCommonT_

align::GlobalVector TrackerGeometryCompare::TrackerCommonT_
private

Definition at line 125 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ type_

int TrackerGeometryCompare::type_
private

Definition at line 155 of file TrackerGeometryCompare.h.

Referenced by Modules.Service::__str__(), Modules.Service::_placeImpl(), Modules.ESSource::_placeImpl(), Modules.ESProducer::_placeImpl(), Modules.ESPrefer::copy(), Mixins._TypedParameterizable::dumpPython(), Modules.ESPrefer::dumpPythonAs(), Modules.Service::dumpSequencePython(), RandomRunSource.RandomRunSource::insertInto(), Modules.Service::insertInto(), ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::moduleLabel_(), Modules.ESProducer::moduleLabel_(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), Modules.ESPrefer::nameInProcessDesc_(), and TrackerGeometryCompare().

◆ useDetId_

int TrackerGeometryCompare::useDetId_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ weightBy_

std::string TrackerGeometryCompare::weightBy_
private

Definition at line 114 of file TrackerGeometryCompare.h.

Referenced by compareGeometries(), diffCommonTrackerSystem(), and TrackerGeometryCompare().

◆ weightById_

bool TrackerGeometryCompare::weightById_
private

Definition at line 118 of file TrackerGeometryCompare.h.

Referenced by compareGeometries(), diffCommonTrackerSystem(), and TrackerGeometryCompare().

◆ weightByIdFile_

std::string TrackerGeometryCompare::weightByIdFile_
private

Definition at line 119 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ weightByIdVector_

std::vector<unsigned int> TrackerGeometryCompare::weightByIdVector_
private

Definition at line 120 of file TrackerGeometryCompare.h.

Referenced by compareGeometries(), diffCommonTrackerSystem(), and TrackerGeometryCompare().

◆ writeToDB_

bool TrackerGeometryCompare::writeToDB_
private

Definition at line 113 of file TrackerGeometryCompare.h.

Referenced by analyze().

◆ xVal_

float TrackerGeometryCompare::xVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ yVal_

float TrackerGeometryCompare::yVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ zVal_

float TrackerGeometryCompare::zVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

Generated for CMSSW Reference Manual by doxygen 1.8.14