MC_SUSY Class Reference

MC validation analysis for SUSY events. More...

Inheritance diagram for MC_SUSY:

[legend]

Collaboration diagram for MC_SUSY:

[legend]

List of all members.

Public Member Functions
	MC_SUSY ()
	Constructor.
AnalysisHandler &	handler () const
	Access the controlling AnalysisHandler object.
void	normalize (AIDA::IHistogram1D *&histo, double norm=1.0)
void	normalize (AIDA::IHistogram2D *&histo, double norm=1.0)
void	scale (AIDA::IHistogram1D *&histo, double scale)
void	scale (AIDA::IHistogram2D *&histo, double scale)
Analysis &	setCrossSection (double xs)
	Set the cross section from the generator.
Analysis methods
void	init ()
void	analyze (const Event &evt)
void	finalize ()
Metadata
Metadata is used for querying from the command line and also for building web pages and the analysis pages in the Rivet manual.
const AnalysisInfo &	info () const
	Get the actual AnalysisInfo object in which all this metadata is stored.
virtual std::string	name () const
	Get the name of the analysis.
virtual std::string	inspireId () const
	Get the Inspire ID code for this analysis.
virtual std::string	spiresId () const
	Get the SPIRES ID code for this analysis (~deprecated).
virtual std::vector< std::string >	authors () const
	Names & emails of paper/analysis authors.
virtual std::string	summary () const
	Get a short description of the analysis.
virtual std::string	description () const
	Get a full description of the analysis.
virtual std::string	runInfo () const
	Information about the events needed as input for this analysis.
virtual std::string	experiment () const
	Experiment which performed and published this analysis.
virtual std::string	collider () const
	Collider on which the experiment ran.
virtual std::string	year () const
	When the original experimental analysis was published.
virtual std::vector< std::string >	references () const
	Journal, and preprint references.
virtual std::string	bibKey () const
	BibTeX citation key for this article.
virtual std::string	bibTeX () const
	BibTeX citation entry for this article.
virtual std::string	status () const
	Whether this analysis is trusted (in any way!).
virtual std::vector< std::string >	todos () const
	Any work to be done on this analysis.
virtual const std::vector < PdgIdPair > &	requiredBeams () const
	Return the allowed pairs of incoming beams required by this analysis.
virtual Analysis &	setRequiredBeams (const std::vector< PdgIdPair > &requiredBeams)
	Declare the allowed pairs of incoming beams required by this analysis.
virtual const std::vector < std::pair< double, double > > &	requiredEnergies () const
	Sets of valid beam energy pairs, in GeV.
virtual Analysis &	setRequiredEnergies (const std::vector< std::pair< double, double > > &requiredEnergies)
	Declare the list of valid beam energy pairs, in GeV.
bool	needsCrossSection () const
	Return true if this analysis needs to know the process cross-section.
Analysis &	setNeedsCrossSection (bool needed=true)
	Declare whether this analysis needs to know the process cross-section from the generator.
Internal metadata modifiying methods
AnalysisInfo &	info ()
	Get the actual AnalysisInfo object in which all this metadata is stored (non-const).
virtual Analysis &	setBeams (PdgId beam1, PdgId beam2)
Run conditions
const ParticlePair &	beams () const
	Incoming beams for this run.
const PdgIdPair	beamIds () const
	Incoming beam IDs for this run.
double	sqrtS () const
	Centre of mass energy for this run.
Analysis / beam compatibility testing
bool	isCompatible (const ParticlePair &beams) const
	Check if analysis is compatible with the provided beam particle IDs and energies.
bool	isCompatible (PdgId beam1, PdgId beam2, double e1, double e2) const
	Check if analysis is compatible with the provided beam particle IDs and energies.
bool	isCompatible (const PdgIdPair &beams, const std::pair< double, double > &energies) const
	Check if analysis is compatible with the provided beam particle IDs and energies.
Projection "getting" functions
std::set< ConstProjectionPtr >	getProjections () const
	Get the contained projections, including recursion.
template<typename PROJ >
const PROJ &	getProjection (const std::string &name) const
	Get the named projection, specifying return type via a template argument.
const Projection &	getProjection (const std::string &name) const
Projection applying functions
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const PROJ &proj) const
	Apply the supplied projection on event.
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const Projection &proj) const
	Apply the supplied projection on event.
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const std::string &name) const
	Apply the named projection on event.
Protected Member Functions
Log &	getLog () const
	Get a Log object based on the name() property of the calling analysis object.
double	crossSection () const
	Get the process cross-section in pb. Throws if this hasn't been set.
double	crossSectionPerEvent () const
size_t	numEvents () const
double	sumOfWeights () const
ProjectionHandler &	getProjHandler () const
	Get a reference to the ProjectionHandler for this thread.
AIDA analysis infrastructure.
AIDA::IAnalysisFactory &	analysisFactory ()
	Access the AIDA analysis factory of the controlling AnalysisHandler object.
AIDA::ITree &	tree ()
	Access the AIDA tree of the controlling AnalysisHandler object.
AIDA::IHistogramFactory &	histogramFactory ()
	Access the AIDA histogram factory of the controlling AnalysisHandler object.
AIDA::IDataPointSetFactory &	datapointsetFactory ()
	Access the AIDA histogram factory of the controlling AnalysisHandler object.
const std::string	histoDir () const
	Get the canonical histogram "directory" path for this analysis.
const std::string	histoPath (const std::string &hname) const
	Get the canonical histogram path for the named histogram in this analysis.
const std::string	makeAxisCode (size_t datasetId, size_t xAxisId, size_t yAxisId) const
	Get the internal histogram name for given d, x and y (cf. HepData).
Internal histogram booking (for use by Analysis sub-classes).
const BinEdges &	binEdges (const std::string &hname) const
	Get bin edges for a named histo (using ref AIDA caching).
const BinEdges &	binEdges (size_t datasetId, size_t xAxisId, size_t yAxisId) const
	Get bin edges for a numbered histo (using ref AIDA caching).
BinEdges	logBinEdges (size_t nbins, double lower, double upper)
	Get bin edges with logarithmic widths.
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram2D *	bookHistogram2D (const std::string &name, size_t nxbins, double xlower, double xupper, size_t nybins, double ylower, double yupper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="", const std::string &ztitle="")
AIDA::IHistogram2D *	bookHistogram2D (const std::string &name, const std::vector< double > &xbinedges, const std::vector< double > &ybinedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="", const std::string &ztitle="")
Internal profile histogram booking (for use by Analysis sub-classes).
AIDA::IProfile1D *	bookProfile1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D *	bookProfile1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D *	bookProfile1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D *	bookProfile1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
Internal data point set booking (for use by Analysis sub-classes).
AIDA::IDataPointSet *	bookDataPointSet (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IDataPointSet *	bookDataPointSet (const std::string &name, size_t npts, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IDataPointSet *	bookDataPointSet (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
Projection registration functions
template<typename PROJ >
const PROJ &	addProjection (const PROJ &proj, const std::string &name)
const Projection &	_addProjection (const Projection &proj, const std::string &name)
	Untemplated function to do the work...
Protected Attributes
string	_defaultname
	Name passed to constructor (used to find .info analysis data file, and as a fallback).
shared_ptr< AnalysisInfo >	_info
	Pointer to analysis metadata object.
bool	_allowProjReg
	Flag to forbid projection registration in analyses until the init phase.
Private Attributes
AIDA::IHistogram1D *	_hist_n_trk
AIDA::IHistogram1D *	_hist_phi_trk
AIDA::IHistogram1D *	_hist_eta_trk
AIDA::IHistogram1D *	_hist_pt_trk
AIDA::IHistogram1D *	_hist_n_jet
AIDA::IHistogram1D *	_hist_phi_jet
AIDA::IHistogram1D *	_hist_eta_jet
AIDA::IHistogram1D *	_hist_pt_jet
AIDA::IHistogram1D *	_hist_n_e
AIDA::IHistogram1D *	_hist_phi_e
AIDA::IHistogram1D *	_hist_eta_e
AIDA::IHistogram1D *	_hist_pt_e
AIDA::IHistogram1D *	_hist_n_mu
AIDA::IHistogram1D *	_hist_phi_mu
AIDA::IHistogram1D *	_hist_eta_mu
AIDA::IHistogram1D *	_hist_pt_mu
AIDA::IHistogram1D *	_hist_n_gamma
AIDA::IHistogram1D *	_hist_phi_gamma
AIDA::IHistogram1D *	_hist_eta_gamma
AIDA::IHistogram1D *	_hist_pt_gamma
AIDA::IHistogram1D *	_hist_n_gammaiso
AIDA::IHistogram1D *	_hist_phi_gammaiso
AIDA::IHistogram1D *	_hist_eta_gammaiso
AIDA::IHistogram1D *	_hist_pt_gammaiso
AIDA::IHistogram1D *	_hist_met
AIDA::IHistogram1D *	_hist_mll_2_ossf_ee
AIDA::IHistogram1D *	_hist_mll_2_ossf_mumu
AIDA::IHistogram1D *	_hist_mll_2_osof_emu
AIDA::IHistogram1D *	_hist_mll_ossf_ee
AIDA::IHistogram1D *	_hist_mll_ossf_mumu
AIDA::IHistogram1D *	_hist_mll_osof_emu
AIDA::IHistogram1D *	_hist_mll_all_ossf_ee
AIDA::IHistogram1D *	_hist_mll_all_ossf_mumu
AIDA::IHistogram1D *	_hist_mll_all_osof_emu

---

Detailed Description

MC validation analysis for SUSY events.

Author:

Andy Buckley

Definition at line 18 of file MC_SUSY.cc.

---

Constructor & Destructor Documentation

MC_SUSY

(

)

[inline]

Constructor.

Definition at line 22 of file MC_SUSY.cc.

      : Analysis("MC_SUSY")
    {    }

---

Member Function Documentation

const Projection & _addProjection	(	const Projection &	proj,
		const std::string &	name
	)			[protected, inherited]

Untemplated function to do the work...

Definition at line 33 of file ProjectionApplier.cc.

References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().

Referenced by ProjectionApplier::addProjection().

                                                                             {
    if (!_allowProjReg) {
      cerr << "Trying to register projection '"
           << proj.name() << "' before init phase in '" << this->name() << "'." << endl;
      exit(2);
    }
    const Projection& reg = getProjHandler().registerProjection(*this, proj, name);
    return reg;
  }

const PROJ& addProjection	(	const PROJ &	proj,
		const std::string &	name
	)			[inline, protected, inherited]

Register a contained projection. The type of the argument is used to instantiate a new projection internally: this new object is applied to events rather than the argument object. Hence you are advised to only use locally-scoped Projection objects in your Projection and Analysis constructors, and to avoid polymorphism (e.g. handling ConcreteProjection via a pointer or reference to type Projection) since this will screw up the internal type management.

Definition at line 113 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

Referenced by ZFinder::_init(), WFinder::_init(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ClusteredPhotons::ClusteredPhotons(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), OPAL_1993_S2692198::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_VH2BB::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_LEADINGJETS::init(), MC_JETS::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), LHCB_2010_S8758301::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2003_WUD_03_11::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), DELPHI_1995_S3137023::init(), D0_2010_S8821313::init(), D0_2010_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CMS_2011_S9120041::init(), CMS_2011_S9088458::init(), CMS_2011_S9086218::init(), CMS_2011_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::init(), CMS_2011_S8950903::init(), CMS_2011_S8884919::init(), CMS_2010_S8656010::init(), CMS_2010_S8547297::init(), CDF_2010_S8591881_QCD::init(), CDF_2010_S8591881_DY::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6217184::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2011_S9225137::init(), ATLAS_2011_S9212183::init(), ATLAS_2011_S9131140::init(), ATLAS_2011_S9128077::init(), ATLAS_2011_S9126244::init(), ATLAS_2011_S9120807::init(), ATLAS_2011_S9108483::init(), ATLAS_2011_S9041966::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S9002537::init(), ATLAS_2011_S8994773::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_S8971293::init(), ATLAS_2011_S8924791::init(), ATLAS_2011_I925932::init(), ATLAS_2011_I919017::init(), ATLAS_2011_CONF_2011_098::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2010_S8919674::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8914702::init(), ATLAS_2010_S8894728::init(), ATLAS_2010_S8817804::init(), ATLAS_2010_S8591806::init(), ATLAS_2010_CONF_2010_049::init(), ALICE_2011_S8945144::init(), ALICE_2011_S8909580::init(), ALICE_2010_S8706239::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), ALEPH_1991_S2435284::init(), IsolationProjection< PROJ1, PROJ2, EST >::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LeptonClusters::LeptonClusters(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), Spherocity::Spherocity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), and VisibleFinalState::VisibleFinalState().

                                                                       {
      const Projection& reg = _addProjection(proj, name);
      const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
      return rtn;
    }

