#include <OuterTrackerMonitorTTStub.h>

Inheritance diagram for OuterTrackerMonitorTTStub:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	OuterTrackerMonitorTTStub (const edm::ParameterSet &)

	~OuterTrackerMonitorTTStub () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Public Attributes
MonitorElement *	Stub_Barrel = 0

MonitorElement *	Stub_Barrel_O = 0

MonitorElement *	Stub_Barrel_W = 0

MonitorElement *	Stub_Barrel_XY = 0

MonitorElement *	Stub_Endcap_Bw_XY = 0

MonitorElement *	Stub_Endcap_Disc = 0

MonitorElement *	Stub_Endcap_Disc_Bw = 0

MonitorElement *	Stub_Endcap_Disc_Fw = 0

MonitorElement *	Stub_Endcap_Disc_O = 0

MonitorElement *	Stub_Endcap_Disc_W = 0

MonitorElement *	Stub_Endcap_Fw_XY = 0

MonitorElement *	Stub_Endcap_Ring = 0

MonitorElement *	Stub_Endcap_Ring_Bw [5] = { 0 , 0 , 0 , 0 , 0 }

MonitorElement *	Stub_Endcap_Ring_Fw [5] = { 0 , 0 , 0 , 0 , 0 }

MonitorElement *	Stub_Endcap_Ring_O = 0

MonitorElement *	Stub_Endcap_Ring_O_Bw [5] = { 0 , 0 , 0 , 0 , 0 }

MonitorElement *	Stub_Endcap_Ring_O_Fw [5] = { 0 , 0 , 0 , 0 , 0 }

MonitorElement *	Stub_Endcap_Ring_W = 0

MonitorElement *	Stub_Endcap_Ring_W_Bw [5]

MonitorElement *	Stub_Endcap_Ring_W_Fw [5]

MonitorElement *	Stub_Eta = 0

MonitorElement *	Stub_Phi = 0

MonitorElement *	Stub_R = 0

MonitorElement *	Stub_RZ = 0

Private Attributes
edm::ParameterSet	conf_

edm::EDGetTokenT< edmNew::DetSetVector< TTStub< Ref_Phase2TrackerDigi_ > > >	tagTTStubsToken_

std::string	topFolderName_

Detailed Description

Definition at line 20 of file OuterTrackerMonitorTTStub.h.

Constructor & Destructor Documentation

OuterTrackerMonitorTTStub::OuterTrackerMonitorTTStub ( const edm::ParameterSet & iConfig )

explicit

Definition at line 40 of file OuterTrackerMonitorTTStub.cc.

References conf_, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, tagTTStubsToken_, and topFolderName_.

                                                                                    : conf_(iConfig) {
   // now do what ever initialization is needed
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
   tagTTStubsToken_ =
       consumes<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>>(conf_.getParameter<edm::InputTag>("TTStubs"));
 }

OuterTrackerMonitorTTStub::~OuterTrackerMonitorTTStub ( )

override

