#include <SectorProcessor.h>

Public Types
typedef unsigned long long	EventNumber_t

typedef PatternRecognition::pattern_ref_t	pattern_ref_t

Public Member Functions
void	configure (const GeometryTranslator tp_geom, const ConditionHelper cond, const SectorProcessorLUT lut, PtAssignmentEngine pt_assign_engine, int verbose, int endcap, int sector, int minBX, int maxBX, int bxWindow, int bxShiftCSC, int bxShiftRPC, int bxShiftGEM, std::string era, const std::vector< int > &zoneBoundaries, int zoneOverlap, bool includeNeighbor, bool duplicateTheta, bool fixZonePhi, bool useNewZones, bool fixME11Edges, const std::vector< std::string > &pattDefinitions, const std::vector< std::string > &symPattDefinitions, bool useSymPatterns, int thetaWindow, int thetaWindowZone0, bool useRPC, bool useSingleHits, bool bugSt2PhDiff, bool bugME11Dupes, bool bugAmbigThetaWin, bool twoStationSameBX, int maxRoadsPerZone, int maxTracks, bool useSecondEarliest, bool bugSameSectorPt0, bool readPtLUTFile, bool fixMode15HighPt, bool bug9BitDPhi, bool bugMode7CLCT, bool bugNegPt, bool bugGMTPhi, bool promoteMode7, int modeQualVer)

void	configure_by_fw_version (unsigned fw_version)

void	process (EventNumber_t ievent, const TriggerPrimitiveCollection &muon_primitives, EMTFHitCollection &out_hits, EMTFTrackCollection &out_tracks) const

void	process_single_bx (int bx, const TriggerPrimitiveCollection &muon_primitives, EMTFHitCollection &out_hits, EMTFTrackCollection &out_tracks, std::deque< EMTFHitCollection > &extended_conv_hits, std::deque< EMTFTrackCollection > &extended_best_track_cands, std::map< pattern_ref_t, int > &patt_lifetime_map) const

	SectorProcessor ()

	~SectorProcessor ()

Private Attributes
bool	bug9BitDPhi_

bool	bugAmbigThetaWin_

bool	bugGMTPhi_

bool	bugME11Dupes_

bool	bugMode7CLCT_

bool	bugNegPt_

bool	bugSameSectorPt0_

bool	bugSt2PhDiff_

int	bxShiftCSC_

int	bxShiftGEM_

int	bxShiftRPC_

int	bxWindow_

const ConditionHelper *	cond_

bool	duplicateTheta_

int	endcap_

std::string	era_

bool	fixME11Edges_

bool	fixMode15HighPt_

bool	fixZonePhi_

bool	includeNeighbor_

const SectorProcessorLUT *	lut_

int	maxBX_

int	maxRoadsPerZone_

int	maxTracks_

int	minBX_

int	modeQualVer_

std::vector< std::string >	pattDefinitions_

bool	promoteMode7_

PtAssignmentEngine *	pt_assign_engine_

bool	readPtLUTFile_

int	sector_

std::vector< std::string >	symPattDefinitions_

int	thetaWindow_

int	thetaWindowZone0_

const GeometryTranslator *	tp_geom_

bool	twoStationSameBX_

bool	useCPPF_

bool	useNewZones_

bool	useRPC_

bool	useSecondEarliest_

bool	useSingleHits_

bool	useSymPatterns_

int	verbose_

std::vector< int >	zoneBoundaries_

int	zoneOverlap_

Detailed Description

Definition at line 29 of file SectorProcessor.h.

Member Typedef Documentation

typedef unsigned long long SectorProcessor::EventNumber_t

Definition at line 34 of file SectorProcessor.h.

typedef PatternRecognition::pattern_ref_t SectorProcessor::pattern_ref_t

Definition at line 35 of file SectorProcessor.h.

Constructor & Destructor Documentation

SectorProcessor::SectorProcessor ( )

explicit

Definition at line 4 of file SectorProcessor.cc.

                                  {
 
 }

SectorProcessor::~SectorProcessor ( )

Definition at line 8 of file SectorProcessor.cc.

                                   {
 
 }

Member Function Documentation

