Skip to content
Snippets Groups Projects

Analysis: dis_analysis_ecal, dis_analysis_hcal, dis_analysis_raw

Merged Wouter Deconinck requested to merge dis_analysis_ecal_hcal_raw into master
4 files
+ 634
3
Compare changes
  • Side-by-side
  • Inline
Files
4
+ 60
0
 
#include "common_bench/benchmark.h"
 
#include "common_bench/mt.h"
 
#include "common_bench/util.h"
 
#include "common_bench/plot.h"
 
 
#include <cmath>
 
#include <fstream>
 
#include <iostream>
 
#include <string>
 
#include <vector>
 
#include <algorithm>
 
#include <utility>
 
 
#include "ROOT/RDataFrame.hxx"
 
#include <TH1D.h>
 
#include <TFitResult.h>
 
#include <TRandom3.h>
 
#include <TCanvas.h>
 
 
#include "fmt/color.h"
 
#include "fmt/core.h"
 
 
#include "nlohmann/json.hpp"
 
#include "eicd/RawCalorimeterHit.h"
 
 
int dis_electrons_ecal(const std::string& config_name)
 
{
 
// read our configuration
 
std::ifstream config_file{config_name};
 
nlohmann::json config;
 
config_file >> config;
 
 
const std::string rec_file = config["rec_file"];
 
const std::string detector = config["detector"];
 
const std::string output_prefix = config["output_prefix"];
 
const std::string test_tag = config["test_tag"];
 
 
fmt::print(fmt::emphasis::bold | fg(fmt::color::forest_green),
 
"Running DIS electron analysis...\n");
 
fmt::print(" - Detector package: {}\n", detector);
 
fmt::print(" - input file: {}\n", rec_file);
 
fmt::print(" - output prefix: {}\n", output_prefix);
 
fmt::print(" - test tag: {}\n", test_tag);
 
 
// create our test definition
 
// test_tag
 
common_bench::Test dis_Q2_resolution{
 
{{"name", fmt::format("{}_Q2_resolution", test_tag)},
 
{"title", "DIS Q2 resolution"},
 
{"description",
 
fmt::format("DIS Q2 resolution with {}, estimated using a Gaussian fit.", detector)},
 
{"quantity", "resolution (in %)"},
 
{"target", "0.1"}}};
 
 
// Run this in multi-threaded mode if desired
 
ROOT::EnableImplicitMT(kNumThreads);
 
ROOT::RDataFrame d("events", rec_file);
 
 
return 0;
 
}
Loading