Skip to content
Snippets Groups Projects
Commit 23dfa62e authored by Sylvester Joosten's avatar Sylvester Joosten
Browse files

Improved replay scripts so they support different modes linked to def-files... 

Improved replay scripts so they support different modes linked to def-files defined in a database file
parent 7ccac04e
Branches
Tags
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
DEF-files/definitions.json 0 → 100644
+ 43
− 0
View file @ 23dfa62e
{
"shms": {
"default": {
"odef": "DEF-files/SHMS/PRODUCTION/pstackana_production.def",
"cuts": "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_production_cuts.def"
},
"all": {
"odef": "DEF-files/SHMS/PRODUCTION/pstackana_production_all.def",
"cuts": "DEF-files/SHMS/PRODUCTION/CUTS/pstackana_production_cuts.def"
},
"scalers": {
"odef": "DEF-files/SHMS/EPICS/epics_short.def",
"cuts": "DEF-files/SHMS/EPICS/pscaler_cuts.def"
},
"raster": {
"odef": "DEF-files/SHMS/RASTER/shms_raster_simple.def",
"cuts": "DEF-files/SHMS/RASTER/shms_raster_simple_cuts.def"
}
},
"hms": {
"default": {
"odef": "DEF-files/HMS/PRODUCTION/hstackana_production.def",
"cuts": "DEF-files/HMS/PRODUCTION/CUTS/hstackana_production_cuts.def"
},
"all": {
"odef": "DEF-files/HMS/PRODUCTION/hstackana_production_all.def",
"cuts": "DEF-files/HMS/PRODUCTION/CUTS/hstackana_production_cuts.def"
},
"scalers": {
"odef": "DEF-files/HMS/EPICS/epics_short.def",
"cuts": "DEF-files/HMS/EPICS/pscaler_cuts.def"
},
"raster": {
"odef": "DEF-files/HMS/RASTER/hms_raster_simple.def",
"cuts": "DEF-files/HMS/RASTER/hms_raster_simple_cuts.def"
}
}
SCRIPTS/include/nlohmann/json.hpp 0 → 100644
+ 20891
− 0
View file @ 23dfa62e
This diff is collapsed.
SCRIPTS/PRODUCTION/replay_hms.cxx → SCRIPTS/src/replay_hms.cxx
+ 32
− 12
View file @ 23dfa62e
......@@ -9,6 +9,8 @@ using namespace std;
// Better formatting
#include <fmt/format.h>
// json file parsing
#include "nlohmann/json.hpp"
// analyzer includes
#include "THaCut.h"
......@@ -41,18 +43,34 @@ std::string coda_file_pattern(bool do_coin) {
return fmt::format("{}_all_{{:05d}}.dat", do_coin ? "coin" : "hms");
}
std::string output_file_pattern(string_view path, string_view content, string_view extension,
bool do_coin, bool do_all) {
return fmt::format("{}/hms{}_{}{}_{{:05d}}_{{}}.{}", path, do_coin ? "_coin" : "", content,
do_all ? "_all" : "", extension);
const std::string& mode, const bool do_coin) {
return fmt::format("{}/hms{}_{}_{}_{{:05d}}_{{}}.{}", path, do_coin ? "_coin" : "", content, mode,
extension);
}
int replay_hms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0,
const bool do_all = false, const bool do_coin = false) {
const std::string mode = "default", const bool do_coin = false) {
// ===========================================================================
// Setup logging
spdlog::set_level(spdlog::level::warn);
spdlog::flush_every(std::chrono::seconds(5));
// ==========================================================================================
// Load the DEF-file info database
// ==========================================================================================
using nlohmann::json;
nlohmann::json def_db;
{
std::ifstream json_input_file("DEF-files/definitions.json");
try {
json_input_file >> def_db;
} catch (json::parse_error) {
std::cerr << "error: json file, " << run_list_fname
<< ", is incomplete or has broken syntax.\n";
return -127;
}
}
// ===========================================================================
// Get RunNumber and MaxEvent if not provided.
if (RunNumber == 0) {
......@@ -78,15 +96,13 @@ int replay_hms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0
vector<TString> pathList;
pathList.push_back(".");
pathList.push_back("./DATA/raw");
pathList.push_back("./raw_coda");
pathList.push_back("./raw");
pathList.push_back("./raw.copiedtotape");
pathList.push_back("./raw/../raw.copiedtotape");
pathList.push_back("./cache");
// 2. Output files
const auto ROOTFileNamePattern =
output_file_pattern("ROOTfiles", "replay_production", "root", do_coin, do_all);
output_file_pattern("ROOTfiles", "replay_production", "root", mode, do_coin);
// Load global parameters
gHcParms->Define("gen_run_number", "Run Number", RunNumber);
......@@ -201,6 +217,11 @@ int replay_hms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0
// Analyzer
// -----------------------------------------------------------
//
// Set up the analyzer - we use the standard one,
// but this could be an experiment-specific one as well.
// The Analyzer controls the reading of the data, executes
// tests/cuts, loops over Acpparatus's and PhysicsModules,
// and executes the output routines.
THcAnalyzer* analyzer = new THcAnalyzer;
analyzer->SetCodaVersion(2);
// analyzer->EnableBenchmarks(true);
......@@ -239,13 +260,12 @@ int replay_hms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0
analyzer->SetOutFile(ROOTFileName.c_str());
// Define DEF-file+
// analyzer->SetOdefFile("DEF-files/HMS/PRODUCTION/pstackana_production_all.def");
analyzer->SetOdefFile(do_all ? "DEF-files/HMS/PRODUCTION/hstackana_production_all.def"
: "DEF-files/HMS/PRODUCTION/hstackana_production.def");
analyzer->SetOdefFile(def_db["hms"][mode]["odef"].c_str());
// Define cuts file
analyzer->SetCutFile("DEF-files/HMS/PRODUCTION/CUTS/hstackana_production_cuts.def"); // optional
analyzer->SetCutFile(def_db["hms"][mode]["cut"].c_str()); // optional
// File to record accounting information for cuts
analyzer->SetSummaryFile(fmt::format(output_file_pattern("REPORT_OUTPUT/PRODUCTION", "summary",
"report", do_coin, do_all),
"report", mode, do_coin),
RunNumber, MaxEvent)
.c_str());
// Start the actual analysis.
......@@ -254,7 +274,7 @@ int replay_hms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0
analyzer->PrintReport(
"TEMPLATES/HMS/PRODUCTION/hstackana_production.template",
fmt::format(output_file_pattern("REPORT_OUTPUT/PRODUCTION", "replay_production", "report",
do_coin, do_all),
do_coin, mode),
RunNumber, MaxEvent)
.c_str());
......
SCRIPTS/PRODUCTION/replay_shms.cxx → SCRIPTS/src/replay_shms.cxx
+ 32
− 12
View file @ 23dfa62e
......@@ -9,6 +9,8 @@ using namespace std;
// Better formatting
#include <fmt/format.h>
// json file parsing
#include "nlohmann/json.hpp"
// analyzer includes
#include "THaCut.h"
......@@ -41,18 +43,34 @@ std::string coda_file_pattern(bool do_coin) {
return fmt::format("{}_all_{{:05d}}.dat", do_coin ? "coin" : "shms");
}
std::string output_file_pattern(string_view path, string_view content, string_view extension,
bool do_coin, bool do_all) {
return fmt::format("{}/shms{}_{}{}_{{:05d}}_{{}}.{}", path, do_coin ? "_coin" : "", content,
do_all ? "_all" : "", extension);
const std::string& mode, const bool do_coin) {
return fmt::format("{}/shms{}_{}_{}_{{:05d}}_{{}}.{}", path, do_coin ? "_coin" : "", content,
mode, extension);
}
int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent = 0,
const bool do_all = false, const bool do_coin = false) {
const std::string mode = "default", const bool do_coin = false) {
// ===========================================================================
// Setup logging
spdlog::set_level(spdlog::level::warn);
spdlog::flush_every(std::chrono::seconds(5));
// ==========================================================================================
// Load the DEF-file info database
// ==========================================================================================
using nlohmann::json;
nlohmann::json def_db;
{
std::ifstream json_input_file("DEF-files/definitions.json");
try {
json_input_file >> def_db;
} catch (json::parse_error) {
std::cerr << "error: json file, " << run_list_fname
<< ", is incomplete or has broken syntax.\n";
return -127;
}
}
// ===========================================================================
// Get RunNumber and MaxEvent if not provided.
if (RunNumber == 0) {
......@@ -78,15 +96,13 @@ int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent =
vector<TString> pathList;
pathList.push_back(".");
pathList.push_back("./DATA/raw");
pathList.push_back("./raw_coda");
pathList.push_back("./raw");
pathList.push_back("./raw.copiedtotape");
pathList.push_back("./raw/../raw.copiedtotape");
pathList.push_back("./cache");
// 2. Output files
const auto ROOTFileNamePattern =
output_file_pattern("ROOTfiles", "replay_production", "root", do_coin, do_all);
output_file_pattern("ROOTfiles", "replay_production", "root", mode, do_coin);
// Load global parameters
gHcParms->Define("gen_run_number", "Run Number", RunNumber);
......@@ -155,7 +171,7 @@ int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent =
THcTrigDet* shms = new THcTrigDet("shms", "SHMS Trigger Information");
shms->SetSpectName("P");
TRG->AddDetector(shms);
THcHelicity* helicity = new THcHelicity("helicity","Helicity Detector");
THcHelicity* helicity = new THcHelicity("helicity", "Helicity Detector");
TRG->AddDetector(helicity);
// ===========================================================================
......@@ -207,6 +223,11 @@ int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent =
// Analyzer
// -----------------------------------------------------------
//
// Set up the analyzer - we use the standard one,
// but this could be an experiment-specific one as well.
// The Analyzer controls the reading of the data, executes
// tests/cuts, loops over Acpparatus's and PhysicsModules,
// and executes the output routines.
THcAnalyzer* analyzer = new THcAnalyzer;
analyzer->SetCodaVersion(2);
// analyzer->EnableBenchmarks(true);
......@@ -245,10 +266,9 @@ int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent =
analyzer->SetOutFile(ROOTFileName.c_str());
// Define DEF-file+
// analyzer->SetOdefFile("DEF-files/SHMS/PRODUCTION/pstackana_production_all.def");
analyzer->SetOdefFile(do_all ? "DEF-files/SHMS/PRODUCTION/pstackana_production_all.def"
: "DEF-files/SHMS/PRODUCTION/pstackana_production.def");
analyzer->SetOdefFile(def_db["shms"][mode]["odef"].c_str());
// Define cuts file
analyzer->SetCutFile("DEF-files/SHMS/PRODUCTION/CUTS/pstackana_production_cuts.def"); // optional
analyzer->SetCutFile(def_db["shms"][mode]["cut"].c_str()); // optional
// File to record accounting information for cuts
analyzer->SetSummaryFile(fmt::format(output_file_pattern("REPORT_OUTPUT/PRODUCTION", "summary",
"report", do_coin, do_all),
......@@ -260,7 +280,7 @@ int replay_shms(Int_t RunNumber = 7160, Int_t MaxEvent = -1, Int_t FirstEvent =
analyzer->PrintReport(
"TEMPLATES/SHMS/PRODUCTION/pstackana_production.template",
fmt::format(output_file_pattern("REPORT_OUTPUT/PRODUCTION", "replay_production", "report",
do_coin, do_all),
mode, do_coin),
RunNumber, MaxEvent)
.c_str());
......
rootlogon.C
+ 1
− 2
View file @ 23dfa62e
void rootlogon() {
gSystem->SetBuildDir("$HOME/.root_build_dir");
gROOT->ProcessLine(".include /usr/local/include/hcana");
gROOT->ProcessLine(".include SCRIPTS/include");
gROOT->ProcessLine(".include .");
R__LOAD_LIBRARY(libfmt)
//gPrintViaErrorHandler = kTRUE;
//gErrorIgnoreLevel = kWarning;
gROOT->SetStyle("Plain");
gStyle->SetPalette(1);
gSystem->Load("libTree");
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment