LVL1::CPMTobAlgorithm Class Reference

This is an internal class, used in the EM/Tau trigger. More...

#include <CPMTobAlgorithm.h>

Collaboration diagram for LVL1::CPMTobAlgorithm:

Public Member Functions
	CPMTobAlgorithm (double eta, double phi, const xAOD::CPMTowerMap_t *ttContainer, const TrigConf::L1Menu *l1menu, int slice=-1)
	~CPMTobAlgorithm ()
double	eta ()
	Accessors for TOB data.
double	phi ()
	Returns phi coordinate of RoI, using standard ATLAS convention.
int	EMClusET ()
	Returns EM cluster ET, limited to 8 bits.
int	EMIsolWord ()
	Returns EM isolation word (5 bits)
int	TauClusET ()
	Returns Tau cluster ET, limited to 8 bits.
int	TauIsolWord ()
	Returns Tau isolation word (5 bits)
bool	isEMRoI ()
	Report whether this passed as an EM TOB.
bool	isTauRoI ()
	Report whether this passed as a Tau TOB.
bool	isEtMax ()
	Does this window pass the local ET maximum condition.
int	CoreET ()
	Additional information for reconstruction & performance studies.
int	EMCoreET ()
	Returns EM core ET.
int	HadCoreET ()
	Returns Had core ET (inner isolation sum)
int	EMIsolET ()
	Returns EM isolation ET.
int	HadIsolET ()
	Returns Had isolation ET.
int	EMLUTClus ()
	Sums with the range & precision of isolation LUT inputs.
int	EMLUTEMIsol ()
	Returns EM isolation ET with range and precision of EM LUT input.
int	EMLUTHadVeto ()
	Returns Had veto ET with range and precision of EM LUT input.
int	TauLUTClus ()
	Returns Tau cluster ET with range and precision of LUT input.
int	TauLUTEMIsol ()
	Returns EM isolation ET with range and precision of Tau LUT input.
xAOD::CPMTobRoI *	EMCPMTobRoI ()
	Create CPMTobRoIs and return pointers to them.
xAOD::CPMTobRoI *	TauCPMTobRoI ()
	Returns CPMTobRoI for TAU hypothesis, provided TAU Tob conditions passed.

Private Member Functions
void	setRoICoord (double eta, double phi)
	threshold values
void	testEtMax (const std::vector< unsigned int > &cores)
	Form all 2x2 clusters within window and test centre is a local ET maximum.
void	emAlgorithm ()
	Check trigger condition and set hits if appropriate.
void	tauAlgorithm ()
	Check trigger condition and set hits if appropriate.
xAOD::CPMTobRoI *	createTobRoI (int type)
	Create a pointer to a CPMTobRoI and return it.

Private Attributes
double	m_refEta
double	m_refPhi
const TrigConf::L1Menu *	m_l1menu {nullptr}
double	m_eta
	Algorithm results.
double	m_phi
unsigned int	m_Core
unsigned int	m_EMClus
unsigned int	m_TauClus
unsigned int	m_EMIsolWord
unsigned int	m_TauIsolWord
unsigned int	m_EMCore
unsigned int	m_EMIsol
unsigned int	m_HadCore
unsigned int	m_HadIsol
bool	m_EtMax
bool	m_EMThresh
bool	m_TauThresh

Static Private Attributes
static const unsigned int	m_maxClus = 0xFF
	Algorithm parameters.
static const unsigned int	m_emLUT_ClusterFirstBit = 1
static const unsigned int	m_emLUT_ClusterNBits = 6
static const unsigned int	m_emLUT_EMIsolFirstBit = 0
static const unsigned int	m_emLUT_EMIsolNBits = 4
static const unsigned int	m_emLUT_HadVetoFirstBit = 0
static const unsigned int	m_emLUT_HadVetoNBits = 3
static const unsigned int	m_tauLUT_ClusterFirstBit = 1
static const unsigned int	m_tauLUT_ClusterNBits = 7
static const unsigned int	m_tauLUT_EMIsolFirstBit = 0
static const unsigned int	m_tauLUT_EMIsolNBits = 6
static const unsigned int	m_noIsol = 999
static const double	m_maxEta = 2.5

Detailed Description

This is an internal class, used in the EM/Tau trigger.

The CPMTobAlgorithm:

