PDG_HADRON_MULTIPLICITIES_RATIOS Class Reference

Inheritance diagram for PDG_HADRON_MULTIPLICITIES_RATIOS:

[legend]

Collaboration diagram for PDG_HADRON_MULTIPLICITIES_RATIOS:

[legend]

List of all members.

Detailed Description

Implementation of PDG hadron multiplicities as ratios to pi+- multiplicity.

Author:: Holger Schulz

Definition at line 14 of file PDG_Hadron_Multiplicities_Ratios.cc.


Public Member Functions
	PDG_HADRON_MULTIPLICITIES_RATIOS ()
	Constructor.
virtual bool	isCompatible (const ParticleName &beam1, const ParticleName &beam2) const
	Is this analysis able to run on the supplied pair of beams?
virtual bool	isCompatible (const BeamPair &beams) const
	Is this analysis able to run on the BeamPair beams ?
AnalysisHandler &	handler () const
	Access the controlling AnalysisHandler object.
void	normalize (AIDA::IHistogram1D *&histo, double norm=1.0)
void	scale (AIDA::IHistogram1D *&histo, double scale)
Analysis &	setCrossSection (double xs)
	Set the cross section from the generator.
bool	needsCrossSection () const
	Return true if this analysis needs to know the process cross-section.
Analysis methods
void	analyze (const Event &e)
void	init ()
void	finalize ()
Metadata
Metadata is used for querying from the command line and also for building web pages and the analysis pages in the Rivet manual.
virtual const AnalysisInfo &	info () const
	Get the actual AnalysisInfo object in which all this metadata is stored.
virtual std::string	name () const
	Get the name of the analysis.
virtual std::string	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 const BeamPair	requiredBeams () const
	Return the pair of incoming beams required by this analysis.
virtual const std::vector < std::pair< double, double > > &	energies () const
	Sets of valid beam energy pairs, in GeV.
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	status () const
	Whether this analysis is trusted (in any way!).
Run conditions
const ParticlePair &	beams () const
	Incoming beams for this run.
const BeamPair	beamIds () const
	Incoming beam IDs for this run.
double	sqrtS () const
	Centre of mass energy for this run.
Projection "getting" functions
std::set< ConstProjectionPtr >	getProjections () const
	Get the contained projections, including recursion.
template<typename PROJ>
const PROJ &	getProjection (const std::string &name) const
	Get the named projection, specifying return type via a template argument.
const Projection &	getProjection (const std::string &name) const
Projection applying functions
template<typename PROJ>
const PROJ &	applyProjection (const Event &evt, const PROJ &proj) const
	Apply the supplied projection on event.
template<typename PROJ>
const PROJ &	applyProjection (const Event &evt, const Projection &proj) const
	Apply the supplied projection on event.
template<typename PROJ>
const PROJ &	applyProjection (const Event &evt, const std::string &name) const
	Apply the named projection on event.
Protected Member Functions
Log &	getLog () const
	Get a Log object based on the name() property of the calling analysis object.
double	crossSection () const
	Get the process cross-section in pb. Throws if this hasn't been set.
double	crossSectionPerEvent () const
size_t	numEvents () const
double	sumOfWeights () const
Analysis &	setBeams (const ParticleName &beam1, const ParticleName &beam2)
Analysis &	setNeedsCrossSection (bool needed)
	Declare whether this analysis needs to know the process cross-section from the generator.
ProjectionHandler &	getProjHandler () const
	Get a reference to the ProjectionHandler for this thread.
AIDA analysis infrastructure.
AIDA::IAnalysisFactory &	analysisFactory ()
	Access the AIDA analysis factory of the controlling AnalysisHandler object.
AIDA::ITree &	tree ()
	Access the AIDA tree of the controlling AnalysisHandler object.
AIDA::IHistogramFactory &	histogramFactory ()
	Access the AIDA histogram factory of the controlling AnalysisHandler object.
AIDA::IDataPointSetFactory &	datapointsetFactory ()
	Access the AIDA histogram factory of the controlling AnalysisHandler object.
const std::string	histoDir () const
	Get the canonical histogram "directory" path for this analysis.
const std::string	histoPath (const std::string &hname) const
	Get the canonical histogram path for the named histogram in this analysis.
Internal histogram booking (for use by Analysis sub-classes).
const BinEdges &	binEdges (const std::string &hname) const
	Get bin edges for a named histo (using ref AIDA caching).
const BinEdges &	binEdges (size_t datasetId, size_t xAxisId, size_t yAxisId) const
	Get bin edges for a numbered histo (using ref AIDA caching).
BinEdges	logBinEdges (size_t nbins, double lower, double upper)
	Get bin edges with logarithmic widths.
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D *	bookHistogram1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
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	_weightedTotalNumPiPlus
AIDA::IHistogram1D *	_histMeanMultiPi0
AIDA::IHistogram1D *	_histMeanMultiKPlus
AIDA::IHistogram1D *	_histMeanMultiK0
AIDA::IHistogram1D *	_histMeanMultiEta
AIDA::IHistogram1D *	_histMeanMultiEtaPrime
AIDA::IHistogram1D *	_histMeanMultiDPlus
AIDA::IHistogram1D *	_histMeanMultiD0
AIDA::IHistogram1D *	_histMeanMultiDPlus_s
AIDA::IHistogram1D *	_histMeanMultiBPlus_B0_d
AIDA::IHistogram1D *	_histMeanMultiBPlus_u
AIDA::IHistogram1D *	_histMeanMultiB0_s
AIDA::IHistogram1D *	_histMeanMultiF0_980
AIDA::IHistogram1D *	_histMeanMultiA0_980Plus
AIDA::IHistogram1D *	_histMeanMultiRho770_0
AIDA::IHistogram1D *	_histMeanMultiRho770Plus
AIDA::IHistogram1D *	_histMeanMultiOmega782
AIDA::IHistogram1D *	_histMeanMultiKStar892Plus
AIDA::IHistogram1D *	_histMeanMultiKStar892_0
AIDA::IHistogram1D *	_histMeanMultiPhi1020
AIDA::IHistogram1D *	_histMeanMultiDStar2010Plus
AIDA::IHistogram1D *	_histMeanMultiDStar2007_0
AIDA::IHistogram1D *	_histMeanMultiDStar_s2112Plus
AIDA::IHistogram1D *	_histMeanMultiBStar
AIDA::IHistogram1D *	_histMeanMultiJPsi1S
AIDA::IHistogram1D *	_histMeanMultiPsi2S
AIDA::IHistogram1D *	_histMeanMultiUpsilon1S
AIDA::IHistogram1D *	_histMeanMultiF1_1285
AIDA::IHistogram1D *	_histMeanMultiF1_1420
AIDA::IHistogram1D *	_histMeanMultiChi_c1_3510
AIDA::IHistogram1D *	_histMeanMultiF2_1270
AIDA::IHistogram1D *	_histMeanMultiF2Prime1525
AIDA::IHistogram1D *	_histMeanMultiK2Star1430Plus
AIDA::IHistogram1D *	_histMeanMultiK2Star1430_0
AIDA::IHistogram1D *	_histMeanMultiBStarStar
AIDA::IHistogram1D *	_histMeanMultiDs1Plus
AIDA::IHistogram1D *	_histMeanMultiDs2Plus
AIDA::IHistogram1D *	_histMeanMultiP
AIDA::IHistogram1D *	_histMeanMultiLambda
AIDA::IHistogram1D *	_histMeanMultiSigma0
AIDA::IHistogram1D *	_histMeanMultiSigmaMinus
AIDA::IHistogram1D *	_histMeanMultiSigmaPlus
AIDA::IHistogram1D *	_histMeanMultiSigmaPlusMinus
AIDA::IHistogram1D *	_histMeanMultiXiMinus
AIDA::IHistogram1D *	_histMeanMultiDelta1232PlusPlus
AIDA::IHistogram1D *	_histMeanMultiSigma1385Minus
AIDA::IHistogram1D *	_histMeanMultiSigma1385Plus
AIDA::IHistogram1D *	_histMeanMultiSigma1385PlusMinus
AIDA::IHistogram1D *	_histMeanMultiXi1530_0
AIDA::IHistogram1D *	_histMeanMultiOmegaMinus
AIDA::IHistogram1D *	_histMeanMultiLambda_c_Plus
AIDA::IHistogram1D *	_histMeanMultiLambda_b_0
AIDA::IHistogram1D *	_histMeanMultiSigma_c_PlusPlus_0
AIDA::IHistogram1D *	_histMeanMultiLambda1520
Friends
class	Projectionhandler

