#include <TMTrackProducer.h>

Inheritance diagram for tmtt::TMTrackProducer:

Public Member Functions
	TMTrackProducer (const edm::ParameterSet &, GlobalCacheTMTT const *globalCacheTMTT)

	~TMTrackProducer () override

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< GlobalCacheTMTT > >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	globalEndJob (GlobalCacheTMTT *globalCacheTMTT)

static std::unique_ptr < GlobalCacheTMTT >	initializeGlobalCache (edm::ParameterSet const &iConfig)

Private Types
typedef std::vector< TTTrack < Ref_Phase2TrackerDigi_ > >	TTTrackCollection

Private Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::EDGetTokenT< TTClusterAssMap >	clusterTruthToken_

bool	debug_

std::unique_ptr< DegradeBend >	degradeBend_

std::map< std::string, std::unique_ptr < TrackFitGeneric > >	fitterWorkerMap_

edm::EDGetTokenT < reco::GenJetCollection >	genJetToken_

Histos &	hists_

HTrphi::ErrorMonitor &	htRphiErrMon_

std::list< TrackerModule >	listTrackerModule_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magneticFieldToken_

bool	runRZfilter_

Settings	settings_

const StubAlgorithmOfficial *	stubAlgo_

std::unique_ptr< StubFEWindows >	stubFEWindows_

edm::EDGetTokenT< TTStubDetSetVec >	stubToken_

edm::EDGetTokenT< TTStubAssMap >	stubTruthToken_

StubWindowSuggest &	stubWindowSuggest_

edm::EDGetTokenT < TrackingParticleCollection >	tpToken_

const TrackerGeometry *	trackerGeometry_

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

const TrackerTopology *	trackerTopology_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	trackerTopologyToken_

std::vector< std::string >	trackFitters_

edm::ESGetToken< StubAlgorithm, TTStubAlgorithmRecord >	ttStubAlgoToken_

std::vector< std::string >	useRZfilter_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< GlobalCacheTMTT > >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 42 of file TMTrackProducer.h.

Member Typedef Documentation

typedef std::vector<TTTrack<Ref_Phase2TrackerDigi_> > tmtt::TMTrackProducer::TTTrackCollection

private

Definition at line 52 of file TMTrackProducer.h.

Constructor & Destructor Documentation

tmtt::TMTrackProducer::TMTrackProducer	(	const edm::ParameterSet &	iConfig,
		GlobalCacheTMTT const *	globalCacheTMTT
	)

explicit

Definition at line 38 of file TMTrackProducer.cc.

