gctpf Namespace Reference

Classes
struct	GCTEtaStrip_t
struct	GCTEtaStripPeak_t
struct	GCTint_t
struct	GCTpfcluster_t
struct	PFcluster_t
struct	Region_t

Functions
GCTint_t	bestOf2 (const GCTint_t &t0, const GCTint_t &t1)
float	getEt (float temp[nTowerEtaSLR][nTowerPhiSLR], int eta, int phi)
GCTint_t	getPeakBin (const GCTEtaStripPeak_t &etaStripPeak)
GCTint_t	getPeakOfStrip (const GCTEtaStrip_t &etaStrip)
GCTint_t	getPeakPosition (const Region_t &region)
Region_t	initStructure (float temp[nTowerEtaSLR][nTowerPhiSLR])
PFcluster_t	pfcluster (float temporary[nTowerEtaSLR][nTowerPhiSLR], int etaoffset, int phioffset)
GCTpfcluster_t	recoPfcluster (float temporary[nTowerEtaSLR][nTowerPhiSLR], int etaoffset, int phioffset)
void	RemoveTmp (float temp[nTowerEtaSLR][nTowerPhiSLR], int eta, int phi)

Function Documentation

bestOf2()

GCTint_t gctpf::bestOf2 ( const GCTint_t &  t0, const GCTint_t &  t1  )

inline

Definition at line 57 of file Phase2L1CaloPFClusterEmulator.h.

References FrontierCondition_GT_autoExpress_cfi::t0, RandomServiceHelper::t1, and x.

Referenced by getPeakBin(), and getPeakOfStrip().

                                                                  {
    GCTint_t x;
    x = (t0.et > t1.et) ? t0 : t1;

    return x;
  }

getEt()

float gctpf::getEt ( float  temp[nTowerEtaSLR][nTowerPhiSLR], int  eta, int  phi  )

inline

Definition at line 121 of file Phase2L1CaloPFClusterEmulator.h.

References PVValHelper::eta, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, nTowerEtaSLR, nTowerPhiSLR, and groupFilesInBlocks::temp.

Referenced by ReflectedIterator::calculatePedestal(), HLTGenericFilter< T1 >::hltFilter(), recoPfcluster(), and ReflectedIterator::subtractPedestal().

                                                                               {
    float et_sumEta[3] = {0., 0., 0.};

    for (int i = 0; i < (nTowerEtaSLR - 2); i++) {
      for (int j = 0; j < (nTowerPhiSLR - 2); j++) {
        if (i + 1 == eta && j + 1 == phi) {
          for (int k = 0; k < 3; k++) {
            et_sumEta[k] = temp[i + k][j] + temp[i + k][j + 1] + temp[i + k][j + 2];
          }
        }
      }
    }

    float pfcluster_et = et_sumEta[0] + et_sumEta[1] + et_sumEta[2];

    return pfcluster_et;
  }

getPeakBin()

GCTint_t gctpf::getPeakBin ( const GCTEtaStripPeak_t &  etaStripPeak )

inline

Definition at line 74 of file Phase2L1CaloPFClusterEmulator.h.

References bestOf2(), and gctpf::GCTEtaStripPeak_t::p.

Referenced by getPeakPosition().

                                                                    {
    GCTint_t best01 = bestOf2(etaStripPeak.p[0], etaStripPeak.p[1]);
    GCTint_t best23 = bestOf2(etaStripPeak.p[2], etaStripPeak.p[3]);
    GCTint_t best45 = bestOf2(etaStripPeak.p[4], etaStripPeak.p[5]);
    GCTint_t best67 = bestOf2(etaStripPeak.p[6], etaStripPeak.p[7]);
    GCTint_t best89 = bestOf2(etaStripPeak.p[8], etaStripPeak.p[9]);
    GCTint_t best1011 = bestOf2(etaStripPeak.p[10], etaStripPeak.p[11]);
    GCTint_t best1213 = bestOf2(etaStripPeak.p[12], etaStripPeak.p[13]);
    GCTint_t best1415 = bestOf2(etaStripPeak.p[14], etaStripPeak.p[15]);
    GCTint_t best1617 = bestOf2(etaStripPeak.p[16], etaStripPeak.p[17]);
    GCTint_t best0123 = bestOf2(best01, best23);
    GCTint_t best4567 = bestOf2(best45, best67);
    GCTint_t best891011 = bestOf2(best89, best1011);
    GCTint_t best12131415 = bestOf2(best1213, best1415);
    GCTint_t best01234567 = bestOf2(best0123, best4567);
    GCTint_t best01234567891011 = bestOf2(best01234567, best891011);
    GCTint_t best121314151617 = bestOf2(best12131415, best1617);
    GCTint_t best12131415161718 = bestOf2(best121314151617, etaStripPeak.p[18]);
    GCTint_t bestAll = bestOf2(best01234567891011, best12131415161718);

    return bestAll;
  }

getPeakOfStrip()

GCTint_t gctpf::getPeakOfStrip ( const GCTEtaStrip_t &  etaStrip )

inline

Definition at line 64 of file Phase2L1CaloPFClusterEmulator.h.

References bestOf2(), and gctpf::GCTEtaStrip_t::t.