Constructor & Destructor Documentation

PDG_HADRON_MULTIPLICITIES_RATIOS ( ) [inline]

Constructor.

Definition at line 18 of file PDG_Hadron_Multiplicities_Ratios.cc.

References PDG_HADRON_MULTIPLICITIES_RATIOS::_weightedTotalNumPiPlus, Rivet::ELECTRON, Rivet::POSITRON, and Analysis::setBeams().

00018                                        : Analysis("PDG_HADRON_MULTIPLICITIES_RATIOS")
00019     {
00020       setBeams(ELECTRON, POSITRON);
00021       _weightedTotalNumPiPlus = 0;
00022     }

Member Function Documentation

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 28 of file PDG_Hadron_Multiplicities_Ratios.cc.

00028                                  {
00029       // First, veto on leptonic events by requiring at least 4 charged FS particles
00030       const FinalState& fs = applyProjection<FinalState>(e, "FS");
00031       const size_t numParticles = fs.particles().size();
00032 
00033       // Even if we only generate hadronic events, we still need a cut on numCharged >= 2.
00034       if (numParticles < 2) {
00035         getLog() << Log::DEBUG << "Failed leptonic event cut" << endl;
00036         vetoEvent;
00037       }
00038       getLog() << Log::DEBUG << "Passed leptonic event cut" << endl;
00039 
00040       // Get event weight for histo filling
00041       const double weight = e.weight();
00042 
00043       getLog() << Log::DEBUG << "sqrt(S) = " << sqrtS()/GeV << " GeV" << endl;
00044 
00045       // Final state of unstable particles to get particle spectra
00046       const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(e, "UFS");
00047 
00048       if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
00049         foreach (const Particle& p, ufs.particles()) {
00050           const PdgId id = abs(p.pdgId());
00051           switch (id) {
00052              case 211:
00053                 _weightedTotalNumPiPlus += weight;
00054                 break;
00055              case 111:
00056                 _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight);
00057                 break;
00058              case 321:
00059                 _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight);
00060                 break;
00061              case 130:
00062              case 310:
00063                 _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight);
00064                 break;
00065              case 221:
00066                 _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight);
00067                 break;
00068              case 331:
00069                 _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight);
00070                 break;
00071              case 411:
00072                 _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight);
00073                 break;
00074              case 421:
00075                 _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight);
00076                 break;
00077              case 431:
00078                 _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight);
00079                 break;
00080              case 9010221:
00081                 _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight);
00082                 break;
00083              case 113:
00084                 _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight);
00085                 break;
00086              case 223:
00087                 _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight);
00088                 break;
00089              case 323:
00090                 _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight);
00091                 break;
00092              case 313:
00093                 _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight);
00094                 break;
00095              case 333:
00096                 _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight);
00097                 break;
00098              case 413:
00099                 _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight);
00100                 break;
00101              case 423:
00102                 _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight);
00103                 break;
00104              case 433:
00105                 _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight);
00106                 break;
00107              case 443:
00108                 _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight);
00109                 break;
00110              case 225:
00111                 _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight);
00112                 break;
00113              case 2212:
00114                 _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight);
00115                 break;
00116              case 3122:
00117                 _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight);
00118                 break;
00119              case 3212:
00120                 _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight);
00121                 break;
00122              case 3312:
00123                 _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight);
00124                 break;
00125              case 2224:
00126                 _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight);
00127                 break;
00128              case 3114:
00129                 _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight);
00130                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00131                 break;
00132              case 3224:
00133                 _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight);
00134                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00135                 break;
00136              case 3324:
00137                 _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight);
00138                 break;
00139              case 3334:
00140                 _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight);
00141                 break;
00142              case 4122:
00143                 _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight);
00144                 break;
00145              case 4222:
00146              case 4112:
00147                 _histMeanMultiSigma_c_PlusPlus_0->fill(_histMeanMultiSigma_c_PlusPlus_0->binMean(0), weight);
00148                 break;
00149              case 3124:
00150                 _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight);
00151                 break;
00152           }
00153         }
00154       }
00155 
00156       if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
00157         foreach (const Particle& p, ufs.particles()) {
00158           const PdgId id = abs(p.pdgId());
00159           switch (id) {
00160              case 211:
00161                 _weightedTotalNumPiPlus += weight;
00162                 break;
00163              case 111:
00164                 _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight);
00165                 break;
00166              case 321:
00167                 _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight);
00168                 break;
00169              case 130:
00170              case 310:
00171                 _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight);
00172                 break;
00173              case 221:
00174                 _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight);
00175                 break;
00176              case 331:
00177                 _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight);
00178                 break;
00179              case 411:
00180                 _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight);
00181                 break;
00182              case 421:
00183                 _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight);
00184                 break;
00185              case 431:
00186                 _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight);
00187                 break;
00188              case 9010221:
00189                 _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight);
00190                 break;
00191              case 113:
00192                 _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight);
00193                 break;
00194              case 323:
00195                 _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight);
00196                 break;
00197              case 313:
00198                 _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight);
00199                 break;
00200              case 333:
00201                 _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight);
00202                 break;
00203              case 413:
00204                 _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight);
00205                 break;
00206              case 423:
00207                 _histMeanMultiDStar2007_0->fill(_histMeanMultiDStar2007_0->binMean(0), weight);
00208                 break;
00209              case 225:
00210                 _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight);
00211                 break;
00212              case 325:
00213                 _histMeanMultiK2Star1430Plus->fill(_histMeanMultiK2Star1430Plus->binMean(0), weight);
00214                 break;
00215              case 315:
00216                 _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight);
00217                 break;
00218              case 2212:
00219                 _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight);
00220                 break;
00221              case 3122:
00222                 _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight);
00223                 break;
00224              case 3312:
00225                 _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight);
00226                 break;
00227              case 3114:
00228                 _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight);
00229                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00230                 break;
00231              case 3224:
00232                 _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight);
00233                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00234                 break;
00235              case 3334:
00236                 _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight);
00237                 break;
00238              case 4122:
00239                 _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight);
00240                 break;
00241           }
00242         }
00243       }
00244 
00245       if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
00246         foreach (const Particle& p, ufs.particles()) {
00247           const PdgId id = abs(p.pdgId());
00248           switch (id) {
00249              case 211:
00250                 _weightedTotalNumPiPlus += weight;
00251                 break;
00252              case 111:
00253                 _histMeanMultiPi0->fill(_histMeanMultiPi0->binMean(0), weight);
00254                 break;
00255              case 321:
00256                 _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight);
00257                 break;
00258              case 130:
00259              case 310:
00260                 _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight);
00261                 break;
00262              case 221:
00263                 _histMeanMultiEta->fill(_histMeanMultiEta->binMean(0), weight);
00264                 break;
00265              case 331:
00266                 _histMeanMultiEtaPrime->fill(_histMeanMultiEtaPrime->binMean(0), weight);
00267                 break;
00268              case 411:
00269                 _histMeanMultiDPlus->fill(_histMeanMultiDPlus->binMean(0), weight);
00270                 break;
00271              case 421:
00272                 _histMeanMultiD0->fill(_histMeanMultiD0->binMean(0), weight);
00273                 break;
00274              case 431:
00275                 _histMeanMultiDPlus_s->fill(_histMeanMultiDPlus_s->binMean(0), weight);
00276                 break;
00277              case 511:
00278                 _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight);
00279                 break;
00280              case 521:
00281                 _histMeanMultiBPlus_B0_d->fill(_histMeanMultiBPlus_B0_d->binMean(0), weight);
00282                 _histMeanMultiBPlus_u->fill(_histMeanMultiBPlus_u->binMean(0), weight);
00283                 break;
00284              case 531:
00285                 _histMeanMultiB0_s->fill(_histMeanMultiB0_s->binMean(0), weight);
00286                 break;
00287              case 9010221:
00288                 _histMeanMultiF0_980->fill(_histMeanMultiF0_980->binMean(0), weight);
00289                 break;
00290              case 9000211:
00291                 _histMeanMultiA0_980Plus->fill(_histMeanMultiA0_980Plus->binMean(0), weight);
00292                 break;
00293              case 113:
00294                 _histMeanMultiRho770_0->fill(_histMeanMultiRho770_0->binMean(0), weight);
00295                 break;
00296              case 213:
00297                 _histMeanMultiRho770Plus->fill(_histMeanMultiRho770Plus->binMean(0), weight);
00298                 break;
00299              case 223:
00300                 _histMeanMultiOmega782->fill(_histMeanMultiOmega782->binMean(0), weight);
00301                 break;
00302              case 323:
00303                 _histMeanMultiKStar892Plus->fill(_histMeanMultiKStar892Plus->binMean(0), weight);
00304                 break;
00305              case 313:
00306                 _histMeanMultiKStar892_0->fill(_histMeanMultiKStar892_0->binMean(0), weight);
00307                 break;
00308              case 333:
00309                 _histMeanMultiPhi1020->fill(_histMeanMultiPhi1020->binMean(0), weight);
00310                 break;
00311              case 413:
00312                 _histMeanMultiDStar2010Plus->fill(_histMeanMultiDStar2010Plus->binMean(0), weight);
00313                 break;
00314              case 433:
00315                 _histMeanMultiDStar_s2112Plus->fill(_histMeanMultiDStar_s2112Plus->binMean(0), weight);
00316                 break;
00317              case 513:
00318              case 523:
00319              case 533:
00320                 _histMeanMultiBStar->fill(_histMeanMultiBStar->binMean(0), weight);
00321                 break;
00322              case 443:
00323                 _histMeanMultiJPsi1S->fill(_histMeanMultiJPsi1S->binMean(0), weight);
00324                 break;
00325              case 100443:
00326                 _histMeanMultiPsi2S->fill(_histMeanMultiPsi2S->binMean(0), weight);
00327                 break;
00328              case 553:
00329                 _histMeanMultiUpsilon1S->fill(_histMeanMultiUpsilon1S->binMean(0), weight);
00330                 break;
00331              case 20223:
00332                 _histMeanMultiF1_1285->fill(_histMeanMultiF1_1285->binMean(0), weight);
00333                 break;
00334              case 20333:
00335                 _histMeanMultiF1_1420->fill(_histMeanMultiF1_1420->binMean(0), weight);
00336                 break;
00337              case 445:
00338                 _histMeanMultiChi_c1_3510->fill(_histMeanMultiChi_c1_3510->binMean(0), weight);
00339                 break;
00340              case 225:
00341                 _histMeanMultiF2_1270->fill(_histMeanMultiF2_1270->binMean(0), weight);
00342                 break;
00343              case 335:
00344                 _histMeanMultiF2Prime1525->fill(_histMeanMultiF2Prime1525->binMean(0), weight);
00345                 break;
00346              case 315:
00347                 _histMeanMultiK2Star1430_0->fill(_histMeanMultiK2Star1430_0->binMean(0), weight);
00348                 break;
00349              case 515:
00350              case 525:
00351              case 535:
00352                 _histMeanMultiBStarStar->fill(_histMeanMultiBStarStar->binMean(0), weight);
00353                 break;
00354              case 10433:
00355              case 20433:
00356                 _histMeanMultiDs1Plus->fill(_histMeanMultiDs1Plus->binMean(0), weight);
00357                 break;
00358              case 435:
00359                 _histMeanMultiDs2Plus->fill(_histMeanMultiDs2Plus->binMean(0), weight);
00360                 break;
00361              case 2212:
00362                 _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight);
00363                 break;
00364              case 3122:
00365                 _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight);
00366                 break;
00367              case 3212:
00368                 _histMeanMultiSigma0->fill(_histMeanMultiSigma0->binMean(0), weight);
00369                 break;
00370              case 3112:
00371                 _histMeanMultiSigmaMinus->fill(_histMeanMultiSigmaMinus->binMean(0), weight);
00372                 _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight);
00373                 break;
00374              case 3222:
00375                 _histMeanMultiSigmaPlus->fill(_histMeanMultiSigmaPlus->binMean(0), weight);
00376                 _histMeanMultiSigmaPlusMinus->fill(_histMeanMultiSigmaPlusMinus->binMean(0), weight);
00377                 break;
00378              case 3312:
00379                 _histMeanMultiXiMinus->fill(_histMeanMultiXiMinus->binMean(0), weight);
00380                 break;
00381              case 2224:
00382                 _histMeanMultiDelta1232PlusPlus->fill(_histMeanMultiDelta1232PlusPlus->binMean(0), weight);
00383                 break;
00384              case 3114:
00385                 _histMeanMultiSigma1385Minus->fill(_histMeanMultiSigma1385Minus->binMean(0), weight);
00386                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00387                 break;
00388              case 3224:
00389                 _histMeanMultiSigma1385Plus->fill(_histMeanMultiSigma1385Plus->binMean(0), weight);
00390                 _histMeanMultiSigma1385PlusMinus->fill(_histMeanMultiSigma1385PlusMinus->binMean(0), weight);
00391                 break;
00392              case 3324:
00393                 _histMeanMultiXi1530_0->fill(_histMeanMultiXi1530_0->binMean(0), weight);
00394                 break;
00395              case 3334:
00396                 _histMeanMultiOmegaMinus->fill(_histMeanMultiOmegaMinus->binMean(0), weight);
00397                 break;
00398              case 4122:
00399                 _histMeanMultiLambda_c_Plus->fill(_histMeanMultiLambda_c_Plus->binMean(0), weight);
00400                 break;
00401              case 5122:
00402                 _histMeanMultiLambda_b_0->fill(_histMeanMultiLambda_b_0->binMean(0), weight);
00403                 break;
00404              case 3124:
00405                 _histMeanMultiLambda1520->fill(_histMeanMultiLambda1520->binMean(0), weight);
00406                 break;
00407           }
00408         }
00409       }
00410 
00411       if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
00412         foreach (const Particle& p, ufs.particles()) {
00413           const PdgId id = abs(p.pdgId());
00414           switch (id) {
00415              case 211:
00416                 _weightedTotalNumPiPlus += weight;
00417                 break;
00418              case 321:
00419                 _histMeanMultiKPlus->fill(_histMeanMultiKPlus->binMean(0), weight);
00420                 break;
00421              case 130:
00422              case 310:
00423                 _histMeanMultiK0->fill(_histMeanMultiK0->binMean(0), weight);
00424                 break;
00425              case 2212:
00426                 _histMeanMultiP->fill(_histMeanMultiP->binMean(0), weight);
00427                 break;
00428              case 3122:
00429                 _histMeanMultiLambda->fill(_histMeanMultiLambda->binMean(0), weight);
00430                 break;
00431           }
00432         }
00433       }
00434 
00435     }

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 439 of file PDG_Hadron_Multiplicities_Ratios.cc.

