CaloTPGTranscoderULUT.h

Go to the documentation of this file.

#ifndef CALOTPGTRANSCODERULUT_H
#define CALOTPGTRANSCODERULUT_H 1

#include <memory>
#include <vector>
#include <set>
#include "CalibFormats/CaloTPG/interface/CaloTPGTranscoder.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "CondFormats/HcalObjects/interface/HcalLutMetadata.h"


class HcalTrigTowerGeometry;

class CaloTPGTranscoderULUT : public CaloTPGTranscoder {
public:
  CaloTPGTranscoderULUT(const std::string& compressionFile="",
                        const std::string& decompressionFile="");
  ~CaloTPGTranscoderULUT() override;
  HcalTriggerPrimitiveSample hcalCompress(const HcalTrigTowerDetId& id, unsigned int sample, int fineGrain) const override;
  EcalTriggerPrimitiveSample ecalCompress(const EcalTrigTowerDetId& id, unsigned int sample, bool fineGrain) const override;

  void rctEGammaUncompress(const HcalTrigTowerDetId& hid, const HcalTriggerPrimitiveSample& hc,
                   const EcalTrigTowerDetId& eid, const EcalTriggerPrimitiveSample& ec, 
                   unsigned int& et, bool& egVecto, bool& activity) const override;
  void rctJetUncompress(const HcalTrigTowerDetId& hid, const HcalTriggerPrimitiveSample& hc,
                   const EcalTrigTowerDetId& eid, const EcalTriggerPrimitiveSample& ec, 
                   unsigned int& et) const override;
  double hcaletValue(const int& ieta, const int& iphi, const int& version, const int& compressedValue) const override;
  double hcaletValue(const HcalTrigTowerDetId& hid, const HcalTriggerPrimitiveSample& hc) const override;
  virtual bool HTvalid(const int ieta, const int iphi, const int version) const;
  virtual const std::vector<unsigned int> getCompressionLUT(const HcalTrigTowerDetId& id) const;
  virtual void setup(HcalLutMetadata const&, HcalTrigTowerGeometry const&,
                     int nctScaleShift, int rctScaleShift,
                     double lsbQIE8, double lsbQIE11, bool allLinear);
  virtual int getOutputLUTId(const HcalTrigTowerDetId& id) const;
  virtual int getOutputLUTId(const int ieta, const int iphi, const int version) const;

 private:
  // Constant
  static const int NOUTLUTS = 4176;

  // Two possible linearization scales
  static const unsigned int REDUCE10BIT = 1024;
  static const unsigned int REDUCE11BIT = 2048;

  // Map different QIE to the right linearization
  static const unsigned int QIE8_OUTPUT_LUT_SIZE = REDUCE10BIT;
  static const unsigned int QIE10_OUTPUT_LUT_SIZE = REDUCE11BIT;
  static const unsigned int QIE11_OUTPUT_LUT_SIZE = REDUCE11BIT;
  static const unsigned int OUTPUT_LUT_SIZE = std::max({QIE8_OUTPUT_LUT_SIZE, QIE10_OUTPUT_LUT_SIZE, QIE11_OUTPUT_LUT_SIZE});
  static const unsigned int TPGMAX = 256;

  // Typedef
  typedef uint8_t LUT;
  typedef std::array<float, TPGMAX> RCTdecompression;

  const HcalTopology* theTopology;
  static const bool newHFphi = true;

  // Member functions
  unsigned int getOutputLUTSize(const HcalTrigTowerDetId& id) const;
  bool isOnlyQIE11(const HcalTrigTowerDetId& id) const;
  void loadHCALCompress(HcalLutMetadata const&, HcalTrigTowerGeometry const&) ; //Analytical compression tables

  // Member Variables
  bool allLinear_ = false;
  double nominal_gain_;
  double lsb_factor_;
  double rct_factor_;
  double nct_factor_;
  double lin8_factor_;
  double lin11_factor_;
  const std::string compressionFile_;
  const std::string decompressionFile_;
  std::vector<int> ietal;
  std::vector<int> ietah;
  std::vector<int> ZS;
  std::vector<int> LUTfactor;

  std::vector<std::vector<LUT>> outputLUT_;
  std::vector<RCTdecompression> hcaluncomp_;

  std::set<HcalDetId> plan1_towers_;
};
#endif