File indexing completed on 2021-02-14 13:08:05
0001 #include "Geometry/EcalTestBeam/interface/EcalTBHodoscopeGeometry.h"
0002 #include "SimDataFormats/EcalTestBeam/interface/HodoscopeDetId.h"
0003
0004 #include <vector>
0005 #include <iostream>
0006
0007 int main() {
0008 EcalTBHodoscopeGeometry theTestGeom;
0009
0010 for (int j = 0; j < theTestGeom.getNPlanes(); ++j) {
0011 for (int i = 0; i < 1000; ++i) {
0012 std::cout << "Position " << -17. + 34. / 1000. * i << " Plane " << j << std::endl;
0013 std::vector<int> firedFibres = theTestGeom.getFiredFibresInPlane(-17. + 34. / 1000. * i, j);
0014 for (int firedFibre : firedFibres) {
0015 std::cout << firedFibre << std::endl;
0016
0017 HodoscopeDetId myDetId = HodoscopeDetId(j, (int)firedFibre);
0018 std::cout << myDetId << std::endl;
0019 }
0020 }
0021 }
0022
0023 return 0;
0024 }