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