|
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
|
// -*- C++ -*-
//
// Package: SiPixelPhase1Digis
// Class: SiPixelPhase1Digis
//
// Original Author: Marcel Schneider
// C++ stuff
#include <iostream>
// CMSSW stuff
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
// DQM Stuff
#include "DQMServices/Core/interface/DQMStore.h"
// Input data stuff
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiPixelDigi/interface/PixelDigi.h"
// PixelDQM Framework
#include "DQM/SiPixelPhase1Common/interface/SiPixelPhase1Base.h"
namespace {
class SiPixelPhase1Digis final : public SiPixelPhase1Base {
// List of quantities to be plotted.
enum {
ADC, // digi ADC readouts
NDIGIS, // number of digis per event and module
NDIGISINCLUSIVE, //Total number of digis in BPix and FPix
NDIGIS_FED, // number of digis per event and FED
NDIGIS_FEDtrend, // number of digis per event and FED
EVENT, // event frequency
MAP, // digi hitmap per module
OCCUPANCY, // like map but coarser
MAX_HIST // a sentinel that gives the number of quantities (not a plot).
};
public:
explicit SiPixelPhase1Digis(const edm::ParameterSet& conf);
void analyze(const edm::Event&, const edm::EventSetup&) override;
private:
edm::EDGetTokenT<edm::DetSetVector<PixelDigi>> srcToken_;
};
SiPixelPhase1Digis::SiPixelPhase1Digis(const edm::ParameterSet& iConfig) : SiPixelPhase1Base(iConfig) {
srcToken_ = consumes<edm::DetSetVector<PixelDigi>>(iConfig.getParameter<edm::InputTag>("src"));
}
void SiPixelPhase1Digis::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
if (!checktrigger(iEvent, iSetup, DCS))
return;
edm::Handle<edm::DetSetVector<PixelDigi>> input;
iEvent.getByToken(srcToken_, input);
if (!input.isValid())
return;
bool hasDigis = false;
edm::DetSetVector<PixelDigi>::const_iterator it;
for (it = input->begin(); it != input->end(); ++it) {
for (PixelDigi const& digi : *it) {
hasDigis = true;
histo[ADC].fill((double)digi.adc(), DetId(it->detId()), &iEvent, digi.column(), digi.row());
histo[MAP].fill(DetId(it->detId()), &iEvent, digi.column(), digi.row());
histo[OCCUPANCY].fill(DetId(it->detId()), &iEvent, digi.column(), digi.row());
histo[NDIGIS].fill(DetId(it->detId()), &iEvent); // count
histo[NDIGISINCLUSIVE].fill(DetId(it->detId()), &iEvent); // count
histo[NDIGIS_FED].fill(DetId(it->detId()), &iEvent);
histo[NDIGIS_FEDtrend].fill(DetId(it->detId()), &iEvent);
}
}
if (hasDigis)
histo[EVENT].fill(DetId(0), &iEvent);
histo[NDIGIS].executePerEventHarvesting(&iEvent);
histo[NDIGISINCLUSIVE].executePerEventHarvesting(&iEvent);
histo[NDIGIS_FED].executePerEventHarvesting(&iEvent);
histo[NDIGIS_FEDtrend].executePerEventHarvesting(&iEvent);
}
} //namespace
DEFINE_FWK_MODULE(SiPixelPhase1Digis);
|