IAnalysisFactory & analysisFactory

(

)

[protected, inherited]

Access the AIDA analysis factory of the controlling AnalysisHandler object.

Definition at line 30 of file Analysis.cc.

References AnalysisHandler::analysisFactory(), and Analysis::handler().

                                              {
    return handler().analysisFactory();
  }

void analyze

(

const Event &

event

)

[inline, virtual]

Analyze one event. A concrete class should here apply the necessary projections on the event and fill the relevant histograms. An overridden function must make sure it first calls the base class function.

Todo:

Resum photons around electrons

Todo:

Resum photons around muons

Implements Analysis.

Definition at line 110 of file MC_SUSY.cc.

References MC_SUSY::_hist_eta_e, MC_SUSY::_hist_eta_gamma, MC_SUSY::_hist_eta_gammaiso, MC_SUSY::_hist_eta_jet, MC_SUSY::_hist_eta_mu, MC_SUSY::_hist_eta_trk, MC_SUSY::_hist_met, MC_SUSY::_hist_mll_2_osof_emu, MC_SUSY::_hist_mll_2_ossf_ee, MC_SUSY::_hist_mll_2_ossf_mumu, MC_SUSY::_hist_mll_all_osof_emu, MC_SUSY::_hist_mll_all_ossf_ee, MC_SUSY::_hist_mll_all_ossf_mumu, MC_SUSY::_hist_mll_osof_emu, MC_SUSY::_hist_mll_ossf_ee, MC_SUSY::_hist_mll_ossf_mumu, MC_SUSY::_hist_n_e, MC_SUSY::_hist_n_gamma, MC_SUSY::_hist_n_gammaiso, MC_SUSY::_hist_n_jet, MC_SUSY::_hist_n_mu, MC_SUSY::_hist_n_trk, MC_SUSY::_hist_phi_e, MC_SUSY::_hist_phi_gamma, MC_SUSY::_hist_phi_gammaiso, MC_SUSY::_hist_phi_jet, MC_SUSY::_hist_phi_mu, MC_SUSY::_hist_phi_trk, MC_SUSY::_hist_pt_e, MC_SUSY::_hist_pt_gamma, MC_SUSY::_hist_pt_gammaiso, MC_SUSY::_hist_pt_jet, MC_SUSY::_hist_pt_mu, MC_SUSY::_hist_pt_trk, Rivet::ANTIMUON, Rivet::deltaR(), Rivet::ELECTRON, FourVector::eta(), Rivet::GeV, JetAlg::jetsByPt(), Rivet::mapAngleMPiToPi(), Vector< N >::mod(), Jet::momentum(), Particle::momentum(), MSG_DEBUG, Rivet::MUON, FinalState::particles(), Particle::pdgId(), FourVector::phi(), Rivet::POSITRON, FourMomentum::pT(), FinalState::size(), Rivet::PID::threeCharge(), MissingMomentum::vectorEt(), and vetoEvent.

                                   {
      const FinalState& tracks = applyProjection<FinalState>(evt, "Tracks");
      if (tracks.particles().empty()) {
        MSG_DEBUG("Failed multiplicity cut");
        vetoEvent;
      }

      // Get event weight
      const double weight = evt.weight();

      // Fill track histos
      _hist_n_trk->fill(tracks.size(), weight);
      foreach (const Particle& t, tracks.particles()) {
        const FourMomentum& p = t.momentum();
        _hist_phi_trk->fill(mapAngleMPiToPi(p.phi()), weight);
        _hist_eta_trk->fill(p.eta(), weight);
        _hist_pt_trk->fill(p.pT()/GeV, weight);
      }

      // Get jets and fill jet histos
      const FastJets& jetpro = applyProjection<FastJets>(evt, "Jets");
      const Jets jets = jetpro.jetsByPt();
      MSG_DEBUG("Jet multiplicity = " << jets.size());
      _hist_n_jet->fill(jets.size(), weight);
      foreach (const Jet& j, jets) {
        const FourMomentum& pj = j.momentum();
        _hist_phi_jet->fill(mapAngleMPiToPi(pj.phi()), weight);
        _hist_eta_jet->fill(pj.eta(), weight);
        _hist_pt_jet->fill(pj.pT()/GeV, weight);
      }

      /// @todo Resum photons around electrons

      // Fill final state electron/positron histos
      const FinalState& efs = applyProjection<FinalState>(evt, "Electrons");
      _hist_n_e->fill(efs.size(), weight);
      vector<FourMomentum> epluses, eminuses;
      foreach (const Particle& e, efs.particles()) {
        const FourMomentum& p = e.momentum();
        _hist_phi_e->fill(mapAngleMPiToPi(p.phi()), weight);
        _hist_eta_e->fill(p.eta(), weight);
        _hist_pt_e->fill(p.pT()/GeV, weight);
        // Add sufficiently hard leptons to collections for m_ll histo
        if (p.pT()/GeV > 20) {
          if (PID::threeCharge(e.pdgId()) > 0) epluses += p; else eminuses += p;
        }
      }

      /// @todo Resum photons around muons

      // Fill final state muon/antimuon histos
      const FinalState& mufs = applyProjection<FinalState>(evt, "Muons");
      _hist_n_mu->fill(mufs.size(), weight);
      vector<FourMomentum> mupluses, muminuses;
      foreach (const Particle& mu, mufs.particles()) {
        const FourMomentum& p = mu.momentum();
        _hist_phi_mu->fill(mapAngleMPiToPi(p.phi()), weight);
        _hist_eta_mu->fill(p.eta(), weight);
        _hist_pt_mu->fill(p.pT()/GeV, weight);
        // Add sufficiently hard leptons to collections for m_ll histo
        if (p.pT()/GeV > 20) {
          if (PID::threeCharge(mu.pdgId()) > 0) mupluses += p; else muminuses += p;
        }
      }

      // Fill final state non-isolated photon histos
      const FinalState& allphotonfs = applyProjection<FinalState>(evt, "AllPhotons");
      _hist_n_gamma->fill(allphotonfs.size(), weight);
      ParticleVector isolatedphotons;
      foreach (const Particle& ph, allphotonfs.particles()) {
        const FourMomentum& p = ph.momentum();
        _hist_phi_gamma->fill(mapAngleMPiToPi(p.phi()), weight);
        _hist_eta_gamma->fill(p.eta(), weight);
        _hist_pt_gamma->fill(p.pT()/GeV, weight);
        // Select isolated photons
        bool isolated = true;
        foreach (const Jet& j, jets) {
          if (deltaR(j.momentum(), p) < 0.2) {
            isolated = false;
            break;
          }
        }
        if (isolated) isolatedphotons += ph;
      }


      // Fill final state isolated photon histos
      _hist_n_gammaiso->fill(isolatedphotons.size(), weight);
      foreach (const Particle& ph_iso, isolatedphotons) {
        const FourMomentum& p = ph_iso.momentum();
        _hist_phi_gammaiso->fill(mapAngleMPiToPi(p.phi()), weight);
        _hist_eta_gammaiso->fill(p.eta(), weight);
        _hist_pt_gammaiso->fill(p.pT()/GeV, weight);
      }

      // Calculate and fill missing Et histos
      const MissingMomentum& met = applyProjection<MissingMomentum>(evt, "MET");
      _hist_met->fill(met.vectorEt().mod()/GeV);

      // Choose highest-pT leptons of each sign and flavour for dilepton mass edges
      const FinalState& lpfs = applyProjection<FinalState>(evt, "LeadingParticles");
      bool eplus_ok(false), eminus_ok(false), muplus_ok(false), muminus_ok(false);
      FourMomentum peplus, peminus, pmuplus, pmuminus;
      foreach (const Particle& p, lpfs.particles()) {
        // Only use leptons above 20 GeV
        if (p.momentum().pT()/GeV < 20) continue;
        // Identify the PID
        const PdgId pid = p.pdgId();
        if (pid == ELECTRON) {
          eminus_ok = true;
          peminus = p.momentum();
        } else if (pid == POSITRON) {
          eplus_ok = true;
          peplus = p.momentum();
        } else if (pid == MUON) {
          muminus_ok = true;
          pmuminus = p.momentum();
        } else if (pid == ANTIMUON) {
          muplus_ok = true;
          pmuplus = p.momentum();
        } else {
          throw Error("Unexpected particle type in leading particles FS!");
        }
      }
      // m_ee
      if (eminus_ok && eplus_ok) {
        const double m_ee = FourMomentum(peplus + peminus).mass();
        _hist_mll_ossf_ee->fill(m_ee/GeV, weight);
        if (epluses.size() == 1 && eminuses.size() == 1)
          _hist_mll_2_ossf_ee->fill(m_ee/GeV, weight);
      }
      // m_mumu
      if (muminus_ok && muplus_ok) {
        const double m_mumu = FourMomentum(pmuplus + pmuminus).mass();
        _hist_mll_ossf_mumu->fill(m_mumu/GeV, weight);
        if (mupluses.size() == 1 && muminuses.size() == 1)
          _hist_mll_2_ossf_mumu->fill(m_mumu/GeV, weight);
      }
      // m_emu (both configurations)
      if (eminus_ok && muplus_ok) {
        const double m_emu = FourMomentum(pmuplus + peminus).mass();
        _hist_mll_osof_emu->fill(m_emu/GeV, weight);
        if (mupluses.size() == 1 && eminuses.size() == 1)
          _hist_mll_2_osof_emu->fill(m_emu/GeV, weight);

      }
      if (muminus_ok && eplus_ok) {
        const double m_mue = FourMomentum(peplus + pmuminus).mass();
        _hist_mll_osof_emu->fill(m_mue/GeV, weight);
        if (epluses.size() == 1 && muminuses.size() == 1)
          _hist_mll_2_osof_emu->fill(m_mue/GeV, weight);
      }


      // m_ll plots using *all* electrons, positrons, muons and antimuons
      // m_ee
      foreach (const FourMomentum& peplus, epluses) {
        foreach (const FourMomentum& peminus, eminuses) {
          const double m_ee = FourMomentum(peplus + peminus).mass();
          _hist_mll_all_ossf_ee->fill(m_ee/GeV, weight);
        }
      }
      // m_mumu
      foreach (const FourMomentum& pmuplus, mupluses) {
        foreach (const FourMomentum& pmuminus, muminuses) {
          const double m_mumu = FourMomentum(pmuplus + pmuminus).mass();
          _hist_mll_all_ossf_mumu->fill(m_mumu/GeV, weight);
        }
      }
      // m_emu (both configurations)
      foreach (const FourMomentum& pmuplus, mupluses) {
        foreach (const FourMomentum& peminus, eminuses) {
          const double m_emu = FourMomentum(pmuplus + peminus).mass();
          _hist_mll_all_osof_emu->fill(m_emu/GeV, weight);
        }
      }
      foreach (const FourMomentum& peplus, epluses) {
        foreach (const FourMomentum& pmuminus, muminuses) {
          const double m_mue = FourMomentum(peplus + pmuminus).mass();
          _hist_mll_all_osof_emu->fill(m_mue/GeV, weight);
        }
      }

    }

