PulseHistory.hpp

Go to the documentation of this file.

#ifndef PULSE_HISTORY_HPP
#define PULSE_HISTORY_HPP

#include "BinaryUtils.hpp"
#include <Rtypes.h>
#include <TString.h>
#include <vector>
 
class TH1D;
class TH2D;
class TProfile;
 
namespace PulseHistory {

const UInt_t kFlagClamped = 0x2000;
const Int_t kNBins = 12;
const Double_t kLogLo = -6.0; 
const Double_t kLogHi =
    -4.0; 
const Int_t kMaxAmpBins = 6;

enum Group {
  kNone = 0,       
  kLongLeft = 1,   
  kLongRight = 2,  
  kShortLeft = 3,  
  kShortRight = 4, 
  kGuard0 = 5,     
  kGuard17 = 6,    
  kNGroups = 7     
};

const char *GroupName(Int_t g);
const char *GroupTag(Int_t g);
Int_t ChainOf(Int_t g);
Bool_t IsLongGroup(Int_t g);

struct Kernel {
  Bool_t ok = kFALSE; 
  Int_t n_amp = 1;    
  Double_t k[kMaxAmpBins]
            [kNBins];        
  Double_t intercept = 0.0;  
  Double_t r2 = 0.0;         
  Double_t rms_before = 0.0; 
  Double_t rms_after = 0.0;  
  Long64_t n = 0;            
  Long64_t n_beam = 0; 

  Kernel();
};

struct Result {
  Kernel kernel[kNGroups];    
  Long64_t n_hits = 0;        
  Long64_t n_seeds = 0;       
  Long64_t n_beam_events = 0; 
  Long64_t n_corrected = 0;   
  Long64_t n_clamped = 0; 
  Long64_t n_clamped_group[kNGroups]; 
  Double_t mean_shift[kNGroups];      
  Bool_t sorted_input = kTRUE; 
  std::vector<Double_t> mode;
  TH2D *dev_vs_pred[kNGroups]; 
  TH1D *dev_before[kNGroups];  
  TH1D *dev_after[kNGroups];   
  TH1D *shift[kNGroups];       
  TProfile
      *dtprev_before[kNGroups]; 
  TProfile *dtprev_after[kNGroups]; 
 
  Result();  
  ~Result(); 
};

Int_t BinOf(Double_t dt_s);
Double_t BinCentreUs(Int_t b);
Int_t AmpBinOf(Double_t e_prev, Double_t mode, Int_t n_amp);

std::vector<Int_t> BuildGroupMap();

Bool_t Measure(std::vector<RawHit> &hits, const std::vector<Int_t> &group_of,
               Result &res, const TString &file_label);

void Apply(std::vector<RawHit> &hits, const std::vector<Int_t> &group_of,
           Result &res);

TString Report(const Result &res, const TString &file_label);
void SavePlots(Result &res, const TString &file_label);
void WriteToEventsFile(const TString &events_subpath, const Result &res);
 
} // namespace PulseHistory
 
#endif