The Rivet MC analysis system: Todo List

The Rivet MC analysis system 1.8.2

Member Analysis::addPlot (AnalysisObjectPtr)

2D versions of scale and normalize... or ditch these completely?

Member Analysis::histoDir () const

This doesn't change: calc and cache at first use!

Member Analysis::needsCrossSection () const

Remove this and require HepMC >= 2.06

Member Analysis::setNeedsCrossSection (bool needed=true)

Remove this and require HepMC >= 2.06

Member AnalysisHandler::addAnalysis (const std::string &analysisname)

Might we want to be able to run an analysis twice, with different params? Requires avoiding histo tree clashes, i.e. storing the histos on the analysis objects.

Member AnalysisLoader::_loadAnalysisPlugins ()

Make sure this is an abs path

Sys-dependent path separator instead of "/"

Member ATLAS_2010_S8894728::_moments_to_stddev (Profile1DPtr moment_profiles[], Scatter2DPtr target_dps)

Assuming unit weights here! Should use N_effective = sumW**2/sumW2?

Member ATLAS_2011_I944826::analyze (const Event &event)

YODA

Member ATLAS_2011_I944826::finalize ()

YODA / Division of histograms to obtain lambdabar/lambda ratios

Member ATLAS_2011_I944826::init ()

YODA

Member ATLAS_2011_I945498::ATLAS_2011_I945498 ()

Set whether your finalize method needs the generator cross section

Member ATLAS_2011_I945498::finalize ()

YODA

Class ATLAS_2011_S9041966

Clean up the debug stuff

Member ATLAS_2011_S9126244::finalize ()