const PROJ& applyProjection	(	const Event &	evt,
		const std::string &	name
	)			const [inline, inherited]

Apply the named projection on event.

Definition at line 81 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                               {
      return pcast<PROJ>(_applyProjection(evt, name));
    }

const PROJ& applyProjection	(	const Event &	evt,
		const Projection &	proj
	)			const [inline, inherited]

Apply the supplied projection on event.

Definition at line 74 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                                {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }

const PROJ& applyProjection	(	const Event &	evt,
		const PROJ &	proj
	)			const [inline, inherited]

Apply the supplied projection on event.

Definition at line 67 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

Referenced by DISFinalState::project().

                                                                          {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }

virtual std::vector<std::string> authors

(

)

const [inline, virtual, inherited]

Names & emails of paper/analysis authors.

Names and email of authors in 'NAME <EMAIL>' format. The first name in the list should be the primary contact person.

Definition at line 136 of file Analysis.hh.

References AnalysisInfo::authors(), and Analysis::info().

                                                 {
      return info().authors();
    }

const PdgIdPair beamIds

(

)

const [inherited]

Incoming beam IDs for this run.

Definition at line 58 of file Analysis.cc.

References AnalysisHandler::beamIds(), and Analysis::handler().

Referenced by UA5_1982_S875503::finalize(), and UA5_1982_S875503::init().

                                          {
    return handler().beamIds();
  }

const ParticlePair & beams

(

)

const [inherited]

Incoming beams for this run.

Definition at line 54 of file Analysis.cc.

References AnalysisHandler::beams(), and Analysis::handler().

Referenced by OPAL_1998_S3780481::analyze(), MC_WPOL::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1996_S3430090::analyze(), DELPHI_1995_S3137023::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_2004_S5765862::analyze(), ALEPH_1996_S3486095::analyze(), MC_WPOL::init(), and Analysis::isCompatible().

                                            {
    return handler().beams();
  }

virtual std::string bibKey

(

)

const [inline, virtual, inherited]

BibTeX citation key for this article.

Definition at line 189 of file Analysis.hh.

References AnalysisInfo::bibKey(), and Analysis::info().

                                     {
      return info().bibKey();
    }

virtual std::string bibTeX

(

)

const [inline, virtual, inherited]

BibTeX citation entry for this article.

Definition at line 194 of file Analysis.hh.

References AnalysisInfo::bibTeX(), and Analysis::info().

                                     {
      return info().bibTeX();
    }

const BinEdges & binEdges	(	size_t	datasetId,
		size_t	xAxisId,
		size_t	yAxisId
	)			const [protected, inherited]

Get bin edges for a numbered histo (using ref AIDA caching).

Definition at line 218 of file Analysis.cc.

References Analysis::binEdges(), and Analysis::makeAxisCode().

                                                                                           {
    const string hname = makeAxisCode(datasetId, xAxisId, yAxisId);
    return binEdges(hname);
  }

const BinEdges & binEdges

(

const std::string &

hname

)

const [protected, inherited]

Get bin edges for a named histo (using ref AIDA caching).

Definition at line 203 of file Analysis.cc.

References Analysis::_cacheBinEdges(), Analysis::_histBinEdges, Analysis::getLog(), MSG_TRACE, Analysis::name(), and Log::TRACE.

Referenced by ATLAS_2011_S8994773::analyze(), ATLAS_2010_S8894728::analyze(), Analysis::binEdges(), Analysis::bookHistogram1D(), Analysis::bookProfile1D(), ATLAS_2011_S9126244::finalize(), D0_2008_S7837160::init(), CMS_2011_S9088458::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9128077::init(), ATLAS_2011_S9002537::init(), ATLAS_2010_S8894728::init(), ALICE_2011_S8909580::init(), and ATLAS_2011_S9126244::initializePlots().

                                                              {
    _cacheBinEdges();
    MSG_TRACE("Using histo bin edges for " << name() << ":" << hname);
    const BinEdges& edges = _histBinEdges.find(hname)->second;
    if (getLog().isActive(Log::TRACE)) {
      stringstream edges_ss;
      foreach (const double be, edges) {
        edges_ss << " " << be;
      }
      MSG_TRACE("Edges:" << edges_ss.str());
    }
    return edges;
  }

