SiStripCommon/plugins/SiStripDetInfoFileWriter.cc

0001 // -*- C++ -*-
0002 //
0003 // Package:    SiStripDetInfoFileWriter
0004 // Class:      SiStripDetInfoFileWriter
0005 //
0006 /**\class SiStripDetInfoFileWriter SiStripDetInfoFileWriter.cc CalibTracker/SiStripCommon/src/SiStripDetInfoFileWriter.cc
0007
0008  Description: <one line class summary>
0009
0010  Implementation:
0011      <Notes on implementation>
0012 */
0013 //
0014 // Original Author:  G. Bruno
0015 //         Created:  Mon Nov 20 10:04:31 CET 2006
0016 //
0017 //
0018
0019 // C++ includes
0020 #include <string>
0021 #include <iostream>
0022 #include <fstream>
0023
0024 // User includes
0025 #include "FWCore/Framework/interface/MakerMacros.h"
0026 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
0027 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0028 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0029 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
0030 #include "Geometry/CommonTopologies/interface/StripTopology.h"
0031 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
0032 #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
0033 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
0034
0035 class SiStripDetInfoFileWriter : public edm::one::EDAnalyzer<edm::one::WatchRuns> {
0036 public:
0037   explicit SiStripDetInfoFileWriter(const edm::ParameterSet&);
0038   ~SiStripDetInfoFileWriter() override;
0039
0040 private:
0041   void beginRun(const edm::Run&, const edm::EventSetup& iSetup) override;
0042   void analyze(const edm::Event&, const edm::EventSetup&) override{};
0043   void endRun(const edm::Run&, const edm::EventSetup& iSetup) override{};
0044
0045 private:
0046   std::ofstream outputFile_;
0047   std::string filePath_;
0048   edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
0049 };
0050
0051 using namespace cms;
0052 using namespace std;
0053
0054 SiStripDetInfoFileWriter::SiStripDetInfoFileWriter(const edm::ParameterSet& iConfig) {
0055   edm::LogInfo("SiStripDetInfoFileWriter::SiStripDetInfoFileWriter");
0056
0057   filePath_ = iConfig.getUntrackedParameter<std::string>("FilePath", std::string("SiStripDetInfo.dat"));
0058   tkGeomToken_ = esConsumes<edm::Transition::BeginRun>();
0059 }
0060
0061 SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter() {
0062   edm::LogInfo("SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter");
0063 }
0064
0065 void SiStripDetInfoFileWriter::beginRun(const edm::Run&, const edm::EventSetup& iSetup) {
0066   outputFile_.open(filePath_.c_str());
0067
0068   if (outputFile_.is_open()) {
0069     const auto& dd = iSetup.getData(tkGeomToken_);
0070
0071     edm::LogInfo("SiStripDetInfoFileWriter::beginRun - got geometry  ") << std::endl;
0072
0073     edm::LogInfo("SiStripDetInfoFileWriter") << " There are " << dd.detUnits().size() << " detectors" << std::endl;
0074
0075     for (const auto& it : dd.detUnits()) {
0076       const StripGeomDetUnit* mit = dynamic_cast<StripGeomDetUnit const*>(it);
0077
0078       if (mit != nullptr) {
0079         uint32_t detid = (mit->geographicalId()).rawId();
0080         double stripLength = mit->specificTopology().stripLength();
0081         unsigned short numberOfAPVs = mit->specificTopology().nstrips() / 128;
0082         float thickness = mit->specificSurface().bounds().thickness();
0083
0084         if (numberOfAPVs < 1 || numberOfAPVs > 6) {
0085           edm::LogError("SiStripDetInfoFileWriter")
0086               << " Problem with Number of strips in detector.. " << mit->specificTopology().nstrips()
0087               << "Will not write this entry to file" << endl;
0088           continue;
0089         }
0090
0091         outputFile_ << detid << " " << numberOfAPVs << " " << stripLength << " " << thickness << "\n";
0092       }
0093     }
0094     outputFile_.close();
0095   } else {
0096     edm::LogError("SiStripDetInfoFileWriter::beginRun - Unable to open file") << endl;
0097     return;
0098   }
0099 }
0100
0101 DEFINE_FWK_MODULE(SiStripDetInfoFileWriter);