CMS 3D CMS Logo
List of all members | Public Types | 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
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > 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
 
bool registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

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_
 
float phiVal_
 
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
 
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
 
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

- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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 ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
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 55 of file TrackerGeometryCompare.cc.

References alignTree_, alphaVal_, badModuleQuality_, betaVal_, looper::cfg, dqmiodatasetharvest::ctr, 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_, 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_.

56  : cpvTokenDDD_(esConsumes()),
57  cpvTokenDD4hep_(esConsumes()),
58  topoToken_(esConsumes()),
59  geomDetToken_(esConsumes()),
60  ptpToken_(esConsumes()),
61  ptitpToken_(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_(nullptr),
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  int ctr = 0;
125  while (!inFile.eof()) {
126  ctr++;
127  unsigned int listId;
128  inFile >> listId;
129  inFile.ignore(256, '\n');
130 
131  weightByIdVector_.push_back(listId);
132  }
133  inFile.close();
134  }
135 
136  //root configuration
137  theFile_ = new TFile(filename_.c_str(), "RECREATE");
138  alignTree_ = new TTree("alignTree",
139  "alignTree"); //,"id:level:mid:mlevel:sublevel:x:y:z:r:phi:a:b:c:dx:dy:dz:dr:dphi:da:db:dc");
140  alignTree_->Branch("id", &id_, "id/I");
141  alignTree_->Branch("badModuleQuality", &badModuleQuality_, "badModuleQuality/I");
142  alignTree_->Branch("inModuleList", &inModuleList_, "inModuleList/I");
143  alignTree_->Branch("level", &level_, "level/I");
144  alignTree_->Branch("mid", &mid_, "mid/I");
145  alignTree_->Branch("mlevel", &mlevel_, "mlevel/I");
146  alignTree_->Branch("sublevel", &sublevel_, "sublevel/I");
147  alignTree_->Branch("x", &xVal_, "x/F");
148  alignTree_->Branch("y", &yVal_, "y/F");
149  alignTree_->Branch("z", &zVal_, "z/F");
150  alignTree_->Branch("r", &rVal_, "r/F");
151  alignTree_->Branch("phi", &phiVal_, "phi/F");
152  alignTree_->Branch("eta", &etaVal_, "eta/F");
153  alignTree_->Branch("alpha", &alphaVal_, "alpha/F");
154  alignTree_->Branch("beta", &betaVal_, "beta/F");
155  alignTree_->Branch("gamma", &gammaVal_, "gamma/F");
156  alignTree_->Branch("dx", &dxVal_, "dx/F");
157  alignTree_->Branch("dy", &dyVal_, "dy/F");
158  alignTree_->Branch("dz", &dzVal_, "dz/F");
159  alignTree_->Branch("dr", &drVal_, "dr/F");
160  alignTree_->Branch("dphi", &dphiVal_, "dphi/F");
161  alignTree_->Branch("dalpha", &dalphaVal_, "dalpha/F");
162  alignTree_->Branch("dbeta", &dbetaVal_, "dbeta/F");
163  alignTree_->Branch("dgamma", &dgammaVal_, "dgamma/F");
164  alignTree_->Branch("du", &duVal_, "du/F");
165  alignTree_->Branch("dv", &dvVal_, "dv/F");
166  alignTree_->Branch("dw", &dwVal_, "dw/F");
167  alignTree_->Branch("da", &daVal_, "da/F");
168  alignTree_->Branch("db", &dbVal_, "db/F");
169  alignTree_->Branch("dg", &dgVal_, "dg/F");
170  alignTree_->Branch("useDetId", &useDetId_, "useDetId/I");
171  alignTree_->Branch("detDim", &detDim_, "detDim/I");
172  alignTree_->Branch("surW", &surWidth_, "surW/F");
173  alignTree_->Branch("surL", &surLength_, "surL/F");
174  alignTree_->Branch("surRot", &surRot_, "surRot[9]/D");
175  alignTree_->Branch("identifiers", &identifiers_, "identifiers[6]/I");
176  alignTree_->Branch("type", &type_, "type/I");
177  alignTree_->Branch("surfDeform", &surfDeform_, "surfDeform[13]/D");
178 
179  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
180  it->setPalette(1);
181  it->addPixel(true);
182  }
183 
184  edm::Service<TFileService> fs;
185  TFileDirectory subDir_All = fs->mkdir("AllSubdetectors");
186  TFileDirectory subDir_PXB = fs->mkdir("PixelBarrel");
187  TFileDirectory subDir_PXF = fs->mkdir("PixelEndcap");
188  for (int ii = 0; ii < 13; ++ii) {
189  std::stringstream histname0;
190  histname0 << "SurfDeform_Par_" << ii;
191  m_h1_[histname0.str()] = subDir_All.make<TH1D>(
192  (histname0.str()).c_str(), (histname0.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
193 
194  std::stringstream histname1;
195  histname1 << "SurfDeform_PixelBarrel_Par_" << ii;
196  m_h1_[histname1.str()] = subDir_PXB.make<TH1D>(
197  (histname1.str()).c_str(), (histname1.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
198 
199  std::stringstream histname2;
200  histname2 << "SurfDeform_PixelEndcap_Par_" << ii;
201  m_h1_[histname2.str()] = subDir_PXF.make<TH1D>(
202  (histname2.str()).c_str(), (histname2.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
203  }
204 }
TrackerGeometryCompare::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: TrackerGeometryCompare.h:91
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:89
TrackerGeometryCompare::levelStrings_
const std::vector< std::string > levelStrings_
Definition: TrackerGeometryCompare.h:105
TrackerGeometryCompare::theSurveyValues
const Alignments * theSurveyValues
Definition: TrackerGeometryCompare.h:101
TrackerGeometryCompare::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: TrackerGeometryCompare.h:92
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:162
TrackerGeometryCompare::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: TrackerGeometryCompare.h:90
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:65
TrackerGeometryCompare::dummyTracker
AlignableTracker * dummyTracker
Definition: TrackerGeometryCompare.h:97
TrackerGeometryCompare::m_h1_
std::map< std::string, TH1D * > m_h1_
Definition: TrackerGeometryCompare.h:164
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:128
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TrackerGeometryCompare::stripQualityToken_
const edm::ESGetToken< SiStripQuality, SiStripQualityRcd > stripQualityToken_
Definition: TrackerGeometryCompare.h:94
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:121
TrackerGeometryCompare::cpvTokenDDD_
const edm::ESGetToken< DDCompactView, IdealGeometryRecord > cpvTokenDDD_
Definition: TrackerGeometryCompare.h:87
TrackerGeometryCompare::weightByIdVector_
std::vector< unsigned int > weightByIdVector_
Definition: TrackerGeometryCompare.h:119
testHGCalDigi_cfg.inFile
inFile
Definition: testHGCalDigi_cfg.py:51
cuy.ii
ii
Definition: cuy.py:589
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:98
TrackerGeometryCompare::cpvTokenDD4hep_
const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord > cpvTokenDD4hep_
Definition: TrackerGeometryCompare.h:88
TrackerGeometryCompare::commonTrackerLevel_
align::StructureType commonTrackerLevel_
Definition: TrackerGeometryCompare.h:122
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:96
TrackerGeometryCompare::theSurveyErrors
const SurveyErrors * theSurveyErrors
Definition: TrackerGeometryCompare.h:102
TrackerGeometryCompare::pixQualityToken_
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
Definition: TrackerGeometryCompare.h:93

Member Function Documentation

◆ addSurveyInfo()

void TrackerGeometryCompare::addSurveyInfo ( Alignable *  ali)
private

Definition at line 958 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().

958  {
959  const auto& comp = ali->components();
960 
961  unsigned int nComp = comp.size();
962 
963  for (unsigned int i = 0; i < nComp; ++i)
964  addSurveyInfo(comp[i]);
965 
966  const SurveyError& error = theSurveyErrors->m_surveyErrors[theSurveyIndex];
967 
968  if (ali->geomDetId().rawId() != error.rawId() || ali->alignableObjectId() != error.structureType()) {
969  throw cms::Exception("DatabaseError") << "Error reading survey info from DB. Mismatched id!";
970  }
971 
972  const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[theSurveyIndex].translation();
973  const CLHEP::HepRotation& rot = theSurveyValues->m_align[theSurveyIndex].rotation();
974 
975  AlignableSurface surf(
976  align::PositionType(pos.x(), pos.y(), pos.z()),
977  align::RotationType(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
978 
979  surf.setWidth(ali->surface().width());
980  surf.setLength(ali->surface().length());
981 
982  ali->setSurvey(new SurveyDet(surf, error.matrix()));
983 
984  ++theSurveyIndex;
985 }
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:101
TrackerGeometryCompare::addSurveyInfo
void addSurveyInfo(Alignable *ali)
Definition: TrackerGeometryCompare.cc:958
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:102
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 226 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_.

226  {
227  if (firstEvent_) {
228  //Retrieve tracker topology from geometry
229  const TrackerTopology* const tTopo = &iSetup.getData(topoToken_);
230 
231  //upload the ROOT geometries
232  createROOTGeometry(iSetup);
233 
234  //setting the levels being used in the geometry comparator
235  edm::LogInfo("TrackerGeometryCompare") << "levels: " << levelStrings_.size();
236  for (const auto& level : levelStrings_) {
237  m_theLevels.push_back(currentTracker->objectIdProvider().stringToId(level));
238  edm::LogInfo("TrackerGeometryCompare") << "level: " << level;
239  edm::LogInfo("TrackerGeometryCompare")
240  << "structure type: " << currentTracker->objectIdProvider().stringToId(level);
241  }
242 
243  //set common tracker system first
244  // if setting the tracker common system
245  if (setCommonTrackerSystem_ != "NONE") {
246  setCommonTrackerSystem();
247  }
248 
249  //compare the goemetries
250  compareGeometries(referenceTracker, currentTracker, tTopo, iSetup);
251  compareSurfaceDeformations(inputTree11_, inputTree12_);
252 
253  //write out ntuple
254  //might be better to do within output module
255 
256  if (writeToDB_) {
257  Alignments* myAlignments = currentTracker->alignments();
258  AlignmentErrorsExtended* myAlignmentErrorsExtended = currentTracker->alignmentErrors();
259 
260  // 2. Store alignment[Error]s to DB
261  edm::Service<cond::service::PoolDBOutputService> poolDbService;
262  // Call service
263  if (!poolDbService.isAvailable()) // Die if not available
264  throw cms::Exception("NotAvailable") << "PoolDBOutputService not available";
265 
266  poolDbService->writeOneIOV<Alignments>(*myAlignments, poolDbService->beginOfTime(), "TrackerAlignmentRcd");
267  poolDbService->writeOneIOV<AlignmentErrorsExtended>(
268  *myAlignmentErrorsExtended, poolDbService->beginOfTime(), "TrackerAlignmentErrorExtendedRcd");
269  }
270 
271  firstEvent_ = false;
272  }
273 }
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:605
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:89
TrackerGeometryCompare::m_theLevels
std::vector< align::StructureType > m_theLevels
Definition: TrackerGeometryCompare.h:58
TrackerGeometryCompare::levelStrings_
const std::vector< std::string > levelStrings_
Definition: TrackerGeometryCompare.h:105
cond::service::PoolDBOutputService::writeOneIOV
Hash writeOneIOV(const T &payload, Time_t time, const std::string &recordName)
Definition: PoolDBOutputService.h:62
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:98
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:96
TrackerGeometryCompare::createROOTGeometry
void createROOTGeometry(const edm::EventSetup &iSetup)
Definition: TrackerGeometryCompare.cc:275
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:453

◆ beginJob()

void TrackerGeometryCompare::beginJob ( )
overridevirtual

Read from DB and print survey info.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 206 of file TrackerGeometryCompare.cc.

References firstEvent_.

206 { firstEvent_ = true; }

◆ compareGeometries()

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

Definition at line 605 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().

608  {
609  using namespace align;
610 
611  const auto& refComp = refAli->components();
612  const auto& curComp = curAli->components();
613 
614  unsigned int nComp = refComp.size();
615  //only perform for designate levels
616  bool useLevel = false;
617  for (unsigned int i = 0; i < m_theLevels.size(); ++i) {
618  if (refAli->alignableObjectId() == m_theLevels[i])
619  useLevel = true;
620  }
621 
622  //another added level for difference between det and detunit
623  //if ((refAli->alignableObjectId()==2)&&(nComp == 1)) useLevel = false;
624 
625  //coordinate matching, etc etc
626  if (useLevel) {
627  DetId detid(refAli->id());
628 
629  CLHEP::Hep3Vector Rtotal, Wtotal, lRtotal, lWtotal;
630  Rtotal.set(0., 0., 0.);
631  Wtotal.set(0., 0., 0.);
632  lRtotal.set(0., 0., 0.);
633  lWtotal.set(0., 0., 0.);
634 
635  bool converged = false;
636 
637  AlgebraicVector diff, check;
638 
639  for (int i = 0; i < 100; i++) {
640  // Get differences between alignments for rotations and translations
641  // both local and global
642  diff = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
643 
644  // 'diffAlignables' returns 'refAli - curAli' for translations and 'curAli - refAli' for rotations.
645  // The plan is to unify this at some point, but a simple change of the sign for one of them was postponed
646  // to do some further checks to understand the rotations better
647  //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.
648  CLHEP::Hep3Vector dR(-diff[0], -diff[1], -diff[2]);
649  CLHEP::Hep3Vector dW(diff[3], diff[4], diff[5]);
650  CLHEP::Hep3Vector dRLocal(-diff[6], -diff[7], -diff[8]);
651  CLHEP::Hep3Vector dWLocal(diff[9], diff[10], diff[11]);
652 
653  // Translations
654  Rtotal += dR;
655  lRtotal += dRLocal;
656 
657  //Rotations
658  CLHEP::HepRotation rot(Wtotal.unit(), Wtotal.mag());
659  CLHEP::HepRotation drot(dW.unit(), dW.mag());
660  rot *= drot;
661  Wtotal.set(rot.axis().x() * rot.delta(), rot.axis().y() * rot.delta(), rot.axis().z() * rot.delta());
662 
663  CLHEP::HepRotation rotLocal(lWtotal.unit(), lWtotal.mag());
664  CLHEP::HepRotation drotLocal(dWLocal.unit(), dWLocal.mag());
665  rotLocal *= drotLocal;
666  lWtotal.set(rotLocal.axis().x() * rotLocal.delta(),
667  rotLocal.axis().y() * rotLocal.delta(),
668  rotLocal.axis().z() * rotLocal.delta());
669 
670  // Move current alignable by shift and check if difference
671  // is smaller than tolerance value
672  // if true, break the loop
673  align::moveAlignable(curAli, diff);
674  float tolerance = 1e-7;
675  check = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
676  align::GlobalVector checkR(check[0], check[1], check[2]);
677  align::GlobalVector checkW(check[3], check[4], check[5]);
678  if ((checkR.mag() < tolerance) && (checkW.mag() < tolerance)) {
679  converged = true;
680  break;
681  }
682  }
683 
684  // give an exception if difference has not fallen below tolerance level
685  // i.e. method has not converged
686  if (!converged) {
687  edm::LogInfo("TrackerGeometryCompare")
688  << "Tolerance Exceeded!(alObjId: " << refAli->alignableObjectId()
689  << ", rawId: " << refAli->geomDetId().rawId() << ", subdetId: " << detid.subdetId() << "): " << diff << check;
690  throw cms::Exception("Tolerance in TrackerGeometryCompare exceeded");
691  }
692 
693  AlgebraicVector TRtot(12);
694  // global
695  TRtot(1) = Rtotal.x();
696  TRtot(2) = Rtotal.y();
697  TRtot(3) = Rtotal.z();
698  TRtot(4) = Wtotal.x();
699  TRtot(5) = Wtotal.y();
700  TRtot(6) = Wtotal.z();
701  // local
702  TRtot(7) = lRtotal.x();
703  TRtot(8) = lRtotal.y();
704  TRtot(9) = lRtotal.z();
705  TRtot(10) = lWtotal.x();
706  TRtot(11) = lWtotal.y();
707  TRtot(12) = lWtotal.z();
708 
709  fillTree(refAli, TRtot, tTopo, iSetup);
710  }
711 
712  // another added level for difference between det and detunit
713  for (unsigned int i = 0; i < nComp; ++i)
714  compareGeometries(refComp[i], curComp[i], tTopo, iSetup);
715 }
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:605
TrackerGeometryCompare::m_theLevels
std::vector< align::StructureType > m_theLevels
Definition: TrackerGeometryCompare.h:58
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:119
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:796
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
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 453 of file TrackerGeometryCompare.cc.

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

Referenced by analyze().

453  {
454  if (inputFilename1_ != "IDEAL" && inputFilename2_ != "IDEAL") {
455  int inputRawid1;
456  int inputRawid2;
457  int inputSubdetid1, inputSubdetid2;
458  int inputDtype1, inputDtype2;
459  std::vector<double> inputDpar1;
460  std::vector<double> inputDpar2;
461  std::vector<double>* p_inputDpar1 = &inputDpar1;
462  std::vector<double>* p_inputDpar2 = &inputDpar2;
463 
464  TTree* refTree = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
465  refTree->SetBranchAddress("irawid", &inputRawid1);
466  refTree->SetBranchAddress("subdetid", &inputSubdetid1);
467  refTree->SetBranchAddress("dtype", &inputDtype1);
468  refTree->SetBranchAddress("dpar", &p_inputDpar1);
469 
470  TTree* curTree = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
471  curTree->SetBranchAddress("irawid", &inputRawid2);
472  curTree->SetBranchAddress("subdetid", &inputSubdetid2);
473  curTree->SetBranchAddress("dtype", &inputDtype2);
474  curTree->SetBranchAddress("dpar", &p_inputDpar2);
475 
476  unsigned int nEntries11 = refTree->GetEntries();
477  unsigned int nEntries12 = curTree->GetEntries();
478 
479  if (nEntries11 != nEntries12) {
480  edm::LogError("TrackerGeometryCompare") << " Surface deformation parameters in two geometries differ!\n";
481  return;
482  }
483 
484  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
485  refTree->GetEntry(iEntry);
486  curTree->GetEntry(iEntry);
487  for (int ii = 0; ii < 13; ++ii) {
488  surfDeform_[ii] = -1.0;
489  }
490  for (int npar = 0; npar < int(inputDpar2.size()); ++npar) {
491  if (inputRawid1 == inputRawid2) {
492  surfDeform_[npar] = inputDpar2.at(npar) - inputDpar1.at(npar);
493  std::stringstream histname0;
494  histname0 << "SurfDeform_Par_" << npar;
495  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
496  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
497  if (inputSubdetid1 == 1 && inputSubdetid2 == 1) {
498  std::stringstream histname1;
499  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
500  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
501  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
502  }
503  if (inputSubdetid1 == 2 && inputSubdetid2 == 2) {
504  std::stringstream histname2;
505  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
506  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
507  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
508  }
509  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
510  }
511  }
512  }
513 
514  } else if (inputFilename1_ == "IDEAL" && inputFilename2_ != "IDEAL") {
515  int inputRawid2;
516  int inputSubdetid2;
517  int inputDtype2;
518  std::vector<double> inputDpar2;
519  std::vector<double>* p_inputDpar2 = &inputDpar2;
520 
521  TTree* curTree = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
522  curTree->SetBranchAddress("irawid", &inputRawid2);
523  curTree->SetBranchAddress("subdetid", &inputSubdetid2);
524  curTree->SetBranchAddress("dtype", &inputDtype2);
525  curTree->SetBranchAddress("dpar", &p_inputDpar2);
526 
527  unsigned int nEntries12 = curTree->GetEntries();
528 
529  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
530  curTree->GetEntry(iEntry);
531  for (int ii = 0; ii < 12; ++ii) {
532  surfDeform_[ii] = -1.0;
533  }
534  for (int npar = 0; npar < int(inputDpar2.size()); ++npar) {
535  surfDeform_[npar] = inputDpar2.at(npar);
536  std::stringstream histname0;
537  histname0 << "SurfDeform_Par_" << npar;
538  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
539  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
540  if (inputSubdetid2 == 1) {
541  std::stringstream histname1;
542  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
543  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
544  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
545  }
546  if (inputSubdetid2 == 2) {
547  std::stringstream histname2;
548  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
549  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
550  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
551  }
552  (m_vtkmap_.at(npar)).fill_current_val(inputRawid2, surfDeform_[npar]);
553  }
554  }
555 
556  } else if (inputFilename1_ != "IDEAL" && inputFilename2_ == "IDEAL") {
557  int inputRawid1;
558  int inputSubdetid1;
559  int inputDtype1;
560  std::vector<double> inputDpar1;
561  std::vector<double>* p_inputDpar1 = &inputDpar1;
562 
563  TTree* refTree = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
564  refTree->SetBranchAddress("irawid", &inputRawid1);
565  refTree->SetBranchAddress("subdetid", &inputSubdetid1);
566  refTree->SetBranchAddress("dtype", &inputDtype1);
567  refTree->SetBranchAddress("dpar", &p_inputDpar1);
568 
569  unsigned int nEntries11 = refTree->GetEntries();
570 
571  for (unsigned int iEntry = 0; iEntry < nEntries11; ++iEntry) {
572  refTree->GetEntry(iEntry);
573  for (int ii = 0; ii < 12; ++ii) {
574  surfDeform_[ii] = -1.0;
575  }
576  for (int npar = 0; npar < int(inputDpar1.size()); ++npar) {
577  surfDeform_[npar] = -inputDpar1.at(npar);
578  std::stringstream histname0;
579  histname0 << "SurfDeform_Par_" << npar;
580  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
581  m_h1_[histname0.str()]->Fill(surfDeform_[npar]);
582  if (inputSubdetid1 == 1) {
583  std::stringstream histname1;
584  histname1 << "SurfDeform_PixelBarrel_Par_" << npar;
585  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
586  m_h1_[histname1.str()]->Fill(surfDeform_[npar]);
587  }
588  if (inputSubdetid1 == 2) {
589  std::stringstream histname2;
590  histname2 << "SurfDeform_PixelEndcap_Par_" << npar;
591  if (TMath::Abs(surfDeform_[npar]) > (m_rangeHigh_ - m_rangeLow_) / (10. * m_nBins_))
592  m_h1_[histname2.str()]->Fill(surfDeform_[npar]);
593  }
594  (m_vtkmap_.at(npar)).fill_current_val(inputRawid1, surfDeform_[npar]);
595  }
596  }
597 
598  } else if (inputFilename1_ == "IDEAL" && inputFilename2_ == "IDEAL") {
599  edm::LogInfo("TrackerGeometryCompare") << ">>>> Comparing IDEAL with IDEAL: nothing to do! <<<<\n";
600  }
601 
602  return;
603 }
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:162
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
TrackerGeometryCompare::m_h1_
std::map< std::string, TH1D * > m_h1_
Definition: TrackerGeometryCompare.h:164
cuy.ii
ii
Definition: cuy.py:589
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125

◆ createROOTGeometry()

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

Definition at line 275 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_, mps_splice::line, Alignments::m_align, AlignmentErrorsExtended::m_alignError, moduleList_, moduleListFile_, moduleListName_, ptitpToken_, ptpToken_, referenceTracker, jetsAK4_CHS_cff::sort, AlCaHLTBitMon_QueryRunRegistry::string, and topoToken_.

Referenced by analyze().

275  {
276  int inputRawId1, inputRawId2;
277  double inputX1, inputY1, inputZ1, inputX2, inputY2, inputZ2;
278  double inputAlpha1, inputBeta1, inputGamma1, inputAlpha2, inputBeta2, inputGamma2;
279 
280  //Retrieve tracker topology from geometry
281  const TrackerTopology* const tTopo = &iSetup.getData(topoToken_);
282 
283  // Fill module IDs from file into a list
284  moduleListFile_.open(moduleListName_);
285  if (moduleListFile_.is_open()) {
286  std::string line;
287  while (!moduleListFile_.eof()) {
288  std::getline(moduleListFile_, line);
289  moduleList_.push_back(std::atoi(line.c_str()));
290  }
291  } else {
292  edm::LogInfo("TrackerGeometryCompare") << "Error: Module list not found! Please verify that given list exists!";
293  }
294 
295  //declare alignments
296  Alignments* alignments1 = new Alignments();
297  AlignmentErrorsExtended* alignmentErrors1 = new AlignmentErrorsExtended();
298  if (inputFilename1_ != "IDEAL") {
299  inputRootFile1_ = new TFile(inputFilename1_.c_str());
300  TTree* inputTree01_ = (TTree*)inputRootFile1_->Get(inputTreenameAlign_.c_str());
301  inputTree01_->SetBranchAddress("rawid", &inputRawId1);
302  inputTree01_->SetBranchAddress("x", &inputX1);
303  inputTree01_->SetBranchAddress("y", &inputY1);
304  inputTree01_->SetBranchAddress("z", &inputZ1);
305  inputTree01_->SetBranchAddress("alpha", &inputAlpha1);
306  inputTree01_->SetBranchAddress("beta", &inputBeta1);
307  inputTree01_->SetBranchAddress("gamma", &inputGamma1);
308 
309  int nEntries1 = inputTree01_->GetEntries();
310  //fill alignments
311  for (int i = 0; i < nEntries1; ++i) {
312  inputTree01_->GetEntry(i);
313  CLHEP::Hep3Vector translation1(inputX1, inputY1, inputZ1);
314  CLHEP::HepEulerAngles eulerangles1(inputAlpha1, inputBeta1, inputGamma1);
315  uint32_t detid1 = inputRawId1;
316  AlignTransform transform1(translation1, eulerangles1, detid1);
317  alignments1->m_align.push_back(transform1);
318 
319  //dummy errors
320  CLHEP::HepSymMatrix clhepSymMatrix(3, 0);
321  AlignTransformErrorExtended transformError(clhepSymMatrix, detid1);
322  alignmentErrors1->m_alignError.push_back(transformError);
323  }
324 
325  // to get the right order, sort by rawId
326  std::sort(alignments1->m_align.begin(), alignments1->m_align.end());
327  std::sort(alignmentErrors1->m_alignError.begin(), alignmentErrors1->m_alignError.end());
328  }
329  //------------------
330  Alignments* alignments2 = new Alignments();
331  AlignmentErrorsExtended* alignmentErrors2 = new AlignmentErrorsExtended();
332  if (inputFilename2_ != "IDEAL") {
333  inputRootFile2_ = new TFile(inputFilename2_.c_str());
334  TTree* inputTree02_ = (TTree*)inputRootFile2_->Get(inputTreenameAlign_.c_str());
335  inputTree02_->SetBranchAddress("rawid", &inputRawId2);
336  inputTree02_->SetBranchAddress("x", &inputX2);
337  inputTree02_->SetBranchAddress("y", &inputY2);
338  inputTree02_->SetBranchAddress("z", &inputZ2);
339  inputTree02_->SetBranchAddress("alpha", &inputAlpha2);
340  inputTree02_->SetBranchAddress("beta", &inputBeta2);
341  inputTree02_->SetBranchAddress("gamma", &inputGamma2);
342 
343  int nEntries2 = inputTree02_->GetEntries();
344  //fill alignments
345  for (int i = 0; i < nEntries2; ++i) {
346  inputTree02_->GetEntry(i);
347  CLHEP::Hep3Vector translation2(inputX2, inputY2, inputZ2);
348  CLHEP::HepEulerAngles eulerangles2(inputAlpha2, inputBeta2, inputGamma2);
349  uint32_t detid2 = inputRawId2;
350  AlignTransform transform2(translation2, eulerangles2, detid2);
351  alignments2->m_align.push_back(transform2);
352 
353  //dummy errors
354  CLHEP::HepSymMatrix clhepSymMatrix(3, 0);
355  AlignTransformErrorExtended transformError(clhepSymMatrix, detid2);
356  alignmentErrors2->m_alignError.push_back(transformError);
357  }
358 
359  // to get the right order, sort by rawId
360  std::sort(alignments2->m_align.begin(), alignments2->m_align.end());
361  std::sort(alignmentErrors2->m_alignError.begin(), alignmentErrors2->m_alignError.end());
362  }
363 
364  //accessing the initial geometry
365  if (!fromDD4hep_) {
366  edm::ESTransientHandle<DDCompactView> cpv = iSetup.getTransientHandle(cpvTokenDDD_);
367  } else {
368  edm::ESTransientHandle<cms::DDCompactView> cpv = iSetup.getTransientHandle(cpvTokenDD4hep_);
369  }
370 
371  const GeometricDet* theGeometricDet = &iSetup.getData(geomDetToken_);
372  const PTrackerParameters* ptp = &iSetup.getData(ptpToken_);
373  const PTrackerAdditionalParametersPerDet* ptitp = &iSetup.getData(ptitpToken_);
374  TrackerGeomBuilderFromGeometricDet trackerBuilder;
375 
376  //reference tracker
377  TrackerGeometry* theRefTracker = trackerBuilder.build(theGeometricDet, ptitp, *ptp, tTopo);
378  if (inputFilename1_ != "IDEAL") {
379  GeometryAligner aligner1;
380  aligner1.applyAlignments<TrackerGeometry>(
381  &(*theRefTracker), &(*alignments1), &(*alignmentErrors1), AlignTransform());
382  }
383  referenceTracker = new AlignableTracker(&(*theRefTracker), tTopo);
384  //referenceTracker->setSurfaceDeformation(surfDef1, true) ;
385 
386  int inputRawid1;
387  int inputRawid2;
388  int inputDtype1, inputDtype2;
389  std::vector<double> inputDpar1;
390  std::vector<double> inputDpar2;
391  std::vector<double>* p_inputDpar1 = &inputDpar1;
392  std::vector<double>* p_inputDpar2 = &inputDpar2;
393 
394  const auto& comp1 = referenceTracker->deepComponents();
395 
396  SurfaceDeformation* surfDef1;
397  if (inputFilename1_ != "IDEAL") {
398  TTree* inputTree11_ = (TTree*)inputRootFile1_->Get(inputTreenameDeform_.c_str());
399  inputTree11_->SetBranchAddress("irawid", &inputRawid1);
400  inputTree11_->SetBranchAddress("dtype", &inputDtype1);
401  inputTree11_->SetBranchAddress("dpar", &p_inputDpar1);
402 
403  unsigned int nEntries11 = inputTree11_->GetEntries();
404  edm::LogInfo("TrackerGeometryCompare") << " nentries11 = " << nEntries11;
405  for (unsigned int iEntry = 0; iEntry < nEntries11; ++iEntry) {
406  inputTree11_->GetEntry(iEntry);
407 
408  surfDef1 = SurfaceDeformationFactory::create(inputDtype1, inputDpar1);
409 
410  if (int(comp1[iEntry]->id()) == inputRawid1) {
411  comp1[iEntry]->setSurfaceDeformation(surfDef1, true);
412  }
413  }
414  }
415 
416  //currernt tracker
417  TrackerGeometry* theCurTracker = trackerBuilder.build(&*theGeometricDet, ptitp, *ptp, tTopo);
418  if (inputFilename2_ != "IDEAL") {
419  GeometryAligner aligner2;
420  aligner2.applyAlignments<TrackerGeometry>(
421  &(*theCurTracker), &(*alignments2), &(*alignmentErrors2), AlignTransform());
422  }
423  currentTracker = new AlignableTracker(&(*theCurTracker), tTopo);
424 
425  const auto& comp2 = currentTracker->deepComponents();
426 
427  SurfaceDeformation* surfDef2;
428  if (inputFilename2_ != "IDEAL") {
429  TTree* inputTree12_ = (TTree*)inputRootFile2_->Get(inputTreenameDeform_.c_str());
430  inputTree12_->SetBranchAddress("irawid", &inputRawid2);
431  inputTree12_->SetBranchAddress("dtype", &inputDtype2);
432  inputTree12_->SetBranchAddress("dpar", &p_inputDpar2);
433 
434  unsigned int nEntries12 = inputTree12_->GetEntries();
435  edm::LogInfo("TrackerGeometryCompare") << " nentries12 = " << nEntries12;
436  for (unsigned int iEntry = 0; iEntry < nEntries12; ++iEntry) {
437  inputTree12_->GetEntry(iEntry);
438 
439  surfDef2 = SurfaceDeformationFactory::create(inputDtype2, inputDpar2);
440 
441  if (int(comp2[iEntry]->id()) == inputRawid2) {
442  comp2[iEntry]->setSurfaceDeformation(surfDef2, true);
443  }
444  }
445  }
446 
447  delete alignments1;
448  delete alignmentErrors1;
449  delete alignments2;
450  delete alignmentErrors2;
451 }
TrackerGeometryCompare::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: TrackerGeometryCompare.h:91
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometryCompare::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: TrackerGeometryCompare.h:89
GeometryAligner
Class to update a given geometry with a set of alignments.
Definition: GeometryAligner.h:33
TrackerGeomBuilderFromGeometricDet::build
TrackerGeometry * build(const GeometricDet *gd, const PTrackerAdditionalParametersPerDet *ptitp, const PTrackerParameters &ptp, const TrackerTopology *tTopo)
Definition: TrackerGeomBuilderFromGeometricDet.cc:44
TrackerGeometryCompare::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: TrackerGeometryCompare.h:92
TrackerGeometryCompare::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: TrackerGeometryCompare.h:90
Alignable::deepComponents
const Alignables & deepComponents() const
Definition: Alignable.h:72
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:128
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
TrackerGeometryCompare::cpvTokenDDD_
const edm::ESGetToken< DDCompactView, IdealGeometryRecord > cpvTokenDDD_
Definition: TrackerGeometryCompare.h:87
mps_splice.line
line
Definition: mps_splice.py:76
jetsAK4_CHS_cff.sort
sort
Definition: jetsAK4_CHS_cff.py:39
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:98
AlignmentErrorsExtended::m_alignError
std::vector< AlignTransformErrorExtended > m_alignError
Definition: AlignmentErrorsExtended.h:19
TrackerGeometryCompare::cpvTokenDD4hep_
const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord > cpvTokenDD4hep_
Definition: TrackerGeometryCompare.h:88
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:96
edm::EventSetup::getTransientHandle
ESTransientHandle< T > getTransientHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:162
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 761 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().

761  {
762  const auto& refComp = refAli->components();
763  const auto& curComp = curAli->components();
764 
765  unsigned int nComp = refComp.size();
766  //only perform for designate levels
767  bool useLevel = false;
768  if (refAli->alignableObjectId() == commonTrackerLevel_)
769  useLevel = true;
770 
771  //useLevel = false;
772  if (useLevel) {
773  CLHEP::Hep3Vector Rtotal, Wtotal;
774  Rtotal.set(0., 0., 0.);
775  Wtotal.set(0., 0., 0.);
776 
777  AlgebraicVector diff = align::diffAlignables(refAli, curAli, weightBy_, weightById_, weightByIdVector_);
778  CLHEP::Hep3Vector dR(diff[0], diff[1], diff[2]);
779  Rtotal += dR;
780  CLHEP::Hep3Vector dW(diff[3], diff[4], diff[5]);
781  CLHEP::HepRotation rot(Wtotal.unit(), Wtotal.mag());
782  CLHEP::HepRotation drot(dW.unit(), dW.mag());
783  rot *= drot;
784  Wtotal.set(rot.axis().x() * rot.delta(), rot.axis().y() * rot.delta(), rot.axis().z() * rot.delta());
785 
786  TrackerCommonT_ = align::GlobalVector(Rtotal.x(), Rtotal.y(), Rtotal.z());
787  TrackerCommonR_ = align::GlobalVector(Wtotal.x(), Wtotal.y(), Wtotal.z());
788  TrackerCommonCM_ = curAli->globalPosition();
789 
790  } else {
791  for (unsigned int i = 0; i < nComp; ++i)
792  diffCommonTrackerSystem(refComp[i], curComp[i]);
793  }
794 }
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometryCompare::TrackerCommonR_
align::GlobalVector TrackerCommonR_
Definition: TrackerGeometryCompare.h:124
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:761
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:119
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:123
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:122
TrackerGeometryCompare::TrackerCommonCM_
align::PositionType TrackerCommonCM_
Definition: TrackerGeometryCompare.h:125

◆ endJob()

void TrackerGeometryCompare::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 208 of file TrackerGeometryCompare.cc.

References alignTree_, m_vtkmap_, AlCaHLTBitMon_QueryRunRegistry::string, surfdir_, and theFile_.

Referenced by o2olib.O2ORunMgr::executeJob().

208  {
209  int iname(0);
210  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
211  std::stringstream mapname;
212  mapname << surfdir_ << "/TkMap_SurfDeform" << iname << ".png";
213  it->save(true, 0, 0, mapname.str());
214  mapname.str(std::string());
215  mapname.clear();
216  mapname << surfdir_ << "/TkMap_SurfDeform" << iname << ".pdf";
217  it->save(true, 0, 0, mapname.str());
218  ++iname;
219  }
220 
221  theFile_->cd();
222  alignTree_->Write();
223  theFile_->Close();
224 }
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:162
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 999 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().

999  {
1000  switch (subdetlevel) {
1001  case 1: {
1002  identifiers_[0] = tTopo->pxbModule(rawid);
1003  identifiers_[1] = tTopo->pxbLadder(rawid);
1004  identifiers_[2] = tTopo->pxbLayer(rawid);
1005  identifiers_[3] = 999;
1006  identifiers_[4] = 999;
1007  identifiers_[5] = 999;
1008  break;
1009  }
1010  case 2: {
1011  identifiers_[0] = tTopo->pxfModule(rawid);
1012  identifiers_[1] = tTopo->pxfPanel(rawid);
1013  identifiers_[2] = tTopo->pxfBlade(rawid);
1014  identifiers_[3] = tTopo->pxfDisk(rawid);
1015  identifiers_[4] = tTopo->pxfSide(rawid);
1016  identifiers_[5] = 999;
1017  break;
1018  }
1019  case 3: {
1020  identifiers_[0] = tTopo->tibModule(rawid);
1021  identifiers_[1] = tTopo->tibStringInfo(rawid)[0];
1022  identifiers_[2] = tTopo->tibStringInfo(rawid)[1];
1023  identifiers_[3] = tTopo->tibStringInfo(rawid)[2];
1024  identifiers_[4] = tTopo->tibLayer(rawid);
1025  identifiers_[5] = 999;
1026  break;
1027  }
1028  case 4: {
1029  identifiers_[0] = tTopo->tidModuleInfo(rawid)[0];
1030  identifiers_[1] = tTopo->tidModuleInfo(rawid)[1];
1031  identifiers_[2] = tTopo->tidRing(rawid);
1032  identifiers_[3] = tTopo->tidWheel(rawid);
1033  identifiers_[4] = tTopo->tidSide(rawid);
1034  identifiers_[5] = 999;
1035  break;
1036  }
1037  case 5: {
1038  identifiers_[0] = tTopo->tobModule(rawid);
1039  identifiers_[1] = tTopo->tobRodInfo(rawid)[0];
1040  identifiers_[2] = tTopo->tobRodInfo(rawid)[1];
1041  identifiers_[3] = tTopo->tobLayer(rawid);
1042  identifiers_[4] = 999;
1043  identifiers_[5] = 999;
1044  break;
1045  }
1046  case 6: {
1047  identifiers_[0] = tTopo->tecModule(rawid);
1048  identifiers_[1] = tTopo->tecRing(rawid);
1049  identifiers_[2] = tTopo->tecPetalInfo(rawid)[0];
1050  identifiers_[3] = tTopo->tecPetalInfo(rawid)[1];
1051  identifiers_[4] = tTopo->tecWheel(rawid);
1052  identifiers_[5] = tTopo->tecSide(rawid);
1053  break;
1054  }
1055  default: {
1056  edm::LogInfo("TrackerGeometryCompare") << "Error: bad subdetid!!";
1057  break;
1058  }
1059  }
1060 }
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:147
TrackerTopology::pxbLayer
unsigned int pxbLayer(const DetId &id) const
Definition: TrackerTopology.h:144
TrackerTopology::tidModuleInfo
std::vector< unsigned int > tidModuleInfo(const DetId &id) const
Definition: TrackerTopology.h:434
TrackerTopology::pxfBlade
unsigned int pxfBlade(const DetId &id) const
Definition: TrackerTopology.h:447
TrackerTopology::tibModule
unsigned int tibModule(const DetId &id) const
Definition: TrackerTopology.h:172
TrackerTopology::tidSide
unsigned int tidSide(const DetId &id) const
Definition: TrackerTopology.h:190
TrackerTopology::pxfModule
unsigned int pxfModule(const DetId &id) const
Definition: TrackerTopology.h:163
TrackerTopology::tecPetalInfo
std::vector< unsigned int > tecPetalInfo(const DetId &id) const
Definition: TrackerTopology.h:407
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:201
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:198
TrackerTopology::pxbLadder
unsigned int pxbLadder(const DetId &id) const
Definition: TrackerTopology.h:155
TrackerTopology::tecRing
unsigned int tecRing(const DetId &id) const
ring id
Definition: TrackerTopology.h:217
TrackerTopology::tecModule
unsigned int tecModule(const DetId &id) const
Definition: TrackerTopology.h:169
TrackerTopology::tecSide
unsigned int tecSide(const DetId &id) const
Definition: TrackerTopology.h:184
TrackerTopology::pxfDisk
unsigned int pxfDisk(const DetId &id) const
Definition: TrackerTopology.h:446
TrackerTopology::pxfPanel
unsigned int pxfPanel(const DetId &id) const
Definition: TrackerTopology.h:450
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
TrackerTopology::pxfSide
unsigned int pxfSide(const DetId &id) const
Definition: TrackerTopology.h:192
TrackerTopology::tibStringInfo
std::vector< unsigned int > tibStringInfo(const DetId &id) const
Definition: TrackerTopology.h:421
TrackerTopology::tobRodInfo
std::vector< unsigned int > tobRodInfo(const DetId &id) const
Definition: TrackerTopology.h:226
TrackerTopology::tidRing
unsigned int tidRing(const DetId &id) const
Definition: TrackerTopology.h:218
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:150
TrackerTopology::tobModule
unsigned int tobModule(const DetId &id) const
Definition: TrackerTopology.h:166
TrackerTopology::pxbModule
unsigned int pxbModule(const DetId &id) const
Definition: TrackerTopology.h:160

◆ fillTree()

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

Definition at line 796 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().

799  {
800  //Get bad modules
801  const SiPixelQuality* SiPixelModules = &iSetup.getData(pixQualityToken_);
802  const SiStripQuality* SiStripModules = &iSetup.getData(stripQualityToken_);
803 
804  id_ = refAli->id();
805 
806  badModuleQuality_ = 0;
807  //check if module has a bad quality tag
808  if (SiPixelModules->IsModuleBad(id_)) {
809  badModuleQuality_ = 1;
810  }
811  if (SiStripModules->IsModuleBad(id_)) {
812  badModuleQuality_ = 1;
813  }
814 
815  //check if module is in a given list of bad/untouched etc. modules
816  inModuleList_ = 0;
817  for (unsigned int i = 0; i < moduleList_.size(); i++) {
818  if (moduleList_[i] == id_) {
819  inModuleList_ = 1;
820  break;
821  }
822  }
823 
824  level_ = refAli->alignableObjectId();
825  //need if ali has no mother
826  if (refAli->mother()) {
827  mid_ = refAli->mother()->geomDetId().rawId();
828  mlevel_ = refAli->mother()->alignableObjectId();
829  } else {
830  mid_ = -1;
831  mlevel_ = -1;
832  }
833  DetId detid(id_);
834  sublevel_ = detid.subdetId();
835  fillIdentifiers(sublevel_, id_, tTopo);
836  xVal_ = refAli->globalPosition().x();
837  yVal_ = refAli->globalPosition().y();
838  zVal_ = refAli->globalPosition().z();
839  align::GlobalVector vec(xVal_, yVal_, zVal_);
840  rVal_ = vec.perp();
841  phiVal_ = vec.phi();
842  etaVal_ = vec.eta();
843  align::RotationType rot = refAli->globalRotation();
844  align::EulerAngles eulerAngles = align::toAngles(rot);
845  alphaVal_ = eulerAngles[0];
846  betaVal_ = eulerAngles[1];
847  gammaVal_ = eulerAngles[2];
848  // global
849  dxVal_ = diff[0];
850  dyVal_ = diff[1];
851  dzVal_ = diff[2];
852  // local
853  duVal_ = diff[6];
854  dvVal_ = diff[7];
855  dwVal_ = diff[8];
856  //...TODO...
857  align::GlobalVector g(dxVal_, dyVal_, dzVal_);
858  //getting dR and dPhi
859  align::GlobalVector vRef(xVal_, yVal_, zVal_);
860  align::GlobalVector vCur(xVal_ + dxVal_, yVal_ + dyVal_, zVal_ + dzVal_);
861  drVal_ = vCur.perp() - vRef.perp();
862  dphiVal_ = vCur.phi() - vRef.phi();
863  // global
864  dalphaVal_ = diff[3];
865  dbetaVal_ = diff[4];
866  dgammaVal_ = diff[5];
867  // local
868  daVal_ = diff[9];
869  dbVal_ = diff[10];
870  dgVal_ = diff[11];
871 
872  //detIdFlag
873  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
874  if (detIdFlag_) {
875  if ((passIdCut(refAli->id())) || (passIdCut(refAli->mother()->id()))) {
876  useDetId_ = 1;
877  } else {
878  useDetId_ = 0;
879  }
880  }
881  }
882  // det module dimension
883  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
884  if (refAli->mother()->alignableObjectId() != align::AlignableDet)
885  detDim_ = 1;
886  else if (refAli->mother()->alignableObjectId() == align::AlignableDet)
887  detDim_ = 2;
888  } else
889  detDim_ = 0;
890 
891  surWidth_ = refAli->surface().width();
892  surLength_ = refAli->surface().length();
893  align::RotationType rt = refAli->globalRotation();
894  surRot_[0] = rt.xx();
895  surRot_[1] = rt.xy();
896  surRot_[2] = rt.xz();
897  surRot_[3] = rt.yx();
898  surRot_[4] = rt.yy();
899  surRot_[5] = rt.yz();
900  surRot_[6] = rt.zx();
901  surRot_[7] = rt.zy();
902  surRot_[8] = rt.zz();
903 
904  //Fill
905  alignTree_->Fill();
906 }
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
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:128
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:94
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
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
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:999
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:93

◆ passIdCut()

bool TrackerGeometryCompare::passIdCut ( uint32_t  id)
private

Definition at line 987 of file TrackerGeometryCompare.cc.

References detIdFlagVector_, and mps_fire::i.

Referenced by fillTree().

987  {
988  bool pass = false;
989  int nEntries = detIdFlagVector_.size();
990 
991  for (int i = 0; i < nEntries; i++) {
992  if (detIdFlagVector_[i] == id)
993  pass = true;
994  }
995 
996  return pass;
997 }
mps_fire.i
i
Definition: mps_fire.py:429
TrackerGeometryCompare::detIdFlagVector_
std::vector< uint32_t > detIdFlagVector_
Definition: TrackerGeometryCompare.h:121

◆ setCommonTrackerSystem()

void TrackerGeometryCompare::setCommonTrackerSystem ( )
private

Definition at line 717 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().

717  {
718  edm::LogInfo("TrackerGeometryCompare") << "Setting Common Tracker System....";
719 
720  // DM_534??AlignableObjectId dummy;
721  // DM_534??_commonTrackerLevel = dummy.nameToType(_setCommonTrackerSystem);
722  commonTrackerLevel_ = currentTracker->objectIdProvider().stringToId(setCommonTrackerSystem_);
723 
724  diffCommonTrackerSystem(referenceTracker, currentTracker);
725 
726  align::EulerAngles dOmega(3);
727  dOmega[0] = TrackerCommonR_.x();
728  dOmega[1] = TrackerCommonR_.y();
729  dOmega[2] = TrackerCommonR_.z();
730  align::RotationType rot = align::toMatrix(dOmega);
731  align::GlobalVector theR = TrackerCommonT_;
732 
733  edm::LogInfo("TrackerGeometryCompare") << "what we get from overlaying the pixels..." << theR << ", " << rot;
734 
735  //transform to the Tracker System
736  align::PositionType trackerCM = currentTracker->globalPosition();
737  align::GlobalVector cmDiff(
738  trackerCM.x() - TrackerCommonCM_.x(), trackerCM.y() - TrackerCommonCM_.y(), trackerCM.z() - TrackerCommonCM_.z());
739 
740  edm::LogInfo("TrackerGeometryCompare") << "Pixel CM: " << TrackerCommonCM_ << ", tracker CM: " << trackerCM;
741 
742  //adjust translational difference factoring in different rotational CM
743  //needed because rotateInGlobalFrame is about CM of alignable, not Tracker
744  const align::GlobalVector::BasicVectorType& lpvgf = cmDiff.basicVector();
745  align::GlobalVector moveV(rot.multiplyInverse(lpvgf) - lpvgf);
746  align::GlobalVector theRprime(theR + moveV);
747 
748  AlgebraicVector TrackerCommonTR(6);
749  TrackerCommonTR(1) = theRprime.x();
750  TrackerCommonTR(2) = theRprime.y();
751  TrackerCommonTR(3) = theRprime.z();
752  TrackerCommonTR(4) = TrackerCommonR_.x();
753  TrackerCommonTR(5) = TrackerCommonR_.y();
754  TrackerCommonTR(6) = TrackerCommonR_.z();
755 
756  edm::LogInfo("TrackerGeometryCompare") << "and after the transformation: " << TrackerCommonTR;
757 
758  align::moveAlignable(currentTracker, TrackerCommonTR);
759 }
TrackerGeometryCompare::TrackerCommonR_
align::GlobalVector TrackerCommonR_
Definition: TrackerGeometryCompare.h:124
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
TrackerGeometryCompare::diffCommonTrackerSystem
void diffCommonTrackerSystem(Alignable *refAli, Alignable *curAli)
Definition: TrackerGeometryCompare.cc:761
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:123
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:13
align::EulerAngles
AlgebraicVector EulerAngles
Definition: Definitions.h:34
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:98
TrackerGeometryCompare::commonTrackerLevel_
align::StructureType commonTrackerLevel_
Definition: TrackerGeometryCompare.h:122
AlignableObjectId::stringToId
align::StructureType stringToId(const char *) const
Definition: AlignableObjectId.cc:284
TrackerGeometryCompare::referenceTracker
AlignableTracker * referenceTracker
Definition: TrackerGeometryCompare.h:96
TrackerGeometryCompare::TrackerCommonCM_
align::PositionType TrackerCommonCM_
Definition: TrackerGeometryCompare.h:125
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 908 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(), jetsAK4_CHS_cff::sort, AlignableTracker::TIDGeomDets(), and HcalDetIdTransform::transform().

910  {
911  //getting the right alignables for the alignment record
912  auto detPB = ali->pixelHalfBarrelGeomDets();
913  auto detPEC = ali->pixelEndcapGeomDets();
914  auto detTIB = ali->innerBarrelGeomDets();
915  auto detTID = ali->TIDGeomDets();
916  auto detTOB = ali->outerBarrelGeomDets();
917  auto detTEC = ali->endcapGeomDets();
918 
919  align::Alignables allGeomDets;
920  std::copy(detPB.begin(), detPB.end(), std::back_inserter(allGeomDets));
921  std::copy(detPEC.begin(), detPEC.end(), std::back_inserter(allGeomDets));
922  std::copy(detTIB.begin(), detTIB.end(), std::back_inserter(allGeomDets));
923  std::copy(detTID.begin(), detTID.end(), std::back_inserter(allGeomDets));
924  std::copy(detTOB.begin(), detTOB.end(), std::back_inserter(allGeomDets));
925  std::copy(detTEC.begin(), detTEC.end(), std::back_inserter(allGeomDets));
926 
927  align::Alignables rcdAlis;
928  for (const auto& i : allGeomDets) {
929  if (i->components().size() == 1) {
930  rcdAlis.push_back(i);
931  } else if (i->components().size() > 1) {
932  rcdAlis.push_back(i);
933  const auto& comp = i->components();
934  for (const auto& j : comp)
935  rcdAlis.push_back(j);
936  }
937  }
938 
939  //turning them into alignments
940  for (const auto& k : rcdAlis) {
941  const SurveyDet* surveyInfo = k->survey();
942  const align::PositionType& pos(surveyInfo->position());
943  align::RotationType rot(surveyInfo->rotation());
944  CLHEP::Hep3Vector clhepVector(pos.x(), pos.y(), pos.z());
945  CLHEP::HepRotation clhepRotation(
946  CLHEP::HepRep3x3(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
947  AlignTransform transform(clhepVector, clhepRotation, k->id());
948  AlignTransformErrorExtended transformError(CLHEP::HepSymMatrix(3, 1), k->id());
949  alignVals->m_align.push_back(transform);
950  alignErrors->m_alignError.push_back(transformError);
951  }
952 
953  // to get the right order, sort by rawId
954  std::sort(alignVals->m_align.begin(), alignVals->m_align.end());
955  std::sort(alignErrors->m_alignError.begin(), alignErrors->m_alignError.end());
956 }
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
jetsAK4_CHS_cff.sort
sort
Definition: jetsAK4_CHS_cff.py:39
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 135 of file TrackerGeometryCompare.h.

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

◆ alphaVal_

float TrackerGeometryCompare::alphaVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ badModuleQuality_

int TrackerGeometryCompare::badModuleQuality_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ betaVal_

float TrackerGeometryCompare::betaVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ commonTrackerLevel_

align::StructureType TrackerGeometryCompare::commonTrackerLevel_
private

Definition at line 122 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ cpvTokenDD4hep_

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

Definition at line 88 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ cpvTokenDDD_

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

Definition at line 87 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ currentTracker

AlignableTracker* TrackerGeometryCompare::currentTracker
private

Definition at line 98 of file TrackerGeometryCompare.h.

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

◆ dalphaVal_

float TrackerGeometryCompare::dalphaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ daVal_

float TrackerGeometryCompare::daVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dbetaVal_

float TrackerGeometryCompare::dbetaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dbVal_

float TrackerGeometryCompare::dbVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detDim_

int TrackerGeometryCompare::detDim_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detIdFlag_

bool TrackerGeometryCompare::detIdFlag_
private

Definition at line 115 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ detIdFlagFile_

std::string TrackerGeometryCompare::detIdFlagFile_
private

Definition at line 116 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ detIdFlagVector_

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

Definition at line 121 of file TrackerGeometryCompare.h.

Referenced by passIdCut(), and TrackerGeometryCompare().

◆ dgammaVal_

float TrackerGeometryCompare::dgammaVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dgVal_

float TrackerGeometryCompare::dgVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dphiVal_

float TrackerGeometryCompare::dphiVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ drVal_

float TrackerGeometryCompare::drVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dummyTracker

AlignableTracker* TrackerGeometryCompare::dummyTracker
private

Definition at line 97 of file TrackerGeometryCompare.h.

◆ duVal_

float TrackerGeometryCompare::duVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dvVal_

float TrackerGeometryCompare::dvVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dwVal_

float TrackerGeometryCompare::dwVal_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dxVal_

float TrackerGeometryCompare::dxVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dyVal_

float TrackerGeometryCompare::dyVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ dzVal_

float TrackerGeometryCompare::dzVal_
private

Definition at line 147 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ etaVal_

float TrackerGeometryCompare::etaVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ filename_

std::string TrackerGeometryCompare::filename_
private

Definition at line 132 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ firstEvent_

bool TrackerGeometryCompare::firstEvent_
private

Definition at line 160 of file TrackerGeometryCompare.h.

Referenced by analyze(), and beginJob().

◆ fromDD4hep_

bool TrackerGeometryCompare::fromDD4hep_
private

Definition at line 111 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ gammaVal_

float TrackerGeometryCompare::gammaVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ geomDetToken_

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

Definition at line 90 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ id_

int TrackerGeometryCompare::id_
private

\ Tree variables

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ identifiers_

uint32_t TrackerGeometryCompare::identifiers_[6]
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillIdentifiers(), and TrackerGeometryCompare().

◆ inModuleList_

int TrackerGeometryCompare::inModuleList_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ inputFilename1_

std::string TrackerGeometryCompare::inputFilename1_
private

Definition at line 107 of file TrackerGeometryCompare.h.

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

◆ inputFilename2_

std::string TrackerGeometryCompare::inputFilename2_
private

Definition at line 108 of file TrackerGeometryCompare.h.

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

◆ inputRootFile1_

TFile* TrackerGeometryCompare::inputRootFile1_
private

Definition at line 136 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

◆ inputRootFile2_

TFile* TrackerGeometryCompare::inputRootFile2_
private

Definition at line 137 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

◆ inputTree01_

TTree* TrackerGeometryCompare::inputTree01_
private

Definition at line 138 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ inputTree02_

TTree* TrackerGeometryCompare::inputTree02_
private

Definition at line 139 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ inputTree11_

TTree* TrackerGeometryCompare::inputTree11_
private

Definition at line 140 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ inputTree12_

TTree* TrackerGeometryCompare::inputTree12_
private

Definition at line 141 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ inputTreenameAlign_

std::string TrackerGeometryCompare::inputTreenameAlign_
private

Definition at line 109 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

◆ inputTreenameDeform_

std::string TrackerGeometryCompare::inputTreenameDeform_
private

Definition at line 110 of file TrackerGeometryCompare.h.

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

◆ level_

int TrackerGeometryCompare::level_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ levelStrings_

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

Definition at line 105 of file TrackerGeometryCompare.h.

Referenced by analyze().

◆ m_h1_

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

Definition at line 164 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_nBins_

int TrackerGeometryCompare::m_nBins_
private

Definition at line 156 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_rangeHigh_

double TrackerGeometryCompare::m_rangeHigh_
private

Definition at line 158 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_rangeLow_

double TrackerGeometryCompare::m_rangeLow_
private

Definition at line 157 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ m_theLevels

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

Definition at line 58 of file TrackerGeometryCompare.h.

Referenced by analyze(), and compareGeometries().

◆ m_vtkmap_

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

Definition at line 162 of file TrackerGeometryCompare.h.

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

◆ mid_

int TrackerGeometryCompare::mid_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ mlevel_

int TrackerGeometryCompare::mlevel_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ moduleInList_

int TrackerGeometryCompare::moduleInList_
private

Definition at line 129 of file TrackerGeometryCompare.h.

◆ moduleList_

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

Definition at line 128 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and fillTree().

◆ moduleListFile_

std::ifstream TrackerGeometryCompare::moduleListFile_
private

Definition at line 127 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ moduleListName_

std::string TrackerGeometryCompare::moduleListName_
private

Definition at line 106 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

◆ phiVal_

float TrackerGeometryCompare::phiVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ pixQualityToken_

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

Definition at line 93 of file TrackerGeometryCompare.h.

Referenced by fillTree().

◆ ptitpToken_

const edm::ESGetToken<PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd> TrackerGeometryCompare::ptitpToken_
private

Definition at line 92 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ ptpToken_

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

Definition at line 91 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

◆ referenceTracker

AlignableTracker* TrackerGeometryCompare::referenceTracker
private

Definition at line 96 of file TrackerGeometryCompare.h.

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

◆ rVal_

float TrackerGeometryCompare::rVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ setCommonTrackerSystem_

std::string TrackerGeometryCompare::setCommonTrackerSystem_
private

Definition at line 114 of file TrackerGeometryCompare.h.

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

◆ stripQualityToken_

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

Definition at line 94 of file TrackerGeometryCompare.h.

Referenced by fillTree().

◆ sublevel_

int TrackerGeometryCompare::sublevel_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surfDeform_

double TrackerGeometryCompare::surfDeform_[13]
private

Definition at line 154 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

◆ surfdir_

std::string TrackerGeometryCompare::surfdir_
private

Definition at line 133 of file TrackerGeometryCompare.h.

Referenced by endJob(), and TrackerGeometryCompare().

◆ surLength_

float TrackerGeometryCompare::surLength_
private

Definition at line 150 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surRot_

double TrackerGeometryCompare::surRot_[9]
private

Definition at line 152 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ surWidth_

float TrackerGeometryCompare::surWidth_
private

Definition at line 150 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ theFile_

TFile* TrackerGeometryCompare::theFile_
private

Definition at line 134 of file TrackerGeometryCompare.h.

Referenced by endJob(), and TrackerGeometryCompare().

◆ theSurveyErrors

const SurveyErrors* TrackerGeometryCompare::theSurveyErrors
private

Definition at line 102 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ theSurveyIndex

unsigned int TrackerGeometryCompare::theSurveyIndex
private

Definition at line 100 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ theSurveyValues

const Alignments* TrackerGeometryCompare::theSurveyValues
private

Definition at line 101 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

◆ topoToken_

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

Definition at line 89 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

◆ TrackerCommonCM_

align::PositionType TrackerGeometryCompare::TrackerCommonCM_
private

Definition at line 125 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ TrackerCommonR_

align::GlobalVector TrackerGeometryCompare::TrackerCommonR_
private

Definition at line 124 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ TrackerCommonT_

align::GlobalVector TrackerGeometryCompare::TrackerCommonT_
private

Definition at line 123 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

◆ type_

int TrackerGeometryCompare::type_
private

Definition at line 153 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 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ weightBy_

std::string TrackerGeometryCompare::weightBy_
private

Definition at line 113 of file TrackerGeometryCompare.h.

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

◆ weightById_

bool TrackerGeometryCompare::weightById_
private

Definition at line 117 of file TrackerGeometryCompare.h.

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

◆ weightByIdFile_

std::string TrackerGeometryCompare::weightByIdFile_
private

Definition at line 118 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

◆ weightByIdVector_

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

Definition at line 119 of file TrackerGeometryCompare.h.

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

◆ writeToDB_

bool TrackerGeometryCompare::writeToDB_
private

Definition at line 112 of file TrackerGeometryCompare.h.

Referenced by analyze().

◆ xVal_

float TrackerGeometryCompare::xVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ yVal_

float TrackerGeometryCompare::yVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

◆ zVal_

float TrackerGeometryCompare::zVal_
private

Definition at line 145 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

Generated for CMSSW Reference Manual by doxygen 1.8.14