#include <TrackFinder.h>

Public Member Functions
void	process (const edm::Event &iEvent, const edm::EventSetup &iSetup, EMTFHitCollection &out_hits, EMTFTrackCollection &out_tracks)

	TrackFinder (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iConsumes)

	~TrackFinder ()

Private Attributes
ConditionHelper	condition_helper_

const edm::ParameterSet	config_

std::string	era_

unsigned int	fw_version_

GeometryTranslator	geometry_translator_

unsigned int	pc_lut_version_

PtAssignmentEngine *	pt_assign_engine_

std::unique_ptr< PtAssignmentEngine >	pt_assign_engine_2016_

std::unique_ptr< PtAssignmentEngine >	pt_assign_engine_2017_

unsigned int	pt_lut_version_

SectorProcessorLUT	sector_processor_lut_

emtf::sector_array< SectorProcessor >	sector_processors_

const edm::EDGetToken	tokenCSC_

const edm::EDGetToken	tokenGEM_

const edm::EDGetToken	tokenRPC_

bool	useCSC_

bool	useGEM_

bool	useRPC_

int	verbose_

Detailed Description

Definition at line 17 of file TrackFinder.h.

Constructor & Destructor Documentation

TrackFinder::TrackFinder	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iConsumes
	)

explicit

Definition at line 11 of file TrackFinder.cc.

