#include <OMTFReconstruction.h>

Public Member Functions
void	beginJob ()

void	beginRun (edm::Run const &, edm::EventSetup const &)

void	endJob ()

	OMTFReconstruction (const edm::ParameterSet &, edm::ConsumesCollector &&)

std::unique_ptr< l1t::RegionalMuonCandBxCollection >	reconstruct (const edm::Event &, const edm::EventSetup &)

	~OMTFReconstruction ()

Private Member Functions
void	getProcessorCandidates (unsigned int iProcessor, l1t::tftype mtfType, int bx, l1t::RegionalMuonCandBxCollection &myCandidates)

void	loadAndFilterDigis (const edm::Event &)

void	writeResultToXML (unsigned int iProcessor, l1t::tftype mtfType, const OMTFinput &myInput, const std::vector< OMTFProcessor::resultsMap > &myResults, const std::vector< l1t::RegionalMuonCand > &candMuons)

Private Attributes
xercesc::DOMElement *	aTopElement

int	bxMax

int	bxMin

edm::Handle< CSCCorrelatedLCTDigiCollection >	cscDigis

edm::Handle< L1MuDTChambPhContainer >	dtPhDigis

edm::Handle< L1MuDTChambThContainer >	dtThDigis

bool	dumpDetailedResultToXML

bool	dumpResultToXML

edm::ESGetToken< L1TMuonOverlapParams, L1TMuonOverlapParamsRcd >	l1TMuonOverlapParamsToken_

edm::ParameterSet	m_Config

std::unique_ptr< IGhostBuster >	m_GhostBuster

OMTFinputMaker	m_InputMaker

OMTFProcessor *	m_OMTF

OMTFConfiguration *	m_OMTFConfig
	OMTF objects. More...

OMTFConfigMaker *	m_OMTFConfigMaker

OMTFSorter	m_Sorter

XMLConfigWriter *	m_Writer

edm::Handle< RPCDigiCollection >	rpcDigis

Detailed Description

Definition at line 37 of file OMTFReconstruction.h.

Constructor & Destructor Documentation

◆ OMTFReconstruction()

