CMS 3D CMS Logo

List of all members | Classes | Public Types | Public Member Functions | Static Public Attributes | Protected Types | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | Static Private Attributes | Friends
L1GctJetFinderBase Class Referenceabstract

Base class to allow implementation of jetFinder algorithms. More...

#include <L1GctJetFinderBase.h>

Inheritance diagram for L1GctJetFinderBase:
L1GctProcessor L1GctHardwareJetFinder L1GctNullJetFinder L1GctSimpleJetFinder L1GctTdrJetFinder

Classes

struct  hfTowerSumsType
 

Public Types

typedef L1GctTwosComplement< L1GctInternEtSum::kJetMissEtNBitsetCompInternJfType
 
typedef L1GctUnsignedInt< L1GctInternEtSum::kTotEtOrHtNBitsetHadType
 
typedef L1GctUnsignedInt< L1GctInternEtSum::kTotEtOrHtNBitsetTotalType
 
typedef L1GctTwosComplement< L1GctInternHtMiss::kJetMissHtNBitshtCompInternJfType
 
typedef std::vector< L1GctJetCandJetVector
 
typedef L1GctJet::lutPtr lutPtr
 
typedef std::vector< lutPtrlutPtrVector
 
enum  maxValues { etTotalMaxValue = L1GctInternEtSum::kTotEtOrHtMaxValue, htTotalMaxValue = L1GctInternEtSum::kTotEtOrHtMaxValue }
 
typedef Pipeline< L1GctJetRawJetPipeline
 
typedef std::vector< L1GctJetRawJetVector
 
typedef std::vector< L1GctRegionRegionsVector
 
typedef unsigned long int ULong
 
typedef unsigned short int UShort
 

Public Member Functions

void fetchInput () override=0
 get input data from sources; to be filled in by derived jetFinders More...
 
unsigned getCenJetSeed () const
 
unsigned getEtaBoundry () const
 
etTotalType getEtSum () const
 Get the scalar sum of Et summed over the input regions. More...
 
etCompInternJfType getExSum () const
 Get the x component of vector Et summed over the input regions. More...
 
etCompInternJfType getEySum () const
 Get the y component of vector Et summed over the input regions. More...
 
unsigned getFwdJetSeed () const
 
hfTowerSumsType getHfSums () const
 Get the Hf tower Et sums and tower-over-threshold counts. More...
 
unsigned getHtmSumJetThreshold () const
 
etHadType getHtSum () const
 Get the scalar sum of Ht summed over jets above threshold. More...
 
unsigned getHttSumJetThreshold () const
 
htCompInternJfType getHxSum () const
 Get the x component of vector Ht summed over jets above threshold. More...
 
htCompInternJfType getHySum () const
 Get the y component of vector Ht summed over jets above threshold. More...
 
RegionsVector getInputRegions () const
 Return input data. More...
 
std::vector< L1GctInternEtSumgetInternalEtSums () const
 get et sums in raw format - to be stored in the event More...
 
std::vector< L1GctInternHtMissgetInternalHtMiss () const
 
std::vector< L1GctInternJetDatagetInternalJets () const
 get output jets in raw format - to be stored in the event More...
 
const lutPtrVector getJetEtCalLuts () const
 Return pointers to calibration LUTs. More...
 
JetVector getJets () const
 Get the located jets. More...
 
RegionsVector getKeptProtoJets () const
 get protoJets kept More...
 
RawJetVector getRawJets () const
 get output jets in raw format More...
 
RegionsVector getRcvdProtoJets () const
 get protoJets received from neighbour More...
 
RegionsVector getSentProtoJets () const
 get protoJets sent to neighbour More...
 
unsigned getTauIsolationThreshold () const
 
unsigned getTauJetSeed () const
 
 L1GctJetFinderBase (int id)
 id is 0-8 for -ve Eta jetfinders, 9-17 for +ve Eta, for increasing Phi. More...
 
void process () override=0
 process the data, fill output buffers; to be filled in by derived jetFinders More...
 
void setEnergySumMasks (const L1GctChannelMask *chmask)
 Set masks for energy summing. More...
 
void setInputRegion (const L1CaloRegion &region)
 Set input data. More...
 
void setJetEtCalibrationLuts (const lutPtrVector &jfluts)
 Set pointer to calibration Lut - needed to complete the setup. More...
 
void setJetFinderParams (const L1GctJetFinderParams *jfpars)
 Set pointer to parameters - needed to complete the setup. More...
 
void setNeighbourJetFinders (const std::vector< L1GctJetFinderBase * > &neighbours)
 Set pointers to neighbours - needed to complete the setup. More...
 
bool setupOk () const
 Check setup is Ok. More...
 
void setupTauAlgo (const bool useImprovedAlgo, const bool ignoreVetoBitsForIsolation)
 Setup the tau algorithm parameters. More...
 
 ~L1GctJetFinderBase () override
 
- Public Member Functions inherited from L1GctProcessor
 L1GctProcessor ()
 
void reset ()
 complete reset of processor More...
 
void setBxRange (const int firstBx, const int numberOfBx)
 define the bunch crossing range to process More...
 
void setNextBx (const int bxnum)
 clear input data buffers and process a new bunch crossing More...
 
void setTerse ()
 
bool setupOk () const
 Method to check the setup for this processor. Returns true by default. More...
 
void setVerbose ()
 control output messages More...
 
virtual ~L1GctProcessor ()
 

Static Public Attributes

static const unsigned int COL_OFFSET = L1GctJetFinderParams::NUMBER_ETA_VALUES + N_EXTRA_REGIONS_ETA00
 The index offset between columns. More...
 
static const unsigned int MAX_JETS_OUT = 6
 Max of 6 jets found per jetfinder in a 2*11 search area. More...
 
static const unsigned int N_EXTRA_REGIONS_ETA00 = 2
 Number of additional regions to process on the "wrong" side of eta=0 (determines COL_OFFSET) More...
 
static const unsigned int N_JF_PER_WHEEL = ((L1CaloRegionDetId::N_PHI) / 2)
 No of jetFinders per Wheel. More...
 

Protected Types

enum  fetchType { TOP, BOT, TOPBOT }
 different ways of getting the neighbour data More...
 

Protected Member Functions

etTotalType calcEtStrip (const UShort strip) const
 Calculates total (raw) energy in a phi strip. More...
 
hfTowerSumsType calcHfSums () const
 Calculates Et sum and number of towers over threshold in Hf. More...
 
etTotalType calcHtStrip (const UShort strip) const
 Calculates total calibrated energy in jets (Ht) sum. More...
 
virtual unsigned centralCol0 () const
 
void doEnergySums ()
 Fill the Et strip sums and Ht sum. All jetFinders should call this in process(). More...
 
void doEtSums ()
 Calculates scalar and vector sum of Et over input regions. More...
 
void doHtSums ()
 Calculates scalar and vector sum of Ht over calibrated jets. More...
 
void fetchProtoJetsFromNeighbour (const fetchType ft)
 fetch the protoJets from neighbour jetFinder More...
 
virtual unsigned maxRegionsIn () const
 
virtual unsigned nCols () const
 
void resetPipelines () override
 
void resetProcessor () override
 Separate reset methods for the processor itself and any data stored in pipelines. More...
 
void setupObjects () override
 Initialise inputs with null objects for the correct bunch crossing if required. More...
 
void sortJets ()
 Sort the found jets. All jetFinders should call this in process(). More...
 
- Protected Member Functions inherited from L1GctProcessor
int bxAbs () const
 
int bxMax () const
 
int bxMin () const
 Support for multiple beam crossing operation. More...
 
int bxRel () const
 
int numOfBx () const
 

Protected Attributes

unsigned m_CenJetSeed
 
unsigned m_EtaBoundry
 
bool m_EtmMask [11]
 
bool m_EttMask [11]
 
unsigned m_FwdJetSeed
 
bool m_gotChannelMask
 Remember whether channel mask have been stored. More...
 
bool m_gotJetEtCalLuts
 Remember whether jet Et calibration Lut pointers have been stored. More...
 
bool m_gotJetFinderParams
 Remember whether jetfinder parameters have been stored. More...
 
bool m_gotNeighbourPointers
 Remember whether the neighbour pointers have been stored. More...
 
bool m_HtmMask [11]
 
unsigned m_HtmSumJetThreshold
 
bool m_HttMask [11]
 
unsigned m_HttSumJetThreshold
 
int m_id
 algo ID More...
 
bool m_idInRange
 Remember whether range check on the input ID was ok. More...
 
bool m_ignoreTauVetoBitsForIsolation
 
RegionsVector m_inputRegions
 
lutPtrVector m_jetEtCalLuts
 Jet Et Conversion LUT pointer. More...
 
RegionsVector m_keptProtoJets
 List of pre-clustered jets retained locally as input to the final clustering. More...
 
unsigned m_minColThisJf
 parameter to determine which Regions belong in our acceptance More...
 
std::vector< L1GctJetFinderBase * > m_neighbourJetFinders
 Store neighbour pointers. More...
 
etTotalType m_outputEtSum
 output Et strip sums and Ht - refactored More...
 
etCompInternJfType m_outputExSum
 
etCompInternJfType m_outputEySum
 
hfTowerSumsType m_outputHfSums
 
etHadType m_outputHtSum
 
htCompInternJfType m_outputHxSum
 
htCompInternJfType m_outputHySum
 
RawJetVector m_outputJets
 output jets More...
 
bool m_positiveEtaWheel
 
RegionsVector m_rcvdProtoJets
 List of pre-clustered jets received from neighbour before the final stage of clustering. More...
 
RegionsVector m_sentProtoJets
 List of pre-clustered jets to be sent to neighbour after the first stage of clustering. More...
 
JetVector m_sortedJets
 
unsigned m_tauIsolationThreshold
 
