EcalFenixTcp.h

Go to the documentation of this file.

#ifndef ECAL_FENIX_TCP_H
#define ECAL_FENIX_TCP_H

#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixBypassLin.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixEtTot.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixFgvbEB.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixMaxof2.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixTcpFgvbEE.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixTcpFormat.h>
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixTcpsFgvbEB.h>

#include <DataFormats/EcalDigi/interface/EBDataFrame.h>
#include <DataFormats/EcalDigi/interface/EEDataFrame.h>
#include <DataFormats/EcalDigi/interface/EcalTriggerPrimitiveSample.h>

#include "FWCore/Framework/interface/EventSetup.h"

#include <iostream>
#include <vector>

class EcalTPGFineGrainEBGroup;
class EcalTPGLutGroup;
class EcalTPGLutIdMap;
class EcalTPGFineGrainEBIdMap;
class EcalTPGFineGrainTowerEE;
class EcalTrigTowerDetId;
class EcalTPGTowerStatus;

class EcalFenixTcp {
private:
  bool debug_;

  int nbMaxStrips_;

  EcalFenixMaxof2 *maxOf2_;
  std::vector<EcalFenixBypassLin *> bypasslin_;
  EcalFenixEtTot *adder_;
  EcalFenixFgvbEB *fgvbEB_;
  EcalFenixTcpFgvbEE *fgvbEE_;
  EcalFenixTcpsFgvbEB *sfgvbEB_;

  EcalFenixTcpFormat *formatter_;

  // permanent data structures
  std::vector<std::vector<int>> bypasslin_out_;
  std::vector<int> adder_out_;
  std::vector<int> maxOf2_out_;
  std::vector<int> fgvb_out_;
  std::vector<int> strip_fgvb_out_;

public:
  // temporary, for timing tests
  void setPointers(const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup,
                   const EcalTPGLutGroup *ecaltpgLutGroup,
                   const EcalTPGLutIdMap *ecaltpgLut,
                   const EcalTPGFineGrainEBIdMap *ecaltpgFineGrainEB,
                   const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE,
                   const EcalTPGTowerStatus *ecaltpgBadTT,
                   const EcalTPGSpike *ecaltpgSpike) {
    ecaltpgFgEBGroup_ = ecaltpgFgEBGroup;
    ecaltpgLutGroup_ = ecaltpgLutGroup;
    ecaltpgLut_ = ecaltpgLut;
    ecaltpgFineGrainEB_ = ecaltpgFineGrainEB;
    ecaltpgFineGrainTowerEE_ = ecaltpgFineGrainTowerEE;
    ecaltpgBadTT_ = ecaltpgBadTT;
    ecaltpgSpike_ = ecaltpgSpike;
  }
  // end temporary, for timing tests

  EcalFenixTcp(const edm::EventSetup &setup,
               bool tcpFormat,
               bool debug,
               bool famos,
               int binOfMax,
               int maxNrSamples,
               int nbMaxStrips);
  virtual ~EcalFenixTcp();

  void process(const edm::EventSetup &setup,
               std::vector<EBDataFrame> &bid,  // dummy argument for template call
               std::vector<std::vector<int>> &tpframetow,
               int nStr,
               std::vector<EcalTriggerPrimitiveSample> &tptow,
               std::vector<EcalTriggerPrimitiveSample> &tptow2,
               bool isInInnerRings,
               EcalTrigTowerDetId thisTower);
  void process(const edm::EventSetup &setup,
               std::vector<EEDataFrame> &bid,  // dummy argument for template call
               std::vector<std::vector<int>> &tpframetow,
               int nStr,
               std::vector<EcalTriggerPrimitiveSample> &tptow,
               std::vector<EcalTriggerPrimitiveSample> &tptow2,
               bool isInInnerRings,
               EcalTrigTowerDetId thisTower);

  void process_part1(std::vector<std::vector<int>> &tpframetow, int nStr, int bitMask);

  void process_part2_barrel(std::vector<std::vector<int>> &,
                            int nStr,
                            int bitMask,
                            const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup,
                            const EcalTPGLutGroup *ecaltpgLutGroup,
                            const EcalTPGLutIdMap *ecaltpgLut,
                            const EcalTPGFineGrainEBIdMap *ecaltpgFineGrainEB,
                            const EcalTPGTowerStatus *ecaltpgBadTT,
                            const EcalTPGSpike *ecaltpgSpike,
                            std::vector<EcalTriggerPrimitiveSample> &tptow,
                            std::vector<EcalTriggerPrimitiveSample> &tptow2,
                            EcalTrigTowerDetId towid);

  void process_part2_endcap(std::vector<std::vector<int>> &,
                            int nStr,
                            int bitMask,
                            const EcalTPGLutGroup *ecaltpgLutGroup,
                            const EcalTPGLutIdMap *ecaltpgLut,
                            const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE,
                            const EcalTPGTowerStatus *ecaltpgBadTT,
                            std::vector<EcalTriggerPrimitiveSample> &tptow,
                            std::vector<EcalTriggerPrimitiveSample> &tptow2,
                            bool isInInnerRings,
                            EcalTrigTowerDetId towid);

  EcalFenixBypassLin *getBypasslin(int i) const { return bypasslin_[i]; }
  EcalFenixEtTot *getAdder() const { return adder_; }
  EcalFenixMaxof2 *getMaxOf2() const { return maxOf2_; }
  EcalFenixTcpFormat *getFormatter() const { return formatter_; }
  EcalFenixFgvbEB *getFGVBEB() const { return fgvbEB_; }
  EcalFenixTcpFgvbEE *getFGVBEE() const { return fgvbEE_; }
  EcalFenixTcpsFgvbEB *getsFGVBEB() const { return sfgvbEB_; }

  const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup_;
  const EcalTPGLutGroup *ecaltpgLutGroup_;
  const EcalTPGLutIdMap *ecaltpgLut_;
  const EcalTPGFineGrainEBIdMap *ecaltpgFineGrainEB_;
  const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE_;
  const EcalTPGTowerStatus *ecaltpgBadTT_;
  const EcalTPGSpike *ecaltpgSpike_;
};

#endif