• identifies towers used in a particular trigger window (4x4 towers)
  
• forms sums required for trigger algorithm
  
• identifies whether passes the RoI condition (local ET maximum)
  
• returns sums (with appropriate precision)
  
• returns hit bits and RoI word (if all requirements met)

Definition at line 43 of file CPMTobAlgorithm.h.

Constructor & Destructor Documentation

CPMTobAlgorithm()

LVL1::CPMTobAlgorithm::CPMTobAlgorithm	(	double	eta,
		double	phi,
		const xAOD::CPMTowerMap_t *	ttContainer,
		const TrigConf::L1Menu *	l1menu,
		int	slice = -1 )

RoI coordinates are centre of window, while key classes are designed for TT coordinates - differ by 0.5* TT_size. Using wrong coordinate convention can lead to precision errors in key generation.

Offsetting by TT_size/4 will put coordinates inside reference TT whether input was reference TT coordinate or RoI centre coordinate. To be extra safe, the code below will then compute the centre of the reference tower.

This should all ensure consistent & safe key generation however the initial coordinate was specified.

Now loop over towers within window and form clusters.

Interesting one is the 9 RoI cores. The logic below fills these as an a 9-element std::vector. The ordering of the elements is as follows: 2 5 8 1 4 7 0 3 6 So the RoI "ET maximum" condition is that 4 >= 0-3, 4 > 5-8 If you picture the cores in a 3x3 array, cores[iphi][ieta], then the index in the vector above = iphi + 3*ieta.

Definition at line 41 of file CPMTobAlgorithm.cxx.

                                                                                  :
  m_l1menu(l1menu),
  m_Core(0),
  m_EMClus(0),
  m_TauClus(0),
  m_EMIsolWord(0),
  m_TauIsolWord(0),
  m_EMCore(0),
  m_EMIsol(0),
  m_HadCore(0),
  m_HadIsol(0),
  m_EtMax(false),
  m_EMThresh(false),
  m_TauThresh(false)
  //m_debug(false)
{
 
  TriggerTowerKey get(phi-M_PI/128.,eta-0.025);
  Coordinate refCoord = get.coord();
  m_refEta = refCoord.eta();
  m_refPhi = refCoord.phi();
 
  // Set coordinate of centre of RoI, starting from reference tower coordinate
  setRoICoord(m_refEta, m_refPhi);
  
 
  std::vector<unsigned int> emEt(4);
  std::vector<unsigned int> cores(9);
  for (int etaOffset = -1; etaOffset <= 2; etaOffset++) {
    double tempEta = m_refEta + etaOffset*0.1;
    for (int phiOffset = -1; phiOffset <= 2; phiOffset++) {
      double tempPhi = m_refPhi + phiOffset*M_PI/32;
      int key = get.ttKey(tempPhi, tempEta);
      xAOD::CPMTowerMap_t::const_iterator tt = ttContainer->find(key);
      if (tt != ttContainer->end() && fabs(tempEta) < m_maxEta) {
        // Get the TT ET values once here, rather than repeat function calls
        int emTT = 0;
        int hadTT = 0;
        if (slice < 0) {  // Default to using peak slice
          emTT  = (tt->second)->emEnergy();
          hadTT = (tt->second)->hadEnergy();
        }
        else {            // Use user-specified slice
          emTT  = (tt->second)->emSliceEnergy(slice);
          hadTT = (tt->second)->hadSliceEnergy(slice);
        }
        // EM Trigger Clusters and Core clusters
        if (etaOffset >= 0 && etaOffset <= 1 && phiOffset >= 0 && phiOffset <= 1) {
          emEt[phiOffset+2*etaOffset] = emTT;
          m_EMCore += emTT;
          m_HadCore += hadTT;
        }
        // Isolation sums
        else {
          m_EMIsol  += emTT;
          m_HadIsol += hadTT;
        }
        // Each tower is part of up to 4 RoI core clusters
        if (etaOffset >= 0) {
          if (phiOffset >= 0) cores[phiOffset+3*etaOffset] += emTT + hadTT;
          if (phiOffset < 2)  cores[phiOffset+3*etaOffset+1] += emTT + hadTT;
        }
        if (etaOffset < 2) {
          if (phiOffset >= 0) cores[phiOffset+3*etaOffset+3] += emTT + hadTT;
          if (phiOffset < 2)  cores[phiOffset+3*etaOffset+4] += emTT + hadTT;
        }
 
      }
    } // phi offset loop
  }   // eta offset loop
 
  // Core of this window is cores[4]
  m_Core = cores[4];
 
  // Form most energetic 2-tower EM cluster
  m_EMClus = ( (emEt[0] > emEt[3]) ? emEt[0] : emEt[3] )
           + ( (emEt[1] > emEt[2]) ? emEt[1] : emEt[2] );
 
  // Tau cluster
  m_TauClus = m_EMClus + m_HadCore;
  
  // Check whether RoI condition is met.
  testEtMax(cores);
 
  // test cluster and isolation conditions
  emAlgorithm();
  tauAlgorithm();
  
}

