ALEPH_1996_S3486095 Class Reference

ALEPH QCD study with event shapes and identified particles. More...

Inheritance diagram for ALEPH_1996_S3486095:
Inheritance graph
[legend]
Collaboration diagram for ALEPH_1996_S3486095:
Collaboration graph
[legend]

List of all members.

Public Member Functions

 ALEPH_1996_S3486095 ()
 Constructor.
AnalysisHandlerhandler () 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)
AnalysissetCrossSection (double xs)
 Set the cross section from the generator.
Analysis methods

void init ()
void analyze (const Event &e)
void finalize ()
 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 AnalysisInfoinfo () 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 spiresId () const
 Get a the SPIRES/Inspire ID code for this analysis.
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 AnalysissetRequiredBeams (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 AnalysissetRequiredEnergies (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.
AnalysissetNeedsCrossSection (bool needed=true)
 Declare whether this analysis needs to know the process cross-section from the generator.
Internal metadata modifiying methods

AnalysisInfoinfo ()
 Get the actual AnalysisInfo object in which all this metadata is stored (non-const).
virtual AnalysissetBeams (PdgId beam1, PdgId beam2)
Run conditions

const ParticlePairbeams () 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< ConstProjectionPtrgetProjections () 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 ProjectiongetProjection (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

LoggetLog () 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
ProjectionHandlergetProjHandler () 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.
Internal histogram booking (for use by Analysis sub-classes).

const BinEdgesbinEdges (const std::string &hname) const
 Get bin edges for a named histo (using ref AIDA caching).
const BinEdgesbinEdges (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

double _weightedTotalPartNum
double _weightedTotalNumPiPlus
double _weightedTotalNumKPlus
double _weightedTotalNumP
double _weightedTotalNumPhoton
double _weightedTotalNumPi0
double _weightedTotalNumEta
double _weightedTotalNumEtaPrime
double _weightedTotalNumK0
double _weightedTotalNumLambda0
double _weightedTotalNumXiMinus
double _weightedTotalNumSigma1385Plus
double _weightedTotalNumXi1530_0
double _weightedTotalNumRho
double _weightedTotalNumOmega782
double _weightedTotalNumKStar892_0
double _weightedTotalNumPhi
double _weightedTotalNumKStar892Plus
double _numChParticles
Histograms

AIDA::IHistogram1D * _histSphericity
AIDA::IHistogram1D * _histAplanarity
AIDA::IHistogram1D * _hist1MinusT
AIDA::IHistogram1D * _histTMinor
AIDA::IHistogram1D * _histY3
AIDA::IHistogram1D * _histHeavyJetMass
AIDA::IHistogram1D * _histCParam
AIDA::IHistogram1D * _histOblateness
AIDA::IHistogram1D * _histScaledMom
AIDA::IHistogram1D * _histRapidityT
AIDA::IHistogram1D * _histPtSIn
AIDA::IHistogram1D * _histPtSOut
AIDA::IHistogram1D * _histJetRate2Durham
AIDA::IHistogram1D * _histJetRate3Durham
AIDA::IHistogram1D * _histJetRate4Durham
AIDA::IHistogram1D * _histJetRate5Durham
AIDA::IHistogram1D * _histLogScaledMom
AIDA::IHistogram1D * _histChMult
AIDA::IHistogram1D * _histMultiPiPlus
AIDA::IHistogram1D * _histMultiKPlus
AIDA::IHistogram1D * _histMultiP
AIDA::IHistogram1D * _histMultiPhoton
AIDA::IHistogram1D * _histMultiPi0
AIDA::IHistogram1D * _histMultiEta
AIDA::IHistogram1D * _histMultiEtaPrime
AIDA::IHistogram1D * _histMultiK0
AIDA::IHistogram1D * _histMultiLambda0
AIDA::IHistogram1D * _histMultiXiMinus
AIDA::IHistogram1D * _histMultiSigma1385Plus
AIDA::IHistogram1D * _histMultiXi1530_0
AIDA::IHistogram1D * _histMultiRho
AIDA::IHistogram1D * _histMultiOmega782
AIDA::IHistogram1D * _histMultiKStar892_0
AIDA::IHistogram1D * _histMultiPhi
AIDA::IHistogram1D * _histMultiKStar892Plus
AIDA::IHistogram1D * _histMeanChMult
AIDA::IHistogram1D * _histMeanChMultRapt05
AIDA::IHistogram1D * _histMeanChMultRapt10
AIDA::IHistogram1D * _histMeanChMultRapt15
AIDA::IHistogram1D * _histMeanChMultRapt20
AIDA::IHistogram1D * _histMeanMultiPi0
AIDA::IHistogram1D * _histMeanMultiEta
AIDA::IHistogram1D * _histMeanMultiEtaPrime
AIDA::IHistogram1D * _histMeanMultiK0
AIDA::IHistogram1D * _histMeanMultiRho
AIDA::IHistogram1D * _histMeanMultiOmega782
AIDA::IHistogram1D * _histMeanMultiPhi
AIDA::IHistogram1D * _histMeanMultiKStar892Plus
AIDA::IHistogram1D * _histMeanMultiKStar892_0
AIDA::IHistogram1D * _histMeanMultiLambda0
AIDA::IHistogram1D * _histMeanMultiSigma0
AIDA::IHistogram1D * _histMeanMultiXiMinus
AIDA::IHistogram1D * _histMeanMultiSigma1385Plus
AIDA::IHistogram1D * _histMeanMultiXi1530_0
AIDA::IHistogram1D * _histMeanMultiOmegaOmegaBar

Detailed Description

ALEPH QCD study with event shapes and identified particles.

Author:
Holger Schulz

Definition at line 20 of file ALEPH_1996_S3486095.cc.


Constructor & Destructor Documentation

ALEPH_1996_S3486095 (  )  [inline]

Constructor.

Definition at line 25 of file ALEPH_1996_S3486095.cc.

References ALEPH_1996_S3486095::_numChParticles, ALEPH_1996_S3486095::_weightedTotalNumEta, ALEPH_1996_S3486095::_weightedTotalNumEtaPrime, ALEPH_1996_S3486095::_weightedTotalNumK0, ALEPH_1996_S3486095::_weightedTotalNumKPlus, ALEPH_1996_S3486095::_weightedTotalNumKStar892_0, ALEPH_1996_S3486095::_weightedTotalNumKStar892Plus, ALEPH_1996_S3486095::_weightedTotalNumLambda0, ALEPH_1996_S3486095::_weightedTotalNumOmega782, ALEPH_1996_S3486095::_weightedTotalNumP, ALEPH_1996_S3486095::_weightedTotalNumPhi, ALEPH_1996_S3486095::_weightedTotalNumPhoton, ALEPH_1996_S3486095::_weightedTotalNumPi0, ALEPH_1996_S3486095::_weightedTotalNumPiPlus, ALEPH_1996_S3486095::_weightedTotalNumRho, ALEPH_1996_S3486095::_weightedTotalNumSigma1385Plus, ALEPH_1996_S3486095::_weightedTotalNumXi1530_0, ALEPH_1996_S3486095::_weightedTotalNumXiMinus, ALEPH_1996_S3486095::_weightedTotalPartNum, Rivet::ELECTRON, Rivet::POSITRON, and Analysis::setBeams().


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().

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

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(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), FinalStateHCM::FinalStateHCM(), 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_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_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::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_2002_S4796047::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_S9120807::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_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_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(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), Spherocity::Spherocity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), and VisibleFinalState::VisibleFinalState().

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

IAnalysisFactory & analysisFactory (  )  [protected, inherited]

Access the AIDA analysis factory of the controlling AnalysisHandler object.

Definition at line 50 of file Analysis.cc.

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

00050                                               {
00051     return handler().analysisFactory();
00052   }

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.

Implements Analysis.

Definition at line 135 of file ALEPH_1996_S3486095.cc.

References ALEPH_1996_S3486095::_hist1MinusT, ALEPH_1996_S3486095::_histAplanarity, ALEPH_1996_S3486095::_histChMult, ALEPH_1996_S3486095::_histCParam, ALEPH_1996_S3486095::_histHeavyJetMass, ALEPH_1996_S3486095::_histLogScaledMom, ALEPH_1996_S3486095::_histMeanChMult, ALEPH_1996_S3486095::_histMeanChMultRapt05, ALEPH_1996_S3486095::_histMeanChMultRapt10, ALEPH_1996_S3486095::_histMeanChMultRapt15, ALEPH_1996_S3486095::_histMeanChMultRapt20, ALEPH_1996_S3486095::_histMeanMultiEta, ALEPH_1996_S3486095::_histMeanMultiEtaPrime, ALEPH_1996_S3486095::_histMeanMultiK0, ALEPH_1996_S3486095::_histMeanMultiKStar892_0, ALEPH_1996_S3486095::_histMeanMultiKStar892Plus, ALEPH_1996_S3486095::_histMeanMultiLambda0, ALEPH_1996_S3486095::_histMeanMultiOmega782, ALEPH_1996_S3486095::_histMeanMultiOmegaOmegaBar, ALEPH_1996_S3486095::_histMeanMultiPhi, ALEPH_1996_S3486095::_histMeanMultiPi0, ALEPH_1996_S3486095::_histMeanMultiRho, ALEPH_1996_S3486095::_histMeanMultiSigma0, ALEPH_1996_S3486095::_histMeanMultiSigma1385Plus, ALEPH_1996_S3486095::_histMeanMultiXi1530_0, ALEPH_1996_S3486095::_histMeanMultiXiMinus, ALEPH_1996_S3486095::_histMultiEta, ALEPH_1996_S3486095::_histMultiEtaPrime, ALEPH_1996_S3486095::_histMultiK0, ALEPH_1996_S3486095::_histMultiKPlus, ALEPH_1996_S3486095::_histMultiKStar892_0, ALEPH_1996_S3486095::_histMultiKStar892Plus, ALEPH_1996_S3486095::_histMultiLambda0, ALEPH_1996_S3486095::_histMultiOmega782, ALEPH_1996_S3486095::_histMultiP, ALEPH_1996_S3486095::_histMultiPhi, ALEPH_1996_S3486095::_histMultiPhoton, ALEPH_1996_S3486095::_histMultiPi0, ALEPH_1996_S3486095::_histMultiPiPlus, ALEPH_1996_S3486095::_histMultiRho, ALEPH_1996_S3486095::_histMultiSigma1385Plus, ALEPH_1996_S3486095::_histMultiXi1530_0, ALEPH_1996_S3486095::_histMultiXiMinus, ALEPH_1996_S3486095::_histOblateness, ALEPH_1996_S3486095::_histPtSIn, ALEPH_1996_S3486095::_histPtSOut, ALEPH_1996_S3486095::_histRapidityT, ALEPH_1996_S3486095::_histScaledMom, ALEPH_1996_S3486095::_histSphericity, ALEPH_1996_S3486095::_histTMinor, ALEPH_1996_S3486095::_histY3, ALEPH_1996_S3486095::_numChParticles, ALEPH_1996_S3486095::_weightedTotalNumEta, ALEPH_1996_S3486095::_weightedTotalNumEtaPrime, ALEPH_1996_S3486095::_weightedTotalNumK0, ALEPH_1996_S3486095::_weightedTotalNumKPlus, ALEPH_1996_S3486095::_weightedTotalNumKStar892_0, ALEPH_1996_S3486095::_weightedTotalNumKStar892Plus, ALEPH_1996_S3486095::_weightedTotalNumLambda0, ALEPH_1996_S3486095::_weightedTotalNumOmega782, ALEPH_1996_S3486095::_weightedTotalNumP, ALEPH_1996_S3486095::_weightedTotalNumPhi, ALEPH_1996_S3486095::_weightedTotalNumPhoton, ALEPH_1996_S3486095::_weightedTotalNumPi0, ALEPH_1996_S3486095::_weightedTotalNumPiPlus, ALEPH_1996_S3486095::_weightedTotalNumRho, ALEPH_1996_S3486095::_weightedTotalNumSigma1385Plus, ALEPH_1996_S3486095::_weightedTotalNumXi1530_0, ALEPH_1996_S3486095::_weightedTotalNumXiMinus, ALEPH_1996_S3486095::_weightedTotalPartNum, Sphericity::aplanarity(), Analysis::beams(), ParisiTensor::C(), FastJets::clusterSeq(), Log::DEBUG, Rivet::dot(), FourMomentum::E(), Analysis::getLog(), Rivet::GeV, Vector< N >::mod(), Particle::momentum(), Thrust::oblateness(), UnstableFinalState::particles(), FinalState::particles(), Particle::pdgId(), Hemispheres::scaledM2high(), Sphericity::sphericity(), Sphericity::sphericityMajorAxis(), Sphericity::sphericityMinorAxis(), Thrust::thrust(), Thrust::thrustAxis(), Thrust::thrustMinor(), FourVector::vector3(), and vetoEvent.

00135                                  {
00136       // First, veto on leptonic events by requiring at least 4 charged FS particles
00137       const FinalState& fs = applyProjection<FinalState>(e, "FS");
00138       const size_t numParticles = fs.particles().size();
00139 
00140       // Even if we only generate hadronic events, we still need a cut on numCharged >= 2.
00141       if (numParticles < 2) {
00142         getLog() << Log::DEBUG << "Failed leptonic event cut" << endl;
00143         vetoEvent;
00144       }
00145       getLog() << Log::DEBUG << "Passed leptonic event cut" << endl;
00146 
00147       // Get event weight for histo filling
00148       const double weight = e.weight();
00149       _weightedTotalPartNum += numParticles * weight;
00150 
00151       // Get beams and average beam momentum
00152       const ParticlePair& beams = applyProjection<Beam>(e, "Beams").beams();
00153       const double meanBeamMom = ( beams.first.momentum().vector3().mod() +
00154                                    beams.second.momentum().vector3().mod() ) / 2.0;
00155       getLog() << Log::DEBUG << "Avg beam momentum = " << meanBeamMom << endl;
00156 
00157       // Thrusts
00158       getLog() << Log::DEBUG << "Calculating thrust" << endl;
00159       const Thrust& thrust = applyProjection<Thrust>(e, "Thrust");
00160       _hist1MinusT->fill(1 - thrust.thrust(), weight);
00161       _histTMinor->fill(thrust.thrustMinor(), weight);
00162       _histOblateness->fill(thrust.oblateness(), weight);
00163 
00164       // Jets
00165       getLog() << Log::DEBUG << "Calculating differential jet rate plots:" << endl;
00166       const FastJets& durjet = applyProjection<FastJets>(e, "DurhamJets");
00167       if (durjet.clusterSeq()) {
00168         double y3 = durjet.clusterSeq()->exclusive_ymerge_max(2);
00169         _histY3->fill(-1. * std::log(y3), weight);
00170       }
00171 
00172       // Sphericities
00173       getLog() << Log::DEBUG << "Calculating sphericity" << endl;
00174       const Sphericity& sphericity = applyProjection<Sphericity>(e, "Sphericity");
00175       _histSphericity->fill(sphericity.sphericity(), weight);
00176       _histAplanarity->fill(sphericity.aplanarity(), weight);
00177 
00178       // C param
00179       getLog() << Log::DEBUG << "Calculating Parisi params" << endl;
00180       const ParisiTensor& parisi = applyProjection<ParisiTensor>(e, "Parisi");
00181       _histCParam->fill(parisi.C(), weight);
00182 
00183       // Hemispheres
00184       getLog() << Log::DEBUG << "Calculating hemisphere variables" << endl;
00185       const Hemispheres& hemi = applyProjection<Hemispheres>(e, "Hemispheres");
00186       _histHeavyJetMass->fill(hemi.scaledM2high(), weight);
00187 
00188       // Iterate over all the charged final state particles.
00189       double Evis = 0.0;
00190       double rapt05 = 0.;
00191       double rapt10 = 0.;
00192       double rapt15 = 0.;
00193       double rapt20 = 0.;
00194       //int numChParticles = 0;
00195       getLog() << Log::DEBUG << "About to iterate over charged FS particles" << endl;
00196       foreach (const Particle& p, fs.particles()) {
00197         // Get momentum and energy of each particle.
00198         const Vector3 mom3 = p.momentum().vector3();
00199         const double energy = p.momentum().E();
00200         Evis += energy;
00201         _numChParticles += weight;
00202 
00203         // Scaled momenta.
00204         const double mom = mom3.mod();
00205         const double scaledMom = mom/meanBeamMom;
00206         const double logInvScaledMom = -std::log(scaledMom);
00207         _histLogScaledMom->fill(logInvScaledMom, weight);
00208         _histScaledMom->fill(scaledMom, weight);
00209 
00210         // Get momenta components w.r.t. thrust and sphericity.
00211         const double momT = dot(thrust.thrustAxis(), mom3);
00212         const double pTinS = dot(mom3, sphericity.sphericityMajorAxis());
00213         const double pToutS = dot(mom3, sphericity.sphericityMinorAxis());
00214         _histPtSIn->fill(fabs(pTinS/GeV), weight);
00215         _histPtSOut->fill(fabs(pToutS/GeV), weight);
00216 
00217         // Calculate rapidities w.r.t. thrust.
00218         const double rapidityT = 0.5 * std::log((energy + momT) / (energy - momT));
00219         _histRapidityT->fill(rapidityT, weight);
00220         if (std::fabs(rapidityT) <= 0.5)  {
00221             rapt05 += 1.0;
00222         }
00223         if (std::fabs(rapidityT) <= 1.0)  {
00224             rapt10 += 1.0;
00225         }
00226         if (std::fabs(rapidityT) <= 1.5) {
00227             rapt15 += 1.0;
00228         }
00229         if (std::fabs(rapidityT) <= 2.0)  {
00230             rapt20 += 1.0;
00231         }
00232 
00233       }
00234 
00235       _histChMult->fill(numParticles, weight);
00236 
00237       _histMeanChMultRapt05->fill(_histMeanChMultRapt05->binMean(0), rapt05 * weight);
00238       _histMeanChMultRapt10->fill(_histMeanChMultRapt10->binMean(0), rapt10 * weight);
00239       _histMeanChMultRapt15->fill(_histMeanChMultRapt15->binMean(0), rapt15 * weight);
00240       _histMeanChMultRapt20->fill(_histMeanChMultRapt20->binMean(0), rapt20 * weight);
00241       _histMeanChMult->fill(_histMeanChMult->binMean(0), numParticles*weight);
00242 
00243 
00244       //// Final state of unstable particles to get particle spectra
00245       const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(e, "UFS");
00246       for (ParticleVector::const_iterator p = ufs.particles().begin(); p != ufs.particles().end(); ++p) {
00247         const Vector3 mom3 = p->momentum().vector3();
00248         int id = abs(p->pdgId());
00249         const double mom = mom3.mod();
00250         const double energy = p->momentum().E();
00251         const double scaledMom = mom/meanBeamMom;
00252         const double scaledEnergy = energy/meanBeamMom;  // meanBeamMom is approximately beam energy
00253         switch (id) {
00254            case 22:
00255               _histMultiPhoton->fill(-1.*std::log(scaledMom), weight);
00256               _weightedTotalNumPhoton += weight;
00257               break;
00258            case -321:
00259            case 321:
00260               _weightedTotalNumKPlus += weight;
00261               _histMultiKPlus->fill(scaledMom, weight);
00262               break;
00263            case 211:
00264            case -211:
00265               _histMultiPiPlus->fill(scaledMom, weight);
00266               _weightedTotalNumPiPlus += weight;
00267               break;
00268            case 2212:
00269            case -2212:
00270               _histMultiP->fill(scaledMom, weight);
00271               _weightedTotalNumP += weight;
00272               break;
00273            case 111:
00274               _histMultiPi0->fill(scaledMom, weight);
00275               _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight);
00276               _weightedTotalNumPi0 += weight;
00277               break;
00278            case 221:
00279               if (scaledMom >= 0.1) {
00280                 _histMultiEta->fill(scaledEnergy, weight);
00281                 _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight);
00282                 _weightedTotalNumEta += weight;
00283               }
00284               break;
00285            case 331:
00286               if (scaledMom >= 0.1) {
00287                 _histMultiEtaPrime->fill(scaledEnergy, weight);
00288                 _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight);
00289                 _weightedTotalNumEtaPrime += weight;
00290               }
00291               break;
00292            case 130: //klong
00293            case 310: //kshort
00294               _histMultiK0->fill(scaledMom, weight);
00295               _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight);
00296               _weightedTotalNumK0 += weight;
00297               break;
00298            case 113:
00299               _histMultiRho->fill(scaledMom, weight);
00300               _histMeanMultiRho->fill(_histMeanMultiRho->binMean(0), weight);
00301               _weightedTotalNumRho += weight;
00302               break;
00303            case 223:
00304               _histMultiOmega782->fill(scaledMom, weight);
00305               _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight);
00306               _weightedTotalNumOmega782 += weight;
00307               break;
00308            case 333:
00309               _histMultiPhi->fill(scaledMom, weight);
00310               _histMeanMultiPhi->fill(_histMeanMultiPhi->binMean(0), weight);
00311               _weightedTotalNumPhi += weight;
00312               break;
00313            case 313:
00314            case -313:
00315               _histMultiKStar892_0->fill(scaledMom, weight);
00316               _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight);
00317               _weightedTotalNumKStar892_0 += weight;
00318               break;
00319            case 323:
00320            case -323:
00321               _histMultiKStar892Plus->fill(scaledEnergy, weight);
00322               _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight);
00323               _weightedTotalNumKStar892Plus += weight;
00324               break;
00325            case 3122:
00326            case -3122:
00327               _histMultiLambda0->fill(scaledMom, weight);
00328               _histMeanMultiLambda0->fill(_histMeanMultiLambda0->binMean(0), weight);
00329               _weightedTotalNumLambda0 += weight;
00330               break;
00331            case 3212:
00332            case -3212:
00333               _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight);
00334               break;
00335            case 3312:
00336            case -3312:
00337               _histMultiXiMinus->fill(scaledEnergy, weight);
00338               _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight);
00339               _weightedTotalNumXiMinus += weight;
00340               break;
00341            case 3114:
00342            case -3114:
00343            case 3224:
00344            case -3224:
00345               _histMultiSigma1385Plus->fill(scaledEnergy, weight);
00346               _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight);
00347               _weightedTotalNumSigma1385Plus += weight;
00348               break;
00349            case 3324:
00350            case -3324:
00351               _histMultiXi1530_0->fill(scaledEnergy, weight);
00352               _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight);
00353               _weightedTotalNumXi1530_0 += weight;
00354               break;
00355            case 3334:
00356               _histMeanMultiOmegaOmegaBar->fill(_histMeanMultiOmegaOmegaBar->binMean(0), weight);
00357               break;
00358         }
00359       }
00360 
00361     }

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().

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

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().

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

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 HadronicFinalState::project(), and FinalStateHCM::project().

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

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 125 of file Analysis.hh.

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

