File indexing completed on 2023-03-17 10:39:53
0001 #include "Alignment/MuonAlignmentAlgorithms/interface/ReadPGInfo.h"
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047
0048
0049
0050
0051
0052
0053
0054
0055
0056
0057
0058
0059
0060
0061
0062
0063
0064
0065
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079
0080
0081
0082
0083
0084
0085
0086
0087
0088
0089
0090
0091
0092
0093
0094
0095
0096
0097
0098
0099
0100
0101
0102
0103
0104
0105
0106
0107
0108
0109
0110
0111
0112
0113
0114
0115
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0126
0127
0128
0129
0130
0131
0132
0133
0134
0135
0136
0137
0138
0139
0140
0141
0142
0143
0144
0145
0146
0147
0148
0149
0150
0151
0152
0153
0154
0155
0156
0157
0158
0159
0160
0161
0162
0163
0164
0165
0166
0167
0168
0169
0170
0171
0172
0173
0174
0175
0176
0177
0178
0179
0180
0181
0182
0183
0184
0185
0186
0187
0188
0189
0190
0191
0192
0193
0194
0195
0196
0197
0198
0199
0200
0201
0202
0203
0204
0205
0206
0207
0208
0209
0210
0211
0212
0213
0214
0215
0216
0217
0218
0219
0220
0221
0222
0223
0224
0225
0226
0227
0228
0229
0230
0231
0232
0233
0234
0235
0236
0237
0238
0239
0240
0241
0242
0243 static char chambers[TOTALCHAMBERS][20] = {
0244 "351010100100009", "351010100100016", "351010100100014", "351010100100017", "351010100100011", "351010100100038",
0245 "351010100100035", "351010100100023", "351010100100036", "351010100100072", "351010100200013", "351010100200015",
0246 "351010100200020", "351010100200019", "351010100200010", "351010100200025", "351010100200037", "351010100200028",
0247 "351010100200034", "351010100200022", "351010100200070", "351010100300002", "351010100300061", "351030100500025",
0248 "351030100500019", "351030100500017", "351030100500021", "351030100500023", "351030100500038", "351030100500039",
0249 "351030100500047", "351030100500040", "351030100500001", "351030100600020", "351030100600026", "351030100600006",
0250 "351030100600024", "351030100600022", "351030100600018", "351030100600036", "351030100600035", "351030100600049",
0251 "351030100600046", "351030100600048", "351030100700014", "351030100710065", "351030100800009", "351020100900010",
0252 "351020100900014", "351020100900028", "351020100900030", "351020100900024", "351020100900022", "351020100900034",
0253 "351020100900040", "351020100900046", "351020101000018", "351020101000016", "351020101000008", "351020101000006",
0254 "351020101000012", "351020101000004", "351020101000032", "351020101000020", "351020101000026", "351020101000042",
0255 "351020101200001", "351010101300051", "351030101700054", "351020102100002", "351040102500003", "351040102500006",
0256 "351040102600004", "351040102600010", "351040102600030", "351040102700002", "351040102700005", "351040102700031",
0257 "351040102800001", "351040102800007", "351020103000004", "351020103000003", "351020103000009", "351040103100004",
0258 "351040103200001", "351040103300003", "351040103400002", "351010103600075", "351030103700033", "351030103800034",
0259 "351030103900031", "351030104000030", "351020104200001", "351020104200002", "351010100100018", "351010100100012",
0260 "351010100100007", "351010100100021", "351010100100067", "351010100200008", "351010100200053", "351010100200044",
0261 "351010100200069", "351010100200071", "351010100300005", "351010100300001", "351010100300024", "351010100300039",
0262 "351010100300045", "351010100300042", "351010100300056", "351010100300058", "351010100300059", "351010100300062",
0263 "351010100300065", "351010100300068", "351010100300073", "351010100300057", "351010100400003", "351010100400004",
0264 "351010100400006", "351010100400027", "351010100400026", "351010100400043", "351010100400041", "351010100400040",
0265 "351010100400055", "351010100400064", "351010100400060", "351010100400063", "351010100400054", "351010100400066",
0266 "351010100400074", "351030100500007", "351030100500005", "351030100500002", "351030100500004", "351030100500045",
0267 "351030100600008", "351030100600037", "351030100600003", "351030100600064", "351030100600063", "351030100700016",
0268 "351030100700028", "351030100700013", "351030100700010", "351030100700011", "351030100700042", "351030100700043",
0269 "351030100700050", "351030100700052", "351030100700061", "351030100700067", "351030100700068", "351030100700069",
0270 "351030100700065", "351030100800015", "351030100800027", "351030100800012", "351030100800041", "351030100800044",
0271 "351030100800051", "351030100800059", "351030100800060", "351030100800062", "351030100800066", "351030100800072",
0272 "351030100800073", "351030100800071", "351030100800074", "351030100800070", "351020100900002", "351020100900038",
0273 "351020100900048", "351020100900052", "351020100900060", "351020101000044", "351020101000036", "351020101000050",
0274 "351020101000056", "351020101000058", "351020101100005", "351020101100007", "351020101100011", "351020101100017",
0275 "351020101100013", "351020101100025", "351020101100027", "351020101100031", "351020101100035", "351020101100041",
0276 "351020101100045", "351020101100047", "351020101100051", "351020101100057", "351020101200009", "351020101200003",
0277 "351020101200021", "351020101200015", "351020101200023", "351020101200019", "351020101200033", "351020101200029",
0278 "351020101200037", "351020101200039", "351020101200049", "351020101200053", "351020101200055", "351020101200043",
0279 "351020101200059", "351010101500052", "351030101900053", "351020102300001", "351040102500016", "351040102500023",
0280 "351040102500021", "351040102500025", "351040102600008", "351040102600011", "351040102600014", "351040102600018",
0281 "351040102600019", "351040102600028", "351040102700013", "351040102700015", "351040102700026", "351040102700020",
0282 "351040102700027", "351040102800009", "351040102800012", "351040102800017", "351040102800024", "351040102800022",
0283 "351040102800029", "351020102900002", "351020102900001", "351020102900005", "351020102900006", "351020103000008",
0284 "351020103000007", "351040103100006", "351040103100010", "351040103200008", "351040103300005", "351040103300009",
0285 "351040103400007", "351010103500030", "351010103500029", "351010103500048", "351010103500050", "351010103500047",
0286 "351010103600031", "351010103600032", "351010103600033", "351010103600049", "351010103600046", "351030103700058",
0287 "351030103800032", "351030103800057", "351030103900056", "351030104000029", "351030104000055", "351040104100001"};
0288
0289 static int position[TOTALCHAMBERS][3] = {
0290 {2, 3, 1}, {2, 5, 1}, {2, 9, 1}, {2, 11, 1}, {1, 2, 1}, {1, 6, 1}, {1, 12, 1}, {1, 10, 1}, {1, 8, 1},
0291 {0, 0, 0}, {2, 6, 1}, {2, 10, 1}, {2, 12, 1}, {1, 5, 1}, {2, 8, 1}, {1, 11, 1}, {1, 3, 1}, {2, 2, 1},
0292 {2, 4, 1}, {1, 9, 1}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {1, 8, 2}, {2, 5, 2}, {2, 11, 2}, {1, 12, 2},
0293 {2, 9, 2}, {1, 6, 2}, {1, 2, 2}, {2, 3, 2}, {1, 10, 2}, {0, 0, 0}, {1, 9, 2}, {2, 12, 2}, {0, 0, 0},
0294 {2, 10, 2}, {2, 8, 2}, {1, 11, 2}, {1, 5, 2}, {1, 3, 2}, {2, 6, 2}, {2, 2, 2}, {2, 4, 2}, {0, 0, 0},
0295 {0, 0, 0}, {0, 0, 0}, {1, 10, 3}, {2, 5, 3}, {1, 6, 3}, {2, 3, 3}, {2, 11, 3}, {1, 2, 3}, {2, 9, 3},
0296 {1, 8, 3}, {1, 12, 3}, {2, 2, 3}, {2, 12, 3}, {2, 8, 3}, {1, 5, 3}, {2, 10, 3}, {2, 6, 3}, {2, 4, 3},
0297 {1, 11, 3}, {1, 9, 3}, {1, 3, 3}, {0, 0, 0}, {1, 4, 1}, {1, 4, 2}, {1, 4, 3}, {2, 5, 4}, {1, 6, 4},
0298 {2, 6, 4}, {1, 5, 4}, {0, 0, 0}, {2, 3, 4}, {1, 2, 4}, {0, 0, 0}, {2, 2, 4}, {1, 3, 4}, {2, 4, 4},
0299 {2, 13, 4}, {0, 0, 0}, {1, 8, 4}, {2, 8, 4}, {1, 12, 4}, {2, 12, 4}, {0, 0, 0}, {1, 10, 4}, {2, 10, 4},
0300 {1, 14, 4}, {2, 14, 4}, {1, 13, 4}, {1, 4, 4}, {2, 7, 1}, {0, 11, 1}, {2, 1, 1}, {0, 3, 1}, {0, 7, 1},
0301 {0, 2, 1}, {0, 10, 1}, {0, 6, 1}, {1, 1, 1}, {1, 7, 1}, {0, 4, 1}, {-2, 8, 1}, {0, 8, 1}, {0, 12, 1},
0302 {-1, 11, 1}, {-2, 10, 1}, {-1, 9, 1}, {-2, 2, 1}, {-1, 5, 1}, {-2, 4, 1}, {-2, 6, 1}, {-1, 1, 1}, {-1, 7, 1},
0303 {-2, 12, 1}, {-1, 6, 1}, {-1, 4, 1}, {-1, 2, 1}, {-1, 8, 1}, {-2, 9, 1}, {-1, 10, 1}, {-2, 11, 1}, {-1, 12, 1},
0304 {0, 9, 1}, {0, 5, 1}, {-2, 7, 1}, {-2, 5, 1}, {-2, 3, 1}, {0, 1, 1}, {-2, 1, 1}, {0, 3, 2}, {0, 7, 2},
0305 {2, 7, 2}, {2, 1, 2}, {0, 11, 2}, {0, 2, 2}, {0, 10, 2}, {0, 6, 2}, {1, 7, 2}, {1, 1, 2}, {0, 4, 2},
0306 {-1, 11, 2}, {-2, 10, 2}, {-1, 5, 2}, {-2, 8, 2}, {0, 8, 2}, {0, 12, 2}, {-1, 9, 2}, {-2, 12, 2}, {-2, 2, 2},
0307 {-2, 6, 2}, {-1, 1, 2}, {-1, 7, 2}, {-2, 4, 2}, {-1, 10, 2}, {-2, 11, 2}, {-1, 2, 2}, {0, 5, 2}, {0, 9, 2},
0308 {-1, 8, 2}, {-1, 4, 2}, {-1, 12, 2}, {-1, 6, 2}, {-2, 5, 2}, {-2, 1, 2}, {-2, 7, 2}, {-2, 3, 2}, {0, 1, 2},
0309 {-2, 9, 2}, {2, 1, 3}, {0, 3, 3}, {0, 11, 3}, {2, 7, 3}, {0, 7, 3}, {0, 10, 3}, {0, 2, 3}, {0, 6, 3},
0310 {1, 1, 3}, {1, 7, 3}, {-1, 5, 3}, {-2, 8, 3}, {-1, 9, 3}, {0, 4, 3}, {-1, 1, 3}, {0, 8, 3}, {-2, 10, 3},
0311 {-1, 11, 3}, {0, 12, 3}, {-2, 12, 3}, {-2, 2, 3}, {-2, 4, 3}, {-2, 6, 3}, {-1, 7, 3}, {-1, 8, 3}, {-1, 4, 3},
0312 {-1, 2, 3}, {-1, 6, 3}, {0, 5, 3}, {-1, 12, 3}, {-1, 10, 3}, {-2, 11, 3}, {0, 9, 3}, {-2, 9, 3}, {-2, 5, 3},
0313 {0, 1, 3}, {-2, 1, 3}, {-2, 3, 3}, {-2, 7, 3}, {-1, 3, 1}, {-1, 3, 2}, {-1, 3, 3}, {2, 7, 4}, {-2, 2, 4},
0314 {0, 7, 4}, {-1, 1, 4}, {-2, 3, 4}, {0, 6, 4}, {-1, 2, 4}, {1, 7, 4}, {0, 1, 4}, {-2, 1, 4}, {0, 3, 4},
0315 {2, 1, 4}, {-2, 6, 4}, {-1, 5, 4}, {-1, 7, 4}, {0, 2, 4}, {0, 5, 4}, {1, 1, 4}, {-2, 5, 4}, {-1, 6, 4},
0316 {-2, 7, 4}, {0, 13, 4}, {0, 4, 4}, {-2, 13, 4}, {-2, 4, 4}, {-1, 13, 4}, {-1, 4, 4}, {0, 12, 4}, {-2, 12, 4},
0317 {-1, 12, 4}, {0, 8, 4}, {-2, 8, 4}, {-1, 8, 4}, {2, 11, 4}, {2, 9, 4}, {-1, 11, 4}, {0, 11, 4}, {-1, 9, 4},
0318 {-2, 9, 4}, {1, 9, 4}, {1, 11, 4}, {-2, 11, 4}, {0, 9, 4}, {-2, 14, 4}, {0, 10, 4}, {-1, 14, 4}, {-2, 10, 4},
0319 {0, 14, 4}, {-1, 10, 4}, {-1, 3, 4}};
0320
0321 ReadPGInfo::ReadPGInfo(const char *name) {
0322 TDirectory *dirSave = gDirectory;
0323 rootFile = new TFile(name);
0324 dirSave->cd();
0325 }
0326
0327 ReadPGInfo::~ReadPGInfo() { delete rootFile; }
0328
0329 char *ReadPGInfo::getId(int wheel, int station, int sector) {
0330 for (int counter = 0; counter < TOTALCHAMBERS; ++counter) {
0331 if (wheel == position[counter][0] && sector == position[counter][1] && station == position[counter][2])
0332 return chambers[counter];
0333 }
0334 return nullptr;
0335 }
0336
0337 TMatrixD ReadPGInfo::giveR(int wheel, int station, int sector) {
0338 TMatrixD *empty = new TMatrixD(0, 0);
0339 char *id = getId(wheel, station, sector);
0340 if (id == nullptr)
0341 return *empty;
0342 TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
0343 TDirectoryFile *myR = (TDirectoryFile *)myDir->Get("R");
0344 TMatrixD *R = (TMatrixD *)myR->Get("matrix");
0345 return *R;
0346 }
0347
0348 TMatrixD ReadPGInfo::giveQCCal(int wheel, int station, int sector) {
0349 TMatrixD *mat = new TMatrixD(0, 0);
0350 TMatrixD qc = giveQC(wheel, station, sector);
0351 if (qc.GetNrows() == 0)
0352 return *mat;
0353 mat->ResizeTo(12, 2);
0354 int maxCount = 12;
0355 if (station == 4)
0356 maxCount = 8;
0357 for (int c = 0; c < maxCount; ++c) {
0358 float error;
0359 if (qc(c, 1) == 0 || qc(c, 3) == 0) {
0360 (*mat)(c, 0) = (qc(c, 0) + qc(c, 2)) / 2.0;
0361 (*mat)(c, 1) = 500;
0362 } else {
0363 error = 1.0 / (1.0 / (qc(c, 1) * qc(c, 1)) + 1.0 / (qc(c, 3) * qc(c, 3)));
0364 (*mat)(c, 0) = (qc(c, 0) / (qc(c, 1) * qc(c, 1)) + qc(c, 2) / (qc(c, 3) * qc(c, 3))) * error;
0365 (*mat)(c, 1) = TMath::Sqrt(error);
0366 }
0367 }
0368 return *mat;
0369 }
0370
0371 TMatrixD ReadPGInfo::giveQC(int wheel, int station, int sector) {
0372 TMatrixD *empty = new TMatrixD(0, 0);
0373 char *id = getId(wheel, station, sector);
0374 if (id == nullptr)
0375 return *empty;
0376 TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
0377 TDirectoryFile *myQC = (TDirectoryFile *)myDir->Get("QCW");
0378 TMatrixD *QC;
0379 if (myQC == nullptr) {
0380 QC = new TMatrixD(0, 0);
0381 } else {
0382 QC = (TMatrixD *)myQC->Get("matrix");
0383 }
0384 return *QC;
0385 }
0386
0387 TMatrixD ReadPGInfo::giveSurvey(int wheel, int station, int sector) {
0388 TMatrixD *empty = new TMatrixD(0, 0);
0389 char *id = getId(wheel, station, sector);
0390 if (id == nullptr)
0391 return *empty;
0392 TDirectoryFile *myDir = (TDirectoryFile *)rootFile->Get(id);
0393 TDirectoryFile *mySur = (TDirectoryFile *)myDir->Get("Survey");
0394 TMatrixD *Survey = (TMatrixD *)mySur->Get("matrix_layer");
0395 return *Survey;
0396 }