L1MuGMTPSB.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: L1MuGMTPSB
//
//   Description: Pipelined Synchronising Buffer module
//
//
//
//   Author :
//   N. Neumeister            CERN EP
//
//   Migrated to CMSSW:
//   I. Mikulec
//
//--------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------

#include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTPSB.h"

//---------------
// C++ Headers --
//---------------

#include <iostream>
#include <iomanip>
#include <vector>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------

#include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTConfig.h"
#include "CondFormats/L1TObjects/interface/L1MuTriggerScales.h"
#include "CondFormats/L1TObjects/interface/L1MuTriggerPtScale.h"
#include "CondFormats/L1TObjects/interface/L1MuGMTChannelMask.h"

#include "L1Trigger/GlobalMuonTrigger/interface/L1MuGlobalMuonTrigger.h"
#include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"

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

// --------------------------------
//       class L1MuGMTPSB
//---------------------------------

//----------------
// Constructors --
//----------------
L1MuGMTPSB::L1MuGMTPSB(const L1MuGlobalMuonTrigger& gmt, edm::ConsumesCollector&& iC)
    : m_gmt(gmt),
      m_RpcMuons(L1MuGMTConfig::MAXRPC),
      m_DtbxMuons(L1MuGMTConfig::MAXDTBX),
      m_CscMuons(L1MuGMTConfig::MAXCSC),
      m_Isol(14, 18, false),
      m_Mip(14, 18, false) {
  m_RpcMuons.reserve(L1MuGMTConfig::MAXRPC);
  m_DtbxMuons.reserve(L1MuGMTConfig::MAXDTBX);
  m_CscMuons.reserve(L1MuGMTConfig::MAXCSC);
  iC.consumes<std::vector<L1MuRegionalCand> >(L1MuGMTConfig::getDTInputTag());
  iC.consumes<std::vector<L1MuRegionalCand> >(L1MuGMTConfig::getCSCInputTag());
  iC.consumes<std::vector<L1MuRegionalCand> >(L1MuGMTConfig::getRPCbInputTag());
  iC.consumes<std::vector<L1MuRegionalCand> >(L1MuGMTConfig::getRPCfInputTag());
  iC.consumes<L1CaloRegionCollection>(L1MuGMTConfig::getMipIsoInputTag());
}

//--------------
// Destructor --
//--------------
L1MuGMTPSB::~L1MuGMTPSB() {
  reset();
  m_RpcMuons.clear();
  m_DtbxMuons.clear();
  m_CscMuons.clear();
}

//--------------
// Operations --
//--------------