00125                                                  {
00126       return info().authors();
00127     }

const PdgIdPair beamIds (  )  const [inherited]

Incoming beam IDs for this run.

Definition at line 78 of file Analysis.cc.

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

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

00078                                           {
00079     return handler().beamIds();
00080   }

const ParticlePair & beams (  )  const [inherited]
virtual std::string bibKey (  )  const [inline, virtual, inherited]

BibTeX citation key for this article.

Definition at line 178 of file Analysis.hh.

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

00178                                      {
00179       return info().bibKey();
00180     }

virtual std::string bibTeX (  )  const [inline, virtual, inherited]

BibTeX citation entry for this article.

Definition at line 183 of file Analysis.hh.

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

00183                                      {
00184       return info().bibTeX();
00185     }

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 223 of file Analysis.cc.

References Analysis::binEdges().

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

const BinEdges & binEdges ( const std::string &  hname  )  const [protected, inherited]

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

Definition at line 208 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(), D0_2008_S7837160::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9002537::init(), and ATLAS_2010_S8894728::init().

00208                                                               {
00209     _cacheBinEdges();
00210     MSG_TRACE("Using histo bin edges for " << name() << ":" << hname);
00211     const BinEdges& edges = _histBinEdges.find(hname)->second;
00212     if (getLog().isActive(Log::TRACE)) {
00213       stringstream edges_ss;
00214       foreach (const double be, edges) {
00215         edges_ss << " " << be;
00216       }
00217       MSG_TRACE("Edges:" << edges_ss.str());
00218     }
00219     return edges;
00220   }

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 418 of file Analysis.cc.

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

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

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 400 of file Analysis.cc.