References tmtt::Settings::clusterTruthInputTag(), clusterTruthToken_, tmtt::trackFitFactory::create(), tmtt::Settings::enableMCtruth(), tmtt::Settings::enableOutputIntermediateTTTracks(), fitterWorkerMap_, tmtt::Settings::genJetInputTag(), genJetToken_, tmtt::Settings::magneticFieldInputTag(), magneticFieldToken_, runRZfilter_, settings_, AlCaHLTBitMon_QueryRunRegistry::string, tmtt::Settings::stubInputTag(), stubToken_, tmtt::Settings::stubTruthInputTag(), stubTruthToken_, tmtt::Settings::tpInputTag(), tpToken_, tmtt::Settings::trackerGeometryInputTag(), trackerGeometryToken_, tmtt::Settings::trackerTopologyInputTag(), trackerTopologyToken_, tmtt::Settings::trackFitters(), trackFitters_, tmtt::Settings::ttStubAlgoInputTag(), ttStubAlgoToken_, tmtt::Settings::useRZfilter(), and useRZfilter_.

       : settings_(iConfig),                                        // Set configuration parameters
         stubWindowSuggest_(globalCacheTMTT->stubWindowSuggest()),  // For tuning FE stub window sizes
         hists_(globalCacheTMTT->hists()),                          // Initialize histograms
         htRphiErrMon_(globalCacheTMTT->htRphiErrMon()),            // rphi HT error monitoring
         debug_(true)                                               // Debug printout
   {
     using namespace edm;
 
     // Get tokens for ES data access.
     magneticFieldToken_ =
         esConsumes<MagneticField, IdealMagneticFieldRecord, Transition::BeginRun>(settings_.magneticFieldInputTag());
     trackerGeometryToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, Transition::BeginRun>(
         settings_.trackerGeometryInputTag());
     trackerTopologyToken_ =
         esConsumes<TrackerTopology, TrackerTopologyRcd, Transition::BeginRun>(settings_.trackerTopologyInputTag());
     ttStubAlgoToken_ =
         esConsumes<StubAlgorithm, TTStubAlgorithmRecord, Transition::BeginRun>(settings_.ttStubAlgoInputTag());
 
     // Get tokens for ED data access.
     stubToken_ = consumes<TTStubDetSetVec>(settings_.stubInputTag());
     if (settings_.enableMCtruth()) {
       // These lines use lots of CPU, even if no use of truth info is made later.
       tpToken_ = consumes<TrackingParticleCollection>(settings_.tpInputTag());
       stubTruthToken_ = consumes<TTStubAssMap>(settings_.stubTruthInputTag());
       clusterTruthToken_ = consumes<TTClusterAssMap>(settings_.clusterTruthInputTag());
       genJetToken_ = consumes<reco::GenJetCollection>(settings_.genJetInputTag());
     }
 
     trackFitters_ = settings_.trackFitters();
     useRZfilter_ = settings_.useRZfilter();
     runRZfilter_ = (not useRZfilter_.empty());  // Do any fitters require an r-z track filter to be run?
 
     // Book histograms.
     //hists_.book();
 
     // Create track fitting algorithm
     for (const string& fitterName : trackFitters_) {
       fitterWorkerMap_[fitterName] = trackFitFactory::create(fitterName, &settings_);
     }
 
     //--- Define EDM output to be written to file (if required)
 
     if (settings_.enableOutputIntermediateTTTracks()) {
       // L1 tracks found by Hough Transform
       produces<TTTrackCollection>("TML1TracksHT").setBranchAlias("TML1TracksHT");
       // L1 tracks found by r-z track filter.
       if (runRZfilter_)
         produces<TTTrackCollection>("TML1TracksRZ").setBranchAlias("TML1TracksRZ");
     }
     // L1 tracks after track fit by each of the fitting algorithms under study
     for (const string& fitterName : trackFitters_) {
       string edmName = string("TML1Tracks") + fitterName;
       produces<TTTrackCollection>(edmName).setBranchAlias(edmName);
     }
   }

tmtt::TMTrackProducer::~TMTrackProducer ( )

inlineoverride

Definition at line 45 of file TMTrackProducer.h.

45 {}

Member Function Documentation

void tmtt::TMTrackProducer::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 97 of file TMTrackProducer.cc.