unsigned m_TauJetSeed
 
bool m_useImprovedTauAlgo
 Setup parameters for the tau jet algorithm. More...
 
- Protected Attributes inherited from L1GctProcessor
bool m_verbose
 Flag to control output messages. More...
 

Private Member Functions

template<int kBitsInput, int kBitsOutput>
L1GctTwosComplement< kBitsOutput > etComponentForJetFinder (const L1GctUnsignedInt< kBitsInput > &etStrip0, const unsigned &fact0, const L1GctUnsignedInt< kBitsInput > &etStrip1, const unsigned &fact1)
 Private method for calculating MEt and MHt components. More...
 

Private Attributes

Pipeline< etTotalTypem_outputEtSumPipe
 "Pipeline memories" for energy sums More...
 
Pipeline< etCompInternJfTypem_outputExSumPipe
 
Pipeline< etCompInternJfTypem_outputEySumPipe
 
Pipeline< etHadTypem_outputHtSumPipe
 
Pipeline< htCompInternJfTypem_outputHxSumPipe
 
Pipeline< htCompInternJfTypem_outputHySumPipe
 
RawJetPipeline m_outputJetsPipe
 Output jets "pipeline memory" for checking. More...
 

Static Private Attributes

static const unsigned int CENTRAL_COL0 = 0
 
static const unsigned int MAX_REGIONS_IN = L1GctJetFinderBase::COL_OFFSET * L1GctJetFinderBase::N_COLS
 The real jetFinders must define these constants. More...
 
static const unsigned int N_COLS = 2
 

Friends

std::ostream & operator<< (std::ostream &os, const L1GctJetFinderBase &algo)
 Overload << operator. More...
 

Detailed Description

Base class to allow implementation of jetFinder algorithms.

The base class defines the reset() method, setXxx() and getXxx() methods. Individual jetFinders must define the fetchInput() and process() methods, using protected methods of the base class where necessary.

The jetFinder looks for jets over a 2x11 search area. Its input region are pushed in from the appropriate (phi, eta) range, including across the eta=0 boundary between Wheels. The input regions are copied into a vector of dimension N_COLS*COL_OFFSET.

The array of input regions is filled in a certain order with respect to the index i:

The jetFinder can also pull in "proto-jets" from adjacent jetFinders. If required by the algorithm, these must be calculated in the fetchInput() method;

Definition at line 48 of file L1GctJetFinderBase.h.

Member Typedef Documentation

Definition at line 59 of file L1GctJetFinderBase.h.

Definition at line 58 of file L1GctJetFinderBase.h.

Definition at line 57 of file L1GctJetFinderBase.h.

Definition at line 60 of file L1GctJetFinderBase.h.

Definition at line 55 of file L1GctJetFinderBase.h.

Definition at line 107 of file L1GctJetFinderBase.h.

Definition at line 108 of file L1GctJetFinderBase.h.

Definition at line 56 of file L1GctJetFinderBase.h.

Definition at line 54 of file L1GctJetFinderBase.h.

Definition at line 53 of file L1GctJetFinderBase.h.

typedef unsigned long int L1GctJetFinderBase::ULong

Definition at line 51 of file L1GctJetFinderBase.h.

typedef unsigned short int L1GctJetFinderBase::UShort

Definition at line 52 of file L1GctJetFinderBase.h.

Member Enumeration Documentation

different ways of getting the neighbour data

Enumerator
TOP 
BOT 
TOPBOT 

Definition at line 223 of file L1GctJetFinderBase.h.

Constructor & Destructor Documentation

L1GctJetFinderBase::L1GctJetFinderBase ( int  id)

id is 0-8 for -ve Eta jetfinders, 9-17 for +ve Eta, for increasing Phi.

Definition at line 23 of file L1GctJetFinderBase.cc.

References m_id, m_idInRange, L1GctProcessor::m_verbose, L1CaloRegionDetId::N_PHI, and L1GctProcessor::reset().

24  : L1GctProcessor(),
25  m_id(id),
27  m_idInRange(false),
29  m_gotJetFinderParams(false),
30  m_gotJetEtCalLuts(false),
31  m_gotChannelMask(false),
33  m_minColThisJf(0),
34  m_CenJetSeed(0),
35  m_FwdJetSeed(0),
36  m_TauJetSeed(0),
37  m_EtaBoundry(0),
39  m_useImprovedTauAlgo(false),
44  m_EttMask(),
45  m_EtmMask(),
46  m_HttMask(),
47  m_HtmMask(),
62  // Call reset to initialise vectors for input and output
63  this->reset();
64  //Check jetfinder setup
65  if (m_id < 0 || m_id >= static_cast<int>(L1CaloRegionDetId::N_PHI)) {
66  if (m_verbose) {
67  edm::LogWarning("L1GctSetupError") << "L1GctJetFinderBase::L1GctJetFinderBase() : Jet Finder ID " << m_id
68  << " has been incorrectly constructed!\n"
69  << "ID number should be between the range of 0 to "
70  << L1CaloRegionDetId::N_PHI - 1 << "\n";
71  }
72  } else {
73  m_idInRange = true;
74  }
75 }
bool m_idInRange
Remember whether range check on the input ID was ok.
void reset()
complete reset of processor
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
RegionsVector m_sentProtoJets
List of pre-clustered jets to be sent to neighbour after the first stage of clustering.
bool m_verbose
Flag to control output messages.
bool m_gotChannelMask
Remember whether channel mask have been stored.
Pipeline< etTotalType > m_outputEtSumPipe
"Pipeline memories" for energy sums
bool m_gotNeighbourPointers
Remember whether the neighbour pointers have been stored.
Pipeline< etHadType > m_outputHtSumPipe
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
std::vector< L1GctJetFinderBase * > m_neighbourJetFinders
Store neighbour pointers.
RegionsVector m_keptProtoJets
List of pre-clustered jets retained locally as input to the final clustering.
static const unsigned int MAX_REGIONS_IN
The real jetFinders must define these constants.
RawJetVector m_outputJets
output jets
RegionsVector m_inputRegions
RawJetPipeline m_outputJetsPipe
Output jets "pipeline memory" for checking.
Pipeline< etCompInternJfType > m_outputExSumPipe
bool m_gotJetFinderParams
Remember whether jetfinder parameters have been stored.
Pipeline< htCompInternJfType > m_outputHySumPipe
static const unsigned N_PHI
unsigned m_minColThisJf
parameter to determine which Regions belong in our acceptance
RegionsVector m_rcvdProtoJets
List of pre-clustered jets received from neighbour before the final stage of clustering.
bool m_gotJetEtCalLuts
Remember whether jet Et calibration Lut pointers have been stored.
bool m_useImprovedTauAlgo
Setup parameters for the tau jet algorithm.
Pipeline< etCompInternJfType > m_outputEySumPipe
hfTowerSumsType m_outputHfSums
Pipeline< htCompInternJfType > m_outputHxSumPipe
L1GctJetFinderBase::~L1GctJetFinderBase ( )
override

Definition at line 77 of file L1GctJetFinderBase.cc.

77 {}

Member Function Documentation

etTotalType L1GctJetFinderBase::calcEtStrip ( const UShort  strip) const
protected

Calculates total (raw) energy in a phi strip.

Referenced by nCols().

L1GctJetFinderBase::hfTowerSumsType L1GctJetFinderBase::calcHfSums ( ) const
protected

Calculates Et sum and number of towers over threshold in Hf.

Definition at line 463 of file L1GctJetFinderBase.cc.

References centralCol0(), COL_OFFSET, EgHLTOffHistBins_cfi::et, mps_fire::i, m_inputRegions, nt, hltrates_dqm_sourceclient-live_cfg::offset, and groupFilesInBlocks::temp.

Referenced by doEnergySums(), and nCols().

463  {
464  static const UShort NUMBER_OF_FRWRD_RINGS = 4;
465  static const UShort NUMBER_OF_INNER_RINGS = 2;
466  std::vector<unsigned> et(NUMBER_OF_INNER_RINGS, 0);
467  std::vector<bool> of(NUMBER_OF_INNER_RINGS, false);
468  std::vector<unsigned> nt(NUMBER_OF_INNER_RINGS, 0);
469 
470  UShort offset = COL_OFFSET * (centralCol0() + 1);
471  for (UShort i = 0; i < NUMBER_OF_FRWRD_RINGS; ++i) {
472  offset--;
473 
474  // Sum HF Et and count jets above threshold over "inner rings"
475  if (i < NUMBER_OF_INNER_RINGS) {
476  et.at(i) += m_inputRegions.at(offset).et();
477  of.at(i) = of.at(i) || m_inputRegions.at(offset).overFlow();
478 
479  et.at(i) += m_inputRegions.at(offset + COL_OFFSET).et();
480  of.at(i) = of.at(i) || m_inputRegions.at(offset + COL_OFFSET).overFlow();
481 
482  if (m_inputRegions.at(offset).fineGrain())
483  nt.at(i)++;
484  if (m_inputRegions.at(offset + COL_OFFSET).fineGrain())
485  nt.at(i)++;
486  }
487  }
488  hfTowerSumsType temp(et.at(0), et.at(1), nt.at(0), nt.at(1));
489  temp.etSum0.setOverFlow(temp.etSum0.overFlow() || of.at(0));
490  temp.etSum1.setOverFlow(temp.etSum1.overFlow() || of.at(1));
491  return temp;
492 }
virtual unsigned centralCol0() const
unsigned short int UShort
static const unsigned int COL_OFFSET
The index offset between columns.
RegionsVector m_inputRegions
int nt
Definition: AMPTWrapper.h:42
etTotalType L1GctJetFinderBase::calcHtStrip ( const UShort  strip) const
protected

Calculates total calibrated energy in jets (Ht) sum.

Referenced by nCols().