References Analysis::bookDataPointSet().

00403                                                                                         {
00404     IDataPointSet* dps = bookDataPointSet(hname, title, xtitle, ytitle);
00405     for (size_t pt = 0; pt < npts; ++pt) {
00406       const double binwidth = (upper-lower)/npts;
00407       const double bincentre = lower + (pt + 0.5) * binwidth;
00408       dps->addPoint();
00409       IMeasurement* meas = dps->point(pt)->coordinate(0);
00410       meas->setValue(bincentre);
00411       meas->setErrorPlus(binwidth/2.0);
00412       meas->setErrorMinus(binwidth/2.0);
00413     }
00414     return dps;
00415   }

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 388 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(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2005_S6217184::init(), CDF_1996_S3418421::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9002537::init(), ATLAS_2010_S8894728::init(), and ALEPH_2004_S5765862::init().

00389                                                                                         {
00390     _makeHistoDir();
00391     const string path = histoPath(hname);
00392     IDataPointSet* dps = datapointsetFactory().create(path, title, 2);
00393     MSG_TRACE("Made data point set " << hname <<  " for " << name());
00394     dps->setXTitle(xtitle);
00395     dps->setYTitle(ytitle);
00396     return dps;
00397   }

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 242 of file Analysis.cc.

References Analysis::bookHistogram1D().

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

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 251 of file Analysis.cc.

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

00253   {
00254     // Get the bin edges (only read the AIDA file once)
00255     const BinEdges edges = binEdges(hname);
00256     _makeHistoDir();
00257     const string path = histoPath(hname);
00258     IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, edges);
00259     MSG_TRACE("Made histogram " << hname <<  " for " << name());
00260     hist->setXTitle(xtitle);
00261     hist->setYTitle(ytitle);
00262     return hist;
00263   }

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 280 of file Analysis.cc.

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

00283                                                                                       {
00284     _makeHistoDir();
00285     const string path = histoPath(hname);
00286     IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, binedges);
00287     MSG_TRACE("Made histogram " << hname <<  " for " << name());
00288     hist->setXTitle(xtitle);
00289     hist->setYTitle(ytitle);
00290     return hist;
00291   }

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 266 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_WWJETS::init(), MC_WPOL::init(), MC_WJETS::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_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::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_2002_S4796047::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_S9120807::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S9002537::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_S8971293::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_2010_S8706239::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), and ALEPH_1991_S2435284::init().

00269                                                                                       {
00270     _makeHistoDir();
00271     const string path = histoPath(hname);
00272     IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, nbins, lower, upper);
00273     MSG_TRACE("Made histogram " << hname <<  " for " << name());
00274     hist->setXTitle(xtitle);
00275     hist->setYTitle(ytitle);
00276     return hist;
00277   }

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 313 of file Analysis.cc.

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

00317                                                             {
00318     _makeHistoDir();
00319     const string path = histoPath(hname);
00320     IHistogram2D* hist =
00321       histogramFactory().createHistogram2D(path, title, xbinedges, ybinedges);
00322     MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
00323     hist->setXTitle(xtitle);
00324     hist->setYTitle(ytitle);
00325     hist->setZTitle(ztitle);
00326     return hist;
00327   }

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 294 of file Analysis.cc.

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