IDataPointSet * bookDataPointSet	(	size_t	datasetId,
		size_t	xAxisId,
		size_t	yAxisId,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 2-dimensional data point set based on the corresponding AIDA data file. The binnings (x-errors) will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property. Book a 2-dimensional data point set based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings (x-errors) will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.

Definition at line 440 of file Analysis.cc.

References Analysis::_cacheXAxisData(), Analysis::_dpsData, Analysis::bookDataPointSet(), Analysis::makeAxisCode(), MSG_TRACE, and Analysis::name().

                                                                                        {
    // Get the bin edges (only read the AIDA file once)
    _cacheXAxisData();
    // Build the axis code
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    //const map<string, vector<DPSXPoint> > xpoints = getDPSXValsErrs(papername);
    MSG_TRACE("Using DPS x-positions for " << name() << ":" << axisCode);
    IDataPointSet* dps = bookDataPointSet(axisCode, title, xtitle, ytitle);
    const vector<DPSXPoint> xpts = _dpsData.find(axisCode)->second;
    for (size_t pt = 0; pt < xpts.size(); ++pt) {
      dps->addPoint();
      IMeasurement* meas = dps->point(pt)->coordinate(0);
      meas->setValue(xpts[pt].val);
      meas->setErrorPlus(xpts[pt].errplus);
      meas->setErrorMinus(xpts[pt].errminus);
    }
    MSG_TRACE("Made DPS " << axisCode <<  " for " << name());
    return dps;
  }

IDataPointSet * bookDataPointSet	(	const std::string &	name,
		size_t	npts,
		double	lower,
		double	upper,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 2-dimensional data point set with equally spaced points in a range. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 422 of file Analysis.cc.

References Analysis::bookDataPointSet().

                                                                                        {
    IDataPointSet* dps = bookDataPointSet(hname, title, xtitle, ytitle);
    for (size_t pt = 0; pt < npts; ++pt) {
      const double binwidth = (upper-lower)/npts;
      const double bincentre = lower + (pt + 0.5) * binwidth;
      dps->addPoint();
      IMeasurement* meas = dps->point(pt)->coordinate(0);
      meas->setValue(bincentre);
      meas->setErrorPlus(binwidth/2.0);
      meas->setErrorMinus(binwidth/2.0);
    }
    return dps;
  }

IDataPointSet * bookDataPointSet	(	const std::string &	name,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 2-dimensional data point set. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 407 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::datapointsetFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by Analysis::bookDataPointSet(), ALEPH_2004_S5765862::finalize(), UA5_1988_S1867512::init(), STAR_2006_S6860818::init(), OPAL_1993_S2692198::init(), MC_XS::init(), MC_JetAnalysis::init(), JADE_OPAL_2000_S4300807::init(), D0_2001_S4674421::init(), CMS_2011_S9088458::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2005_S6217184::init(), CDF_1996_S3418421::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9128077::init(), ATLAS_2011_S9002537::init(), ATLAS_2010_S8894728::init(), ALEPH_2004_S5765862::init(), and ATLAS_2011_S9126244::initializePlots().

                                                                                        {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IDataPointSet* dps = datapointsetFactory().create(path, title, 2);
    MSG_TRACE("Made data point set " << hname <<  " for " << name());
    dps->setXTitle(xtitle);
    dps->setYTitle(ytitle);
    return dps;
  }

IHistogram1D * bookHistogram1D	(	size_t	datasetId,
		size_t	xAxisId,
		size_t	yAxisId,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D histogram based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.

Definition at line 237 of file Analysis.cc.

References Analysis::bookHistogram1D(), and Analysis::makeAxisCode().

  {
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    return bookHistogram1D(axisCode, title, xtitle, ytitle);
  }

IHistogram1D * bookHistogram1D	(	const std::string &	name,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D histogram based on the name in the corresponding AIDA file. The binnings will be obtained by reading the bundled AIDA data record file with the same filename as the analysis' name() property.

Definition at line 246 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

  {
    // Get the bin edges (only read the AIDA file once)
    const BinEdges edges = binEdges(hname);
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, edges);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setXTitle(xtitle);
    hist->setYTitle(ytitle);
    return hist;
  }

IHistogram1D * bookHistogram1D	(	const std::string &	name,
		const std::vector< double > &	binedges,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D histogram with non-uniform bins defined by the vector of bin edges binedges . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 281 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                                                      {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, binedges);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setXTitle(xtitle);
    hist->setYTitle(ytitle);
    return hist;
  }

IHistogram1D * bookHistogram1D	(	const std::string &	name,
		size_t	nbins,
		double	lower,
		double	upper,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D histogram with nbins uniformly distributed across the range lower - upper . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 264 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by Analysis::bookHistogram1D(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_XS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_VH2BB::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), LHCB_2010_S8758301::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2003_WUD_03_11::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), DELPHI_1995_S3137023::init(), D0_2010_S8821313::init(), D0_2010_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CMS_2011_S9120041::init(), CMS_2011_S9088458::init(), CMS_2011_S9086218::init(), CMS_2011_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::init(), CMS_2011_S8950903::init(), CMS_2011_S8884919::init(), CMS_2010_S8656010::init(), CMS_2010_S8547297::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2011_S9225137::init(), ATLAS_2011_S9212183::init(), ATLAS_2011_S9131140::init(), ATLAS_2011_S9128077::init(), ATLAS_2011_S9120807::init(), ATLAS_2011_S9108483::init(), ATLAS_2011_S9041966::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S9002537::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_S8971293::init(), ATLAS_2011_I925932::init(), ATLAS_2011_I919017::init(), ATLAS_2011_CONF_2011_098::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2010_S8919674::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8914702::init(), ATLAS_2010_S8817804::init(), ATLAS_2010_S8591806::init(), ATLAS_2010_CONF_2010_049::init(), ALICE_2011_S8945144::init(), ALICE_2011_S8909580::init(), ALICE_2010_S8706239::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), ALEPH_1991_S2435284::init(), and ATLAS_2011_S9126244::initializePlots().

                                                                                      {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, nbins, lower, upper);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setXTitle(xtitle);
    hist->setYTitle(ytitle);
    return hist;
  }

IHistogram2D * bookHistogram2D	(	const std::string &	name,
		const std::vector< double > &	xbinedges,
		const std::vector< double > &	ybinedges,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = "",
		const std::string &	ztitle = ""
	)			[protected, inherited]

Book a 2D histogram with non-uniform bins defined by the vectorx of bin edges xbinedges and ybinedges. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 320 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                            {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IHistogram2D* hist =
      histogramFactory().createHistogram2D(path, title, xbinedges, ybinedges);
    MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
    hist->setXTitle(xtitle);
    hist->setYTitle(ytitle);
    hist->setZTitle(ztitle);
    return hist;
  }

IHistogram2D * bookHistogram2D	(	const std::string &	name,
		size_t	nxbins,
		double	xlower,
		double	xupper,
		size_t	nybins,
		double	ylower,
		double	yupper,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = "",
		const std::string &	ztitle = ""
	)			[protected, inherited]

Book a 2D histogram with nxbins and nybins uniformly distributed across the ranges xlower - xupper and ylower - yupper respectively along the x- and y-axis. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 298 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                            {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IHistogram2D* hist =
      histogramFactory().createHistogram2D(path, title, nxbins, xlower, xupper,
                       nybins, ylower, yupper);
    MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
    hist->setXTitle(xtitle);
    hist->setYTitle(ytitle);
    hist->setZTitle(ztitle);
    return hist;
  }

IProfile1D * bookProfile1D	(	size_t	datasetId,
		size_t	xAxisId,
		size_t	yAxisId,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D profile histogram based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.

Definition at line 343 of file Analysis.cc.

References Analysis::bookProfile1D(), and Analysis::makeAxisCode().

                                                                                  {
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    return bookProfile1D(axisCode, title, xtitle, ytitle);
  }

IProfile1D * bookProfile1D	(	const std::string &	name,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D profile histogram based on the name in the corresponding AIDA file. The binnings will be obtained by reading the bundled AIDA data record file with the same filename as the analysis' name() property.

Definition at line 351 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

  {
    // Get the bin edges (only read the AIDA file once)
    const BinEdges edges = binEdges(hname);
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IProfile1D* prof = histogramFactory().createProfile1D(path, title, edges);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setXTitle(xtitle);
    prof->setYTitle(ytitle);
    return prof;
  }

IProfile1D * bookProfile1D	(	const std::string &	name,
		const std::vector< double > &	binedges,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D profile histogram with non-uniform bins defined by the vector of bin edges binedges . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 386 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                                                  {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IProfile1D* prof = histogramFactory().createProfile1D(path, title, binedges);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setXTitle(xtitle);
    prof->setYTitle(ytitle);
    return prof;
  }

IProfile1D * bookProfile1D	(	const std::string &	name,
		size_t	nbins,
		double	lower,
		double	upper,
		const std::string &	title = "",
		const std::string &	xtitle = "",
		const std::string &	ytitle = ""
	)			[protected, inherited]

Book a 1D profile histogram with nbins uniformly distributed across the range lower - upper . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)

Definition at line 369 of file Analysis.cc.

References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by Analysis::bookProfile1D(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2006_S6860818::init(), MC_WPOL::init(), MC_PHOTONJETUE::init(), MC_LEADINGJETS::init(), MC_GENERIC::init(), H1_2000_S4129130::init(), H1_1994_S2919893::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), D0_1996_S3324664::init(), CMS_2011_S9120041::init(), CMS_2011_S8884919::init(), CDF_2010_S8591881_QCD::init(), CDF_2010_S8591881_DY::init(), CDF_2009_S8233977::init(), CDF_2008_S7782535::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2005_S6217184::init(), CDF_2004_S5839831::init(), CDF_2001_S4751469::init(), ATLAS_2011_S8994773::init(), ATLAS_2011_S8924791::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8894728::init(), ATLAS_2010_S8591806::init(), ALICE_2011_S8945144::init(), ALICE_2010_S8706239::init(), and ATLAS_2011_S9126244::initializePlots().

                                                                                  {
    _makeHistoDir();
    const string path = histoPath(hname);
    if (path.find(" ") != string::npos) {
      throw Error("Histogram path '" + path + "' is invalid: spaces are not permitted in paths");
    }
    IProfile1D* prof = histogramFactory().createProfile1D(path, title, nbins, lower, upper);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setXTitle(xtitle);
    prof->setYTitle(ytitle);
    return prof;
  }

virtual std::string collider

(

)

const [inline, virtual, inherited]

Collider on which the experiment ran.

Definition at line 174 of file Analysis.hh.

References AnalysisInfo::collider(), and Analysis::info().

                                       {
      return info().collider();
    }

double crossSection

(

)

const [protected, inherited]

Get the process cross-section in pb. Throws if this hasn't been set.

Definition at line 166 of file Analysis.cc.

References Analysis::_crossSection, Analysis::_gotCrossSection, and Analysis::name().

Referenced by STAR_2006_S6870392::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_XS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_VH2BB::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_DIPHOTON::finalize(), LHCB_2010_S8758301::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S9086218::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7541902::finalize(), CDF_2008_S7540469::finalize(), CDF_2006_S6653332::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2001_S4517016::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_S9108483::finalize(), ATLAS_2011_S9041966::finalize(), ATLAS_2011_I919017::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), and ATLAS_2010_CONF_2010_049::finalize().

                                      {
    if (!_gotCrossSection || std::isnan(_crossSection)) {
      string errMsg = "You did not set the cross section for the analysis " + name();
      throw Error(errMsg);
    }
    return _crossSection;
  }

double crossSectionPerEvent

(

)

const [protected, inherited]

Get the process cross-section per generated event in pb. Throws if this hasn't been set.

Definition at line 174 of file Analysis.cc.

References Analysis::_crossSection, and Analysis::sumOfWeights().

Referenced by UA1_1990_S2044935::finalize(), MC_WPOL::finalize(), MC_PHOTONJETS::finalize(), D0_2010_S8671338::finalize(), D0_2009_S8349509::finalize(), D0_2001_S4674421::finalize(), CDF_2007_S7057202::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2011_S9128077::finalize(), and ATLAS_2010_S8817804::finalize().

                                              {
    const double sumW = sumOfWeights();
    assert(sumW != 0.0);
    return _crossSection / sumW;
  }

IDataPointSetFactory & datapointsetFactory

(

)

[protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 45 of file Analysis.cc.

References AnalysisHandler::datapointsetFactory(), and Analysis::handler().

Referenced by Analysis::bookDataPointSet(), and Analysis::scale().

                                                      {
    return handler().datapointsetFactory();
  }

virtual std::string description

(

)

const [inline, virtual, inherited]

Get a full description of the analysis.

Full textual description of this analysis, what it is useful for, what experimental techniques are applied, etc. Should be treated as a chunk of restructuredText (http://docutils.sourceforge.net/rst.html), with equations to be rendered as LaTeX with amsmath operators.

Definition at line 155 of file Analysis.hh.

References AnalysisInfo::description(), and Analysis::info().

                                          {
      return info().description();
    }

virtual std::string experiment

(

)

const [inline, virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 169 of file Analysis.hh.

References AnalysisInfo::experiment(), and Analysis::info().

                                         {
      return info().experiment();
    }

void finalize

(

)

[inline, virtual]

Finalize this analysis object. A concrete class should here make all necessary operations on the histograms. Writing the histograms to a file is, however, done by the Rivet class. An overridden function must make sure it first calls the base class function.

Todo:

Normalisations

Implements Analysis.

Definition at line 296 of file MC_SUSY.cc.

                    {
      /// @todo Normalisations
    }

Log & getLog

(

)

const [protected, inherited]

Get a Log object based on the name() property of the calling analysis object.

Reimplemented from ProjectionApplier.

Definition at line 97 of file Analysis.cc.

References Analysis::name().

Referenced by Analysis::binEdges().

                              {
    string logname = "Rivet.Analysis." + name();
    return Log::getLog(logname);
  }

const Projection& getProjection

(

const std::string &

name

)

const [inline, inherited]

Get the named projection (non-templated, so returns as a reference to a Projection base class).

Definition at line 57 of file ProjectionApplier.hh.

References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler().

                                                                 {
      return getProjHandler().getProjection(*this, name);
    }

const PROJ& getProjection

(

const std::string &

name

)

const [inline, inherited]

Get the named projection, specifying return type via a template argument.

Definition at line 49 of file ProjectionApplier.hh.

References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler().

Referenced by ProjectionApplier::_applyProjection(), Rivet::pcmp(), and Hemispheres::project().

                                                           {
      const Projection& p = getProjHandler().getProjection(*this, name);
      return pcast<PROJ>(p);
    }

std::set<ConstProjectionPtr> getProjections

(

)

const [inline, inherited]

Get the contained projections, including recursion.

Definition at line 42 of file ProjectionApplier.hh.

References ProjectionHandler::DEEP, ProjectionHandler::getChildProjections(), and ProjectionApplier::getProjHandler().

Referenced by Projection::beamPairs().

                                                      {
      return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
    }

ProjectionHandler& getProjHandler

(

)

const [inline, protected, inherited]

Get a reference to the ProjectionHandler for this thread.

Definition at line 94 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

Referenced by ProjectionApplier::_addProjection(), ProjectionApplier::getProjection(), ProjectionApplier::getProjections(), and ProjectionApplier::~ProjectionApplier().

                                              {
      return _projhandler;
    }

AnalysisHandler& handler

(

)

const [inline, inherited]

Access the controlling AnalysisHandler object.

Definition at line 296 of file Analysis.hh.

References Analysis::_analysishandler.

Referenced by Analysis::analysisFactory(), Analysis::beamIds(), Analysis::beams(), Analysis::datapointsetFactory(), Analysis::histoDir(), Analysis::histogramFactory(), Analysis::numEvents(), Analysis::sqrtS(), Analysis::sumOfWeights(), and Analysis::tree().

{ return *_analysishandler; }

const string histoDir

(

)

const [protected, inherited]

Get the canonical histogram "directory" path for this analysis.

Todo:

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

Definition at line 63 of file Analysis.cc.

References Analysis::handler(), Analysis::name(), and AnalysisHandler::runName().

Referenced by Analysis::_makeHistoDir(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), MC_WJETS::finalize(), D0_2008_S7837160::finalize(), D0_2008_S7719523::finalize(), CMS_2011_S8978280::finalize(), ATLAS_2011_S9128077::finalize(), and Analysis::histoPath().

                                        {
    /// @todo This doesn't change: calc and cache at first use!
    string path = "/" + name();
    if (handler().runName().length() > 0) {
      path = "/" + handler().runName() + path;
    }
    while (find_first(path, "//")) {
      replace_all(path, "//", "/");
    }
    return path;
  }

IHistogramFactory & histogramFactory

(

)

[protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 40 of file Analysis.cc.

References Analysis::handler(), and AnalysisHandler::histogramFactory().

Referenced by Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_GENERIC::finalize(), H1_1995_S3167097::finalize(), D0_2008_S7837160::finalize(), D0_2008_S7719523::finalize(), CMS_2011_S9088458::finalize(), CMS_2011_S8978280::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2011_S9126244::finalize(), ATLAS_2011_S9002537::finalize(), ALICE_2011_S8909580::finalize(), and ATLAS_2011_S9126244::finalizeQ0GapFraction().

                                                {
    return handler().histogramFactory();
  }

const string histoPath

(

const std::string &

hname

)

const [protected, inherited]

Get the canonical histogram path for the named histogram in this analysis.

Definition at line 76 of file Analysis.cc.

References Analysis::histoDir().

Referenced by Analysis::bookDataPointSet(), Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), MC_JetAnalysis::finalize(), MC_GENERIC::finalize(), ATLAS_2011_S9126244::finalize(), and ALICE_2011_S8909580::finalize().

                                                            {
    const string path = histoDir() + "/" + hname;
    return path;
  }

AnalysisInfo& info

(

)

[inline, inherited]

Get the actual AnalysisInfo object in which all this metadata is stored (non-const).

Definition at line 248 of file Analysis.hh.

References Analysis::_info.

                         {
      assert(_info.get() != 0 && "No AnalysisInfo object :O");
      return *_info;
    }

const AnalysisInfo& info

(

)

const [inline, inherited]

Get the actual AnalysisInfo object in which all this metadata is stored.

Definition at line 108 of file Analysis.hh.

References Analysis::_info.

Referenced by Analysis::authors(), Analysis::bibKey(), Analysis::bibTeX(), Analysis::collider(), Analysis::description(), Analysis::experiment(), Analysis::inspireId(), Analysis::name(), Analysis::needsCrossSection(), Analysis::references(), Analysis::requiredBeams(), Analysis::requiredEnergies(), Analysis::runInfo(), Analysis::setNeedsCrossSection(), Analysis::setRequiredBeams(), Analysis::setRequiredEnergies(), Analysis::spiresId(), Analysis::status(), Analysis::summary(), Analysis::todos(), and Analysis::year().

                                     {
      assert(_info.get() != 0 && "No AnalysisInfo object :O");
      return *_info;
    }

void init

(

)

[inline, virtual]

Initialize this analysis object. A concrete class should here book all necessary histograms. An overridden function must make sure it first calls the base class function.

Todo:

LSP eta, pT, phi, mass: no reliable cross-scenario LSP PID but maybe plot for all of chi^0_1, gravitino, sneutrino, gluino, ... or identify the LSP as any PID::isSUSY (?) particle with status = 1?

Implements Analysis.

Definition at line 31 of file MC_SUSY.cc.

References MC_SUSY::_hist_eta_e, MC_SUSY::_hist_eta_gamma, MC_SUSY::_hist_eta_gammaiso, MC_SUSY::_hist_eta_jet, MC_SUSY::_hist_eta_mu, MC_SUSY::_hist_eta_trk, MC_SUSY::_hist_met, MC_SUSY::_hist_mll_2_osof_emu, MC_SUSY::_hist_mll_2_ossf_ee, MC_SUSY::_hist_mll_2_ossf_mumu, MC_SUSY::_hist_mll_all_osof_emu, MC_SUSY::_hist_mll_all_ossf_ee, MC_SUSY::_hist_mll_all_ossf_mumu, MC_SUSY::_hist_mll_osof_emu, MC_SUSY::_hist_mll_ossf_ee, MC_SUSY::_hist_mll_ossf_mumu, MC_SUSY::_hist_n_e, MC_SUSY::_hist_n_gamma, MC_SUSY::_hist_n_gammaiso, MC_SUSY::_hist_n_jet, MC_SUSY::_hist_n_mu, MC_SUSY::_hist_n_trk, MC_SUSY::_hist_phi_e, MC_SUSY::_hist_phi_gamma, MC_SUSY::_hist_phi_gammaiso, MC_SUSY::_hist_phi_jet, MC_SUSY::_hist_phi_mu, MC_SUSY::_hist_phi_trk, MC_SUSY::_hist_pt_e, MC_SUSY::_hist_pt_gamma, MC_SUSY::_hist_pt_gammaiso, MC_SUSY::_hist_pt_jet, MC_SUSY::_hist_pt_mu, MC_SUSY::_hist_pt_trk, IdentifiedFinalState::acceptId(), IdentifiedFinalState::acceptIdPair(), LeadingParticlesFinalState::addParticleIdPair(), ProjectionApplier::addProjection(), FastJets::ANTIKT, Analysis::bookHistogram1D(), Rivet::ELECTRON, Rivet::GeV, Rivet::MUON, Rivet::PHOTON, and Rivet::PI.

                {
      // Basic final state
      const FinalState fs(-4.0, 4.0, 10*GeV);

      // Tracks and jets
      addProjection(ChargedFinalState(fs), "Tracks");
      addProjection(FastJets(fs, FastJets::ANTIKT, 0.7), "Jets");

      IdentifiedFinalState photonfs(fs);
      photonfs.acceptId(PHOTON);
      addProjection(photonfs, "AllPhotons");

      IdentifiedFinalState efs(fs);
      efs.acceptIdPair(ELECTRON);
      addProjection(efs, "Electrons");

      IdentifiedFinalState mufs(fs);
      mufs.acceptIdPair(MUON);
      addProjection(mufs, "Muons");

      MissingMomentum missing(fs);
      addProjection(missing, "MET");

      LeadingParticlesFinalState lpfs(fs);
      lpfs.addParticleIdPair(ELECTRON);
      lpfs.addParticleIdPair(MUON);
      addProjection(lpfs, "LeadingParticles");

      _hist_n_trk   = bookHistogram1D("n-trk", 50, 0.5, 300.5);
      _hist_phi_trk = bookHistogram1D("phi-trk", 50, -PI, PI);
      _hist_eta_trk = bookHistogram1D("eta-trk", 50, -4, 4);
      _hist_pt_trk  = bookHistogram1D("pt-trk", 100, 0.0, 1500);

      _hist_n_jet   = bookHistogram1D("n-jet", 21, -0.5, 20.5);
      _hist_phi_jet = bookHistogram1D("phi-jet", 50, -PI, PI);
      _hist_eta_jet = bookHistogram1D("eta-jet", 50, -4, 4);
      _hist_pt_jet  = bookHistogram1D("pt-jet", 100, 0.0, 1500);

      _hist_n_e   = bookHistogram1D("n-e", 11, -0.5, 10.5);
      _hist_phi_e = bookHistogram1D("phi-e", 50, -PI, PI);
      _hist_eta_e = bookHistogram1D("eta-e", 50, -4, 4);
      _hist_pt_e  = bookHistogram1D("pt-e", 100, 0.0, 500);

      _hist_n_mu   = bookHistogram1D("n-mu", 11, -0.5, 10.5);
      _hist_phi_mu = bookHistogram1D("phi-mu", 50, -PI, PI);
      _hist_eta_mu = bookHistogram1D("eta-mu", 50, -4, 4);
      _hist_pt_mu  = bookHistogram1D("pt-mu", 100, 0.0, 500);

      _hist_n_gamma   = bookHistogram1D("n-gamma", 11, -0.5, 10.5);
      _hist_phi_gamma = bookHistogram1D("phi-gamma", 50, -PI, PI);
      _hist_eta_gamma = bookHistogram1D("eta-gamma", 50, -4, 4);
      _hist_pt_gamma  = bookHistogram1D("pt-gamma", 100, 0.0, 500);

      _hist_n_gammaiso   = bookHistogram1D("n-gamma-iso", 11, -0.5, 10.5);
      _hist_phi_gammaiso = bookHistogram1D("phi-gamma-iso", 50, -PI, PI);
      _hist_eta_gammaiso = bookHistogram1D("eta-gamma-iso", 50, -4, 4);
      _hist_pt_gammaiso  = bookHistogram1D("pt-gamma-iso", 100, 0.0, 500);

      _hist_met = bookHistogram1D("Etmiss", 100, 0.0, 1500);

      _hist_mll_ossf_ee   = bookHistogram1D("mll-ossf-ee", 50, 0.0, 500);
      _hist_mll_ossf_mumu = bookHistogram1D("mll-ossf-mumu", 50, 0.0, 500);
      _hist_mll_osof_emu  = bookHistogram1D("mll-osof-emu", 50, 0.0, 500);

      _hist_mll_all_ossf_ee   = bookHistogram1D("mll-all-ossf-ee", 50, 0.0, 500);
      _hist_mll_all_ossf_mumu = bookHistogram1D("mll-all-ossf-mumu", 50, 0.0, 500);
      _hist_mll_all_osof_emu  = bookHistogram1D("mll-all-osof-emu", 50, 0.0, 500);

      _hist_mll_2_ossf_ee   = bookHistogram1D("mll-2-ossf-ee", 50, 0.0, 500);
      _hist_mll_2_ossf_mumu = bookHistogram1D("mll-2-ossf-mumu", 50, 0.0, 500);
      _hist_mll_2_osof_emu  = bookHistogram1D("mll-2-osof-emu", 50, 0.0, 500);

      /// @todo LSP eta, pT, phi, mass: no reliable cross-scenario LSP PID but
      /// maybe plot for all of chi^0_1, gravitino, sneutrino, gluino, ... or
      /// identify the LSP as any PID::isSUSY (?) particle with status = 1?
    }

virtual std::string inspireId

(

)

const [inline, virtual, inherited]

Get the Inspire ID code for this analysis.

Definition at line 123 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::inspireId().

                                        {
      return info().inspireId();
    }

bool isCompatible	(	const PdgIdPair &	beams,
		const std::pair< double, double > &	energies
	)			const [inherited]

Check if analysis is compatible with the provided beam particle IDs and energies.

bool isCompatible	(	PdgId	beam1,
		PdgId	beam2,
		double	e1,
		double	e2
	)			const [inherited]

Check if analysis is compatible with the provided beam particle IDs and energies.

Definition at line 122 of file Analysis.cc.

References Analysis::beams(), and Analysis::isCompatible().

                                                                                  {
    PdgIdPair beams(beam1, beam2);
    pair<double,double> energies(e1, e2);
    return isCompatible(beams, energies);
  }

bool isCompatible

(

const ParticlePair &

beams

)

const [inherited]

Check if analysis is compatible with the provided beam particle IDs and energies.

Definition at line 116 of file Analysis.cc.

Referenced by Analysis::isCompatible().

                                                             {
    return isCompatible(beams.first.pdgId(),  beams.second.pdgId(),
                        beams.first.energy(), beams.second.energy());
  }

BinEdges logBinEdges	(	size_t	nbins,
		double	lower,
		double	upper
	)			[protected, inherited]

Get bin edges with logarithmic widths.

Definition at line 224 of file Analysis.cc.

Referenced by MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), and MC_DIPHOTON::init().

                                                                         {
    assert(lower>0.0);
    assert(upper>lower);
    double loglower=log10(lower);
    double logupper=log10(upper);
    vector<double> binedges;
    double stepwidth=(logupper-loglower)/double(nbins);
    for (size_t i=0; i<=nbins; ++i) {
      binedges.push_back(pow(10.0, loglower+double(i)*stepwidth));
    }
    return binedges;
  }

const string makeAxisCode	(	size_t	datasetId,
		size_t	xAxisId,
		size_t	yAxisId
	)			const [protected, inherited]

Get the internal histogram name for given d, x and y (cf. HepData).

Definition at line 82 of file Analysis.cc.

Referenced by Analysis::binEdges(), Analysis::bookDataPointSet(), Analysis::bookHistogram1D(), Analysis::bookProfile1D(), and ATLAS_2011_S9126244::finalize().

                                                                                            {
    stringstream axisCode;
    axisCode << "d";
    if (datasetId < 10) axisCode << 0;
    axisCode << datasetId;
    axisCode << "-x";
    if (xAxisId < 10) axisCode << 0;
    axisCode << xAxisId;
    axisCode << "-y";
    if (yAxisId < 10) axisCode << 0;
    axisCode << yAxisId;
    return axisCode.str();
  }

virtual std::string name

(

)

const [inline, virtual, inherited]

Get the name of the analysis.

By default this is computed by combining the results of the experiment, year and Spires ID metadata methods and you should only override it if there's a good reason why those won't work.

Implements ProjectionApplier.

Definition at line 118 of file Analysis.hh.

References Analysis::_defaultname, Analysis::info(), and AnalysisInfo::name().

Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::bookDataPointSet(), Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), Analysis::crossSection(), Analysis::getLog(), Analysis::histoDir(), Analysis::normalize(), HistoHandler::registerAnalysisObject(), and Analysis::scale().

                                   {
      return (info().name().empty()) ? _defaultname : info().name();
    }

bool needsCrossSection

(

)

const [inline, inherited]

Return true if this analysis needs to know the process cross-section.

Definition at line 232 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::needsCrossSection().

                                   {
      return info().needsCrossSection();
    }

void normalize	(	AIDA::IHistogram2D *&	histo,
		double	norm = 1.0
	)			[inherited]

Normalize the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path. It has the same effect as scale(histo, norm/sumOfWeights).

Parameters:

	histo	The histogram to be normalised.
	norm	The new area of the histogram.

Warning:

The old histogram will be deleted, and its pointer set to zero.

Definition at line 552 of file Analysis.cc.

References MSG_ERROR, MSG_TRACE, MSG_WARNING, Analysis::name(), Analysis::scale(), and Analysis::tree().

                                                                {
    if (!histo) {
      MSG_ERROR("Failed to normalize histo=NULL in analysis "
                << name() << " (norm=" << norm << ")");
      return;
    }
    const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
    MSG_TRACE("Normalizing histo " << hpath << " to " << norm);

    double oldintg = 0.0;
    int nxBins = histo->xAxis().bins();
    int nyBins = histo->yAxis().bins();
    for (int ixBin = 0; ixBin != nxBins; ++ixBin)
      for (int iyBin = 0; iyBin != nyBins; ++iyBin) {
      // Leaving out factor of binWidth because AIDA's "height"
      // already includes a width factor.
    oldintg += histo->binHeight(ixBin, iyBin); // * histo->axis().binWidth(iBin);
    }
    if (oldintg == 0.0) {
      MSG_WARNING("Histo " << hpath << " has null integral during normalization");
      return;
    }

    // Scale by the normalisation factor.
    scale(histo, norm/oldintg);
  }

void normalize	(	AIDA::IHistogram1D *&	histo,
		double	norm = 1.0
	)			[inherited]

Normalize the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path. It has the same effect as scale(histo, norm/sumOfWeights).

Parameters:

	histo	The histogram to be normalised.
	norm	The new area of the histogram.

Warning:

The old histogram will be deleted, and its pointer set to zero.

Definition at line 482 of file Analysis.cc.

References MSG_ERROR, MSG_TRACE, MSG_WARNING, Analysis::name(), Analysis::scale(), and Analysis::tree().

Referenced by TASSO_1990_S2148048::finalize(), STAR_2008_S7869363::finalize(), OPAL_2001_S4553896::finalize(), OPAL_1998_S3780481::finalize(), JADE_1998_S3612880::finalize(), H1_1994_S2919893::finalize(), ExampleAnalysis::finalize(), DELPHI_2003_WUD_03_11::finalize(), DELPHI_2002_069_CONF_603::finalize(), DELPHI_1996_S3430090::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8821313::finalize(), D0_2010_S8671338::finalize(), D0_2009_S8320160::finalize(), D0_2008_S7554427::finalize(), D0_2008_S6879055::finalize(), D0_2007_S7075677::finalize(), D0_2004_S5992206::finalize(), D0_2001_S4674421::finalize(), D0_1996_S3324664::finalize(), D0_1996_S3214044::finalize(), CMS_2011_S9120041::finalize(), CMS_2011_S8968497::finalize(), CMS_2011_S8957746::finalize(), CMS_2011_S8950903::finalize(), CMS_2011_S8884919::finalize(), CDF_2009_NOTE_9936::finalize(), CDF_2004_S5839831::finalize(), CDF_2001_S4751469::finalize(), CDF_1997_S3541940::finalize(), CDF_1996_S3418421::finalize(), CDF_1996_S3349578::finalize(), CDF_1996_S3108457::finalize(), CDF_1994_S2952106::finalize(), CDF_1993_S2742446::finalize(), ATLAS_2011_S8971293::finalize(), ALICE_2010_S8625980::finalize(), ALICE_2010_S8624100::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), and ATLAS_2011_I919017::safeinvscale().

                                                                {
    if (!histo) {
      MSG_ERROR("Failed to normalize histo=NULL in analysis "
                << name() << " (norm=" << norm << ")");
      return;
    }
    const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
    MSG_TRACE("Normalizing histo " << hpath << " to " << norm);

    double oldintg = 0.0;
    int nBins = histo->axis().bins();
    for (int iBin = 0; iBin != nBins; ++iBin) {
      // Leaving out factor of binWidth because AIDA's "height" already includes a width factor.
      oldintg += histo->binHeight(iBin); // * histo->axis().binWidth(iBin);
    }
    if (oldintg == 0.0) {
      MSG_WARNING("Histo " << hpath << " has null integral during normalization");
      return;
    }

    // Scale by the normalisation factor.
    scale(histo, norm/oldintg);
  }

size_t numEvents

(

)

const [protected, inherited]

Get the number of events seen (via the analysis handler). Use in the finalize phase only.

Definition at line 103 of file Analysis.cc.

References Analysis::handler(), and AnalysisHandler::numEvents().

Referenced by LHCB_2010_S8758301::finalize().

                                   {
    return handler().numEvents();
  }

virtual std::vector<std::string> references

(

)

const [inline, virtual, inherited]

Journal, and preprint references.

Definition at line 184 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::references().

                                                    {
      return info().references();
    }

virtual const std::vector<PdgIdPair>& requiredBeams

(

)

const [inline, virtual, inherited]

Return the allowed pairs of incoming beams required by this analysis.

Definition at line 210 of file Analysis.hh.

References AnalysisInfo::beams(), and Analysis::info().

                                                              {
      return info().beams();
    }

virtual const std::vector<std::pair<double, double> >& requiredEnergies

(

)

const [inline, virtual, inherited]

Sets of valid beam energy pairs, in GeV.

Definition at line 221 of file Analysis.hh.

References AnalysisInfo::energies(), and Analysis::info().

Referenced by Analysis::setRequiredEnergies().

                                                                                {
      return info().energies();
    }

virtual std::string runInfo

(

)

const [inline, virtual, inherited]

Information about the events needed as input for this analysis.

Event types, energies, kinematic cuts, particles to be considered stable, etc. etc. Should be treated as a restructuredText bullet list (http://docutils.sourceforge.net/rst.html)

Definition at line 164 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::runInfo().

                                      {
      return info().runInfo();
    }

void scale	(	AIDA::IHistogram2D *&	histo,
		double	scale
	)			[inherited]

Multiplicatively scale the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path.

Parameters:

	histo	The histogram to be scaled.
	scale	The factor used to multiply the histogram bin heights.

Warning:

The old histogram will be deleted, and its pointer set to zero.

Definition at line 580 of file Analysis.cc.

References Analysis::datapointsetFactory(), MSG_ERROR, MSG_TRACE, Analysis::name(), and Analysis::tree().

                                                             {
    if (!histo) {
      MSG_ERROR("Failed to scale histo=NULL in analysis "
                << name() << " (scale=" << scale << ")");
      return;
    }
    const string hpath =
      tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
    MSG_TRACE("Scaling histo " << hpath);

    vector<double> x, y, z, ex, ey, ez;
    for (size_t ix = 0, Nx = histo->xAxis().bins(); ix < Nx; ++ix)
      for (size_t iy = 0, Ny = histo->yAxis().bins(); iy < Ny; ++iy) {
    x.push_back(0.5 * (histo->xAxis().binLowerEdge(ix) +
               histo->xAxis().binUpperEdge(ix)));
    ex.push_back(histo->xAxis().binWidth(ix)*0.5);
    y.push_back(0.5 * (histo->yAxis().binLowerEdge(iy) +
               histo->yAxis().binUpperEdge(iy)));
    ey.push_back(histo->yAxis().binWidth(iy)*0.5);

    // "Bin height" is a misnomer in the AIDA spec: width is neglected.
    // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
    z.push_back(histo->binHeight(ix, iy)*scale/
            (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy)));
    // "Bin error" is a misnomer in the AIDA spec: width is neglected.
    // We'd like to do this: ey.push_back(histo->binError(i) * scale);
    ez.push_back(histo->binError(ix, iy)*scale/
             (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy)));
    }

    string title = histo->title();
    string xtitle = histo->xtitle();
    string ytitle = histo->ytitle();
    string ztitle = histo->ztitle();

    tree().mkdir("/tmpnormalize");
    tree().mv(hpath, "/tmpnormalize");

    if (hpath.find(" ") != string::npos) {
      throw Error("Histogram path '" + hpath + "' is invalid: spaces are not permitted in paths");
    }
    AIDA::IDataPointSet* dps =
      datapointsetFactory().createXYZ(hpath, title, x, y, z, ex, ey, ez);
    dps->setXTitle(xtitle);
    dps->setYTitle(ytitle);
    dps->setZTitle(ztitle);

    tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
    tree().rmdir("/tmpnormalize");

    // Set histo pointer to null - it can no longer be used.
    histo = 0;
  }

void scale	(	AIDA::IHistogram1D *&	histo,
		double	scale
	)			[inherited]

Multiplicatively scale the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path.

Parameters:

	histo	The histogram to be scaled.
	scale	The factor used to multiply the histogram bin heights.

Warning:

The old histogram will be deleted, and its pointer set to zero.

Definition at line 507 of file Analysis.cc.

References Analysis::datapointsetFactory(), MSG_ERROR, MSG_TRACE, Analysis::name(), and Analysis::tree().

Referenced by UA5_1989_S1926373::finalize(), UA5_1987_S1640666::finalize(), UA5_1986_S1583476::finalize(), UA5_1982_S875503::finalize(), UA1_1990_S2044935::finalize(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), SFM_1984_S1178091::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_2004_S6132243::finalize(), OPAL_1998_S3780481::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_XS::finalize(), MC_WWJETS::finalize(), MC_WPOL::finalize(), MC_WJETS::finalize(), MC_VH2BB::finalize(), MC_TTBAR::finalize(), MC_PHOTONJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), LHCB_2010_S8758301::finalize(), JADE_OPAL_2000_S4300807::finalize(), H1_2000_S4129130::finalize(), H1_1994_S2919893::finalize(), E735_1998_S3905616::finalize(), DELPHI_1996_S3430090::finalize(), D0_2010_S8671338::finalize(), D0_2010_S8570965::finalize(), D0_2009_S8349509::finalize(), D0_2009_S8202443::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S8978280::finalize(), CMS_2010_S8656010::finalize(), CMS_2010_S8547297::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7540469::finalize(), CDF_2007_S7057202::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), CDF_1990_S2089246::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2011_S9131140::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_S9108483::finalize(), ATLAS_2011_S9041966::finalize(), ATLAS_2011_I925932::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2010_S8591806::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ALICE_2011_S8945144::finalize(), ALICE_2011_S8909580::finalize(), ALICE_2010_S8706239::finalize(), ALICE_2010_S8625980::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), ALEPH_1991_S2435284::finalize(), Analysis::normalize(), ATLAS_2011_I919017::safeinvscale(), and BinnedHistogram< T >::scale().

                                                             {
    if (!histo) {
      MSG_ERROR("Failed to scale histo=NULL in analysis "
                << name() << " (scale=" << scale << ")");
      return;
    }
    const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
    MSG_TRACE("Scaling histo " << hpath);

    vector<double> x, y, ex, ey;
    for (size_t i = 0, N = histo->axis().bins(); i < N; ++i) {
      x.push_back(0.5 * (histo->axis().binLowerEdge(i) + histo->axis().binUpperEdge(i)));
      ex.push_back(histo->axis().binWidth(i)*0.5);

      // "Bin height" is a misnomer in the AIDA spec: width is neglected.
      // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
      y.push_back(histo->binHeight(i)*scale/histo->axis().binWidth(i));

      // "Bin error" is a misnomer in the AIDA spec: width is neglected.
      // We'd like to do this: ey.push_back(histo->binError(i) * scale);
      ey.push_back(histo->binError(i)*scale/histo->axis().binWidth(i));
    }

    string title = histo->title();
    string xtitle = histo->xtitle();
    string ytitle = histo->ytitle();

    tree().mkdir("/tmpnormalize");
    tree().mv(hpath, "/tmpnormalize");

    if (hpath.find(" ") != string::npos) {
      throw Error("Histogram path '" + hpath + "' is invalid: spaces are not permitted in paths");
    }
    AIDA::IDataPointSet* dps = datapointsetFactory().createXY(hpath, title, x, y, ex, ey);
    dps->setXTitle(xtitle);
    dps->setYTitle(ytitle);

    tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
    tree().rmdir("/tmpnormalize");

    // Set histo pointer to null - it can no longer be used.
    histo = 0;
  }

