CastorSimpleRecAlgo.h

Go to the documentation of this file.

#ifndef CASTORSIMPLERECALGO_H
#define CASTORSIMPLERECALGO_H 1

#include "DataFormats/HcalDigi/interface/CastorDataFrame.h"
#include "DataFormats/HcalRecHit/interface/CastorRecHit.h"
#include "CalibFormats/CastorObjects/interface/CastorCoder.h"
#include "CalibFormats/CastorObjects/interface/CastorCalibrations.h"
#include "CalibCalorimetry/CastorCalib/interface/CastorPulseContainmentCorrection.h"
#include <memory>

class CastorSimpleRecAlgo {
public:
  CastorSimpleRecAlgo(
      int firstSample, int samplesToAdd, bool correctForTimeslew, bool correctForContainment, float fixedPhaseNs);
  CastorSimpleRecAlgo(int firstSample, int samplesToAdd);

  CastorRecHit reconstruct(const CastorDataFrame& digi,
                           const CastorCoder& coder,
                           const CastorCalibrations& calibs) const;

  // sets rechit saturation status bit on if ADC count is >= maxADCvalue
  void checkADCSaturation(CastorRecHit& rechit, const CastorDataFrame& digi, const int& maxADCvalue) const;

  //++++ Saturation Correction +++++
  // recover pulse shape if ADC count is >= masADCvalue
  void recoverADCSaturation(CastorRecHit& rechit,
                            const CastorCoder& coder,
                            const CastorCalibrations& calibs,
                            const CastorDataFrame& digi,
                            const int& maxADCvalue,
                            const double& satCorrConst) const;

  void resetTimeSamples(int f, int t) {
    firstSample_ = f;
    samplesToAdd_ = t;
  }

private:
  int firstSample_, samplesToAdd_;
  bool correctForTimeslew_;
  std::unique_ptr<CastorPulseContainmentCorrection> pulseCorr_;
};

#endif