|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
|
// -*- C++ -*-
//
// Package: MuonAlignment
// Class : MuonAlignmentInputSurveyDB
//
// Implementation:
// <Notes on implementation>
//
// Original Author: Jim Pivarski
// Created: Thu Mar 6 17:30:46 CST 2008
// $Id: MuonAlignmentInputSurveyDB.cc,v 1.1 2008/03/15 20:26:46 pivarski Exp $
//
// system include files
#include "FWCore/Framework/interface/ESHandle.h"
// user include files
#include "Alignment/MuonAlignment/interface/MuonAlignmentInputSurveyDB.h"
#include "CondFormats/AlignmentRecord/interface/DTSurveyRcd.h"
#include "CondFormats/AlignmentRecord/interface/DTSurveyErrorExtendedRcd.h"
#include "CondFormats/AlignmentRecord/interface/CSCSurveyRcd.h"
#include "CondFormats/AlignmentRecord/interface/CSCSurveyErrorExtendedRcd.h"
#include "Alignment/CommonAlignment/interface/SurveyDet.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
//
// constants, enums and typedefs
//
//
// static data member definitions
//
//
// constructors and destructor
//
MuonAlignmentInputSurveyDB::MuonAlignmentInputSurveyDB(const DTGeometry* dtGeometry,
const CSCGeometry* cscGeometry,
const GEMGeometry* gemGeometry,
const Alignments* dtSurvey,
const Alignments* cscSurvey,
const SurveyErrors* dtSurveyError,
const SurveyErrors* cscSurveyError)
: dtGeometry_(dtGeometry),
cscGeometry_(cscGeometry),
gemGeometry_(gemGeometry),
dtSurvey_(dtSurvey),
cscSurvey_(cscSurvey),
dtSurveyError_(dtSurveyError),
cscSurveyError_(cscSurveyError) {}
// MuonAlignmentInputSurveyDB::MuonAlignmentInputSurveyDB(const MuonAlignmentInputSurveyDB& rhs)
// {
// // do actual copying here;
//
MuonAlignmentInputSurveyDB::~MuonAlignmentInputSurveyDB() {}
//
// assignment operators
//
// const MuonAlignmentInputSurveyDB& MuonAlignmentInputSurveyDB::operator=(const MuonAlignmentInputSurveyDB& rhs)
// {
// //An exception safe implementation is
// MuonAlignmentInputSurveyDB temp(rhs);
// swap(rhs);
//
// return *this;
// }
//
// member functions
//
AlignableMuon* MuonAlignmentInputSurveyDB::newAlignableMuon() const {
AlignableMuon* output = new AlignableMuon(dtGeometry_, cscGeometry_, gemGeometry_);
unsigned int theSurveyIndex = 0;
const Alignments* theSurveyValues = dtSurvey_;
const SurveyErrors* theSurveyErrors = dtSurveyError_;
const auto& barrels = output->DTBarrel();
for (const auto& iter : barrels) {
addSurveyInfo_(iter, &theSurveyIndex, theSurveyValues, theSurveyErrors);
}
theSurveyIndex = 0;
theSurveyValues = cscSurvey_;
theSurveyErrors = cscSurveyError_;
const auto& endcaps = output->CSCEndcaps();
for (const auto& iter : endcaps) {
addSurveyInfo_(iter, &theSurveyIndex, theSurveyValues, theSurveyErrors);
}
return output;
}
// This function was copied (with minimal modifications) from
// Alignment/CommonAlignmentProducer/plugins/AlignmentProducer.cc
// (version CMSSW_2_0_0_pre1), guaranteed to work the same way
// unless AlignmentProducer.cc's version changes!
void MuonAlignmentInputSurveyDB::addSurveyInfo_(Alignable* ali,
unsigned int* theSurveyIndex,
const Alignments* theSurveyValues,
const SurveyErrors* theSurveyErrors) const {
const auto& comp = ali->components();
unsigned int nComp = comp.size();
for (unsigned int i = 0; i < nComp; ++i)
addSurveyInfo_(comp[i], theSurveyIndex, theSurveyValues, theSurveyErrors);
const SurveyError& error = theSurveyErrors->m_surveyErrors[*theSurveyIndex];
if (ali->geomDetId().rawId() != error.rawId() || ali->alignableObjectId() != error.structureType()) {
throw cms::Exception("DatabaseError") << "Error reading survey info from DB. Mismatched id!";
}
const CLHEP::Hep3Vector& pos = theSurveyValues->m_align[*theSurveyIndex].translation();
const CLHEP::HepRotation& rot = theSurveyValues->m_align[*theSurveyIndex].rotation();
AlignableSurface surf(
align::PositionType(pos.x(), pos.y(), pos.z()),
align::RotationType(rot.xx(), rot.xy(), rot.xz(), rot.yx(), rot.yy(), rot.yz(), rot.zx(), rot.zy(), rot.zz()));
surf.setWidth(ali->surface().width());
surf.setLength(ali->surface().length());
ali->setSurvey(new SurveyDet(surf, error.matrix()));
(*theSurveyIndex)++;
}
//
// const member functions
//
//
// static member functions
//
|