00439                 {
00440       addProjection(ChargedFinalState(), "FS");
00441       addProjection(UnstableFinalState(), "UFS");
00442 
00443       if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
00444         _histMeanMultiPi0                = bookHistogram1D( 2, 1, 1);
00445         _histMeanMultiKPlus              = bookHistogram1D( 3, 1, 1);
00446         _histMeanMultiK0                 = bookHistogram1D( 4, 1, 1);
00447         _histMeanMultiEta                = bookHistogram1D( 5, 1, 1);
00448         _histMeanMultiEtaPrime           = bookHistogram1D( 6, 1, 1);
00449         _histMeanMultiDPlus              = bookHistogram1D( 7, 1, 1);
00450         _histMeanMultiD0                 = bookHistogram1D( 8, 1, 1);
00451         _histMeanMultiDPlus_s            = bookHistogram1D( 9, 1, 1);
00452         _histMeanMultiF0_980             = bookHistogram1D(13, 1, 1);
00453         _histMeanMultiRho770_0           = bookHistogram1D(15, 1, 1);
00454         _histMeanMultiOmega782           = bookHistogram1D(17, 1, 1);
00455         _histMeanMultiKStar892Plus       = bookHistogram1D(18, 1, 1);
00456         _histMeanMultiKStar892_0         = bookHistogram1D(19, 1, 1);
00457         _histMeanMultiPhi1020            = bookHistogram1D(20, 1, 1);
00458         _histMeanMultiDStar2010Plus      = bookHistogram1D(21, 1, 1);
00459         _histMeanMultiDStar2007_0        = bookHistogram1D(22, 1, 1);
00460         _histMeanMultiDStar_s2112Plus    = bookHistogram1D(23, 1, 1);
00461         _histMeanMultiJPsi1S             = bookHistogram1D(25, 1, 1);
00462         _histMeanMultiF2_1270            = bookHistogram1D(31, 1, 1);
00463         _histMeanMultiP                  = bookHistogram1D(38, 1, 1);
00464         _histMeanMultiLambda             = bookHistogram1D(39, 1, 1);
00465         _histMeanMultiSigma0             = bookHistogram1D(40, 1, 1);
00466         _histMeanMultiXiMinus            = bookHistogram1D(44, 1, 1);
00467         _histMeanMultiDelta1232PlusPlus  = bookHistogram1D(45, 1, 1);
00468         _histMeanMultiSigma1385Minus     = bookHistogram1D(46, 1, 1);
00469         _histMeanMultiSigma1385Plus      = bookHistogram1D(47, 1, 1);
00470         _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 1);
00471         _histMeanMultiXi1530_0           = bookHistogram1D(49, 1, 1);
00472         _histMeanMultiOmegaMinus         = bookHistogram1D(50, 1, 1);
00473         _histMeanMultiLambda_c_Plus      = bookHistogram1D(51, 1, 1);
00474         _histMeanMultiSigma_c_PlusPlus_0 = bookHistogram1D(53, 1, 1);
00475         _histMeanMultiLambda1520         = bookHistogram1D(54, 1, 1);
00476       }
00477 
00478       if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
00479         _histMeanMultiPi0                = bookHistogram1D( 2, 1, 2);
00480         _histMeanMultiKPlus              = bookHistogram1D( 3, 1, 2);
00481         _histMeanMultiK0                 = bookHistogram1D( 4, 1, 2);
00482         _histMeanMultiEta                = bookHistogram1D( 5, 1, 2);
00483         _histMeanMultiEtaPrime           = bookHistogram1D( 6, 1, 2);
00484         _histMeanMultiDPlus              = bookHistogram1D( 7, 1, 2);
00485         _histMeanMultiD0                 = bookHistogram1D( 8, 1, 2);
00486         _histMeanMultiDPlus_s            = bookHistogram1D( 9, 1, 2);
00487         _histMeanMultiF0_980             = bookHistogram1D(13, 1, 2);
00488         _histMeanMultiRho770_0           = bookHistogram1D(15, 1, 2);
00489         _histMeanMultiKStar892Plus       = bookHistogram1D(18, 1, 2);
00490         _histMeanMultiKStar892_0         = bookHistogram1D(19, 1, 2);
00491         _histMeanMultiPhi1020            = bookHistogram1D(20, 1, 2);
00492         _histMeanMultiDStar2010Plus      = bookHistogram1D(21, 1, 2);
00493         _histMeanMultiDStar2007_0        = bookHistogram1D(22, 1, 2);
00494         _histMeanMultiF2_1270            = bookHistogram1D(31, 1, 2);
00495         _histMeanMultiK2Star1430Plus     = bookHistogram1D(33, 1, 1);
00496         _histMeanMultiK2Star1430_0       = bookHistogram1D(34, 1, 1);
00497         _histMeanMultiP                  = bookHistogram1D(38, 1, 2);
00498         _histMeanMultiLambda             = bookHistogram1D(39, 1, 2);
00499         _histMeanMultiXiMinus            = bookHistogram1D(44, 1, 2);
00500         _histMeanMultiSigma1385Minus     = bookHistogram1D(46, 1, 2);
00501         _histMeanMultiSigma1385Plus      = bookHistogram1D(47, 1, 2);
00502         _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 2);
00503         _histMeanMultiOmegaMinus         = bookHistogram1D(50, 1, 2);
00504         _histMeanMultiLambda_c_Plus      = bookHistogram1D(51, 1, 2);
00505       }
00506 
00507       if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
00508         _histMeanMultiPi0                = bookHistogram1D( 2, 1, 3);
00509         _histMeanMultiKPlus              = bookHistogram1D( 3, 1, 3);
00510         _histMeanMultiK0                 = bookHistogram1D( 4, 1, 3);
00511         _histMeanMultiEta                = bookHistogram1D( 5, 1, 3);
00512         _histMeanMultiEtaPrime           = bookHistogram1D( 6, 1, 3);
00513         _histMeanMultiDPlus              = bookHistogram1D( 7, 1, 3);
00514         _histMeanMultiD0                 = bookHistogram1D( 8, 1, 3);
00515         _histMeanMultiDPlus_s            = bookHistogram1D( 9, 1, 3);
00516         _histMeanMultiBPlus_B0_d         = bookHistogram1D(10, 1, 1);
00517         _histMeanMultiBPlus_u            = bookHistogram1D(11, 1, 1);
00518         _histMeanMultiB0_s               = bookHistogram1D(12, 1, 1);
00519         _histMeanMultiF0_980             = bookHistogram1D(13, 1, 3);
00520         _histMeanMultiA0_980Plus         = bookHistogram1D(14, 1, 1);
00521         _histMeanMultiRho770_0           = bookHistogram1D(15, 1, 3);
00522         _histMeanMultiRho770Plus         = bookHistogram1D(16, 1, 1);
00523         _histMeanMultiOmega782           = bookHistogram1D(17, 1, 2);
00524         _histMeanMultiKStar892Plus       = bookHistogram1D(18, 1, 3);
00525         _histMeanMultiKStar892_0         = bookHistogram1D(19, 1, 3);
00526         _histMeanMultiPhi1020            = bookHistogram1D(20, 1, 3);
00527         _histMeanMultiDStar2010Plus      = bookHistogram1D(21, 1, 3);
00528         _histMeanMultiDStar_s2112Plus    = bookHistogram1D(23, 1, 2);
00529         _histMeanMultiBStar              = bookHistogram1D(24, 1, 1);
00530         _histMeanMultiJPsi1S             = bookHistogram1D(25, 1, 2);
00531         _histMeanMultiPsi2S              = bookHistogram1D(26, 1, 1);
00532         _histMeanMultiUpsilon1S          = bookHistogram1D(27, 1, 1);
00533         _histMeanMultiF1_1285            = bookHistogram1D(28, 1, 1);
00534         _histMeanMultiF1_1420            = bookHistogram1D(29, 1, 1);
00535         _histMeanMultiChi_c1_3510        = bookHistogram1D(30, 1, 1);
00536         _histMeanMultiF2_1270            = bookHistogram1D(31, 1, 3);
00537         _histMeanMultiF2Prime1525        = bookHistogram1D(32, 1, 1);
00538         _histMeanMultiK2Star1430_0       = bookHistogram1D(34, 1, 2);
00539         _histMeanMultiBStarStar          = bookHistogram1D(35, 1, 1);
00540         _histMeanMultiDs1Plus            = bookHistogram1D(36, 1, 1);
00541         _histMeanMultiDs2Plus            = bookHistogram1D(37, 1, 1);
00542         _histMeanMultiP                  = bookHistogram1D(38, 1, 3);
00543         _histMeanMultiLambda             = bookHistogram1D(39, 1, 3);
00544         _histMeanMultiSigma0             = bookHistogram1D(40, 1, 2);
00545         _histMeanMultiSigmaMinus         = bookHistogram1D(41, 1, 1);
00546         _histMeanMultiSigmaPlus          = bookHistogram1D(42, 1, 1);
00547         _histMeanMultiSigmaPlusMinus     = bookHistogram1D(43, 1, 1);
00548         _histMeanMultiXiMinus            = bookHistogram1D(44, 1, 3);
00549         _histMeanMultiDelta1232PlusPlus  = bookHistogram1D(45, 1, 2);
00550         _histMeanMultiSigma1385Minus     = bookHistogram1D(46, 1, 3);
00551         _histMeanMultiSigma1385Plus      = bookHistogram1D(47, 1, 3);
00552         _histMeanMultiSigma1385PlusMinus = bookHistogram1D(48, 1, 3);
00553         _histMeanMultiXi1530_0           = bookHistogram1D(49, 1, 2);
00554         _histMeanMultiOmegaMinus         = bookHistogram1D(50, 1, 3);
00555         _histMeanMultiLambda_c_Plus      = bookHistogram1D(51, 1, 3);
00556         _histMeanMultiLambda_b_0         = bookHistogram1D(52, 1, 1);
00557         _histMeanMultiLambda1520         = bookHistogram1D(54, 1, 2);
00558       }
00559 
00560       if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
00561         _histMeanMultiKPlus             = bookHistogram1D( 3, 1, 4);
00562         _histMeanMultiK0                = bookHistogram1D( 4, 1, 4);
00563         _histMeanMultiP                 = bookHistogram1D(38, 1, 4);
00564         _histMeanMultiLambda            = bookHistogram1D(39, 1, 4);
00565       }
00566     }