virtual unsigned L1GctJetFinderBase::centralCol0 ( ) const
inlineprotectedvirtual

Reimplemented in L1GctTdrJetFinder, L1GctHardwareJetFinder, L1GctSimpleJetFinder, and L1GctNullJetFinder.

Definition at line 334 of file L1GctJetFinderBase.h.

References CENTRAL_COL0.

Referenced by calcHfSums(), and doEtSums().

334 { return CENTRAL_COL0; }
static const unsigned int CENTRAL_COL0
void L1GctJetFinderBase::doEnergySums ( )
protected

Fill the Et strip sums and Ht sum. All jetFinders should call this in process().

Definition at line 353 of file L1GctJetFinderBase.cc.

References calcHfSums(), doEtSums(), doHtSums(), and m_outputHfSums.

Referenced by nCols(), L1GctNullJetFinder::process(), L1GctHardwareJetFinder::process(), and L1GctTdrJetFinder::process().

353  {
354  // Refactored energy sums code - find scalar and vector sums
355  // of Et and Ht instead of strip stums
356  doEtSums();
357  doHtSums();
358 
359  //calculate the Hf tower Et sums and tower-over-threshold counts
361 
362  return;
363 }
void doEtSums()
Calculates scalar and vector sum of Et over input regions.
hfTowerSumsType calcHfSums() const
Calculates Et sum and number of towers over threshold in Hf.
void doHtSums()
Calculates scalar and vector sum of Ht over calibrated jets.
hfTowerSumsType m_outputHfSums
void L1GctJetFinderBase::doEtSums ( )
protected

Calculates scalar and vector sum of Et over input regions.

Definition at line 366 of file L1GctJetFinderBase.cc.

References L1GctProcessor::bxRel(), centralCol0(), COL_OFFSET, etTotalMaxValue, mps_fire::i, LEDCalibrationChannels::ieta, m_EttMask, m_id, m_inputRegions, m_outputEtSum, m_outputEtSumPipe, m_outputExSum, m_outputExSumPipe, m_outputEySum, m_outputEySumPipe, N_EXTRA_REGIONS_ETA00, hltrates_dqm_sourceclient-live_cfg::offset, L1GctUnsignedInt< nBits >::overFlow(), L1GctUnsignedInt< nBits >::setOverFlow(), and L1GctUnsignedInt< nBits >::setValue().

Referenced by doEnergySums(), and nCols().

366  {
367  unsigned et0 = 0;
368  unsigned et1 = 0;
369  bool of = false;
370 
371  // Add the Et values from regions 2 to 12 for strip 0,
372  // the Et values from regions 15 to 25 for strip 1.
373  unsigned offset = COL_OFFSET * centralCol0();
374  unsigned ieta = 0;
375  for (UShort i = offset + N_EXTRA_REGIONS_ETA00; i < offset + COL_OFFSET; ++i, ++ieta) {
376  if (!m_EttMask[ieta]) {
377  et0 += m_inputRegions.at(i).et();
378  of |= m_inputRegions.at(i).overFlow();
379  et1 += m_inputRegions.at(i + COL_OFFSET).et();
380  of |= m_inputRegions.at(i + COL_OFFSET).overFlow();
381  }
382  }
383 
384  etTotalType etStrip0(et0);
385  etTotalType etStrip1(et1);
386  etStrip0.setOverFlow(etStrip0.overFlow() || of);
387  etStrip1.setOverFlow(etStrip1.overFlow() || of);
388  unsigned xfact0 = (4 * m_id + 6) % 36;
389  unsigned xfact1 = (4 * m_id + 8) % 36;
390  unsigned yfact0 = (4 * m_id + 15) % 36;
391  unsigned yfact1 = (4 * m_id + 17) % 36;
392  m_outputEtSum = etStrip0 + etStrip1;
393  if (m_outputEtSum.overFlow())
395  m_outputExSum = etComponentForJetFinder<L1GctInternEtSum::kTotEtOrHtNBits, L1GctInternEtSum::kJetMissEtNBits>(
396  etStrip0, xfact0, etStrip1, xfact1);
397  m_outputEySum = etComponentForJetFinder<L1GctInternEtSum::kTotEtOrHtNBits, L1GctInternEtSum::kJetMissEtNBits>(
398  etStrip0, yfact0, etStrip1, yfact1);
399 
402  m_outputEySumPipe.store(m_outputEySum, bxRel());
403 }
etTotalType m_outputEtSum
output Et strip sums and Ht - refactored
int bxRel() const
Pipeline< etTotalType > m_outputEtSumPipe
"Pipeline memories" for energy sums
virtual unsigned centralCol0() const
static const unsigned int N_EXTRA_REGIONS_ETA00
Number of additional regions to process on the "wrong" side of eta=0 (determines COL_OFFSET) ...
unsigned short int UShort
static const unsigned int COL_OFFSET
The index offset between columns.
RegionsVector m_inputRegions
L1GctUnsignedInt< L1GctInternEtSum::kTotEtOrHtNBits > etTotalType
etCompInternJfType m_outputEySum
Pipeline< etCompInternJfType > m_outputExSumPipe
void setValue(unsigned value)
Set value from unsigned.
bool overFlow() const
access overflow
etCompInternJfType m_outputExSum
Pipeline< etCompInternJfType > m_outputEySumPipe
void L1GctJetFinderBase::doHtSums ( )
protected

Calculates scalar and vector sum of Ht over calibrated jets.

Definition at line 406 of file L1GctJetFinderBase.cc.

References L1GctProcessor::bxRel(), htTotalMaxValue, mps_fire::i, LEDCalibrationChannels::ieta, m_HtmMask, m_HtmSumJetThreshold, m_HttMask, m_HttSumJetThreshold, m_id, m_jetEtCalLuts, m_outputHtSum, m_outputHtSumPipe, m_outputHxSum, m_outputHxSumPipe, m_outputHySum, m_outputHySumPipe, m_outputJets, MAX_JETS_OUT, L1GctUnsignedInt< nBits >::overFlow(), L1GctTwosComplement< nBits >::overFlow(), L1GctUnsignedInt< nBits >::setOverFlow(), L1GctTwosComplement< nBits >::setOverFlow(), and L1GctUnsignedInt< nBits >::setValue().

Referenced by doEnergySums(), and nCols().

406  {
407  unsigned htt = 0;
408  unsigned ht0 = 0;
409  unsigned ht1 = 0;
410  bool of = false;
411 
412  for (UShort i = 0; i < MAX_JETS_OUT; ++i) {
413  // Only sum Ht for valid jets
414  if (!m_outputJets.at(i).isNullJet()) {
415  unsigned ieta = m_outputJets.at(i).rctEta();
416  unsigned htJet = m_outputJets.at(i).calibratedEt(m_jetEtCalLuts.at(ieta));
417  // Scalar sum of Htt, with associated threshold
418  if (htJet >= m_HttSumJetThreshold && !m_HttMask[ieta]) {
419  htt += htJet;
420  }
421  // Strip sums, for input to Htm calculation, with associated threshold
422  if (htJet >= m_HtmSumJetThreshold && !m_HtmMask[ieta]) {
423  if (m_outputJets.at(i).rctPhi() == 0) {
424  ht0 += htJet;
425  }
426  if (m_outputJets.at(i).rctPhi() == 1) {
427  ht1 += htJet;
428  }
429  of |= m_outputJets.at(i).overFlow();
430  }
431  }
432  }
433 
434  etHadType httTotal(htt);
435  etHadType htStrip0(ht0);
436  etHadType htStrip1(ht1);
437  httTotal.setOverFlow(httTotal.overFlow() || of);
438  if (httTotal.overFlow())
439  httTotal.setValue(htTotalMaxValue);
440  htStrip0.setOverFlow(htStrip0.overFlow() || of);
441  htStrip1.setOverFlow(htStrip1.overFlow() || of);
442  unsigned xfact0 = (4 * m_id + 10) % 36;
443  unsigned xfact1 = (4 * m_id + 4) % 36;
444  unsigned yfact0 = (4 * m_id + 19) % 36;
445  unsigned yfact1 = (4 * m_id + 13) % 36;
446  m_outputHtSum = httTotal;
447  m_outputHxSum = etComponentForJetFinder<L1GctInternEtSum::kTotEtOrHtNBits, L1GctInternHtMiss::kJetMissHtNBits>(
448  htStrip0, xfact0, htStrip1, xfact1);
449  m_outputHySum = etComponentForJetFinder<L1GctInternEtSum::kTotEtOrHtNBits, L1GctInternHtMiss::kJetMissHtNBits>(
450  htStrip0, yfact0, htStrip1, yfact1);
451 
452  // Common overflow for Ht components
453  bool htmOverFlow = m_outputHxSum.overFlow() || m_outputHySum.overFlow();
454  m_outputHxSum.setOverFlow(htmOverFlow);
455  m_outputHySum.setOverFlow(htmOverFlow);
456 
459  m_outputHySumPipe.store(m_outputHySum, bxRel());
460 }
bool overFlow() const
access overflow
void setOverFlow(bool oflow)
set the overflow bit
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
htCompInternJfType m_outputHySum
int bxRel() const
htCompInternJfType m_outputHxSum
Pipeline< etHadType > m_outputHtSumPipe
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
unsigned short int UShort
L1GctUnsignedInt< L1GctInternEtSum::kTotEtOrHtNBits > etHadType
RawJetVector m_outputJets
output jets
Pipeline< htCompInternJfType > m_outputHySumPipe
Pipeline< htCompInternJfType > m_outputHxSumPipe
template<int kBitsInput, int kBitsOutput>
L1GctTwosComplement< kBitsOutput > L1GctJetFinderBase::etComponentForJetFinder ( const L1GctUnsignedInt< kBitsInput > &  etStrip0,
const unsigned &  fact0,
const L1GctUnsignedInt< kBitsInput > &  etStrip1,
const unsigned &  fact1 
)
private