References tmtt::TrackerModule::ModuleTypeCfg::barrelVsType, tmtt::Settings::barrelVsType(), ecalTB2006H4_GenSimDigiReco_cfg::bField, degradeBend_, edm::ESHandleBase::description(), edm::EventSetup::getData(), edm::EventSetup::getHandle(), edm::getParameterSet(), MagneticField::inTesla(), TrackerTopology::isLower(), listTrackerModule_, magneticFieldToken_, edm::eventsetup::ComponentDescription::pid_, tmtt::TrackerModule::ModuleTypeCfg::pitchVsType, tmtt::Settings::pitchVsType(), tmtt::TrackerModule::ModuleTypeCfg::psVsType, tmtt::Settings::psVsType(), tmtt::StubWindowSuggest::setFEWindows(), tmtt::Settings::setMagneticField(), settings_, tmtt::TrackerModule::ModuleTypeCfg::spaceVsType, tmtt::Settings::spaceVsType(), stubAlgo_, stubFEWindows_, stubWindowSuggest_, DetId::subdetId(), submitPVValidationJobs::t, runonSM::text, StripSubdetector::TID, tmtt::TrackerModule::ModuleTypeCfg::tiltedVsType, tmtt::Settings::tiltedVsType(), StripSubdetector::TOB, trackerGeometry_, trackerGeometryToken_, trackerTopology_, trackerTopologyToken_, ttStubAlgoToken_, and PV3DBase< T, PVType, FrameType >::z().

                                                                                 {
     // Get the B-field and store its value in the Settings class.
     const MagneticField* theMagneticField = &(iSetup.getData(magneticFieldToken_));
     float bField = theMagneticField->inTesla(GlobalPoint(0, 0, 0)).z();  // B field in Tesla.
     settings_.setMagneticField(bField);
 
     // Set also B field in GlobalCacheTMTT (used only for Histogramming)
     globalCache()->settings().setMagneticField(bField);
 
     std::stringstream text;
     text << "\n--- B field = " << bField << " Tesla ---\n";
     std::call_once(
         printOnce, [](string t) { PrintL1trk() << t; }, text.str());
 
     // Get tracker geometry
     trackerGeometry_ = &(iSetup.getData(trackerGeometryToken_));
     trackerTopology_ = &(iSetup.getData(trackerTopologyToken_));
 
     // Loop over tracker modules to get module info.
 
     // Identifies tracker module type for firmware.
     TrackerModule::ModuleTypeCfg moduleTypeCfg;
     moduleTypeCfg.pitchVsType = settings_.pitchVsType();
     moduleTypeCfg.spaceVsType = settings_.spaceVsType();
     moduleTypeCfg.barrelVsType = settings_.barrelVsType();
     moduleTypeCfg.psVsType = settings_.psVsType();
     moduleTypeCfg.tiltedVsType = settings_.tiltedVsType();
 
     listTrackerModule_.clear();
     for (const GeomDet* gd : trackerGeometry_->dets()) {
       DetId detId = gd->geographicalId();
       // Phase 2 Outer Tracker uses TOB for entire barrel & TID for entire endcap.
       if (detId.subdetId() != StripSubdetector::TOB && detId.subdetId() != StripSubdetector::TID)
         continue;
       if (trackerTopology_->isLower(detId)) {  // Select only lower of the two sensors in a module.
         // Store info about this tracker module.
         listTrackerModule_.emplace_back(trackerGeometry_, trackerTopology_, moduleTypeCfg, detId);
       }
     }
 
     // Takes one copy of this to GlobalCacheTMTT for later histogramming.
     globalCache()->setListTrackerModule(listTrackerModule_);
 
     // Get TTStubProducerAlgorithm algorithm, to adjust stub bend FE encoding.
     stubAlgo_ = dynamic_cast<const StubAlgorithmOfficial*>(&iSetup.getData(ttStubAlgoToken_));
     // Get FE stub window size from TTStub producer configuration
     const edm::ESHandle<StubAlgorithm> stubAlgoHandle = iSetup.getHandle(ttStubAlgoToken_);
     const edm::ParameterSet& pSetStubAlgo = getParameterSet(stubAlgoHandle.description()->pid_);
     stubFEWindows_ = std::make_unique<StubFEWindows>(pSetStubAlgo);
     // Initialize utilities needing FE window size.
     stubWindowSuggest_.setFEWindows(stubFEWindows_.get());
     degradeBend_ = std::make_unique<DegradeBend>(trackerTopology_, stubFEWindows_.get(), stubAlgo_);
   }

void tmtt::TMTrackProducer::globalEndJob ( GlobalCacheTMTT * globalCacheTMTT )

static

Definition at line 402 of file TMTrackProducer.cc.