OMTFReconstruction::OMTFReconstruction	(	const edm::ParameterSet &	theConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 21 of file OMTFReconstruction.cc.

     : m_Config(theConfig),
       l1TMuonOverlapParamsToken_(
           iC.esConsumes<L1TMuonOverlapParams, L1TMuonOverlapParamsRcd, edm::Transition::BeginRun>()),
       m_OMTFConfig(nullptr),
       m_InputMaker(iC, false),
       m_OMTF(nullptr),
       aTopElement(nullptr),
       m_OMTFConfigMaker(nullptr),
       m_Writer(nullptr) {
   dumpResultToXML = m_Config.getParameter<bool>("dumpResultToXML");
   dumpDetailedResultToXML = m_Config.getParameter<bool>("dumpDetailedResultToXML");
   //m_Config.getParameter<std::string>("XMLDumpFileName");
   bxMin = m_Config.exists("bxMin") ? m_Config.getParameter<int>("bxMin") : 0;
   bxMax = m_Config.exists("bxMax") ? m_Config.getParameter<int>("bxMax") : 0;
 }

References bxMax, bxMin, dumpDetailedResultToXML, dumpResultToXML, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and m_Config.

◆ ~OMTFReconstruction()

OMTFReconstruction::~OMTFReconstruction ( )

Definition at line 39 of file OMTFReconstruction.cc.

                                         {
   delete m_OMTFConfig;
   delete m_OMTF;
  
   if (m_Writer)
     delete m_Writer;
 }

References m_OMTF, m_OMTFConfig, and m_Writer.

Member Function Documentation

◆ beginJob()

void OMTFReconstruction::beginJob ( void )

Definition at line 49 of file OMTFReconstruction.cc.

                                   {
   m_OMTFConfig = new OMTFConfiguration();
   m_OMTF = new OMTFProcessor();
 }

References m_OMTF, and m_OMTFConfig.

Referenced by L1TMuonOverlapTrackProducer::beginJob().

◆ beginRun()

void OMTFReconstruction::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

Definition at line 63 of file OMTFReconstruction.cc.

                                                                             {
   const L1TMuonOverlapParams* omtfParams = &iSetup.getData(l1TMuonOverlapParamsToken_);
  
   if (!omtfParams) {
     edm::LogError("L1TMuonOverlapTrackProducer") << "Could not retrieve parameters from Event Setup" << std::endl;
   }
  
   m_OMTFConfig->configure(omtfParams);
   m_OMTF->configure(m_OMTFConfig, omtfParams);
   //m_GhostBuster.setNphiBins(m_OMTFConfig->nPhiBins());
  
   if (m_OMTFConfig->fwVersion() >= 5) {
     //  if(m_Config.exists("ghostBusterType") ) {
     //    if(m_Config.getParameter<std::string>("ghostBusterType") == "GhostBusterPreferRefDt")
     m_GhostBuster = std::make_unique<GhostBusterPreferRefDt>(m_OMTFConfig);
   } else {
     m_GhostBuster = std::make_unique<GhostBuster>();
   }
  
   m_Sorter.initialize(m_OMTFConfig);
   m_Sorter.setNphiBins(m_OMTFConfig->nPhiBins());
  
   m_InputMaker.initialize(iSetup, m_OMTFConfig);
  
   if (dumpResultToXML) {
     m_Writer = new XMLConfigWriter(m_OMTFConfig);
     std::string fName = "OMTF";
     m_Writer->initialiseXMLDocument(fName);
   }
 }

References OMTFProcessor::configure(), OMTFConfiguration::configure(), dumpResultToXML, MainPageGenerator::fName, OMTFConfiguration::fwVersion(), edm::EventSetup::getData(), XMLConfigWriter::initialiseXMLDocument(), OMTFSorter::initialize(), OMTFinputMaker::initialize(), l1TMuonOverlapParamsToken_, m_GhostBuster, m_InputMaker, m_OMTF, m_OMTFConfig, m_Sorter, m_Writer, OMTFConfiguration::nPhiBins(), fakeOmtfParams_cff::omtfParams, OMTFSorter::setNphiBins(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by L1TMuonOverlapTrackProducer::beginRun().

◆ endJob()

void OMTFReconstruction::endJob ( void )

Definition at line 55 of file OMTFReconstruction.cc.

                                 {
   if (dumpResultToXML) {
     std::string fName = m_Config.getParameter<std::string>("XMLDumpFileName");
     m_Writer->finaliseXMLDocument(fName);
   }
 }

References dumpResultToXML, XMLConfigWriter::finaliseXMLDocument(), MainPageGenerator::fName, edm::ParameterSet::getParameter(), m_Config, m_Writer, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by L1TMuonOverlapTrackProducer::endJob().

◆ getProcessorCandidates()

void OMTFReconstruction::getProcessorCandidates	(	unsigned int	iProcessor,
		l1t::tftype	mtfType,
		int	bx,
		l1t::RegionalMuonCandBxCollection &	myCandidates
	)

private

Definition at line 137 of file OMTFReconstruction.cc.

                                                                                                  {
   m_InputMaker.setFlag(0);
   OMTFinput input = m_InputMaker.buildInputForProcessor(
       dtPhDigis.product(), dtThDigis.product(), cscDigis.product(), rpcDigis.product(), iProcessor, mtfType, bx);
   int flag = m_InputMaker.getFlag();
  
   const std::vector<OMTFProcessor::resultsMap>& results = m_OMTF->processInput(iProcessor, input);
  
   std::vector<AlgoMuon> algoCandidates;
  
   m_Sorter.sortRefHitResults(results, algoCandidates);
  
   // perform GB
   std::vector<AlgoMuon> gbCandidates = m_GhostBuster->select(algoCandidates);
  
   // fill RegionalMuonCand colleciton
   std::vector<l1t::RegionalMuonCand> candMuons = m_Sorter.candidates(iProcessor, mtfType, gbCandidates);
  
   //fill outgoing collection
   for (auto& candMuon : candMuons) {
     candMuon.setHwQual(candMuon.hwQual() | flag);  //FIXME temporary debug fix
     omtfCandidates.push_back(bx, candMuon);
   }
  
   //dump to XML
   if (bx == 0)
     writeResultToXML(iProcessor, mtfType, input, results, candMuons);
 }

References OMTFinputMaker::buildInputForProcessor(), l1GtPatternGenerator_cfi::bx, OMTFSorter::candidates(), cscDigis, dtPhDigis, dtThDigis, RemoveAddSevLevel::flag, OMTFinputMaker::getFlag(), input, m_GhostBuster, m_InputMaker, m_OMTF, m_Sorter, OMTFProcessor::processInput(), edm::Handle< T >::product(), BXVector< T >::push_back(), bookConverter::results, rpcDigis, OMTFinputMaker::setFlag(), OMTFSorter::sortRefHitResults(), and writeResultToXML().

Referenced by reconstruct().

◆ loadAndFilterDigis()

void OMTFReconstruction::loadAndFilterDigis ( const edm::Event & iEvent )

private

Definition at line 123 of file OMTFReconstruction.cc.

                                                                   {
   // Filter digis by dropping digis from selected (by cfg.py) subsystems
   if (!m_Config.getParameter<bool>("dropDTPrimitives")) {
     iEvent.getByLabel(m_Config.getParameter<edm::InputTag>("srcDTPh"), dtPhDigis);
     iEvent.getByLabel(m_Config.getParameter<edm::InputTag>("srcDTTh"), dtThDigis);
   }
   if (!m_Config.getParameter<bool>("dropRPCPrimitives"))
     iEvent.getByLabel(m_Config.getParameter<edm::InputTag>("srcRPC"), rpcDigis);
   if (!m_Config.getParameter<bool>("dropCSCPrimitives"))
     iEvent.getByLabel(m_Config.getParameter<edm::InputTag>("srcCSC"), cscDigis);
 }

References cscDigis, dtPhDigis, dtThDigis, edm::ParameterSet::getParameter(), iEvent, m_Config, and rpcDigis.

Referenced by reconstruct().

◆ reconstruct()

std::unique_ptr< l1t::RegionalMuonCandBxCollection > OMTFReconstruction::reconstruct	(	const edm::Event &	iEvent,
		const edm::EventSetup &
	)

The order is important: first put omtf_pos candidates, then omtf_neg.

Definition at line 95 of file OMTFReconstruction.cc.

                                                                                                          {
   loadAndFilterDigis(iEvent);
  
   if (dumpResultToXML)
     aTopElement = m_Writer->writeEventHeader(iEvent.id().event());
  
   std::unique_ptr<l1t::RegionalMuonCandBxCollection> candidates(new l1t::RegionalMuonCandBxCollection);
   candidates->setBXRange(bxMin, bxMax);
  
   for (int bx = bxMin; bx <= bxMax; bx++) {
     for (unsigned int iProcessor = 0; iProcessor < m_OMTFConfig->nProcessors(); ++iProcessor)
       getProcessorCandidates(iProcessor, l1t::tftype::omtf_pos, bx, *candidates);
  
     for (unsigned int iProcessor = 0; iProcessor < m_OMTFConfig->nProcessors(); ++iProcessor)
       getProcessorCandidates(iProcessor, l1t::tftype::omtf_neg, bx, *candidates);
  
     edm::LogInfo("OMTFReconstruction") << "OMTF:  Number of candidates in BX=" << bx << ": " << candidates->size(bx)
                                        << std::endl;
     ;
   }
  
   return candidates;
 }

References aTopElement, l1GtPatternGenerator_cfi::bx, bxMax, bxMin, HLT_FULL_cff::candidates, dumpResultToXML, getProcessorCandidates(), iEvent, loadAndFilterDigis(), m_OMTFConfig, m_Writer, OMTFConfiguration::nProcessors(), l1t::omtf_neg, l1t::omtf_pos, and XMLConfigWriter::writeEventHeader().

Referenced by L1TMuonOverlapTrackProducer::produce().

◆ writeResultToXML()

void OMTFReconstruction::writeResultToXML	(	unsigned int	iProcessor,
		l1t::tftype	mtfType,
		const OMTFinput &	myInput,
		const std::vector< OMTFProcessor::resultsMap > &	myResults,
		const std::vector< l1t::RegionalMuonCand > &	candMuons
	)

private

Dump only regions, where a candidate was found

Definition at line 170 of file OMTFReconstruction.cc.

                                                                                            {
   int endcap = (mtfType == l1t::omtf_neg) ? -1 : ((mtfType == l1t::omtf_pos) ? +1 : 0);
   OmtfName board(iProcessor, endcap);
  
   //Write data to XML file
   if (dumpResultToXML) {
     xercesc::DOMElement* aProcElement = m_Writer->writeEventData(aTopElement, board, input);
     for (unsigned int iRefHit = 0; iRefHit < m_OMTFConfig->nTestRefHits(); ++iRefHit) {
       AlgoMuon algoMuon = m_Sorter.sortRefHitResults(results[iRefHit], 0);  //charge=0 means ignore charge
       if (algoMuon.getPt()) {
         m_Writer->writeAlgoMuon(aProcElement, iRefHit, algoMuon);
         if (dumpDetailedResultToXML) {
           for (auto& itKey : results[iRefHit])
             m_Writer->writeResultsData(aProcElement, iRefHit, itKey.first, itKey.second);
         }
       }
     }
     for (auto& candMuon : candMuons)
       m_Writer->writeCandMuon(aProcElement, candMuon);
   }
 }