bar.hh

Go to the documentation of this file.

 
#ifndef BAR_HH
#define BAR_HH
 
#include <iostream>
#include <fstream>
#include <vector>
#include <cfloat>
#include <string>
#include <regex>
 
#include <TH3D.h>
#include <TRandom3.h>
#include <TTree.h>
#include <TMath.h>
#include <TFile.h>
 
#include "globals.hh"
#include "daq.hh"
 
class BarLYSO
{
public:
    BarLYSO(const char* inputFilename, Int_t threadID);
    ~BarLYSO();
 
    void SetParsDistro();
    void InitializeBaselines(Int_t event);
    void ClearContainers();
    void SetSamplingTimes();
    void SetFrontWaveform(Int_t channel, Double_t start);
    void SetBackWaveform(Int_t channel, Double_t start);
    void SaveEvent();
    void SaveBar();
 
    inline void SetSigmaNoise(Double_t newSigmaNoise) { fSigmaNoise = newSigmaNoise; } 
    inline void SetInputFilename(std::string newInputFilename) { fInputFilename = newInputFilename; } 
    inline void SetChargeCuts(Double_t min, Double_t max)
    {
        fChargeCuts[0] = min;
        fChargeCuts[1] = max;
    }
    inline void SetHisto_A(Double_t nbins, Double_t min, Double_t max)
    {
        fHisto_A[0] = TMath::Nint(nbins);
        fHisto_A[1] = min;
        fHisto_A[2] = max;
    }
    inline void SetHisto_Tau_rise(Double_t nbins, Double_t min, Double_t max)
    {
        fHisto_Tau_rise[0] = TMath::Nint(nbins);
        fHisto_Tau_rise[1] = min;
        fHisto_Tau_rise[2] = max;
    }
    inline void SetHisto_Tau_dec(Double_t nbins, Double_t min, Double_t max)
    {
        fHisto_Tau_dec[0] = TMath::Nint(nbins);
        fHisto_Tau_dec[1] = min;
        fHisto_Tau_dec[2] = max;
    }
 
    inline void SetEvents(Int_t events) { EVENTS = events; } 
    inline Int_t GetEvents() const { return EVENTS; } 
    inline Int_t GetID() const { return fID; } 
    inline DAQ *GetDAQ() const { return fDAQ; }
 
private:
    Int_t fEvent; 
    Int_t fID; 
    Int_t EVENTS;
    Int_t fThreadID;
 
    std::vector<std::vector<Double_t>> fFront; 
    std::vector<std::vector<Double_t>> fBack;  
    std::vector<std::vector<Double_t>> fTimes_F;
    std::vector<std::vector<Double_t>> fTimes_B;
 
    TH3D *hPars; 
    TRandom3 *fRandPars; 
    TRandom3 *fRandNoise; 
    Double_t fSigmaNoise; 
    std::string fInputFilename; 
    Double_t fChargeCuts[2]; 
    Double_t fHisto_A[3]; 
    Double_t fHisto_Tau_rise[3]; 
    Double_t fHisto_Tau_dec[3]; 
    DAQ *fDAQ;
 
    std::string GenerateOutputFilename(const char *inputFilename);
 
    TFile *fOutFile = nullptr;
    TTree *fOutTree = nullptr;
 
    inline Double_t Add_Noise() { return fRandNoise->Gaus(BASELINE, fDAQ->fSigmaNoise); };
    Double_t Wave_OnePhel(Double_t t, Double_t A, Double_t tau_rise, Double_t tau_dec, Double_t timePhel);
 
};
 
#endif  // BAR_HH