l1t::MicroGMTIsolationUnit Class Reference

#include <MicroGMTIsolationUnit.h>

Public Member Functions
void	calculate5by1Sums (const MicroGMTConfiguration::CaloInputCollection &, int bx)
int	calculate5by5Sum (unsigned index) const
void	extrapolateMuons (MicroGMTConfiguration::InterMuonList &) const
int	getCaloIndex (MicroGMTConfiguration::InterMuon &) const
void	isolate (MicroGMTConfiguration::InterMuonList &) const
void	isolatePreSummed (MicroGMTConfiguration::InterMuonList &muons) const
	MicroGMTIsolationUnit (const edm::ParameterSet &)
void	setTowerSums (const MicroGMTConfiguration::CaloInputCollection &inputs, int bx)
virtual	~MicroGMTIsolationUnit ()

Private Attributes
std::vector< int >	m_5by1TowerSums
MicroGMTAbsoluteIsolationCheckLUT	m_AbsIsoCheckMem
MicroGMTExtrapolationLUT	m_BEtaExtrapolation
MicroGMTExtrapolationLUT	m_BPhiExtrapolation
std::map< tftype, MicroGMTExtrapolationLUT * >	m_etaExtrapolationLUTs
MicroGMTExtrapolationLUT	m_FEtaExtrapolation
MicroGMTExtrapolationLUT	m_FPhiExtrapolation
MicroGMTCaloIndexSelectionLUT	m_IdxSelMemEta
MicroGMTCaloIndexSelectionLUT	m_IdxSelMemPhi
bool	m_initialSums
MicroGMTExtrapolationLUT	m_OEtaExtrapolation
MicroGMTExtrapolationLUT	m_OPhiExtrapolation
std::map< tftype, MicroGMTExtrapolationLUT * >	m_phiExtrapolationLUTs
MicroGMTRelativeIsolationCheckLUT	m_RelIsoCheckMem
std::map< int, int >	m_towerEnergies

Detailed Description

Definition at line 12 of file MicroGMTIsolationUnit.h.

Constructor & Destructor Documentation

l1t::MicroGMTIsolationUnit::MicroGMTIsolationUnit ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 9 of file MicroGMTIsolationUnit.cc.

References l1t::bmtf, l1t::emtf_neg, l1t::emtf_pos, m_BEtaExtrapolation, m_BPhiExtrapolation, m_etaExtrapolationLUTs, m_FEtaExtrapolation, m_FPhiExtrapolation, m_OEtaExtrapolation, m_OPhiExtrapolation, m_phiExtrapolationLUTs, l1t::omtf_neg, and l1t::omtf_pos.

                                                                               :
  m_BEtaExtrapolation(iConfig, "BEtaExtrapolationLUTSettings", 0), m_BPhiExtrapolation(iConfig, "BPhiExtrapolationLUTSettings", 1), m_OEtaExtrapolation(iConfig, "OEtaExtrapolationLUTSettings", 0),
  m_OPhiExtrapolation(iConfig, "OPhiExtrapolationLUTSettings", 1), m_FEtaExtrapolation(iConfig, "FEtaExtrapolationLUTSettings", 0), m_FPhiExtrapolation(iConfig, "FPhiExtrapolationLUTSettings", 1),
  m_IdxSelMemEta(iConfig, "IdxSelMemEtaLUTSettings", 0), m_IdxSelMemPhi(iConfig, "IdxSelMemPhiLUTSettings", 1), m_RelIsoCheckMem(iConfig, "RelIsoCheckMemLUTSettings"),
  m_AbsIsoCheckMem(iConfig, "AbsIsoCheckMemLUTSettings"), m_initialSums(false)
{
  m_etaExtrapolationLUTs[tftype::bmtf] = &m_BEtaExtrapolation;
  m_phiExtrapolationLUTs[tftype::bmtf] = &m_BPhiExtrapolation;
  m_etaExtrapolationLUTs[tftype::omtf_pos] = &m_OEtaExtrapolation;
  m_etaExtrapolationLUTs[tftype::omtf_neg] = &m_OEtaExtrapolation;
  m_phiExtrapolationLUTs[tftype::omtf_pos] = &m_OPhiExtrapolation;
  m_phiExtrapolationLUTs[tftype::omtf_neg] = &m_OPhiExtrapolation;
  m_etaExtrapolationLUTs[tftype::emtf_pos] = &m_FEtaExtrapolation;
  m_etaExtrapolationLUTs[tftype::emtf_neg] = &m_FEtaExtrapolation;
  m_phiExtrapolationLUTs[tftype::emtf_pos] = &m_FPhiExtrapolation;
  m_phiExtrapolationLUTs[tftype::emtf_neg] = &m_FPhiExtrapolation;
}

l1t::MicroGMTIsolationUnit::~MicroGMTIsolationUnit ( )

virtual

Definition at line 27 of file MicroGMTIsolationUnit.cc.

{
}

Member Function Documentation

