L1GctHardwareJetFinder.h

Go to the documentation of this file.

#ifndef L1GCTHARDWAREJETFINDER_H_
#define L1GCTHARDWAREJETFINDER_H_

#include "L1Trigger/GlobalCaloTrigger/interface/L1GctJetFinderBase.h"

#include <boost/cstdint.hpp> //for uint16_t
#include <vector>

/*
 * \author Greg Heath
 * \date June 2006
 */



class L1GctHardwareJetFinder : public L1GctJetFinderBase
{
 public:

  L1GctHardwareJetFinder(int id);
                 
  ~L1GctHardwareJetFinder();
   
  friend std::ostream& operator << (std::ostream& os, const L1GctHardwareJetFinder& algo);

  virtual void reset();

  virtual void fetchInput();

  virtual void process();

 protected:

  // Each jetFinder must define the constants as private and copy the
  // function definitions below.
  virtual unsigned maxRegionsIn() const { return MAX_REGIONS_IN; }
  virtual unsigned centralCol0() const { return CENTRAL_COL0; }
  virtual unsigned nCols() const { return N_COLS; }

private:

  static const unsigned int MAX_REGIONS_IN; 
  static const unsigned int N_COLS;
  static const unsigned int CENTRAL_COL0;

  bool m_positiveEtaWheel;

  RegionsVector m_localMaxima;
  RegionsVector m_clusters;

  unsigned m_numberOfClusters;

  // Additional clusters to avoid double counting of jets across eta=0
  RegionsVector m_localMax00;
  RegionsVector  m_cluster00;
  
  void findProtoJets();  
  L1GctRegion makeProtoJet(L1GctRegion localMax);
  void findJets();  

  void findLocalMaxima();
  void findProtoClusters();
  void findFinalClusters();

  void convertClustersToProtoJets();
  void convertClustersToOutputJets();

};

std::ostream& operator << (std::ostream& os, const L1GctHardwareJetFinder& algo);

#endif /*L1GCTHARDWAREJETFINDER_H_*/