~CPMTobAlgorithm()

LVL1::CPMTobAlgorithm::~CPMTobAlgorithm

(

)

Definition at line 157 of file CPMTobAlgorithm.cxx.

                                     {
}

Member Function Documentation

CoreET()

int LVL1::CPMTobAlgorithm::CoreET

(

)

Additional information for reconstruction & performance studies.

Returns RoI Core ET.

Definition at line 374 of file CPMTobAlgorithm.cxx.

                                {
  return m_Core;
}

createTobRoI()

xAOD::CPMTobRoI * LVL1::CPMTobAlgorithm::createTobRoI ( int type )

private

Create a pointer to a CPMTobRoI and return it.

This function does not check whether a valid RoI exists - that is left to the calling functions (but as this is a private method that should be enough)

Definition at line 494 of file CPMTobAlgorithm.cxx.

                                                         {
  
  xAOD::CPMTobRoI* pointer = 0;
  
  int et = 0;
  int isol = 0;
  
  if (type == TrigT1CaloDefs::emTobType) {
    et   = EMClusET();
    isol = EMIsolWord();
  }
  else if (type == TrigT1CaloDefs::tauTobType) {
    et   = TauClusET();
    isol = TauIsolWord();
  }
  else return pointer;
  
  // Need to calculate hardware coordinate
  // Object coordinate is RoI centre, but this works better if use reference tower coord
  float eta = m_eta - 0.05;
  float phi = m_phi - M_PI/64;
  if (phi < 0) phi += 2*M_PI; 
  Coordinate coord(phi, eta);
  CoordToHardware convertor;
  
  int crate = convertor.cpCrate(coord);
  int module = convertor.cpModule(coord);
  int chip = convertor.cpModuleFPGA(coord);
  int location = convertor.cpModuleLocalRoI(coord);
  
  /*
  float phi = ( (m_phi >= 0) ?  m_phi : 2.*M_PI + m_phi );
  int crate = 2*phi/M_PI;
  int module = (m_eta + 2.8)/0.4 + 1;
  int iPhi = (phi - crate*M_PI/2)*32/M_PI;
  int chip = iPhi/2;
  int iEta = (m_eta + 2.8 - module*0.4)/0.1;
  int iSide = iEta/2;
  int location = (iSide << 2) + ((iPhi&1) << 1) + (iEta&1);
  */
  
  pointer = new xAOD::CPMTobRoI;
  pointer->makePrivateStore();
  pointer->initialize(crate, module, chip, location, type, et, isol);
  
  return pointer;
  
}

emAlgorithm()

void LVL1::CPMTobAlgorithm::emAlgorithm ( )

private

Check trigger condition and set hits if appropriate.

Disable isolation by default, then set cut if defined and in range. Expect parameters to be on 100 MeV scale whatever the digit scale, so need to rescale if digit scale differs. Also must rescale ET sums to same scale for comparison