virtual Analysis& setBeams	(	PdgId	beam1,
		PdgId	beam2
	)			[inline, virtual, inherited]

Set the required beams

Deprecated:

To be removed in 2.0.0. Use .info file and AnalysisInfo class instead

Todo:

Print out a warning to use setRequiredBeams() instead (and really to use .info files)

Definition at line 255 of file Analysis.hh.

References Analysis::setRequiredBeams().

                                                         {
      /// @todo Print out a warning to use setRequiredBeams() instead (and really to use .info files)
      return setRequiredBeams(std::vector<PdgIdPair>(1, make_pair(beam1, beam2)));
    }

Analysis & setCrossSection

(

double

xs

)

[inherited]

Set the cross section from the generator.

Definition at line 160 of file Analysis.cc.

References Analysis::_crossSection, and Analysis::_gotCrossSection.

                                               {
    _crossSection = xs;
    _gotCrossSection = true;
    return *this;
  }

Analysis& setNeedsCrossSection

(

bool

needed = true

)

[inline, inherited]

Declare whether this analysis needs to know the process cross-section from the generator.

Definition at line 236 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::setNeedsCrossSection().

Referenced by ATLAS_2011_S9108483::ATLAS_2011_S9108483(), LHCB_2010_S8758301::LHCB_2010_S8758301(), and MC_JetAnalysis::MC_JetAnalysis().

                                                     {
      info().setNeedsCrossSection(needed);
      return *this;
    }

