calibration.h

#ifndef CALIBRATION_H_INCLUDED
#define CALIBRATION_H_INCLUDED
 
/********************************
 * ReadoutBoard calibration:
 * - convert adc, time bins in 
 *   mV and nanoseconds.
 *******************************/
 
#include <vector>
#include <string>
#include <map>
 
#include "tof_typedefs.h"
#include "events.h"
 
class RBEvent;
 
 
void spike_cleaning_drs4(Vec<Vec<f32>> &wf, u16 tCell, i32 spikes[]);
 
void spike_cleaning_simple(Vec<Vec<f32>> &voltages, bool calibrated = true);
 
void spike_cleaning_all(Vec<Vec<f32>> &voltages, bool calibrated = true);
 
struct RBCalibration {
  static constexpr u16 HEAD = 0xAAAA;
  static constexpr u16 TAIL = 0x5555;
  static bool serialize_event_data;
 
  u8 rb_id;
  f32 d_v;
  u32 timestamp;
  Vec<Vec<f32>> v_offsets;
  Vec<Vec<f32>> v_dips;
  Vec<Vec<f32>> v_incs;
  Vec<Vec<f32>> t_bin;
  // data used to calculate calibration constants
  Vec<RBEvent> noi_data;
  Vec<RBEvent> vcal_data;
  Vec<RBEvent> tcal_data;
 
  RBCalibration();
 
  Vec<Vec<f32>> voltages    (const RBEvent &event, bool spike_cleaning = false) const;
  Vec<Vec<f32>> nanoseconds (const RBEvent &event) const;
  
  Vec<f32> voltages   (const RBEvent &event, const u8 channel) const;
  Vec<f32> nanoseconds(const RBEvent &event, const u8 channel) const;
 
  static auto from_bytestream(const Vec<u8> &bytestream, u64 &pos, 
                              bool discard_events = true)
     -> RBCalibration;
 
  [[deprecated("Replaced by from_file, we will no longer write human readable calibration files")]]
  static RBCalibration from_txtfile(const String &filename);
 
 
  static auto from_file(const String &filename, bool discard_events = true)
    -> RBCalibration;
  auto to_string() const -> std::string;
 
  static void disable_eventdata();
 
  private:
 
    bool channel_check(u8 channel) const;
};
 
// FIXME - we want to put everything in this namespace,
// start with the new stuff
namespace Gaps {
  auto load_tof_calibrations(std::string const &pathname) -> std::map<u8, RBCalibration>;
}
 
std::ostream& operator<<(std::ostream& os, const RBCalibration& pck);
 
typedef std::map<u8, RBCalibration> RBCalibrationMap;
 
#endif