2 #ifndef RIVET_RivetHandler_HH 3 #define RIVET_RivetHandler_HH 5 #include "Rivet/Config/RivetCommon.hh" 6 #include "Rivet/Particle.hh" 7 #include "Rivet/AnalysisLoader.hh" 8 #include "Rivet/Tools/RivetYODA.hh" 15 typedef std::shared_ptr<Analysis> AnaHandle;
53 double sumW()
const {
return _eventCounter->sumW(); }
55 double sumW2()
const {
return _eventCounter->sumW2(); }
58 const vector<string>&
weightNames()
const {
return _weightNames; }
64 size_t numWeights()
const {
return _weightNames.size(); }
76 void setWeightCap(
const double maxWeight) { _weightCap = maxWeight; }
103 void setCrossSection(
const pair<double, double>& xsec,
bool isUserSupplied=
false);
111 _xs.get()->setActiveWeightIdx(_rivetDefaultWeightIdx);
112 const YODA::Scatter1D::Points& ps = _xs->points();
113 if (ps.size() != 1) {
114 string errMsg =
"value missing when requesting nominal cross-section";
117 double xs = ps[0].x();
118 _xs.get()->unsetActiveWeight();
131 MSG_DEBUG(
"Setting run beams = " << beams <<
" @ " <<
sqrtS()/GeV <<
" GeV");
136 const ParticlePair&
beams()
const {
return _beams; }
144 double sqrtS()
const;
172 std::vector<AnaHandle> rtn;
173 rtn.reserve(_analyses.size());
174 for (
const auto& apair : _analyses) rtn.push_back(apair.second);
179 AnaHandle
analysis(
const std::string& analysisname) {
180 if ( _analyses.find(analysisname) == _analyses.end() )
181 throw LookupError(
"No analysis named '" + analysisname +
"' registered in AnalysisHandler");
183 return _analyses[analysisname];
185 throw LookupError(
"No analysis named '" + analysisname +
"' registered in AnalysisHandler");
224 void init(
const GenEvent& event);
230 void analyze(
const GenEvent& event);
236 void analyze(
const GenEvent* event);
249 void readData(
const std::string& filename);
252 vector<YODA::AnalysisObjectPtr>
getYodaAOs(
bool includeraw=
false)
const;
256 const YODA::AnalysisObjectPtr
getPreload(
string path)
const {
257 auto it = _preloads.find(path);
258 if ( it == _preloads.end() )
return nullptr;
263 void writeData(
const string& filename)
const;
271 dump(dumpfile, period);
279 void dump(
const string& dumpfile,
int period) {
280 _dumpPeriod = period;
281 _dumpFile = dumpfile;
298 void mergeYodas(
const vector<string>& aofiles,
299 const vector<string>& delopts=vector<string>(),
300 const vector<string>& addopts=vector<string>(),
311 enum class Stage { OTHER, INIT, FINALIZE };
328 vector<MultiweightAOPtr> getRivetAOs()
const;
331 void stripOptions(YODA::AnalysisObjectPtr ao,
const vector<string>& delopts)
const;
335 void pushToPersistent();
343 Stage _stage = Stage::OTHER;
346 std::map<std::string, AnaHandle> _analyses;
351 map<string,YODA::AnalysisObjectPtr> _preloads;
354 vector<YODA::AnalysisObjectPtr> _finalizedAOs;
361 std::vector<std::string> _weightNames;
362 std::vector<std::valarray<double> > _subEventWeights;
366 std::vector<size_t> _weightIndices;
369 std::string _runname;
372 mutable CounterPtr _eventCounter;
378 std::pair<double,double> _userxs;
393 std::string _matchWeightNames;
396 std::string _unmatchWeightNames;
399 std::string _nominalWeightName;
413 size_t _defaultWeightIdx;
416 size_t _rivetDefaultWeightIdx;
Definition: MC_Cent_pPb.hh:10
void deselectMultiWeights(std::string patterns="")
Setter for _unmatchWeightNames.
const ParticlePair & beams() const
Get the beam particles for this run, usually determined from the first event.
Definition: AnalysisHandler.hh:136
double nominalCrossSection() const
Get the nominal cross-section.
Definition: AnalysisHandler.hh:110
AnalysisHandler & addAnalysis(Analysis *analysis)
Add an analysis to the run list by object.
bool haveNamedWeights() const
Are any of the weights non-numeric?
vector< YODA::AnalysisObjectPtr > getYodaAOs(bool includeraw=false) const
Get all YODA analysis objects (across all weights, optionally including RAW)
void writeData(const string &filename) const
Write all analyses' plots (via getData) to the named file.
void analyze(const GenEvent &event)
Analyze the given event by reference.
void init(const GenEvent &event)
Initialize a run, with the run beams taken from the example event.
double sumW2() const
Access to the sum of squared-weights.
Definition: AnalysisHandler.hh:55
void setWeightCap(const double maxWeight)
Set the weight cap.
Definition: AnalysisHandler.hh:76
void skipMultiWeights(bool ignore=false)
Setter for _skipWeights.
std::vector< AnaHandle > analyses() const
Get the collection of currently registered analyses.
Definition: AnalysisHandler.hh:171
std::vector< std::string > analysisNames() const
Get a list of the currently registered analyses' names.
void setCrossSection(double xsec, double xsecerr, bool isUserSupplied=false)
Set the cross-section for the process being generated (alternative signature)
Definition: AnalysisHandler.hh:105
Scatter1DPtr crossSection() const
Get the cross-section known to the handler.
Definition: AnalysisHandler.hh:100
PdgIdPair beamIds() const
#define MSG_DEBUG(x)
Debug messaging, not enabled by default, using MSG_LVL.
Definition: Logging.hh:195
void setNoAODump()
Configure the AnalysisObject dump rate and destination.
Definition: AnalysisHandler.hh:274
AnalysisHandler(const string &runname="")
Preferred constructor, with optional run name.
Stage
Definition: AnalysisHandler.hh:311
void readData(const std::string &filename)
Read analysis plots into the histo collection (via addData) from the named file.
void selectMultiWeights(std::string patterns="")
Setter for _matchWeightNames.
Logging system for controlled & formatted writing to stdout.
Definition: Logging.hh:10
void setWeightNames(const GenEvent &ge)
Set the weight names from a GenEvent.
AnalysisHandler & removeAnalyses(const std::vector< std::string > &analysisnames)
Remove analyses from the run list using their names.
void setNominalWeightName(std::string name="")
Setter for _nominalWeightName.
~AnalysisHandler()
The destructor is not virtual, as this class should not be inherited from.
Generic runtime Rivet error.
Definition: Exceptions.hh:12
AnalysisHandler & addAnalyses(const std::vector< std::string > &analysisnames)
Add analyses to the run list using their names.
void mergeYodas(const vector< string > &aofiles, const vector< string > &delopts=vector< string >(), const vector< string > &addopts=vector< string >(), bool equiv=false)
Merge the vector of YODA files, using the cross-section and weight information provided in each...
void dump(const string &dumpfile, int period)
Definition: AnalysisHandler.hh:279
This is the base class of all analysis classes in Rivet.
Definition: Analysis.hh:64
size_t numWeights() const
Indices of the weights in the original weight matrix.
Definition: AnalysisHandler.hh:64
AnalysisHandler & removeAnalysis(const std::string &analysisname)
Remove an analysis from the run list using its name.
void checkBeams(bool check=true)
Definition: AnalysisHandler.hh:148
const vector< string > & weightNames() const
Names of event weight categories.
Definition: AnalysisHandler.hh:58
Stage stage() const
Return the current processing stage.
Definition: AnalysisHandler.hh:314
Error relating to looking up analysis objects in the register.
Definition: Exceptions.hh:61
AnalysisHandler & operator=(const AnalysisHandler &)=delete
The assignment operator is deleted, so it can never be called.
void setNLOSmearing(double frac)
Set the relative width of the NLO smearing window.
Definition: AnalysisHandler.hh:79
The key class for coordination of Analysis objects and the event loop.
Definition: AnalysisHandler.hh:23
const std::map< std::string, AnaHandle > & analysesMap() const
Get the collection of currently registered analyses.
Definition: AnalysisHandler.hh:166
double sumW() const
Access the sum of the event weights seen.
Definition: AnalysisHandler.hh:53
const YODA::AnalysisObjectPtr getPreload(string path) const
Definition: AnalysisHandler.hh:256
string runName() const
Get the name of this run.
void setAODump(const string &dumpfile, int period)
Configure the AnalysisObject dump rate and destination.
Definition: AnalysisHandler.hh:270
void setCrossSection(const pair< double, double > &xsec, bool isUserSupplied=false)
Set the cross-section for the process being generated.
size_t defaultWeightIndex() const
Get the index of the nominal weight-stream.
Definition: AnalysisHandler.hh:73
std::vector< std::string > stdAnalysisNames() const
Get a list of the official analysis names for this release.
void setIgnoreBeams(bool ignore=true)
AnaHandle analysis(const std::string &analysisname)
Get a registered analysis by name.
Definition: AnalysisHandler.hh:179
AnalysisHandler & setRunBeams(const ParticlePair &beams)
Set the beam particles for this run.
Definition: AnalysisHandler.hh:129