virtual Analysis& setRequiredBeams

(

const std::vector< PdgIdPair > &

requiredBeams

)

[inline, virtual, inherited]

Declare the allowed pairs of incoming beams required by this analysis.

Definition at line 214 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::setBeams().

Referenced by Analysis::setBeams().

                                                                                  {
      info().setBeams(requiredBeams);
      return *this;
    }

virtual Analysis& setRequiredEnergies

(

const std::vector< std::pair< double, double > > &

requiredEnergies

)

[inline, virtual, inherited]

Declare the list of valid beam energy pairs, in GeV.

Definition at line 225 of file Analysis.hh.

References Analysis::info(), Analysis::requiredEnergies(), and AnalysisInfo::setEnergies().

                                                                                                       {
      info().setEnergies(requiredEnergies);
      return *this;
    }

virtual std::string spiresId

(

)

const [inline, virtual, inherited]

Get the SPIRES ID code for this analysis (~deprecated).

Definition at line 128 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::spiresId().

                                       {
      return info().spiresId();
    }

double sqrtS

(

)

const [inherited]

Centre of mass energy for this run.

Definition at line 50 of file Analysis.cc.

References Analysis::handler(), and AnalysisHandler::sqrtS().

Referenced by UA1_1990_S2044935::analyze(), TASSO_1990_S2148048::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), MC_TTBAR::analyze(), JADE_1998_S3612880::analyze(), CMS_2011_S9120041::analyze(), CDF_2004_S5839831::analyze(), ATLAS_2010_S8918562::analyze(), ATLAS_2010_S8894728::analyze(), ALICE_2010_S8625980::analyze(), ALEPH_2004_S5765862::analyze(), UA1_1990_S2044935::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), JADE_1998_S3612880::finalize(), CMS_2011_S8978280::finalize(), CDF_2004_S5839831::finalize(), ATLAS_2010_S8918562::finalize(), ALICE_2010_S8625980::finalize(), ALEPH_2004_S5765862::finalize(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1986_S1583476::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), MC_DIPHOTON::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), CMS_2011_S8978280::init(), CMS_2011_S8884919::init(), CMS_2010_S8547297::init(), CDF_2004_S5839831::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), ATLAS_2011_S8994773::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8894728::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), and ALEPH_2004_S5765862::init().

                               {
    return handler().sqrtS();
  }

