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