void SectorProcessor::configure	(	const GeometryTranslator *	tp_geom,
		const ConditionHelper *	cond,
		const SectorProcessorLUT *	lut,
		PtAssignmentEngine *	pt_assign_engine,
		int	verbose,
		int	endcap,
		int	sector,
		int	minBX,
		int	maxBX,
		int	bxWindow,
		int	bxShiftCSC,
		int	bxShiftRPC,
		int	bxShiftGEM,
		std::string	era,
		const std::vector< int > &	zoneBoundaries,
		int	zoneOverlap,
		bool	includeNeighbor,
		bool	duplicateTheta,
		bool	fixZonePhi,
		bool	useNewZones,
		bool	fixME11Edges,
		const std::vector< std::string > &	pattDefinitions,
		const std::vector< std::string > &	symPattDefinitions,
		bool	useSymPatterns,
		int	thetaWindow,
		int	thetaWindowZone0,
		bool	useRPC,
		bool	useSingleHits,
		bool	bugSt2PhDiff,
		bool	bugME11Dupes,
		bool	bugAmbigThetaWin,
		bool	twoStationSameBX,
		int	maxRoadsPerZone,
		int	maxTracks,
		bool	useSecondEarliest,
		bool	bugSameSectorPt0,
		bool	readPtLUTFile,
		bool	fixMode15HighPt,
		bool	bug9BitDPhi,
		bool	bugMode7CLCT,
		bool	bugNegPt,
		bool	bugGMTPhi,
		bool	promoteMode7,
		int	modeQualVer
	)

Definition at line 12 of file SectorProcessor.cc.

References bug9BitDPhi_, bugAmbigThetaWin_, bugGMTPhi_, bugME11Dupes_, bugMode7CLCT_, bugNegPt_, bugSameSectorPt0_, bugSt2PhDiff_, bxShiftCSC_, bxShiftGEM_, bxShiftRPC_, bxWindow_, cond_, duplicateTheta_, makeMuonMisalignmentScenario::endcap, endcap_, year_2016_cff::era, era_, fixME11Edges_, fixMode15HighPt_, fixZonePhi_, includeNeighbor_, lumiPlot::lut, lut_, emtf::MAX_ENDCAP, emtf::MAX_TRIGSECTOR, simTrackMatching_cfi::maxBX, maxBX_, maxRoadsPerZone_, RecoTauCombinatoricProducer_cfi::maxTracks, maxTracks_, emtf::MIN_ENDCAP, emtf::MIN_TRIGSECTOR, simTrackMatching_cfi::minBX, minBX_, modeQualVer_, pattDefinitions_, promoteMode7_, pt_assign_engine_, readPtLUTFile_, sector_, symPattDefinitions_, thetaWindow_, thetaWindowZone0_, tp_geom_, twoStationSameBX_, useNewZones_, useRPC_, useSecondEarliest_, useSingleHits_, useSymPatterns_, verbose, verbose_, zoneBoundaries_, and zoneOverlap_.

   {
   if (not(emtf::MIN_ENDCAP <= endcap && endcap <= emtf::MAX_ENDCAP))
     { edm::LogError("L1T") << "emtf::MIN_ENDCAP = " << emtf::MIN_ENDCAP 
                << ", emtf::MAX_ENDCAP = " << emtf::MAX_ENDCAP
                << ", endcap = " << endcap; return; }
   if (not(emtf::MIN_TRIGSECTOR <= sector && sector <= emtf::MAX_TRIGSECTOR))
     { edm::LogError("L1T") << "emtf::MIN_TRIGSECTOR = " << emtf::MIN_TRIGSECTOR 
                << ", emtf::MAX_TRIGSECTOR = " << emtf::MAX_TRIGSECTOR
                << ", endcap = " << sector; return; }
   if (not(tp_geom != nullptr))
     { edm::LogError("L1T") << "tp_geom = nullptr"; return; }
   if (not(cond != nullptr))
     { edm::LogError("L1T") << "cond = nullptr"; return; }
   if (not(lut != nullptr))
     { edm::LogError("L1T") << "lut = nullptr"; return; }
   if (not(pt_assign_engine != nullptr))
     { edm::LogError("L1T") << "pt_assign_engine = nullptr"; return; }
 
   tp_geom_          = tp_geom;
   cond_             = cond;
   lut_              = lut;
   pt_assign_engine_ = pt_assign_engine;
 
   verbose_  = verbose;
   endcap_   = endcap;
   sector_   = sector;
 
   minBX_       = minBX;
   maxBX_       = maxBX;
   bxWindow_    = bxWindow;
   bxShiftCSC_  = bxShiftCSC;
   bxShiftRPC_  = bxShiftRPC;
   bxShiftGEM_  = bxShiftGEM;
 
   era_ = era;
 
   zoneBoundaries_     = zoneBoundaries;
   zoneOverlap_        = zoneOverlap;
   includeNeighbor_    = includeNeighbor;
   duplicateTheta_     = duplicateTheta;
   fixZonePhi_         = fixZonePhi;
   useNewZones_        = useNewZones;
   fixME11Edges_       = fixME11Edges;
 
   pattDefinitions_    = pattDefinitions;
   symPattDefinitions_ = symPattDefinitions;
   useSymPatterns_     = useSymPatterns;
 
   thetaWindow_        = thetaWindow;
   thetaWindowZone0_   = thetaWindowZone0;
   useRPC_             = useRPC;
   useSingleHits_      = useSingleHits;
   bugSt2PhDiff_       = bugSt2PhDiff;
   bugME11Dupes_       = bugME11Dupes;
   bugAmbigThetaWin_   = bugAmbigThetaWin;
   twoStationSameBX_   = twoStationSameBX;
 
   maxRoadsPerZone_    = maxRoadsPerZone;
   maxTracks_          = maxTracks;
   useSecondEarliest_  = useSecondEarliest;
   bugSameSectorPt0_   = bugSameSectorPt0;
 
   readPtLUTFile_      = readPtLUTFile;
   fixMode15HighPt_    = fixMode15HighPt;
   bug9BitDPhi_        = bug9BitDPhi;
   bugMode7CLCT_       = bugMode7CLCT;
   bugNegPt_           = bugNegPt;
   bugGMTPhi_          = bugGMTPhi;
   promoteMode7_       = promoteMode7;
   modeQualVer_        = modeQualVer;
 }

