Normalization.hpp

Go to the documentation of this file.

#ifndef NORMALIZATION_HPP
#define NORMALIZATION_HPP
 
#include "Constants.hpp"
#include <Rtypes.h>
#include <TTree.h>

struct EnergyView {
  UShort_t leftdE_adc[16];  
  UShort_t rightdE_adc[16]; 
  UShort_t strip0_adc;      
  UShort_t strip17_adc;     
  UShort_t hits_adc[36];    
  Short_t cathode_adc;      
  Short_t grid_adc;         

  Float_t gain_left[16];  
  Float_t gain_right[16]; 
  Float_t offset_left[16];
  Float_t offset_right[16];
  Float_t gain_strip0;  
  Float_t gain_strip17; 
  Float_t gain_cathode; 
  Float_t gain_grid;
  Float_t strip_factor[18];
  Float_t strip_offset[18];
  Bool_t is_normed;

  Double_t left[16];  
  Double_t right[16]; 
  Double_t strip0;    
  Double_t strip17;   
  Double_t cathode;   
  Double_t grid;      
 
  TTree *tree_;       
  Int_t loaded_tree_; 

  EnergyView()
      : strip0_adc(0), strip17_adc(0), cathode_adc(0), grid_adc(0),
        gain_strip0(0.0f), gain_strip17(0.0f), gain_cathode(0.0f),
        gain_grid(0.0f), is_normed(kFALSE), strip0(0.0), strip17(0.0),
        cathode(0.0), grid(0.0), tree_(nullptr), loaded_tree_(-1) {
    for (Int_t k = 0; k < 16; k++) {
      leftdE_adc[k] = 0;
      rightdE_adc[k] = 0;
      gain_left[k] = 0.0f;
      gain_right[k] = 0.0f;
      offset_left[k] = 0.0f;
      offset_right[k] = 0.0f;
      left[k] = 0.0;
      right[k] = 0.0;
    }
    for (Int_t s = 0; s < 18; s++) {
      strip_factor[s] = 1.0f;
      strip_offset[s] = 0.0f;
    }
  }

  Bool_t Attach(TTree *t);

  void Decode();

  void LoadGains();

  Double_t Total(Int_t strip) const {
    if (strip <= 0)
      return strip0;
    if (strip >= 17)
      return strip17;
    return left[strip - 1] + right[strip - 1];
  }

  Double_t Axis(Int_t axis) const {
    return axis == GATE_AXIS_GRID ? grid : Total(axis);
  }

  void Totals(Double_t *out) const {
    for (Int_t s = 0; s < 18; s++)
      out[s] = Total(s);
  }
};
 
#endif