SectorProcessor Class Reference

#include <SectorProcessor.h>

Public Types
typedef PatternRecognition::pattern_ref_t	pattern_ref_t

Public Member Functions
void	configure (const EMTFSetup *setup, int verbose, int endcap, int sector)
void	process (const edm::EventID &event_id, 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
int	endcap_
int	sector_
const EMTFSetup *	setup_
int	verbose_

Detailed Description

Definition at line 23 of file SectorProcessor.h.

Member Typedef Documentation

pattern_ref_t

typedef PatternRecognition::pattern_ref_t SectorProcessor::pattern_ref_t

Definition at line 28 of file SectorProcessor.h.

Constructor & Destructor Documentation

SectorProcessor()

SectorProcessor::SectorProcessor ( )

explicit

Definition at line 3 of file SectorProcessor.cc.

{}

~SectorProcessor()

SectorProcessor::~SectorProcessor

(

)

Definition at line 5 of file SectorProcessor.cc.

{}

Member Function Documentation

configure()

void SectorProcessor::configure	(	const EMTFSetup *	setup,
		int	verbose,
		int	endcap,
		int	sector
	)

Definition at line 7 of file SectorProcessor.cc.

References emtf_assert, makeMuonMisalignmentScenario::endcap, endcap_, emtf::MAX_ENDCAP, emtf::MAX_TRIGSECTOR, emtf::MIN_ENDCAP, emtf::MIN_TRIGSECTOR, sector_, singleTopDQM_cfi::setup, setup_, verbose, and verbose_.

                                                                                           {
  emtf_assert(setup != nullptr);
  emtf_assert(emtf::MIN_ENDCAP <= endcap && endcap <= emtf::MAX_ENDCAP);
  emtf_assert(emtf::MIN_TRIGSECTOR <= sector && sector <= emtf::MAX_TRIGSECTOR);

  setup_ = setup;
  verbose_ = verbose;
  endcap_ = endcap;
  sector_ = sector;
}

process()

void SectorProcessor::process	(	const edm::EventID &	event_id,
		const TriggerPrimitiveCollection &	muon_primitives,
		EMTFHitCollection &	out_hits,
		EMTFTrackCollection &	out_tracks
	)		const

Definition at line 18 of file SectorProcessor.cc.

References simKBmtfDigis_cfi::bx, looper::cfg, gather_cfg::cout, endcap_, EMTFSetup::getVersionControl(), dqmiodumpmetadata::n, process_single_bx(), sector_, setup_, findQualityFiles::size, and verbose_.

                                                                     {
  auto cfg = setup_->getVersionControl();

  // 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 = cfg.bxWindow_ - 1;

  for (int bx = cfg.minBX_; bx <= cfg.maxBX_ + delayBX; ++bx) {
    if (verbose_ > 0) {  // debug
      std::cout << "Event: " << event_id << " Endcap: " << endcap_ << " Sector: " << sector_ << " 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 >= cfg.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;
}

process_single_bx()

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 61 of file SectorProcessor.cc.

References simKBmtfDigis_cfi::bx, looper::cfg, PrimitiveSelection::configure(), AngleCalculation::configure(), BestTrackSelection::configure(), SingleHitTrack::configure(), PrimitiveConversion::configure(), PrimitiveMatching::configure(), PtAssignment::configure(), PatternRecognition::configure(), endcap_, EMTFSetup::getGeometryTranslator(), EMTFSetup::getPtAssignmentEngine(), EMTFSetup::getPtAssignmentEngineDxy(), EMTFSetup::getSectorProcessorLUT(), EMTFSetup::getVersionControl(), PrimitiveSelection::merge(), PrimitiveSelection::merge_no_truncate(), SingleHitTrack::process(), BestTrackSelection::process(), AngleCalculation::process(), PrimitiveMatching::process(), PrimitiveSelection::process(), PatternRecognition::process(), PtAssignment::process(), PrimitiveConversion::process(), sector_, setup_, and verbose_.

Referenced by process().

                                                                                             {
  auto cfg = setup_->getVersionControl();

  auto tp_geom_ = &(setup_->getGeometryTranslator());
  auto pc_lut_ = &(setup_->getSectorProcessorLUT());
  auto pt_assign_engine_ = setup_->getPtAssignmentEngine();
  auto pt_assign_engine_dxy_ = setup_->getPtAssignmentEngineDxy();

  // ___________________________________________________________________________
  // Configure

  PrimitiveSelection prim_sel;
  prim_sel.configure(verbose_,
                     endcap_,
                     sector_,
                     bx,
                     cfg.bxShiftCSC_,
                     cfg.bxShiftRPC_,
                     cfg.bxShiftGEM_,
                     cfg.bxShiftME0_,
                     cfg.includeNeighbor_,
                     cfg.duplicateTheta_,
                     cfg.bugME11Dupes_);

  PrimitiveConversion prim_conv;
  prim_conv.configure(tp_geom_,
                      pc_lut_,
                      verbose_,
                      endcap_,
                      sector_,
                      bx,
                      cfg.bxShiftCSC_,
                      cfg.bxShiftRPC_,
                      cfg.bxShiftGEM_,
                      cfg.bxShiftME0_,
                      cfg.zoneBoundaries_,
                      cfg.zoneOverlap_,
                      cfg.duplicateTheta_,
                      cfg.fixZonePhi_,
                      cfg.useNewZones_,
                      cfg.fixME11Edges_,
                      cfg.bugME11Dupes_,
                      cfg.useRun3CCLUT_OTMB_,
                      cfg.useRun3CCLUT_TMB_);

  PatternRecognition patt_recog;
  patt_recog.configure(verbose_,
                       endcap_,
                       sector_,
                       bx,
                       cfg.bxWindow_,
                       cfg.pattDefinitions_,
                       cfg.symPattDefinitions_,
                       cfg.useSymPatterns_,
                       cfg.maxRoadsPerZone_,
                       cfg.useSecondEarliest_);

  PrimitiveMatching prim_match;
  prim_match.configure(
      verbose_, endcap_, sector_, bx, cfg.fixZonePhi_, cfg.useNewZones_, cfg.bugSt2PhDiff_, cfg.bugME11Dupes_);

  AngleCalculation angle_calc;
  angle_calc.configure(verbose_,
                       endcap_,
                       sector_,
                       bx,
                       cfg.bxWindow_,
                       cfg.thetaWindow_,
                       cfg.thetaWindowZone0_,
                       cfg.bugME11Dupes_,
                       cfg.bugAmbigThetaWin_,
                       cfg.twoStationSameBX_);

  BestTrackSelection btrack_sel;
  btrack_sel.configure(verbose_,
                       endcap_,
                       sector_,
                       bx,
                       cfg.bxWindow_,
                       cfg.maxRoadsPerZone_,
                       cfg.maxTracks_,
                       cfg.useSecondEarliest_,
                       cfg.bugSameSectorPt0_);

  SingleHitTrack single_hit;
  single_hit.configure(verbose_, endcap_, sector_, bx, cfg.maxTracks_, cfg.useSingleHits_);

  PtAssignment pt_assign;
  pt_assign.configure(pt_assign_engine_,
                      pt_assign_engine_dxy_,
                      verbose_,
                      endcap_,
                      sector_,
                      bx,
                      cfg.readPtLUTFile_,
                      cfg.fixMode15HighPt_,
                      cfg.bug9BitDPhi_,
                      cfg.bugMode7CLCT_,
                      cfg.bugNegPt_,
                      cfg.bugGMTPhi_,
                      cfg.promoteMode7_,
                      cfg.modeQualVer_,
                      cfg.pbFileName_);

  std::map<int, TriggerPrimitiveCollection> selected_dt_map;
  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_me0_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(emtf::DTTag(), muon_primitives, selected_dt_map);
  prim_sel.process(emtf::CSCTag(), muon_primitives, selected_csc_map);
  prim_sel.process(emtf::RPCTag(), muon_primitives, selected_rpc_map);
  prim_sel.process(emtf::GEMTag(), muon_primitives, selected_gem_map);
  prim_sel.process(emtf::ME0Tag(), muon_primitives, selected_me0_map);
  prim_sel.merge(
      selected_dt_map, selected_csc_map, selected_rpc_map, selected_gem_map, selected_me0_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_dt_map,
                               selected_csc_map,
                               selected_rpc_map,
                               selected_gem_map,
                               selected_me0_map,
                               inclusive_selected_prim_map);
    prim_conv.process(inclusive_selected_prim_map, inclusive_conv_hits);

    // Clear the input maps to save memory
    selected_dt_map.clear();
    selected_csc_map.clear();
    selected_rpc_map.clear();
    selected_gem_map.clear();
    selected_me0_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;
}

Member Data Documentation

endcap_

int SectorProcessor::endcap_

private

Definition at line 55 of file SectorProcessor.h.

Referenced by configure(), process(), and process_single_bx().

sector_

int SectorProcessor::sector_

private

Definition at line 55 of file SectorProcessor.h.

Referenced by configure(), process(), and process_single_bx().

setup_

const EMTFSetup* SectorProcessor::setup_

private

Definition at line 53 of file SectorProcessor.h.

Referenced by configure(), process(), process_single_bx(), and upgradeWorkflowComponents.UpgradeWorkflow::setup().

verbose_

int SectorProcessor::verbose_

private

Definition at line 55 of file SectorProcessor.h.

Referenced by configure(), process(), and process_single_bx().