WCSimReader.h

Go to the documentation of this file.

#ifndef WCSimReader_H
#define WCSimReader_H

#include <string>
#include <iostream>
#include <sstream>

#include "WCSimRootOptions.hh"
#include "WCSimRootEvent.hh"
#include "WCSimRootGeom.hh"

#include "TChain.h"

#include "Tool.h"
#include "TFile.h"
#include "TVectorT.h"
#include "Stopwatch.h"

class WCSimReader: public Tool {


 public:

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


 private:
  //methods used in Initialise

  bool ReadTree(TChain * chain);
  bool AddTreeToChain(const char * fname, TChain * chain);
  bool CompareTree(TChain * chain, int mode);
  template <typename T> bool CompareVariable(T v1, T v2, const char * tag);

  //methods used in Execute

  SubSample GetDigits();

  TChain * m_chain_opt;
  TChain * m_chain_event;
  TChain * m_chain_geom;

  WCSimRootOptions * m_wcsim_opt;
  WCSimRootEvent   * m_wcsim_event_ID;
  WCSimRootEvent   * m_wcsim_event_OD;
  WCSimRootGeom    * m_wcsim_geom;
  WCSimRootTrigger * m_wcsim_trigger;

  long int m_current_event_num;
  long int m_first_event_num;
  long int m_n_events;

  // interpose an interval between events
  // this allows to convert simulated data, which always start at t = 0,
  // to look like real data in which events are separated over time
  // have to be used with care: dark noise should be generated in WCSim for a given time window, not around physics hits
  // units: ns
  float m_interpose_interval;

  std::string m_input_filename;
  std::string m_input_filelist;

  util::Stopwatch * m_stopwatch;
  std::string m_stopwatch_file;

  int m_verbose;

  std::stringstream m_ss;

  void StreamToLog(int level) {
    Log(m_ss.str(), level, m_verbose);
    m_ss.str("");
  }

  enum LogLevel {FATAL=-1, ERROR=0, WARN=1, INFO=2, DEBUG1=3, DEBUG2=4, DEBUG3=5};
};


#endif