|
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
|
// -*- C++ -*-
//
// Package: SiPixelPhase1Clusters
// Class: SiPixelPhase1Clusters
//
// Original Author: Marcel Schneider
#include "DQM/SiPixelPhase1Common/interface/SiPixelPhase1Base.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
#include "Geometry/CommonTopologies/interface/PixelTopology.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
namespace {
class SiPixelPhase1Clusters final : public SiPixelPhase1Base {
enum {
CHARGE,
BIGPIXELCHARGE,
NOTBIGPIXELCHARGE,
SIZE,
SIZEX,
SIZEY,
NCLUSTERS,
NCLUSTERSINCLUSIVE,
EVENTRATE,
POSITION_B,
POSITION_F,
POSITION_XZ,
POSITION_YZ,
SIZE_VS_ETA,
READOUT_CHARGE,
READOUT_NCLUSTERS,
PIXEL_TO_STRIP_RATIO
};
public:
explicit SiPixelPhase1Clusters(const edm::ParameterSet& conf);
void analyze(const edm::Event&, const edm::EventSetup&) override;
private:
edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster>> pixelSrcToken_;
edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster>> stripSrcToken_;
edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> trackerGeometryToken_;
};
SiPixelPhase1Clusters::SiPixelPhase1Clusters(const edm::ParameterSet& iConfig)
: SiPixelPhase1Base(iConfig), trackerGeometryToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()) {
pixelSrcToken_ = consumes<edmNew::DetSetVector<SiPixelCluster>>(iConfig.getParameter<edm::InputTag>("pixelSrc"));
stripSrcToken_ = consumes<edmNew::DetSetVector<SiStripCluster>>(iConfig.getParameter<edm::InputTag>("stripSrc"));
}
void SiPixelPhase1Clusters::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
if (!checktrigger(iEvent, iSetup, DCS))
return;
edm::Handle<edmNew::DetSetVector<SiPixelCluster>> inputPixel;
iEvent.getByToken(pixelSrcToken_, inputPixel);
if (!inputPixel.isValid())
return;
edm::Handle<edmNew::DetSetVector<SiStripCluster>> inputStrip;
iEvent.getByToken(stripSrcToken_, inputStrip);
if (inputStrip.isValid()) {
if (!inputStrip.product()->data().empty()) {
histo[PIXEL_TO_STRIP_RATIO].fill(
(double)inputPixel.product()->data().size() / (double)inputStrip.product()->data().size(),
DetId(0),
&iEvent);
}
}
bool hasClusters = false;
const TrackerGeometry& tracker = iSetup.getData(trackerGeometryToken_);
edmNew::DetSetVector<SiPixelCluster>::const_iterator it;
for (it = inputPixel->begin(); it != inputPixel->end(); ++it) {
auto id = DetId(it->detId());
const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(tracker.idToDet(id));
const PixelTopology& topol = theGeomDet->specificTopology();
for (SiPixelCluster const& cluster : *it) {
int row = cluster.x() - 0.5, col = cluster.y() - 0.5;
const std::vector<SiPixelCluster::Pixel> pixelsVec = cluster.pixels();
for (unsigned int i = 0; i < pixelsVec.size(); ++i) {
float pixx = pixelsVec[i].x; // index as float=iteger, row index
float pixy = pixelsVec[i].y; // same, col index
bool bigInX = topol.isItBigPixelInX(int(pixx));
bool bigInY = topol.isItBigPixelInY(int(pixy));
float pixel_charge = pixelsVec[i].adc;
if (bigInX == true || bigInY == true) {
histo[BIGPIXELCHARGE].fill(pixel_charge, id, &iEvent, col, row);
}
else {
histo[NOTBIGPIXELCHARGE].fill(pixel_charge, id, &iEvent, col, row);
}
}
histo[READOUT_CHARGE].fill(double(cluster.charge()), id, &iEvent, col, row);
histo[CHARGE].fill(double(cluster.charge()), id, &iEvent, col, row);
histo[SIZE].fill(double(cluster.size()), id, &iEvent, col, row);
histo[SIZEX].fill(double(cluster.sizeX()), id, &iEvent, col, row);
histo[SIZEY].fill(double(cluster.sizeY()), id, &iEvent, col, row);
histo[NCLUSTERS].fill(id, &iEvent, col, row);
histo[NCLUSTERSINCLUSIVE].fill(id, &iEvent);
hasClusters = true;
if (cluster.size() > 1) {
histo[READOUT_NCLUSTERS].fill(id, &iEvent);
}
LocalPoint clustlp = topol.localPosition(MeasurementPoint(cluster.x(), cluster.y()));
GlobalPoint clustgp = theGeomDet->surface().toGlobal(clustlp);
histo[POSITION_B].fill(clustgp.z(), clustgp.phi(), id, &iEvent);
histo[POSITION_F].fill(clustgp.x(), clustgp.y(), id, &iEvent);
histo[POSITION_XZ].fill(clustgp.x(), clustgp.z(), id, &iEvent);
histo[POSITION_YZ].fill(clustgp.y(), clustgp.z(), id, &iEvent);
histo[SIZE_VS_ETA].fill(clustgp.eta(), cluster.sizeY(), id, &iEvent);
}
}
if (hasClusters)
histo[EVENTRATE].fill(DetId(0), &iEvent);
histo[NCLUSTERS].executePerEventHarvesting(&iEvent);
histo[READOUT_NCLUSTERS].executePerEventHarvesting(&iEvent);
histo[NCLUSTERSINCLUSIVE].executePerEventHarvesting(&iEvent);
}
} // namespace
DEFINE_FWK_MODULE(SiPixelPhase1Clusters);
|