CalibrateBeam.hpp

Go to the documentation of this file.

#ifndef CALIBRATE_BEAM_HPP
#define CALIBRATE_BEAM_HPP
 
#include "BeamFit2D.hpp"
#include "Constants.hpp"
#include "FileSet.hpp"
#include "IOUtils.hpp"
#include "InitUtils.hpp"
#include "Normalization.hpp"
#include "Paths.hpp"
#include "PlottingUtils.hpp"
#include <Rtypes.h>
#include <TCanvas.h>
#include <TDirectory.h>
#include <TEllipse.h>
#include <TF1.h>
#include <TF2.h>
#include <TFile.h>
#include <TFitResult.h>
#include <TGraph.h>
#include <TGraphErrors.h>
#include <TH1.h>
#include <TH1D.h>
#include <TH1F.h>
#include <TH2.h>
#include <TH2F.h>
#include <TKey.h>
#include <TLegend.h>
#include <TLine.h>
#include <TList.h>
#include <TMath.h>
#include <TROOT.h>
#include <TSpectrum.h>
#include <TString.h>
#include <TSystem.h>
#include <TTree.h>
#include <algorithm>
#include <cmath>
#include <fstream>
#include <iostream>
#include <map>
#include <mutex>
#include <queue>
#include <set>
#include <thread>
#include <toml++/toml.hpp>
#include <vector>

struct ChannelCal {
  TString name;                 
  Char_t side;                  
  Int_t strip;                  
  Double_t fit_adc = 0.0;       
  Double_t fit_sigma_adc = 0.0; 
  Long64_t n_samples = 0.0;     
  Double_t gain = -1.0;
  Double_t ridge_ratio = 0.0;
};

struct StripAlignmentResult {
  Bool_t ok = kFALSE; 
  Double_t beam_e_min =
      0.0; 
  Double_t beam_e_max = 0.0; 
  Double_t factors[18] = {};
  Double_t centroids[18] = {};
};

class CalibrateBeam {
public:
  static std::vector<ChannelCal> BuildChannels();

  static void
  CalibrateBeamOneSubfile(const FileSpec &spec,
                          const std::vector<ChannelCal> &chans_template);

  static void
  SaveCalibSampleOverlay(const std::vector<ChannelCal> &chans,
                         const std::vector<std::vector<Float_t>> &samples,
                         const TString &plot_subdir, const TString &file_label);

  static void AggregateRidgeRatiosForRun(Int_t run,
                                         const std::vector<FileSpec> &specs);
  static void AggregateEresTomlForRun(Int_t run,
                                      const std::vector<FileSpec> &specs);

  static void Run(const TString &file_label = "");
};
 
#endif