Definition at line 205 of file CPMTobAlgorithm.cxx.

                                      {
  
  // Initialise to correct defaults, to be sure
  m_EMThresh = false;
  m_EMIsolWord = 0;
  
  // Don't waste time if it isn't an RoI candidate
  if (!m_EtMax) return;
  
  // Does this pass min TOB pT cut?
  // Need to worry about digit scale not being GeV
  float emEnergyScale{1}; 
  unsigned int tobPtMinEM{0};
  emEnergyScale = m_l1menu->thrExtraInfo().EM().emScale();
  tobPtMinEM = m_l1menu->thrExtraInfo().EM().ptMinToTopo();
 
  unsigned int thresh = tobPtMinEM * emEnergyScale;
 
  if (m_EMClus <= thresh) return;
  
  // Passed, so there is a TOB here
  m_EMThresh = true;
 
  // Now calculate isolation word
  // Define cluster and isolation values with appropriate scales and range for isolation LUT  
  unsigned int clusMask    = ( (1<<m_emLUT_ClusterNBits)-1 )<<m_emLUT_ClusterFirstBit;
  unsigned int emIsolMask  = ( (1<<m_emLUT_EMIsolNBits) -1 )<<m_emLUT_EMIsolFirstBit;
  unsigned int hadIsolMask = ( (1<<m_emLUT_HadVetoNBits)-1 )<<m_emLUT_HadVetoFirstBit;
 
  // If cluster overflows, set all threshold bits
  if (m_EMClus >= clusMask) {
    m_EMIsolWord = 0x1F;
    return;
  }
  
  // For consistency with LUT filling and menu parameters, convert ET sums to 100 MeV units
  int clus    = (( (m_EMClus < clusMask)     ? (m_EMClus  & clusMask)    : clusMask ) * 10)/emEnergyScale;
  int emIsol  = (( (m_EMIsol < emIsolMask)   ? (m_EMIsol  & emIsolMask)  : emIsolMask ) * 10)/emEnergyScale;
  int hadVeto = (( (m_HadCore < hadIsolMask) ? (m_HadCore & hadIsolMask) : hadIsolMask ) * 10)/emEnergyScale;

  std::vector<int> emisolcuts(5,m_noIsol);
  std::vector<int> hadisolcuts(5,m_noIsol);
  for(size_t bit = 1; bit<=5; ++bit) {
    const TrigConf::IsolationLegacy & iso = m_l1menu->thrExtraInfo().EM().isolation("EMIsoForEMthr", bit);
    if(! iso.isDefined() ) {
      continue;
    }
    int offset = iso.offset();
    int slope  = iso.slope();
    int mincut = iso.mincut();
    int upperlimit = iso.upperlimit()*10;
    if (clus < upperlimit && bit >= 1 && bit <= 5) {
      // As "slope" = 10* slope, rescale clus too to get fraction right.
      int cut = offset + (slope != 0 ? 10*clus/slope : 0);
      if (cut < mincut) cut = mincut;
      emisolcuts[bit-1] = cut;
    }
  }
  for(size_t bit = 1; bit<=5; ++bit) {
    const TrigConf::IsolationLegacy & iso = m_l1menu->thrExtraInfo().EM().isolation("HAIsoForEMthr", bit);
    if(! iso.isDefined() ) {
      continue;
    }
    int offset = iso.offset();
    int slope  = iso.slope();
    int mincut = iso.mincut();
    int upperlimit = iso.upperlimit()*10;
    if (clus < upperlimit && bit >= 1 && bit <= 5) {
      // As "slope" = 10* slope, rescale clus too to get fraction right.
      int cut = offset + (slope != 0 ? 10*clus/slope : 0);
      if (cut < mincut) cut = mincut;
      hadisolcuts[bit-1] = cut;
    }
  }
 
  // Set isolation bits
  for (unsigned int isol = 0; isol < 5; ++isol) {
    if (emIsol <= emisolcuts[isol] && hadVeto <= hadisolcuts[isol]) m_EMIsolWord += (1 << isol);
  }
    
}

EMClusET()

int LVL1::CPMTobAlgorithm::EMClusET

(

)

Returns EM cluster ET, limited to 8 bits.

Definition at line 379 of file CPMTobAlgorithm.cxx.

                                  {
  return ( (m_EMClus < m_maxClus) ? m_EMClus : m_maxClus );
}

EMCoreET()

int LVL1::CPMTobAlgorithm::EMCoreET

(

)

Returns EM core ET.

Definition at line 409 of file CPMTobAlgorithm.cxx.

                                  {
  return m_EMCore;
}

EMCPMTobRoI()

xAOD::CPMTobRoI * LVL1::CPMTobAlgorithm::EMCPMTobRoI

(

)

Create CPMTobRoIs and return pointers to them.

Returns CPMTobRoI for EM hypothesis, provided EM Tob conditions passed.

Will return a null pointer if object does not pass hypothesis. It is the user's responsibility to check