void SectorProcessor::configure_by_fw_version ( unsigned fw_version )

Settings for 2018 (by default just use settings in simEmtfDigis_cfi.py) ///

Settings for 2017 (by default just use settings in simEmtfDigis_cfi.py) ///

Settings for all of 2016 (following order in simEmtfDigis_cfi.py) ///

Definition at line 99 of file SectorProcessor.cc.

References bug9BitDPhi_, bugAmbigThetaWin_, bugGMTPhi_, bugME11Dupes_, bugMode7CLCT_, bugNegPt_, bugSameSectorPt0_, bugSt2PhDiff_, bxWindow_, gather_cfg::cout, era_, fixME11Edges_, fixMode15HighPt_, fixZonePhi_, maxBX_, modeQualVer_, promoteMode7_, thetaWindow_, thetaWindowZone0_, twoStationSameBX_, useRPC_, useSecondEarliest_, useSymPatterns_, and verbose_.

                                                                  {
   if (verbose_ > 0) {
     std::cout << "Configure SectorProcessor with fw_version: " << fw_version << std::endl;
   }
   
   if (fw_version == 0 || fw_version == 123456)  // fw_version '123456' is from the fake conditions
     return;
 
   // For now, no switches later than FW version 47864 (end-of-year 2016)
   // Beggining in late 2016, "fw_version" in O2O populated with timestamp, rather than FW version
   // tm fw_time = gmtime(fw_version);  (See https://linux.die.net/man/3/gmtime, https://www.epochconverter.com)
 
   if (fw_version >= 1514764800) {  // January 1, 2018
 
     // Settings for all of 2018 (following order in simEmtfDigis_cfi.py)
     // BXWindow(2) and BugAmbigThetaWin(F) deployed sometime before stable beams on March 20, not quite sure when - AWB 26.04.18
     // TwoStationSameBX(T), ThetaWindowZone0(4), and ModeQualVer(2) to be deployed sometime between May 17 and May 31 - AWB 14.05.18
 
     // Global parameters
     // Defaults : CSCEnable(T), RPCEnable(T), GEMEnable(F), Era("Run2_2018"), MinBX(-3), MaxBX(+3), BXWindow(2)
     // --------------------------------------------------------------------------------------------------------
     era_      = "Run2_2018";  // Era for CMSSW customization
     bxWindow_ = 2;            // Number of BX whose primitives can be included in the same track
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : PrimConvLUT(1), ZoneBoundaries(0,41,49,87,127), ZoneOverlap(2), IncludeNeighbor(T),
     //            DuplicateThteta(T), FixZonePhi(T), UseNewZones(F), FixME11Edges(T)
     // ------------------------------------------------------------------------------
 
     // spPRParams16 : Sector processor pattern-recognition parameters
     // Defaults : PatternDefinitions(...), SymPatternDefinitions(...), UseSymmetricalPatterns(T)
     // -----------------------------------------------------------------------------------------
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : ThetaWindow(8), ThetaWindowZone0(4), UseSingleHits(F), BugSt2PhDiff(F),
     //            BugME11Dupes(F), BugAmbigThetaWin(F), TwoStationSameBX(T)
     // ----------------------------------------------------------------------------------
     thetaWindow_      = 8;     // Maximum dTheta between primitives in the same track
     thetaWindowZone0_ = 4;     // Maximum dTheta between primitives in the same track in Zone 0 (ring 1)
     bugAmbigThetaWin_ = false; // Can allow dThetas outside window when there are 2 LCTs in the same chamber
     twoStationSameBX_ = true;  // Requires the hits in two-station tracks to have the same BX
 
     // spGCParams16 : Sector processor ghost-cancellation parameters
     // Defaults : MaxRoadsPerZone(3), MaxTracks(3), UseSecondEarliest(T), BugSameSectorPt0(F)
     // --------------------------------------------------------------------------------------
 
     // spPAParams16 : Sector processor pt-assignment parameters
     // Defaults : ReadPtLUTFile(F), FixMode15HighPt(T), Bug9BitDPhi(F), BugMode7CLCT(F),
     //            BugNegPt(F), BugGMTPhi(F), PromoteMode7(F), ModeQualVer(2)
     // ---------------------------------------------------------------------------------
     modeQualVer_ = 2;      // Version 2 contains modified mode-quality mapping for 2018
     promoteMode7_ = false; // Assign station 2-3-4 tracks with |eta| > 1.6 SingleMu quality
 
 
     // ___________________________________________________________________________
     // Versions in 2018 - no external documentation
     // As of the beginning of 2018 EMTF O2O was broken, not updating the database with online conditions
     // Firmware version reported for runs before 318841 (June 28) is 1504018578 (Aug. 29, 2017) even though
     //   updates occured in February and March of 2018.  May need to re-write history in the database. - AWB 30.06.18
 
     if (fw_version < 1529520380) {  // June 20, 2018
       thetaWindowZone0_ = 8;     // Maximum dTheta between primitives in the same track in Zone 0 (ring 1)
       twoStationSameBX_ = false; // Requires the hits in two-station tracks to have the same BX
       modeQualVer_ = 1;          // Version 2 contains modified mode-quality mapping for 2018
       promoteMode7_ = true;      // Assign station 2-3-4 tracks with |eta| > 1.6 SingleMu quality
     }
 
     return;
   }
 
   else if (fw_version >= 50000) {
 
     // Settings for all of 2017 (following order in simEmtfDigis_cfi.py)
 
     // Global parameters
     // Defaults : CSCEnable(T), RPCEnable(T), GEMEnable(F), Era("Run2_2018"), MinBX(-3), MaxBX(+3), BXWindow(2)
     // --------------------------------------------------------------------------------------------------------
     era_      = "Run2_2017";  // Era for CMSSW customization
     bxWindow_ = 3;            // Number of BX whose primitives can be included in the same track
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : PrimConvLUT(1), ZoneBoundaries(0,41,49,87,127), ZoneOverlap(2), IncludeNeighbor(T),
     //            DuplicateThteta(T), FixZonePhi(T), UseNewZones(F), FixME11Edges(T)
     // ------------------------------------------------------------------------------
 
     // spPRParams16 : Sector processor pattern-recognition parameters
     // Defaults : PatternDefinitions(...), SymPatternDefinitions(...), UseSymmetricalPatterns(T)
     // -----------------------------------------------------------------------------------------
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : ThetaWindow(8), ThetaWindowZone0(4), UseSingleHits(F), BugSt2PhDiff(F),
     //            BugME11Dupes(F), BugAmbigThetaWin(F), TwoStationSameBX(T)
     // ----------------------------------------------------------------------------------
     thetaWindow_      = 8;     // Maximum dTheta between primitives in the same track
     thetaWindowZone0_ = 8;     // Maximum dTheta between primitives in the same track in Zone 0 (ring 1)
     bugAmbigThetaWin_ = true;  // Can allow dThetas outside window when there are 2 LCTs in the same chamber
     twoStationSameBX_ = false; // Requires the hits in two-station tracks to have the same BX
 
     // spGCParams16 : Sector processor ghost-cancellation parameters
     // Defaults : MaxRoadsPerZone(3), MaxTracks(3), UseSecondEarliest(T), BugSameSectorPt0(F)
     // --------------------------------------------------------------------------------------
 
     // spPAParams16 : Sector processor pt-assignment parameters
     // Defaults : ReadPtLUTFile(F), FixMode15HighPt(T), Bug9BitDPhi(F), BugMode7CLCT(F),
     //            BugNegPt(F), BugGMTPhi(F), PromoteMode7(F)
     // ---------------------------------------------------------------------------------
     modeQualVer_ = 1; // Version 2 contains modified mode-quality mapping for 2018
 
     // ___________________________________________________________________________
     // Versions in 2017 - no full documentation, can refer to https://twiki.cern.ch/twiki/bin/viewauth/CMS/L1KnownIssues
 
     // Before July 9th (runs < 298653), all mode 7 tracks (station 2-3-4) assigned quality 11
     // July 9th - 29th (runs 298653 - 300087), mode 7 tracks with |eta| > 1.6 in sector -6 assigned quality 12
     // After July 29th (runs >= 300088), mode 7 track promotion applied in all sectors
     // For some reason, the FW version in the database is 1496792995, at least for runs >= 298034 (July 4),
     //   which is towards the end of run 2017B (could not check earlier runs).  This corresponds to the date "June 6th",
     //   which is a month earlier than the first firmware update to apply this promotion.  So something's screwey.
     // Since July 18 is in the middle of the single-sector-fix period, would like to use a firmware version with
     //   roughly that date.  But this may require an intervention in the database. - AWB 04.08.17
     // Last firmware version in 2017 was 1504018578 (Aug. 29, 2017).
     if (fw_version >= 1496792995)
       promoteMode7_ = true;  // Assign station 2-3-4 tracks with |eta| > 1.6 SingleMu quality
 
     return;
   }
 
   else {
 
     // Global parameters
     // Defaults : CSCEnable(T), RPCEnable(T), GEMEnable(F), Era("Run2_2018"), MinBX(-3), MaxBX(+3), BXWindow(2)
     // --------------------------------------------------------------------------------------------------------
     useRPC_   = false;        // Use clustered RPC hits from CPPF in track-building
     era_      = "Run2_2016";  // Era for CMSSW customization
     // maxBX_                 // Depends on FW version, see below
     bxWindow_ = 3;            // Number of BX whose primitives can be included in the same track
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : PrimConvLUT(1), ZoneBoundaries(0,41,49,87,127), ZoneOverlap(2), IncludeNeighbor(T),
     //            DuplicateThteta(T), FixZonePhi(T), UseNewZones(F), FixME11Edges(T)
     // ------------------------------------------------------------------------------
     // primConvLUT_         // Should be 0 for 2016, set using get_pc_lut_version() from ConditionsHelper.cc
     // fixZonePhi_          // Depends on FW version, see below
     fixME11Edges_ = false;  // Improved small fraction of buggy LCT coordinate transformations
 
     // spPRParams16 : Sector processor pattern-recognition parameters
     // Defaults : PatternDefinitions(...), SymPatternDefinitions(...), UseSymmetricalPatterns(T)
     // -----------------------------------------------------------------------------------------
     // useSymPatterns_  // Depends on FW version, see below
 
     // spTBParams16 : Sector processor track-building parameters
     // Defaults : ThetaWindow(8), ThetaWindowZone0(4), UseSingleHits(F), BugSt2PhDiff(F),
     //            BugME11Dupes(F), BugAmbigThetaWin(F), TwoStationSameBX(T)
     // ----------------------------------------------------------------------------------
     thetaWindow_      = 4;     // Maximum dTheta between primitives in the same track
     thetaWindowZone0_ = 4;     // Maximum dTheta between primitives in the same track in Zone 0 (ring 1)
     // bugSt2PhDiff_           // Depends on FW version, see below
     // bugME11Dupes_           // Depends on FW version, see below
     bugAmbigThetaWin_ = true;  // Can allow dThetas outside window when there are 2 LCTs in the same chamber
     twoStationSameBX_ = false; // Requires the hits in two-station tracks to have the same BX
 
     // spGCParams16 : Sector processor ghost-cancellation parameters
     // Defaults : MaxRoadsPerZone(3), MaxTracks(3), UseSecondEarliest(T), BugSameSectorPt0(F)
     // --------------------------------------------------------------------------------------
     // useSecondEarliest_  // Depends on FW version, see below
     // bugSameSectorPt0_   // Depends on FW version, see below
 
     // spPAParams16 : Sector processor pt-assignment parameters
     // Defaults : ReadPtLUTFile(F), FixMode15HighPt(T), Bug9BitDPhi(F), BugMode7CLCT(F),
     //            BugNegPt(F), BugGMTPhi(F), PromoteMode7(F)
     // ---------------------------------------------------------------------------------
     // fixMode15HighPt_   // Depends on FW version, see below
     // bug9BitDPhi_       // Depends on FW version, see below
     // bugMode7CLCT_      // Depends on FW version, see below
     // bugNegPt_          // Depends on FW version, see below
     bugGMTPhi_   = true;  // Some drift in uGMT phi conversion, off by up to a few degrees
     modeQualVer_ = 1;     // Version 2 contains modified mode-quality mapping for 2018
 
   } // End default settings for 2016
 
 
   // ___________________________________________________________________________
   // Versions in 2016 - refer to docs/EMTF_FW_LUT_versions_2016_draft2.xlsx
 
   // 1st_LCT_BX / 2nd_LCT_BX  (should also make unpacker configurable - AWB 21.07.17)
   // FW: Before: Earliest LCT used to assign BX, tracks only cancel within same BX
   //     After:  Second-earliest LCT used to assign BX, tracks cancel over 3 BX, improved LCT recovery
   useSecondEarliest_  = (fw_version < 46773) ? false : true;  // Changed Sept. 5
 
   // 8_BX_readout / 7_BX_readout
   // SW: DAQ readout changed from to [-3, +4] BX to [-3, +3] BX
   maxBX_              = (fw_version < 47109) ? +4 : +3;       // Changed Sept. 28
 
   // Asymm_patterns / Symm_patterns
   // FW: Changed from 9 asymmetric patterns to 5 symmetric patterns for track building
   useSymPatterns_     = (fw_version < 47214) ? false : true;  // Changed Oct. 6
 
   // HiPt_outlier
   // LUT: High-pT fix puts outlier LCTs in mode 15 tracks back in a straight line
   fixMode15HighPt_    = (fw_version < 46650) ? false : true;  // Changed July 25
 
   // Link_monitor (unpacker only)
   // FW: Added MPC link monitoring
 
   // ___________________________________________________________________________
   // Bugs
 
   // DAQ_ID (unpacker only; should make configurable - AWB 21.07.17)
   // FW: DAQ ME with output CSC ID range 0 - 8 instead of 1 - 9
   //     SP output ME2_ID, 3_ID, and 4_ID filled with 4, 5, or 6 when they should have been 7, 8, or 9.
 
   // ME_ID_FR
   // FW: Incorrect ME_ID fields in DAQ, wrong FR bits and some dPhi wrap-around in pT LUT address
   // - Unpacker only, or not worth emulating
 
   // DAQ_miss_LCT (unpacker only)
   // FW: LCTs only output if there was a track in the sector
 
   // Sector_pT_0
   // FW: Only highest-quality track in a sector assigned pT; others assigned pT = 0
   bugSameSectorPt0_   = (fw_version < 46650) ? true : false;  // Fixed July 22
 
   // Sector_bad_pT
   // FW: Tracks sometimes assigned pT of track in previous BX
   // - This is an ongoing (very rare) bug which occurs when 2 tracks try to access the same "bank" in the pT LUT
   //   It would be very difficult to emulate exactly, but the logic from Alex Madorsky is below
   // ## macro for detecting same bank address
   // ## bank and chip must match, and valid flags must be set
   // ## a and b are indexes 0,1,2
   // ## [X:Y] are bit portions from ptlut address words
   // `define sb(a,b) (ptlut_addr[a][29:26] == ptlut_addr[b][29:26] && ptlut_addr[a][5:2] == ptlut_addr[b][5:2] && ptlut_addr_val[a] && ptlut_addr_val[b])
   // ## This macro is used like this:
   // if (`sb(0,2) || `sb(1,2)) {disable PT readout for track 2}
 
 
   // DAQ_BX_3_LCT (unpacker only)
   // SW: LCTs in BX -3 only reported if there was a track in the sector
   // - not applicable
 
   // DAQ_BX_23_LCT (unpacker only)
   // SW: LCTs in BX -2 and -3 only reported if there was a track in the sector
   // - not applicable
 
   // pT_dPhi_bits
   // FW: dPhi wrap-around in modes 3, 5, 6, 9, 10, 12
   bug9BitDPhi_        = (fw_version < 47214) ? true : false;  // Fixed Oct. 6
 
   // Pattern_phi / ME1_neigh_phi
   // FW: Pattern phi slightly offset from true LCT phi; also ME3/4 pattern width off
   //     Pattern phi of neighbor hits in ME1 miscalculated
   fixZonePhi_         = (fw_version < 47214) ? false : true;  // Fixed Oct. 6
 
   // LCT_station_2
   // FW: Reduced LCT matching window in station 2, resulting in demoted tracks and inefficiency
   bugSt2PhDiff_       = (47109 <= fw_version && fw_version < 47249) ? true : false;  // Bug introduced Oct. 6, fixed Oct. 19
 
   // LCT_theta_dup
   // FW: LCTs matched to track may take theta value from other LCT in the same chamber
   bugME11Dupes_       = (fw_version < 47423) ? true : false;  // Fixed Nov. 1
 
   // LCT_7_10_neg_pT (E)
   // LUT: Written with incorrect values for mode 7 CLCT, mode 10 random offset, all modes negative (1/pT) set to 3 instead of 511
   bugMode7CLCT_       = (fw_version < 47864) ? true : false;  // Fixed sometime after Nov. 1
   bugNegPt_           = (fw_version < 47864) ? true : false;  // Fixed sometime after Nov. 1
 
 
 }

