|
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
// system include files
#include <memory>
// user include files
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
#include "CondFormats/SiStripObjects/interface/SiStripDetVOff.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.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 SiStripDetVOffFakeBuilder : public edm::one::EDAnalyzer<> {
public:
explicit SiStripDetVOffFakeBuilder(const edm::ParameterSet& iConfig);
~SiStripDetVOffFakeBuilder() override;
virtual void initialize(const edm::EventSetup&);
void analyze(const edm::Event&, const edm::EventSetup&) override;
private:
bool printdebug_;
const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkGeomToken_;
std::vector<uint32_t> detids;
};
using namespace std;
using namespace cms;
SiStripDetVOffFakeBuilder::SiStripDetVOffFakeBuilder(const edm::ParameterSet& iConfig)
: printdebug_(iConfig.getUntrackedParameter<bool>("printDebug", false)), tkGeomToken_(esConsumes()) {}
SiStripDetVOffFakeBuilder::~SiStripDetVOffFakeBuilder() = default;
void SiStripDetVOffFakeBuilder::initialize(const edm::EventSetup& iSetup) {
const auto& tkGeom = iSetup.getData(tkGeomToken_);
edm::LogInfo("SiStripDetVOffFakeBuilder") << " There are " << tkGeom.detUnits().size() << " detectors" << std::endl;
for (const auto& it : tkGeom.detUnits()) {
if (dynamic_cast<StripGeomDetUnit const*>(it) != nullptr) {
uint32_t detid = (it->geographicalId()).rawId();
const StripTopology& p = dynamic_cast<StripGeomDetUnit const*>(it)->specificTopology();
unsigned short Nstrips = p.nstrips();
if (Nstrips < 1 || Nstrips > 768) {
edm::LogError("SiStripDetVOffFakeBuilder")
<< " Problem with Number of strips in detector.. " << p.nstrips() << " Exiting program" << endl;
exit(1);
}
detids.push_back(detid);
if (printdebug_)
edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid " << detid;
}
}
}
void SiStripDetVOffFakeBuilder::analyze(const edm::Event& evt, const edm::EventSetup& iSetup) {
initialize(iSetup);
unsigned int run = evt.id().run();
edm::LogInfo("SiStripDetVOffFakeBuilder")
<< "... creating dummy SiStripDetVOff Data for Run " << run << "\n " << std::endl;
SiStripDetVOff SiStripDetVOff_;
// std::vector<uint32_t> TheDetIdHVVector;
for (std::vector<uint32_t>::const_iterator it = detids.begin(); it != detids.end(); it++) {
//Generate HV and LV for each channel, if at least one of the two is off fill the value
int hv = rand() % 20;
int lv = rand() % 20;
if (hv <= 2) {
edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " << *it << " HV\t OFF" << std::endl;
SiStripDetVOff_.put(*it, 1, -1);
// TheDetIdHVVector.push_back(*it);
}
if (lv <= 2) {
edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " << *it << " LV\t OFF" << std::endl;
SiStripDetVOff_.put(*it, -1, 1);
// TheDetIdHVVector.push_back(*it);
}
if (lv <= 2 || hv <= 2)
edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " << *it << " V\t OFF" << std::endl;
}
// SiStripDetVOff_.put(TheDetIdHVVector);
//End now write DetVOff data in DB
edm::Service<cond::service::PoolDBOutputService> mydbservice;
if (mydbservice.isAvailable()) {
try {
if (mydbservice->isNewTagRequest("SiStripDetVOffRcd")) {
mydbservice->createOneIOV<SiStripDetVOff>(SiStripDetVOff_, mydbservice->beginOfTime(), "SiStripDetVOffRcd");
} else {
mydbservice->appendOneIOV<SiStripDetVOff>(SiStripDetVOff_, mydbservice->currentTime(), "SiStripDetVOffRcd");
}
} catch (const cond::Exception& er) {
edm::LogError("SiStripDetVOffFakeBuilder") << er.what() << std::endl;
} catch (const std::exception& er) {
edm::LogError("SiStripDetVOffFakeBuilder") << "caught std::exception " << er.what() << std::endl;
}
} else {
edm::LogError("SiStripDetVOffFakeBuilder") << "Service is unavailable" << std::endl;
}
}
#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(SiStripDetVOffFakeBuilder);
|