Definition at line 469 of file CPMTobAlgorithm.cxx.

                                                {
  if (isEMRoI()) {
    xAOD::CPMTobRoI* pointer = createTobRoI(TrigT1CaloDefs::emTobType); 
    return pointer;
  }
  
  return 0;
 
}

EMIsolET()

int LVL1::CPMTobAlgorithm::EMIsolET

(

)

Returns EM isolation ET.

Definition at line 399 of file CPMTobAlgorithm.cxx.

                                  {
  return m_EMIsol;
}

EMIsolWord()

int LVL1::CPMTobAlgorithm::EMIsolWord

(

)

Returns EM isolation word (5 bits)

Definition at line 389 of file CPMTobAlgorithm.cxx.

                                    {
  return m_EMIsolWord;
}

EMLUTClus()

int LVL1::CPMTobAlgorithm::EMLUTClus

(

)

Sums with the range & precision of isolation LUT inputs.

Returns EM cluster ET with range and precision of LUT input.

Definition at line 424 of file CPMTobAlgorithm.cxx.

                                   {
  int ET = (m_EMClus >> m_emLUT_ClusterFirstBit);
  if (ET > (1<< m_emLUT_ClusterNBits) - 1) ET = (1<<m_emLUT_ClusterNBits) - 1;
  return (ET << m_emLUT_ClusterFirstBit);
}

EMLUTEMIsol()

int LVL1::CPMTobAlgorithm::EMLUTEMIsol

(

)

Returns EM isolation ET with range and precision of EM LUT input.

Definition at line 438 of file CPMTobAlgorithm.cxx.

                                     {
  int ET = (m_EMIsol >> m_emLUT_EMIsolFirstBit);
  if (ET > (1<<m_emLUT_EMIsolNBits) - 1) ET = (1<<m_emLUT_EMIsolNBits) - 1;
  return (ET << m_emLUT_EMIsolFirstBit);
}

EMLUTHadVeto()

int LVL1::CPMTobAlgorithm::EMLUTHadVeto

(

)

Returns Had veto ET with range and precision of EM LUT input.

Definition at line 445 of file CPMTobAlgorithm.cxx.

                                      {
  int ET = (m_HadCore >> m_emLUT_HadVetoFirstBit);
  if (ET > (1<<m_emLUT_HadVetoNBits) - 1) ET = (1<<m_emLUT_HadVetoNBits) - 1;
  return (ET << m_emLUT_HadVetoFirstBit);
}

eta()

double LVL1::CPMTobAlgorithm::eta

(

)

Accessors for TOB data.

Returns eta coordinate of RoI.

Definition at line 459 of file CPMTobAlgorithm.cxx.

                                {
  return m_eta;
}

HadCoreET()

int LVL1::CPMTobAlgorithm::HadCoreET

(

)

Returns Had core ET (inner isolation sum)

Definition at line 414 of file CPMTobAlgorithm.cxx.

                                   {
  return m_HadCore;
}

HadIsolET()

int LVL1::CPMTobAlgorithm::HadIsolET

(

)

Returns Had isolation ET.

Definition at line 404 of file CPMTobAlgorithm.cxx.

                                   {
  return m_HadIsol;
}

isEMRoI()

bool LVL1::CPMTobAlgorithm::isEMRoI

(

)

Report whether this passed as an EM TOB.

Definition at line 356 of file CPMTobAlgorithm.cxx.

                                 {
  
  bool passed = ( m_EMThresh && m_EtMax );
 
  return passed;
}

isEtMax()

bool LVL1::CPMTobAlgorithm::isEtMax

(

)

Does this window pass the local ET maximum condition.

Definition at line 419 of file CPMTobAlgorithm.cxx.

                                  {
  return m_EtMax;
}

isTauRoI()

bool LVL1::CPMTobAlgorithm::isTauRoI

(

)

Report whether this passed as a Tau TOB.

Definition at line 364 of file CPMTobAlgorithm.cxx.

                                  {
 
  bool passed = ( m_TauThresh && m_EtMax );
 
  return passed;
}

phi()

double LVL1::CPMTobAlgorithm::phi

(

)

Returns phi coordinate of RoI, using standard ATLAS convention.

