OMTFReconstruction.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/L1TMuon/interface/RegionalMuonCand.h"
#include "DataFormats/L1TMuon/interface/RegionalMuonCandFwd.h"

#include "CondFormats/DataRecord/interface/L1TMuonOverlapParamsRcd.h"
#include "CondFormats/L1TObjects/interface/L1TMuonOverlapParams.h"

#include "L1Trigger/L1TMuonOverlap/interface/OMTFProcessor.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFinput.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFReconstruction.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFConfiguration.h"
#include "L1Trigger/L1TMuonOverlap/interface/XMLConfigWriter.h"
#include "L1Trigger/L1TMuonOverlap/interface/OmtfName.h"
#include "L1Trigger/L1TMuonOverlap/interface/GhostBusterPreferRefDt.h"

#include "L1Trigger/RPCTrigger/interface/RPCConst.h"

OMTFReconstruction::OMTFReconstruction() :
  m_OMTFConfig(nullptr), m_OMTF(nullptr), aTopElement(nullptr), m_OMTFConfigMaker(nullptr), m_Writer(nullptr){}
OMTFReconstruction::OMTFReconstruction(const edm::ParameterSet& theConfig) :
  m_Config(theConfig), m_OMTFConfig(nullptr), 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;
}
OMTFReconstruction::~OMTFReconstruction(){
  
  delete m_OMTFConfig;
  delete m_OMTF;  

  if (m_Writer) delete m_Writer;
}

void OMTFReconstruction::beginJob() {
  
    m_OMTFConfig = new OMTFConfiguration();
    m_OMTF = new OMTFProcessor();

}
void OMTFReconstruction::endJob(){

  if(dumpResultToXML){
    std::string fName = m_Config.getParameter<std::string>("XMLDumpFileName");
    m_Writer->finaliseXMLDocument(fName);
  } 
}
void OMTFReconstruction::beginRun(edm::Run const& run, edm::EventSetup const& iSetup) {

  const L1TMuonOverlapParamsRcd& omtfRcd = iSetup.get<L1TMuonOverlapParamsRcd>();
  
  edm::ESHandle<L1TMuonOverlapParams> omtfParamsHandle;
  omtfRcd.get(omtfParamsHandle);

  const L1TMuonOverlapParams* omtfParams = omtfParamsHandle.product();

  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.reset(new GhostBusterPreferRefDt(m_OMTFConfig) );
  }
  else {
      m_GhostBuster.reset(new 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);
  }
}
std::unique_ptr<l1t::RegionalMuonCandBxCollection> OMTFReconstruction::reconstruct(const edm::Event& iEvent, const edm::EventSetup& evSetup) {

  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;
}

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

  // 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);
  
}

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

  
  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);
}
void OMTFReconstruction::writeResultToXML(unsigned int iProcessor,  l1t::tftype mtfType,  const OMTFinput &input, 
               const std::vector<OMTFProcessor::resultsMap> & results,
               const std::vector<l1t::RegionalMuonCand> & candMuons ){

  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);
  }
}