References tmtt::Histos::endJobAnalysis(), tmtt::GlobalCacheTMTT::hists(), tmtt::GlobalCacheTMTT::htRphiErrMon(), tmtt::GlobalCacheTMTT::listTrackerModule(), tmtt::Settings::numEtaRegions(), tmtt::Settings::numPhiSectors(), tmtt::StubWindowSuggest::printResults(), tmtt::Settings::printStubWindows(), tmtt::GlobalCacheTMTT::settings(), tmtt::GlobalCacheTMTT::stubWindowSuggest(), and tmtt::Histos::trackerGeometryAnalysis().

                                                                      {
     const Settings& settings = globalCacheTMTT->settings();
 
     // Print stub window sizes that TMTT recommends CMS uses in FE chips.
     if (settings.printStubWindows())
       globalCacheTMTT->stubWindowSuggest().printResults();
 
     // Print (once) info about tracker geometry.
     globalCacheTMTT->hists().trackerGeometryAnalysis(globalCacheTMTT->listTrackerModule());
 
     PrintL1trk() << "\n Number of (eta,phi) sectors used = (" << settings.numEtaRegions() << ","
                  << settings.numPhiSectors() << ")";
 
     // Print job summary
     globalCacheTMTT->hists().endJobAnalysis(&(globalCacheTMTT->htRphiErrMon()));
   }

std::unique_ptr< GlobalCacheTMTT > tmtt::TMTrackProducer::initializeGlobalCache ( edm::ParameterSet const & iConfig )

static

Definition at line 34 of file TMTrackProducer.cc.

References iConfig.

                                                                                                       {
     return std::make_unique<GlobalCacheTMTT>(iConfig);
   }

void tmtt::TMTrackProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 153 of file TMTrackProducer.cc.

