HcaluLUTTPGCoder.h

Go to the documentation of this file.

#ifndef CALIBCALORIMETRY_HCALTPGALGOS_HCALNOMINALTPGCODER_H
#define CALIBCALORIMETRY_HCALTPGALGOS_HCALNOMINALTPGCODER_H 1

#include "CalibFormats/HcalObjects/interface/HcalTPGCoder.h"
#include "CalibFormats/HcalObjects/interface/HcalNominalCoder.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include <vector>

class HcalDbService;

class HcaluLUTTPGCoder : public HcalTPGCoder {
public:

  HcaluLUTTPGCoder(const HcalTopology* topo);
  virtual ~HcaluLUTTPGCoder();
  virtual void adc2Linear(const HBHEDataFrame& df, IntegerCaloSamples& ics) const;
  virtual void adc2Linear(const HFDataFrame& df, IntegerCaloSamples& ics) const;
  virtual void compress(const IntegerCaloSamples& ics, const std::vector<bool>& featureBits, HcalTriggerPrimitiveDigi& tp) const;
  virtual unsigned short adc2Linear(HcalQIESample sample,HcalDetId id) const;
  virtual float getLUTPedestal(HcalDetId id) const;
  virtual float getLUTGain(HcalDetId id) const;

  void update(const HcalDbService& conditions);
  void update(const char* filename, bool appendMSB = false);
  void updateXML(const char* filename);
  void setLUTGenerationMode(bool gen){ LUTGenerationMode_ = gen; };
  void setMaskBit(int bit){ bitToMask_ = bit; };
  std::vector<unsigned short> getLinearizationLUTWithMSB(const HcalDetId& id) const;
  void lookupMSB(const HBHEDataFrame& df, std::vector<bool>& msb) const;
  bool getMSB(const HcalDetId& id, int adc) const;
  int getLUTId(HcalSubdetector id, int ieta, int iphi, int depth) const;
  int getLUTId(uint32_t rawid) const;
  int getLUTId(const HcalDetId& detid) const;

private:
  // typedef
  typedef unsigned short LutElement;
  typedef std::vector<LutElement> Lut;

  // constants
  static const size_t INPUT_LUT_SIZE = 128;
  static const int    nFi_ = 72;
  static const float  lsb_;
  
  // member variables
  const HcalTopology* topo_;
  bool LUTGenerationMode_;
  int  bitToMask_;
  int  firstHBEta_, lastHBEta_, nHBEta_, maxDepthHB_, sizeHB_;
  int  firstHEEta_, lastHEEta_, nHEEta_, maxDepthHE_, sizeHE_;
  int  firstHFEta_, lastHFEta_, nHFEta_, maxDepthHF_, sizeHF_;
  std::vector< Lut > inputLUT_;
  std::vector<float> gain_;
  std::vector<float> ped_;
};

#endif