Definition at line 464 of file CPMTobAlgorithm.cxx.

                                {
  return ( (m_phi <= M_PI) ? m_phi : m_phi - 2.*M_PI);
}

setRoICoord()

void LVL1::CPMTobAlgorithm::setRoICoord ( double eta, double phi )

private

threshold values

Compute RoI coordinate.

Input coordinate should be inside "reference tower" of window. Computes RoI coordinate as centre of 2x2 tower RoI core

Definition at line 162 of file CPMTobAlgorithm.cxx.

                                                            {
 
  // TriggerTowerKey::coord() returns centre of TriggerTower
  TriggerTowerKey keyLL(phi,eta);
  Coordinate lowerLeft = keyLL.coord();
 
  // Get centre of tower at eta+1, phi+1, i.e. opposite corner of RoI core
  double dEta = keyLL.dEta(lowerLeft);
  double dPhi = keyLL.dPhi(lowerLeft);  
 
  TriggerTowerKey keyUR(phi+dPhi, eta + dEta); 
  Coordinate upperRight = keyUR.coord();
 
  // CoordinateRange object will compute centre, correcting for wrap-around
  CoordinateRange roi(lowerLeft.phi(),upperRight.phi(),lowerLeft.eta(),upperRight.eta());
  m_eta = roi.eta();
  m_phi = roi.phi();
 
  // Range check (shouldn't be needed, but floating point precision errors possible)
  if       (m_eta > 2.5) m_eta = 2.5;
  else if (m_eta < -2.5) m_eta = -2.5;
}

tauAlgorithm()

void LVL1::CPMTobAlgorithm::tauAlgorithm ( )

private

Check trigger condition and set hits if appropriate.

Definition at line 293 of file CPMTobAlgorithm.cxx.

                                       {
  
  // Initialise to correct defaults, to be sure
  m_TauThresh = false;
  m_TauIsolWord = 0;
  
  // Don't waste time if it isn't an RoI candidate
  if (!m_EtMax) return;
  
  // Does this pass min TOB pT cut?
  // Need to worry about digit scale not being GeV
  float emEnergyScale{1}; 
  unsigned int tobPtMinTau{0};
  emEnergyScale = m_l1menu->thrExtraInfo().EM().emScale();
  tobPtMinTau = m_l1menu->thrExtraInfo().TAU().ptMinToTopo();
  unsigned int thresh = tobPtMinTau * emEnergyScale;
 
  if (m_TauClus <= thresh) return;
  
  // Passed, so there is a TOB here
  m_TauThresh = true;
  
  // Now calculate isolation word
  // Define cluster and isolation values with appropriate scales and range for isolation LUT
  unsigned int clusMask    = ( (1<<m_tauLUT_ClusterNBits)-1 )<<m_tauLUT_ClusterFirstBit;
  unsigned int emIsolMask  = ( (1<<m_tauLUT_EMIsolNBits) -1 )<<m_tauLUT_EMIsolFirstBit;
 
  // If cluster overflows, set all threshold bits
  if (m_TauClus >= clusMask) {
     m_TauIsolWord = 0x1F;
     return;
  }
                 
  // Convert to 100 MeV units to match parameters
  int clus    = (( (m_TauClus < clusMask)    ? (m_TauClus & clusMask)    : clusMask ) * 10)/emEnergyScale;
  int emIsol  = (( (m_EMIsol < emIsolMask)   ? (m_EMIsol  & emIsolMask)  : emIsolMask ) * 10)/emEnergyScale;
  
  // Get isolation values from menu (placeholder code - example logic)
  std::vector<int> emisolcuts(5,m_noIsol);
  for(size_t bit = 1; bit<=5; ++bit) {
    const TrigConf::IsolationLegacy & iso = m_l1menu->thrExtraInfo().TAU().isolation("EMIsoForTAUthr", bit);
    if(! iso.isDefined() ) {
      continue;
    }
    int offset = iso.offset();
    int slope  = iso.slope();
    int mincut = iso.mincut();
    int upperlimit = iso.upperlimit()*10;
    if (clus < upperlimit && bit >= 1 && bit <= 5) {
      // As "slope" = 10* slope, rescale clus too to get fraction right.
      int cut = offset + (slope != 0 ? 10*clus/slope : 0);
      if (cut < mincut) cut = mincut;
      emisolcuts[bit-1] = cut;
    }
  }
 
  // Set isolation bits
  for (unsigned int isol = 0; isol < 5; ++isol) 
    if (emIsol <= emisolcuts[isol]) m_TauIsolWord += (1 << isol);
 
}