Private method for calculating MEt and MHt components.

Definition at line 499 of file L1GctJetFinderBase.cc.

References L1GctProcessor::m_verbose, heavyionUCCDQM_cfi::maxEt, L1GctUnsignedInt< nBits >::overFlow(), groupFilesInBlocks::temp, and L1GctUnsignedInt< nBits >::value().

503  {
504  // typedefs and constants
505  typedef L1GctTwosComplement<kBitsOutput> OutputType;
506 
507  // The sin(phi), cos(phi) factors are represented in 15 bits,
508  // as numbers in the range -2^14 to 2^14.
509  // We multiply each input strip Et by the required factor
510  // then shift, to divide by 2^13. This gives an extra bit
511  // of precision on the LSB of the output values.
512  // It's important to avoid systematically biasing the Ex, Ey
513  // component values because this results in an asymmetric
514  // distribution in phi for the final MEt.
515  // The extra LSB is required because one of the factors is 0.5.
516  // Applying this factor without the extra LSB corrects odd values
517  // systematically down by 0.5; or all values by 0.25
518  // on average, giving a shift of -2 units in Ex.
519 
520  static const int internalComponentSize = 15;
521  static const int maxEt = 1 << internalComponentSize;
522 
523  static const int kBitsFactor = internalComponentSize + kBitsInput + 1;
524  static const int maxFactor = 1 << kBitsFactor;
525 
526  static const int bitsToShift = internalComponentSize - 2;
527  static const int halfInputLsb = 1 << (bitsToShift - 1);
528 
529  // These factors correspond to the sine of angles from -90 degrees to
530  // 90 degrees in 10 degree steps, multiplied by 16383 and written
531  // as a <kBitsFactor>-bit 2s-complement number.
532  const int factors[19] = {maxFactor - 16383,
533  maxFactor - 16134,
534  maxFactor - 15395,
535  maxFactor - 14188,
536  maxFactor - 12550,
537  maxFactor - 10531,
538  maxFactor - 8192,
539  maxFactor - 5603,
540  maxFactor - 2845,
541  0,
542  2845,
543  5603,
544  8192,
545  10531,
546  12550,
547  14188,
548  15395,
549  16134,
550  16383};
551 
552  int rotatedValue0, rotatedValue1, myFact;
553  int etComponentSum = 0;
554 
555  if (fact0 >= 36 || fact1 >= 36) {
556  if (m_verbose) {
557  edm::LogError("L1GctProcessingError") << "L1GctJetLeafCard::rotateEtValue() has been called with factor numbers "
558  << fact0 << " and " << fact1 << "; should be less than 36 \n";
559  }
560  } else {
561  // First strip - choose the required multiplication factor
562  if (fact0 > 18) {
563  myFact = factors[(36 - fact0)];
564  } else {
565  myFact = factors[fact0];
566  }
567 
568  // Multiply the Et value by the factor.
569  rotatedValue0 = static_cast<int>(etStrip0.value()) * myFact;
570 
571  // Second strip - choose the required multiplication factor
572  if (fact1 > 18) {
573  myFact = factors[(36 - fact1)];
574  } else {
575  myFact = factors[fact1];
576  }
577 
578  // Multiply the Et value by the factor.
579  rotatedValue1 = static_cast<int>(etStrip1.value()) * myFact;
580 
581  // Add the two scaled values together, with full resolution including
582  // fractional parts from the sin(phi), cos(phi) scaling.
583  // Adjust the value to avoid truncation errors since these
584  // accumulate and cause problems for the missing Et measurement.
585  // Then discard the 13 LSB and interpret the result as
586  // a 15-bit twos complement integer.
587  etComponentSum = ((rotatedValue0 + rotatedValue1) + halfInputLsb) >> bitsToShift;
588 
589  etComponentSum = etComponentSum & (maxEt - 1);
590  if (etComponentSum >= (maxEt / 2)) {
591  etComponentSum = etComponentSum - maxEt;
592  }
593  }
594 
595  // Store as a TwosComplement format integer and return
596  OutputType temp(etComponentSum);
597  temp.setOverFlow(temp.overFlow() || etStrip0.overFlow() || etStrip1.overFlow());
598  return temp;
599 }
bool m_verbose
Flag to control output messages.
unsigned value() const
access value as unsigned
Definition of signed integer types with overflow.
bool overFlow() const
access overflow
void L1GctJetFinderBase::fetchInput ( )
overridepure virtual

get input data from sources; to be filled in by derived jetFinders

Implements L1GctProcessor.

Implemented in L1GctTdrJetFinder, L1GctHardwareJetFinder, L1GctSimpleJetFinder, and L1GctNullJetFinder.

Referenced by L1GctJetLeafCard::fetchInput(), and setupOk().

void L1GctJetFinderBase::fetchProtoJetsFromNeighbour ( const fetchType  ft)
protected

fetch the protoJets from neighbour jetFinder

Helper functions for the fetchInput() and process() methods fetch the protoJets from neighbour jetFinder

Definition at line 312 of file L1GctJetFinderBase.cc.

References BOT, dqmiolumiharvest::j, m_neighbourJetFinders, m_rcvdProtoJets, MAX_JETS_OUT, groupFilesInBlocks::temp, TOP, and TOPBOT.

Referenced by nCols(), and L1GctHardwareJetFinder::process().

312  {
313  switch (ft) {
314  case TOP:
315  m_rcvdProtoJets = m_neighbourJetFinders.at(0)->getSentProtoJets();
316  break;
317  case BOT:
318  m_rcvdProtoJets = m_neighbourJetFinders.at(1)->getSentProtoJets();
319  break;
320  case TOPBOT:
321  // Copy half the jets from each neighbour
322  static const unsigned int MAX_TOPBOT_JETS = MAX_JETS_OUT / 2;
323  unsigned j = 0;
325  temp = m_neighbourJetFinders.at(0)->getSentProtoJets();
326  for (; j < MAX_TOPBOT_JETS; ++j) {
327  m_rcvdProtoJets.at(j) = temp.at(j);
328  }
329  temp = m_neighbourJetFinders.at(1)->getSentProtoJets();
330  for (; j < MAX_JETS_OUT; ++j) {
331  m_rcvdProtoJets.at(j) = temp.at(j);
332  }
333  break;
334  }
335 }
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
std::vector< L1GctJetFinderBase * > m_neighbourJetFinders
Store neighbour pointers.
RegionsVector m_rcvdProtoJets
List of pre-clustered jets received from neighbour before the final stage of clustering.
std::vector< L1GctRegion > RegionsVector
unsigned L1GctJetFinderBase::getCenJetSeed ( ) const
inline

Definition at line 205 of file L1GctJetFinderBase.h.

References m_CenJetSeed.

205 { return m_CenJetSeed; }
unsigned L1GctJetFinderBase::getEtaBoundry ( ) const
inline

Definition at line 208 of file L1GctJetFinderBase.h.

References m_EtaBoundry.

208 { return m_EtaBoundry; }
etTotalType L1GctJetFinderBase::getEtSum ( ) const
inline

Get the scalar sum of Et summed over the input regions.

Definition at line 185 of file L1GctJetFinderBase.h.

References m_outputEtSum.

Referenced by L1GctJetLeafCard::process().

etCompInternJfType L1GctJetFinderBase::getExSum ( ) const
inline

Get the x component of vector Et summed over the input regions.

Definition at line 186 of file L1GctJetFinderBase.h.

References m_outputExSum.

Referenced by L1GctJetLeafCard::process().

etCompInternJfType L1GctJetFinderBase::getEySum ( ) const
inline

Get the y component of vector Et summed over the input regions.

Definition at line 189 of file L1GctJetFinderBase.h.

References m_outputEySum.

Referenced by L1GctJetLeafCard::process().

unsigned L1GctJetFinderBase::getFwdJetSeed ( ) const
inline

Definition at line 206 of file L1GctJetFinderBase.h.

References m_FwdJetSeed.

206 { return m_FwdJetSeed; }
hfTowerSumsType L1GctJetFinderBase::getHfSums ( ) const
inline

Get the Hf tower Et sums and tower-over-threshold counts.

Definition at line 200 of file L1GctJetFinderBase.h.

References m_outputHfSums.

Referenced by L1GctJetLeafCard::process().

unsigned L1GctJetFinderBase::getHtmSumJetThreshold ( ) const
inline
etHadType L1GctJetFinderBase::getHtSum ( ) const
inline

Get the scalar sum of Ht summed over jets above threshold.

Definition at line 192 of file L1GctJetFinderBase.h.

References m_outputHtSum.

Referenced by L1GctJetLeafCard::process().

unsigned L1GctJetFinderBase::getHttSumJetThreshold ( ) const
inline

Definition at line 210 of file L1GctJetFinderBase.h.

References m_HttSumJetThreshold.

210 { return m_HttSumJetThreshold; }
htCompInternJfType L1GctJetFinderBase::getHxSum ( ) const
inline

Get the x component of vector Ht summed over jets above threshold.

Definition at line 193 of file L1GctJetFinderBase.h.

References m_outputHxSum.

Referenced by L1GctJetLeafCard::process().

htCompInternJfType L1GctJetFinderBase::getHySum ( ) const
inline

Get the y component of vector Ht summed over jets above threshold.

Definition at line 196 of file L1GctJetFinderBase.h.

References m_outputHySum.

Referenced by L1GctJetLeafCard::process().

RegionsVector L1GctJetFinderBase::getInputRegions ( ) const
inline

Return input data.

Definition at line 158 of file L1GctJetFinderBase.h.

References m_inputRegions.

158 { return m_inputRegions; }
RegionsVector m_inputRegions
std::vector< L1GctInternEtSum > L1GctJetFinderBase::getInternalEtSums ( ) const

get et sums in raw format - to be stored in the event

Definition at line 280 of file L1GctJetFinderBase.cc.

References l1GtPatternGenerator_cfi::bx, L1GctProcessor::bxMin(), L1GctInternEtSum::fromEmulatorJetMissEt(), L1GctInternEtSum::fromEmulatorJetTotEt(), L1GctInternEtSum::fromEmulatorJetTotHt(), m_outputEtSumPipe, m_outputExSumPipe, m_outputEySumPipe, m_outputHtSumPipe, L1GctProcessor::numOfBx(), and mps_fire::result.

Referenced by getRawJets().

280  {
281  std::vector<L1GctInternEtSum> result;
282  for (int bx = 0; bx < numOfBx(); bx++) {
283  result.push_back(L1GctInternEtSum::fromEmulatorJetTotEt(m_outputEtSumPipe.contents.at(bx).value(),
284  m_outputEtSumPipe.contents.at(bx).overFlow(),
285  static_cast<int16_t>(bx - bxMin())));
286  result.push_back(L1GctInternEtSum::fromEmulatorJetMissEt(m_outputExSumPipe.contents.at(bx).value(),
287  m_outputExSumPipe.contents.at(bx).overFlow(),
288  static_cast<int16_t>(bx - bxMin())));
289  result.push_back(L1GctInternEtSum::fromEmulatorJetMissEt(m_outputEySumPipe.contents.at(bx).value(),
290  m_outputEySumPipe.contents.at(bx).overFlow(),
291  static_cast<int16_t>(bx - bxMin())));
292  result.push_back(L1GctInternEtSum::fromEmulatorJetTotHt(m_outputHtSumPipe.contents.at(bx).value(),
293  m_outputHtSumPipe.contents.at(bx).overFlow(),
294  static_cast<int16_t>(bx - bxMin())));
295  }
296  return result;
297 }
static L1GctInternEtSum fromEmulatorJetTotEt(unsigned totEt, bool overFlow, int16_t bx)
Emulator constructors.
int bxMin() const
Support for multiple beam crossing operation.
int numOfBx() const
Pipeline< etTotalType > m_outputEtSumPipe
"Pipeline memories" for energy sums
static L1GctInternEtSum fromEmulatorJetTotHt(unsigned totHt, bool overFlow, int16_t bx)
Pipeline< etHadType > m_outputHtSumPipe
Pipeline< etCompInternJfType > m_outputExSumPipe
static L1GctInternEtSum fromEmulatorJetMissEt(int missEtxOrEty, bool overFlow, int16_t bx)
Pipeline< etCompInternJfType > m_outputEySumPipe
std::vector< L1GctInternHtMiss > L1GctJetFinderBase::getInternalHtMiss ( ) const

Definition at line 299 of file L1GctJetFinderBase.cc.

References l1GtPatternGenerator_cfi::bx, L1GctProcessor::bxMin(), L1GctInternHtMiss::emulatorJetMissHt(), m_outputHxSumPipe, m_outputHySumPipe, L1GctProcessor::numOfBx(), and mps_fire::result.

Referenced by getRawJets().

299  {
300  std::vector<L1GctInternHtMiss> result;
301  for (int bx = 0; bx < numOfBx(); bx++) {
302  result.push_back(L1GctInternHtMiss::emulatorJetMissHt(m_outputHxSumPipe.contents.at(bx).value(),
303  m_outputHySumPipe.contents.at(bx).value(),
304  m_outputHxSumPipe.contents.at(bx).overFlow(),
305  static_cast<int16_t>(bx - bxMin())));
306  }
307  return result;
308 }
int bxMin() const
Support for multiple beam crossing operation.
int numOfBx() const
static L1GctInternHtMiss emulatorJetMissHt(const int htx, const int hty, const bool overFlow, const int16_t bx)
Named ctor for making missing Ht x & y components object from emulator (jetFinder output)...
Pipeline< htCompInternJfType > m_outputHySumPipe
Pipeline< htCompInternJfType > m_outputHxSumPipe
std::vector< L1GctInternJetData > L1GctJetFinderBase::getInternalJets ( ) const

get output jets in raw format - to be stored in the event

Definition at line 263 of file L1GctJetFinderBase.cc.

References L1GctProcessor::Pipeline< T >::contents, L1GctInternJetData::fromEmulator(), metsig::jet, m_jetEtCalLuts, m_outputJetsPipe, and mps_fire::result.

Referenced by getRawJets().

263  {
264  std::vector<L1GctInternJetData> result;
265  for (RawJetVector::const_iterator jet = m_outputJetsPipe.contents.begin(); jet != m_outputJetsPipe.contents.end();
266  jet++) {
267  result.push_back(L1GctInternJetData::fromEmulator(jet->id(),
268  jet->bx(),
269  jet->calibratedEt(m_jetEtCalLuts.at(jet->rctEta())),
270  jet->overFlow(),
271  jet->tauVeto(),
272  jet->hwEta(),
273  jet->hwPhi(),
274  jet->rank(m_jetEtCalLuts.at(jet->rctEta()))));
275  }
276  return result;
277 }
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
RawJetPipeline m_outputJetsPipe
Output jets "pipeline memory" for checking.
static L1GctInternJetData fromEmulator(L1CaloRegionDetId rgn, int16_t bx, uint16_t et, bool overFlow, bool tauVeto, uint8_t eta, uint8_t phi, uint16_t rank)
std::vector< T > contents
const lutPtrVector L1GctJetFinderBase::getJetEtCalLuts ( ) const
inline

Return pointers to calibration LUTs.

Definition at line 180 of file L1GctJetFinderBase.h.

References m_jetEtCalLuts.

180 { return m_jetEtCalLuts; }
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
JetVector L1GctJetFinderBase::getJets ( ) const
inline
RegionsVector L1GctJetFinderBase::getKeptProtoJets ( ) const
inline

get protoJets kept

Definition at line 167 of file L1GctJetFinderBase.h.

References m_keptProtoJets.

167 { return m_keptProtoJets; }
RegionsVector m_keptProtoJets
List of pre-clustered jets retained locally as input to the final clustering.
RawJetVector L1GctJetFinderBase::getRawJets ( ) const
inline

get output jets in raw format

Definition at line 170 of file L1GctJetFinderBase.h.

References L1GctProcessor::Pipeline< T >::contents, getInternalEtSums(), getInternalHtMiss(), getInternalJets(), and m_outputJetsPipe.

170 { return m_outputJetsPipe.contents; }
RawJetPipeline m_outputJetsPipe
Output jets "pipeline memory" for checking.
std::vector< T > contents
RegionsVector L1GctJetFinderBase::getRcvdProtoJets ( ) const
inline

get protoJets received from neighbour

Definition at line 164 of file L1GctJetFinderBase.h.

References m_rcvdProtoJets.

164 { return m_rcvdProtoJets; }
RegionsVector m_rcvdProtoJets
List of pre-clustered jets received from neighbour before the final stage of clustering.
RegionsVector L1GctJetFinderBase::getSentProtoJets ( ) const
inline

get protoJets sent to neighbour

Definition at line 161 of file L1GctJetFinderBase.h.

References m_sentProtoJets.

161 { return m_sentProtoJets; }
RegionsVector m_sentProtoJets
List of pre-clustered jets to be sent to neighbour after the first stage of clustering.
unsigned L1GctJetFinderBase::getTauIsolationThreshold ( ) const
inline

Definition at line 209 of file L1GctJetFinderBase.h.

References m_tauIsolationThreshold.

209 { return m_tauIsolationThreshold; }
unsigned L1GctJetFinderBase::getTauJetSeed ( ) const
inline

Definition at line 207 of file L1GctJetFinderBase.h.

References m_TauJetSeed.

207 { return m_TauJetSeed; }
virtual unsigned L1GctJetFinderBase::maxRegionsIn ( ) const
inlineprotectedvirtual

*** End of event data ***


Reimplemented in L1GctTdrJetFinder, L1GctHardwareJetFinder, L1GctSimpleJetFinder, and L1GctNullJetFinder.

Definition at line 333 of file L1GctJetFinderBase.h.

References MAX_REGIONS_IN.

Referenced by resetProcessor(), and setupObjects().

333 { return MAX_REGIONS_IN; }
static const unsigned int MAX_REGIONS_IN
The real jetFinders must define these constants.
virtual unsigned L1GctJetFinderBase::nCols ( ) const
inlineprotectedvirtual
void L1GctJetFinderBase::process ( )
overridepure virtual

process the data, fill output buffers; to be filled in by derived jetFinders

Implements L1GctProcessor.

Implemented in L1GctTdrJetFinder, L1GctHardwareJetFinder, L1GctSimpleJetFinder, and L1GctNullJetFinder.

Referenced by L1GctJetLeafCard::process(), and setupOk().

void L1GctJetFinderBase::resetPipelines ( )
overrideprotectedvirtual

Implements L1GctProcessor.

Definition at line 204 of file L1GctJetFinderBase.cc.

References m_outputEtSumPipe, m_outputExSumPipe, m_outputEySumPipe, m_outputHtSumPipe, m_outputHxSumPipe, m_outputHySumPipe, m_outputJetsPipe, L1GctProcessor::numOfBx(), and L1GctProcessor::Pipeline< T >::reset().

Referenced by getHtmSumJetThreshold().