References condition_helper_, config_, makeMuonMisalignmentScenario::endcap, fw_version_, geometry_translator_, edm::ParameterSet::getParameter(), emtf::MAX_ENDCAP, emtf::MAX_TRIGSECTOR, simTrackMatching_cfi::maxBX, RecoTauCombinatoricProducer_cfi::maxTracks, emtf::MIN_ENDCAP, emtf::MIN_TRIGSECTOR, simTrackMatching_cfi::minBX, pc_lut_version_, pt_assign_engine_, pt_assign_engine_2016_, pt_assign_engine_2017_, pt_lut_version_, sector_processor_lut_, sector_processors_, and verbose_.

                                                                                          :
     geometry_translator_(),
     condition_helper_(),
     sector_processor_lut_(),
     pt_assign_engine_(),
     sector_processors_(),
     config_(iConfig),
     tokenCSC_(iConsumes.consumes<CSCTag::digi_collection>(iConfig.getParameter<edm::InputTag>("CSCInput"))),
     tokenRPC_(iConsumes.consumes<RPCTag::digi_collection>(iConfig.getParameter<edm::InputTag>("RPCInput"))),
     tokenGEM_(iConsumes.consumes<GEMTag::digi_collection>(iConfig.getParameter<edm::InputTag>("GEMInput"))),
     verbose_(iConfig.getUntrackedParameter<int>("verbosity")),
     useCSC_(iConfig.getParameter<bool>("CSCEnable")),
     useRPC_(iConfig.getParameter<bool>("RPCEnable")),
     useGEM_(iConfig.getParameter<bool>("GEMEnable"))
 {
 
   auto minBX       = iConfig.getParameter<int>("MinBX");
   auto maxBX       = iConfig.getParameter<int>("MaxBX");
   auto bxWindow    = iConfig.getParameter<int>("BXWindow");
   auto bxShiftCSC  = iConfig.getParameter<int>("CSCInputBXShift");
   auto bxShiftRPC  = iConfig.getParameter<int>("RPCInputBXShift");
   auto bxShiftGEM  = iConfig.getParameter<int>("GEMInputBXShift");
 
   const auto& spPCParams16 = config_.getParameter<edm::ParameterSet>("spPCParams16");
   auto zoneBoundaries     = spPCParams16.getParameter<std::vector<int> >("ZoneBoundaries");
   auto zoneOverlap        = spPCParams16.getParameter<int>("ZoneOverlap");
   auto includeNeighbor    = spPCParams16.getParameter<bool>("IncludeNeighbor");
   auto duplicateTheta     = spPCParams16.getParameter<bool>("DuplicateTheta");
   auto fixZonePhi         = spPCParams16.getParameter<bool>("FixZonePhi");
   auto useNewZones        = spPCParams16.getParameter<bool>("UseNewZones");
   auto fixME11Edges       = spPCParams16.getParameter<bool>("FixME11Edges");
 
   const auto& spPRParams16 = config_.getParameter<edm::ParameterSet>("spPRParams16");
   auto pattDefinitions    = spPRParams16.getParameter<std::vector<std::string> >("PatternDefinitions");
   auto symPattDefinitions = spPRParams16.getParameter<std::vector<std::string> >("SymPatternDefinitions");
   auto useSymPatterns     = spPRParams16.getParameter<bool>("UseSymmetricalPatterns");
 
   const auto& spTBParams16 = config_.getParameter<edm::ParameterSet>("spTBParams16");
   auto thetaWindow        = spTBParams16.getParameter<int>("ThetaWindow");
   auto thetaWindowRPC     = spTBParams16.getParameter<int>("ThetaWindowRPC");
   auto useSingleHits      = spTBParams16.getParameter<bool>("UseSingleHits");
   auto bugSt2PhDiff       = spTBParams16.getParameter<bool>("BugSt2PhDiff");
   auto bugME11Dupes       = spTBParams16.getParameter<bool>("BugME11Dupes");
 
   const auto& spGCParams16 = config_.getParameter<edm::ParameterSet>("spGCParams16");
   auto maxRoadsPerZone    = spGCParams16.getParameter<int>("MaxRoadsPerZone");
   auto maxTracks          = spGCParams16.getParameter<int>("MaxTracks");
   auto useSecondEarliest  = spGCParams16.getParameter<bool>("UseSecondEarliest");
   auto bugSameSectorPt0   = spGCParams16.getParameter<bool>("BugSameSectorPt0");
 
   const auto& spPAParams16 = config_.getParameter<edm::ParameterSet>("spPAParams16");
   auto ptLUTVersion       = spPAParams16.getParameter<int>("PtLUTVersion");
   auto readPtLUTFile      = spPAParams16.getParameter<bool>("ReadPtLUTFile");
   auto fixMode15HighPt    = spPAParams16.getParameter<bool>("FixMode15HighPt");
   auto bug9BitDPhi        = spPAParams16.getParameter<bool>("Bug9BitDPhi");
   auto bugMode7CLCT       = spPAParams16.getParameter<bool>("BugMode7CLCT");
   auto bugNegPt           = spPAParams16.getParameter<bool>("BugNegPt");
   auto bugGMTPhi          = spPAParams16.getParameter<bool>("BugGMTPhi");
   auto promoteMode7       = spPAParams16.getParameter<bool>("PromoteMode7");
 
   pt_assign_engine_2016_.reset(new PtAssignmentEngine2016());
   pt_assign_engine_2017_.reset(new PtAssignmentEngine2017());
 
   fw_version_     = 999999;
   pt_lut_version_ = 888888;
   pc_lut_version_ = 777777;
 
   // Configure sector processors
   for (int endcap = emtf::MIN_ENDCAP; endcap <= emtf::MAX_ENDCAP; ++endcap) {
     for (int sector = emtf::MIN_TRIGSECTOR; sector <= emtf::MAX_TRIGSECTOR; ++sector) {
       const int es = (endcap - emtf::MIN_ENDCAP) * (emtf::MAX_TRIGSECTOR - emtf::MIN_TRIGSECTOR + 1) + (sector - emtf::MIN_TRIGSECTOR);
       
       sector_processors_.at(es).configure(
           &geometry_translator_,
           &condition_helper_,
           &sector_processor_lut_,
           &pt_assign_engine_,
           verbose_, endcap, sector,
           minBX, maxBX, bxWindow, bxShiftCSC, bxShiftRPC, bxShiftGEM,
           zoneBoundaries, zoneOverlap,
           includeNeighbor, duplicateTheta, fixZonePhi, useNewZones, fixME11Edges,
           pattDefinitions, symPattDefinitions, useSymPatterns,
           thetaWindow, thetaWindowRPC, useSingleHits, bugSt2PhDiff, bugME11Dupes,
           maxRoadsPerZone, maxTracks, useSecondEarliest, bugSameSectorPt0,
           ptLUTVersion, readPtLUTFile, fixMode15HighPt, bug9BitDPhi, bugMode7CLCT, bugNegPt, bugGMTPhi, promoteMode7
       );
     }
   }
 
 } // End constructor: TrackFinder::TrackFinder()

