28 for(
int stationCounter = 0; stationCounter < 4; stationCounter++) {
29 for(
int sectorCounter = 0; sectorCounter < 14; sectorCounter++) {
30 if(
chambers[stationCounter][sectorCounter]->getNumberPoints() > 2) {
44 for (
int cont_stat = 0; cont_stat < 4; cont_stat++) {
46 for(
int cont_sect = 0; cont_sect < 14; cont_sect++) {
56 double x,
y,
z,
rms, dx, dy, dz;
57 file >> code >> x >> y >> z >> rms >> dx >> dy >> dz;
59 x = x/10.0; y=y/10.0; z=z/10.0; dx=dx/10.0; dy=dy/10.0; dz=dz/10.0;rms=rms/10.0;
60 station = code/10000 - 1;
61 sector = (code-(station+1)*10000)/100 - 1;
66 disp(0,0) = dx; disp(1,0) = dy; disp(2,0) = dz;
71 GlobalPoint rt(
r(0,0)-disp(0,0),
r(1,0)-disp(1,0),
r(2,0)-disp(2,0));
77 rLocal(0,0) = rl.
x(); rLocal(1,0) = rl.
y(); rLocal(2,0) = rl.
z();
79 rTeo(0,0) = rtl.
x(); rTeo(1,0) = rtl.
y(); rTeo(2,0) = rtl.
z();
80 TMatrixD
diff = rLocal-rTeo;
113 float zOffset, deltax, deltay, deltaz,
alpha,
beta, gamma;
114 wheeltowheel >> zOffset >> deltax >> deltay >> deltaz >> alpha >> beta >> gamma;
115 wheeltowheel.close();
121 delta(0,0) = deltax/10.0;
122 delta(1,0) = deltay/10.0;
123 delta(2,0) = deltaz/10.0;
127 TMatrixD alpha_m(3,3);
128 TMatrixD beta_m(3,3);
129 TMatrixD gamma_m(3,3);
133 for(
int k = 0;
k < 3;
k++) {
141 alpha_m(1,1) =
cos(alpha);
142 alpha_m(1,2) =
sin(alpha);
143 alpha_m(2,1) = -
sin(alpha);
144 alpha_m(2,2) =
cos(alpha);
145 beta_m(0,0) =
cos(beta);
146 beta_m(0,2) = -
sin(beta);
147 beta_m(2,0) =
sin(beta);
148 beta_m(2,2) =
cos(beta);
149 gamma_m(0,0) =
cos(gamma);
150 gamma_m(0,1) =
sin(gamma);
151 gamma_m(1,0) = -
sin(gamma);
152 gamma_m(1,1) =
cos(gamma);
153 Rot = alpha_m*beta_m*gamma_m;
159 for(
int stationCounter = 0; stationCounter < 4; stationCounter++) {
160 for(
int sectorCounter = 0; sectorCounter < 14; sectorCounter++) {
void ReadChambers(edm::ESHandle< DTGeometry >)
Sin< T >::type sin(const T &t)
LocalPoint toLocal(const GlobalPoint &gp) const
Conversion to the R.F. of the GeomDet.
int getNumberPoints() const
std::ostream & operator<<(std::ostream &out, const ALILine &li)
std::string nameOfWheelInfoFile
uint32_t rawId() const
get the raw id
tuple obj
Example code starts here #.
const DTSurveyChamber * getChamber(int, int) const
Cos< T >::type cos(const T &t)
DTSurveyChamber *** chambers
void addPoint(int, const TMatrixD &, const TMatrixD &, const TMatrixD &)
std::string nameOfChamberInfoFile
DTSurvey(const std::string &, const std::string &, int)