Calibrate.h

Go to the documentation of this file.

#ifndef Calibrate_H
#define Calibrate_H

#include <string>
#include <iostream>
#include <algorithm>

#include "Tool.h"
#include "Stopwatch.h"


class Calibrate: public Tool {


 public:

  static bool validCalibrationDataFound;

  Calibrate();
  bool Initialise(std::string configfile,DataModel &data);
  bool Execute(); 
  bool Finalise(); 


 private:
  void ApplyCalibration(std::vector<SubSample> &samples);
  void ApplyCalibration(const int& PMTid, float& charge, TimeDelta::short_time_t& time,bool& masked);
  void DumpHits(std::vector<SubSample> &samples);
  void DumpHit(const int& PMTid, float& charge, TimeDelta::short_time_t& time,bool& masked);
  void readCalibration(std::string calibrationFileName);
  util::Stopwatch * m_stopwatch;
  std::string m_stopwatch_file_name{""};
  std::string m_calibration_file_name{""};
  
  enum LogLevel {FATAL=-1, ERROR=0, WARN=1, INFO=2, DEBUG1=3, DEBUG2=4, DEBUG3=5};
  
  int m_verbose{ERROR};

  std::stringstream m_ss;

  void StreamToLog(int level) {
    Log(m_ss.str(), level, m_verbose);
    m_ss.str("");
  }
  struct calibrationData
  {
    bool masked; //< mask to define whether a particular PMT should be used or not
    float chargeScaling; //< Factor to multiply charge by
  };
  typedef int PMTChannelID ; //< Unique identifier for a PMT
  std::map<PMTChannelID,calibrationData> calibrationInformation;

};

namespace ranges
{
  template<typename Range, typename Function>
  Function for_each(Range& range, Function f)
  {
    return std::for_each(begin(range),end(range),f);
  }
}
#endif