204  {
206  m_outputEtSumPipe.reset(numOfBx());
207  m_outputExSumPipe.reset(numOfBx());
208  m_outputEySumPipe.reset(numOfBx());
209  m_outputHtSumPipe.reset(numOfBx());
210  m_outputHxSumPipe.reset(numOfBx());
211  m_outputHySumPipe.reset(numOfBx());
212 }
int numOfBx() const
Pipeline< etTotalType > m_outputEtSumPipe
"Pipeline memories" for energy sums
Pipeline< etHadType > m_outputHtSumPipe
RawJetPipeline m_outputJetsPipe
Output jets "pipeline memory" for checking.
void reset(const unsigned nBx)
Pipeline< etCompInternJfType > m_outputExSumPipe
Pipeline< htCompInternJfType > m_outputHySumPipe
Pipeline< etCompInternJfType > m_outputEySumPipe
Pipeline< htCompInternJfType > m_outputHxSumPipe
void L1GctJetFinderBase::resetProcessor ( )
overrideprotectedvirtual

Separate reset methods for the processor itself and any data stored in pipelines.

Implements L1GctProcessor.

Definition at line 179 of file L1GctJetFinderBase.cc.

References m_inputRegions, m_keptProtoJets, m_outputEtSum, m_outputExSum, m_outputEySum, m_outputHfSums, m_outputHtSum, m_outputHxSum, m_outputHySum, m_outputJets, m_rcvdProtoJets, m_sentProtoJets, m_sortedJets, MAX_JETS_OUT, maxRegionsIn(), and L1GctJetFinderBase::hfTowerSumsType::reset().

Referenced by L1GctNullJetFinder::fetchInput(), and getHtmSumJetThreshold().

179  {
180  m_inputRegions.clear();
181  m_inputRegions.resize(this->maxRegionsIn());
182  m_outputJets.clear();
183  m_outputJets.resize(MAX_JETS_OUT);
184  m_sortedJets.clear();
185  m_sortedJets.resize(MAX_JETS_OUT);
186 
187  m_sentProtoJets.clear();
189  m_rcvdProtoJets.clear();
191  m_keptProtoJets.clear();
193 
194  m_outputEtSum = 0;
195  m_outputExSum = 0;
196  m_outputEySum = 0;
197  m_outputHtSum = 0;
198  m_outputHxSum = 0;
199  m_outputHySum = 0;
200 
202 }
etTotalType m_outputEtSum
output Et strip sums and Ht - refactored
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
RegionsVector m_sentProtoJets
List of pre-clustered jets to be sent to neighbour after the first stage of clustering.
htCompInternJfType m_outputHySum
htCompInternJfType m_outputHxSum
RegionsVector m_keptProtoJets
List of pre-clustered jets retained locally as input to the final clustering.
RawJetVector m_outputJets
output jets
RegionsVector m_inputRegions
etCompInternJfType m_outputEySum
virtual unsigned maxRegionsIn() const
RegionsVector m_rcvdProtoJets
List of pre-clustered jets received from neighbour before the final stage of clustering.
etCompInternJfType m_outputExSum
hfTowerSumsType m_outputHfSums
void L1GctJetFinderBase::setEnergySumMasks ( const L1GctChannelMask chmask)

Set masks for energy summing.

Set et sum masks from ChannelMask object - needed to complete the setup.

Definition at line 130 of file L1GctJetFinderBase.cc.

References LEDCalibrationChannels::ieta, m_EtmMask, m_EttMask, m_gotChannelMask, m_HtmMask, m_HttMask, m_id, m_positiveEtaWheel, L1GctChannelMask::missingEtMask(), L1GctChannelMask::missingHtMask(), L1GctJetFinderParams::NUMBER_ETA_VALUES, L1GctChannelMask::totalEtMask(), and L1GctChannelMask::totalHtMask().

130  {
131  bool matchCheckEttAndEtm = true;
132  if (chmask != nullptr) {
133  static const unsigned N_ETA = L1GctJetFinderParams::NUMBER_ETA_VALUES;
134  for (unsigned ieta = 0; ieta < N_ETA; ++ieta) {
135  unsigned globalEta = (m_positiveEtaWheel ? N_ETA + ieta : N_ETA - (ieta + 1));
136  m_EttMask[ieta] = chmask->totalEtMask(globalEta);
137  m_EtmMask[ieta] = chmask->missingEtMask(globalEta);
138  m_HttMask[ieta] = chmask->totalHtMask(globalEta);
139  m_HtmMask[ieta] = chmask->missingHtMask(globalEta);
140 
141  matchCheckEttAndEtm &= (m_EttMask[ieta] == m_EtmMask[ieta]);
142  }
143  if (!matchCheckEttAndEtm)
144  edm::LogWarning("L1GctSetupError")
145  << "L1GctJetFinderBase::setEnergySumMasks() : In Jet Finder ID " << m_id
146  << " setting eta-dependent masks for Et sums: you cannot have different masks for total and missing Et\n";
147  m_gotChannelMask = true;
148  }
149 }
bool m_gotChannelMask
Remember whether channel mask have been stored.
bool totalHtMask(unsigned ieta) const
bool missingHtMask(unsigned ieta) const
static const unsigned NUMBER_ETA_VALUES
Number of eta bins used in correction.
bool totalEtMask(unsigned ieta) const
bool missingEtMask(unsigned ieta) const
void L1GctJetFinderBase::setInputRegion ( const L1CaloRegion region)

Set input data.

Definition at line 237 of file L1GctJetFinderBase.cc.

References COL_OFFSET, mps_fire::i, m_id, m_inputRegions, m_minColThisJf, L1GctRegion::makeJfInputRegion(), N_EXTRA_REGIONS_ETA00, N_JF_PER_WHEEL, L1CaloRegionDetId::N_PHI, nCols(), L1CaloRegion::rctCrate(), L1CaloRegion::rctEta(), and L1CaloRegion::rctPhi().

Referenced by setupOk().

237  {
238  static const unsigned NPHI = L1CaloRegionDetId::N_PHI;
239  static const unsigned N_00 = N_EXTRA_REGIONS_ETA00;
240  unsigned crate = region.rctCrate();
241  // Find the column for this region in a global (eta,phi) array
242  // Note the column numbers here are not the same as region->gctPhi()
243  // because the RCT crates are not numbered from phi=0.
244  unsigned colAbsolute = (crate + 1) * 2 + region.rctPhi();
245  unsigned colRelative = ((colAbsolute + NPHI) - m_minColThisJf) % NPHI;
246  if (colRelative < this->nCols()) {
247  // We are in the right range in phi
248  // Now check we are in the right wheel (positive or negative eta)
249  if ((crate / N_JF_PER_WHEEL) == (m_id / N_JF_PER_WHEEL)) {
250  unsigned i = colRelative * COL_OFFSET + N_00 + region.rctEta();
252  } else {
253  // Accept neighbouring regions from the other wheel
254  if (region.rctEta() < N_00) {
255  unsigned i = colRelative * COL_OFFSET + N_00 - (region.rctEta() + 1);
257  }
258  }
259  }
260 }
unsigned rctEta() const
get local eta index (within RCT crate)
Definition: L1CaloRegion.h:147
unsigned rctCrate() const
get RCT crate ID
Definition: L1CaloRegion.h:138
static L1GctRegion makeJfInputRegion(const L1CaloRegion &)
Definition: L1GctRegion.cc:7
unsigned rctPhi() const
get local phi index (within RCT crate)
Definition: L1CaloRegion.h:150
static const unsigned int N_EXTRA_REGIONS_ETA00
Number of additional regions to process on the "wrong" side of eta=0 (determines COL_OFFSET) ...
static const unsigned int COL_OFFSET
The index offset between columns.
RegionsVector m_inputRegions
virtual unsigned nCols() const
static const unsigned N_PHI
unsigned m_minColThisJf
parameter to determine which Regions belong in our acceptance
static const unsigned int N_JF_PER_WHEEL
No of jetFinders per Wheel.
void L1GctJetFinderBase::setJetEtCalibrationLuts ( const lutPtrVector jfluts)

Set pointer to calibration Lut - needed to complete the setup.

Definition at line 124 of file L1GctJetFinderBase.cc.

References m_gotJetEtCalLuts, m_jetEtCalLuts, and L1GctJetFinderParams::NUMBER_ETA_VALUES.

124  {
125  m_jetEtCalLuts = jfluts;
127 }
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
static const unsigned NUMBER_ETA_VALUES
Number of eta bins used in correction.
bool m_gotJetEtCalLuts
Remember whether jet Et calibration Lut pointers have been stored.
void L1GctJetFinderBase::setJetFinderParams ( const L1GctJetFinderParams jfpars)

Set pointer to parameters - needed to complete the setup.

Definition at line 112 of file L1GctJetFinderBase.cc.

References L1GctJetFinderParams::getCenForJetEtaBoundary(), L1GctJetFinderParams::getCenJetEtSeedGct(), L1GctJetFinderParams::getForJetEtSeedGct(), L1GctJetFinderParams::getHtJetEtThresholdGct(), L1GctJetFinderParams::getMHtJetEtThresholdGct(), L1GctJetFinderParams::getTauIsoEtThresholdGct(), L1GctJetFinderParams::getTauJetEtSeedGct(), m_CenJetSeed, m_EtaBoundry, m_FwdJetSeed, m_gotJetFinderParams, m_HtmSumJetThreshold, m_HttSumJetThreshold, m_tauIsolationThreshold, and m_TauJetSeed.

112  {
113  m_CenJetSeed = jfpars->getCenJetEtSeedGct();
114  m_FwdJetSeed = jfpars->getForJetEtSeedGct();
115  m_TauJetSeed = jfpars->getTauJetEtSeedGct();
120  m_gotJetFinderParams = true;
121 }
unsigned getCenJetEtSeedGct() const
unsigned getTauJetEtSeedGct() const
unsigned getTauIsoEtThresholdGct() const
unsigned getMHtJetEtThresholdGct() const
unsigned getCenForJetEtaBoundary() const
unsigned getHtJetEtThresholdGct() const
unsigned getForJetEtSeedGct() const
bool m_gotJetFinderParams
Remember whether jetfinder parameters have been stored.
void L1GctJetFinderBase::setNeighbourJetFinders ( const std::vector< L1GctJetFinderBase * > &  neighbours)