00298                                                             {
00299     _makeHistoDir();
00300     const string path = histoPath(hname);
00301     IHistogram2D* hist =
00302       histogramFactory().createHistogram2D(path, title, nxbins, xlower, xupper,
00303                        nybins, ylower, yupper);
00304     MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
00305     hist->setXTitle(xtitle);
00306     hist->setYTitle(ytitle);
00307     hist->setZTitle(ztitle);
00308     return hist;
00309   }

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 333 of file Analysis.cc.

References Analysis::bookProfile1D().

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

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 341 of file Analysis.cc.

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

00343   {
00344     // Get the bin edges (only read the AIDA file once)
00345     const BinEdges edges = binEdges(hname);
00346     _makeHistoDir();
00347     const string path = histoPath(hname);
00348     IProfile1D* prof = histogramFactory().createProfile1D(path, title, edges);
00349     MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
00350     prof->setXTitle(xtitle);
00351     prof->setYTitle(ytitle);
00352     return prof;
00353   }

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 370 of file Analysis.cc.

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

00373                                                                                   {
00374     _makeHistoDir();
00375     const string path = histoPath(hname);
00376     IProfile1D* prof = histogramFactory().createProfile1D(path, title, binedges);
00377     MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
00378     prof->setXTitle(xtitle);
00379     prof->setYTitle(ytitle);
00380     return prof;
00381   }

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 356 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_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_2002_S4796047::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(), and ALICE_2010_S8706239::init().

00359                                                                                   {
00360     _makeHistoDir();
00361     const string path = histoPath(hname);
00362     IProfile1D* prof = histogramFactory().createProfile1D(path, title, nbins, lower, upper);
00363     MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
00364     prof->setXTitle(xtitle);
00365     prof->setYTitle(ytitle);
00366     return prof;
00367   }

virtual std::string collider (  )  const [inline, virtual, inherited]

Collider on which the experiment ran.

Definition at line 163 of file Analysis.hh.

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

00163                                        {
00164       return info().collider();
00165     }

double crossSection (  )  const [protected, inherited]

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

Definition at line 171 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_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_DIPHOTON::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(), 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_S9120807::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), and ATLAS_2010_CONF_2010_049::finalize().

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

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 179 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(), and ATLAS_2010_S8817804::finalize().

00179                                               {
00180     const double sumW = sumOfWeights();
00181     assert(sumW > 0);
00182     return _crossSection / sumW;
00183   }

IDataPointSetFactory & datapointsetFactory (  )  [protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 65 of file Analysis.cc.

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

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

00065                                                       {
00066     return handler().datapointsetFactory();
00067   }

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 144 of file Analysis.hh.

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

00144                                           {
00145       return info().description();
00146     }

virtual std::string experiment (  )  const [inline, virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 158 of file Analysis.hh.

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

00158                                          {
00159       return info().experiment();
00160     }

void finalize (  )  [inline, virtual]

Finalize.

Implements Analysis.

Definition at line 366 of file ALEPH_1996_S3486095.cc.

References ALEPH_1996_S3486095::_hist1MinusT, ALEPH_1996_S3486095::_histAplanarity, ALEPH_1996_S3486095::_histChMult, ALEPH_1996_S3486095::_histCParam, ALEPH_1996_S3486095::_histHeavyJetMass, ALEPH_1996_S3486095::_histLogScaledMom, ALEPH_1996_S3486095::_histMeanChMult, ALEPH_1996_S3486095::_histMeanChMultRapt05, ALEPH_1996_S3486095::_histMeanChMultRapt10, ALEPH_1996_S3486095::_histMeanChMultRapt15, ALEPH_1996_S3486095::_histMeanChMultRapt20, ALEPH_1996_S3486095::_histMeanMultiEta, ALEPH_1996_S3486095::_histMeanMultiEtaPrime, ALEPH_1996_S3486095::_histMeanMultiK0, ALEPH_1996_S3486095::_histMeanMultiKStar892_0, ALEPH_1996_S3486095::_histMeanMultiKStar892Plus, ALEPH_1996_S3486095::_histMeanMultiLambda0, ALEPH_1996_S3486095::_histMeanMultiOmega782, ALEPH_1996_S3486095::_histMeanMultiOmegaOmegaBar, ALEPH_1996_S3486095::_histMeanMultiPhi, ALEPH_1996_S3486095::_histMeanMultiPi0, ALEPH_1996_S3486095::_histMeanMultiRho, ALEPH_1996_S3486095::_histMeanMultiSigma0, ALEPH_1996_S3486095::_histMeanMultiSigma1385Plus, ALEPH_1996_S3486095::_histMeanMultiXi1530_0, ALEPH_1996_S3486095::_histMeanMultiXiMinus, ALEPH_1996_S3486095::_histMultiEta, ALEPH_1996_S3486095::_histMultiEtaPrime, ALEPH_1996_S3486095::_histMultiK0, ALEPH_1996_S3486095::_histMultiKPlus, ALEPH_1996_S3486095::_histMultiKStar892_0, ALEPH_1996_S3486095::_histMultiKStar892Plus, ALEPH_1996_S3486095::_histMultiLambda0, ALEPH_1996_S3486095::_histMultiOmega782, ALEPH_1996_S3486095::_histMultiP, ALEPH_1996_S3486095::_histMultiPhi, ALEPH_1996_S3486095::_histMultiPhoton, ALEPH_1996_S3486095::_histMultiPi0, ALEPH_1996_S3486095::_histMultiPiPlus, ALEPH_1996_S3486095::_histMultiRho, ALEPH_1996_S3486095::_histMultiSigma1385Plus, ALEPH_1996_S3486095::_histMultiXi1530_0, ALEPH_1996_S3486095::_histMultiXiMinus, ALEPH_1996_S3486095::_histOblateness, ALEPH_1996_S3486095::_histPtSIn, ALEPH_1996_S3486095::_histPtSOut, ALEPH_1996_S3486095::_histRapidityT, ALEPH_1996_S3486095::_histScaledMom, ALEPH_1996_S3486095::_histSphericity, ALEPH_1996_S3486095::_histTMinor, ALEPH_1996_S3486095::_histY3, ALEPH_1996_S3486095::_weightedTotalPartNum, Analysis::normalize(), Analysis::scale(), and Analysis::sumOfWeights().

00366                     {
00367       // Normalize inclusive single particle distributions to the average number
00368       // of charged particles per event.
00369       const double avgNumParts = _weightedTotalPartNum / sumOfWeights();
00370 
00371       normalize(_histPtSIn, avgNumParts);
00372       normalize(_histPtSOut, avgNumParts);
00373 
00374       normalize(_histRapidityT, avgNumParts);
00375       normalize(_histY3);
00376 
00377       normalize(_histLogScaledMom, avgNumParts);
00378       normalize(_histScaledMom, avgNumParts);
00379 
00380       // particle spectra
00381       scale(_histMultiPiPlus        ,1./sumOfWeights());
00382       scale(_histMultiKPlus         ,1./sumOfWeights());
00383       scale(_histMultiP             ,1./sumOfWeights());
00384       scale(_histMultiPhoton        ,1./sumOfWeights());
00385       scale(_histMultiPi0           ,1./sumOfWeights());
00386       scale(_histMultiEta           ,1./sumOfWeights());
00387       scale(_histMultiEtaPrime      ,1./sumOfWeights());
00388       scale(_histMultiK0            ,1./sumOfWeights());
00389       scale(_histMultiLambda0       ,1./sumOfWeights());
00390       scale(_histMultiXiMinus       ,1./sumOfWeights());
00391       scale(_histMultiSigma1385Plus ,1./sumOfWeights());
00392       scale(_histMultiXi1530_0      ,1./sumOfWeights());
00393       scale(_histMultiRho           ,1./sumOfWeights());
00394       scale(_histMultiOmega782      ,1./sumOfWeights());
00395       scale(_histMultiKStar892_0    ,1./sumOfWeights());
00396       scale(_histMultiPhi           ,1./sumOfWeights());
00397 
00398       scale(_histMultiKStar892Plus  ,1./sumOfWeights());
00399 
00400       //normalize(_histMultiPiPlus        ,_weightedTotalNumPiPlus / sumOfWeights());
00401       //normalize(_histMultiKPlus         ,_weightedTotalNumKPlus/sumOfWeights());
00402       //normalize(_histMultiP             ,_weightedTotalNumP/sumOfWeights());
00403       //normalize(_histMultiPhoton            ,_weightedTotalNumPhoton/sumOfWeights());
00404       //normalize(_histMultiPi0           ,_weightedTotalNumPi0/sumOfWeights());
00405       //normalize(_histMultiEta           ,_weightedTotalNumEta/sumOfWeights());
00406       //normalize(_histMultiEtaPrime      ,_weightedTotalNumEtaPrime/sumOfWeights());
00407       //normalize(_histMultiK0            ,_weightedTotalNumK0/sumOfWeights());
00408       //normalize(_histMultiLambda0       ,_weightedTotalNumLambda0/sumOfWeights());
00409       //normalize(_histMultiXiMinus       ,_weightedTotalNumXiMinus/sumOfWeights());
00410       //normalize(_histMultiSigma1385Plus ,_weightedTotalNumSigma1385Plus/sumOfWeights());
00411       //normalize(_histMultiXi1530_0      ,_weightedTotalNumXi1530_0 /sumOfWeights());
00412       //normalize(_histMultiRho           ,_weightedTotalNumRho/sumOfWeights());
00413       //normalize(_histMultiOmegaMinus    ,_weightedTotalNumOmegaMinus/sumOfWeights());
00414       //normalize(_histMultiKStar892_0    ,_weightedTotalNumKStar892_0/sumOfWeights());
00415       //normalize(_histMultiPhi           ,_weightedTotalNumPhi/sumOfWeights());
00416 
00417       //normalize(_histMultiKStar892Plus  ,_weightedTotalNumKStar892Plus/sumOfWeights());
00418 
00419       // event shape
00420       normalize(_hist1MinusT);
00421       normalize(_histTMinor);
00422       normalize(_histOblateness);
00423 
00424       normalize(_histSphericity);
00425       normalize(_histAplanarity);
00426       normalize(_histHeavyJetMass);
00427       normalize(_histCParam);
00428 
00429 
00430       // mean multiplicities
00431       scale(_histChMult              , 2.0/sumOfWeights()); // taking into account the binwidth of 2
00432       scale(_histMeanChMult          , 1.0/sumOfWeights());
00433       scale(_histMeanChMultRapt05    , 1.0/sumOfWeights());
00434       scale(_histMeanChMultRapt10    , 1.0/sumOfWeights());
00435       scale(_histMeanChMultRapt15    , 1.0/sumOfWeights());
00436       scale(_histMeanChMultRapt20    , 1.0/sumOfWeights());
00437 
00438 
00439       scale(_histMeanMultiPi0          , 1.0/sumOfWeights());
00440       scale(_histMeanMultiEta          , 1.0/sumOfWeights());
00441       scale(_histMeanMultiEtaPrime     , 1.0/sumOfWeights());
00442       scale(_histMeanMultiK0           , 1.0/sumOfWeights());
00443       scale(_histMeanMultiRho          , 1.0/sumOfWeights());
00444       scale(_histMeanMultiOmega782     , 1.0/sumOfWeights());
00445       scale(_histMeanMultiPhi          , 1.0/sumOfWeights());
00446       scale(_histMeanMultiKStar892Plus , 1.0/sumOfWeights());
00447       scale(_histMeanMultiKStar892_0   , 1.0/sumOfWeights());
00448       scale(_histMeanMultiLambda0      , 1.0/sumOfWeights());
00449       scale(_histMeanMultiSigma0       , 1.0/sumOfWeights());
00450       scale(_histMeanMultiXiMinus      , 1.0/sumOfWeights());
00451       scale(_histMeanMultiSigma1385Plus, 1.0/sumOfWeights());
00452       scale(_histMeanMultiXi1530_0     , 1.0/sumOfWeights());
00453       scale(_histMeanMultiOmegaOmegaBar, 1.0/sumOfWeights());
00454     }

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 102 of file Analysis.cc.

References Analysis::name().

Referenced by CDF_2004_S5839831::_calcTransCones(), CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), D0_1996_S3214044::_safeMass(), CDF_1997_S3541940::_safeMass(), CDF_1996_S3349578::_safeMass(), CDF_1996_S3349578::_threeJetAnalysis(), ZEUS_2001_S4815815::analyze(), UA5_1986_S1583476::analyze(), TASSO_1990_S2148048::analyze(), STAR_2009_UE_HELEN::analyze(), STAR_2008_S7993412::analyze(), STAR_2006_S6870392::analyze(), STAR_2006_S6860818::analyze(), STAR_2006_S6500200::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), OPAL_1998_S3780481::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_LEADINGJETS::analyze(), MC_JetAnalysis::analyze(), JADE_1998_S3612880::analyze(), H1_1994_S2919893::analyze(), DELPHI_2003_WUD_03_11::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1996_S3430090::analyze(), DELPHI_1995_S3137023::analyze(), D0_2009_S8349509::analyze(), D0_2009_S8202443::analyze(), D0_2008_S7863608::analyze(), D0_2008_S7837160::analyze(), D0_2008_S7719523::analyze(), D0_2008_S7662670::analyze(), D0_2008_S7554427::analyze(), D0_2007_S7075677::analyze(), D0_2004_S5992206::analyze(), CMS_2011_S8884919::analyze(), CDF_2010_S8591881_QCD::analyze(), CDF_2010_S8591881_DY::analyze(), CDF_2009_S8383952::analyze(), CDF_2008_S8095620::analyze(), CDF_2008_S7782535::analyze(), CDF_2008_S7540469::analyze(), CDF_2008_NOTE_9351::analyze(), CDF_2008_LEADINGJETS::analyze(), CDF_2006_S6653332::analyze(), CDF_2004_S5839831::analyze(), CDF_2001_S4751469::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_1991_S2435284::analyze(), Analysis::binEdges(), CDF_2009_S8057893::CDF_2009_S8057893::analyze(), UA5_1986_S1583476::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), DELPHI_2003_WUD_03_11::finalize(), D0_2001_S4674421::finalize(), CDF_2006_S6653332::finalize(), ZEUS_2001_S4815815::init(), and ALEPH_2004_S5765862::init().

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

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().

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

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().

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

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().

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

