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), m_Mip(14,18) {

  m_RpcMuons.reserve(L1MuGMTConfig::MAXRPC);
  m_DtbxMuons.reserve(L1MuGMTConfig::MAXDTBX);
  m_CscMuons.reserve(L1MuGMTConfig::MAXCSC);
  m_Isol.init(false);
  m_Mip.init(false);
  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.init(false);
  m_Mip.init(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;
    }
  }

}