L1TCaloLayer1/test/testUCTGeometry.cpp

0001 #include <iostream>
0002 #include <stdlib.h>
0003 #include <stdint.h>
0004
0005 using namespace std;
0006
0007 #include "L1Trigger/L1TCaloLayer1/src/UCTGeometry.hh"
0008
0009 int main(int argc, char** argv) {
0010   UCTGeometry g;
0011   for (int caloPhi = 1; caloPhi <= l1tcalo::MaxCaloPhi; caloPhi++) {
0012     for (int caloEta = -l1tcalo::MaxCaloEta; caloEta <= l1tcalo::MaxCaloEta; caloEta++) {
0013       if (caloEta == 0 || abs(caloEta) == 29)
0014         continue;
0015       uint32_t crt = g.getCrate(caloEta, caloPhi);
0016       uint32_t crd = g.getCard(caloEta, caloPhi);
0017       uint32_t rgn = g.getRegion(caloEta, caloPhi);
0018       uint32_t eta = g.getiEta(caloEta);
0019       uint32_t phi = g.getiPhi(caloPhi);
0020       int cEta = g.getCaloEtaIndex((caloEta < 0), rgn, eta);
0021       int cPhi = g.getCaloPhiIndex(crt, crd, rgn, phi);
0022       if (cEta != caloEta || cPhi != caloPhi) {
0023         cout << "(caloEta, caloPhi) = (" << caloEta << ", " << caloPhi << ") "
0024              << "(crt,crd,rgn,eta,phi) = (" << crt << ", " << crd << ", " << rgn << ", " << eta << ", " << phi << ")"
0025              << endl;
0026         cerr << "(caloEta, caloPhi) = (" << cEta << ", " << cPhi << ") "
0027              << "Obtained instead   " << endl;
0028       }
0029     }
0030   }
0031   return 0;
0032 }