SiStripDetInfoFileWriter.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripDetInfoFileWriter
// Class:      SiStripDetInfoFileWriter
//
//
// Original Author:  G. Bruno
//         Created:  Mon Nov 20 10:04:31 CET 2006
//
//

// C++ includes
#include <string>
#include <iostream>
#include <fstream>

// User includes
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

class SiStripDetInfoFileWriter : public edm::one::EDAnalyzer<edm::one::WatchRuns> {
public:
  explicit SiStripDetInfoFileWriter(const edm::ParameterSet&);
  ~SiStripDetInfoFileWriter() override;

private:
  void beginRun(const edm::Run&, const edm::EventSetup& iSetup) override;
  void analyze(const edm::Event&, const edm::EventSetup&) override{};
  void endRun(const edm::Run&, const edm::EventSetup& iSetup) override{};

private:
  std::ofstream outputFile_;
  std::string filePath_;
  edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
};

using namespace cms;
using namespace std;

SiStripDetInfoFileWriter::SiStripDetInfoFileWriter(const edm::ParameterSet& iConfig) {
  edm::LogInfo("SiStripDetInfoFileWriter::SiStripDetInfoFileWriter");

  filePath_ = iConfig.getUntrackedParameter<std::string>("FilePath", std::string("SiStripDetInfo.dat"));
  tkGeomToken_ = esConsumes<edm::Transition::BeginRun>();
}

SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter() {
  edm::LogInfo("SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter");
}

void SiStripDetInfoFileWriter::beginRun(const edm::Run&, const edm::EventSetup& iSetup) {
  outputFile_.open(filePath_.c_str());

  if (outputFile_.is_open()) {
    const auto& dd = iSetup.getData(tkGeomToken_);

    edm::LogInfo("SiStripDetInfoFileWriter::beginRun - got geometry  ") << std::endl;

    edm::LogInfo("SiStripDetInfoFileWriter") << " There are " << dd.detUnits().size() << " detectors" << std::endl;

    for (const auto& it : dd.detUnits()) {
      const StripGeomDetUnit* mit = dynamic_cast<StripGeomDetUnit const*>(it);

      if (mit != nullptr) {
        uint32_t detid = (mit->geographicalId()).rawId();
        double stripLength = mit->specificTopology().stripLength();
        unsigned short numberOfAPVs = mit->specificTopology().nstrips() / 128;
        float thickness = mit->specificSurface().bounds().thickness();

        if (numberOfAPVs < 1 || numberOfAPVs > 6) {
          edm::LogError("SiStripDetInfoFileWriter")
              << " Problem with Number of strips in detector.. " << mit->specificTopology().nstrips()
              << "Will not write this entry to file" << endl;
          continue;
        }

        outputFile_ << detid << " " << numberOfAPVs << " " << stripLength << " " << thickness << "\n";
      }
    }
    outputFile_.close();
  } else {
    edm::LogError("SiStripDetInfoFileWriter::beginRun - Unable to open file") << endl;
    return;
  }
}

DEFINE_FWK_MODULE(SiStripDetInfoFileWriter);