|
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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
|
#include "CalibTracker/SiStripDCS/test/Synchronization/SyncDCSO2O.h"
#include "CoralBase/TimeStamp.h"
#include "CondFormats/Common/interface/Time.h"
#include "CondFormats/Common/interface/TimeConversions.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include <iostream>
#include <algorithm>
#include "TFile.h"
#include "TCanvas.h"
SyncDCSO2O::SyncDCSO2O(const edm::ParameterSet& iConfig) : dcsToken_(esConsumes()) {
// get all digi collections
digiProducersList_ = iConfig.getParameter<Parameters>("DigiProducersList");
}
SyncDCSO2O::~SyncDCSO2O() {
std::cout << "Analyzed events with digis = " << timeInfo_.size() << std::endl;
// First sort the timeInfo vector by time
std::sort(timeInfo_.begin(), timeInfo_.end(), SortByTime());
TFile* outputFile = new TFile("digisAndHVvsTime.root", "RECREATE");
outputFile->cd();
TH1F* digis = new TH1F("digis", "digis", timeInfo_.size(), 0, timeInfo_.size());
TH1F* digisWithMasking = new TH1F("digisWithMasking", "digisWithMasking", timeInfo_.size(), 0, timeInfo_.size());
TH1F* HVoff = new TH1F("HVoff", "HVoff", timeInfo_.size(), 0, timeInfo_.size());
TH1F* time = new TH1F("time", "time", timeInfo_.size(), 0, timeInfo_.size());
std::vector<TimeInfo>::const_iterator it = timeInfo_.begin();
// Float_t * timeArray = new Float_t[timeInfo_.size()];
unsigned int i = 1;
for (; it != timeInfo_.end(); ++it, ++i) {
digis->SetBinContent(i, it->digiOccupancy);
digisWithMasking->SetBinContent(i, it->digiOccupancyWithMasking);
HVoff->SetBinContent(i, it->HVoff);
// Store only with seconds precision
coral::TimeStamp coralTime(cond::time::to_boost(it->time));
// N.B. we add 1 hour to the coralTime because it is the conversion from posix_time which is non-adjusted.
// The shift of +1 gives the CERN time zone.
TDatime date1(coralTime.year(),
coralTime.month(),
coralTime.day(),
coralTime.hour() + 1,
coralTime.minute(),
coralTime.second());
// timeArray[i-1] = date1.Convert();
time->SetBinContent(i, date1.Convert());
}
// TCanvas * canvas = new TCanvas("digiAndHVvsTimeCanvas", "digi and HVoff vs time", 1000, 800);
// canvas->Draw();
// canvas->cd();
digis->Draw();
digis->Write();
digisWithMasking->Draw("same");
digisWithMasking->Write();
HVoff->Draw("same");
HVoff->SetLineColor(2);
HVoff->Write();
time->Draw("same");
time->Write();
// TGraph * digisGraph = buildGraph(digis, timeArray);
// digisGraph->Write("digisGraph");
// TGraph * digisWithMaskingGraph = buildGraph(digisWithMasking, timeArray);
// digisWithMaskingGraph->Write("digisWithMaskingGraph");
// TGraph * HVoffGraph = buildGraph(HVoff, timeArray);
// HVoffGraph->Write("HVoffGraph");
outputFile->Write();
outputFile->Close();
// delete[] timeArray;
}
// ------------ method called to for each event ------------
void SyncDCSO2O::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
using namespace edm;
const SiStripDetVOff* detVOff = &iSetup.getData(dcsToken_);
std::vector<uint32_t> detIds;
detVOff->getDetIds(detIds);
std::cout << "Number of DetIds with HV off = " << detIds.size() << std::endl;
// stringstream ss;
// detVOff->printSummary(ss);
// std::cout << ss.str() << std::endl;
// double seconds = (iEvent.time().value() >> 32);
std::cout << "event time = " << iEvent.time().value() << std::endl;
// std::cout << "event time in seconds = " << seconds << std::endl;
coral::TimeStamp coralTime(cond::time::to_boost(iEvent.time().value()));
std::cout << "year = " << coralTime.year() << ", month = " << coralTime.month() << ", day = " << coralTime.day();
// N.B. we add 1 hour to the coralTime because it is the conversion from posix_time which is non-adjusted.
// The shift of +1 gives the CERN time zone.
std::cout << ", hour = " << coralTime.hour() + 1 << ", minute = " << coralTime.minute()
<< ", second = " << coralTime.second();
std::cout << ", nanosecond = " << coralTime.nanosecond() << std::endl;
getDigis(iEvent);
if (!(digiDetsetVector_[0].isValid()))
std::cout << "NOT VALID DIGI COLLECTION 0" << std::endl;
else {
edm::DetSetVector<SiStripDigi>::const_iterator it = digiDetsetVector_[0]->begin();
unsigned int totDigis = 0;
unsigned int totDigisWithMasking = 0;
for (; it != digiDetsetVector_[0]->end(); ++it) {
totDigis += it->size();
// Compute also number of digis masking detIds off according to DetVOff
if (!(detVOff->IsModuleHVOff(it->detId()))) {
totDigisWithMasking += it->size();
}
}
std::cout << "digis = " << totDigis << std::endl;
timeInfo_.push_back(TimeInfo(iEvent.time().value(), totDigis, totDigisWithMasking, detVOff->getHVoffCounts()));
}
}
void SyncDCSO2O::getDigis(const edm::Event& iEvent) {
using namespace edm;
int icoll = 0;
Parameters::iterator itDigiProducersList = digiProducersList_.begin();
for (; itDigiProducersList != digiProducersList_.end(); ++itDigiProducersList) {
std::string digiProducer = itDigiProducersList->getParameter<std::string>("DigiProducer");
std::string digiLabel = itDigiProducersList->getParameter<std::string>("DigiLabel");
// std::cout << "Reading digi for " << digiProducer << " with label: " << digiLabel << std::endl;
iEvent.getByLabel(digiProducer, digiLabel, digiDetsetVector_[icoll]);
icoll++;
}
}
/// Build TGraphs with quantity vs time
TGraph* SyncDCSO2O::buildGraph(TH1F* histo, Float_t* timeArray) {
unsigned int arraySize = histo->GetNbinsX();
// Note that the array reported has [0] = underflow and [Nbins+1] = overflow
Float_t* valueArray = (histo->GetArray()) + 1;
TGraph* graph = new TGraph(arraySize, valueArray, timeArray);
graph->Draw("A*");
graph->GetXaxis()->SetTimeDisplay(1);
graph->GetXaxis()->SetLabelOffset(0.02);
graph->GetXaxis()->SetTimeFormat("#splitline{ %d}{%H:%M}");
graph->GetXaxis()->SetTimeOffset(0, "gmt");
graph->GetYaxis()->SetRangeUser(0, 16000);
graph->GetXaxis()->SetTitle("day/hour");
graph->GetXaxis()->SetTitleSize(0.03);
graph->GetXaxis()->SetTitleColor(kBlack);
graph->GetXaxis()->SetTitleOffset(1.80);
graph->GetYaxis()->SetTitle("number of digis");
graph->GetYaxis()->SetTitleSize(0.03);
graph->GetYaxis()->SetTitleColor(kBlack);
graph->GetYaxis()->SetTitleOffset(1.80);
graph->SetTitle();
return graph;
}
// ------------ method called once each job just before starting event loop ------------
void SyncDCSO2O::beginJob() {}
// ------------ method called once each job just after ending the event loop ------------
void SyncDCSO2O::endJob() {}
|