virtual std::string status

(

)

const [inline, virtual, inherited]

Whether this analysis is trusted (in any way!).

Definition at line 199 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::status().

                                     {
      return (info().status().empty()) ? "UNVALIDATED" : info().status();
    }

virtual std::string summary

(

)

const [inline, virtual, inherited]

Get a short description of the analysis.

Short (one sentence) description used as an index entry. Use description() to provide full descriptive paragraphs of analysis details.

Definition at line 145 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::summary().

                                      {
      return info().summary();
    }

double sumOfWeights

(

)

const [protected, inherited]

Get the sum of event weights seen (via the analysis handler). Use in the finalize phase only.

Definition at line 108 of file Analysis.cc.

References Analysis::handler(), and AnalysisHandler::sumOfWeights().

Referenced by Analysis::crossSectionPerEvent(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_1998_S3780481::finalize(), OPAL_1993_S2692198::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_XS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_VH2BB::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), LHCB_2010_S8758301::finalize(), JADE_OPAL_2000_S4300807::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S9086218::finalize(), CMS_2011_S8978280::finalize(), CMS_2010_S8656010::finalize(), CMS_2010_S8547297::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7541902::finalize(), CDF_2008_S7540469::finalize(), CDF_2006_S6653332::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2001_S4517016::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_S9108483::finalize(), ATLAS_2011_S9041966::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ALICE_2011_S8945144::finalize(), ALICE_2011_S8909580::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), and ALEPH_1991_S2435284::finalize().

                                      {
    return handler().sumOfWeights();
  }