References tmtt::L1fittedTrack::accepted(), clusterTruthToken_, submitPVResolutionJobs::count, debug_, degradeBend_, tmtt::L1fittedTrack::digitizeTrack(), tmtt::Settings::enableDigitize(), tmtt::Settings::enableOutputIntermediateTTTracks(), dataset::end, tmtt::Sector::etaMax(), tmtt::Sector::etaMin(), tmtt::MiniHTstage::exec(), tmtt::MuxHToutputs::exec(), tmtt::Histos::fill(), tmtt::DupFitTrkKiller::filter(), fitterWorkerMap_, genJetToken_, tmtt::InputData::getTPs(), tmtt::Stub::GP, hists_, tmtt::Stub::HT, htRphiErrMon_, listTrackerModule_, tmtt::ConverterToTTTrack::makeTTTrack(), tmtt::Settings::miniHTstage(), eostools::move(), tmtt::Settings::muxOutputsHT(), tmtt::Settings::numEtaRegions(), tmtt::Settings::numPhiSectors(), AlCaHLTBitMon_ParallelJobs::p, tmtt::Sector::phiCentre(), edm::Event::put(), tmtt::Make3Dtracks::run(), runRZfilter_, alignCSCRings::s, settings_, tmtt::HTrphi::store(), AlCaHLTBitMon_QueryRunRegistry::string, stubToken_, stubTruthToken_, stubWindowSuggest_, tmtt::Stub::TF, tpToken_, tmtt::HTbase::trackCands2D(), tmtt::Make3Dtracks::trackCands3D(), trackerGeometry_, trackerTopology_, trackFitters_, and useRZfilter_.

                                                                              {
     // Note useful info about MC truth particles and about reconstructed stubs .
     InputData inputData(iEvent,
                         iSetup,
                         &settings_,
                         &stubWindowSuggest_,
                         degradeBend_.get(),
                         trackerGeometry_,
                         trackerTopology_,
                         listTrackerModule_,
                         tpToken_,
                         stubToken_,
                         stubTruthToken_,
                         clusterTruthToken_,
                         genJetToken_);
 
     const list<TP>& vTPs = inputData.getTPs();
     const list<Stub*>& vStubs = inputData.stubs();
 
     // Creates matrix of Sector objects, which decide which stubs are in which (eta,phi) sector
     Array2D<unique_ptr<Sector>> mSectors(settings_.numPhiSectors(), settings_.numEtaRegions());
     // Create matrix of r-phi Hough-Transform arrays, with one-to-one correspondence to sectors.
     Array2D<unique_ptr<HTrphi>> mHtRphis(settings_.numPhiSectors(), settings_.numEtaRegions());
     // Create matrix of Make3Dtracks objects, to run optional r-z track filter, with one-to-one correspondence to sectors.
     Array2D<unique_ptr<Make3Dtracks>> mMake3Dtrks(settings_.numPhiSectors(), settings_.numEtaRegions());
     // Create matrix of tracks from each fitter in each sector
     Array2D<map<string, std::list<L1fittedTrack>>> mapmFitTrks(settings_.numPhiSectors(), settings_.numEtaRegions());
     // Final tracks after duplicate removal from each track fitter in entire tracker.
     map<string, list<const L1fittedTrack*>> mapFinalTracks;
 
     //=== Initialization
     // Create utility for converting L1 tracks from our private format to official CMSSW EDM format.
     const ConverterToTTTrack converter(&settings_);
 
     // Pointers to TTTrack collections for ED output.
     auto htTTTracksForOutput = std::make_unique<TTTrackCollection>();
     auto rzTTTracksForOutput = std::make_unique<TTTrackCollection>();
     map<string, unique_ptr<TTTrackCollection>> allFitTTTracksForOutput;
     for (const string& fitterName : trackFitters_) {
       auto fitTTTracksForOutput = std::make_unique<TTTrackCollection>();
       allFitTTTracksForOutput[fitterName] = std::move(fitTTTracksForOutput);
     }
 
     //=== Do tracking in the r-phi Hough transform within each sector.
 
     // Fill Hough-Transform arrays with stubs.
     for (unsigned int iPhiSec = 0; iPhiSec < settings_.numPhiSectors(); iPhiSec++) {
       for (unsigned int iEtaReg = 0; iEtaReg < settings_.numEtaRegions(); iEtaReg++) {
         // Initialize constants for this sector.
         mSectors(iPhiSec, iEtaReg) = std::make_unique<Sector>(&settings_, iPhiSec, iEtaReg);
         Sector* sector = mSectors(iPhiSec, iEtaReg).get();
 
         mHtRphis(iPhiSec, iEtaReg) = std::make_unique<HTrphi>(
             &settings_, iPhiSec, iEtaReg, sector->etaMin(), sector->etaMax(), sector->phiCentre(), &htRphiErrMon_);
         HTrphi* htRphi = mHtRphis(iPhiSec, iEtaReg).get();
 
         // Check sector is enabled (always true, except if user disabled some for special studies).
         if (settings_.isHTRPhiEtaRegWhitelisted(iEtaReg)) {
           for (Stub* stub : vStubs) {
             // Digitize stub as would be at input to GP. This doesn't need the nonant number, since we assumed an integer number of
             // phi digitisation  bins inside an nonant. N.B. This changes the coordinates & bend stored in the stub.
 
             if (settings_.enableDigitize())
               stub->digitize(iPhiSec, Stub::DigiStage::GP);
 
             // Check if stub is inside this sector
             bool inside = sector->inside(stub);
 
             if (inside) {
               // Check which eta subsectors within the sector the stub is compatible with (if subsectors being used).
               const vector<bool> inEtaSubSecs = sector->insideEtaSubSecs(stub);
 
               // Digitize stub if as would be at input to HT, which slightly degrades its coord. & bend resolution, affecting the HT performance.
               if (settings_.enableDigitize())
                 stub->digitize(iPhiSec, Stub::DigiStage::HT);
 
               // Store stub in Hough transform array for this sector, indicating its compatibility with eta subsectors with sector.
               htRphi->store(stub, inEtaSubSecs);
             }
           }
         }
 
         // Find tracks in r-phi HT array.
         htRphi->end();  // Calls htArrayRphi_.end() -> HTBase::end()
       }
     }
 
     if (settings_.muxOutputsHT() > 0) {
       // Multiplex outputs of several HT onto one pair of output opto-links.
       // This only affects tracking performance if option busySectorKill is enabled, so that tracks that
       // can't be sent down the link within the time-multiplexed period are killed.
       MuxHToutputs muxHT(&settings_);
       muxHT.exec(mHtRphis);
     }
 
     // Optionally, run 2nd stage mini HT -- WITHOUT TRUNCATION ???
     if (settings_.miniHTstage()) {
       MiniHTstage miniHTstage(&settings_);
       miniHTstage.exec(mHtRphis);
     }
 
     //=== Make 3D tracks, optionally running r-z track filters (such as Seed Filter) & duplicate track removal.
 
     for (unsigned int iPhiSec = 0; iPhiSec < settings_.numPhiSectors(); iPhiSec++) {
       for (unsigned int iEtaReg = 0; iEtaReg < settings_.numEtaRegions(); iEtaReg++) {
         const Sector* sector = mSectors(iPhiSec, iEtaReg).get();
 
         // Get tracks found by r-phi HT.
         const HTrphi* htRphi = mHtRphis(iPhiSec, iEtaReg).get();
         const list<L1track2D>& vecTracksRphi = htRphi->trackCands2D();
 
         // Initialize utility for making 3D tracks from 2D ones.
         mMake3Dtrks(iPhiSec, iEtaReg) = std::make_unique<Make3Dtracks>(
             &settings_, iPhiSec, iEtaReg, sector->etaMin(), sector->etaMax(), sector->phiCentre());
         Make3Dtracks* make3Dtrk = mMake3Dtrks(iPhiSec, iEtaReg).get();
 
         // Convert 2D tracks found by HT to 3D tracks (optionally by running r-z filters & duplicate track removal)
         make3Dtrk->run(vecTracksRphi);
 
         if (settings_.enableOutputIntermediateTTTracks()) {
           // Convert these tracks to EDM format for output (used for collaborative work outside TMTT group).
           // Do this for tracks output by HT & optionally also for those output by r-z track filter.
           const list<L1track3D>& vecTrk3D_ht = make3Dtrk->trackCands3D(false);
           for (const L1track3D& trk : vecTrk3D_ht) {
             TTTrack<Ref_Phase2TrackerDigi_> htTTTrack = converter.makeTTTrack(&trk, iPhiSec, iEtaReg);
             htTTTracksForOutput->push_back(htTTTrack);
           }
 
           if (runRZfilter_) {
             const list<L1track3D>& vecTrk3D_rz = make3Dtrk->trackCands3D(true);
             for (const L1track3D& trk : vecTrk3D_rz) {
               TTTrack<Ref_Phase2TrackerDigi_> rzTTTrack = converter.makeTTTrack(&trk, iPhiSec, iEtaReg);
               rzTTTracksForOutput->push_back(rzTTTrack);
             }
           }
         }
       }
     }
 
     //=== Do a helix fit to all the track candidates.
 
     // Loop over all the fitting algorithms we are trying.
     for (const string& fitterName : trackFitters_) {
       for (unsigned int iPhiSec = 0; iPhiSec < settings_.numPhiSectors(); iPhiSec++) {
         for (unsigned int iEtaReg = 0; iEtaReg < settings_.numEtaRegions(); iEtaReg++) {
           const Make3Dtracks* make3Dtrk = mMake3Dtrks(iPhiSec, iEtaReg).get();
 
           // Does this fitter require r-z track filter to be run before it?
           bool useRZfilt = (std::count(useRZfilter_.begin(), useRZfilter_.end(), fitterName) > 0);
 
           // Get 3D track candidates found by Hough transform (plus optional r-z filters/duplicate removal) in this sector.
           const list<L1track3D>& vecTrk3D = make3Dtrk->trackCands3D(useRZfilt);
 
           // Find list where fitted tracks will be stored.
           list<L1fittedTrack>& fitTrksInSec = mapmFitTrks(iPhiSec, iEtaReg)[fitterName];
 
           // Fit all tracks in this sector
           for (const L1track3D& trk : vecTrk3D) {
             // Ensure stubs assigned to this track is digitized with respect to the phi sector the track is in.
             if (settings_.enableDigitize()) {
               const vector<Stub*>& stubsOnTrk = trk.stubs();
               for (Stub* s : stubsOnTrk) {
                 // Also digitize stub in way this specific track fitter uses it.
                 s->digitize(iPhiSec, Stub::DigiStage::TF);
               }
             }
 
             L1fittedTrack fitTrk = fitterWorkerMap_[fitterName]->fit(trk);
 
             if (fitTrk.accepted()) {  // If fitter accepted track, then store it.
               // Optionally digitize fitted track, degrading slightly resolution.
               if (settings_.enableDigitize())
                 fitTrk.digitizeTrack(fitterName);
               // Store fitted tracks, such that there is one fittedTracks corresponding to each HT tracks.
               fitTrksInSec.push_back(fitTrk);
             }
           }
         }
       }
     }
 
     // Run duplicate track removal on the fitted tracks if requested.
 
     // Initialize the duplicate track removal algorithm that can optionally be run after the track fit.
     DupFitTrkKiller killDupFitTrks(&settings_);
 
     // Loop over all the fitting algorithms we used.
     for (const string& fitterName : trackFitters_) {
       for (unsigned int iPhiSec = 0; iPhiSec < settings_.numPhiSectors(); iPhiSec++) {
         for (unsigned int iEtaReg = 0; iEtaReg < settings_.numEtaRegions(); iEtaReg++) {
           // Get fitted tracks in sector
           const list<L1fittedTrack>& fitTrksInSec = mapmFitTrks(iPhiSec, iEtaReg)[fitterName];
 
           // Run duplicate removal
           list<const L1fittedTrack*> filteredFitTrksInSec = killDupFitTrks.filter(fitTrksInSec);
 
           // Prepare TTTrack collection.
           for (const L1fittedTrack* fitTrk : filteredFitTrksInSec) {
             // Convert these fitted tracks to EDM format for output (used for collaborative work outside TMTT group).
             TTTrack<Ref_Phase2TrackerDigi_> fitTTTrack = converter.makeTTTrack(fitTrk, iPhiSec, iEtaReg);
             allFitTTTracksForOutput[fitterName]->push_back(fitTTTrack);
           }
 
           // Store fitted tracks from entire tracker.
           mapFinalTracks[fitterName].insert(
               mapFinalTracks[fitterName].end(), filteredFitTrksInSec.begin(), filteredFitTrksInSec.end());
         }
       }
     }
 
     // Debug printout
     if (debug_) {
       PrintL1trk() << "INPUT #TPs = " << vTPs.size() << " #STUBs = " << vStubs.size();
       unsigned int numHTtracks = 0;
       for (unsigned int iPhiSec = 0; iPhiSec < settings_.numPhiSectors(); iPhiSec++) {
         for (unsigned int iEtaReg = 0; iEtaReg < settings_.numEtaRegions(); iEtaReg++) {
           const Make3Dtracks* make3Dtrk = mMake3Dtrks(iPhiSec, iEtaReg).get();
           numHTtracks += make3Dtrk->trackCands3D(false).size();
         }
       }
       PrintL1trk() << "Number of tracks after HT = " << numHTtracks;
       for (const auto& p : mapFinalTracks) {
         const string& fitName = p.first;
         const list<const L1fittedTrack*> fittedTracks = p.second;
         PrintL1trk() << "Number of tracks after " << fitName << " track helix fit = " << fittedTracks.size();
       }
     }
 
     // Allow histogramming to plot undigitized variables.
     for (Stub* stub : vStubs) {
       if (settings_.enableDigitize())
         stub->setDigitizeWarningsOn(false);
     }
 
     // Fill histograms to monitor input data & tracking performance.
     hists_.fill(inputData, mSectors, mHtRphis, mMake3Dtrks, mapFinalTracks);
 
     //=== Store output EDM track and hardware stub collections.
     if (settings_.enableOutputIntermediateTTTracks()) {
       iEvent.put(std::move(htTTTracksForOutput), "TML1TracksHT");
       if (runRZfilter_)
         iEvent.put(std::move(rzTTTracksForOutput), "TML1TracksRZ");
     }
     for (const string& fitterName : trackFitters_) {
       string edmName = string("TML1Tracks") + fitterName;
       iEvent.put(std::move(allFitTTTracksForOutput[fitterName]), edmName);
     }
   }