Belle II Software release-09-00-04
PXDMappingLookup.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9#include <pxd/unpacking/PXDMappingLookup.h>
10#include <framework/logging/Logger.h>
11
12#include <stdio.h>
13#include <vector>
14
15using namespace std;
16using namespace Belle2;
17using namespace Belle2::PXD;
18
20int PXDMappingLookup::getDCDID(const int u, const int v, const VxdID sensorID)
21{
22
23 //const SensorInfo& info = dynamic_cast<const SensorInfo&>(VXD::GeoCache::get(
24 // sensorID));
25 int ret = 0;
26 int iLayer = sensorID.getLayerNumber();
27 int iSensor = sensorID.getSensorNumber();
28 int chUMap;
29 chUMap = u;
30 int Direction = 1;
31 if (((iLayer == 1) && (iSensor == 1)) || ((iLayer == 2) && (iSensor == 2))) {
32 // inverted counting in U direction:
33 chUMap = 249 - u;
34 Direction = -1;
35 }
36 float fID;
37 fID = (float)chUMap / 62.5;
38 if ((chUMap == 62) || (chUMap == 187)) {
39 if ((v % 4) >= 2) fID = fID + Direction;
40 }
41 ret = (int)fID;
42 return ret + 1;
43} //getDCDID
44
46int PXDMappingLookup::getSWBID(const int v)
47{
48 return (int)(v / 128) + 1;
49} //getSWBID
50
52void PXDMappingLookup::map_rc_to_uv_IF_OB(unsigned int& v_cellID, unsigned int& u_cellID, const unsigned int dhp_id,
53 const unsigned int dhe_ID)
54{
55 unsigned int DCD_channel = 0;
56 unsigned int Drain = 0;
57 unsigned int row = 0;
58 static std::vector <unsigned int> LUT_IF_OB = {
59 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 1023, 1023, 1023, 1023, 1023, 1023, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239,
60 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207,
61 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175,
62 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143,
63 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
64 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79,
65 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
66 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
67 490, 491, 492, 493, 494, 495, 496, 497, 498, 499, 1023, 1023, 1023, 1023, 1023, 1023, 474, 475, 476, 477, 478, 479, 480, 481, 482, 483, 484, 485, 486, 487, 488, 489,
68 458, 459, 460, 461, 462, 463, 464, 465, 466, 467, 468, 469, 470, 471, 472, 473, 442, 443, 444, 445, 446, 447, 448, 449, 450, 451, 452, 453, 454, 455, 456, 457,
69 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, 441, 410, 411, 412, 413, 414, 415, 416, 417, 418, 419, 420, 421, 422, 423, 424, 425,
70 394, 395, 396, 397, 398, 399, 400, 401, 402, 403, 404, 405, 406, 407, 408, 409, 378, 379, 380, 381, 382, 383, 384, 385, 386, 387, 388, 389, 390, 391, 392, 393,
71 362, 363, 364, 365, 366, 367, 368, 369, 370, 371, 372, 373, 374, 375, 376, 377, 346, 347, 348, 349, 350, 351, 352, 353, 354, 355, 356, 357, 358, 359, 360, 361,
72 330, 331, 332, 333, 334, 335, 336, 337, 338, 339, 340, 341, 342, 343, 344, 345, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323, 324, 325, 326, 327, 328, 329,
73 298, 299, 300, 301, 302, 303, 304, 305, 306, 307, 308, 309, 310, 311, 312, 313, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297,
74 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265,
75 740, 741, 742, 743, 744, 745, 746, 747, 748, 749, 1023, 1023, 1023, 1023, 1023, 1023, 724, 725, 726, 727, 728, 729, 730, 731, 732, 733, 734, 735, 736, 737, 738, 739,
76 708, 709, 710, 711, 712, 713, 714, 715, 716, 717, 718, 719, 720, 721, 722, 723, 692, 693, 694, 695, 696, 697, 698, 699, 700, 701, 702, 703, 704, 705, 706, 707,
77 676, 677, 678, 679, 680, 681, 682, 683, 684, 685, 686, 687, 688, 689, 690, 691, 660, 661, 662, 663, 664, 665, 666, 667, 668, 669, 670, 671, 672, 673, 674, 675,
78 644, 645, 646, 647, 648, 649, 650, 651, 652, 653, 654, 655, 656, 657, 658, 659, 628, 629, 630, 631, 632, 633, 634, 635, 636, 637, 638, 639, 640, 641, 642, 643,
79 612, 613, 614, 615, 616, 617, 618, 619, 620, 621, 622, 623, 624, 625, 626, 627, 596, 597, 598, 599, 600, 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611,
80 580, 581, 582, 583, 584, 585, 586, 587, 588, 589, 590, 591, 592, 593, 594, 595, 564, 565, 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 576, 577, 578, 579,
81 548, 549, 550, 551, 552, 553, 554, 555, 556, 557, 558, 559, 560, 561, 562, 563, 532, 533, 534, 535, 536, 537, 538, 539, 540, 541, 542, 543, 544, 545, 546, 547,
82 516, 517, 518, 519, 520, 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 500, 501, 502, 503, 504, 505, 506, 507, 508, 509, 510, 511, 512, 513, 514, 515,
83 990, 991, 992, 993, 994, 995, 996, 997, 998, 999, 1023, 1023, 1023, 1023, 1023, 1023, 974, 975, 976, 977, 978, 979, 980, 981, 982, 983, 984, 985, 986, 987, 988, 989,
84 958, 959, 960, 961, 962, 963, 964, 965, 966, 967, 968, 969, 970, 971, 972, 973, 942, 943, 944, 945, 946, 947, 948, 949, 950, 951, 952, 953, 954, 955, 956, 957,
85 926, 927, 928, 929, 930, 931, 932, 933, 934, 935, 936, 937, 938, 939, 940, 941, 910, 911, 912, 913, 914, 915, 916, 917, 918, 919, 920, 921, 922, 923, 924, 925,
86 894, 895, 896, 897, 898, 899, 900, 901, 902, 903, 904, 905, 906, 907, 908, 909, 878, 879, 880, 881, 882, 883, 884, 885, 886, 887, 888, 889, 890, 891, 892, 893,
87 862, 863, 864, 865, 866, 867, 868, 869, 870, 871, 872, 873, 874, 875, 876, 877, 846, 847, 848, 849, 850, 851, 852, 853, 854, 855, 856, 857, 858, 859, 860, 861,
88 830, 831, 832, 833, 834, 835, 836, 837, 838, 839, 840, 841, 842, 843, 844, 845, 814, 815, 816, 817, 818, 819, 820, 821, 822, 823, 824, 825, 826, 827, 828, 829,
89 798, 799, 800, 801, 802, 803, 804, 805, 806, 807, 808, 809, 810, 811, 812, 813, 782, 783, 784, 785, 786, 787, 788, 789, 790, 791, 792, 793, 794, 795, 796, 797,
90 766, 767, 768, 769, 770, 771, 772, 773, 774, 775, 776, 777, 778, 779, 780, 781, 750, 751, 752, 753, 754, 755, 756, 757, 758, 759, 760, 761, 762, 763, 764, 765
91 };
92
93 // B2INFO("Remapped :: From COL $" << u_cellID << " ROW $" << v_cellID);
94 DCD_channel = (u_cellID << 2) + (v_cellID & 0x3) + 256 * dhp_id;
95 Drain = LUT_IF_OB[DCD_channel];
96 u_cellID = Drain >> 2;
97 row = (v_cellID & ~0x3) + (Drain & 0x3); // no ~ bei drain
98 row = row % 768;
99 if ((dhe_ID & 0x20) == 0) {//if inner module
100 v_cellID = 768 - 1 - row ;
101
102 } else {//if outer module
103 v_cellID = row ;
104 }
105// B2DEBUG(99,"Remapped ::To COL COL $" << u_cellID << " ROW $" << v_cellID);
106}
107
108void PXDMappingLookup::map_uv_to_rc_IF_OB(unsigned int& v_cellID, unsigned int& u_cellID,
109 __attribute__((unused)) unsigned int& dhp_id,
110 const unsigned int dhe_ID)
111{
112 B2FATAL("Code to be written");
113 unsigned int row;
114 if ((dhe_ID & 0x20) == 0) { //if inner module
115 row = 768 - 1 - v_cellID ;
116 } else { //if outer module
117 row = v_cellID;
118 }
119 v_cellID = row;
120 u_cellID = 0;
121}
122
124void PXDMappingLookup::map_rc_to_uv_IB_OF(unsigned int& v_cellID, unsigned int& u_cellID, const unsigned int dhp_id,
125 const unsigned int dhe_ID)
126{
127 unsigned int DCD_channel = 0;
128 unsigned int Drain = 0;
129 unsigned int row = 0;
130 static std::vector <unsigned int> LUT_IB_OF = {
131 243, 242, 241, 240, 247, 246, 245, 244, 251, 250, 1023, 1023, 1023, 1023, 1023, 1023, 227, 226, 225, 224, 231, 230, 229, 228, 235, 234, 233, 232, 239, 238, 237, 236,
132 211, 210, 209, 208, 215, 214, 213, 212, 219, 218, 217, 216, 223, 222, 221, 220, 195, 194, 193, 192, 199, 198, 197, 196, 203, 202, 201, 200, 207, 206, 205, 204,
133 179, 178, 177, 176, 183, 182, 181, 180, 187, 186, 185, 184, 191, 190, 189, 188, 163, 162, 161, 160, 167, 166, 165, 164, 171, 170, 169, 168, 175, 174, 173, 172,
134 147, 146, 145, 144, 151, 150, 149, 148, 155, 154, 153, 152, 159, 158, 157, 156, 131, 130, 129, 128, 135, 134, 133, 132, 139, 138, 137, 136, 143, 142, 141, 140,
135 115, 114, 113, 112, 119, 118, 117, 116, 123, 122, 121, 120, 127, 126, 125, 124, 99, 98, 97, 96, 103, 102, 101, 100, 107, 106, 105, 104, 111, 110, 109, 108,
136 83, 82, 81, 80, 87, 86, 85, 84, 91, 90, 89, 88, 95, 94, 93, 92, 67, 66, 65, 64, 71, 70, 69, 68, 75, 74, 73, 72, 79, 78, 77, 76,
137 51, 50, 49, 48, 55, 54, 53, 52, 59, 58, 57, 56, 63, 62, 61, 60, 35, 34, 33, 32, 39, 38, 37, 36, 43, 42, 41, 40, 47, 46, 45, 44,
138 19, 18, 17, 16, 23, 22, 21, 20, 27, 26, 25, 24, 31, 30, 29, 28, 3, 2, 1, 0, 7, 6, 5, 4, 11, 10, 9, 8, 15, 14, 13, 12,
139 489, 488, 495, 494, 493, 492, 499, 498, 497, 496, 1023, 1023, 1023, 1023, 1023, 1023, 473, 472, 479, 478, 477, 476, 483, 482, 481, 480, 487, 486, 485, 484, 491, 490,
140 457, 456, 463, 462, 461, 460, 467, 466, 465, 464, 471, 470, 469, 468, 475, 474, 441, 440, 447, 446, 445, 444, 451, 450, 449, 448, 455, 454, 453, 452, 459, 458,
141 425, 424, 431, 430, 429, 428, 435, 434, 433, 432, 439, 438, 437, 436, 443, 442, 409, 408, 415, 414, 413, 412, 419, 418, 417, 416, 423, 422, 421, 420, 427, 426,
142 393, 392, 399, 398, 397, 396, 403, 402, 401, 400, 407, 406, 405, 404, 411, 410, 377, 376, 383, 382, 381, 380, 387, 386, 385, 384, 391, 390, 389, 388, 395, 394,
143 361, 360, 367, 366, 365, 364, 371, 370, 369, 368, 375, 374, 373, 372, 379, 378, 345, 344, 351, 350, 349, 348, 355, 354, 353, 352, 359, 358, 357, 356, 363, 362,
144 329, 328, 335, 334, 333, 332, 339, 338, 337, 336, 343, 342, 341, 340, 347, 346, 313, 312, 319, 318, 317, 316, 323, 322, 321, 320, 327, 326, 325, 324, 331, 330,
145 297, 296, 303, 302, 301, 300, 307, 306, 305, 304, 311, 310, 309, 308, 315, 314, 281, 280, 287, 286, 285, 284, 291, 290, 289, 288, 295, 294, 293, 292, 299, 298,
146 265, 264, 271, 270, 269, 268, 275, 274, 273, 272, 279, 278, 277, 276, 283, 282, 249, 248, 255, 254, 253, 252, 259, 258, 257, 256, 263, 262, 261, 260, 267, 266,
147 743, 742, 741, 740, 747, 746, 745, 744, 751, 750, 1023, 1023, 1023, 1023, 1023, 1023, 727, 726, 725, 724, 731, 730, 729, 728, 735, 734, 733, 732, 739, 738, 737, 736,
148 711, 710, 709, 708, 715, 714, 713, 712, 719, 718, 717, 716, 723, 722, 721, 720, 695, 694, 693, 692, 699, 698, 697, 696, 703, 702, 701, 700, 707, 706, 705, 704,
149 679, 678, 677, 676, 683, 682, 681, 680, 687, 686, 685, 684, 691, 690, 689, 688, 663, 662, 661, 660, 667, 666, 665, 664, 671, 670, 669, 668, 675, 674, 673, 672,
150 647, 646, 645, 644, 651, 650, 649, 648, 655, 654, 653, 652, 659, 658, 657, 656, 631, 630, 629, 628, 635, 634, 633, 632, 639, 638, 637, 636, 643, 642, 641, 640,
151 615, 614, 613, 612, 619, 618, 617, 616, 623, 622, 621, 620, 627, 626, 625, 624, 599, 598, 597, 596, 603, 602, 601, 600, 607, 606, 605, 604, 611, 610, 609, 608,
152 583, 582, 581, 580, 587, 586, 585, 584, 591, 590, 589, 588, 595, 594, 593, 592, 567, 566, 565, 564, 571, 570, 569, 568, 575, 574, 573, 572, 579, 578, 577, 576,
153 551, 550, 549, 548, 555, 554, 553, 552, 559, 558, 557, 556, 563, 562, 561, 560, 535, 534, 533, 532, 539, 538, 537, 536, 543, 542, 541, 540, 547, 546, 545, 544,
154 519, 518, 517, 516, 523, 522, 521, 520, 527, 526, 525, 524, 531, 530, 529, 528, 503, 502, 501, 500, 507, 506, 505, 504, 511, 510, 509, 508, 515, 514, 513, 512,
155 989, 988, 995, 994, 993, 992, 999, 998, 997, 996, 1023, 1023, 1023, 1023, 1023, 1023, 973, 972, 979, 978, 977, 976, 983, 982, 981, 980, 987, 986, 985, 984, 991, 990,
156 957, 956, 963, 962, 961, 960, 967, 966, 965, 964, 971, 970, 969, 968, 975, 974, 941, 940, 947, 946, 945, 944, 951, 950, 949, 948, 955, 954, 953, 952, 959, 958,
157 925, 924, 931, 930, 929, 928, 935, 934, 933, 932, 939, 938, 937, 936, 943, 942, 909, 908, 915, 914, 913, 912, 919, 918, 917, 916, 923, 922, 921, 920, 927, 926,
158 893, 892, 899, 898, 897, 896, 903, 902, 901, 900, 907, 906, 905, 904, 911, 910, 877, 876, 883, 882, 881, 880, 887, 886, 885, 884, 891, 890, 889, 888, 895, 894,
159 861, 860, 867, 866, 865, 864, 871, 870, 869, 868, 875, 874, 873, 872, 879, 878, 845, 844, 851, 850, 849, 848, 855, 854, 853, 852, 859, 858, 857, 856, 863, 862,
160 829, 828, 835, 834, 833, 832, 839, 838, 837, 836, 843, 842, 841, 840, 847, 846, 813, 812, 819, 818, 817, 816, 823, 822, 821, 820, 827, 826, 825, 824, 831, 830,
161 797, 796, 803, 802, 801, 800, 807, 806, 805, 804, 811, 810, 809, 808, 815, 814, 781, 780, 787, 786, 785, 784, 791, 790, 789, 788, 795, 794, 793, 792, 799, 798,
162 765, 764, 771, 770, 769, 768, 775, 774, 773, 772, 779, 778, 777, 776, 783, 782, 749, 748, 755, 754, 753, 752, 759, 758, 757, 756, 763, 762, 761, 760, 767, 766
163 };
164
165 DCD_channel = (u_cellID << 2) + (v_cellID & 0x3) + 256 * dhp_id;
166 Drain = LUT_IB_OF[DCD_channel];
167 if (Drain >> 2 >= 250) u_cellID = 255;
168 else u_cellID = 250 - 1 - (Drain >> 2);
169// row = (v_cellID / 4) * 4 + Drain % 4;
170 row = (v_cellID & ~0x3) + ((~Drain) & 0x3); // ~ bei drain
171 row = row % 768;
172 if ((dhe_ID & 0x20) == 0) { //if inner module
173 v_cellID = 768 - 1 - row ;
174 } else { //if outer module
175 v_cellID = row ;
176 }
177// B2DEBUG(99,"Remapped ::To COL COL $" << u_cellID << " ROW $" << v_cellID);
178}
179
180void PXDMappingLookup::map_uv_to_rc_IB_OF(unsigned int& v_cellID, unsigned int& u_cellID,
181 __attribute__((unused)) unsigned int& dhp_id,
182 const unsigned int dhe_ID)
183{
184 B2FATAL("Code to be written");
185 // slow way until we have tables
186 unsigned int row;
187 if ((dhe_ID & 0x20) == 0) { //if inner module
188 row = 768 - 1 - v_cellID ;
189 } else { //if outer module
190 row = v_cellID;
191 }
192 v_cellID = row;
193 u_cellID = 0;
194}
195
196void PXDMappingLookup::write_mapping_to_file(void)
197{
198 FILE* file = fopen("lut.csv", "wt+");
199 if (file) {
200 for (unsigned int row = 0; row < 4; ++row) {
201 for (unsigned int dhp_id = 0; dhp_id < 4; ++dhp_id) {
202 for (unsigned int col = 0; col < 64; ++col) {
203 unsigned int coli;
204 unsigned int rowi;
205 coli = col;
206 rowi = row;
207 map_rc_to_uv_IF_OB(rowi, coli, dhp_id, 0x00);
208 fprintf(file, "%u; %u; %u; %u; %u; ", row, dhp_id, col, rowi, coli);
209 coli = col;
210 rowi = row;
211 map_rc_to_uv_IB_OF(rowi, coli, dhp_id, 0x00);
212 fprintf(file, "%u; %u\n", rowi, coli);
213 }
214 }
215 }
216
217 for (unsigned int row = 0; row < 4; ++row) {
218 for (unsigned int dhp_id = 0; dhp_id < 4; ++dhp_id) {
219 for (unsigned int col = 0; col < 64; ++col) {
220 unsigned int coli;
221 unsigned int rowi;
222 coli = col;
223 rowi = row;
224 map_rc_to_uv_IF_OB(rowi, coli, dhp_id, 0x01);
225 fprintf(file, "%u; %u; %u; %u; %u; ", row, dhp_id, col, rowi, coli);
226 coli = col;
227 rowi = row;
228 map_rc_to_uv_IB_OF(rowi, coli, dhp_id, 0x01);
229 fprintf(file, "%u; %u\n", rowi, coli);
230 }
231 }
232 }
233
234
235 for (unsigned int row = 0; row < 4; ++row) {
236 for (unsigned int dhp_id = 0; dhp_id < 4; ++dhp_id) {
237 for (unsigned int col = 0; col < 64; ++col) {
238 unsigned int coli;
239 unsigned int rowi;
240 coli = col;
241 rowi = row;
242 map_rc_to_uv_IF_OB(rowi, coli, dhp_id, 0x20);
243 fprintf(file, "%u; %u; %u; %u; %u; ", row, dhp_id, col, rowi, coli);
244 coli = col;
245 rowi = row;
246 map_rc_to_uv_IB_OF(rowi, coli, dhp_id, 0x20);
247 fprintf(file, "%u; %u\n", rowi, coli);
248 }
249 }
250 }
251
252
253 for (unsigned int row = 0; row < 4; ++row) {
254 for (unsigned int dhp_id = 0; dhp_id < 4; ++dhp_id) {
255 for (unsigned int col = 0; col < 64; ++col) {
256 unsigned int coli;
257 unsigned int rowi;
258 coli = col;
259 rowi = row;
260 map_rc_to_uv_IF_OB(rowi, coli, dhp_id, 0x21);
261 fprintf(file, "%u; %u; %u; %u; %u; ", row, dhp_id, col, rowi, coli);
262 coli = col;
263 rowi = row;
264 map_rc_to_uv_IB_OF(rowi, coli, dhp_id, 0x21);
265 fprintf(file, "%u; %u\n", rowi, coli);
266 }
267 }
268 }
269
270 fclose(file);
271 }
272}
273
274void PXDMappingLookup::write_inversemapping_to_file(void)
275{
276 B2FATAL("Code to be written");
277}
278
279
280void PXDMappingLookup::check(void)
281{
282 const unsigned int dhe_ID[4] = {0x02, 0x03, 0x22, 0x23};// one of each kind, IB, OB, IF, OF
283 for (int i = 0; i < 4; i++) {
284 unsigned int dhe = dhe_ID[i];
285 for (unsigned int u_org = 0; i < 768; i++) {
286 for (unsigned int v_org = 0; v_org < 250; v_org++) {
287 unsigned int r, c, u, v;
288 unsigned int dhp_id = 0;
289
290 // if we call IF-OB with wrong DHE ID, we do not expect the correct result _BUT_ still it should be consistent
291 c = u_org;
292 map_uv_to_rc_IF_OB(r, c, dhp_id, dhe);
293 v = r;
294 u = c;
295 map_rc_to_uv_IF_OB(v, u, dhp_id, dhe);
296
297 if (u != u_org || v != v_org) {
298 B2ERROR("Mapping failed (IF-OB)! DHE $%02X V/U %u %u -> R/C %u %u DHP (%u) -> V/U %u %u" <<
299 dhe << v_org << u_org << r << c << dhp_id << v << u);
300 }
301
302 // if we call IB-OF with wrong DHE ID, we do not expect the correct result _BUT_ still it should be consistent
303 c = u_org;
304 r = v_org;
305 map_uv_to_rc_IB_OF(r, c, dhp_id, dhe);
306 v = r;
307 u = c;
308 map_rc_to_uv_IB_OF(v, u, dhp_id, dhe);
309
310 if (u != u_org || v != v_org) {
311 B2ERROR("Mapping failed (IB-OF)! DHE $%02X V/U %u %u -> R/C %u %u DHP (%u) -> V/U %u %u" <<
312 dhe << v_org << u_org << r << c << dhp_id << v << u);
313 }
314 }
315 }
316 }
317}
Belle2::PXD::PXDMappingLookup::getDCDID
static int getDCDID(const int u, const int v, const VxdID sensorID)
get ID of DCD for giving pixel, range: 1..4.
Definition: PXDMappingLookup.cc:20
Belle2::PXD::PXDMappingLookup::check
static void check(void)
consistency check
Definition: PXDMappingLookup.cc:280
Belle2::PXD::PXDMappingLookup::map_rc_to_uv_IF_OB
static void map_rc_to_uv_IF_OB(unsigned int &row_u, unsigned int &col_v, const unsigned int dhp_id, const unsigned int dhe_ID)
Maps row/col of inner forward (IF) and outer backward (OB) modules of the PXD to U/V cell.
Definition: PXDMappingLookup.cc:52
Belle2::PXD::PXDMappingLookup::getSWBID
static int getSWBID(const int v)
get ID of SWB for giving pixel, range: 1..6.
Definition: PXDMappingLookup.cc:46
Belle2::PXD::PXDMappingLookup::map_rc_to_uv_IB_OF
static void map_rc_to_uv_IB_OF(unsigned int &row_u, unsigned int &col_v, const unsigned int dhp_id, const unsigned int dhe_ID)
Maps row/cols of inner backward (IB) and outer forward (OF) modules of the PXD to U/V cell.
Definition: PXDMappingLookup.cc:124
Belle2::PXD::PXDMappingLookup::map_uv_to_rc_IF_OB
static void map_uv_to_rc_IF_OB(unsigned int &row_u, unsigned int &col_v, unsigned int &dhp_id, const unsigned int dhe_ID)
Maps U/V cell of inner forward (IF) and outer backward (OB) modules of the PXD to row/col/dhpid.
Definition: PXDMappingLookup.cc:108
Belle2::PXD::PXDMappingLookup::map_uv_to_rc_IB_OF
static void map_uv_to_rc_IB_OF(unsigned int &row_u, unsigned int &col_v, unsigned int &dhp_id, const unsigned int dhe_ID)
Maps U/V cell of inner backward (IB) and outer forward (OF) modules of the PXD to row/col/dhpid.
Definition: PXDMappingLookup.cc:180
Belle2::PXD::PXDMappingLookup::write_mapping_to_file
static void write_mapping_to_file(void)
write out LUT to file for cross-check
Definition: PXDMappingLookup.cc:196
Belle2::PXD::PXDMappingLookup::write_inversemapping_to_file
static void write_inversemapping_to_file(void)
write out inverse LUT to file for cross-check
Definition: PXDMappingLookup.cc:274
Belle2::VxdID
Class to uniquely identify a any structure of the PXD and SVD.
Definition: VxdID.h:33
Belle2::VxdID::getSensorNumber
baseType getSensorNumber() const
Get the sensor id.
Definition: VxdID.h:100
Belle2::VxdID::getLayerNumber
baseType getLayerNumber() const
Get the layer id.
Definition: VxdID.h:96
Belle2::PXD
Namespace to encapsulate code needed for simulation and reconstrucion of the PXD.
Definition: PXDCalibrationUtilities.h:18
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17
std
STL namespace.