TauClusET()

int LVL1::CPMTobAlgorithm::TauClusET

(

)

Returns Tau cluster ET, limited to 8 bits.

Definition at line 384 of file CPMTobAlgorithm.cxx.

                                   {
  return ( (m_TauClus < m_maxClus) ? m_TauClus : m_maxClus );
}

TauCPMTobRoI()

xAOD::CPMTobRoI * LVL1::CPMTobAlgorithm::TauCPMTobRoI

(

)

Returns CPMTobRoI for TAU hypothesis, provided TAU Tob conditions passed.

Definition at line 480 of file CPMTobAlgorithm.cxx.

                                                 {
  
  xAOD::CPMTobRoI* pointer = 0;
  
  if (isTauRoI()) pointer = createTobRoI(TrigT1CaloDefs::tauTobType);
  
  return pointer;
 
}

TauIsolWord()

int LVL1::CPMTobAlgorithm::TauIsolWord

(

)

Returns Tau isolation word (5 bits)

Definition at line 394 of file CPMTobAlgorithm.cxx.

                                     {
  return m_TauIsolWord;
}

TauLUTClus()

int LVL1::CPMTobAlgorithm::TauLUTClus

(

)

Returns Tau cluster ET with range and precision of LUT input.

Definition at line 431 of file CPMTobAlgorithm.cxx.

                                    {
  int ET = (m_TauClus >> m_tauLUT_ClusterFirstBit);
  if (ET > (1<<m_tauLUT_ClusterNBits) - 1) ET = (1<<m_tauLUT_ClusterNBits) - 1;
  return (ET << m_tauLUT_ClusterFirstBit);
}

TauLUTEMIsol()

int LVL1::CPMTobAlgorithm::TauLUTEMIsol

(

)

Returns EM isolation ET with range and precision of Tau LUT input.

Definition at line 452 of file CPMTobAlgorithm.cxx.

                                      {
  int ET = (m_EMIsol >> m_tauLUT_EMIsolFirstBit);
  if (ET > (1<<m_tauLUT_EMIsolNBits) - 1) ET = (1<<m_tauLUT_EMIsolNBits) - 1;
  return (ET << m_tauLUT_EMIsolFirstBit);
}

testEtMax()

void LVL1::CPMTobAlgorithm::testEtMax ( const std::vector< unsigned int > & cores )

private

Form all 2x2 clusters within window and test centre is a local ET maximum.

Clusters form a list, arranged as
2 5 8
1 4 7
0 3 6
Then 4 = this window's RoI core, and ET max condition is
4 >= 0, 1, 2, 3 && 4 < 5, 6, 7, 8

Definition at line 186 of file CPMTobAlgorithm.cxx.

                                                                        {
  
  
  // RoI condition test
  m_EtMax = true;
  for (int i = 0; i < 4; i++) if (m_Core < cores[i])  m_EtMax = false;
  for (int i = 5; i < 9; i++) if (m_Core <= cores[i]) m_EtMax = false;
  
}

Member Data Documentation

m_Core

unsigned int LVL1::CPMTobAlgorithm::m_Core

private

Definition at line 91 of file CPMTobAlgorithm.h.

m_EMClus

unsigned int LVL1::CPMTobAlgorithm::m_EMClus

private

Definition at line 92 of file CPMTobAlgorithm.h.

m_EMCore

unsigned int LVL1::CPMTobAlgorithm::m_EMCore

private

Definition at line 96 of file CPMTobAlgorithm.h.

m_EMIsol

unsigned int LVL1::CPMTobAlgorithm::m_EMIsol

private

Definition at line 97 of file CPMTobAlgorithm.h.

m_EMIsolWord

unsigned int LVL1::CPMTobAlgorithm::m_EMIsolWord

private

Definition at line 94 of file CPMTobAlgorithm.h.

m_emLUT_ClusterFirstBit

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_ClusterFirstBit = 1

staticprivate

Definition at line 109 of file CPMTobAlgorithm.h.

m_emLUT_ClusterNBits

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_ClusterNBits = 6

staticprivate

