 |
DAS
3.0
Das Analysis System
|
#include <boost/test/included/unit_test.hpp>
#include <boost/exception/all.hpp>
#include <boost/property_tree/ptree.hpp>
#include <boost/property_tree/info_parser.hpp>
#include <boost/property_tree/json_parser.hpp>
#include <boost/property_tree/xml_parser.hpp>
#include "test.h"
#include "exceptions.h"
#include "Options.h"
#include "UserInfo.h"
#include "MetaInfo.h"
#include "example01.cc"
#include "example02.cc"
#include "example03.cc"
#include "example04.cc"
#include "../bin/getMetaInfo.cc"
#include "../bin/printMetaInfo.cc"
#include "../bin/printEntries.cc"
#include "../bin/printBranches.cc"
#include "../bin/forceMetaInfo.cc"
|
template<bool DARWIN_TEST_EXCEPTIONS = false> |
void | example02 (const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0}) |
|
template<bool DARWIN_TEST_EXCEPTIONS = false> |
void | example03 (const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0}) |
|
void | example04 (const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0}) |
|
| BOOST_AUTO_TEST_CASE (output) |
|
| BOOST_AUTO_TEST_CASE (getters_empty_input) |
|
| BOOST_AUTO_TEST_CASE (getters_broken_input) |
|
| BOOST_AUTO_TEST_CASE (getters_wrong_objects) |
|
| BOOST_AUTO_TEST_CASE (ntuple_producer) |
|
| BOOST_AUTO_TEST_CASE (ntuple_weight_modifier) |
|
| BOOST_AUTO_TEST_CASE (ntuple_JEC_modifier) |
|
| BOOST_AUTO_TEST_CASE (projection) |
|
| BOOST_AUTO_TEST_CASE (all_steps) |
|
| BOOST_AUTO_TEST_CASE (getBool) |
|
| BOOST_AUTO_TEST_CASE (forceMetaInfo) |
|
| BOOST_AUTO_TEST_CASE (GetFirstTreeLocation) |
|
| BOOST_AUTO_TEST_CASE (printEntries) |
|
| BOOST_AUTO_TEST_CASE (printBranches) |
|
| BOOST_AUTO_TEST_CASE (printMetaInfo) |
|
| BOOST_AUTO_TEST_CASE (closure) |
|
| BOOST_AUTO_TEST_CASE (exceptions) |
|
◆ BOOST_TEST_MODULE
#define BOOST_TEST_MODULE test |
◆ DARWIN_GIT_REPO
#define DARWIN_GIT_REPO DARWIN |
◆ DOXYGEN_SHOULD_SKIP_THIS
#define DOXYGEN_SHOULD_SKIP_THIS |
◆ BOOST_AUTO_TEST_CASE() [1/17]
BOOST_AUTO_TEST_CASE |
( |
all_steps |
| ) |
|
276 BOOST_REQUIRE_NO_THROW( pt::read_info(
"example.info",
config) );
279 BOOST_REQUIRE_NO_THROW( DP::example01(
"ntuple01.root",
config, steering) );
280 BOOST_REQUIRE_NO_THROW(
DP::example02(
"ntuple01.root",
"ntuple02.root",
config, steering) );
281 BOOST_REQUIRE_NO_THROW(
DP::example03(
"ntuple02.root",
"ntuple03.root",
config, steering) );
282 BOOST_REQUIRE_NO_THROW(
DP::example04(
"ntuple03.root",
"hists04.root",
config, steering) );
◆ BOOST_AUTO_TEST_CASE() [2/17]
BOOST_AUTO_TEST_CASE |
( |
closure |
| ) |
|
349 unique_ptr<TFile>
f(TFile::Open(
"ntuple02.root"));
353 vector<string> exts {
"info",
"xml",
"json"};
354 for (
string ext: exts) {
355 fs::path file =
"out1." + ext;
356 vector<fs::path>
inputs = {
"ntuple02.root"};
360 auto getHash = [](
const fs::path& ntuple) {
361 auto rootfile = unique_ptr<TFile>(TFile::Open(ntuple.c_str(),
"READ"));
362 auto events = unique_ptr<TTree>(rootfile->Get<TTree>(
"events"));
367 size_t Hash = getHash(
"ntuple02.root");
370 auto chain = [&
config,&getHash](
const fs::path& fOut,
bool isMC) {
371 fs::remove(
"ntuple02.root");
372 fs::remove(
"ntuple03.root");
376 DP::example01(
"ntuple01.root",
config, steering);
380 vector<fs::path>
inputs = {
"hists04.root"};
382 size_t Hash = getHash(
"ntuple02.root");
383 fs::remove(
"ntuple01.root");
384 fs::rename(
"hists04.root", Form(
"hists04_%s_%s.root", &(fOut.extension().c_str()[1]),
isMC ?
"true" :
"false"));
388 pt::read_info(
"out1.info",
config);
389 int oldseed =
config.get<
int>(
"preseed");
390 config.put<
int>(
"preseed", 2*oldseed);
391 BOOST_TEST( Hash != chain(
"out2.info",
true ) );
392 config.put<
int>(
"preseed", oldseed);
393 BOOST_TEST( Hash == chain(
"out2.info",
true ) );
394 BOOST_TEST( Hash != chain(
"out2.info",
false) );
396 pt::read_json(
"out1.json",
config);
397 BOOST_TEST( Hash == chain(
"out2.json",
true ) );
398 BOOST_TEST( Hash != chain(
"out2.json",
false) );
400 pt::read_xml(
"out1.xml",
config);
402 BOOST_TEST( Hash == chain(
"out2.xml",
true ) );
403 BOOST_TEST( Hash != chain(
"out2.xml",
false) );
405 vector<fs::path>
inputs = {
"hists04_info_true.root" ,
"hists04_json_true.root" ,
"hists04_xml_true.root" };
407 inputs = {
"hists04_info_false.root",
"hists04_json_false.root",
"hists04_xml_false.root"};
409 inputs = {
"hists04_info_true.root" ,
"hists04_info_false.root"};
412 for (
string name: {
"true",
"false",
"failure"})
413 fs::remove(
name +
".info");
415 for (
string ext: exts) {
416 fs::remove(
"out1." + ext);
417 fs::remove(
"out2." + ext);
418 for (
string b: {
"true",
"false"})
419 fs::remove(
"hists04_" + ext +
"_" + b +
".root");
◆ BOOST_AUTO_TEST_CASE() [3/17]
BOOST_AUTO_TEST_CASE |
( |
exceptions |
| ) |
|
426 auto f = unique_ptr<TFile>(TFile::Open(
"ntuple02.root",
"READ"));
427 auto t = unique_ptr<TTree>(
f->Get<TTree>(
"events"));
429 cout <<
DE::BadInput(
"Here is a bad input", metainfo).what() << endl;
431 cout <<
DE::AnomalousEvent(
"Here is an anomally (just for the example)", t).what() << endl;
435 BOOST_REQUIRE_NO_THROW( pt::read_info(DARWIN
"/test/example.info",
config) );
438 DP::example01(
"ntuple01.root",
config, steering);
439 BOOST_REQUIRE_THROW( DP::example02<true>(
"ntuple01.root",
"ntuple02.root",
config, steering), wrapexcept<DE::AnomalousEvent> );
440 BOOST_REQUIRE_THROW( DP::example01<true>(
"ntuple01.root",
config, steering), wrapexcept<DE::BadInput > );
442 fs::remove(
"inexistent.out");
443 fs::remove(
"ntuple01.root");
◆ BOOST_AUTO_TEST_CASE() [4/17]
BOOST_AUTO_TEST_CASE |
( |
forceMetaInfo |
| ) |
|
297 vector<fs::path>
input = {
"ntuple02.root"};
300 input = {
"ntuple01.root"};
302 unique_ptr<TFile>
f(TFile::Open(
"ntuple01.root",
"READ"));
303 auto events = unique_ptr<TTree>(
f->Get<TTree>(
"events"));
306 BOOST_TEST( !metainfo.Get<
bool>(
"git",
"reproducible") );
◆ BOOST_AUTO_TEST_CASE() [5/17]
BOOST_AUTO_TEST_CASE |
( |
getBool |
| ) |
|
◆ BOOST_AUTO_TEST_CASE() [6/17]
BOOST_AUTO_TEST_CASE |
( |
GetFirstTreeLocation |
| ) |
|
314 TFile::Open(
"GetFirstTreeLocation.root",
"RECREATE")->Close();
◆ BOOST_AUTO_TEST_CASE() [7/17]
BOOST_AUTO_TEST_CASE |
( |
getters_broken_input |
| ) |
|
89 vector<fs::path>
files{
"/this/dir/does/not/exist"};
◆ BOOST_AUTO_TEST_CASE() [8/17]
BOOST_AUTO_TEST_CASE |
( |
getters_empty_input |
| ) |
|
80 vector<fs::path>
files;
82 BOOST_TEST_MESSAGE(
"Testing exceptions for empty list of input files" );
83 BOOST_REQUIRE_THROW( flow.GetInputHist(
"h") , wrapexcept<runtime_error> );
84 BOOST_REQUIRE_THROW( flow.GetInputTree({1, 0},
"t"), wrapexcept<runtime_error> );
◆ BOOST_AUTO_TEST_CASE() [9/17]
BOOST_AUTO_TEST_CASE |
( |
getters_wrong_objects |
| ) |
|
95 vector<fs::path>
files;
96 for (
int i = 1; i < 4; ++i) {
97 fs::path
name = Form(
"getters_%d.root", i);
100 unique_ptr<TFile>
f(TFile::Open(
name.c_str(),
"RECREATE"));
101 auto h = make_unique<TH1F>(
"h",
"h", 1, 0., 1.);
102 auto t = make_unique<TTree>(
"t",
"t");
103 h->SetDirectory(
f.get());
104 t->SetDirectory(
f.get());
110 BOOST_TEST_MESSAGE(
"Testing exceptions for inexisting histogram" );
111 BOOST_REQUIRE_THROW( flow.GetInputHist(
"thisHistDoesNotExist"), wrapexcept<DE::BadInput> );
112 BOOST_TEST_MESSAGE(
"Testing exceptions for empty tree" );
113 BOOST_REQUIRE_THROW( flow.GetInputTree({1, 0},
"t" ), wrapexcept<invalid_argument> );
114 BOOST_TEST_MESSAGE(
"Testing exceptions for inexisting tree" );
115 BOOST_REQUIRE_THROW( flow.GetInputTree({1, 0},
"thisTreeDoesNotExist"), wrapexcept<DE::BadInput> );
117 BOOST_REQUIRE_NO_THROW( flow.GetInputHist(
"h") );
◆ BOOST_AUTO_TEST_CASE() [10/17]
BOOST_AUTO_TEST_CASE |
( |
ntuple_JEC_modifier |
| ) |
|
193 fs::path
input, output;
198 .output(
"output", &output,
"output ROOT file");
200 for (
string mode: {
"fill",
"Friend"}) {
201 string cmd =
"example03 ntuple02_" + mode +
".root ntuple03_" + mode +
".root -j 3 -k 1 -" + mode.at(0) +
" -s";
208 const auto& slice =
options.slice();
209 const int steering =
options.steering();
214 vector<DP::RecJet> * recs_f =
nullptr,
219 flow_F.GetInputTree();
220 BOOST_REQUIRE_NO_THROW( recs_f = flow_f.GetBranchReadOnly<vector<DP::RecJet>>(
"recJets") );
221 BOOST_REQUIRE_NO_THROW( recs_F = flow_F.GetBranchReadOnly<vector<DP::RecJet>>(
"recJets") );
223 string active_branches_f = flow_f.DumpActiveBranches(),
224 active_branches_F = flow_F.DumpActiveBranches();
225 cout <<
"Without friends: " << active_branches_f << endl;
226 cout <<
"With friends: " << active_branches_F << endl;
227 BOOST_TEST( active_branches_f == active_branches_F );
229 auto tIn_f = flow_f.GetInputTree(),
230 tIn_F = flow_F.GetInputTree();
231 BOOST_TEST( tIn_f->GetEntries() == tIn_F->GetEntries() );
233 for (
long long i = 0; i < tIn_f->GetEntries(); ++i) {
236 BOOST_TEST( recs_f->size() == recs_F->size() );
237 for (
size_t j = 0; j < recs_f->size(); ++j)
238 BOOST_TEST( recs_f->at(j).CorrPt() == recs_F->at(j).CorrPt() );
241 fs::remove(
"ntuple02_fill.root");
242 fs::remove(
"ntuple03_fill.root");
◆ BOOST_AUTO_TEST_CASE() [11/17]
BOOST_AUTO_TEST_CASE |
( |
ntuple_producer |
| ) |
|
129 options.output(
"output" , &output,
"output ROOT file")
130 .arg<
bool> (
"isMC" ,
"flags.isMC" ,
"flag" )
131 .arg<float> (
"R" ,
"flags.R" ,
"R parameter in jet clustering algorithm")
132 .arg<
int > (
"year" ,
"flags.year" ,
"year (20xx)" );
134 auto const args =
tokenize(
"example01 ntuple01.root -c example.info -f");
138 const int steering =
options.steering();
140 BOOST_REQUIRE_NO_THROW( DP::example01(output,
config, steering) );
◆ BOOST_AUTO_TEST_CASE() [12/17]
BOOST_AUTO_TEST_CASE |
( |
ntuple_weight_modifier |
| ) |
|
145 fs::path
input, output;
150 .output(
"output", &output,
"output ROOT file")
151 .arg<fs::path>(
"filters",
"corrections.filters",
"location of list of event filters to apply");
153 for (
string mode: {
"fill",
"Friend"}) {
154 string cmd =
"example02 ntuple01.root ntuple02_" + mode +
".root -c example.info -v -j 2 -k 1 -" + mode.at(0);
160 const auto& slice =
options.slice();
161 const int steering =
options.steering();
173 tIn_F = flow_F.GetInputTree();
174 BOOST_TEST( tIn_f->GetEntries() == tIn_F->GetEntries() );
176 BOOST_TEST_MESSAGE(
"Trying to load an inexistant branch" );
180 BOOST_TEST_MESSAGE(
"Loading the analog branch in the friend and no-friend n-tuples." );
181 BOOST_REQUIRE_NO_THROW( event_f = flow_f.GetBranchReadOnly<
DP::RecEvent>(
"recEvent") );
182 BOOST_REQUIRE_NO_THROW( event_F = flow_F.GetBranchReadOnly<
DP::RecEvent>(
"recEvent") );
184 for (
long long i = 0; i < tIn_f->GetEntries(); ++i) {
◆ BOOST_AUTO_TEST_CASE() [13/17]
BOOST_AUTO_TEST_CASE |
( |
output |
| ) |
|
74 BOOST_REQUIRE_THROW(
DT::GetOutputFile(
"/this/file/does/not/exists.root"), wrapexcept<fs::filesystem_error> );
◆ BOOST_AUTO_TEST_CASE() [14/17]
BOOST_AUTO_TEST_CASE |
( |
printBranches |
| ) |
|
334 vector<fs::path>
inputs {
"ntuple01.root"};
◆ BOOST_AUTO_TEST_CASE() [15/17]
BOOST_AUTO_TEST_CASE |
( |
printEntries |
| ) |
|
320 auto f = TFile::Open(
"printEntries.root",
"RECREATE");
322 auto d =
f->mkdir(
"dir2"); d->cd();
323 auto dd = d->mkdir(
"subdir"); dd->cd();
324 auto t =
new TTree(
"one_tree",
"");
329 fs::remove(
"printEntries.root");
◆ BOOST_AUTO_TEST_CASE() [16/17]
BOOST_AUTO_TEST_CASE |
( |
printMetaInfo |
| ) |
|
340 vector<fs::path>
inputs {
"ntuple01.root"};
342 tree.put<
string>(
"path",
"preseed");
◆ BOOST_AUTO_TEST_CASE() [17/17]
BOOST_AUTO_TEST_CASE |
( |
projection |
| ) |
|
247 fs::path
input, output;
251 .output(
"output", &output,
"output ROOT file");
253 BOOST_TEST( fs::exists(
"ntuple03_Friend.root") );
254 auto const args =
tokenize(
"example04 ntuple03_Friend.root hists04.root -v");
258 const auto& slice =
options.slice();
259 const int steering =
options.steering();
263 fs::remove(
"ntuple01.root");
264 fs::remove(
"ntuple02_Friend.root");
265 fs::remove(
"ntuple03_Friend.root");
266 fs::remove(
"hists04.root");
name
Definition: DYToLL_M-50_13TeV_pythia8_cff_GEN_SIM_RECOBEFMIX_DIGI_L1_DIGI2RAW_L1Reco_RECO.py:48
options
Definition: DYToLL_M-50_13TeV_pythia8_cff_GEN_SIM_RECOBEFMIX_DIGI_L1_DIGI2RAW_L1Reco_RECO.py:41
void example03(const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
Definition: main.cc:51
static bool verbose
Definition: Step.h:40
args
Definition: Ntupliser_cfg.py:11
BOOST_TEST(gendijet.Rapidity()==1.3573785791881385)
Hash from TTree obtained by combination of different elements.
Definition: UserInfo.h:286
f
Definition: Ntupliser_cfg.py:289
void example02(const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
Definition: main.cc:43
void example04(const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
Definition: main.cc:58
files
Definition: compareFiles.py:125
cmd
Definition: Core-cfgcmd.txt:1
config
Definition: Ntupliser_cfg.py:297
input
Definition: DYToLL_M-50_13TeV_pythia8_cff_GEN_SIM_RECOBEFMIX_DIGI_L1_DIGI2RAW_L1Reco_RECO.py:35
Generic detector-level event.
Definition: Event.h:64
string isMC
Definition: Ntupliser_cfg.py:59
def inputs
Definition: jercExample.py:118
Generic exception for problematic event (during event loop).
Definition: exceptions.h:63
Weights weights
e.g. cross section normalisation
Definition: Event.h:26