ProjectionHandler& getProjHandler (  )  const [inline, protected, inherited]
AnalysisHandler& handler (  )  const [inline, inherited]
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 83 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(), and Analysis::histoPath().

00083                                         {
00084     /// @todo This doesn't change: calc and cache at first use!
00085     string path = "/" + name();
00086     if (handler().runName().length() > 0) {
00087       path = "/" + handler().runName() + path;
00088     }
00089     while (find_first(path, "//")) {
00090       replace_all(path, "//", "/");
00091     }
00092     return path;
00093   }

IHistogramFactory & histogramFactory (  )  [protected, inherited]
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 96 of file Analysis.cc.

References Analysis::histoDir().

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

00096                                                             {
00097     const string path = histoDir() + "/" + hname;
00098     return path;
00099   }

AnalysisInfo& info (  )  [inline, inherited]

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

Definition at line 237 of file Analysis.hh.

References Analysis::_info.

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

const AnalysisInfo& info (  )  const [inline, inherited]
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.

Implements Analysis.

Definition at line 54 of file ALEPH_1996_S3486095.cc.

References ALEPH_1996_S3486095::_hist1MinusT, ALEPH_1996_S3486095::_histAplanarity, ALEPH_1996_S3486095::_histChMult, ALEPH_1996_S3486095::_histCParam, ALEPH_1996_S3486095::_histHeavyJetMass, ALEPH_1996_S3486095::_histLogScaledMom, ALEPH_1996_S3486095::_histMeanChMult, ALEPH_1996_S3486095::_histMeanChMultRapt05, ALEPH_1996_S3486095::_histMeanChMultRapt10, ALEPH_1996_S3486095::_histMeanChMultRapt15, ALEPH_1996_S3486095::_histMeanChMultRapt20, ALEPH_1996_S3486095::_histMeanMultiEta, ALEPH_1996_S3486095::_histMeanMultiEtaPrime, ALEPH_1996_S3486095::_histMeanMultiK0, ALEPH_1996_S3486095::_histMeanMultiKStar892_0, ALEPH_1996_S3486095::_histMeanMultiKStar892Plus, ALEPH_1996_S3486095::_histMeanMultiLambda0, ALEPH_1996_S3486095::_histMeanMultiOmega782, ALEPH_1996_S3486095::_histMeanMultiOmegaOmegaBar, ALEPH_1996_S3486095::_histMeanMultiPhi, ALEPH_1996_S3486095::_histMeanMultiPi0, ALEPH_1996_S3486095::_histMeanMultiRho, ALEPH_1996_S3486095::_histMeanMultiSigma0, ALEPH_1996_S3486095::_histMeanMultiSigma1385Plus, ALEPH_1996_S3486095::_histMeanMultiXi1530_0, ALEPH_1996_S3486095::_histMeanMultiXiMinus, ALEPH_1996_S3486095::_histMultiEta, ALEPH_1996_S3486095::_histMultiEtaPrime, ALEPH_1996_S3486095::_histMultiK0, ALEPH_1996_S3486095::_histMultiKPlus, ALEPH_1996_S3486095::_histMultiKStar892_0, ALEPH_1996_S3486095::_histMultiKStar892Plus, ALEPH_1996_S3486095::_histMultiLambda0, ALEPH_1996_S3486095::_histMultiOmega782, ALEPH_1996_S3486095::_histMultiP, ALEPH_1996_S3486095::_histMultiPhi, ALEPH_1996_S3486095::_histMultiPhoton, ALEPH_1996_S3486095::_histMultiPi0, ALEPH_1996_S3486095::_histMultiPiPlus, ALEPH_1996_S3486095::_histMultiRho, ALEPH_1996_S3486095::_histMultiSigma1385Plus, ALEPH_1996_S3486095::_histMultiXi1530_0, ALEPH_1996_S3486095::_histMultiXiMinus, ALEPH_1996_S3486095::_histOblateness, ALEPH_1996_S3486095::_histPtSIn, ALEPH_1996_S3486095::_histPtSOut, ALEPH_1996_S3486095::_histRapidityT, ALEPH_1996_S3486095::_histScaledMom, ALEPH_1996_S3486095::_histSphericity, ALEPH_1996_S3486095::_histTMinor, ALEPH_1996_S3486095::_histY3, ProjectionApplier::addProjection(), Analysis::bookHistogram1D(), and FastJets::DURHAM.