void finalize ( ) [inline, virtual]

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

Implements Analysis.

Definition at line 571 of file PDG_Hadron_Multiplicities_Ratios.cc.

00571                     {
00572       if (sqrtS()/GeV >= 9.5 && sqrtS()/GeV <= 10.5) {
00573         scale(_histMeanMultiPi0               , 1.0/_weightedTotalNumPiPlus);
00574         scale(_histMeanMultiKPlus             , 1.0/_weightedTotalNumPiPlus);
00575         scale(_histMeanMultiK0                , 1.0/_weightedTotalNumPiPlus);
00576         scale(_histMeanMultiEta               , 1.0/_weightedTotalNumPiPlus);
00577         scale(_histMeanMultiEtaPrime          , 1.0/_weightedTotalNumPiPlus);
00578         scale(_histMeanMultiDPlus             , 1.0/_weightedTotalNumPiPlus);
00579         scale(_histMeanMultiD0                , 1.0/_weightedTotalNumPiPlus);
00580         scale(_histMeanMultiDPlus_s           , 1.0/_weightedTotalNumPiPlus);
00581         scale(_histMeanMultiF0_980            , 1.0/_weightedTotalNumPiPlus);
00582         scale(_histMeanMultiRho770_0          , 1.0/_weightedTotalNumPiPlus);
00583         scale(_histMeanMultiOmega782          , 1.0/_weightedTotalNumPiPlus);
00584         scale(_histMeanMultiKStar892Plus      , 1.0/_weightedTotalNumPiPlus);
00585         scale(_histMeanMultiKStar892_0        , 1.0/_weightedTotalNumPiPlus);
00586         scale(_histMeanMultiPhi1020           , 1.0/_weightedTotalNumPiPlus);
00587         scale(_histMeanMultiDStar2010Plus     , 1.0/_weightedTotalNumPiPlus);
00588         scale(_histMeanMultiDStar2007_0       , 1.0/_weightedTotalNumPiPlus);
00589         scale(_histMeanMultiDStar_s2112Plus   , 1.0/_weightedTotalNumPiPlus);
00590         scale(_histMeanMultiJPsi1S            , 1.0/_weightedTotalNumPiPlus);
00591         scale(_histMeanMultiF2_1270           , 1.0/_weightedTotalNumPiPlus);
00592         scale(_histMeanMultiP                 , 1.0/_weightedTotalNumPiPlus);
00593         scale(_histMeanMultiLambda            , 1.0/_weightedTotalNumPiPlus);
00594         scale(_histMeanMultiSigma0            , 1.0/_weightedTotalNumPiPlus);
00595         scale(_histMeanMultiXiMinus           , 1.0/_weightedTotalNumPiPlus);
00596         scale(_histMeanMultiDelta1232PlusPlus , 1.0/_weightedTotalNumPiPlus);
00597         scale(_histMeanMultiSigma1385Minus    , 1.0/_weightedTotalNumPiPlus);
00598         scale(_histMeanMultiSigma1385Plus     , 1.0/_weightedTotalNumPiPlus);
00599         scale(_histMeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus);
00600         scale(_histMeanMultiXi1530_0          , 1.0/_weightedTotalNumPiPlus);
00601         scale(_histMeanMultiOmegaMinus        , 1.0/_weightedTotalNumPiPlus);
00602         scale(_histMeanMultiLambda_c_Plus     , 1.0/_weightedTotalNumPiPlus);
00603         scale(_histMeanMultiSigma_c_PlusPlus_0, 1.0/_weightedTotalNumPiPlus);
00604         scale(_histMeanMultiLambda1520        , 1.0/_weightedTotalNumPiPlus);
00605       }
00606 
00607       if (sqrtS()/GeV >= 29 && sqrtS()/GeV <= 35) {
00608         scale(_histMeanMultiPi0               , 1.0/_weightedTotalNumPiPlus);
00609         scale(_histMeanMultiKPlus             , 1.0/_weightedTotalNumPiPlus);
00610         scale(_histMeanMultiK0                , 1.0/_weightedTotalNumPiPlus);
00611         scale(_histMeanMultiEta               , 1.0/_weightedTotalNumPiPlus);
00612         scale(_histMeanMultiEtaPrime          , 1.0/_weightedTotalNumPiPlus);
00613         scale(_histMeanMultiDPlus             , 1.0/_weightedTotalNumPiPlus);
00614         scale(_histMeanMultiD0                , 1.0/_weightedTotalNumPiPlus);
00615         scale(_histMeanMultiDPlus_s           , 1.0/_weightedTotalNumPiPlus);
00616         scale(_histMeanMultiF0_980            , 1.0/_weightedTotalNumPiPlus);
00617         scale(_histMeanMultiRho770_0          , 1.0/_weightedTotalNumPiPlus);
00618         scale(_histMeanMultiKStar892Plus      , 1.0/_weightedTotalNumPiPlus);
00619         scale(_histMeanMultiKStar892_0        , 1.0/_weightedTotalNumPiPlus);
00620         scale(_histMeanMultiPhi1020           , 1.0/_weightedTotalNumPiPlus);
00621         scale(_histMeanMultiDStar2010Plus     , 1.0/_weightedTotalNumPiPlus);
00622         scale(_histMeanMultiDStar2007_0       , 1.0/_weightedTotalNumPiPlus);
00623         scale(_histMeanMultiF2_1270           , 1.0/_weightedTotalNumPiPlus);
00624         scale(_histMeanMultiK2Star1430Plus    , 1.0/_weightedTotalNumPiPlus);
00625         scale(_histMeanMultiK2Star1430_0      , 1.0/_weightedTotalNumPiPlus);
00626         scale(_histMeanMultiP                 , 1.0/_weightedTotalNumPiPlus);
00627         scale(_histMeanMultiLambda            , 1.0/_weightedTotalNumPiPlus);
00628         scale(_histMeanMultiXiMinus           , 1.0/_weightedTotalNumPiPlus);
00629         scale(_histMeanMultiSigma1385Minus    , 1.0/_weightedTotalNumPiPlus);
00630         scale(_histMeanMultiSigma1385Plus     , 1.0/_weightedTotalNumPiPlus);
00631         scale(_histMeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus);
00632         scale(_histMeanMultiOmegaMinus        , 1.0/_weightedTotalNumPiPlus);
00633         scale(_histMeanMultiLambda_c_Plus     , 1.0/_weightedTotalNumPiPlus);
00634       }
00635 
00636       if (sqrtS()/GeV >= 89.5 && sqrtS()/GeV <= 91.8) {
00637         scale(_histMeanMultiPi0               , 1.0/_weightedTotalNumPiPlus);
00638         scale(_histMeanMultiKPlus             , 1.0/_weightedTotalNumPiPlus);
00639         scale(_histMeanMultiK0                , 1.0/_weightedTotalNumPiPlus);
00640         scale(_histMeanMultiEta               , 1.0/_weightedTotalNumPiPlus);
00641         scale(_histMeanMultiEtaPrime          , 1.0/_weightedTotalNumPiPlus);
00642         scale(_histMeanMultiDPlus             , 1.0/_weightedTotalNumPiPlus);
00643         scale(_histMeanMultiD0                , 1.0/_weightedTotalNumPiPlus);
00644         scale(_histMeanMultiDPlus_s           , 1.0/_weightedTotalNumPiPlus);
00645         scale(_histMeanMultiBPlus_B0_d        , 1.0/_weightedTotalNumPiPlus);
00646         scale(_histMeanMultiBPlus_u           , 1.0/_weightedTotalNumPiPlus);
00647         scale(_histMeanMultiB0_s              , 1.0/_weightedTotalNumPiPlus);
00648         scale(_histMeanMultiF0_980            , 1.0/_weightedTotalNumPiPlus);
00649         scale(_histMeanMultiA0_980Plus        , 1.0/_weightedTotalNumPiPlus);
00650         scale(_histMeanMultiRho770_0          , 1.0/_weightedTotalNumPiPlus);
00651         scale(_histMeanMultiRho770Plus        , 1.0/_weightedTotalNumPiPlus);
00652         scale(_histMeanMultiOmega782          , 1.0/_weightedTotalNumPiPlus);
00653         scale(_histMeanMultiKStar892Plus      , 1.0/_weightedTotalNumPiPlus);
00654         scale(_histMeanMultiKStar892_0        , 1.0/_weightedTotalNumPiPlus);
00655         scale(_histMeanMultiPhi1020           , 1.0/_weightedTotalNumPiPlus);
00656         scale(_histMeanMultiDStar2010Plus     , 1.0/_weightedTotalNumPiPlus);
00657         scale(_histMeanMultiDStar_s2112Plus   , 1.0/_weightedTotalNumPiPlus);
00658         scale(_histMeanMultiBStar             , 1.0/_weightedTotalNumPiPlus);
00659         scale(_histMeanMultiJPsi1S            , 1.0/_weightedTotalNumPiPlus);
00660         scale(_histMeanMultiPsi2S             , 1.0/_weightedTotalNumPiPlus);
00661         scale(_histMeanMultiUpsilon1S         , 1.0/_weightedTotalNumPiPlus);
00662         scale(_histMeanMultiF1_1285           , 1.0/_weightedTotalNumPiPlus);
00663         scale(_histMeanMultiF1_1420           , 1.0/_weightedTotalNumPiPlus);
00664         scale(_histMeanMultiChi_c1_3510       , 1.0/_weightedTotalNumPiPlus);
00665         scale(_histMeanMultiF2_1270           , 1.0/_weightedTotalNumPiPlus);
00666         scale(_histMeanMultiF2Prime1525       , 1.0/_weightedTotalNumPiPlus);
00667         scale(_histMeanMultiK2Star1430_0      , 1.0/_weightedTotalNumPiPlus);
00668         scale(_histMeanMultiBStarStar         , 1.0/_weightedTotalNumPiPlus);
00669         scale(_histMeanMultiDs1Plus           , 1.0/_weightedTotalNumPiPlus);
00670         scale(_histMeanMultiDs2Plus           , 1.0/_weightedTotalNumPiPlus);
00671         scale(_histMeanMultiP                 , 1.0/_weightedTotalNumPiPlus);
00672         scale(_histMeanMultiLambda            , 1.0/_weightedTotalNumPiPlus);
00673         scale(_histMeanMultiSigma0            , 1.0/_weightedTotalNumPiPlus);
00674         scale(_histMeanMultiSigmaMinus        , 1.0/_weightedTotalNumPiPlus);
00675         scale(_histMeanMultiSigmaPlus         , 1.0/_weightedTotalNumPiPlus);
00676         scale(_histMeanMultiSigmaPlusMinus    , 1.0/_weightedTotalNumPiPlus);
00677         scale(_histMeanMultiXiMinus           , 1.0/_weightedTotalNumPiPlus);
00678         scale(_histMeanMultiDelta1232PlusPlus , 1.0/_weightedTotalNumPiPlus);
00679         scale(_histMeanMultiSigma1385Minus    , 1.0/_weightedTotalNumPiPlus);
00680         scale(_histMeanMultiSigma1385Plus     , 1.0/_weightedTotalNumPiPlus);
00681         scale(_histMeanMultiSigma1385PlusMinus, 1.0/_weightedTotalNumPiPlus);
00682         scale(_histMeanMultiXi1530_0          , 1.0/_weightedTotalNumPiPlus);
00683         scale(_histMeanMultiOmegaMinus        , 1.0/_weightedTotalNumPiPlus);
00684         scale(_histMeanMultiLambda_c_Plus     , 1.0/_weightedTotalNumPiPlus);
00685         scale(_histMeanMultiLambda_b_0        , 1.0/_weightedTotalNumPiPlus);
00686         scale(_histMeanMultiLambda1520        , 1.0/_weightedTotalNumPiPlus);
00687       }
00688 
00689       if (sqrtS()/GeV >= 130 && sqrtS()/GeV <= 200) {
00690         scale(_histMeanMultiKPlus            , 1.0/_weightedTotalNumPiPlus);
00691         scale(_histMeanMultiK0               , 1.0/_weightedTotalNumPiPlus);
00692         scale(_histMeanMultiP                , 1.0/_weightedTotalNumPiPlus);
00693         scale(_histMeanMultiLambda           , 1.0/_weightedTotalNumPiPlus);
00694       }
00695     }

