|
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
|
/**_________________________________________________________________
class: AlcaPCCEventProducer.cc
authors: Sam Higginbotham (shigginb@cern.ch), Chris Palmer (capalmer@cern.ch), Attila Radl (attila.radl@cern.ch)
________________________________________________________________**/
// C++ standard
#include <string>
// CMS
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Luminosity/interface/PixelClusterCountsInEvent.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
#include "FWCore/Framework/interface/ConsumesCollector.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/LuminosityBlock.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "TMath.h"
//The class
class AlcaPCCEventProducer : public edm::global::EDProducer<> {
public:
explicit AlcaPCCEventProducer(const edm::ParameterSet&);
~AlcaPCCEventProducer() override = default;
void produce(edm::StreamID id, edm::Event& e, edm::EventSetup const& c) const final;
static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
private:
const edm::InputTag pixelClusterLabel_;
const std::string trigstring_; //specifies the trigger Rand or ZeroBias
const bool savePerROCInfo_; // save per ROC data (important for the special fills)
const edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > pixelToken_;
static constexpr int rowsperroc = 52;
static constexpr int colsperroc = 80;
static constexpr int nROCcolumns = 8;
};
//--------------------------------------------------------------------------------------------------
AlcaPCCEventProducer::AlcaPCCEventProducer(const edm::ParameterSet& iConfig)
: pixelClusterLabel_(iConfig.getParameter<edm::InputTag>("pixelClusterLabel")),
trigstring_(iConfig.getUntrackedParameter<std::string>("trigstring", "alcaPCCEvent")),
savePerROCInfo_(iConfig.getParameter<bool>("savePerROCInfo")),
pixelToken_(consumes<edmNew::DetSetVector<SiPixelCluster> >(pixelClusterLabel_)) {
produces<reco::PixelClusterCountsInEvent, edm::Transition::Event>(trigstring_);
}
//--------------------------------------------------------------------------------------------------
void AlcaPCCEventProducer::produce(edm::StreamID id, edm::Event& iEvent, edm::EventSetup const& iSetup) const {
std::unique_ptr<reco::PixelClusterCountsInEvent> thePCCob = std::make_unique<reco::PixelClusterCountsInEvent>();
unsigned int bx = iEvent.bunchCrossing();
//Looping over the clusters and adding the counts up
edm::Handle<edmNew::DetSetVector<SiPixelCluster> > hClusterColl;
iEvent.getByToken(pixelToken_, hClusterColl);
const edmNew::DetSetVector<SiPixelCluster>& clustColl = *(hClusterColl.product());
// ----------------------------------------------------------------------
// -- Clusters without tracks
for (auto const& mod : clustColl) {
if (mod.empty()) {
continue;
}
DetId detId = mod.id();
if (savePerROCInfo_) {
// Iterate over Clusters in module to fill per ROC histogram
for (auto const& cluster : mod) {
for (int i = 0; i < cluster.size(); ++i) {
const auto pix = cluster.pixel(i);
// TODO: add roc threshold to config if(di.adc > fRocThreshold_) {
if (pix.adc > 0) {
int irow = pix.x / rowsperroc; /* constant column direction is along x-axis */
int icol = pix.y / colsperroc; /* constant row direction is along y-axis */
/* generate the folling roc index that is going to map with ROC id as
8 9 10 11 12 13 14 15
0 1 2 3 4 5 6 7 */
int key = icol + irow * nROCcolumns;
thePCCob->incrementRoc(((detId << 7) + key), 1);
}
}
}
}
int nCluster = mod.size();
thePCCob->increment(detId(), nCluster);
thePCCob->setbxID(bx);
}
iEvent.put(std::move(thePCCob), std::string(trigstring_));
}
//--------------------------------------------------------------------------------------------------
void AlcaPCCEventProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription evtParamDesc;
evtParamDesc.add<edm::InputTag>("pixelClusterLabel", edm::InputTag("siPixelClustersForLumi"));
evtParamDesc.addUntracked<std::string>("trigstring", "alcaPCCEvent");
evtParamDesc.add<bool>("savePerROCInfo", true);
descriptions.add("alcaPCCEventProducer", evtParamDesc);
}
DEFINE_FWK_MODULE(AlcaPCCEventProducer);
|