void SectorProcessor::process	(	EventNumber_t	ievent,
		const TriggerPrimitiveCollection &	muon_primitives,
		EMTFHitCollection &	out_hits,
		EMTFTrackCollection &	out_tracks
	)		const

Definition at line 375 of file SectorProcessor.cc.

References bxWindow_, gather_cfg::cout, endcap_, maxBX_, minBX_, gen::n, process_single_bx(), sector_, findQualityFiles::size, and verbose_.

         {
 
   // if (endcap_ == 1 && sector_ == 1) {
   //   std::cout << "\nConfigured with era " << era_ << ", thetaWindowZone0 = " << thetaWindowZone0_ << ", bugAmbigThetaWin = "
   //          << bugAmbigThetaWin_ << ", twoStationSameBX = " << twoStationSameBX_ << ", promoteMode7_ = " << promoteMode7_ << std::endl;
   // }
 
   // List of converted hits, extended from previous BXs
   // deque (double-ended queue) is similar to a vector, but allows insertion or deletion of elements at both beginning and end
   std::deque<EMTFHitCollection> extended_conv_hits;
 
   // List of best track candidates, extended from previous BXs
   std::deque<EMTFTrackCollection> extended_best_track_cands;
 
   // Map of pattern detector --> lifetime, tracked across BXs
   std::map<pattern_ref_t, int> patt_lifetime_map;
 
   // ___________________________________________________________________________
   // Run each sector processor for every BX, taking into account the BX window
 
   int delayBX = bxWindow_ - 1;
 
   for (int bx = minBX_; bx <= maxBX_ + delayBX; ++bx) {
     if (verbose_ > 0) {  // debug
       std::cout << "Endcap: " << endcap_ << " Sector: " << sector_ << " Event: " << ievent << " BX: " << bx << std::endl;
     }
 
     process_single_bx(
         bx,
         muon_primitives,
         out_hits,
         out_tracks,
         extended_conv_hits,
         extended_best_track_cands,
         patt_lifetime_map
     );
 
     // Drop earliest BX outside of BX window
     if (bx >= minBX_ + delayBX) {
       extended_conv_hits.pop_front();
 
       int n = emtf::zone_array<int>().size();
       extended_best_track_cands.erase(extended_best_track_cands.end()-n, extended_best_track_cands.end());  // pop_back
     }
   }  // end loop over bx
 
   return;
 }