Set pointers to neighbours - needed to complete the setup.

Set pointers to neighbours.

Definition at line 80 of file L1GctJetFinderBase.cc.

References m_gotNeighbourPointers, m_id, m_neighbourJetFinders, and L1GctProcessor::m_verbose.

Referenced by L1GctJetLeafCard::setNeighbourLeafCards().

80  {
82  if (neighbours.size() == 2) {
83  m_neighbourJetFinders = neighbours;
84  } else {
85  m_gotNeighbourPointers = false;
86  if (m_verbose) {
87  edm::LogWarning("L1GctSetupError") << "L1GctJetFinderBase::setNeighbourJetFinders() : In Jet Finder ID " << m_id
88  << " size of input vector should be 2, but is in fact " << neighbours.size()
89  << "\n";
90  }
91  }
92  if (m_neighbourJetFinders.at(0) == nullptr) {
93  m_gotNeighbourPointers = false;
94  if (m_verbose) {
95  edm::LogWarning("L1GctSetupError") << "L1GctJetFinderBase::setNeighbourJetFinders() : In Jet Finder ID " << m_id
96  << " first neighbour pointer is set to zero\n";
97  }
98  }
99  if (m_neighbourJetFinders.at(1) == nullptr) {
100  m_gotNeighbourPointers = false;
101  if (m_verbose) {
102  edm::LogWarning("L1GctSetupError") << "L1GctJetFinderBase::setNeighbourJetFinders() : In Jet Finder ID " << m_id
103  << " second neighbour pointer is set to zero\n";
104  }
105  }
107  edm::LogError("L1GctSetupError") << "Jet Finder ID " << m_id << " has incorrect assignment of neighbour pointers";
108  }
109 }
bool m_verbose
Flag to control output messages.
bool m_gotNeighbourPointers
Remember whether the neighbour pointers have been stored.
std::vector< L1GctJetFinderBase * > m_neighbourJetFinders
Store neighbour pointers.
void L1GctJetFinderBase::setupObjects ( )
overrideprotectedvirtual

Initialise inputs with null objects for the correct bunch crossing if required.

Initialise inputs with null objects for the correct bunch crossing If no other input candidates "arrive", we have the correct bunch crossing to propagate through the processing.

Create a null input region with the right bunch crossing, and fill the input candidates with copies of this.

The same for the lists of pre-clustered jets passed between neighbour jetFinders

The same for the lists of output jets

Implements L1GctProcessor.

Definition at line 217 of file L1GctJetFinderBase.cc.

References L1GctProcessor::bxAbs(), m_inputRegions, m_keptProtoJets, m_outputJets, m_rcvdProtoJets, m_sentProtoJets, MAX_JETS_OUT, maxRegionsIn(), L1GctJet::setBx(), and L1CaloRegion::setBx().

Referenced by L1GctNullJetFinder::fetchInput(), and getHtmSumJetThreshold().

217  {
220  L1GctRegion tempRgn;
221  tempRgn.setBx(bxAbs());
222  m_inputRegions.assign(this->maxRegionsIn(), tempRgn);
223 
226  m_sentProtoJets.assign(MAX_JETS_OUT, tempRgn);
227  m_rcvdProtoJets.assign(MAX_JETS_OUT, tempRgn);
228  m_keptProtoJets.assign(MAX_JETS_OUT, tempRgn);
229 
231  L1GctJet tempJet;
232  tempJet.setBx(bxAbs());
233  m_outputJets.assign(MAX_JETS_OUT, tempJet);
234 }
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
RegionsVector m_sentProtoJets
List of pre-clustered jets to be sent to neighbour after the first stage of clustering.
Gct version of a calorimeter region, used within GCT emulation.
Definition: L1GctRegion.h:17
void setBx(int16_t bx)
set bx
A Level-1 jet candidate, used within GCT emulation.
Definition: L1GctJet.h:26
void setBx(const int16_t bx)
Definition: L1GctJet.h:53
RegionsVector m_keptProtoJets
List of pre-clustered jets retained locally as input to the final clustering.
RawJetVector m_outputJets
output jets
RegionsVector m_inputRegions
int bxAbs() const
virtual unsigned maxRegionsIn() const
RegionsVector m_rcvdProtoJets
List of pre-clustered jets received from neighbour before the final stage of clustering.
bool L1GctJetFinderBase::setupOk ( ) const
inline

Check setup is Ok.

Definition at line 141 of file L1GctJetFinderBase.h.

References ecalcalib_dqm_sourceclient-live_cfg::algo, fetchInput(), m_gotChannelMask, m_gotJetEtCalLuts, m_gotJetFinderParams, m_gotNeighbourPointers, m_idInRange, operator<<, process(), HLT_2018_cff::region, and setInputRegion().

Referenced by L1GctHardwareJetFinder::fetchInput(), L1GctNullJetFinder::process(), L1GctHardwareJetFinder::process(), L1GctTdrJetFinder::process(), and L1GctJetLeafCard::setupOk().

141  {
143  }
bool m_idInRange
Remember whether range check on the input ID was ok.
bool m_gotChannelMask
Remember whether channel mask have been stored.
bool m_gotNeighbourPointers
Remember whether the neighbour pointers have been stored.
bool m_gotJetFinderParams
Remember whether jetfinder parameters have been stored.
bool m_gotJetEtCalLuts
Remember whether jet Et calibration Lut pointers have been stored.
void L1GctJetFinderBase::setupTauAlgo ( const bool  useImprovedAlgo,
const bool  ignoreVetoBitsForIsolation 
)
inline

Setup the tau algorithm parameters.

Definition at line 135 of file L1GctJetFinderBase.h.

References m_ignoreTauVetoBitsForIsolation, and m_useImprovedTauAlgo.

135  {
136  m_useImprovedTauAlgo = useImprovedAlgo;
137  m_ignoreTauVetoBitsForIsolation = ignoreVetoBitsForIsolation;
138  }
bool m_useImprovedTauAlgo
Setup parameters for the tau jet algorithm.
void L1GctJetFinderBase::sortJets ( )
protected

Sort the found jets. All jetFinders should call this in process().

Definition at line 338 of file L1GctJetFinderBase.cc.

References L1GctProcessor::bxRel(), L1GctJetSorter::getSortedJets(), dqmiolumiharvest::j, m_jetEtCalLuts, m_outputJets, m_outputJetsPipe, m_sortedJets, MAX_JETS_OUT, and L1GctProcessor::Pipeline< T >::store().

Referenced by nCols(), L1GctHardwareJetFinder::process(), and L1GctTdrJetFinder::process().

338  {
339  JetVector tempJets(MAX_JETS_OUT);
340  for (unsigned j = 0; j < MAX_JETS_OUT; j++) {
341  tempJets.at(j) = m_outputJets.at(j).jetCand(m_jetEtCalLuts);
342  }
343 
344  // Sort the jets
345  L1GctJetSorter jSorter(tempJets);
346  m_sortedJets = jSorter.getSortedJets();
347 
348  //store jets in "pipeline memory" for checking
350 }
static const unsigned int MAX_JETS_OUT
Max of 6 jets found per jetfinder in a 2*11 search area.
std::vector< L1GctJetCand > JetVector
void store(const T &thisBx, const int bxNum)
int bxRel() const
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
RawJetVector m_outputJets
output jets
RawJetPipeline m_outputJetsPipe
Output jets "pipeline memory" for checking.

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  os,
const L1GctJetFinderBase algo 
)
friend

Overload << operator.

Definition at line 151 of file L1GctJetFinderBase.cc.

Referenced by setupOk().

151  {
152  using std::endl;
153  os << "ID = " << algo.m_id << endl;
154  os << "Calibration lut pointers stored for " << algo.m_jetEtCalLuts.size() << " eta bins" << endl;
155  for (unsigned ieta = 0; ieta < algo.m_jetEtCalLuts.size(); ieta++) {
156  os << "Eta bin " << ieta << ", JetEtCalibrationLut* = " << algo.m_jetEtCalLuts.at(ieta) << endl;
157  }
158  os << "No of input regions " << algo.m_inputRegions.size() << endl;
159  // for(unsigned i=0; i < algo.m_inputRegions.size(); ++i)
160  // {
161  // os << algo.m_inputRegions.at(i);
162  // }
163  os << "No of output jets " << algo.m_outputJets.size() << endl;
164  // for(unsigned i=0; i < algo.m_outputJets.size(); ++i)
165  // {
166  // os << algo.m_outputJets.at(i);
167  // }
168  os << "Output total scalar Et " << algo.m_outputEtSum << endl;
169  os << "Output vector Et x component " << algo.m_outputExSum << endl;
170  os << "Output vector Et y component " << algo.m_outputEySum << endl;
171  os << "Output total scalar Ht " << algo.m_outputHtSum << endl;
172  os << "Output vector Ht x component " << algo.m_outputHxSum << endl;
173  os << "Output vector Ht y component " << algo.m_outputHySum << endl;
174  os << endl;
175 
176  return os;
177 }
etTotalType m_outputEtSum
output Et strip sums and Ht - refactored
htCompInternJfType m_outputHySum
htCompInternJfType m_outputHxSum
lutPtrVector m_jetEtCalLuts
Jet Et Conversion LUT pointer.
RawJetVector m_outputJets
output jets
RegionsVector m_inputRegions
etCompInternJfType m_outputEySum
etCompInternJfType m_outputExSum

Member Data Documentation

const unsigned int L1GctJetFinderBase::CENTRAL_COL0 = 0
staticprivate

Definition at line 364 of file L1GctJetFinderBase.h.

Referenced by centralCol0().