TrackFinder::~TrackFinder ( )

Definition at line 102 of file TrackFinder.cc.

                           {
 
 }

Member Function Documentation

void TrackFinder::process	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		EMTFHitCollection &	out_hits,
		EMTFTrackCollection &	out_tracks
	)

Definition at line 106 of file TrackFinder.cc.

References relativeConstraints::chamber, ConditionHelper::checkAndUpdateConditions(), L1TMuon::GeometryTranslator::checkAndUpdateGeometry(), condition_helper_, gather_cfg::cout, makeMuonMisalignmentScenario::endcap, edm::EventID::event(), EMTFSubsystemCollector::extractPrimitives(), fw_version_, geometry_translator_, ConditionHelper::get_fw_version(), ConditionHelper::get_pc_lut_version(), ConditionHelper::get_pt_lut_version(), ConditionHelper::getForest(), edm::EventBase::id(), edm::EventBase::isRealData(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kRPC, PtAssignmentEngine::load(), emtf::MAX_ENDCAP, emtf::MAX_TRIGSECTOR, emtf::MIN_ENDCAP, emtf::MIN_TRIGSECTOR, AlCaHLTBitMon_ParallelJobs::p, pc_lut_version_, pt_assign_engine_, pt_assign_engine_2016_, pt_assign_engine_2017_, pt_lut_version_, SectorProcessorLUT::read(), sector_processor_lut_, sector_processors_, PtAssignmentEngine::set_ptLUTVersion(), relativeConstraints::station, digitizers_cfi::strip, lumiQTWidget::t, tokenCSC_, tokenGEM_, tokenRPC_, useCSC_, useGEM_, useRPC_, and verbose_.

   {
 
   // Clear output collections
   out_hits.clear();
   out_tracks.clear();
 
   // Get the geometry for TP conversions
   geometry_translator_.checkAndUpdateGeometry(iSetup);
 
   // ___________________________________________________________________________
   // Get the conditions: firmware version, PC LUT version, pT BDT forests
   bool new_conditions = condition_helper_.checkAndUpdateConditions(iEvent, iSetup);
   // std::cout << "\nDo we have new conditions? Yes (1) or no (0)? -- " << new_conditions << std::endl;
   // std::cout << "  * Before updating conditions, fw_version_ = " << fw_version_ << ", pt_lut_version_ = " << pt_lut_version_ << ", pc_lut_version_ = " << pc_lut_version_ << std::endl;
 
   if (new_conditions) {
 
     fw_version_ = condition_helper_.get_fw_version();
     // No RPC or GEM hits in 2016
     if (fw_version_ != 0 && fw_version_ < 50000) {
       useRPC_ = false;
       useGEM_ = false;
     }
 
     pt_lut_version_ = condition_helper_.get_pt_lut_version();
     pc_lut_version_ = condition_helper_.get_pc_lut_version();
 
     // Reload primitive conversion LUTs if necessary
     // std::cout << "Configured with pc_lut_version_ = " << pc_lut_version_ << std::endl;
     sector_processor_lut_.read( pc_lut_version_ );
 
     // std::cout << "Configured with pt_lut_version_ = " << pt_lut_version_ << std::endl;
     if ( pt_lut_version_ <= 5 ) {
       pt_assign_engine_ = pt_assign_engine_2016_.get();
       pt_assign_engine_->set_ptLUTVersion( pt_lut_version_ );
     } else {
       pt_assign_engine_ = pt_assign_engine_2017_.get();
       pt_assign_engine_->set_ptLUTVersion( pt_lut_version_ );
     }
     // Reload pT LUT if necessary
     pt_assign_engine_->load( &(condition_helper_.getForest()) );
   }
 
   // std::cout << "  * After changes, fw_version_ = " << fw_version_ << ", pt_lut_version_ = " << pt_lut_version_ << ", pc_lut_version_ = " << pc_lut_version_ << std::endl;
 
   // ___________________________________________________________________________
   // Extract all trigger primitives
 
   TriggerPrimitiveCollection muon_primitives;
 
   EMTFSubsystemCollector collector;
   if (useCSC_)
     collector.extractPrimitives(CSCTag(), iEvent, tokenCSC_, muon_primitives);
   if (useRPC_)
     collector.extractPrimitives(RPCTag(), iEvent, tokenRPC_, muon_primitives);
   if (useGEM_)
     collector.extractPrimitives(GEMTag(), iEvent, tokenGEM_, muon_primitives);
 
   // Check trigger primitives
   if (verbose_ > 2) {  // debug
     std::cout << "Num of TriggerPrimitive: " << muon_primitives.size() << std::endl;
     for (const auto& p : muon_primitives) {
       p.print(std::cout);
     }
   }
 
   // ___________________________________________________________________________
   // Run each sector processor
 
   // MIN/MAX ENDCAP and TRIGSECTOR set in interface/Common.h
   for (int endcap = emtf::MIN_ENDCAP; endcap <= emtf::MAX_ENDCAP; ++endcap) {
     for (int sector = emtf::MIN_TRIGSECTOR; sector <= emtf::MAX_TRIGSECTOR; ++sector) {
       const int es = (endcap - emtf::MIN_ENDCAP) * (emtf::MAX_TRIGSECTOR - emtf::MIN_TRIGSECTOR + 1) + (sector - emtf::MIN_TRIGSECTOR);
 
       // Run-dependent configure. This overwrites many of the configurables passed by the python config file.
       if (new_conditions) {
     if (iEvent.isRealData()) {
       sector_processors_.at(es).configure_by_fw_version( fw_version_ );
     }
     sector_processors_.at(es).set_pt_lut_version( pt_lut_version_ );
       }
 
       // Process
       sector_processors_.at(es).process(
           iEvent.id().event(),
           muon_primitives,
           out_hits,
           out_tracks
       );
     }
   }
 
 
   // ___________________________________________________________________________
   // Check emulator input and output. They are printed in a way that is friendly
   // for comparison with the firmware simulator.
 
   if (verbose_ > 0) {  // debug
 
     for (int endcap = emtf::MIN_ENDCAP; endcap <= emtf::MAX_ENDCAP; ++endcap) {
       for (int sector = emtf::MIN_TRIGSECTOR; sector <= emtf::MAX_TRIGSECTOR; ++sector) {
         const int es = (endcap - emtf::MIN_ENDCAP) * (emtf::MAX_TRIGSECTOR - emtf::MIN_TRIGSECTOR + 1) + (sector - emtf::MIN_TRIGSECTOR);
 
         // _____________________________________________________________________
         // This prints the hits as raw text input to the firmware simulator
         // "12345" is the BX separator
 
         std::cout << "==== Endcap " << endcap << " Sector " << sector << " Hits ====" << std::endl;
         std::cout << "bx e s ss st vf ql cp wg id bd hs" << std::endl;
 
         bool empty_sector = true;
         for (const auto& h : out_hits) {
           if (h.Sector_idx() != es)  continue;
           empty_sector = false;
         }
 
         for (int ibx = -3-5; (ibx < +3+5+5) && !empty_sector; ++ibx) {
 
           for (const auto& h : out_hits) {
             if (h.Subsystem() == TriggerPrimitive::kCSC) {
               if (h.Sector_idx() != es)  continue;
               if (h.BX() != ibx)  continue;
 
               int bx        = 1;
               int endcap    = (h.Endcap() == 1) ? 1 : 2;
               int sector    = h.PC_sector();
               int station   = (h.PC_station() == 0 && h.Subsector() == 1) ? 1 : h.PC_station();
               int chamber   = h.PC_chamber() + 1;
               int strip     = (h.Station() == 1 && h.Ring() == 4) ? h.Strip() + 128 : h.Strip();  // ME1/1a
               int wire      = h.Wire();
               int valid     = 1;
               std::cout << bx << " " << endcap << " " << sector << " " << h.Subsector() << " "
                   << station << " " << valid << " " << h.Quality() << " " << h.Pattern() << " "
                   << wire << " " << chamber << " " << h.Bend() << " " << strip << std::endl;
 
             } else if (h.Subsystem() == TriggerPrimitive::kRPC) {
               if (h.Sector_idx() != es)  continue;
               if (h.BX()+5 != ibx)  continue;  // RPC hits should be supplied 5 BX later relative to CSC hits
 
               // Assign RPC link index. Code taken from src/PrimitiveSelection.cc
               int rpc_sub = -1;
               int rpc_chm = -1;
               if (!h.Neighbor()) {
                 rpc_sub = h.Subsector() - 1;
               } else {
                 rpc_sub = 6;
               }
               if (h.Station() <= 2) {
                 rpc_chm = (h.Station() - 1);
               } else {
                 rpc_chm = 2 + (h.Station() - 3)*2 + (h.Ring() - 2);
               }
 
               int bx        = 1;
               int endcap    = (h.Endcap() == 1) ? 1 : 2;
               int sector    = h.PC_sector();
               int station   = rpc_sub;
               int chamber   = rpc_chm + 1;
               int strip     = (h.Phi_fp() >> 2);
               int wire      = (h.Theta_fp() >> 2);
               int valid     = 2;  // this marks RPC stub
               std::cout << bx << " " << endcap << " " << sector << " " << 0 << " "
                   << station << " " << valid << " " << 0 << " " << 0 << " "
                   << wire << " " << chamber << " " << 0 << " " << strip << std::endl;
             }
           }  // end loop over hits
 
           std::cout << "12345" << std::endl;
         }  // end loop over bx
 
         // _____________________________________________________________________
         // This prints the tracks as raw text output from the firmware simulator
 
         std::cout << "==== Endcap " << endcap << " Sector " << sector << " Tracks ====" << std::endl;
         std::cout << "bx e s a mo et ph cr q pt" << std::endl;
 
         for (const auto& t : out_tracks) {
           if (t.Sector_idx() != es)  continue;
 
           std::cout << t.BX() << " " << (t.Endcap() == 1 ? 1 : 2) << " " << t.Sector() << " " << t.PtLUT().address << " " << t.Mode() << " "
               << (t.GMT_eta() >= 0 ? t.GMT_eta() : t.GMT_eta()+512)<< " " << t.GMT_phi() << " "
               << t.GMT_charge() << " " << t.GMT_quality() << " " << t.Pt() << std::endl;
         }  // end loop over tracks
 
       }  // end loop over sector
     }  // end loop over endcap
   }  // end debug
 
   return;
 }

Member Data Documentation

ConditionHelper TrackFinder::condition_helper_

private

Definition at line 33 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

const edm::ParameterSet TrackFinder::config_

private

Definition at line 46 of file TrackFinder.h.

Referenced by TrackFinder().

std::string TrackFinder::era_

private

Definition at line 54 of file TrackFinder.h.

unsigned int TrackFinder::fw_version_

private

Definition at line 34 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

GeometryTranslator TrackFinder::geometry_translator_

private

Definition at line 31 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

unsigned int TrackFinder::pc_lut_version_

private

Definition at line 36 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

PtAssignmentEngine* TrackFinder::pt_assign_engine_

private

Definition at line 40 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

std::unique_ptr<PtAssignmentEngine> TrackFinder::pt_assign_engine_2016_

private

Definition at line 41 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

std::unique_ptr<PtAssignmentEngine> TrackFinder::pt_assign_engine_2017_

private

Definition at line 42 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

unsigned int TrackFinder::pt_lut_version_

private

Definition at line 35 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

SectorProcessorLUT TrackFinder::sector_processor_lut_

private

Definition at line 38 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

emtf::sector_array<SectorProcessor> TrackFinder::sector_processors_

private

Definition at line 44 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

const edm::EDGetToken TrackFinder::tokenCSC_

private

Definition at line 48 of file TrackFinder.h.

Referenced by process().

const edm::EDGetToken TrackFinder::tokenGEM_

private

Definition at line 48 of file TrackFinder.h.

Referenced by process().

const edm::EDGetToken TrackFinder::tokenRPC_

private

Definition at line 48 of file TrackFinder.h.

Referenced by process().

bool TrackFinder::useCSC_

private

Definition at line 52 of file TrackFinder.h.

Referenced by process().

bool TrackFinder::useGEM_

private

Definition at line 52 of file TrackFinder.h.

Referenced by process().

bool TrackFinder::useRPC_

private

Definition at line 52 of file TrackFinder.h.

Referenced by process().

int TrackFinder::verbose_

private

Definition at line 50 of file TrackFinder.h.

Referenced by process(), and TrackFinder().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation