Page MenuHomeHEPForge
Files F9501567

No OneTemporary
Actions

View Options

diff --git a/src/Core/AnalysisInfo.cc b/src/Core/AnalysisInfo.cc
--- a/src/Core/AnalysisInfo.cc
+++ b/src/Core/AnalysisInfo.cc
@@ -1,194 +1,197 @@
#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/AnalysisInfo.hh"
#include "Rivet/Tools/RivetPaths.hh"
#include "Rivet/Tools/Utils.hh"
#include "Rivet/Tools/Logging.hh"
#include "yaml-cpp/yaml.h"
#include <iostream>
#include <fstream>
#include <unistd.h>
#ifdef YAML_NAMESPACE
#define YAML YAML_NAMESPACE
#endif
namespace Rivet {
namespace {
Log& getLog() {
return Log::getLog("Rivet.AnalysisInfo");
}
}
/// Static factory method
unique_ptr<AnalysisInfo> AnalysisInfo::make(const std::string& ananame) {
// Returned AI, in semi-null state
unique_ptr<AnalysisInfo> ai( new AnalysisInfo );
ai->_beams += make_pair(PID::ANY, PID::ANY);
ai->_name = ananame;
/// If no ana data file found, return null AI
const string datapath = findAnalysisInfoFile(ananame + ".info");
if (datapath.empty()) {
MSG_DEBUG("No datafile " << ananame + ".info found");
return ai;
}
// Read data from YAML document
MSG_DEBUG("Reading analysis data from " << datapath);
YAML::Node doc;
try {
doc = YAML::LoadFile(datapath);
} catch (const YAML::ParserException& ex) {
MSG_ERROR("Parse error when reading analysis data from " << datapath << " (" << ex.what() << ")");
return ai;
}
#define THROW_INFOERR(KEY) throw InfoError("Problem in info parsing while accessing key " + string(KEY) + " in file " + datapath)
// Simple scalars (test for nullness before casting)
#define TRY_GETINFO(KEY, VAR) try { if (doc[KEY] && !doc[KEY].IsNull()) ai->_ ## VAR = doc[KEY].as<string>(); } catch (...) { THROW_INFOERR(KEY); }
- #define TRY_GETINFO_DBL(KEY, VAR) try { if (doc[KEY] && !doc[KEY].IsNull()) ai->_ ## VAR = doc[KEY].as<double>(); } catch (...) { THROW_INFOERR(KEY); }
-
+ #define TRY_GETINFO_DEFAULT(KEY, VAR, DEFAULT) try { if (doc[KEY] && !doc[KEY].IsNull()) ai->_ ## VAR = doc[KEY].as<string>(); } catch (...) { ai->_ ## VAR = DEFAULT; }
+ #define TRY_GETINFO_DBL(KEY, VAR, DEFAULT) try { if (doc[KEY] && !doc[KEY].IsNull()) ai->_ ## VAR = doc[KEY].as<double>(); } catch (...) { THROW_INFOERR(KEY); }
+ #define TRY_GETINFO_DBL_DEFAULT(KEY, VAR, DEFAULT) try { if (doc[KEY] && !doc[KEY].IsNull()) ai->_ ## VAR = doc[KEY].as<double>(); } catch (...) { ai->_ ## VAR = DEFAULT; }
TRY_GETINFO("Name", name);
TRY_GETINFO("Summary", summary);
TRY_GETINFO("Status", status);
TRY_GETINFO("RunInfo", runInfo);
TRY_GETINFO("Description", description);
TRY_GETINFO("Experiment", experiment);
TRY_GETINFO("Collider", collider);
TRY_GETINFO("Year", year);
TRY_GETINFO("SpiresID", spiresId);
TRY_GETINFO("InspireID", inspireId);
TRY_GETINFO("BibKey", bibKey);
TRY_GETINFO("BibTeX", bibTeX);
- TRY_GETINFO_DBL("Luminosity_fb", luminosityfb);
+ TRY_GETINFO_DBL_DEFAULT("Luminosity_fb", luminosityfb, -1);
#undef TRY_GETINFO
+ #undef TRY_GETINFO_DEFAULT
#undef TRY_GETINFO_DBL
+ #undef TRY_GETINFO_DBL_DEFAULT
// Normalise the status info to upper-case
ai->_status = toUpper(ai->_status);
// Sequences (test the seq *and* each entry for nullness before casting)
#define TRY_GETINFO_SEQ(KEY, VAR) try { \
if (doc[KEY] && !doc[KEY].IsNull()) { \
const YAML::Node& VAR = doc[KEY]; \
for (size_t i = 0; i < VAR.size(); ++i) \
if (!VAR[i].IsNull()) ai->_ ## VAR += VAR[i].as<string>(); \
} } catch (...) { THROW_INFOERR(KEY); }
TRY_GETINFO_SEQ("Authors", authors);
TRY_GETINFO_SEQ("References", references);
TRY_GETINFO_SEQ("ToDo", todos);
TRY_GETINFO_SEQ("Keywords", keywords);
TRY_GETINFO_SEQ("Options", options);
TRY_GETINFO_SEQ("ReleaseTests", validation);
#undef TRY_GETINFO_SEQ
// Build the option map
ai->buildOptionMap();
// A boolean with some name flexibility
try {
if (doc["NeedsCrossSection"]) ai->_needsCrossSection = doc["NeedsCrossSection"].as<bool>();
else if (doc["NeedCrossSection"]) ai->_needsCrossSection = doc["NeedCrossSection"].as<bool>();
} catch (...) {
THROW_INFOERR("NeedsCrossSection|NeedCrossSection");
}
// Check if reentrant
if ( ai->statuscheck("REENTRANT") && !ai->statuscheck("NOTREENTRY") )
ai->_reentrant = true;
// Beam particle identities
try {
if (doc["Beams"]) {
const YAML::Node& beams = doc["Beams"];
vector<PdgIdPair> beam_pairs;
if (beams.size() == 2 && beams[0].IsScalar() && beams[0].IsScalar()) {
beam_pairs += PID::make_pdgid_pair(beams[0].as<string>(), beams[1].as<string>());
} else {
for (size_t i = 0; i < beams.size(); ++i) {
const YAML::Node& bp = beams[i];
if (bp.size() != 2 || !bp[0].IsScalar() || !bp[0].IsScalar())
throw InfoError("Beam ID pairs have to be either a 2-tuple or a list of 2-tuples of particle names");
beam_pairs += PID::make_pdgid_pair(bp[0].as<string>(), bp[1].as<string>());
}
}
ai->_beams = beam_pairs;
}
} catch (...) { THROW_INFOERR("Beams"); }
// Beam energies
try {
if (doc["Energies"]) {
vector< pair<double,double> > beam_energy_pairs;
for (size_t i = 0; i < doc["Energies"].size(); ++i) {
const YAML::Node& be = doc["Energies"][i];
if (be.IsScalar()) {
// If beam energy is a scalar, then assume symmetric beams each with half that energy
beam_energy_pairs += make_pair(be.as<double>()/2.0, be.as<double>()/2.0);
} else if (be.IsSequence()) {
if (be.size() != 2)
throw InfoError("Beam energies have to be a list of either numbers or pairs of numbers");
beam_energy_pairs += make_pair(be[0].as<double>(), be[1].as<double>());
} else {
throw InfoError("Beam energies have to be a list of either numbers or pairs of numbers");
}
}
ai->_energies = beam_energy_pairs;
}
} catch (...) { THROW_INFOERR("Energies"); }
#undef THROW_INFOERR
MSG_TRACE("AnalysisInfo pointer = " << ai.get());
return ai;
}
string toString(const AnalysisInfo& ai) {
std::stringstream ss;
ss << ai.name();
ss << " - " << ai.summary();
// ss << " - " << ai.beams();
// ss << " - " << ai.energies();
ss << " (" << ai.status() << ")";
return ss.str();
}
void AnalysisInfo::buildOptionMap() {
_optionmap.clear();
for ( auto opttag : _options ) {
std::vector<std::string> optv = split(opttag, "=");
std::string optname = optv[0];
for ( auto opt : split(optv[1], ",") )
_optionmap[optname].insert(opt);
}
}
bool AnalysisInfo::validOption(std::string key, std::string val) const {
auto opt = _optionmap.find(key);
// The option is required to be defined in the .info file.
if ( opt == _optionmap.end() ) return false;
// If the selection option is among the range of given options,
// we are fine.
if ( opt->second.find(val) != opt->second.end() ) return true;
// Wildcard selection option for value types is #.
if ( opt->second.size() == 1 && *opt->second.begin() == "#" ) {
std::istringstream ss(val);
double test;
if ( ss >> test ) return true;
}
// Wildcard selection option for any type is *.
if ( opt->second.size() == 1 && *opt->second.begin() == "*" )
return true;
return false;
}
}

File Metadata

Mime Type
text/x-diff
Expires
Sun, Feb 23, 2:42 PM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4486682
Default Alt Text
(7 KB)

Event Timeline