Referenced by getPeakPosition().

                                                                {
    GCTint_t best12 = bestOf2(etaStrip.t[1], etaStrip.t[2]);
    GCTint_t best34 = bestOf2(etaStrip.t[3], etaStrip.t[4]);
    GCTint_t best56 = bestOf2(etaStrip.t[5], etaStrip.t[6]);
    GCTint_t best1234 = bestOf2(best12, best34);
    GCTint_t bestAll = bestOf2(best1234, best56);

    return bestAll;
  }

getPeakPosition()

GCTint_t gctpf::getPeakPosition ( const Region_t &  region )

inline

Definition at line 97 of file Phase2L1CaloPFClusterEmulator.h.

References getPeakBin(), getPeakOfStrip(), mps_fire::i, SiStripPI::max, nTowerEtaSLR, gctpf::GCTEtaStripPeak_t::p, and nano_mu_digi_cff::region.

Referenced by recoPfcluster().

                                                          {
    GCTEtaStripPeak_t etaPeak;
    for (int i = 0; i < nTowerEtaSLR - 2; i++) {
      etaPeak.p[i] = getPeakOfStrip(region.s[i + 1]);
    }
    GCTint_t max = getPeakBin(etaPeak);

    return max;
  }

initStructure()

Region_t gctpf::initStructure ( float  temp[nTowerEtaSLR][nTowerPhiSLR] )

inline

Definition at line 107 of file Phase2L1CaloPFClusterEmulator.h.

References mps_fire::i, dqmiolumiharvest::j, nTowerEtaSLR, nTowerPhiSLR, and groupFilesInBlocks::temp.

Referenced by recoPfcluster().

                                                                        {
    Region_t r;

    for (int i = 0; i < nTowerPhiSLR; i++) {
      for (int j = 0; j < nTowerEtaSLR; j++) {
        r.s[j].t[i].et = temp[j][i];
        r.s[j].t[i].eta = j;
        r.s[j].t[i].phi = i;
      }
    }

    return r;
  }

pfcluster()

PFcluster_t gctpf::pfcluster ( float  temporary[nTowerEtaSLR][nTowerPhiSLR], int  etaoffset, int  phioffset  )

inline

Definition at line 177 of file Phase2L1CaloPFClusterEmulator.h.

References gctpf::GCTpfcluster_t::et, gctpf::GCTpfcluster_t::eta, gctpf::PFcluster_t::GCTpfclusters, mps_fire::i, nPFClusterSLR, gctpf::GCTpfcluster_t::phi, and recoPfcluster().

Referenced by cms::PFClusterMETProducer::produce(), and Phase2L1CaloPFClusterEmulator::produce().

                                                                                                          {
    GCTpfcluster_t pfcluster[nPFClusterSLR];

    for (int i = 0; i < nPFClusterSLR; i++) {
      pfcluster[i] = recoPfcluster(temporary, etaoffset, phioffset);
    }

    PFcluster_t GCTPfclusters;

    for (int i = 0; i < nPFClusterSLR; i++) {
      GCTPfclusters.GCTpfclusters[i].et = pfcluster[i].et;
      GCTPfclusters.GCTpfclusters[i].eta = pfcluster[i].eta;
      GCTPfclusters.GCTpfclusters[i].phi = pfcluster[i].phi;
    }

    return GCTPfclusters;
  }

recoPfcluster()

GCTpfcluster_t gctpf::recoPfcluster ( float  temporary[nTowerEtaSLR][nTowerPhiSLR], int  etaoffset, int  phioffset  )

inline

Definition at line 153 of file Phase2L1CaloPFClusterEmulator.h.

References gctpf::GCTpfcluster_t::et, gctpf::GCTpfcluster_t::eta, gctpf::GCTint_t::eta, getEt(), getPeakPosition(), initStructure(), nTowerEtaSLR, nTowerPhiSLR, gctpf::GCTpfcluster_t::phi, gctpf::GCTint_t::phi, nano_mu_digi_cff::region, and RemoveTmp().

Referenced by pfcluster().

                                                                                                                 {
    GCTpfcluster_t pfclusterReturn;

    Region_t region;

    region = initStructure(temporary);

    GCTint_t regionMax = getPeakPosition(region);

    float pfcluster_et = getEt(temporary, regionMax.eta, regionMax.phi);

    RemoveTmp(temporary, regionMax.eta, regionMax.phi);

    if (!(regionMax.eta >= 2 && regionMax.eta < (nTowerEtaSLR - 2) && regionMax.phi >= 2 &&
          regionMax.phi < (nTowerPhiSLR - 2)))
      pfcluster_et = 0;

    pfclusterReturn.et = pfcluster_et;
    pfclusterReturn.eta = regionMax.eta - 2 + etaoffset;
    pfclusterReturn.phi = regionMax.phi - 2 + phioffset;

    return pfclusterReturn;
  }

RemoveTmp()

void gctpf::RemoveTmp ( float  temp[nTowerEtaSLR][nTowerPhiSLR], int  eta, int  phi  )

inline

Definition at line 139 of file Phase2L1CaloPFClusterEmulator.h.

References PVValHelper::eta, mps_fire::i, dqmiolumiharvest::j, nTowerEtaSLR, nTowerPhiSLR, and groupFilesInBlocks::temp.

Referenced by recoPfcluster().

                                                                                  {
    for (int i = 0; i < nTowerEtaSLR; i++) {
      if (i + 1 >= eta && i <= eta + 1) {
        for (int j = 0; j < nTowerPhiSLR; j++) {
          if (j + 1 >= phi && j <= phi + 1) {
            temp[i][j] = 0;
          }
        }
      }
    }

    return;
  }