const AnalysisInfo & info ( ) const [virtual, inherited]

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

Definition at line 125 of file Analysis.cc.

References Analysis::_info.

Referenced by Analysis::energies(), and Analysis::requiredBeams().

00125                                            {
00126     assert(_info.get() != 0);
00127     return *_info;
00128   }

string name ( ) const [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 130 of file Analysis.cc.

References Analysis::_defaultname, and Analysis::_info.

Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), Analysis::_makeHistoDir(), AnalysisHandler::analysisNames(), AnalysisHandler::analyze(), Analysis::crossSection(), Analysis::getLog(), Analysis::histoDir(), AnalysisHandler::init(), AnalysisBuilderBase::name(), Analysis::normalize(), HistoHandler::registerAnalysisObject(), AnalysisHandler::removeIncompatibleAnalyses(), and Analysis::scale().

00130                               {
00131     if (_info && !_info->name().empty()) return _info->name();
00132     return _defaultname;
00133   }

string spiresId ( ) const [virtual, inherited]

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

Definition at line 135 of file Analysis.cc.

References Analysis::_info.

00135                                   {
00136     if (!_info) return "NONE";
00137     return _info->spiresId();
00138   }

vector< string > authors ( ) const [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 140 of file Analysis.cc.

References Analysis::_info.

00140                                          {
00141     if (!_info) return std::vector<std::string>();
00142     return _info->authors();
00143   }

string summary ( ) const [virtual, inherited]

Get a short description of the analysis.

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

Definition at line 145 of file Analysis.cc.

References Analysis::_info.

00145                                  {
00146     if (!_info) return "NONE";
00147     return _info->summary();
00148   }

string description ( ) const [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 150 of file Analysis.cc.

References Analysis::_info.

00150                                      {
00151     if (!_info) return "NONE";
00152     return _info->description();
00153   }

string runInfo ( ) const [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 155 of file Analysis.cc.

References Analysis::_info.

00155                                  {
00156     if (!_info) return "NONE";
00157     return _info->runInfo();
00158   }

string experiment ( ) const [virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 164 of file Analysis.cc.

References Analysis::_info.

00164                                     {
00165     if (!_info) return "NONE";
00166     return _info->experiment();
00167   }

string collider ( ) const [virtual, inherited]

Collider on which the experiment ran.

Definition at line 169 of file Analysis.cc.

References Analysis::_info.

00169                                   {
00170     if (!_info) return "NONE";
00171     return _info->collider();
00172   }

const BeamPair requiredBeams ( ) const [virtual, inherited]

Return the pair of incoming beams required by this analysis.

Definition at line 189 of file Analysis.cc.

References Analysis::beams(), Analysis::info(), and Rivet::make_pdgid_pair().

Referenced by Analysis::isCompatible().

00189                                                {
00190     return make_pdgid_pair(info().beams());
00191   }

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

Sets of valid beam energy pairs, in GeV.

Definition at line 160 of file Analysis.cc.

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

00160                                                                      {
00161     return info().energies();
00162   }

string year ( ) const [virtual, inherited]

When the original experimental analysis was published.

When the refereed paper on which this is based was published, according to SPIRES.

Definition at line 174 of file Analysis.cc.

References Analysis::_info.

00174                               {
00175     if (!_info) return "NONE";
00176     return _info->year();
00177   }

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

Journal, and preprint references.

Definition at line 179 of file Analysis.cc.

References Analysis::_info.

00179                                             {
00180     if (!_info) return vector<string>();
00181     return _info->references();
00182   }

string status ( ) const [virtual, inherited]

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

Definition at line 184 of file Analysis.cc.

References Analysis::_info.

Referenced by AnalysisHandler::init().

00184                                 {
00185     if (!_info) return "UNVALIDATED";
00186     return _info->status();
00187   }

const ParticlePair & beams ( ) const [inherited]

Incoming beams for this run.

Definition at line 79 of file Analysis.cc.

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

Referenced by OPAL_1998_S3780481::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1995_S3137023::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_1996_S3486095::analyze(), Analysis::isCompatible(), and Analysis::requiredBeams().

00079                                             {
00080     return handler().beams();
00081   }

const BeamPair beamIds ( ) const [inherited]

Incoming beam IDs for this run.

Definition at line 83 of file Analysis.cc.

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

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

00083                                          {
00084     return handler().beamIds();
00085   }

double sqrtS ( ) const [inherited]

Centre of mass energy for this run.

Definition at line 75 of file Analysis.cc.

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

Referenced by UA1_1990_S2044935::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), CDF_2004_S5839831::analyze(), ALEPH_2004_S5765862::analyze(), UA1_1990_S2044935::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), CDF_2004_S5839831::finalize(), CDF_2002_S4796047::finalize(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1986_S1583476::init(), UA1_1990_S2044935::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), MC_ZJETS::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_DIPHOTON::init(), JADE_OPAL_2000_S4300807::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), and ALEPH_2004_S5765862::init().

00075                                {
00076     return handler().sqrtS();
00077   }

bool isCompatible	(	const ParticleName &	beam1,
		const ParticleName &	beam2
	)			const `[virtual, inherited]`

Is this analysis able to run on the supplied pair of beams?

Todo:: Need to also check internal consistency of the analysis' beam requirements with those of the projections it uses.

Definition at line 200 of file Analysis.cc.

References Analysis::beams(), Rivet::compatible(), and Analysis::requiredBeams().

Referenced by AnalysisHandler::removeIncompatibleAnalyses().

00200                                                                                         {
00201     BeamPair beams(beam1, beam2);
00202     return compatible(beams, requiredBeams());
00203     /// @todo Need to also check internal consistency of the analysis'
00204     /// beam requirements with those of the projections it uses.
00205   }

bool isCompatible ( const BeamPair & beams ) const [virtual, inherited]

Is this analysis able to run on the BeamPair beams ?

Todo:: Need to also check internal consistency of the analysis' beam requirements with those of the projections it uses.

Definition at line 207 of file Analysis.cc.

References Rivet::compatible(), and Analysis::requiredBeams().

00207                                                          {
00208     return compatible(beams, requiredBeams());
00209     /// @todo Need to also check internal consistency of the analysis'
00210     /// beam requirements with those of the projections it uses.
00211   }

AnalysisHandler & handler ( ) const [inherited]

Access the controlling AnalysisHandler object.

Definition at line 244 of file Analysis.cc.

References Analysis::_analysishandler.

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

00244                                            {
00245     return *_analysishandler;
00246   }

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

References Log::ERROR, Analysis::getLog(), Analysis::name(), Analysis::scale(), Log::TRACE, Analysis::tree(), and Log::WARN.

Referenced by OPAL_1998_S3780481::finalize(), H1_1994_S2919893::finalize(), ExampleAnalysis::finalize(), DELPHI_2002_069_CONF_603::finalize(), DELPHI_1995_S3137023::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(), CDF_2007_S7057202::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(), ALEPH_2004_S5765862::finalize(), and ALEPH_1996_S3486095::finalize().

00487                                                                 {
00488     if (!histo) {
00489       getLog() << Log::ERROR << "Failed to normalise histo=NULL in analysis "
00490                << name() << "(norm=" << norm << ")" << endl;
00491       return;
00492     }
00493     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00494     getLog() << Log::TRACE << "Normalizing histo " << hpath << " to " << norm << endl;
00495  
00496     double oldintg = 0.0;
00497     int nBins = histo->axis().bins();
00498     for (int iBin = 0; iBin != nBins; ++iBin) {
00499       // Leaving out factor of binWidth because AIDA's "height" already includes a width factor.
00500       oldintg += histo->binHeight(iBin); // * histo->axis().binWidth(iBin);
00501     }
00502     if (oldintg == 0.0) {
00503       getLog() << Log::WARN << "Histo " << hpath << " has null integral during normalisation" << endl;
00504       return;
00505     }
00506 
00507     // Scale by the normalisation factor.
00508     scale(histo, norm/oldintg);
00509   }

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

References Analysis::datapointsetFactory(), Log::ERROR, Analysis::getLog(), Analysis::name(), Log::TRACE, and Analysis::tree().

00512                                                              {
00513     if (!histo) {
00514       getLog() << Log::ERROR << "Failed to scale histo=NULL in analysis "
00515           << name() << "(scale=" << scale << ")" << endl;
00516       return;
00517     }
00518     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00519     getLog() << Log::TRACE << "Scaling histo " << hpath << endl;
00520  
00521     vector<double> x, y, ex, ey;
00522     for (size_t i = 0, N = histo->axis().bins(); i < N; ++i) {
00523       x.push_back(0.5 * (histo->axis().binLowerEdge(i) + histo->axis().binUpperEdge(i)));
00524       ex.push_back(histo->axis().binWidth(i)*0.5);
00525 
00526       // "Bin height" is a misnomer in the AIDA spec: width is neglected.
00527       // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
00528       y.push_back(histo->binHeight(i)*scale/histo->axis().binWidth(i));
00529 
00530       // "Bin error" is a misnomer in the AIDA spec: width is neglected.
00531       // We'd like to do this: ey.push_back(histo->binError(i) * scale);
00532       ey.push_back(histo->binError(i)*scale/(0.5*histo->axis().binWidth(i)));
00533     }
00534  
00535     string title = histo->title();
00536     string xtitle = histo->xtitle();
00537     string ytitle = histo->ytitle();
00538 
00539     tree().mkdir("/tmpnormalize");
00540     tree().mv(hpath, "/tmpnormalize");
00541  
00542     AIDA::IDataPointSet* dps = datapointsetFactory().createXY(hpath, title, x, y, ex, ey);
00543     dps->setXTitle(xtitle);
00544     dps->setYTitle(ytitle);
00545  
00546     tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
00547     tree().rmdir("/tmpnormalize");
00548  
00549     // Set histo pointer to null - it can no longer be used.
00550     histo = 0;
00551   }

Analysis & setCrossSection ( double xs ) [inherited]

Set the cross section from the generator.

Definition at line 214 of file Analysis.cc.

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

Referenced by AnalysisHandler::setCrossSection().

00214                                                {
00215     _crossSection = xs;
00216     _gotCrossSection = true;
00217     return *this;
00218   }

bool needsCrossSection ( ) const [inherited]

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

Definition at line 220 of file Analysis.cc.

References Analysis::_needsCrossSection.

Referenced by AnalysisHandler::needCrossSection().

00220                                          {
00221     return _needsCrossSection;
00222   }

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

References Log::getLog(), and Analysis::name().

Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), 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(), UA1_1990_S2044935::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(), PDG_HADRON_MULTIPLICITIES::analyze(), OPAL_1998_S3780481::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_LEADINGJETS::analyze(), JADE_OPAL_2000_S4300807::analyze(), H1_1994_S2919893::analyze(), ExampleAnalysis::analyze(), DELPHI_2002_069_CONF_603::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_2008_S6879055::analyze(), D0_2007_S7075677::analyze(), D0_2006_S6438750::analyze(), D0_2004_S5992206::analyze(), D0_2001_S4674421::analyze(), D0_1998_S3711838::analyze(), CDF_2009_S8383952::analyze(), CDF_2009_S8233977::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_2005_S6217184::analyze(), CDF_2004_S5839831::analyze(), CDF_2001_S4751469::analyze(), CDF_2000_S4155203::analyze(), CDF_1994_S2952106::analyze(), CDF_1991_S2313472::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_1991_S2435284::analyze(), CDF_2009_S8057893::CDF_2009_S8057893::analyze(), UA5_1986_S1583476::finalize(), UA1_1990_S2044935::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), D0_2001_S4674421::finalize(), CDF_2009_S8233977::finalize(), CDF_2007_S7057202::finalize(), CDF_2006_S6653332::finalize(), ZEUS_2001_S4815815::init(), JADE_OPAL_2000_S4300807::init(), ALEPH_2004_S5765862::init(), Analysis::normalize(), and Analysis::scale().

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

