|
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
122
123
124
125
|
#ifndef _CSCBADWIRESCONDITIONS_H
#define _CSCBADWIRESCONDITIONS_H
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/SourceFactory.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <cmath>
#include <memory>
#include "CondFormats/CSCObjects/interface/CSCBadWires.h"
#include "CondFormats/DataRecord/interface/CSCBadWiresRcd.h"
#include <DataFormats/MuonDetId/interface/CSCDetId.h>
class CSCBadWiresConditions : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
CSCBadWiresConditions(const edm::ParameterSet &);
~CSCBadWiresConditions() override;
inline static CSCBadWires *prefillBadWires();
typedef std::unique_ptr<CSCBadWires> ReturnType;
ReturnType produceBadWires(const CSCBadWiresRcd &);
private:
// ----------member data ---------------------------
void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
const edm::IOVSyncValue &,
edm::ValidityInterval &) override;
CSCBadWires *cndbBadWires;
};
#include <fstream>
#include <iostream>
#include <vector>
// to workaround plugin library
inline CSCBadWires *CSCBadWiresConditions::prefillBadWires() {
CSCBadWires *cndbbadwires = new CSCBadWires();
int new_index, new_chan;
int new_layer, new_channel, new_flag1, new_flag2, new_flag3, new_pointer;
std::vector<int> new_index_id;
std::vector<short int> new_layer_id;
std::vector<short int> new_chan_id;
std::vector<int> new_badchannels;
std::vector<short int> new_flag1_id;
std::vector<short int> new_flag2_id;
std::vector<short int> new_flag3_id;
std::vector<int> new_cham_id;
std::vector<int> new_point;
// when bad channels exist use the following
int new_nrlines1 = 0;
int new_nrlines2 = 0;
// when empty
// int new_nrlines1=-1;
// int new_nrlines2=-1;
std::ifstream newdata1;
std::ifstream newdata2;
newdata1.open("badwires1.dat", std::ios::in);
if (!newdata1) {
std::cerr << "Error: badwires1.dat -> no such file!" << std::endl;
exit(1);
}
newdata2.open("badwires2.dat", std::ios::in); // channelcontainer
if (!newdata2) {
std::cerr << "Error: badwires2.dat -> no such file!" << std::endl;
exit(1);
}
while (!newdata1.eof()) {
newdata1 >> new_index >> new_pointer >> new_chan;
new_index_id.push_back(new_index);
new_point.push_back(new_pointer);
new_badchannels.push_back(new_chan);
new_nrlines1++;
}
newdata1.close();
while (!newdata2.eof()) {
newdata2 >> new_layer >> new_channel >> new_flag1 >> new_flag2 >> new_flag3;
new_layer_id.push_back(new_layer);
new_chan_id.push_back(new_channel);
new_flag1_id.push_back(new_flag1);
new_flag2_id.push_back(new_flag2);
new_flag3_id.push_back(new_flag3);
new_nrlines2++;
}
newdata2.close();
CSCBadWires::BadChamberContainer &itemvector1 = cndbbadwires->chambers;
itemvector1.resize(new_nrlines1);
CSCBadWires::BadChannelContainer &itemvector2 = cndbbadwires->channels;
itemvector2.resize(new_nrlines2);
cndbbadwires->numberOfBadChannels = new_nrlines2;
for (int i = 0; i < new_nrlines1; i++) {
itemvector1[i].chamber_index = new_index_id[i];
itemvector1[i].pointer = new_point[i];
itemvector1[i].bad_channels = new_badchannels[i];
}
for (int j = 0; j < new_nrlines2; ++j) {
itemvector2[j].layer = new_layer_id[j];
itemvector2[j].channel = new_chan_id[j];
itemvector2[j].flag1 = new_flag1_id[j];
itemvector2[j].flag2 = new_flag2_id[j];
itemvector2[j].flag3 = new_flag3_id[j];
}
return cndbbadwires;
}
#endif
|