void SectorProcessor::process_single_bx	(	int	bx,
		const TriggerPrimitiveCollection &	muon_primitives,
		EMTFHitCollection &	out_hits,
		EMTFTrackCollection &	out_tracks,
		std::deque< EMTFHitCollection > &	extended_conv_hits,
		std::deque< EMTFTrackCollection > &	extended_best_track_cands,
		std::map< pattern_ref_t, int > &	patt_lifetime_map
	)		const

Definition at line 429 of file SectorProcessor.cc.

References bug9BitDPhi_, bugAmbigThetaWin_, bugGMTPhi_, bugME11Dupes_, bugMode7CLCT_, bugNegPt_, bugSameSectorPt0_, bugSt2PhDiff_, bxShiftCSC_, bxShiftGEM_, bxShiftRPC_, bxWindow_, BestTrackSelection::configure(), SingleHitTrack::configure(), PrimitiveSelection::configure(), AngleCalculation::configure(), PrimitiveConversion::configure(), PtAssignment::configure(), PrimitiveMatching::configure(), PatternRecognition::configure(), duplicateTheta_, endcap_, fixME11Edges_, fixMode15HighPt_, fixZonePhi_, includeNeighbor_, lut_, maxRoadsPerZone_, maxTracks_, PrimitiveSelection::merge(), PrimitiveSelection::merge_no_truncate(), modeQualVer_, pattDefinitions_, SingleHitTrack::process(), AngleCalculation::process(), BestTrackSelection::process(), PrimitiveSelection::process(), PrimitiveMatching::process(), PtAssignment::process(), PrimitiveConversion::process(), PatternRecognition::process(), promoteMode7_, pt_assign_engine_, readPtLUTFile_, sector_, symPattDefinitions_, thetaWindow_, thetaWindowZone0_, tp_geom_, twoStationSameBX_, useNewZones_, useRPC_, useSecondEarliest_, useSingleHits_, useSymPatterns_, verbose_, zoneBoundaries_, and zoneOverlap_.