void l1t::MicroGMTIsolationUnit::calculate5by1Sums	(	const MicroGMTConfiguration::CaloInputCollection &	inputs,
		int	bx
	)

Definition at line 75 of file MicroGMTIsolationUnit.cc.

References BXVector< T >::at(), and BXVector< T >::size().

{
  m_5by1TowerSums.clear();
  if (inputs.size(bx) == 0) return;

  for (int iphi = 0; iphi < 36; ++iphi) {
    int iphiIndexOffset = iphi*28;
    m_5by1TowerSums.push_back(inputs.at(bx, iphiIndexOffset).etBits()+inputs.at(bx, iphiIndexOffset+1).etBits()+inputs.at(bx, iphiIndexOffset+2).etBits());//ieta = 0 (tower -28)
    m_5by1TowerSums.push_back(inputs.at(bx, iphiIndexOffset-1).etBits()+inputs.at(bx, iphiIndexOffset).etBits()+inputs.at(bx, iphiIndexOffset+1).etBits()+inputs.at(bx, iphiIndexOffset+2).etBits()); //
    for (int ieta = 2; ieta < 26; ++ieta) {
      int sum = 0;
      for (int dIEta = -2; dIEta <= 2; ++dIEta) {
        sum += inputs.at(bx, iphiIndexOffset+dIEta).etBits();
      }
      m_5by1TowerSums.push_back(sum);
    }
    m_5by1TowerSums.push_back(inputs.at(bx, iphiIndexOffset+1).etBits()+inputs.at(bx, iphiIndexOffset).etBits()+inputs.at(bx, iphiIndexOffset-1).etBits()+inputs.at(bx, iphiIndexOffset-2).etBits());
    m_5by1TowerSums.push_back(inputs.at(bx, iphiIndexOffset).etBits()+inputs.at(bx, iphiIndexOffset-1).etBits()+inputs.at(bx, iphiIndexOffset-2).etBits());//ieta = 0 (tower 28)
  }

  m_initialSums = true;
}

int l1t::MicroGMTIsolationUnit::calculate5by5Sum

(

unsigned

index

)

const

Definition at line 100 of file MicroGMTIsolationUnit.cc.

References min().

{
  if (index > m_5by1TowerSums.size()) {
    edm::LogWarning("energysum out of bounds!");
    return 0;
  }
  // phi wrap around:
  int returnSum = 0;
  for (int dIPhi = -2; dIPhi <= 2; ++dIPhi) {
    int currIndex = (index + dIPhi*28)%1008; // wrap-around at top
    if (currIndex < 0) currIndex = 1008+currIndex;
    if ((unsigned)currIndex < m_5by1TowerSums.size()) {
      returnSum += m_5by1TowerSums[currIndex];
    } else {
      edm::LogWarning("energysum out of bounds!");
    }
  }
  return std::min(31, returnSum);
}

void l1t::MicroGMTIsolationUnit::extrapolateMuons

(

MicroGMTConfiguration::InterMuonList &

inputmuons

)

const

Definition at line 51 of file MicroGMTIsolationUnit.cc.

References funct::abs(), HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, RPCpg::mu, and jetcorrextractor::sign().

Referenced by L1TMicroGMTProducer::produce().

                                                                                               {
  for (auto &mu : inputmuons) {
    // only use 6 LSBs of pt:
    int ptRed = mu->hwPt() & 0b111111;
    // here we drop the two LSBs and masking the MSB
    int etaAbsRed = (std::abs(mu->hwEta()) >> 2) & ((1 << 6) - 1);

    int deltaPhi = 0;
    int deltaEta = 0;

    if (mu->hwPt() < 64) { // extrapolation only for "low" pT muons
      int sign = 1;
      if (mu->hwSign() == 0) {
        sign = -1;
      }
      deltaPhi = (m_phiExtrapolationLUTs.at(mu->trackFinderType())->lookup(etaAbsRed, ptRed) << 3) * sign;
      deltaEta = (m_etaExtrapolationLUTs.at(mu->trackFinderType())->lookup(etaAbsRed, ptRed) << 3);
    }

    mu->setExtrapolation(deltaEta, deltaPhi);
  }
}

int l1t::MicroGMTIsolationUnit::getCaloIndex

(

MicroGMTConfiguration::InterMuon &

mu

)

const

Definition at line 32 of file MicroGMTIsolationUnit.cc.

References eta, l1t::MicroGMTConfiguration::getTwosComp(), l1t::GMTInternalMuon::hwDEta(), l1t::GMTInternalMuon::hwDPhi(), l1t::GMTInternalMuon::hwEta(), l1t::GMTInternalMuon::hwGlobalPhi(), phi(), l1t::GMTInternalMuon::setHwCaloEta(), and l1t::GMTInternalMuon::setHwCaloPhi().