Definition at line 110 of file CPMTobAlgorithm.h.

m_emLUT_EMIsolFirstBit

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_EMIsolFirstBit = 0

staticprivate

Definition at line 111 of file CPMTobAlgorithm.h.

m_emLUT_EMIsolNBits

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_EMIsolNBits = 4

staticprivate

Definition at line 112 of file CPMTobAlgorithm.h.

m_emLUT_HadVetoFirstBit

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_HadVetoFirstBit = 0

staticprivate

Definition at line 113 of file CPMTobAlgorithm.h.

m_emLUT_HadVetoNBits

const unsigned int LVL1::CPMTobAlgorithm::m_emLUT_HadVetoNBits = 3

staticprivate

Definition at line 114 of file CPMTobAlgorithm.h.

m_EMThresh

bool LVL1::CPMTobAlgorithm::m_EMThresh

private

Definition at line 101 of file CPMTobAlgorithm.h.

m_eta

double LVL1::CPMTobAlgorithm::m_eta

private

Algorithm results.

Definition at line 89 of file CPMTobAlgorithm.h.

m_EtMax

bool LVL1::CPMTobAlgorithm::m_EtMax

private

Definition at line 100 of file CPMTobAlgorithm.h.

m_HadCore

unsigned int LVL1::CPMTobAlgorithm::m_HadCore

private

Definition at line 98 of file CPMTobAlgorithm.h.

m_HadIsol

unsigned int LVL1::CPMTobAlgorithm::m_HadIsol

private

Definition at line 99 of file CPMTobAlgorithm.h.

m_l1menu

const TrigConf::L1Menu* LVL1::CPMTobAlgorithm::m_l1menu {nullptr}

private

Definition at line 86 of file CPMTobAlgorithm.h.

{nullptr};

m_maxClus

const unsigned int LVL1::CPMTobAlgorithm::m_maxClus = 0xFF

staticprivate

Algorithm parameters.

Definition at line 107 of file CPMTobAlgorithm.h.

m_maxEta

const double LVL1::CPMTobAlgorithm::m_maxEta = 2.5

staticprivate

Definition at line 123 of file CPMTobAlgorithm.h.

m_noIsol

const unsigned int LVL1::CPMTobAlgorithm::m_noIsol = 999

staticprivate

Definition at line 121 of file CPMTobAlgorithm.h.

m_phi

double LVL1::CPMTobAlgorithm::m_phi

private

Definition at line 90 of file CPMTobAlgorithm.h.

m_refEta

double LVL1::CPMTobAlgorithm::m_refEta

private

Definition at line 84 of file CPMTobAlgorithm.h.

m_refPhi

double LVL1::CPMTobAlgorithm::m_refPhi

private

Definition at line 85 of file CPMTobAlgorithm.h.

m_TauClus

unsigned int LVL1::CPMTobAlgorithm::m_TauClus

private

Definition at line 93 of file CPMTobAlgorithm.h.

m_TauIsolWord

unsigned int LVL1::CPMTobAlgorithm::m_TauIsolWord

private

Definition at line 95 of file CPMTobAlgorithm.h.

m_tauLUT_ClusterFirstBit

const unsigned int LVL1::CPMTobAlgorithm::m_tauLUT_ClusterFirstBit = 1

staticprivate

Definition at line 116 of file CPMTobAlgorithm.h.

m_tauLUT_ClusterNBits

const unsigned int LVL1::CPMTobAlgorithm::m_tauLUT_ClusterNBits = 7

staticprivate

Definition at line 117 of file CPMTobAlgorithm.h.

m_tauLUT_EMIsolFirstBit

const unsigned int LVL1::CPMTobAlgorithm::m_tauLUT_EMIsolFirstBit = 0

staticprivate

Definition at line 118 of file CPMTobAlgorithm.h.

m_tauLUT_EMIsolNBits

const unsigned int LVL1::CPMTobAlgorithm::m_tauLUT_EMIsolNBits = 6

staticprivate

Definition at line 119 of file CPMTobAlgorithm.h.

m_TauThresh

bool LVL1::CPMTobAlgorithm::m_TauThresh

private

Definition at line 102 of file CPMTobAlgorithm.h.

---

The documentation for this class was generated from the following files:

• CPMTobAlgorithm.h
• CPMTobAlgorithm.cxx