CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
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]
 
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 () noexcept
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag) noexcept
 
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

typedef AlignTransform TrackerGeometryCompare::SurveyValue

Definition at line 43 of file TrackerGeometryCompare.h.

typedef Alignments TrackerGeometryCompare::SurveyValues

Definition at line 44 of file TrackerGeometryCompare.h.

Constructor & Destructor Documentation

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

Do nothing. Required by framework.

Definition at line 55 of file TrackerGeometryCompare.cc.

References alignTree_, alphaVal_, badModuleQuality_, betaVal_, dqmiodatasetharvest::ctr, dalphaVal_, daVal_, dbetaVal_, dbVal_, detDim_, detIdFlag_, detIdFlagFile_, detIdFlagVector_, dgammaVal_, dgVal_, dphiVal_, drVal_, duVal_, dvVal_, dwVal_, dxVal_, dyVal_, dzVal_, etaVal_, filename_, groupFilesInBlocks::fin, gammaVal_, edm::ParameterSet::getUntrackedParameter(), gpuClustering::id, id_, identifiers_, cuy::ii, inModuleList_, inputFilename1_, inputFilename2_, inputTreenameAlign_, inputTreenameDeform_, level_, m_h1_, m_nBins_, m_rangeHigh_, m_rangeLow_, m_vtkmap_, TFileDirectory::make(), mid_, TFileService::mkdir(), mlevel_, moduleListName_, phiVal_, rVal_, setCommonTrackerSystem_, AlCaHLTBitMon_QueryRunRegistry::string, sublevel_, surfDeform_, 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 
98  weightBy_ = cfg.getUntrackedParameter<std::string>("weightBy");
99  setCommonTrackerSystem_ = cfg.getUntrackedParameter<std::string>("setCommonTrackerSystem");
100  detIdFlag_ = cfg.getUntrackedParameter<bool>("detIdFlag");
101  detIdFlagFile_ = cfg.getUntrackedParameter<std::string>("detIdFlagFile");
102  weightById_ = cfg.getUntrackedParameter<bool>("weightById");
103  weightByIdFile_ = cfg.getUntrackedParameter<std::string>("weightByIdFile");
104 
105  // if want to use, make id cut list
106  if (detIdFlag_) {
107  std::ifstream fin;
108  fin.open(detIdFlagFile_.c_str());
109 
110  while (!fin.eof() && fin.good()) {
111  uint32_t id;
112  fin >> id;
113  detIdFlagVector_.push_back(id);
114  }
115  fin.close();
116  }
117 
118  // turn weightByIdFile into weightByIdVector
119  if (weightById_) {
120  std::ifstream inFile;
121  inFile.open(weightByIdFile_.c_str());
122  int ctr = 0;
123  while (!inFile.eof()) {
124  ctr++;
125  unsigned int listId;
126  inFile >> listId;
127  inFile.ignore(256, '\n');
128 
129  weightByIdVector_.push_back(listId);
130  }
131  inFile.close();
132  }
133 
134  //root configuration
135  theFile_ = new TFile(filename_.c_str(), "RECREATE");
136  alignTree_ = new TTree("alignTree",
137  "alignTree"); //,"id:level:mid:mlevel:sublevel:x:y:z:r:phi:a:b:c:dx:dy:dz:dr:dphi:da:db:dc");
138  alignTree_->Branch("id", &id_, "id/I");
139  alignTree_->Branch("badModuleQuality", &badModuleQuality_, "badModuleQuality/I");
140  alignTree_->Branch("inModuleList", &inModuleList_, "inModuleList/I");
141  alignTree_->Branch("level", &level_, "level/I");
142  alignTree_->Branch("mid", &mid_, "mid/I");
143  alignTree_->Branch("mlevel", &mlevel_, "mlevel/I");
144  alignTree_->Branch("sublevel", &sublevel_, "sublevel/I");
145  alignTree_->Branch("x", &xVal_, "x/F");
146  alignTree_->Branch("y", &yVal_, "y/F");
147  alignTree_->Branch("z", &zVal_, "z/F");
148  alignTree_->Branch("r", &rVal_, "r/F");
149  alignTree_->Branch("phi", &phiVal_, "phi/F");
150  alignTree_->Branch("eta", &etaVal_, "eta/F");
151  alignTree_->Branch("alpha", &alphaVal_, "alpha/F");
152  alignTree_->Branch("beta", &betaVal_, "beta/F");
153  alignTree_->Branch("gamma", &gammaVal_, "gamma/F");
154  alignTree_->Branch("dx", &dxVal_, "dx/F");
155  alignTree_->Branch("dy", &dyVal_, "dy/F");
156  alignTree_->Branch("dz", &dzVal_, "dz/F");
157  alignTree_->Branch("dr", &drVal_, "dr/F");
158  alignTree_->Branch("dphi", &dphiVal_, "dphi/F");
159  alignTree_->Branch("dalpha", &dalphaVal_, "dalpha/F");
160  alignTree_->Branch("dbeta", &dbetaVal_, "dbeta/F");
161  alignTree_->Branch("dgamma", &dgammaVal_, "dgamma/F");
162  alignTree_->Branch("du", &duVal_, "du/F");
163  alignTree_->Branch("dv", &dvVal_, "dv/F");
164  alignTree_->Branch("dw", &dwVal_, "dw/F");
165  alignTree_->Branch("da", &daVal_, "da/F");
166  alignTree_->Branch("db", &dbVal_, "db/F");
167  alignTree_->Branch("dg", &dgVal_, "dg/F");
168  alignTree_->Branch("useDetId", &useDetId_, "useDetId/I");
169  alignTree_->Branch("detDim", &detDim_, "detDim/I");
170  alignTree_->Branch("surW", &surWidth_, "surW/F");
171  alignTree_->Branch("surL", &surLength_, "surL/F");
172  alignTree_->Branch("surRot", &surRot_, "surRot[9]/D");
173  alignTree_->Branch("identifiers", &identifiers_, "identifiers[6]/I");
174  alignTree_->Branch("type", &type_, "type/I");
175  alignTree_->Branch("surfDeform", &surfDeform_, "surfDeform[13]/D");
176 
177  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
178  it->setPalette(1);
179  it->addPixel(true);
180  }
181 
182  edm::Service<TFileService> fs;
183  TFileDirectory subDir_All = fs->mkdir("AllSubdetectors");
184  TFileDirectory subDir_PXB = fs->mkdir("PixelBarrel");
185  TFileDirectory subDir_PXF = fs->mkdir("PixelEndcap");
186  for (int ii = 0; ii < 13; ++ii) {
187  std::stringstream histname0;
188  histname0 << "SurfDeform_Par_" << ii;
189  m_h1_[histname0.str()] = subDir_All.make<TH1D>(
190  (histname0.str()).c_str(), (histname0.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
191 
192  std::stringstream histname1;
193  histname1 << "SurfDeform_PixelBarrel_Par_" << ii;
194  m_h1_[histname1.str()] = subDir_PXB.make<TH1D>(
195  (histname1.str()).c_str(), (histname1.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
196 
197  std::stringstream histname2;
198  histname2 << "SurfDeform_PixelEndcap_Par_" << ii;
199  m_h1_[histname2.str()] = subDir_PXF.make<TH1D>(
200  (histname2.str()).c_str(), (histname2.str()).c_str(), m_nBins_, m_rangeLow_, m_rangeHigh_);
201  }
202 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
TrackerGeometryCompare::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: TrackerGeometryCompare.h:91
gpuClustering::id
uint16_t *__restrict__ id
Definition: gpuClusterChargeCut.h:20
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:161
TrackerGeometryCompare::geomDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > geomDetToken_
Definition: TrackerGeometryCompare.h:90
cuy.ii
int ii
Definition: cuy.py:589
TrackerGeometryCompare::dummyTracker
AlignableTracker * dummyTracker
Definition: TrackerGeometryCompare.h:97
TrackerGeometryCompare::m_h1_
std::map< std::string, TH1D * > m_h1_
Definition: TrackerGeometryCompare.h:163
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
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
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
TrackerGeometryCompare::currentTracker
AlignableTracker * currentTracker
Definition: TrackerGeometryCompare.h:98
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
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

void TrackerGeometryCompare::addSurveyInfo ( Alignable *  ali)
private

Definition at line 956 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, SurveyError::matrix(), SurveyError::rawId(), DetId::rawId(), makeMuonMisalignmentScenario::rot, AlignableSurface::setLength(), Alignable::setSurvey(), AlignableSurface::setWidth(), SurveyError::structureType(), Alignable::surface(), theSurveyErrors, theSurveyIndex, theSurveyValues, and AlignableSurface::width().

956  {
957  const auto& comp = ali->components();
958 
959  unsigned int nComp = comp.size();
960 
961  for (unsigned int i = 0; i < nComp; ++i)
962  addSurveyInfo(comp[i]);
963 
964  const SurveyError& error = theSurveyErrors->m_surveyErrors[theSurveyIndex];
965 
966  if (ali->geomDetId().rawId() != error.rawId() || ali->alignableObjectId() != error.structureType()) {
967  throw cms::Exception("DatabaseError") << "Error reading survey info from DB. Mismatched id!";
968  }
969 
970  const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[theSurveyIndex].translation();
971  const CLHEP::HepRotation& rot = theSurveyValues->m_align[theSurveyIndex].rotation();
972 
973  AlignableSurface surf(
974  align::PositionType(pos.x(), pos.y(), pos.z()),
975  align::RotationType(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
976 
977  surf.setWidth(ali->surface().width());
978  surf.setLength(ali->surface().length());
979 
980  ali->setSurvey(new SurveyDet(surf, error.matrix()));
981 
982  ++theSurveyIndex;
983 }
AlignableSurface::width
align::Scalar width() const
Definition: AlignableSurface.h:30
mps_fire.i
i
Definition: mps_fire.py:428
SurveyError::matrix
ErrorMatrix matrix() const
Definition: SurveyError.h:63
TrackerGeometryCompare::theSurveyValues
const Alignments * theSurveyValues
Definition: TrackerGeometryCompare.h:101
TrackerGeometryCompare::addSurveyInfo
void addSurveyInfo(Alignable *ali)
Definition: TrackerGeometryCompare.cc:956
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:249
SurveyError::structureType
uint8_t structureType() const
Definition: SurveyError.h:59
Alignments::m_align
std::vector< AlignTransform > m_align
Definition: Alignments.h:19
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
SurveyError::rawId
align::ID rawId() const
Definition: SurveyError.h:61
Alignable::components
virtual const Alignables & components() const =0
Return vector of all direct components.
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:132
AlignableSurface::length
align::Scalar length() const
Definition: AlignableSurface.h:32
Alignable::setSurvey
void setSurvey(const SurveyDet *)
Set survey info.
Definition: Alignable.cc:266
TrackerGeometryCompare::theSurveyErrors
const SurveyErrors * theSurveyErrors
Definition: TrackerGeometryCompare.h:102
SurveyErrors::m_surveyErrors
std::vector< SurveyError > m_surveyErrors
Definition: SurveyErrors.h:22
Alignable::geomDetId
const DetId & geomDetId() const
Definition: Alignable.h:177
void TrackerGeometryCompare::analyze ( const edm::Event &  ,
const edm::EventSetup &  iSetup 
)
overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 224 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(), testEve_cfg::level, levelStrings_, m_theLevels, AlignableTracker::objectIdProvider(), referenceTracker, setCommonTrackerSystem(), setCommonTrackerSystem_, AlignableObjectId::stringToId(), topoToken_, cond::service::PoolDBOutputService::writeOneIOV(), and writeToDB_.

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

Read from DB and print survey info.

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 204 of file TrackerGeometryCompare.cc.

References firstEvent_.

204 { firstEvent_ = true; }
void TrackerGeometryCompare::compareGeometries ( Alignable *  refAli,
Alignable *  curAli,
const TrackerTopology *  tTopo,
const edm::EventSetup &  iSetup 
)
private

Definition at line 603 of file TrackerGeometryCompare.cc.

References Alignable::alignableObjectId(), trackerTree::check(), Alignable::components(), change_name::diff, align::diffAlignables(), alignCSCRings::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().

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

Definition at line 451 of file TrackerGeometryCompare.cc.

References Abs(), cuy::ii, inputFilename1_, inputFilename2_, inputRootFile1_, inputRootFile2_, inputTreenameDeform_, m_h1_, m_nBins_, m_rangeHigh_, m_rangeLow_, m_vtkmap_, and surfDeform_.

Referenced by analyze().

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

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

Referenced by analyze().

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

Definition at line 759 of file TrackerGeometryCompare.cc.

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

Referenced by setCommonTrackerSystem().

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

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 206 of file TrackerGeometryCompare.cc.

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

Referenced by o2olib.O2ORunMgr::executeJob().

206  {
207  int iname(0);
208  for (std::vector<TrackerMap>::iterator it = m_vtkmap_.begin(); it != m_vtkmap_.end(); ++it) {
209  std::stringstream mapname;
210  mapname << "TkMap_SurfDeform" << iname << ".png";
211  it->save(true, 0, 0, mapname.str());
212  mapname.str(std::string());
213  mapname.clear();
214  mapname << "TkMap_SurfDeform" << iname << ".pdf";
215  it->save(true, 0, 0, mapname.str());
216  ++iname;
217  }
218 
219  theFile_->cd();
220  alignTree_->Write();
221  theFile_->Close();
222 }
TrackerGeometryCompare::m_vtkmap_
std::vector< TrackerMap > m_vtkmap_
Definition: TrackerGeometryCompare.h:161
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
void TrackerGeometryCompare::fillIdentifiers ( int  subdetlevel,
int  rawid,
const TrackerTopology *  tTopo 
)
private

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

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

Definition at line 794 of file TrackerGeometryCompare.cc.

References align::AlignableDet, align::AlignableDetUnit, Alignable::alignableObjectId(), alignTree_, alphaVal_, badModuleQuality_, betaVal_, dalphaVal_, daVal_, dbetaVal_, dbVal_, detDim_, detIdFlag_, dgammaVal_, dgVal_, 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_, TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), PV3DBase< T, PVType, FrameType >::y(), yVal_, TkRotation< T >::yx(), TkRotation< T >::yy(), TkRotation< T >::yz(), PV3DBase< T, PVType, FrameType >::z(), zVal_, TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

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

797  {
798  //Get bad modules
799  const SiPixelQuality* SiPixelModules = &iSetup.getData(pixQualityToken_);
800  const SiStripQuality* SiStripModules = &iSetup.getData(stripQualityToken_);
801 
802  id_ = refAli->id();
803 
804  badModuleQuality_ = 0;
805  //check if module has a bad quality tag
806  if (SiPixelModules->IsModuleBad(id_)) {
807  badModuleQuality_ = 1;
808  }
809  if (SiStripModules->IsModuleBad(id_)) {
810  badModuleQuality_ = 1;
811  }
812 
813  //check if module is in a given list of bad/untouched etc. modules
814  inModuleList_ = 0;
815  for (unsigned int i = 0; i < moduleList_.size(); i++) {
816  if (moduleList_[i] == id_) {
817  inModuleList_ = 1;
818  break;
819  }
820  }
821 
822  level_ = refAli->alignableObjectId();
823  //need if ali has no mother
824  if (refAli->mother()) {
825  mid_ = refAli->mother()->geomDetId().rawId();
826  mlevel_ = refAli->mother()->alignableObjectId();
827  } else {
828  mid_ = -1;
829  mlevel_ = -1;
830  }
831  DetId detid(id_);
832  sublevel_ = detid.subdetId();
833  fillIdentifiers(sublevel_, id_, tTopo);
834  xVal_ = refAli->globalPosition().x();
835  yVal_ = refAli->globalPosition().y();
836  zVal_ = refAli->globalPosition().z();
837  align::GlobalVector vec(xVal_, yVal_, zVal_);
838  rVal_ = vec.perp();
839  phiVal_ = vec.phi();
840  etaVal_ = vec.eta();
841  align::RotationType rot = refAli->globalRotation();
842  align::EulerAngles eulerAngles = align::toAngles(rot);
843  alphaVal_ = eulerAngles[0];
844  betaVal_ = eulerAngles[1];
845  gammaVal_ = eulerAngles[2];
846  // global
847  dxVal_ = diff[0];
848  dyVal_ = diff[1];
849  dzVal_ = diff[2];
850  // local
851  duVal_ = diff[6];
852  dvVal_ = diff[7];
853  dwVal_ = diff[8];
854  //...TODO...
855  align::GlobalVector g(dxVal_, dyVal_, dzVal_);
856  //getting dR and dPhi
857  align::GlobalVector vRef(xVal_, yVal_, zVal_);
858  align::GlobalVector vCur(xVal_ + dxVal_, yVal_ + dyVal_, zVal_ + dzVal_);
859  drVal_ = vCur.perp() - vRef.perp();
860  dphiVal_ = vCur.phi() - vRef.phi();
861  // global
862  dalphaVal_ = diff[3];
863  dbetaVal_ = diff[4];
864  dgammaVal_ = diff[5];
865  // local
866  daVal_ = diff[9];
867  dbVal_ = diff[10];
868  dgVal_ = diff[11];
869 
870  //detIdFlag
871  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
872  if (detIdFlag_) {
873  if ((passIdCut(refAli->id())) || (passIdCut(refAli->mother()->id()))) {
874  useDetId_ = 1;
875  } else {
876  useDetId_ = 0;
877  }
878  }
879  }
880  // det module dimension
881  if (refAli->alignableObjectId() == align::AlignableDetUnit) {
882  if (refAli->mother()->alignableObjectId() != align::AlignableDet)
883  detDim_ = 1;
884  else if (refAli->mother()->alignableObjectId() == align::AlignableDet)
885  detDim_ = 2;
886  } else
887  detDim_ = 0;
888 
889  surWidth_ = refAli->surface().width();
890  surLength_ = refAli->surface().length();
891  align::RotationType rt = refAli->globalRotation();
892  surRot_[0] = rt.xx();
893  surRot_[1] = rt.xy();
894  surRot_[2] = rt.xz();
895  surRot_[3] = rt.yx();
896  surRot_[4] = rt.yy();
897  surRot_[5] = rt.yz();
898  surRot_[6] = rt.zx();
899  surRot_[7] = rt.zy();
900  surRot_[8] = rt.zz();
901 
902  //Fill
903  alignTree_->Fill();
904 }
TkRotation::xx
T xx() const
Definition: extTkRotation.h:193
AlignableSurface::width
align::Scalar width() const
Definition: AlignableSurface.h:30
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
mps_fire.i
i
Definition: mps_fire.py:428
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
TkRotation::yx
T yx() const
Definition: extTkRotation.h:196
Alignable::globalRotation
const RotationType & globalRotation() const
Return the global orientation of the object.
Definition: Alignable.h:138
TrackerGeometryCompare::moduleList_
std::vector< int > moduleList_
Definition: TrackerGeometryCompare.h:128
hcal_runs.rt
list rt
Definition: hcal_runs.py:76
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
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:128
TkRotation::zx
T zx() const
Definition: extTkRotation.h:199
TkRotation::xy
T xy() const
Definition: extTkRotation.h:194
TkRotation::zz
T zz() const
Definition: extTkRotation.h:201
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
change_name.diff
list diff
Definition: change_name.py:13
Alignable::alignableObjectId
virtual StructureType alignableObjectId() const =0
Return the alignable type identifier.
TkRotation::zy
T zy() const
Definition: extTkRotation.h:200
SiStripQuality::IsModuleBad
bool IsModuleBad(const uint32_t &detid) const
Definition: SiStripQuality.cc:558
align::toAngles
EulerAngles toAngles(const RotationType &)
Convert rotation matrix to angles about x-, y-, z-axes (frame rotation).
Definition: Utilities.cc:8
TkRotation::yy
T yy() const
Definition: extTkRotation.h:197
SiPixelQuality::IsModuleBad
bool IsModuleBad(const uint32_t &detid) const
Definition: SiPixelQuality.cc:93
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:132
DetId
Definition: DetId.h:17
align::EulerAngles
AlgebraicVector EulerAngles
Definition: Definitions.h:34
AlignableSurface::length
align::Scalar length() const
Definition: AlignableSurface.h:32
TrackerGeometryCompare::fillIdentifiers
void fillIdentifiers(int subdetlevel, int rawid, const TrackerTopology *tTopo)
Definition: TrackerGeometryCompare.cc:997
TkRotation::xz
T xz() const
Definition: extTkRotation.h:195
Alignable::globalPosition
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:135
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
Alignable::mother
Alignable * mother() const
Return pointer to container alignable (if any)
Definition: Alignable.h:91
TkRotation::yz
T yz() const
Definition: extTkRotation.h:198
Alignable::geomDetId
const DetId & geomDetId() const
Definition: Alignable.h:177
TrackerGeometryCompare::pixQualityToken_
const edm::ESGetToken< SiPixelQuality, SiPixelQualityRcd > pixQualityToken_
Definition: TrackerGeometryCompare.h:93
bool TrackerGeometryCompare::passIdCut ( uint32_t  id)
private

Definition at line 985 of file TrackerGeometryCompare.cc.

References detIdFlagVector_, and mps_fire::i.

Referenced by fillTree().

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

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

715  {
716  edm::LogInfo("TrackerGeometryCompare") << "Setting Common Tracker System....";
717 
718  // DM_534??AlignableObjectId dummy;
719  // DM_534??_commonTrackerLevel = dummy.nameToType(_setCommonTrackerSystem);
720  commonTrackerLevel_ = currentTracker->objectIdProvider().stringToId(setCommonTrackerSystem_);
721 
722  diffCommonTrackerSystem(referenceTracker, currentTracker);
723 
724  align::EulerAngles dOmega(3);
725  dOmega[0] = TrackerCommonR_.x();
726  dOmega[1] = TrackerCommonR_.y();
727  dOmega[2] = TrackerCommonR_.z();
728  align::RotationType rot = align::toMatrix(dOmega);
729  align::GlobalVector theR = TrackerCommonT_;
730 
731  edm::LogInfo("TrackerGeometryCompare") << "what we get from overlaying the pixels..." << theR << ", " << rot;
732 
733  //transform to the Tracker System
734  align::PositionType trackerCM = currentTracker->globalPosition();
735  align::GlobalVector cmDiff(
736  trackerCM.x() - TrackerCommonCM_.x(), trackerCM.y() - TrackerCommonCM_.y(), trackerCM.z() - TrackerCommonCM_.z());
737 
738  edm::LogInfo("TrackerGeometryCompare") << "Pixel CM: " << TrackerCommonCM_ << ", tracker CM: " << trackerCM;
739 
740  //adjust translational difference factoring in different rotational CM
741  //needed because rotateInGlobalFrame is about CM of alignable, not Tracker
742  const align::GlobalVector::BasicVectorType& lpvgf = cmDiff.basicVector();
743  align::GlobalVector moveV(rot.multiplyInverse(lpvgf) - lpvgf);
744  align::GlobalVector theRprime(theR + moveV);
745 
746  AlgebraicVector TrackerCommonTR(6);
747  TrackerCommonTR(1) = theRprime.x();
748  TrackerCommonTR(2) = theRprime.y();
749  TrackerCommonTR(3) = theRprime.z();
750  TrackerCommonTR(4) = TrackerCommonR_.x();
751  TrackerCommonTR(5) = TrackerCommonR_.y();
752  TrackerCommonTR(6) = TrackerCommonR_.z();
753 
754  edm::LogInfo("TrackerGeometryCompare") << "and after the transformation: " << TrackerCommonTR;
755 
756  align::moveAlignable(currentTracker, TrackerCommonTR);
757 }
AlignableObjectId::stringToId
align::StructureType stringToId(const char *) const
Definition: AlignableObjectId.cc:284
TrackerGeometryCompare::TrackerCommonR_
align::GlobalVector TrackerCommonR_
Definition: TrackerGeometryCompare.h:124
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
TrackerGeometryCompare::diffCommonTrackerSystem
void diffCommonTrackerSystem(Alignable *refAli, Alignable *curAli)
Definition: TrackerGeometryCompare.cc:759
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
TrackerGeometryCompare::TrackerCommonT_
align::GlobalVector TrackerCommonT_
Definition: TrackerGeometryCompare.h:123
AlignableTracker::objectIdProvider
const AlignableObjectId & objectIdProvider() const
Return tracker alignable object ID provider derived from the tracker&#39;s geometry.
Definition: AlignableTracker.h:110
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
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
Alignable::globalPosition
const PositionType & globalPosition() const
Return the global position of the object.
Definition: Alignable.h:135
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
align::moveAlignable
void moveAlignable(Alignable *ali, AlgebraicVector diff)
Moves the alignable by the AlgebraicVector.
Definition: AlignTools.cc:84
void TrackerGeometryCompare::surveyToTracker ( AlignableTracker *  ali,
Alignments *  alignVals,
AlignmentErrorsExtended *  alignErrors 
)
private

Definition at line 906 of file TrackerGeometryCompare.cc.

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

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

TTree* TrackerGeometryCompare::alignTree_
private

Definition at line 134 of file TrackerGeometryCompare.h.

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

float TrackerGeometryCompare::alphaVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

int TrackerGeometryCompare::badModuleQuality_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::betaVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

align::StructureType TrackerGeometryCompare::commonTrackerLevel_
private

Definition at line 122 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

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

Definition at line 88 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

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

Definition at line 87 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

AlignableTracker* TrackerGeometryCompare::currentTracker
private

Definition at line 98 of file TrackerGeometryCompare.h.

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

float TrackerGeometryCompare::dalphaVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::daVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dbetaVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dbVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

int TrackerGeometryCompare::detDim_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

bool TrackerGeometryCompare::detIdFlag_
private

Definition at line 115 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::detIdFlagFile_
private

Definition at line 116 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

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

Definition at line 121 of file TrackerGeometryCompare.h.

Referenced by passIdCut(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dgammaVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dgVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dphiVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::drVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

AlignableTracker* TrackerGeometryCompare::dummyTracker
private

Definition at line 97 of file TrackerGeometryCompare.h.

float TrackerGeometryCompare::duVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dvVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dwVal_
private

Definition at line 148 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dxVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dyVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::dzVal_
private

Definition at line 146 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::etaVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::filename_
private

Definition at line 132 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

bool TrackerGeometryCompare::firstEvent_
private

Definition at line 159 of file TrackerGeometryCompare.h.

Referenced by analyze(), and beginJob().

bool TrackerGeometryCompare::fromDD4hep_
private

Definition at line 111 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

float TrackerGeometryCompare::gammaVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

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

Definition at line 90 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

int TrackerGeometryCompare::id_
private

\ Tree variables

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

uint32_t TrackerGeometryCompare::identifiers_[6]
private

Definition at line 150 of file TrackerGeometryCompare.h.

Referenced by fillIdentifiers(), and TrackerGeometryCompare().

int TrackerGeometryCompare::inModuleList_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::inputFilename1_
private

Definition at line 107 of file TrackerGeometryCompare.h.

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

std::string TrackerGeometryCompare::inputFilename2_
private

Definition at line 108 of file TrackerGeometryCompare.h.

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

TFile* TrackerGeometryCompare::inputRootFile1_
private

Definition at line 135 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

TFile* TrackerGeometryCompare::inputRootFile2_
private

Definition at line 136 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and createROOTGeometry().

TTree* TrackerGeometryCompare::inputTree01_
private

Definition at line 137 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

TTree* TrackerGeometryCompare::inputTree02_
private

Definition at line 138 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

TTree* TrackerGeometryCompare::inputTree11_
private

Definition at line 139 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

TTree* TrackerGeometryCompare::inputTree12_
private

Definition at line 140 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

std::string TrackerGeometryCompare::inputTreenameAlign_
private

Definition at line 109 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::inputTreenameDeform_
private

Definition at line 110 of file TrackerGeometryCompare.h.

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

int TrackerGeometryCompare::level_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

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

Definition at line 105 of file TrackerGeometryCompare.h.

Referenced by analyze().

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

Definition at line 163 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

int TrackerGeometryCompare::m_nBins_
private

Definition at line 155 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

double TrackerGeometryCompare::m_rangeHigh_
private

Definition at line 157 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

double TrackerGeometryCompare::m_rangeLow_
private

Definition at line 156 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

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

Definition at line 58 of file TrackerGeometryCompare.h.

Referenced by analyze(), and compareGeometries().

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

Definition at line 161 of file TrackerGeometryCompare.h.

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

int TrackerGeometryCompare::mid_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

int TrackerGeometryCompare::mlevel_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

int TrackerGeometryCompare::moduleInList_
private

Definition at line 129 of file TrackerGeometryCompare.h.

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

Definition at line 128 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and fillTree().

std::ifstream TrackerGeometryCompare::moduleListFile_
private

Definition at line 127 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

std::string TrackerGeometryCompare::moduleListName_
private

Definition at line 106 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry(), and TrackerGeometryCompare().

float TrackerGeometryCompare::phiVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

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

Definition at line 93 of file TrackerGeometryCompare.h.

Referenced by fillTree().

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

Definition at line 92 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

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

Definition at line 91 of file TrackerGeometryCompare.h.

Referenced by createROOTGeometry().

AlignableTracker* TrackerGeometryCompare::referenceTracker
private

Definition at line 96 of file TrackerGeometryCompare.h.

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

float TrackerGeometryCompare::rVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::setCommonTrackerSystem_
private

Definition at line 114 of file TrackerGeometryCompare.h.

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

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

Definition at line 94 of file TrackerGeometryCompare.h.

Referenced by fillTree().

int TrackerGeometryCompare::sublevel_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

double TrackerGeometryCompare::surfDeform_[13]
private

Definition at line 153 of file TrackerGeometryCompare.h.

Referenced by compareSurfaceDeformations(), and TrackerGeometryCompare().

float TrackerGeometryCompare::surLength_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

double TrackerGeometryCompare::surRot_[9]
private

Definition at line 151 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::surWidth_
private

Definition at line 149 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

TFile* TrackerGeometryCompare::theFile_
private

Definition at line 133 of file TrackerGeometryCompare.h.

Referenced by endJob(), and TrackerGeometryCompare().

const SurveyErrors* TrackerGeometryCompare::theSurveyErrors
private

Definition at line 102 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

unsigned int TrackerGeometryCompare::theSurveyIndex
private

Definition at line 100 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

const Alignments* TrackerGeometryCompare::theSurveyValues
private

Definition at line 101 of file TrackerGeometryCompare.h.

Referenced by addSurveyInfo().

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

Definition at line 89 of file TrackerGeometryCompare.h.

Referenced by analyze(), and createROOTGeometry().

align::PositionType TrackerGeometryCompare::TrackerCommonCM_
private

Definition at line 125 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

align::GlobalVector TrackerGeometryCompare::TrackerCommonR_
private

Definition at line 124 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

align::GlobalVector TrackerGeometryCompare::TrackerCommonT_
private

Definition at line 123 of file TrackerGeometryCompare.h.

Referenced by diffCommonTrackerSystem(), and setCommonTrackerSystem().

int TrackerGeometryCompare::type_
private

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

int TrackerGeometryCompare::useDetId_
private

Definition at line 143 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

std::string TrackerGeometryCompare::weightBy_
private

Definition at line 113 of file TrackerGeometryCompare.h.

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

bool TrackerGeometryCompare::weightById_
private

Definition at line 117 of file TrackerGeometryCompare.h.

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

std::string TrackerGeometryCompare::weightByIdFile_
private

Definition at line 118 of file TrackerGeometryCompare.h.

Referenced by TrackerGeometryCompare().

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

Definition at line 119 of file TrackerGeometryCompare.h.

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

bool TrackerGeometryCompare::writeToDB_
private

Definition at line 112 of file TrackerGeometryCompare.h.

Referenced by analyze().

float TrackerGeometryCompare::xVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::yVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

float TrackerGeometryCompare::zVal_
private

Definition at line 144 of file TrackerGeometryCompare.h.

Referenced by fillTree(), and TrackerGeometryCompare().

Generated for CMSSW Reference Manual by doxygen 1.8.5