|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
|
// -*- C++ -*-
//
// Package: SiStripDetInfoFileWriter
// Class: SiStripDetInfoFileWriter
//
/**\class SiStripDetInfoFileWriter SiStripDetInfoFileWriter.cc CalibTracker/SiStripCommon/src/SiStripDetInfoFileWriter.cc
Description: <one line class summary>
Implementation:
<Notes on implementation>
*/
//
// 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);
|