Referenced by process().

         {
 
   // ___________________________________________________________________________
   // Configure
 
   PrimitiveSelection prim_sel;
   prim_sel.configure(
       verbose_, endcap_, sector_, bx,
       bxShiftCSC_, bxShiftRPC_, bxShiftGEM_,
       includeNeighbor_, duplicateTheta_,
       bugME11Dupes_
   );
 
   PrimitiveConversion prim_conv;
   prim_conv.configure(
       tp_geom_, lut_,
       verbose_, endcap_, sector_, bx,
       bxShiftCSC_, bxShiftRPC_, bxShiftGEM_,
       zoneBoundaries_, zoneOverlap_,
       duplicateTheta_, fixZonePhi_, useNewZones_, fixME11Edges_,
       bugME11Dupes_
   );
 
   PatternRecognition patt_recog;
   patt_recog.configure(
       verbose_, endcap_, sector_, bx,
       bxWindow_,
       pattDefinitions_, symPattDefinitions_, useSymPatterns_,
       maxRoadsPerZone_, useSecondEarliest_
   );
 
   PrimitiveMatching prim_match;
   prim_match.configure(
       verbose_, endcap_, sector_, bx,
       fixZonePhi_, useNewZones_,
       bugSt2PhDiff_, bugME11Dupes_
   );
 
   AngleCalculation angle_calc;
   angle_calc.configure(
       verbose_, endcap_, sector_, bx,
       bxWindow_,
       thetaWindow_, thetaWindowZone0_,
       bugME11Dupes_, bugAmbigThetaWin_, twoStationSameBX_
   );
 
   BestTrackSelection btrack_sel;
   btrack_sel.configure(
       verbose_, endcap_, sector_, bx,
       bxWindow_,
       maxRoadsPerZone_, maxTracks_, useSecondEarliest_,
       bugSameSectorPt0_
   );
 
   SingleHitTrack single_hit;
   single_hit.configure(
       verbose_, endcap_, sector_, bx,
       maxTracks_,
       useSingleHits_
   );
 
   PtAssignment pt_assign;
   pt_assign.configure(
       pt_assign_engine_,
       verbose_, endcap_, sector_, bx,
       readPtLUTFile_, fixMode15HighPt_,
       bug9BitDPhi_, bugMode7CLCT_, bugNegPt_,
       bugGMTPhi_, promoteMode7_, modeQualVer_
   );
 
   std::map<int, TriggerPrimitiveCollection> selected_csc_map;
   std::map<int, TriggerPrimitiveCollection> selected_rpc_map;
   std::map<int, TriggerPrimitiveCollection> selected_gem_map;
   std::map<int, TriggerPrimitiveCollection> selected_prim_map;
   std::map<int, TriggerPrimitiveCollection> inclusive_selected_prim_map;
 
   EMTFHitCollection conv_hits;  // "converted" hits converted by primitive converter
   EMTFHitCollection inclusive_conv_hits;
 
   emtf::zone_array<EMTFRoadCollection> zone_roads;  // each zone has its road collection
 
   emtf::zone_array<EMTFTrackCollection> zone_tracks;  // each zone has its track collection
 
   EMTFTrackCollection best_tracks;  // "best" tracks selected from all the zones
 
   // ___________________________________________________________________________
   // Process
 
   // Select muon primitives that belong to this sector and this BX.
   // Put them into maps with an index that roughly corresponds to
   // each input link.
   // From src/PrimitiveSelection.cc
   prim_sel.process(CSCTag(), muon_primitives, selected_csc_map);
   if (useRPC_) {
     prim_sel.process(RPCTag(), muon_primitives, selected_rpc_map);
   }
   prim_sel.process(GEMTag(), muon_primitives, selected_gem_map);
   prim_sel.merge(selected_csc_map, selected_rpc_map, selected_gem_map, selected_prim_map);
 
   // Convert trigger primitives into "converted" hits
   // A converted hit consists of integer representations of phi, theta, and zones
   // From src/PrimitiveConversion.cc
   prim_conv.process(selected_prim_map, conv_hits);
   extended_conv_hits.push_back(conv_hits);
 
   {
     // Keep all the converted hits for the use of data-emulator comparisons.
     // They include the extra ones that are not used in track building and the subsequent steps.
     prim_sel.merge_no_truncate(selected_csc_map, selected_rpc_map, selected_gem_map, inclusive_selected_prim_map);
     prim_conv.process(inclusive_selected_prim_map, inclusive_conv_hits);
 
     // Clear the input maps to save memory
     selected_csc_map.clear();
     selected_rpc_map.clear();
     selected_gem_map.clear();
   }
 
   // Detect patterns in all zones, find 3 best roads in each zone
   // From src/PatternRecognition.cc
   patt_recog.process(extended_conv_hits, patt_lifetime_map, zone_roads);
 
   // Match the trigger primitives to the roads, create tracks
   // From src/PrimitiveMatching.cc
   prim_match.process(extended_conv_hits, zone_roads, zone_tracks);
 
   // Calculate deflection angles for each track and fill track variables
   // From src/AngleCalculation.cc
   angle_calc.process(zone_tracks);
   extended_best_track_cands.insert(extended_best_track_cands.begin(), zone_tracks.begin(), zone_tracks.end());  // push_front
 
   // Select 3 "best" tracks from all the zones
   // From src/BestTrackSelection.cc
   btrack_sel.process(extended_best_track_cands, best_tracks);
 
   // Insert single LCTs from station 1 as tracks
   // From src/SingleHitTracks.cc
   single_hit.process(conv_hits, best_tracks);
 
   // Construct pT address, assign pT, calculate other GMT quantities
   // From src/PtAssignment.cc
   pt_assign.process(best_tracks);
 
   // ___________________________________________________________________________
   // Output
 
   out_hits.insert(out_hits.end(), inclusive_conv_hits.begin(), inclusive_conv_hits.end());
   out_tracks.insert(out_tracks.end(), best_tracks.begin(), best_tracks.end());
 
   return;
 }