00054                 {
00055       // Set up projections
00056       addProjection(Beam(), "Beams");
00057       const ChargedFinalState cfs;
00058       addProjection(cfs, "FS");
00059       addProjection(UnstableFinalState(), "UFS");
00060       addProjection(FastJets(cfs, FastJets::DURHAM, 0.7), "DurhamJets");
00061       addProjection(Sphericity(cfs), "Sphericity");
00062       addProjection(ParisiTensor(cfs), "Parisi");
00063       const Thrust thrust(cfs);
00064       addProjection(thrust, "Thrust");
00065       addProjection(Hemispheres(thrust), "Hemispheres");
00066 
00067       // Book histograms
00068       _histSphericity   = bookHistogram1D(1, 1, 1);
00069       _histAplanarity   = bookHistogram1D(2, 1, 1);
00070 
00071       _hist1MinusT      = bookHistogram1D(3, 1, 1);
00072       _histTMinor       = bookHistogram1D(4, 1, 1);
00073 
00074       _histY3           = bookHistogram1D(5, 1, 1);
00075       _histHeavyJetMass = bookHistogram1D(6, 1, 1);
00076       _histCParam       = bookHistogram1D(7, 1, 1);
00077       _histOblateness   = bookHistogram1D(8, 1, 1);
00078 
00079       _histScaledMom    = bookHistogram1D(9, 1, 1);
00080       _histRapidityT    = bookHistogram1D(10, 1, 1);
00081 
00082       _histPtSIn        = bookHistogram1D(11, 1, 1);
00083       _histPtSOut       = bookHistogram1D(12, 1, 1);
00084 
00085       _histLogScaledMom = bookHistogram1D(17, 1, 1);
00086 
00087       _histChMult       = bookHistogram1D(18, 1, 1);
00088       _histMeanChMult   = bookHistogram1D(19, 1, 1);
00089 
00090       _histMeanChMultRapt05= bookHistogram1D(20, 1, 1);
00091       _histMeanChMultRapt10= bookHistogram1D(21, 1, 1);
00092       _histMeanChMultRapt15= bookHistogram1D(22, 1, 1);
00093       _histMeanChMultRapt20= bookHistogram1D(23, 1, 1);
00094 
00095 
00096       // Particle spectra
00097       _histMultiPiPlus        = bookHistogram1D(25, 1, 1);
00098       _histMultiKPlus         = bookHistogram1D(26, 1, 1);
00099       _histMultiP             = bookHistogram1D(27, 1, 1);
00100       _histMultiPhoton        = bookHistogram1D(28, 1, 1);
00101       _histMultiPi0           = bookHistogram1D(29, 1, 1);
00102       _histMultiEta           = bookHistogram1D(30, 1, 1);
00103       _histMultiEtaPrime      = bookHistogram1D(31, 1, 1);
00104       _histMultiK0            = bookHistogram1D(32, 1, 1);
00105       _histMultiLambda0       = bookHistogram1D(33, 1, 1);
00106       _histMultiXiMinus       = bookHistogram1D(34, 1, 1);
00107       _histMultiSigma1385Plus = bookHistogram1D(35, 1, 1);
00108       _histMultiXi1530_0      = bookHistogram1D(36, 1, 1);
00109       _histMultiRho           = bookHistogram1D(37, 1, 1);
00110       _histMultiOmega782      = bookHistogram1D(38, 1, 1);
00111       _histMultiKStar892_0    = bookHistogram1D(39, 1, 1);
00112       _histMultiPhi           = bookHistogram1D(40, 1, 1);
00113 
00114       _histMultiKStar892Plus  = bookHistogram1D(43, 1, 1);
00115 
00116       // Mean multiplicities
00117       _histMeanMultiPi0           = bookHistogram1D(44, 1,  2);
00118       _histMeanMultiEta           = bookHistogram1D(44, 1,  3);
00119       _histMeanMultiEtaPrime      = bookHistogram1D(44, 1,  4);
00120       _histMeanMultiK0            = bookHistogram1D(44, 1,  5);
00121       _histMeanMultiRho           = bookHistogram1D(44, 1,  6);
00122       _histMeanMultiOmega782      = bookHistogram1D(44, 1,  7);
00123       _histMeanMultiPhi           = bookHistogram1D(44, 1,  8);
00124       _histMeanMultiKStar892Plus  = bookHistogram1D(44, 1,  9);
00125       _histMeanMultiKStar892_0    = bookHistogram1D(44, 1, 10);
00126       _histMeanMultiLambda0       = bookHistogram1D(44, 1, 11);
00127       _histMeanMultiSigma0        = bookHistogram1D(44, 1, 12);
00128       _histMeanMultiXiMinus       = bookHistogram1D(44, 1, 13);
00129       _histMeanMultiSigma1385Plus = bookHistogram1D(44, 1, 14);
00130       _histMeanMultiXi1530_0      = bookHistogram1D(44, 1, 15);
00131       _histMeanMultiOmegaOmegaBar = bookHistogram1D(44, 1, 16);
00132     }

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 127 of file Analysis.cc.

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

00127                                                                                   {
00128     PdgIdPair beams(beam1, beam2);
00129     pair<double,double> energies(e1, e2);
00130     return isCompatible(beams, energies);
00131   }

bool isCompatible ( const ParticlePair beams  )  const [inherited]

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

Definition at line 121 of file Analysis.cc.

Referenced by Analysis::isCompatible().

00121                                                              {
00122     return isCompatible(beams.first.pdgId(),  beams.second.pdgId(),
00123                         beams.first.energy(), beams.second.energy());
00124   }

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

Get bin edges with logarithmic widths.

Definition at line 229 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().

00229                                                                          {
00230     assert(lower>0.0);
00231     assert(upper>lower);
00232     double loglower=log10(lower);
00233     double logupper=log10(upper);
00234     vector<double> binedges;
00235     double stepwidth=(logupper-loglower)/double(nbins);
00236     for (size_t i=0; i<=nbins; ++i) {
00237       binedges.push_back(pow(10.0, loglower+double(i)*stepwidth));
00238     }
00239     return binedges;
00240   }

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 112 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().

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

bool needsCrossSection (  )  const [inline, inherited]

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

Definition at line 221 of file Analysis.hh.

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

00221                                    {
00222       return info().needsCrossSection();
00223     }

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 527 of file Analysis.cc.

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

00527                                                                 {
00528     if (!histo) {
00529       MSG_ERROR("Failed to normalize histo=NULL in analysis "
00530                 << name() << " (norm=" << norm << ")");
00531       return;
00532     }
00533     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00534     MSG_TRACE("Normalizing histo " << hpath << " to " << norm);
00535 
00536     double oldintg = 0.0;
00537     int nxBins = histo->xAxis().bins();
00538     int nyBins = histo->yAxis().bins();
00539     for (int ixBin = 0; ixBin != nxBins; ++ixBin)
00540       for (int iyBin = 0; iyBin != nyBins; ++iyBin) {
00541       // Leaving out factor of binWidth because AIDA's "height"
00542       // already includes a width factor.
00543     oldintg += histo->binHeight(ixBin, iyBin); // * histo->axis().binWidth(iBin);
00544     }
00545     if (oldintg == 0.0) {
00546       MSG_WARNING("Histo " << hpath << " has null integral during normalization");
00547       return;
00548     }
00549 
00550     // Scale by the normalisation factor.
00551     scale(histo, norm/oldintg);
00552   }

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 460 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_S8968497::finalize(), CMS_2011_S8957746::finalize(), CMS_2011_S8884919::finalize(), CDF_2009_NOTE_9936::finalize(), CDF_2004_S5839831::finalize(), CDF_2002_S4796047::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(), and ALEPH_1996_S3486095::finalize().

00460                                                                 {
00461     if (!histo) {
00462       MSG_ERROR("Failed to normalize histo=NULL in analysis "
00463                 << name() << " (norm=" << norm << ")");
00464       return;
00465     }
00466     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00467     MSG_TRACE("Normalizing histo " << hpath << " to " << norm);
00468 
00469     double oldintg = 0.0;
00470     int nBins = histo->axis().bins();
00471     for (int iBin = 0; iBin != nBins; ++iBin) {
00472       // Leaving out factor of binWidth because AIDA's "height" already includes a width factor.
00473       oldintg += histo->binHeight(iBin); // * histo->axis().binWidth(iBin);
00474     }
00475     if (oldintg == 0.0) {
00476       MSG_WARNING("Histo " << hpath << " has null integral during normalization");
00477       return;
00478     }
00479 
00480     // Scale by the normalisation factor.
00481     scale(histo, norm/oldintg);
00482   }

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 108 of file Analysis.cc.

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

00108                                    {
00109     return handler().numEvents();
00110   }

virtual std::vector<std::string> references (  )  const [inline, virtual, inherited]

Journal, and preprint references.

Definition at line 173 of file Analysis.hh.

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

00173                                                     {
00174       return info().references();
00175     }

virtual const std::vector<PdgIdPair>& requiredBeams (  )  const [inline, virtual, inherited]

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

Definition at line 199 of file Analysis.hh.

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

00199                                                               {
00200       return info().beams();
00201     }

virtual const std::vector<std::pair<double, double> >& requiredEnergies (  )  const [inline, virtual, inherited]

Sets of valid beam energy pairs, in GeV.

Definition at line 210 of file Analysis.hh.

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

Referenced by Analysis::setRequiredEnergies().