//
// receive data
//
void L1MuGMTPSB::receiveData(edm::Event& e, int bx) {

  edm::Handle<std::vector<L1MuRegionalCand> > rc_handle;

  const L1MuGMTChannelMask* theChannelMask = L1MuGMTConfig::getGMTChanMask();
  unsigned mask = theChannelMask->getSubsystemMask();

  if ((L1MuGMTConfig::getDTInputTag()).label() != "none" && !(mask & 1)) {
    e.getByLabel(L1MuGMTConfig::getDTInputTag(), rc_handle);
    if (rc_handle.isValid()) {
      getDTBX(rc_handle.product(), bx);
    } else {
      if (L1MuGMTConfig::Debug(1)) {
        edm::LogWarning("GlobalMuonTrigger")
            << "\nWarning: GlobalMuonTrigger: input tag " << L1MuGMTConfig::getDTInputTag()
            << "\nrequested, but not found in the event." << std::endl;
      }
    }
  }
  if ((L1MuGMTConfig::getCSCInputTag()).label() != "none" && !(mask & 4)) {
    e.getByLabel(L1MuGMTConfig::getCSCInputTag(), rc_handle);
    if (rc_handle.isValid()) {
      getCSC(rc_handle.product(), bx);
    } else {
      if (L1MuGMTConfig::Debug(1)) {
        edm::LogWarning("GlobalMuonTrigger")
            << "\nWarning: GlobalMuonTrigger: input tag " << L1MuGMTConfig::getCSCInputTag()
            << "\nrequested, but not found in the event." << std::endl;
      }
    }
  }
  if ((L1MuGMTConfig::getRPCbInputTag()).label() != "none" && !(mask & 2)) {
    e.getByLabel(L1MuGMTConfig::getRPCbInputTag(), rc_handle);
    if (rc_handle.isValid()) {
      getRPCb(rc_handle.product(), bx);
    } else {
      if (L1MuGMTConfig::Debug(1)) {
        edm::LogWarning("GlobalMuonTrigger")
            << "\nWarning: GlobalMuonTrigger: input tag " << L1MuGMTConfig::getRPCbInputTag()
            << "\nrequested, but not found in the event." << std::endl;
      }
    }
  }
  if ((L1MuGMTConfig::getRPCfInputTag()).label() != "none" && !(mask & 8)) {
    e.getByLabel(L1MuGMTConfig::getRPCfInputTag(), rc_handle);
    if (rc_handle.isValid()) {
      getRPCf(rc_handle.product(), bx);
    } else {
      if (L1MuGMTConfig::Debug(1)) {
        edm::LogWarning("GlobalMuonTrigger")
            << "\nWarning: GlobalMuonTrigger: input tag " << L1MuGMTConfig::getRPCfInputTag()
            << "\nrequested, but not found in the event." << std::endl;
      }
    }
  }


  const L1MuTriggerScales* theTriggerScales = L1MuGMTConfig::getTriggerScales();
  const L1MuTriggerPtScale* theTriggerPtScale = L1MuGMTConfig::getTriggerPtScale();

  // store data in readout record
  for (int i = 0; i < 4; i++) {
    L1MuRegionalCand* cand = &(m_DtbxMuons[i]);
    cand->setPhiValue(theTriggerScales->getPhiScale()->getLowEdge(cand->phi_packed()));
    cand->setEtaValue(theTriggerScales->getRegionalEtaScale(cand->type_idx())->getCenter(cand->eta_packed()));
    cand->setPtValue(theTriggerPtScale->getPtScale()->getLowEdge(cand->pt_packed()));
    // cand->setPtValue( theTriggerScales->getPtScale()->getLowEdge(cand->pt_packed()) );
    m_gmt.currentReadoutRecord()->setInputCand(i, *cand);
  }
  for (int i = 0; i < 4; i++) {
    L1MuRegionalCand* cand = &(m_RpcMuons[i]);
    cand->setPhiValue(theTriggerScales->getPhiScale()->getLowEdge(cand->phi_packed()));
    cand->setEtaValue(theTriggerScales->getRegionalEtaScale(cand->type_idx())->getCenter(cand->eta_packed()));
    cand->setPtValue(theTriggerPtScale->getPtScale()->getLowEdge(cand->pt_packed()));
    // cand->setPtValue( theTriggerScales->getPtScale()->getLowEdge(cand->pt_packed()) );
    m_gmt.currentReadoutRecord()->setInputCand(i + 4, *cand);
  }
  for (int i = 0; i < 4; i++) {
    L1MuRegionalCand* cand = &(m_CscMuons[i]);
    cand->setPhiValue(theTriggerScales->getPhiScale()->getLowEdge(cand->phi_packed()));
    cand->setEtaValue(theTriggerScales->getRegionalEtaScale(cand->type_idx())->getCenter(cand->eta_packed()));
    cand->setPtValue(theTriggerPtScale->getPtScale()->getLowEdge(cand->pt_packed()));
    // cand->setPtValue( theTriggerScales->getPtScale()->getLowEdge(cand->pt_packed()) );
    m_gmt.currentReadoutRecord()->setInputCand(i + 8, *cand);
  }
  for (int i = 0; i < 4; i++) {
    L1MuRegionalCand* cand = &(m_RpcMuons[i + 4]);
    cand->setPhiValue(theTriggerScales->getPhiScale()->getLowEdge(cand->phi_packed()));
    cand->setEtaValue(theTriggerScales->getRegionalEtaScale(cand->type_idx())->getCenter(cand->eta_packed()));
    cand->setPtValue(theTriggerPtScale->getPtScale()->getLowEdge(cand->pt_packed()));
    // cand->setPtValue( theTriggerScales->getPtScale()->getLowEdge(cand->pt_packed()) );
    m_gmt.currentReadoutRecord()->setInputCand(i + 12, *cand);
  }

  // if there is at least one muon start the calorimeter trigger

  if (L1MuGMTConfig::getCaloTrigger() && !empty())
    getCalo(e);
}

