#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 DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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] = { 0 , 0 , 0 , 0 , 0 }

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

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_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 21 of file OuterTrackerMonitorTTStub.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 39 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 48 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 58 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
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 140 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);
   }
 }