virtual std::vector<std::string> todos

(

)

const [inline, virtual, inherited]

Any work to be done on this analysis.

Definition at line 204 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::todos().

                                               {
      return info().todos();
    }

ITree & tree

(

)

[protected, inherited]

Access the AIDA tree of the controlling AnalysisHandler object.

Definition at line 35 of file Analysis.cc.

References Analysis::handler(), and AnalysisHandler::tree().

Referenced by Analysis::_makeHistoDir(), Analysis::normalize(), and Analysis::scale().

                        {
    return handler().tree();
  }

virtual std::string year

(

)

const [inline, virtual, inherited]

When the original experimental analysis was published.

Definition at line 179 of file Analysis.hh.

References Analysis::info(), and AnalysisInfo::year().

                                   {
      return info().year();
    }

---

Member Data Documentation

bool _allowProjReg [protected, inherited]

Flag to forbid projection registration in analyses until the init phase.

Definition at line 140 of file ProjectionApplier.hh.

Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().

string _defaultname [protected, inherited]

Name passed to constructor (used to find .info analysis data file, and as a fallback).

Definition at line 547 of file Analysis.hh.

Referenced by Analysis::Analysis(), and Analysis::name().

AIDA::IHistogram1D * _hist_eta_e [private]

Definition at line 307 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_eta_gamma [private]

Definition at line 309 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_eta_gammaiso [private]

Definition at line 310 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_eta_jet [private]

Definition at line 306 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_eta_mu [private]

Definition at line 308 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_eta_trk [private]

Definition at line 305 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_met [private]

Definition at line 311 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_2_osof_emu [private]

Definition at line 312 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_mll_2_ossf_ee [private]

Definition at line 312 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_2_ossf_mumu [private]

Definition at line 312 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_all_osof_emu [private]

Definition at line 314 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_mll_all_ossf_ee [private]

Definition at line 314 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_all_ossf_mumu [private]

Definition at line 314 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_osof_emu [private]

Definition at line 313 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_mll_ossf_ee [private]

Definition at line 313 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_mll_ossf_mumu [private]

Definition at line 313 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_e [private]

Definition at line 307 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_gamma [private]

Definition at line 309 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_gammaiso [private]

Definition at line 310 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_jet [private]

Definition at line 306 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_mu [private]

Definition at line 308 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D* _hist_n_trk [private]

Definition at line 305 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_e [private]

Definition at line 307 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_gamma [private]

Definition at line 309 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_gammaiso [private]

Definition at line 310 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_jet [private]

Definition at line 306 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_mu [private]

Definition at line 308 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_phi_trk [private]

Definition at line 305 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_e [private]

Definition at line 307 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_gamma [private]

Definition at line 309 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_gammaiso [private]

Definition at line 310 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_jet [private]

Definition at line 306 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_mu [private]

Definition at line 308 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

AIDA::IHistogram1D * _hist_pt_trk [private]

Definition at line 305 of file MC_SUSY.cc.

Referenced by MC_SUSY::analyze(), and MC_SUSY::init().

shared_ptr<AnalysisInfo> _info [protected, inherited]

Pointer to analysis metadata object.

Definition at line 550 of file Analysis.hh.

Referenced by Analysis::Analysis(), and Analysis::info().

---

The documentation for this class was generated from the following file:

• MC_SUSY.cc