double crossSection ( ) const [protected, inherited]

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

Definition at line 229 of file Analysis.cc.

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

00229                                       {
00230     if (!_gotCrossSection || _crossSection < 0) {
00231       string errMsg = "You did not set the cross section for the analysis " + name();
00232       throw Error(errMsg);
00233     }
00234     return _crossSection;
00235   }

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

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

Referenced by UA1_1990_S2044935::finalize(), MC_PHOTONJETS::finalize(), D0_2001_S4674421::finalize(), CDF_2007_S7057202::finalize(), and CDF_1988_S1865951::finalize().

00237                                               {
00238     const double sumW = sumOfWeights();
00239     assert(sumW > 0);
00240     return _crossSection / sumW;
00241   }

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

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

00112                                    {
00113     return handler().numEvents();
00114   }

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 117 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(), MC_ZJETS::finalize(), MC_WJETS::finalize(), MC_TTBAR::finalize(), MC_JetAnalysis::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(), 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_1998_S3618439::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), and ALEPH_1991_S2435284::finalize().

00117                                       {
00118     return handler().sumOfWeights();
00119   }

IAnalysisFactory & analysisFactory ( ) [protected, inherited]

Access the AIDA analysis factory of the controlling AnalysisHandler object.

Definition at line 55 of file Analysis.cc.

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

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