{
  // handle the wrap-around of phi:
  int phi = (mu.hwGlobalPhi() + mu.hwDPhi())%576;
  if (phi < 0) {
    phi = 576+phi;
  }

  int phiIndex = m_IdxSelMemPhi.lookup(phi);
  int eta = mu.hwEta()+mu.hwDEta();
  eta = MicroGMTConfiguration::getTwosComp(eta, 9);
  int etaIndex = m_IdxSelMemEta.lookup(eta);
  mu.setHwCaloEta(etaIndex);
  mu.setHwCaloPhi(phiIndex);

  return phiIndex + etaIndex*36;
}

void l1t::MicroGMTIsolationUnit::isolate

(

MicroGMTConfiguration::InterMuonList &

muons

)

const

Definition at line 121 of file MicroGMTIsolationUnit.cc.

References CastorDataFrameFilter_impl::energySum(), RPCpg::mu, and objects.IsoTrackAnalyzer::relIso.

{
  for (auto& mu : muons) {
    int caloIndex = getCaloIndex(*mu);
    int energySum = calculate5by5Sum(caloIndex);
    mu->setHwIsoSum(energySum);

    int absIso = m_AbsIsoCheckMem.lookup(energySum);
    int relIso = m_RelIsoCheckMem.lookup(energySum, mu->hwPt());

    mu->setHwRelIso(relIso);
    mu->setHwAbsIso(absIso);
  }
}

void l1t::MicroGMTIsolationUnit::isolatePreSummed

(

MicroGMTConfiguration::InterMuonList &

muons

)

const

Definition at line 149 of file MicroGMTIsolationUnit.cc.

References CastorDataFrameFilter_impl::energySum(), RPCpg::mu, and objects.IsoTrackAnalyzer::relIso.

Referenced by L1TMicroGMTProducer::produce().

{
  for (auto mu : muons) {
    int caloIndex = getCaloIndex(*mu);
    int energySum = 0;
    if (m_towerEnergies.count(caloIndex) == 1) {
      energySum = m_towerEnergies.at(caloIndex);
    }

    mu->setHwIsoSum(energySum);

    int absIso = m_AbsIsoCheckMem.lookup(energySum);
    int relIso = m_RelIsoCheckMem.lookup(energySum, mu->hwPt());

    mu->setHwRelIso(relIso);
    mu->setHwAbsIso(absIso);
  }

}

void l1t::MicroGMTIsolationUnit::setTowerSums	(	const MicroGMTConfiguration::CaloInputCollection &	inputs,
		int	bx
	)

Definition at line 136 of file MicroGMTIsolationUnit.cc.

References BXVector< T >::begin(), BXVector< T >::end(), input, and BXVector< T >::size().

Referenced by L1TMicroGMTProducer::produce().

                                                                                                          {
  m_towerEnergies.clear();
  if (inputs.size(bx) == 0) return;
  for (auto input = inputs.begin(bx); input != inputs.end(bx); ++input) {
    if ( input->etBits() != 0 ) {
      m_towerEnergies[input->hwEta()*36+input->hwPhi()] = input->etBits();
    }
  }

  m_initialSums = true;

}

Member Data Documentation

std::vector<int> l1t::MicroGMTIsolationUnit::m_5by1TowerSums

private

Definition at line 50 of file MicroGMTIsolationUnit.h.

MicroGMTAbsoluteIsolationCheckLUT l1t::MicroGMTIsolationUnit::m_AbsIsoCheckMem

private

Definition at line 48 of file MicroGMTIsolationUnit.h.

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_BEtaExtrapolation

private

Definition at line 34 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_BPhiExtrapolation

private

Definition at line 35 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

std::map<tftype, MicroGMTExtrapolationLUT*> l1t::MicroGMTIsolationUnit::m_etaExtrapolationLUTs

private

Definition at line 42 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_FEtaExtrapolation

private

Definition at line 38 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_FPhiExtrapolation

private

Definition at line 39 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTCaloIndexSelectionLUT l1t::MicroGMTIsolationUnit::m_IdxSelMemEta

private

Definition at line 44 of file MicroGMTIsolationUnit.h.

MicroGMTCaloIndexSelectionLUT l1t::MicroGMTIsolationUnit::m_IdxSelMemPhi

private

Definition at line 45 of file MicroGMTIsolationUnit.h.

bool l1t::MicroGMTIsolationUnit::m_initialSums

private

Definition at line 52 of file MicroGMTIsolationUnit.h.

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_OEtaExtrapolation

private

Definition at line 36 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTExtrapolationLUT l1t::MicroGMTIsolationUnit::m_OPhiExtrapolation

private

Definition at line 37 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

std::map<tftype, MicroGMTExtrapolationLUT*> l1t::MicroGMTIsolationUnit::m_phiExtrapolationLUTs

private

Definition at line 41 of file MicroGMTIsolationUnit.h.

Referenced by MicroGMTIsolationUnit().

MicroGMTRelativeIsolationCheckLUT l1t::MicroGMTIsolationUnit::m_RelIsoCheckMem

private

Definition at line 47 of file MicroGMTIsolationUnit.h.

std::map<int, int> l1t::MicroGMTIsolationUnit::m_towerEnergies

private

Definition at line 51 of file MicroGMTIsolationUnit.h.