Definition at line 47 of file OuterTrackerMonitorTTStub.cc.

                                                       {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void OuterTrackerMonitorTTStub::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Track Trigger Stubs

Geometry

Loop over input Stubs

Make reference stub

Get det ID (place of the stub)

Get trigger displacement/offset

Define position stub by position inner cluster

Definition at line 55 of file OuterTrackerMonitorTTStub.cc.

References edmNew::DetSetVector< T >::begin(), TtSemiLepEvtBuilder_cfi::disc, edmNew::DetSetVector< T >::end(), PV3DBase< T, PVType, FrameType >::eta(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), TrackerGeometry::idToDet(), edm::HandleBase::isValid(), TrackerTopology::layer(), Topology::localPosition(), edmNew::makeRefTo(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), edm::ESHandle< T >::product(), relativeConstraints::ring, Stub_Barrel, Stub_Barrel_O, Stub_Barrel_W, Stub_Barrel_XY, Stub_Endcap_Bw_XY, Stub_Endcap_Disc, Stub_Endcap_Disc_Bw, Stub_Endcap_Disc_Fw, Stub_Endcap_Disc_O, Stub_Endcap_Disc_W, Stub_Endcap_Fw_XY, Stub_Endcap_Ring, Stub_Endcap_Ring_Bw, Stub_Endcap_Ring_Fw, Stub_Endcap_Ring_O, Stub_Endcap_Ring_O_Bw, Stub_Endcap_Ring_O_Fw, Stub_Endcap_Ring_W, Stub_Endcap_Ring_W_Bw, Stub_Endcap_Ring_W_Fw, Stub_Eta, Stub_Phi, Stub_R, Stub_RZ, DetId::subdetId(), GeomDet::surface(), tagTTStubsToken_, StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, Surface::toGlobal(), GeomDet::topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                            {
   edm::Handle<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>> Phase2TrackerDigiTTStubHandle;
   iEvent.getByToken(tagTTStubsToken_, Phase2TrackerDigiTTStubHandle);
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   const TrackerTopology *tTopo;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   tTopo = tTopoHandle.product();
 
   edm::ESHandle<TrackerGeometry> tGeometryHandle;
   const TrackerGeometry *theTrackerGeometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(tGeometryHandle);
   theTrackerGeometry = tGeometryHandle.product();
 
   typename edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>::const_iterator inputIter;
   typename edmNew::DetSet<TTStub<Ref_Phase2TrackerDigi_>>::const_iterator contentIter;
   // Adding protection
   if (!Phase2TrackerDigiTTStubHandle.isValid())
     return;
 
   for (inputIter = Phase2TrackerDigiTTStubHandle->begin(); inputIter != Phase2TrackerDigiTTStubHandle->end();
        ++inputIter) {
     for (contentIter = inputIter->begin(); contentIter != inputIter->end(); ++contentIter) {
       edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>, TTStub<Ref_Phase2TrackerDigi_>> tempStubRef =
           edmNew::makeRefTo(Phase2TrackerDigiTTStubHandle, contentIter);
 
       //  tempStubRef->getDetId() gives the stackDetId, not rawId
       DetId detIdStub = theTrackerGeometry->idToDet((tempStubRef->getClusterRef(0))->getDetId())->geographicalId();
 
       double displStub = tempStubRef->getTriggerDisplacement();
       double offsetStub = tempStubRef->getTriggerOffset();
 
       MeasurementPoint mp = (tempStubRef->getClusterRef(0))->findAverageLocalCoordinates();
       const GeomDet *theGeomDet = theTrackerGeometry->idToDet(detIdStub);
       Global3DPoint posStub = theGeomDet->surface().toGlobal(theGeomDet->topology().localPosition(mp));
 
       Stub_Eta->Fill(posStub.eta());
       Stub_Phi->Fill(posStub.phi());
       Stub_R->Fill(posStub.perp());
       Stub_RZ->Fill(posStub.z(), posStub.perp());
 
       if (detIdStub.subdetId() == static_cast<int>(StripSubdetector::TOB)) {  // Phase 2 Outer Tracker Barrel
         Stub_Barrel->Fill(tTopo->layer(detIdStub));
         Stub_Barrel_XY->Fill(posStub.x(), posStub.y());
         Stub_Barrel_W->Fill(tTopo->layer(detIdStub), displStub - offsetStub);
         Stub_Barrel_O->Fill(tTopo->layer(detIdStub), offsetStub);
       } else if (detIdStub.subdetId() == static_cast<int>(StripSubdetector::TID)) {  // Phase 2 Outer Tracker Endcap
         int disc = tTopo->layer(detIdStub);                                          // returns wheel
         int ring = tTopo->tidRing(detIdStub);
         Stub_Endcap_Disc->Fill(disc);
         Stub_Endcap_Ring->Fill(ring);
         Stub_Endcap_Disc_W->Fill(disc, displStub - offsetStub);
         Stub_Endcap_Ring_W->Fill(ring, displStub - offsetStub);
         Stub_Endcap_Disc_O->Fill(disc, offsetStub);
         Stub_Endcap_Ring_O->Fill(ring, offsetStub);
 
         if (posStub.z() > 0) {
           Stub_Endcap_Fw_XY->Fill(posStub.x(), posStub.y());
           Stub_Endcap_Disc_Fw->Fill(disc);
           Stub_Endcap_Ring_Fw[disc - 1]->Fill(ring);
           Stub_Endcap_Ring_W_Fw[disc - 1]->Fill(ring, displStub - offsetStub);
           Stub_Endcap_Ring_O_Fw[disc - 1]->Fill(ring, offsetStub);
         } else {
           Stub_Endcap_Bw_XY->Fill(posStub.x(), posStub.y());
           Stub_Endcap_Disc_Bw->Fill(disc);
           Stub_Endcap_Ring_Bw[disc - 1]->Fill(ring);
           Stub_Endcap_Ring_W_Bw[disc - 1]->Fill(ring, displStub - offsetStub);
           Stub_Endcap_Ring_O_Bw[disc - 1]->Fill(ring, offsetStub);
         }
       }
     }
   }
 }  // end of method

void OuterTrackerMonitorTTStub::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

override

Definition at line 137 of file OuterTrackerMonitorTTStub.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), conf_, DEFINE_FWK_MODULE, edm::ParameterSet::getParameter(), mps_fire::i, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, Stub_Barrel, Stub_Barrel_O, Stub_Barrel_W, Stub_Barrel_XY, Stub_Endcap_Bw_XY, Stub_Endcap_Disc, Stub_Endcap_Disc_Bw, Stub_Endcap_Disc_Fw, Stub_Endcap_Disc_O, Stub_Endcap_Disc_W, Stub_Endcap_Fw_XY, Stub_Endcap_Ring, Stub_Endcap_Ring_Bw, Stub_Endcap_Ring_Fw, Stub_Endcap_Ring_O, Stub_Endcap_Ring_O_Bw, Stub_Endcap_Ring_O_Fw, Stub_Endcap_Ring_W, Stub_Endcap_Ring_W_Bw, Stub_Endcap_Ring_W_Fw, Stub_Eta, Stub_Phi, Stub_R, Stub_RZ, and topFolderName_.

                                                                         {
   std::string HistoName;
   const int numDiscs = 5;
   iBooker.setCurrentFolder(topFolderName_ + "/Stubs/Position");
 
   edm::ParameterSet psTTStub_Barrel_XY = conf_.getParameter<edm::ParameterSet>("TH2TTStub_Position");
   HistoName = "Stub_Barrel_XY";
   Stub_Barrel_XY = iBooker.book2D(HistoName,
                                   HistoName,
                                   psTTStub_Barrel_XY.getParameter<int32_t>("Nbinsx"),
                                   psTTStub_Barrel_XY.getParameter<double>("xmin"),
                                   psTTStub_Barrel_XY.getParameter<double>("xmax"),
                                   psTTStub_Barrel_XY.getParameter<int32_t>("Nbinsy"),
                                   psTTStub_Barrel_XY.getParameter<double>("ymin"),
                                   psTTStub_Barrel_XY.getParameter<double>("ymax"));
   Stub_Barrel_XY->setAxisTitle("L1 Stub Barrel position x [cm]", 1);
   Stub_Barrel_XY->setAxisTitle("L1 Stub Barrel position y [cm]", 2);
 
   edm::ParameterSet psTTStub_Endcap_Fw_XY = conf_.getParameter<edm::ParameterSet>("TH2TTStub_Position");
   HistoName = "Stub_Endcap_Fw_XY";
   Stub_Endcap_Fw_XY = iBooker.book2D(HistoName,
                                      HistoName,
                                      psTTStub_Endcap_Fw_XY.getParameter<int32_t>("Nbinsx"),
                                      psTTStub_Endcap_Fw_XY.getParameter<double>("xmin"),
                                      psTTStub_Endcap_Fw_XY.getParameter<double>("xmax"),
                                      psTTStub_Endcap_Fw_XY.getParameter<int32_t>("Nbinsy"),
                                      psTTStub_Endcap_Fw_XY.getParameter<double>("ymin"),
                                      psTTStub_Endcap_Fw_XY.getParameter<double>("ymax"));
   Stub_Endcap_Fw_XY->setAxisTitle("L1 Stub Endcap position x [cm]", 1);
   Stub_Endcap_Fw_XY->setAxisTitle("L1 Stub Endcap position y [cm]", 2);
 
   edm::ParameterSet psTTStub_Endcap_Bw_XY = conf_.getParameter<edm::ParameterSet>("TH2TTStub_Position");
   HistoName = "Stub_Endcap_Bw_XY";
   Stub_Endcap_Bw_XY = iBooker.book2D(HistoName,
                                      HistoName,
                                      psTTStub_Endcap_Bw_XY.getParameter<int32_t>("Nbinsx"),
                                      psTTStub_Endcap_Bw_XY.getParameter<double>("xmin"),
                                      psTTStub_Endcap_Bw_XY.getParameter<double>("xmax"),
                                      psTTStub_Endcap_Bw_XY.getParameter<int32_t>("Nbinsy"),
                                      psTTStub_Endcap_Bw_XY.getParameter<double>("ymin"),
                                      psTTStub_Endcap_Bw_XY.getParameter<double>("ymax"));
   Stub_Endcap_Bw_XY->setAxisTitle("L1 Stub Endcap position x [cm]", 1);
   Stub_Endcap_Bw_XY->setAxisTitle("L1 Stub Endcap position y [cm]", 2);
 
   // TTStub #rho vs. z
   edm::ParameterSet psTTStub_RZ = conf_.getParameter<edm::ParameterSet>("TH2TTStub_RZ");
   HistoName = "Stub_RZ";
   Stub_RZ = iBooker.book2D(HistoName,
                            HistoName,
                            psTTStub_RZ.getParameter<int32_t>("Nbinsx"),
                            psTTStub_RZ.getParameter<double>("xmin"),
                            psTTStub_RZ.getParameter<double>("xmax"),
                            psTTStub_RZ.getParameter<int32_t>("Nbinsy"),
                            psTTStub_RZ.getParameter<double>("ymin"),
                            psTTStub_RZ.getParameter<double>("ymax"));
   Stub_RZ->setAxisTitle("L1 Stub position z [cm]", 1);
   Stub_RZ->setAxisTitle("L1 Stub position #rho [cm]", 2);
 
   iBooker.setCurrentFolder(topFolderName_ + "/Stubs");
   // TTStub eta
   edm::ParameterSet psTTStub_Eta = conf_.getParameter<edm::ParameterSet>("TH1TTStub_Eta");
   HistoName = "Stub_Eta";
   Stub_Eta = iBooker.book1D(HistoName,
                             HistoName,
                             psTTStub_Eta.getParameter<int32_t>("Nbinsx"),
                             psTTStub_Eta.getParameter<double>("xmin"),
                             psTTStub_Eta.getParameter<double>("xmax"));
   Stub_Eta->setAxisTitle("#eta", 1);
   Stub_Eta->setAxisTitle("# L1 Stubs ", 2);
 
   // TTStub phi
   edm::ParameterSet psTTStub_Phi = conf_.getParameter<edm::ParameterSet>("TH1TTStub_Phi");
   HistoName = "Stub_Phi";
   Stub_Phi = iBooker.book1D(HistoName,
                             HistoName,
                             psTTStub_Phi.getParameter<int32_t>("Nbinsx"),
                             psTTStub_Phi.getParameter<double>("xmin"),
                             psTTStub_Phi.getParameter<double>("xmax"));
   Stub_Phi->setAxisTitle("#phi", 1);
   Stub_Phi->setAxisTitle("# L1 Stubs ", 2);
 
   // TTStub R
   edm::ParameterSet psTTStub_R = conf_.getParameter<edm::ParameterSet>("TH1TTStub_R");
   HistoName = "Stub_R";
   Stub_R = iBooker.book1D(HistoName,
                           HistoName,
                           psTTStub_R.getParameter<int32_t>("Nbinsx"),
                           psTTStub_R.getParameter<double>("xmin"),
                           psTTStub_R.getParameter<double>("xmax"));
   Stub_R->setAxisTitle("R", 1);
   Stub_R->setAxisTitle("# L1 Stubs ", 2);
 
   iBooker.setCurrentFolder(topFolderName_ + "/Stubs/NStubs");
   // TTStub barrel stack
   edm::ParameterSet psTTStub_Barrel = conf_.getParameter<edm::ParameterSet>("TH1TTStub_Layers");
   HistoName = "NStubs_Barrel";
   Stub_Barrel = iBooker.book1D(HistoName,
                                HistoName,
                                psTTStub_Barrel.getParameter<int32_t>("Nbinsx"),
                                psTTStub_Barrel.getParameter<double>("xmin"),
                                psTTStub_Barrel.getParameter<double>("xmax"));
   Stub_Barrel->setAxisTitle("Barrel Layer", 1);
   Stub_Barrel->setAxisTitle("# L1 Stubs ", 2);
 
   // TTStub Endcap stack
   edm::ParameterSet psTTStub_ECDisc = conf_.getParameter<edm::ParameterSet>("TH1TTStub_Discs");
   HistoName = "NStubs_Endcap_Disc";
   Stub_Endcap_Disc = iBooker.book1D(HistoName,
                                     HistoName,
                                     psTTStub_ECDisc.getParameter<int32_t>("Nbinsx"),
                                     psTTStub_ECDisc.getParameter<double>("xmin"),
                                     psTTStub_ECDisc.getParameter<double>("xmax"));
   Stub_Endcap_Disc->setAxisTitle("Endcap Disc", 1);
   Stub_Endcap_Disc->setAxisTitle("# L1 Stubs ", 2);
 
   // TTStub Endcap stack
   HistoName = "NStubs_Endcap_Disc_Fw";
   Stub_Endcap_Disc_Fw = iBooker.book1D(HistoName,
                                        HistoName,
                                        psTTStub_ECDisc.getParameter<int32_t>("Nbinsx"),
                                        psTTStub_ECDisc.getParameter<double>("xmin"),
                                        psTTStub_ECDisc.getParameter<double>("xmax"));
   Stub_Endcap_Disc_Fw->setAxisTitle("Forward Endcap Disc", 1);
   Stub_Endcap_Disc_Fw->setAxisTitle("# L1 Stubs ", 2);
 
   // TTStub Endcap stack
   HistoName = "NStubs_Endcap_Disc_Bw";
   Stub_Endcap_Disc_Bw = iBooker.book1D(HistoName,
                                        HistoName,
                                        psTTStub_ECDisc.getParameter<int32_t>("Nbinsx"),
                                        psTTStub_ECDisc.getParameter<double>("xmin"),
                                        psTTStub_ECDisc.getParameter<double>("xmax"));
   Stub_Endcap_Disc_Bw->setAxisTitle("Backward Endcap Disc", 1);
   Stub_Endcap_Disc_Bw->setAxisTitle("# L1 Stubs ", 2);
 
   edm::ParameterSet psTTStub_ECRing = conf_.getParameter<edm::ParameterSet>("TH1TTStub_Rings");
   HistoName = "NStubs_Endcap_Ring";
   Stub_Endcap_Ring = iBooker.book1D(HistoName,
                                     HistoName,
                                     psTTStub_ECRing.getParameter<int32_t>("Nbinsx"),
                                     psTTStub_ECRing.getParameter<double>("xmin"),
                                     psTTStub_ECRing.getParameter<double>("xmax"));
   Stub_Endcap_Ring->setAxisTitle("Endcap Ring", 1);
   Stub_Endcap_Ring->setAxisTitle("# L1 Stubs ", 2);
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "NStubs_Disc+" + std::to_string(i + 1);
     // TTStub Endcap stack
     Stub_Endcap_Ring_Fw[i] = iBooker.book1D(HistoName,
                                             HistoName,
                                             psTTStub_ECRing.getParameter<int32_t>("Nbinsx"),
                                             psTTStub_ECRing.getParameter<double>("xmin"),
                                             psTTStub_ECRing.getParameter<double>("xmax"));
     Stub_Endcap_Ring_Fw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_Fw[i]->setAxisTitle("# L1 Stubs ", 2);
   }
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "NStubs_Disc-" + std::to_string(i + 1);
     // TTStub Endcap stack
     Stub_Endcap_Ring_Bw[i] = iBooker.book1D(HistoName,
                                             HistoName,
                                             psTTStub_ECRing.getParameter<int32_t>("Nbinsx"),
                                             psTTStub_ECRing.getParameter<double>("xmin"),
                                             psTTStub_ECRing.getParameter<double>("xmax"));
     Stub_Endcap_Ring_Bw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_Bw[i]->setAxisTitle("# L1 Stubs ", 2);
   }
 
   // TTStub displ/offset
   edm::ParameterSet psTTStub_Barrel_2D = conf_.getParameter<edm::ParameterSet>("TH2TTStub_DisOf_Layer");
   edm::ParameterSet psTTStub_ECDisc_2D = conf_.getParameter<edm::ParameterSet>("TH2TTStub_DisOf_Disc");
   edm::ParameterSet psTTStub_ECRing_2D = conf_.getParameter<edm::ParameterSet>("TH2TTStub_DisOf_Ring");
 
   iBooker.setCurrentFolder(topFolderName_ + "/Stubs/Width");
   HistoName = "Stub_Width_Barrel";
   Stub_Barrel_W = iBooker.book2D(HistoName,
                                  HistoName,
                                  psTTStub_Barrel_2D.getParameter<int32_t>("Nbinsx"),
                                  psTTStub_Barrel_2D.getParameter<double>("xmin"),
                                  psTTStub_Barrel_2D.getParameter<double>("xmax"),
                                  psTTStub_Barrel_2D.getParameter<int32_t>("Nbinsy"),
                                  psTTStub_Barrel_2D.getParameter<double>("ymin"),
                                  psTTStub_Barrel_2D.getParameter<double>("ymax"));
   Stub_Barrel_W->setAxisTitle("Barrel Layer", 1);
   Stub_Barrel_W->setAxisTitle("Displacement - Offset", 2);
 
   HistoName = "Stub_Width_Endcap_Disc";
   Stub_Endcap_Disc_W = iBooker.book2D(HistoName,
                                       HistoName,
                                       psTTStub_ECDisc_2D.getParameter<int32_t>("Nbinsx"),
                                       psTTStub_ECDisc_2D.getParameter<double>("xmin"),
                                       psTTStub_ECDisc_2D.getParameter<double>("xmax"),
                                       psTTStub_ECDisc_2D.getParameter<int32_t>("Nbinsy"),
                                       psTTStub_ECDisc_2D.getParameter<double>("ymin"),
                                       psTTStub_ECDisc_2D.getParameter<double>("ymax"));
   Stub_Endcap_Disc_W->setAxisTitle("Endcap Disc", 1);
   Stub_Endcap_Disc_W->setAxisTitle("Displacement - Offset", 2);
 
   HistoName = "Stub_Width_Endcap_Ring";
   Stub_Endcap_Ring_W = iBooker.book2D(HistoName,
                                       HistoName,
                                       psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                       psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                       psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                       psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                       psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                       psTTStub_ECRing_2D.getParameter<double>("ymax"));
   Stub_Endcap_Ring_W->setAxisTitle("Endcap Ring", 1);
   Stub_Endcap_Ring_W->setAxisTitle("Trigger Offset", 2);
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "Stub_Width_Disc+" + std::to_string(i + 1);
     Stub_Endcap_Ring_W_Fw[i] = iBooker.book2D(HistoName,
                                               HistoName,
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymax"));
     Stub_Endcap_Ring_W_Fw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_W_Fw[i]->setAxisTitle("Displacement - Offset", 2);
   }
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "Stub_Width_Disc-" + std::to_string(i + 1);
     Stub_Endcap_Ring_W_Bw[i] = iBooker.book2D(HistoName,
                                               HistoName,
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymax"));
     Stub_Endcap_Ring_W_Bw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_W_Bw[i]->setAxisTitle("Displacement - Offset", 2);
   }
 
   iBooker.setCurrentFolder(topFolderName_ + "/Stubs/Offset");
   HistoName = "Stub_Offset_Barrel";
   Stub_Barrel_O = iBooker.book2D(HistoName,
                                  HistoName,
                                  psTTStub_Barrel_2D.getParameter<int32_t>("Nbinsx"),
                                  psTTStub_Barrel_2D.getParameter<double>("xmin"),
                                  psTTStub_Barrel_2D.getParameter<double>("xmax"),
                                  psTTStub_Barrel_2D.getParameter<int32_t>("Nbinsy"),
                                  psTTStub_Barrel_2D.getParameter<double>("ymin"),
                                  psTTStub_Barrel_2D.getParameter<double>("ymax"));
   Stub_Barrel_O->setAxisTitle("Barrel Layer", 1);
   Stub_Barrel_O->setAxisTitle("Trigger Offset", 2);
 
   HistoName = "Stub_Offset_Endcap_Disc";
   Stub_Endcap_Disc_O = iBooker.book2D(HistoName,
                                       HistoName,
                                       psTTStub_ECDisc_2D.getParameter<int32_t>("Nbinsx"),
                                       psTTStub_ECDisc_2D.getParameter<double>("xmin"),
                                       psTTStub_ECDisc_2D.getParameter<double>("xmax"),
                                       psTTStub_ECDisc_2D.getParameter<int32_t>("Nbinsy"),
                                       psTTStub_ECDisc_2D.getParameter<double>("ymin"),
                                       psTTStub_ECDisc_2D.getParameter<double>("ymax"));
   Stub_Endcap_Disc_O->setAxisTitle("Endcap Disc", 1);
   Stub_Endcap_Disc_O->setAxisTitle("Trigger Offset", 2);
 
   HistoName = "Stub_Offset_Endcap_Ring";
   Stub_Endcap_Ring_O = iBooker.book2D(HistoName,
                                       HistoName,
                                       psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                       psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                       psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                       psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                       psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                       psTTStub_ECRing_2D.getParameter<double>("ymax"));
   Stub_Endcap_Ring_O->setAxisTitle("Endcap Ring", 1);
   Stub_Endcap_Ring_O->setAxisTitle("Trigger Offset", 2);
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "Stub_Offset_Disc+" + std::to_string(i + 1);
     Stub_Endcap_Ring_O_Fw[i] = iBooker.book2D(HistoName,
                                               HistoName,
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymax"));
     Stub_Endcap_Ring_O_Fw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_O_Fw[i]->setAxisTitle("Trigger Offset", 2);
   }
 
   for (int i = 0; i < numDiscs; i++) {
     HistoName = "Stub_Offset_Disc-" + std::to_string(i + 1);
     Stub_Endcap_Ring_O_Bw[i] = iBooker.book2D(HistoName,
                                               HistoName,
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsx"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmin"),
                                               psTTStub_ECRing_2D.getParameter<double>("xmax"),
                                               psTTStub_ECRing_2D.getParameter<int32_t>("Nbinsy"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymin"),
                                               psTTStub_ECRing_2D.getParameter<double>("ymax"));
     Stub_Endcap_Ring_O_Bw[i]->setAxisTitle("Endcap Ring", 1);
     Stub_Endcap_Ring_O_Bw[i]->setAxisTitle("Trigger Offset", 2);
   }
 }