DTGeometryTest.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ESTransientHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"

#include <iostream>
#include <string>

using namespace std;
using namespace cms;
using namespace edm;

class DTGeometryTest : public one::EDAnalyzer<> {
public:
  explicit DTGeometryTest(const ParameterSet&);

  void beginJob() override {}
  void analyze(Event const& iEvent, EventSetup const&) override;
  void endJob() override {}

private:  
  const string m_label;
  const edm::ESGetToken<DTGeometry, MuonGeometryRecord> m_token;
};

DTGeometryTest::DTGeometryTest(const ParameterSet& iConfig)
  : m_label(iConfig.getUntrackedParameter<string>("fromDataLabel", "")),
    m_token(esConsumes<DTGeometry, MuonGeometryRecord>(edm::ESInputTag{"", m_label}))
{}

void
DTGeometryTest::analyze(const Event&, const EventSetup& iEventSetup)
{
  LogVerbatim("Geometry") << "DTGeometryTest::analyze: " << m_label;
  ESTransientHandle<DTGeometry> pDD = iEventSetup.getTransientHandle(m_token);
  
  LogVerbatim("Geometry") << " Geometry node for DTGeom is  " << &(*pDD);   
  LogVerbatim("Geometry") << " I have " << pDD->detTypes().size()    << " detTypes";
  LogVerbatim("Geometry") << " I have " << pDD->detUnits().size()    << " detUnits";
  LogVerbatim("Geometry") << " I have " << pDD->dets().size()        << " dets";
  LogVerbatim("Geometry") << " I have " << pDD->layers().size()      << " layers";
  LogVerbatim("Geometry") << " I have " << pDD->superLayers().size() << " superlayers";
  LogVerbatim("Geometry") << " I have " << pDD->chambers().size()    << " chambers";
  
  // check chamber
  LogVerbatim("Geometry") << "CHAMBERS " << string(120, '-');

  LogVerbatim("Geometry").log([&](auto& log) {
      for(auto det : pDD->chambers()){
    const BoundPlane& surf = det->surface();
    log << "Chamber " << det->id() 
        << " Position " << surf.position()
        << " normVect " << surf.normalVector() 
        << " bounds W/H/L: " << surf.bounds().width() << "/" 
        << surf.bounds().thickness() << "/" << surf.bounds().length() << "\n";
      }
    });
  LogVerbatim("Geometry") << "END " << string(120, '-');
  
  // check superlayers
  LogVerbatim("Geometry") << "SUPERLAYERS " << string(120, '-');
  LogVerbatim("Geometry").log([&](auto& log) {
      for(auto det : pDD->superLayers()) {
    const BoundPlane& surf = det->surface();
    log << "SuperLayer " << det->id() 
        << " chamber " << det->chamber()->id()
        << " Position " << surf.position()
        << " normVect " << surf.normalVector() 
        << " bounds W/H/L: " << surf.bounds().width() << "/" 
        << surf.bounds().thickness() << "/" << surf.bounds().length() << "\n";
      }
    });
  LogVerbatim("Geometry") << "END " << string(120, '-');

  // check layers
  LogVerbatim("Geometry") << "LAYERS " << string(120, '-');

  LogVerbatim("Geometry").log([&](auto& log) {
      for(auto det : pDD->layers()){
    const DTTopology& topo = det->specificTopology();
    const BoundPlane& surf=det->surface();
    log << "Layer " << det->id() << " SL " << det->superLayer()->id() 
        << " chamber " << det->chamber()->id() 
        << " Topology W/H/L: " 
        << topo.cellWidth() << "/" << topo.cellHeight() << "/" << topo.cellLenght() 
        << " first/last/# wire " << topo.firstChannel() << "/" << topo.lastChannel() << "/" << topo.channels()
        << " Position " << surf.position()
        << " normVect " << surf.normalVector() 
        << " bounds W/H/L: " << surf.bounds().width() << "/" 
        << surf.bounds().thickness() << "/" << surf.bounds().length() << "\n";
      }
    });
  LogVerbatim("Geometry") << "END " << string(120, '-');
}

DEFINE_FWK_MODULE(DTGeometryTest);