EventData.hh

Go to the documentation of this file.

#ifndef EVENTDATA_h
#define EVENTDATA_h

#include <vector>
#include <string>
#include <stdint.h>
#include <iomanip>

#include "Rtypes.h" //has the classdef macro
#include "ChannelData.hh"
#include "SumOfIntegral.hh"


//notice: members with comment starting with ! are not saved

class EventData{
  //interface members
public:
  EventData() { Clear(); }
  virtual ~EventData() {} //anything need cleaning up?
  void Clear(); //inlined below
  static const char* GetBranchName(){ return "event"; }
  void Print (int verbosity);
public:
  
  enum STATUS_FLAGS { NORMAL=0, ID_MISMATCH=1, /*enter more here*/};
  
  //data members
  int run_id; 
  int event_id; 
  uint64_t status;             
  int trigger_count; 
  long timestamp; 
  uint64_t dt; 
  uint64_t event_time; 
  int nchans; 
  bool saturated; 
  bool pulses_aligned; 
  std::vector<double> generic; 
  std::vector<ChannelData> channels;  
  std::vector<SumOfIntegral> sum_of_int;  
  std::vector<Roi> roi_sum_of_int;       
  
  // Energy parameters
  bool s1_valid;  
  bool s1_fixed_valid; 
  bool s2_valid;  
  bool s2_fixed_valid; 
  bool s1s2_valid; 
  bool s1s2_fixed_valid; 
  double s1_start_time; 
  double s1_end_time; 
  double s2_start_time; 
  double s2_end_time; 
  
  double drift_time; 
  
  double s1_full;   
  double s2_full;   
  double s1_fixed;   
  double s2_fixed;   
  double max_s1;    
  double max_s2;    
  int max_s1_chan; 
  int max_s2_chan; 
  double f90_full;  
  double f90_fixed; 
  double gatti; 
  double ll_r; 
    
  //position  
  bool position_valid; 
  double x;          
  double y;          
  double z;          
  double bary_valid; 
  double bary_x;     
  double bary_y;     
    
  ChannelData* GetChannelByID(int id){
    std::vector<ChannelData>::iterator it = channels.begin();
    while(it != channels.end() && it->channel_id != id) it++;
    return (it == channels.end() ? 0 : &(*it) );
  }
  
  Pulse* GetPulse(size_t pulsenum, int channel_id = ChannelData::CH_SUM){ 
    ChannelData* ch = GetChannelByID(channel_id);
    if(ch && ch->pulses.size() > pulsenum)
      return &(ch->pulses[pulsenum]);
    return 0;
  }
  Roi* GetROI(size_t region, int channel_id = ChannelData::CH_SUM){
    ChannelData* ch = GetChannelByID(channel_id);
    if(ch && ch->regions.size() > region)
      return &(ch->regions[region]);
    return 0;
  }
  
  ClassDef(EventData,14)
};



inline void EventData::Clear()
{
  run_id = -1;
  event_id = -1;
  status = NORMAL;
  trigger_count = -1;
  timestamp = 0;
  dt = 0;
  event_time = 0;
  nchans = -1;
  saturated = false;
  pulses_aligned = false;
  generic.clear();
  s1_valid = false;
  s2_valid = false;
  s1s2_valid = false;
  s1_start_time = 0;
  s1_end_time = 0;
  s2_start_time = 0;
  s2_end_time = 0;
  s1_full = 0;
  s2_full = 0;
  s1_fixed = 0;
  s2_fixed = 0;
  max_s1 = 0;
  max_s2 = 0;
  max_s1_chan = -1;
  max_s2_chan = -1;
  f90_full = 0;
  f90_fixed = 0;
  gatti = 0;
  ll_r = 0;
  drift_time = 0;
  position_valid = false;
  x=0; 
  y=0;
  z=0;
  bary_valid = false;
  bary_x = 0;
  bary_y = 0;
  channels.clear();
}
#endif