ITree & tree ( ) [protected, inherited]

Access the AIDA tree of the controlling AnalysisHandler object.

Definition at line 60 of file Analysis.cc.

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

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

00060                         {
00061     return handler().tree();
00062   }

IHistogramFactory & histogramFactory ( ) [protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 65 of file Analysis.cc.

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

Referenced by STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), H1_1995_S3167097::finalize(), D0_2008_S7837160::finalize(), D0_2008_S7719523::finalize(), and D0_2001_S4674421::finalize().

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

IDataPointSetFactory & datapointsetFactory ( ) [protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 70 of file Analysis.cc.

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

Referenced by Analysis::scale().

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

const string histoDir ( ) const [protected, inherited]

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

Definition at line 88 of file Analysis.cc.

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

Referenced by Analysis::_makeHistoDir(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), D0_2008_S7837160::finalize(), D0_2008_S7719523::finalize(), and D0_2001_S4674421::finalize().

00088                                         {
00089     string path = "/" + name();
00090     if (handler().runName().length() > 0) {
00091       path = "/" + handler().runName() + path;
00092     }
00093     while (find_first(path, "//")) {
00094       replace_all(path, "//", "/");
00095     }
00096     return path;
00097   }

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

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

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

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

Referenced by Analysis::binEdges(), and D0_2008_S7837160::init().

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

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

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

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

Get bin edges with logarithmic widths.

Definition at line 292 of file Analysis.cc.

Referenced by MC_ZJETS::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), and MC_DIPHOTON::init().