const unsigned int L1GctJetFinderBase::COL_OFFSET = L1GctJetFinderParams::NUMBER_ETA_VALUES + N_EXTRA_REGIONS_ETA00
static
unsigned L1GctJetFinderBase::m_CenJetSeed
protected
unsigned L1GctJetFinderBase::m_EtaBoundry
protected
bool L1GctJetFinderBase::m_EtmMask[11]
protected

Definition at line 290 of file L1GctJetFinderBase.h.

Referenced by setEnergySumMasks().

bool L1GctJetFinderBase::m_EttMask[11]
protected

Definition at line 289 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), and setEnergySumMasks().

unsigned L1GctJetFinderBase::m_FwdJetSeed
protected
bool L1GctJetFinderBase::m_gotChannelMask
protected

Remember whether channel mask have been stored.

Definition at line 244 of file L1GctJetFinderBase.h.

Referenced by setEnergySumMasks(), and setupOk().

bool L1GctJetFinderBase::m_gotJetEtCalLuts
protected

Remember whether jet Et calibration Lut pointers have been stored.

Definition at line 241 of file L1GctJetFinderBase.h.

Referenced by setJetEtCalibrationLuts(), and setupOk().

bool L1GctJetFinderBase::m_gotJetFinderParams
protected

Remember whether jetfinder parameters have been stored.

Definition at line 238 of file L1GctJetFinderBase.h.

Referenced by setJetFinderParams(), and setupOk().

bool L1GctJetFinderBase::m_gotNeighbourPointers
protected

Remember whether the neighbour pointers have been stored.

Definition at line 235 of file L1GctJetFinderBase.h.

Referenced by setNeighbourJetFinders(), and setupOk().

bool L1GctJetFinderBase::m_HtmMask[11]
protected

Definition at line 292 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), and setEnergySumMasks().

unsigned L1GctJetFinderBase::m_HtmSumJetThreshold
protected

Definition at line 286 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getHtmSumJetThreshold(), and setJetFinderParams().

bool L1GctJetFinderBase::m_HttMask[11]
protected

Definition at line 291 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), and setEnergySumMasks().

unsigned L1GctJetFinderBase::m_HttSumJetThreshold
protected

Definition at line 285 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getHttSumJetThreshold(), and setJetFinderParams().

int L1GctJetFinderBase::m_id
protected
bool L1GctJetFinderBase::m_idInRange
protected

Remember whether range check on the input ID was ok.

Definition at line 232 of file L1GctJetFinderBase.h.

Referenced by L1GctJetFinderBase(), and setupOk().

bool L1GctJetFinderBase::m_ignoreTauVetoBitsForIsolation
protected
RegionsVector L1GctJetFinderBase::m_inputRegions
protected
lutPtrVector L1GctJetFinderBase::m_jetEtCalLuts
protected

Jet Et Conversion LUT pointer.

Definition at line 266 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getInternalJets(), getJetEtCalLuts(), operator<<(), setJetEtCalibrationLuts(), and sortJets().

RegionsVector L1GctJetFinderBase::m_keptProtoJets
protected

List of pre-clustered jets retained locally as input to the final clustering.

Definition at line 309 of file L1GctJetFinderBase.h.

Referenced by L1GctHardwareJetFinder::convertClustersToProtoJets(), L1GctHardwareJetFinder::findFinalClusters(), getKeptProtoJets(), resetProcessor(), and setupObjects().

unsigned L1GctJetFinderBase::m_minColThisJf
protected

parameter to determine which Regions belong in our acceptance

Definition at line 253 of file L1GctJetFinderBase.h.

Referenced by L1GctHardwareJetFinder::L1GctHardwareJetFinder(), L1GctTdrJetFinder::L1GctTdrJetFinder(), and setInputRegion().

std::vector<L1GctJetFinderBase*> L1GctJetFinderBase::m_neighbourJetFinders
protected

Store neighbour pointers.

Definition at line 229 of file L1GctJetFinderBase.h.

Referenced by fetchProtoJetsFromNeighbour(), and setNeighbourJetFinders().

etTotalType L1GctJetFinderBase::m_outputEtSum
protected

output Et strip sums and Ht - refactored

Definition at line 316 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getEtSum(), operator<<(), and resetProcessor().

Pipeline<etTotalType> L1GctJetFinderBase::m_outputEtSumPipe
private

"Pipeline memories" for energy sums

Definition at line 370 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getInternalEtSums(), and resetPipelines().

etCompInternJfType L1GctJetFinderBase::m_outputExSum
protected

Definition at line 317 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getExSum(), operator<<(), and resetProcessor().

Pipeline<etCompInternJfType> L1GctJetFinderBase::m_outputExSumPipe
private

Definition at line 371 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getInternalEtSums(), and resetPipelines().

etCompInternJfType L1GctJetFinderBase::m_outputEySum
protected

Definition at line 318 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getEySum(), operator<<(), and resetProcessor().

Pipeline<etCompInternJfType> L1GctJetFinderBase::m_outputEySumPipe
private

Definition at line 372 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), getInternalEtSums(), and resetPipelines().

hfTowerSumsType L1GctJetFinderBase::m_outputHfSums
protected

Definition at line 323 of file L1GctJetFinderBase.h.

Referenced by doEnergySums(), getHfSums(), and resetProcessor().

etHadType L1GctJetFinderBase::m_outputHtSum
protected

Definition at line 319 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getHtSum(), operator<<(), and resetProcessor().

Pipeline<etHadType> L1GctJetFinderBase::m_outputHtSumPipe
private

Definition at line 373 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getInternalEtSums(), and resetPipelines().

htCompInternJfType L1GctJetFinderBase::m_outputHxSum
protected

Definition at line 320 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getHxSum(), operator<<(), and resetProcessor().

Pipeline<htCompInternJfType> L1GctJetFinderBase::m_outputHxSumPipe
private

Definition at line 374 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getInternalHtMiss(), and resetPipelines().

htCompInternJfType L1GctJetFinderBase::m_outputHySum
protected

Definition at line 321 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getHySum(), operator<<(), and resetProcessor().

Pipeline<htCompInternJfType> L1GctJetFinderBase::m_outputHySumPipe
private

Definition at line 375 of file L1GctJetFinderBase.h.

Referenced by doHtSums(), getInternalHtMiss(), and resetPipelines().

RawJetVector L1GctJetFinderBase::m_outputJets
protected
RawJetPipeline L1GctJetFinderBase::m_outputJetsPipe
private

Output jets "pipeline memory" for checking.

Definition at line 367 of file L1GctJetFinderBase.h.

Referenced by getInternalJets(), getRawJets(), resetPipelines(), and sortJets().

bool L1GctJetFinderBase::m_positiveEtaWheel
protected

*** Geometry parameters ***

Positive/negative eta flag

Definition at line 250 of file L1GctJetFinderBase.h.

Referenced by L1GctHardwareJetFinder::findLocalMaxima(), L1GctHardwareJetFinder::makeProtoJet(), and setEnergySumMasks().

RegionsVector L1GctJetFinderBase::m_rcvdProtoJets
protected

List of pre-clustered jets received from neighbour before the final stage of clustering.

Definition at line 307 of file L1GctJetFinderBase.h.

Referenced by fetchProtoJetsFromNeighbour(), L1GctHardwareJetFinder::findFinalClusters(), getRcvdProtoJets(), resetProcessor(), and setupObjects().

RegionsVector L1GctJetFinderBase::m_sentProtoJets
protected

List of pre-clustered jets to be sent to neighbour after the first stage of clustering.

Definition at line 305 of file L1GctJetFinderBase.h.

Referenced by L1GctHardwareJetFinder::convertClustersToProtoJets(), getSentProtoJets(), resetProcessor(), and setupObjects().

JetVector L1GctJetFinderBase::m_sortedJets
protected

Definition at line 313 of file L1GctJetFinderBase.h.

Referenced by getJets(), resetProcessor(), and sortJets().

unsigned L1GctJetFinderBase::m_tauIsolationThreshold
protected
unsigned L1GctJetFinderBase::m_TauJetSeed
protected

Definition at line 262 of file L1GctJetFinderBase.h.

Referenced by getTauJetSeed(), and setJetFinderParams().

bool L1GctJetFinderBase::m_useImprovedTauAlgo
protected

Setup parameters for the tau jet algorithm.

Definition at line 271 of file L1GctJetFinderBase.h.

Referenced by L1GctHardwareJetFinder::findFinalClusters(), L1GctHardwareJetFinder::makeProtoJet(), and setupTauAlgo().

const unsigned int L1GctJetFinderBase::MAX_JETS_OUT = 6
static
const unsigned int L1GctJetFinderBase::MAX_REGIONS_IN = L1GctJetFinderBase::COL_OFFSET * L1GctJetFinderBase::N_COLS
staticprivate

The real jetFinders must define these constants.

Dependent on number of rows and columns.

Definition at line 362 of file L1GctJetFinderBase.h.

Referenced by maxRegionsIn().

const unsigned int L1GctJetFinderBase::N_COLS = 2
staticprivate

Definition at line 363 of file L1GctJetFinderBase.h.

Referenced by nCols().

const unsigned int L1GctJetFinderBase::N_EXTRA_REGIONS_ETA00 = 2
static

Number of additional regions to process on the "wrong" side of eta=0 (determines COL_OFFSET)

Definition at line 115 of file L1GctJetFinderBase.h.

Referenced by doEtSums(), L1GctHardwareJetFinder::findFinalClusters(), L1GctHardwareJetFinder::makeProtoJet(), and setInputRegion().

const unsigned int L1GctJetFinderBase::N_JF_PER_WHEEL = ((L1CaloRegionDetId::N_PHI) / 2)
static

No of jetFinders per Wheel.

Definition at line 113 of file L1GctJetFinderBase.h.

Referenced by setInputRegion().