//
// clear PSB
//
void L1MuGMTPSB::reset() {
  std::vector<L1MuRegionalCand>::iterator iter;
  iter = m_RpcMuons.begin();
  while (iter != m_RpcMuons.end())
    (*(iter++)).reset();

  iter = m_DtbxMuons.begin();
  while (iter != m_DtbxMuons.end())
    (*(iter++)).reset();

  iter = m_CscMuons.begin();
  while (iter != m_CscMuons.end())
    (*(iter++)).reset();

  m_Isol.reset(false);
  m_Mip.reset(false);
}

//
// print muons
//
void L1MuGMTPSB::print() const {
  edm::LogVerbatim("GMT_PSB_info") << " ";
  printDTBX();
  printRPCbarrel();
  printCSC();
  printRPCendcap();
  edm::LogVerbatim("GMT_PSB_info") << " ";
}

//
// return RPC muon
//
const L1MuRegionalCand* L1MuGMTPSB::RPCMuon(int index) const {
  return (index < (int)L1MuGMTConfig::MAXRPC && index >= 0) ? &(m_RpcMuons[index]) : nullptr;
}

//
// return DTBX muon
//
const L1MuRegionalCand* L1MuGMTPSB::DTBXMuon(int index) const {
  return (index < (int)L1MuGMTConfig::MAXDTBX && index >= 0) ? &(m_DtbxMuons[index]) : nullptr;
}

//
// return CSC muon
//
const L1MuRegionalCand* L1MuGMTPSB::CSCMuon(int index) const {
  return (index < (int)L1MuGMTConfig::MAXCSC && index >= 0) ? &(m_CscMuons[index]) : nullptr;
}

//
// count number of non empty RPC muons
//
int L1MuGMTPSB::numberRPC() const {
  int count = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter = m_RpcMuons.begin();
  while (iter != m_RpcMuons.end()) {
    if (!(*iter).empty())
      count++;
    iter++;
  }
  return count;
}

//
// count number of non empty DT muons
//
int L1MuGMTPSB::numberDTBX() const {
  int count = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter = m_DtbxMuons.begin();
  while (iter != m_DtbxMuons.end()) {
    if (!(*iter).empty())
      count++;
    iter++;
  }
  return count;
}

//
// count number of non empty CSC muons
//
int L1MuGMTPSB::numberCSC() const {
  int count = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter = m_CscMuons.begin();
  while (iter != m_CscMuons.end()) {
    if (!(*iter).empty())
      count++;
    iter++;
  }
  return count;
}

//
// are there any data in the PSB
//
bool L1MuGMTPSB::empty() const {
  int number = numberRPC() + numberDTBX() + numberCSC();

  return (number == 0);
}

//
// get muons from RPCb Trigger
//
void L1MuGMTPSB::getRPCb(std::vector<L1MuRegionalCand> const* data, int bx) {
  int irpcb = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter;
  for (iter = data->begin(); iter != data->end(); iter++) {
    if ((*iter).bx() != bx)
      continue;
    if (irpcb < (int)L1MuGMTConfig::MAXRPCbarrel) {
      if (!(*iter).empty())
        m_RpcMuons[irpcb] = (*iter);
      irpcb++;
    }
  }
}

//
// get muons from RPCf Trigger
//
void L1MuGMTPSB::getRPCf(std::vector<L1MuRegionalCand> const* data, int bx) {
  int irpcf = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter;
  for (iter = data->begin(); iter != data->end(); iter++) {
    if ((*iter).bx() != bx)
      continue;
    if (irpcf < (int)L1MuGMTConfig::MAXRPCendcap) {
      if (!(*iter).empty())
        m_RpcMuons[irpcf + 4] = (*iter);
      irpcf++;
    }
  }
}

//
// get muons from barrel Muon Trigger Track Finder
//
void L1MuGMTPSB::getDTBX(std::vector<L1MuRegionalCand> const* data, int bx) {
  // temporary hack with bxoffset - to be removed, trigger bx should be 0
  int bxoffset = 0;
  int idtbx = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter;
  for (iter = data->begin(); iter != data->end(); iter++) {
    if (L1MuGMTConfig::Debug(2))
      edm::LogVerbatim("") << "DTTF BX: " << (*iter).bx() << " my bx: " << bx;
    if ((*iter).bx() > 10)
      bxoffset = 16;
    if ((*iter).bx() != bx + bxoffset)
      continue;
    if (idtbx < (int)L1MuGMTConfig::MAXDTBX) {
      m_DtbxMuons[idtbx] = (*iter);
      m_DtbxMuons[idtbx].setBx(bx);
      idtbx++;
    }
  }
}