00210                                                                                 {
00211       return info().energies();
00212     }

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 153 of file Analysis.hh.

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

00153                                       {
00154       return info().runInfo();
00155     }

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 555 of file Analysis.cc.

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

00555                                                              {
00556     if (!histo) {
00557       MSG_ERROR("Failed to scale histo=NULL in analysis "
00558                 << name() << " (scale=" << scale << ")");
00559       return;
00560     }
00561     const string hpath =
00562       tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00563     MSG_TRACE("Scaling histo " << hpath);
00564 
00565     vector<double> x, y, z, ex, ey, ez;
00566     for (size_t ix = 0, Nx = histo->xAxis().bins(); ix < Nx; ++ix)
00567       for (size_t iy = 0, Ny = histo->yAxis().bins(); iy < Ny; ++iy) {
00568     x.push_back(0.5 * (histo->xAxis().binLowerEdge(ix) +
00569                histo->xAxis().binUpperEdge(ix)));
00570     ex.push_back(histo->xAxis().binWidth(ix)*0.5);
00571     y.push_back(0.5 * (histo->yAxis().binLowerEdge(iy) +
00572                histo->yAxis().binUpperEdge(iy)));
00573     ey.push_back(histo->yAxis().binWidth(iy)*0.5);
00574 
00575     // "Bin height" is a misnomer in the AIDA spec: width is neglected.
00576     // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
00577     z.push_back(histo->binHeight(ix, iy)*scale/
00578             (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy)));
00579     // "Bin error" is a misnomer in the AIDA spec: width is neglected.
00580     // We'd like to do this: ey.push_back(histo->binError(i) * scale);
00581     ez.push_back(histo->binError(ix, iy)*scale/
00582              (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy)));
00583     }
00584 
00585     string title = histo->title();
00586     string xtitle = histo->xtitle();
00587     string ytitle = histo->ytitle();
00588     string ztitle = histo->ztitle();
00589 
00590     tree().mkdir("/tmpnormalize");
00591     tree().mv(hpath, "/tmpnormalize");
00592 
00593     AIDA::IDataPointSet* dps =
00594       datapointsetFactory().createXYZ(hpath, title, x, y, z, ex, ey, ez);
00595     dps->setXTitle(xtitle);
00596     dps->setYTitle(ytitle);
00597     dps->setZTitle(ztitle);
00598 
00599     tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
00600     tree().rmdir("/tmpnormalize");
00601 
00602     // Set histo pointer to null - it can no longer be used.
00603     histo = 0;
00604   }

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 485 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_WWJETS::finalize(), MC_WPOL::finalize(), MC_WJETS::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_S9120807::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_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(), and BinnedHistogram< T >::scale().

00485                                                              {
00486     if (!histo) {
00487       MSG_ERROR("Failed to scale histo=NULL in analysis "
00488                 << name() << " (scale=" << scale << ")");
00489       return;
00490     }
00491     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00492     MSG_TRACE("Scaling histo " << hpath);
00493 
00494     vector<double> x, y, ex, ey;
00495     for (size_t i = 0, N = histo->axis().bins(); i < N; ++i) {
00496       x.push_back(0.5 * (histo->axis().binLowerEdge(i) + histo->axis().binUpperEdge(i)));
00497       ex.push_back(histo->axis().binWidth(i)*0.5);
00498 
00499       // "Bin height" is a misnomer in the AIDA spec: width is neglected.
00500       // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
00501       y.push_back(histo->binHeight(i)*scale/histo->axis().binWidth(i));
00502 
00503       // "Bin error" is a misnomer in the AIDA spec: width is neglected.
00504       // We'd like to do this: ey.push_back(histo->binError(i) * scale);
00505       ey.push_back(histo->binError(i)*scale/histo->axis().binWidth(i));
00506     }
00507 
00508     string title = histo->title();
00509     string xtitle = histo->xtitle();
00510     string ytitle = histo->ytitle();
00511 
00512     tree().mkdir("/tmpnormalize");
00513     tree().mv(hpath, "/tmpnormalize");
00514 
00515     AIDA::IDataPointSet* dps = datapointsetFactory().createXY(hpath, title, x, y, ex, ey);
00516     dps->setXTitle(xtitle);
00517     dps->setYTitle(ytitle);
00518 
00519     tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
00520     tree().rmdir("/tmpnormalize");
00521 
00522     // Set histo pointer to null - it can no longer be used.
00523     histo = 0;
00524   }

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 244 of file Analysis.hh.

References Analysis::setRequiredBeams().

