|
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
|
// -*- C++ -*-
//
// Package: SiStripMonitorRawData
// Class: SiStripMonitorRawData
//
/**\class SiStripMonitorDigi SiStripMonitorDigi.cc
DQM/SiStripMonitorDigi/src/SiStripMonitorDigi.cc
Description: <one line class summary>
Implementation:
<Notes on implementation>
*/
//
// Original Author: Suchandra Dutta
// Created: Fri June 1 17:00:00 CET 2007
//
//
#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
#include "DQM/SiStripMonitorPedestals/interface/SiStripMonitorRawData.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/SiStripDigi/interface/SiStripRawDigi.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
// std
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <numeric>
SiStripMonitorRawData::SiStripMonitorRawData(edm::ParameterSet const &iConfig)
: BadFedNumber(nullptr),
dqmStore_(edm::Service<DQMStore>().operator->()),
conf_(iConfig),
detCablingToken_(esConsumes<edm::Transition::BeginRun>())
{
// retrieve producer name of input StripDigiCollection
std::string digiProducer = conf_.getParameter<std::string>("DigiProducer");
std::string digiType = "VirginRaw";
digiToken_ = consumes<edm::DetSetVector<SiStripRawDigi>>(edm::InputTag(digiProducer, digiType));
edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData "
<< " Constructing....... ";
}
SiStripMonitorRawData::~SiStripMonitorRawData() {
edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData "
<< " Destructing....... ";
}
//
// -- BeginRun
//
void SiStripMonitorRawData::bookHistograms(DQMStore::IBooker &ibooker,
const edm::Run &run,
const edm::EventSetup &eSetup) {
if (BadFedNumber)
BadFedNumber->Reset();
if (detCablingWatcher_.check(eSetup)) {
SelectedDetIds.clear();
const auto &detcabling = eSetup.getData(detCablingToken_);
detcabling.addActiveDetectorsRawIds(SelectedDetIds);
edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::bookHistograms: "
<< " Creating MEs for new Cabling ";
ibooker.setCurrentFolder("Track/GlobalParameter");
if (!BadFedNumber) {
BadFedNumber =
ibooker.book1D("FaultyFedNumberAndChannel", "Faulty Fed Id and Channel and Numbers", 60000, 0.5, 600.5);
BadFedNumber->setAxisTitle("Fed Id and Channel numbers", 1);
}
}
}
// ------------ method called to produce the data ------------
void SiStripMonitorRawData::analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) {
edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::analyze: Run " << iEvent.id().run() << " Event "
<< iEvent.id().event();
const auto &detcabling = iSetup.getData(detCablingToken_);
// get DigiCollection object from Event
edm::Handle<edm::DetSetVector<SiStripRawDigi>> digi_collection;
iEvent.getByToken(digiToken_, digi_collection);
for (std::vector<uint32_t>::const_iterator idetid = SelectedDetIds.begin(), iEnd = SelectedDetIds.end();
idetid != iEnd;
++idetid) {
std::vector<edm::DetSet<SiStripRawDigi>>::const_iterator digis = digi_collection->find((*idetid));
if (digis == digi_collection->end() || digis->data.empty() || digis->data.size() > 768) {
std::vector<const FedChannelConnection *> fed_conns = detcabling.getConnections((*idetid));
for (unsigned int k = 0; k < fed_conns.size(); k++) {
if (fed_conns[k] && fed_conns[k]->isConnected()) {
float fed_id = fed_conns[k]->fedId() + 0.01 * fed_conns[k]->fedCh();
BadFedNumber->Fill(fed_id);
}
}
continue;
}
}
}
//
// -- End Run
//
void SiStripMonitorRawData::dqmEndRun(edm::Run const &run, edm::EventSetup const &eSetup) {
bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
if (outputMEsInRootFile) {
dqmStore_->save(outputFileName);
}
}
//
// -- End Job
//
void SiStripMonitorRawData::endJob(void) {
edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::EndJob: "
<< " Finishing!! ";
}
DEFINE_FWK_MODULE(SiStripMonitorRawData);
|