//
// get muons from CSC Track Finder
//
void L1MuGMTPSB::getCSC(std::vector<L1MuRegionalCand> const* data, int bx) {
  int icsc = 0;
  std::vector<L1MuRegionalCand>::const_iterator iter;
  for (iter = data->begin(); iter != data->end(); iter++) {
    if ((*iter).bx() != bx)
      continue;
    if (icsc < (int)L1MuGMTConfig::MAXCSC) {
      m_CscMuons[icsc] = (*iter);
      icsc++;
    }
  }
}

//
// print barrel RPC muons
//
void L1MuGMTPSB::printRPCbarrel() const {
  edm::LogVerbatim("GMT_PSB_info") << "RPC barrel  muons received by the GMT :";

  for (unsigned i = 0; i < L1MuGMTConfig::MAXRPCbarrel; i++) {
    if (!m_RpcMuons[i].empty())
      m_RpcMuons[i].print();
  }
}

//
// print endcap RPC muons
//
void L1MuGMTPSB::printRPCendcap() const {
  edm::LogVerbatim("GMT_PSB_info") << "RPC endcap  muons received by the GMT :";

  for (unsigned i = 0; i < L1MuGMTConfig::MAXRPCendcap; i++) {
    if (!m_RpcMuons[i + 4].empty())
      m_RpcMuons[i + 4].print();
  }
}

//
// print DTBX muons
//
void L1MuGMTPSB::printDTBX() const {
  edm::LogVerbatim("GMT_PSB_info") << "DTBX muons received by the GMT :";

  for (unsigned i = 0; i < L1MuGMTConfig::MAXDTBX; i++) {
    if (!m_DtbxMuons[i].empty())
      m_DtbxMuons[i].print();
  }
}

//
// print CSC muons
//
void L1MuGMTPSB::printCSC() const {
  edm::LogVerbatim("GMT_PSB_info") << "CSC  muons received by the GMT :";

  for (unsigned i = 0; i < L1MuGMTConfig::MAXCSC; i++) {
    if (!m_CscMuons[i].empty())
      m_CscMuons[i].print();
  }
}

//
// get data from regional calorimeter trigger
//
void L1MuGMTPSB::getCalo(edm::Event& e) {
  edm::Handle<L1CaloRegionCollection> calocoll_h;
  e.getByLabel(L1MuGMTConfig::getMipIsoInputTag(), calocoll_h);
  if (calocoll_h.isValid()) {
    L1CaloRegionCollection const* regions = calocoll_h.product();
    L1CaloRegionCollection::const_iterator iter;

    //  edm::LogVerbatim("GMT_PSB_info") << "MIP/QUIET bits rceived by the GMT :";

    for (iter = regions->begin(); iter != regions->end(); iter++) {
      if ((*iter).id().ieta() < 4 || (*iter).id().ieta() > 17 || (*iter).id().iphi() > 17)
        continue;
      m_Isol.set((*iter).id().ieta() - 4, (*iter).id().iphi(), (*iter).quiet());
      m_Mip.set((*iter).id().ieta() - 4, (*iter).id().iphi(), (*iter).mip());

      if ((*iter).quiet())
        m_gmt.currentReadoutRecord()->setQuietbit((*iter).id().ieta() - 4, (*iter).id().iphi());

      if ((*iter).mip())
        m_gmt.currentReadoutRecord()->setMIPbit((*iter).id().ieta() - 4, (*iter).id().iphi());

      //    edm::LogVerbatim("GMT_PSB_info") << (*iter).id().ieta()-4 << " "
      //                                     << (*iter).id().iphi() << " "
      //                                     << (*iter).quiet() << " "
      //                                     << (*iter).mip();
    }
  } else {
    if (L1MuGMTConfig::Debug(1)) {
      edm::LogWarning("GlobalMuonTrigger")
          << "\nWarning: GlobalMuonTrigger: input tag " << L1MuGMTConfig::getMipIsoInputTag()
          << "\nrequested, but not found in the event." << std::endl;
    }
  }
}