69 theTrackerGeometry = tGeometryHandle.
product();
75 if (!Phase2TrackerDigiTTStubHandle.
isValid())
78 for (inputIter = Phase2TrackerDigiTTStubHandle->begin(); inputIter != Phase2TrackerDigiTTStubHandle->end();
80 for (contentIter = inputIter->
begin(); contentIter != inputIter->
end(); ++contentIter) {
87 DetId detIdStub = theTrackerGeometry->
idToDet((tempStubRef->getClusterRef(0))->getDetId())->geographicalId();
90 double displStub = tempStubRef->getTriggerDisplacement();
91 double offsetStub = tempStubRef->getTriggerOffset();
94 MeasurementPoint mp = (tempStubRef->getClusterRef(0))->findAverageLocalCoordinates();
95 const GeomDet *theGeomDet = theTrackerGeometry->
idToDet(detIdStub);
118 if (posStub.
z() > 0) {
141 const int numDiscs = 5;
145 HistoName =
"Stub_Barrel_XY";
158 HistoName =
"Stub_Endcap_Fw_XY";
171 HistoName =
"Stub_Endcap_Bw_XY";
185 HistoName =
"Stub_RZ";
200 HistoName =
"Stub_Eta";
211 HistoName =
"Stub_Phi";
222 HistoName =
"Stub_R";
234 HistoName =
"NStubs_Barrel";
245 HistoName =
"NStubs_Endcap_Disc";
255 HistoName =
"NStubs_Endcap_Disc_Fw";
265 HistoName =
"NStubs_Endcap_Disc_Bw";
275 HistoName =
"NStubs_Endcap_Ring";
284 for (
int i = 0;
i < numDiscs;
i++) {
285 HistoName =
"NStubs_Disc+" + std::to_string(
i + 1);
296 for (
int i = 0;
i < numDiscs;
i++) {
297 HistoName =
"NStubs_Disc-" + std::to_string(
i + 1);
314 HistoName =
"Stub_Width_Barrel";
326 HistoName =
"Stub_Width_Endcap_Disc";
330 psTTStub_ECDisc_2D.getParameter<
double>(
"xmin"),
331 psTTStub_ECDisc_2D.getParameter<
double>(
"xmax"),
332 psTTStub_ECDisc_2D.getParameter<int32_t>(
"Nbinsy"),
333 psTTStub_ECDisc_2D.getParameter<
double>(
"ymin"),
334 psTTStub_ECDisc_2D.getParameter<
double>(
"ymax"));
338 HistoName =
"Stub_Width_Endcap_Ring";
342 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
343 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
344 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
345 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
346 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
350 for (
int i = 0;
i < numDiscs;
i++) {
351 HistoName =
"Stub_Width_Disc+" + std::to_string(
i + 1);
355 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
356 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
357 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
358 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
359 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
364 for (
int i = 0;
i < numDiscs;
i++) {
365 HistoName =
"Stub_Width_Disc-" + std::to_string(
i + 1);
369 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
370 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
371 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
372 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
373 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
379 HistoName =
"Stub_Offset_Barrel";
391 HistoName =
"Stub_Offset_Endcap_Disc";
395 psTTStub_ECDisc_2D.getParameter<
double>(
"xmin"),
396 psTTStub_ECDisc_2D.getParameter<
double>(
"xmax"),
397 psTTStub_ECDisc_2D.getParameter<int32_t>(
"Nbinsy"),
398 psTTStub_ECDisc_2D.getParameter<
double>(
"ymin"),
399 psTTStub_ECDisc_2D.getParameter<
double>(
"ymax"));
403 HistoName =
"Stub_Offset_Endcap_Ring";
407 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
408 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
409 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
410 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
411 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
415 for (
int i = 0;
i < numDiscs;
i++) {
416 HistoName =
"Stub_Offset_Disc+" + std::to_string(
i + 1);
420 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
421 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
422 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
423 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
424 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
429 for (
int i = 0;
i < numDiscs;
i++) {
430 HistoName =
"Stub_Offset_Disc-" + std::to_string(
i + 1);
434 psTTStub_ECRing_2D.getParameter<
double>(
"xmin"),
435 psTTStub_ECRing_2D.getParameter<
double>(
"xmax"),
436 psTTStub_ECRing_2D.getParameter<int32_t>(
"Nbinsy"),
437 psTTStub_ECRing_2D.getParameter<
double>(
"ymin"),
438 psTTStub_ECRing_2D.getParameter<
double>(
"ymax"));
edm::EDGetTokenT< edmNew::DetSetVector< TTStub< Ref_Phase2TrackerDigi_ > > > tagTTStubsToken_
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
T getParameter(std::string const &) const
edm::Ref< typename HandleT::element_type, typename HandleT::element_type::value_type::value_type > makeRefTo(const HandleT &iHandle, typename HandleT::element_type::value_type::const_iterator itIter)
const_iterator end(bool update=false) const
MonitorElement * Stub_Endcap_Bw_XY
unsigned int tidRing(const DetId &id) const
std::string topFolderName_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
virtual const Topology & topology() const
Geom::Phi< T > phi() const
MonitorElement * Stub_Endcap_Ring_O_Fw[5]
MonitorElement * Stub_Endcap_Fw_XY
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
MonitorElement * Stub_Barrel_O
const Plane & surface() const
The nominal surface of the GeomDet.
MonitorElement * Stub_Barrel_XY
~OuterTrackerMonitorTTStub() override
void analyze(const edm::Event &, const edm::EventSetup &) override
MonitorElement * Stub_Endcap_Ring_Bw[5]
MonitorElement * Stub_Endcap_Disc_Fw
MonitorElement * Stub_Eta
#define DEFINE_FWK_MODULE(type)
OuterTrackerMonitorTTStub(const edm::ParameterSet &)
MonitorElement * Stub_Endcap_Ring_O_Bw[5]
void setCurrentFolder(std::string const &fullpath)
MonitorElement * Stub_Endcap_Disc
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
MonitorElement * book1D(Args &&...args)
MonitorElement * Stub_Endcap_Disc_Bw
Class to store the L1 Track Trigger stubs.
MonitorElement * Stub_Endcap_Ring
MonitorElement * Stub_Endcap_Ring_O
MonitorElement * book2D(Args &&...args)
MonitorElement * Stub_Endcap_Disc_O
virtual LocalPoint localPosition(const MeasurementPoint &) const =0
MonitorElement * Stub_Endcap_Disc_W
unsigned int layer(const DetId &id) const
MonitorElement * Stub_Barrel
MonitorElement * Stub_Endcap_Ring_Fw[5]
MonitorElement * Stub_Phi
const TrackerGeomDet * idToDet(DetId) const override
MonitorElement * Stub_Endcap_Ring_W_Bw[5]
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
T const * product() const
MonitorElement * Stub_Barrel_W
const_iterator begin(bool update=false) const
MonitorElement * Stub_Endcap_Ring_W
MonitorElement * Stub_Endcap_Ring_W_Fw[5]