EcalFenixLinearizer.h

Go to the documentation of this file.

#ifndef SIMCALORIMETRY_ECALTRIGPRIMALGOS_ECALFENIXLINEARIZER_H
#define SIMCALORIMETRY_ECALTRIGPRIMALGOS_ECALFENIXLINEARIZER_H
 
#include <CondFormats/EcalObjects/interface/EcalTPGCrystalStatus.h>
#include <CondFormats/EcalObjects/interface/EcalTPGLinearizationConst.h>
#include <CondFormats/EcalObjects/interface/EcalTPGPedestals.h>
#include <DataFormats/EcalDigi/interface/EcalMGPASample.h>
 
#include <vector>
 
class EcalFenixLinearizer {
private:
  bool famos_;
  int uncorrectedSample_;
  int gainID_;
  int base_;
  int mult_;
  int shift_;
  int strip_;
  bool init_;
 
  const EcalTPGLinearizationConstant *linConsts_;
  const EcalTPGPedestal *peds_;
  const EcalTPGCrystalStatusCode *badXStatus_;
 
  std::vector<const EcalTPGCrystalStatusCode *> vectorbadXStatus_;
 
  int setInput(const EcalMGPASample &RawSam);
  int process();
 
public:
  EcalFenixLinearizer(bool famos);
  virtual ~EcalFenixLinearizer();
 
  template <class T>
  void process(const T &, std::vector<int> &);
  void setParameters(uint32_t raw,
                     const EcalTPGPedestals *ecaltpPed,
                     const EcalTPGLinearizationConst *ecaltpLin,
                     const EcalTPGCrystalStatus *ecaltpBadX);
};
 
template <class T>
void EcalFenixLinearizer::process(const T &df, std::vector<int> &output_percry) {
  // We know a tower numbering is:
  // S1 S2 S3 S4 S5
  //
  // 4  5  14 15 24
  // 3  6  13 16 23
  // 2  7  12 17 22
  // 1  8  11 18 21
  // 0  9  10 19 20
  for (int i = 0; i < df.size(); i++) {
    setInput(df[i]);
    output_percry[i] = process();
  }
 
  return;
}
 
#endif