Referenced by ALEPH_1991_S2435284::ALEPH_1991_S2435284(), ALEPH_1996_S3196992::ALEPH_1996_S3196992(), ALEPH_1996_S3486095::ALEPH_1996_S3486095(), ALEPH_2004_S5765862::ALEPH_2004_S5765862(), ALICE_2010_S8624100::ALICE_2010_S8624100(), ALICE_2010_S8625980::ALICE_2010_S8625980(), ALICE_2010_S8706239::ALICE_2010_S8706239(), ATLAS_2010_S8817804::ATLAS_2010_S8817804(), ATLAS_2011_S8924791::ATLAS_2011_S8924791(), ATLAS_2011_S9120807::ATLAS_2011_S9120807(), BELLE_2006_S6265367::BELLE_2006_S6265367(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1990_S2089246::CDF_1990_S2089246(), CDF_1993_S2742446::CDF_1993_S2742446(), CDF_1994_S2952106::CDF_1994_S2952106(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1996_S3349578::CDF_1996_S3349578(), CDF_1996_S3418421::CDF_1996_S3418421(), CDF_1997_S3541940::CDF_1997_S3541940(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2001_S4751469::CDF_2001_S4751469(), CDF_2002_S4796047::CDF_2002_S4796047(), CDF_2004_S5839831::CDF_2004_S5839831(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2005_S6217184::CDF_2005_S6217184(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_LEADINGJETS::CDF_2008_LEADINGJETS(), CDF_2008_NOTE_9351::CDF_2008_NOTE_9351(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7782535::CDF_2008_S7782535(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2008_S8095620::CDF_2008_S8095620(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8057893::CDF_2009_S8057893::CDF_2009_S8057893(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), CDF_2010_S8591881_DY::CDF_2010_S8591881_DY(), CDF_2010_S8591881_QCD::CDF_2010_S8591881_QCD(), CMS_2011_S8957746::CMS_2011_S8957746(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_2001_S4674421::D0_2001_S4674421(), D0_2004_S5992206::D0_2004_S5992206(), D0_2006_S6438750::D0_2006_S6438750(), D0_2007_S7075677::D0_2007_S7075677(), D0_2008_S6879055::D0_2008_S6879055(), D0_2008_S7554427::D0_2008_S7554427(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7837160::D0_2008_S7837160(), D0_2008_S7863608::D0_2008_S7863608(), D0_2009_S8202443::D0_2009_S8202443(), D0_2009_S8320160::D0_2009_S8320160(), D0_2009_S8349509::D0_2009_S8349509(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), D0_2010_S8671338::D0_2010_S8671338(), DELPHI_1995_S3137023::DELPHI_1995_S3137023(), DELPHI_1996_S3430090::DELPHI_1996_S3430090(), DELPHI_2002_069_CONF_603::DELPHI_2002_069_CONF_603(), E735_1998_S3905616::E735_1998_S3905616(), H1_1994_S2919893::H1_1994_S2919893(), H1_1995_S3167097::H1_1995_S3167097(), H1_2000_S4129130::H1_2000_S4129130(), JADE_1998_S3612880::JADE_1998_S3612880(), JADE_OPAL_2000_S4300807::JADE_OPAL_2000_S4300807(), LHCB_2010_S8758301::LHCB_2010_S8758301(), OPAL_1993_S2692198::OPAL_1993_S2692198(), OPAL_1998_S3780481::OPAL_1998_S3780481(), OPAL_2001_S4553896::OPAL_2001_S4553896(), PDG_HADRON_MULTIPLICITIES::PDG_HADRON_MULTIPLICITIES(), PDG_HADRON_MULTIPLICITIES_RATIOS::PDG_HADRON_MULTIPLICITIES_RATIOS(), SFM_1984_S1178091::SFM_1984_S1178091(), STAR_2006_S6500200::STAR_2006_S6500200(), STAR_2006_S6860818::STAR_2006_S6860818(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), STAR_2008_S7993412::STAR_2008_S7993412(), STAR_2009_UE_HELEN::STAR_2009_UE_HELEN(), TASSO_1990_S2148048::TASSO_1990_S2148048(), UA1_1990_S2044935::UA1_1990_S2044935(), UA5_1986_S1583476::UA5_1986_S1583476(), UA5_1987_S1640666::UA5_1987_S1640666(), UA5_1988_S1867512::UA5_1988_S1867512(), UA5_1989_S1926373::UA5_1989_S1926373(), and ZEUS_2001_S4815815::ZEUS_2001_S4815815().

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

Analysis & setCrossSection ( double  xs  )  [inherited]

Set the cross section from the generator.

Definition at line 165 of file Analysis.cc.

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

00165                                                {
00166     _crossSection = xs;
00167     _gotCrossSection = true;
00168     return *this;
00169   }

Analysis& setNeedsCrossSection ( bool  needed = true  )  [inline, inherited]

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

Definition at line 225 of file Analysis.hh.

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

Referenced by ALICE_2010_S8624100::ALICE_2010_S8624100(), ALICE_2010_S8625980::ALICE_2010_S8625980(), ALICE_2010_S8706239::ALICE_2010_S8706239(), ATLAS_2010_CONF_2010_049::ATLAS_2010_CONF_2010_049(), ATLAS_2010_S8591806::ATLAS_2010_S8591806(), ATLAS_2010_S8817804::ATLAS_2010_S8817804(), ATLAS_2010_S8914702::ATLAS_2010_S8914702(), ATLAS_2010_S8919674::ATLAS_2010_S8919674(), ATLAS_2011_CONF_2011_090::ATLAS_2011_CONF_2011_090(), ATLAS_2011_S8971293::ATLAS_2011_S8971293(), ATLAS_2011_S8983313::ATLAS_2011_S8983313(), ATLAS_2011_S9019561::ATLAS_2011_S9019561(), ATLAS_2011_S9120807::ATLAS_2011_S9120807(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1993_S2742446::CDF_1993_S2742446(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), CMS_2011_S8884919::CMS_2011_S8884919(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_2000_S4480767::D0_2000_S4480767(), D0_2001_S4674421::D0_2001_S4674421(), D0_2006_S6438750::D0_2006_S6438750(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7863608::D0_2008_S7863608(), D0_2009_S8349509::D0_2009_S8349509(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), D0_2010_S8671338::D0_2010_S8671338(), D0_2010_S8821313::D0_2010_S8821313(), LHCB_2010_S8758301::LHCB_2010_S8758301(), MC_DIPHOTON::MC_DIPHOTON(), MC_HJETS::MC_HJETS(), MC_JetAnalysis::MC_JetAnalysis(), MC_JETS::MC_JETS(), MC_PHOTONJETS::MC_PHOTONJETS(), MC_WJETS::MC_WJETS(), MC_WPOL::MC_WPOL(), MC_WWJETS::MC_WWJETS(), MC_XS::MC_XS(), MC_ZJETS::MC_ZJETS(), MC_ZZJETS::MC_ZZJETS(), OPAL_2001_S4553896::OPAL_2001_S4553896(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), and UA1_1990_S2044935::UA1_1990_S2044935().

00225                                                      {
00226       info().setNeedsCrossSection(needed);
00227       return *this;
00228     }

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 203 of file Analysis.hh.

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

Referenced by Analysis::setBeams().

00203                                                                                   {
00204       info().setBeams(requiredBeams);
00205       return *this;
00206     }

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 214 of file Analysis.hh.

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

00214                                                                                                        {
00215       info().setEnergies(requiredEnergies);
00216       return *this;
00217     }

virtual std::string spiresId (  )  const [inline, virtual, inherited]

Get a the SPIRES/Inspire ID code for this analysis.

Definition at line 117 of file Analysis.hh.

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

00117                                        {
00118       return info().spiresId();
00119     }

double sqrtS (  )  const [inherited]

Centre of mass energy for this run.

Definition at line 70 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(), JADE_1998_S3612880::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(), CDF_2002_S4796047::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_2002_S4796047::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().

00070                                {
00071     return handler().sqrtS();
00072   }

virtual std::string status (  )  const [inline, virtual, inherited]

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

Definition at line 188 of file Analysis.hh.

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

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

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 134 of file Analysis.hh.

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

00134                                       {
00135       return info().summary();
00136     }

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 113 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_WWJETS::finalize(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::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_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_S9120807::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), and ALEPH_1991_S2435284::finalize().

00113                                       {
00114     return handler().sumOfWeights();
00115   }

virtual std::vector<std::string> todos (  )  const [inline, virtual, inherited]

Any work to be done on this analysis.

Definition at line 193 of file Analysis.hh.

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

00193                                                {
00194       return info().todos();
00195     }

ITree & tree (  )  [protected, inherited]

Access the AIDA tree of the controlling AnalysisHandler object.

Definition at line 55 of file Analysis.cc.

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

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

00055                         {
00056     return handler().tree();
00057   }

virtual std::string year (  )  const [inline, virtual, inherited]

When the original experimental analysis was published.

Definition at line 168 of file Analysis.hh.

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

00168                                    {
00169       return info().year();
00170     }


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 533 of file Analysis.hh.

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

AIDA::IHistogram1D* _hist1MinusT [private]
AIDA::IHistogram1D* _histAplanarity [private]
AIDA::IHistogram1D* _histChMult [private]
AIDA::IHistogram1D* _histCParam [private]
AIDA::IHistogram1D* _histHeavyJetMass [private]
AIDA::IHistogram1D* _histJetRate2Durham [private]

Definition at line 502 of file ALEPH_1996_S3486095.cc.

AIDA::IHistogram1D* _histJetRate3Durham [private]

Definition at line 503 of file ALEPH_1996_S3486095.cc.

AIDA::IHistogram1D* _histJetRate4Durham [private]

Definition at line 504 of file ALEPH_1996_S3486095.cc.

AIDA::IHistogram1D* _histJetRate5Durham [private]

Definition at line 505 of file ALEPH_1996_S3486095.cc.

AIDA::IHistogram1D* _histLogScaledMom [private]
AIDA::IHistogram1D* _histMeanChMult [private]
AIDA::IHistogram1D* _histMeanChMultRapt05 [private]
AIDA::IHistogram1D* _histMeanChMultRapt10 [private]
AIDA::IHistogram1D* _histMeanChMultRapt15 [private]
AIDA::IHistogram1D* _histMeanChMultRapt20 [private]
AIDA::IHistogram1D* _histMeanMultiEta [private]
AIDA::IHistogram1D* _histMeanMultiEtaPrime [private]
AIDA::IHistogram1D* _histMeanMultiK0 [private]
AIDA::IHistogram1D* _histMeanMultiKStar892_0 [private]
AIDA::IHistogram1D* _histMeanMultiKStar892Plus [private]
AIDA::IHistogram1D* _histMeanMultiLambda0 [private]
AIDA::IHistogram1D* _histMeanMultiOmega782 [private]
AIDA::IHistogram1D* _histMeanMultiOmegaOmegaBar [private]
AIDA::IHistogram1D* _histMeanMultiPhi [private]
AIDA::IHistogram1D* _histMeanMultiPi0 [private]
AIDA::IHistogram1D* _histMeanMultiRho [private]
AIDA::IHistogram1D* _histMeanMultiSigma0 [private]
AIDA::IHistogram1D* _histMeanMultiSigma1385Plus [private]
AIDA::IHistogram1D* _histMeanMultiXi1530_0 [private]
AIDA::IHistogram1D* _histMeanMultiXiMinus [private]
AIDA::IHistogram1D* _histMultiEta [private]
AIDA::IHistogram1D* _histMultiEtaPrime [private]
AIDA::IHistogram1D* _histMultiK0 [private]
AIDA::IHistogram1D* _histMultiKPlus [private]
AIDA::IHistogram1D* _histMultiKStar892_0 [private]
AIDA::IHistogram1D* _histMultiKStar892Plus [private]
AIDA::IHistogram1D* _histMultiLambda0 [private]
AIDA::IHistogram1D* _histMultiOmega782 [private]
AIDA::IHistogram1D* _histMultiP [private]
AIDA::IHistogram1D* _histMultiPhi [private]
AIDA::IHistogram1D* _histMultiPhoton [private]
AIDA::IHistogram1D* _histMultiPi0 [private]
AIDA::IHistogram1D* _histMultiPiPlus [private]
AIDA::IHistogram1D* _histMultiRho [private]
AIDA::IHistogram1D* _histMultiSigma1385Plus [private]
AIDA::IHistogram1D* _histMultiXi1530_0 [private]
AIDA::IHistogram1D* _histMultiXiMinus [private]
AIDA::IHistogram1D* _histOblateness [private]
AIDA::IHistogram1D* _histPtSIn [private]
AIDA::IHistogram1D* _histPtSOut [private]
AIDA::IHistogram1D* _histRapidityT [private]
AIDA::IHistogram1D* _histScaledMom [private]
AIDA::IHistogram1D* _histSphericity [private]
AIDA::IHistogram1D* _histTMinor [private]
AIDA::IHistogram1D* _histY3 [private]
shared_ptr<AnalysisInfo> _info [protected, inherited]

Pointer to analysis metadata object.

Definition at line 536 of file Analysis.hh.

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

double _numChParticles [private]
double _weightedTotalNumEta [private]
double _weightedTotalNumEtaPrime [private]
double _weightedTotalNumK0 [private]
double _weightedTotalNumKPlus [private]
double _weightedTotalNumKStar892_0 [private]
double _weightedTotalNumLambda0 [private]
double _weightedTotalNumOmega782 [private]
double _weightedTotalNumP [private]
double _weightedTotalNumPhi [private]
double _weightedTotalNumPhoton [private]
double _weightedTotalNumPi0 [private]
double _weightedTotalNumPiPlus [private]
double _weightedTotalNumRho [private]
double _weightedTotalNumXi1530_0 [private]
double _weightedTotalNumXiMinus [private]
double _weightedTotalPartNum [private]

Store the weighted sums of numbers of charged / charged+neutral particles - used to calculate average number of particles for the inclusive single particle distributions' normalisations.

Definition at line 463 of file ALEPH_1996_S3486095.cc.

Referenced by ALEPH_1996_S3486095::ALEPH_1996_S3486095(), ALEPH_1996_S3486095::analyze(), and ALEPH_1996_S3486095::finalize().


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