00292                                                                          {
00293     assert(lower>0.0);
00294     assert(upper>lower);
00295     double loglower=log10(lower);
00296     double logupper=log10(upper);
00297     vector<double> binedges;
00298     double stepwidth=(logupper-loglower)/double(nbins);
00299     for (size_t i=0; i<=nbins; ++i) {
00300       binedges.push_back(pow(10.0, loglower+double(i)*stepwidth));
00301     }
00302     return binedges;
00303   }

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 = `""`
	)			`[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!)

Referenced by ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::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_1998_S3780481::init(), MC_ZJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_DIPHOTON::init(), MC_DIJET::init(), JADE_OPAL_2000_S4300807::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1995_S3137023::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_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::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_1991_S2313472::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2010_S8591806::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), and ALEPH_1991_S2435284::init().

AIDA::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!)

AIDA::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.

AIDA::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.

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 = `""`
	)			`[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!)

Referenced by UA5_1988_S1867512::init(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2006_S6860818::init(), MC_PHOTONJETUE::init(), MC_LEADINGJETS::init(), H1_2000_S4129130::init(), H1_1994_S2919893::init(), DELPHI_2002_069_CONF_603::init(), D0_1996_S3324664::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(), and ATLAS_2010_S8591806::init().

AIDA::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!)

AIDA::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.

AIDA::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.

AIDA::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!)

Referenced by STAR_2006_S6860818::init(), MC_JetAnalysis::init(), JADE_OPAL_2000_S4300807::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_1996_S3418421::init(), and ALEPH_2004_S5765862::init().

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 = `""`
	)			`[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!)

AIDA::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.

Analysis & setBeams	(	const ParticleName &	beam1,
		const ParticleName &	beam2
	)			`[protected, inherited]`

Set the colliding beam pair.

Deprecated:: Use .info file and AnalysisInfo class instead

Definition at line 193 of file Analysis.cc.

References Analysis::_info.

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(), BELLE_2006_S6265367::BELLE_2006_S6265367(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1990_S2089246::CDF_1990_S2089246(), CDF_1991_S2313472::CDF_1991_S2313472(), 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_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(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_1998_S3711838::D0_1998_S3711838(), 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(), DELPHI_1995_S3137023::DELPHI_1995_S3137023(), 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_OPAL_2000_S4300807::JADE_OPAL_2000_S4300807(), OPAL_1998_S3780481::OPAL_1998_S3780481(), 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_S7993412::STAR_2008_S7993412(), STAR_2009_UE_HELEN::STAR_2009_UE_HELEN(), UA1_1990_S2044935::UA1_1990_S2044935(), UA5_1986_S1583476::UA5_1986_S1583476(), UA5_1988_S1867512::UA5_1988_S1867512(), UA5_1989_S1926373::UA5_1989_S1926373(), and ZEUS_2001_S4815815::ZEUS_2001_S4815815().

00193                                                                                    {
00194     assert(_info.get() != 0);
00195     _info->_beams = make_pair(beam1, beam2);
00196     return *this;
00197   }

Analysis & setNeedsCrossSection ( bool needed ) [protected, inherited]

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

Definition at line 224 of file Analysis.cc.

References Analysis::_needsCrossSection.

Referenced by ATLAS_2010_S8591806::ATLAS_2010_S8591806(), CDF_1991_S2313472::CDF_1991_S2313472(), 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_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_1998_S3711838::D0_1998_S3711838(), 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_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), MC_DIPHOTON::MC_DIPHOTON(), MC_JetAnalysis::MC_JetAnalysis(), MC_JETS::MC_JETS(), MC_PHOTONJETS::MC_PHOTONJETS(), MC_WJETS::MC_WJETS(), MC_ZJETS::MC_ZJETS(), STAR_2006_S6870392::STAR_2006_S6870392(), and UA1_1990_S2044935::UA1_1990_S2044935().

00224                                                       {
00225     _needsCrossSection = needed;
00226     return *this;
00227   }

std::set<ConstProjectionPtr> getProjections ( ) const [inline, inherited]

Get the contained projections, including recursion.

Definition at line 43 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

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

const PROJ& getProjection ( const std::string & name ) const [inline, inherited]

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

Definition at line 50 of file ProjectionApplier.hh.

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

Referenced by VetoedFinalState::compare(), Rivet::pcmp(), and Hemispheres::project().

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

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 58 of file ProjectionApplier.hh.

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

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

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

Apply the supplied projection on event.

Definition at line 68 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

Referenced by HadronicFinalState::project(), and FinalStateHCM::project().

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

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

Apply the supplied projection on event.

Definition at line 75 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

Apply the named projection on event.

Definition at line 82 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

ProjectionHandler& getProjHandler ( ) const [inline, protected, inherited]

Get a reference to the ProjectionHandler for this thread.

Definition at line 95 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

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

00095                                               {
00096       assert(_projhandler);
00097       return *_projhandler;
00098     }

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 115 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

Referenced by ZFinder::_init(), WFinder::_init(), VetoedFinalState::addVetoOnThisFinalState(), 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(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::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_1998_S3780481::init(), MC_ZJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_LEADINGJETS::init(), MC_JETS::init(), MC_DIPHOTON::init(), MC_DIJET::init(), JADE_OPAL_2000_S4300807::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1995_S3137023::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_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::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_1991_S2313472::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2010_S8591806::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), ALEPH_1991_S2435284::init(), IsolationProjection::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), KtJets::KtJets(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState::LossyFinalState(), MergedFinalState::MergedFinalState(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerUA5::TriggerUA5(), and VetoedFinalState::VetoedFinalState().

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

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, Log::ERROR, ProjectionApplier::getLog(), ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().

Referenced by ProjectionApplier::addProjection().

00034                                                                              {
00035     if (!_allowProjReg) {
00036       getLog() << Log::ERROR << "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   }