YODA divide and binEdges foreach (const ATLAS_2011_S9126244_Plots& plots, m_selectionPlots) { Calculate the gap fraction for each slice for (size_t x=0; x<plots._h_gapVsDeltaYVeto.getHistograms().size(); x++) { histogramFactory().divide(histoPath(makeAxisCode(plots.m_gapFractionDeltaYHistIndex+x, 1, plots.selectionType)), (plots._h_gapVsDeltaYVeto.getHistograms()[x]), (plots._h_gapVsDeltaYInc.getHistograms()[x])); histogramFactory().destroy(plots._h_gapVsDeltaYVeto.getHistograms()[x]); histogramFactory().destroy(plots._h_gapVsDeltaYInc.getHistograms()[x]); }

Member ATLAS_2011_S9126244::m_q0BinEdges

YODA void finalizeQ0GapFraction(double totalWeightSum, Scatter2DPtr gapFractionDP, Histo1DPtr vetoPtHist, int binNumber) { double vetoPtWeightSum = 0.0; for (int x=0; x<binNumber-1; x++) { vetoPtWeightSum += vetoPtHist->bin(x).area();

Member ATLAS_2012_I1083318::finalize ()

YODA / first construct jet multi ratio

Member ATLAS_2012_I1091481::init ()

Initialise and register projections here

Member CDF_1996_S3108457::analyze (const Event &event)

include gaussian jet energy resolution smearing?

Member CDF_1996_S3108457::init ()

Ratios of mass histograms left out: Binning doesn't work out

Member CDF_1996_S3349578::analyze (const Event &event)

include gaussian jet energy resolution smearing?

Member CDF_1996_S3418421::finalize ()

calculate errors while analysing and fill them here as well

Member CDF_2001_S4751469::analyze (const Event &event)

Copy the permanent histos to get the binnings more robustly

Member CDF_2004_S5839831::finalize ()

Take these normalisations from the data histo (it can't come from just the MC)

Member CDF_2006_S6653332::analyze (const Event &event)

Use jet contents rather than accessing quarks directly

Use nicer looping

Use jet contents rather than accessing quarks directly

Use nicer looping

Use jet contents rather than accessing quarks directly

Member CDF_2008_LEADINGJETS::analyze (const Event &e)

Implement Run II min bias trigger cf. CDF_2009?

Member CDF_2008_NOTE_9351::analyze (const Event &e)

Replace with foreach

Replace with PID::isLepton

Member CDF_2008_S7540469::analyze (const Event &event)

Make this neater, using the JetAlg interface and the built-in sorting

Member CDF_2008_S7541902::finalize ()

YODA Replace with autobooking etc. once YODA in place

Class CDF_2009_S8057893

Finish!

Member CDF_2009_S8233977::analyze (const Event &evt)

The pT and sum(ET) distributions look slightly different from Niccolo's Monte Carlo plots. Still waiting for his answer.

Member CDF_2010_S8591881_DY::analyze (const Event &e)

Replace with foreach

Replace with PID::isLepton

Member CDF_2010_S8591881_QCD::analyze (const Event &e)

Implement Run II min bias trigger cf. CDF_2009?

Member CentralEtHCM::project (const Event &e)

Generalise rapidity cut value

Member CMS_2011_S9215166::analyze (const Event &event)

"Diffractive" veto should really also veto dijet events?

Member CMS_2012_I1102908::init ()

Can we manage to only register these as they are "really" created in the finalize()?

Member ConstRandomFilter::operator() (const Particle &)

Use a better RNG

Member D0_1996_S3214044::init ()

Use correct jet algorithm

Member D0_1996_S3324664::finalize ()

Understand this!

Prefer to scale rather than normalize, if possible

Member D0_1996_S3324664::init ()

Use correct jet algorithm

Member D0_2008_S7837160::analyze (const Event &event)

Use pT cut in WFinder

Any ETmiss cut?

Member D0_2008_S7837160::init ()

Use separate pT and ETmiss cuts in WFinder

Member DELPHI_2003_WUD_03_11::finalize ()

Scale rather than normalize?

Member E735_1998_S3905616::init ()

E735 trigger?

Member Event::_geNormAlignment ()

e+ e- convention? B-factories different from LEP?

Member FastJets::_init3 (fastjet::JetDefinition::Plugin *plugin)

Should we be copying the plugin?

Member FastJets::_pseudojetsToJets (const PseudoJets &pjets) const

Cache?

Member FastJets::clusterSeqArea () const

Throw error if no area def? Or just blindly call dynamic_cast?

Member FastJets::reset ()

_cseq = fastjet::ClusterSequence();

Member FinalState::FinalState (const vector< pair< double, double > > &etaRanges, double minpt=0.0 *GeV)

Properly check whether any of these eta ranges (or their combination) are actually open

Member FinalState::project (const Event &e)

In general, we'd like to calculate a restrictive FS based on the most restricted superset FS.

Member FourMomentum::boostVector () const

Be careful about c=1 convention...

Member FoxWolframMoments::FoxWolframMoments (const FinalState &fsp)

Let the user supply any projection they like?

Member FoxWolframMoments::getFoxWolframMoment (unsigned int l) const

What?!?

Member H1_1994_S2919893::analyze (const Event &event)

Use units in w comparisons... what are the units?

Is this profile meant to be filled with 2 weight factors?

Use angle function

Is this profile meant to be filled with 2 weight factors?

Use angle function

Class H1_1995_S3167097

Check this analysis!

Member H1_1995_S3167097::_nbin

Remove statics!

Member H1_1995_S3167097::_nbin

Remove statics!

Member H1_1995_S3167097::_xmax

Remove statics!

Member H1_1995_S3167097::_xmax

Remove statics!

Member H1_1995_S3167097::finalize ()

Automate this sort of thing so that the analysis code is more readable.

Member H1_1995_S3167097::init ()

Automate this sort of thing so that the analysis code is more readable.

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Group Histograms

Add jet masses and d(rap)

Member IS_PARTON_PDGID (id)

Use inline PID functions instead

Member JADE_1998_S3612880::JADE_1998_S3612880 ()

Set whether your finalize method needs the generator cross section

Member JADE_OPAL_2000_S4300807::analyze (const Event &e)

Put this in an index loop?

Member Jet::neutralEnergy () const

Jet::containsMatch(Matcher m) { ... if m(pid) return true; ... }

Member JetAlg::jets (double ptmin=0.0, double ptmax=MAXDOUBLE, double rapmin=-MAXDOUBLE, double rapmax=MAXDOUBLE, RapScheme rapscheme=PSEUDORAPIDITY) const