ALEPH_1996_S3486095 Class Reference

Inheritance diagram for ALEPH_1996_S3486095:

[legend]

Collaboration diagram for ALEPH_1996_S3486095:

[legend]

List of all members.

Detailed Description

ALEPH QCD study with event shapes and identified particles.

Author:: Holger Schulz

Definition at line 20 of file ALEPH_1996_S3486095.cc.


Public Member Functions
	ALEPH_1996_S3486095 ()
	Constructor.
virtual bool	isCompatible (PdgId beam1, PdgId beam2) const
	Is this analysis able to run on the supplied pair of beams?
virtual bool	isCompatible (const PdgIdPair &beams) const
	Is this analysis able to run on the PdgIdPair 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	init ()
void	analyze (const Event &e)
void	finalize ()
	Finalize.
Metadata
Metadata is used for querying from the command line and also for building web pages and the analysis pages in the Rivet manual.
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 std::vector < PdgIdPair > &	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	bibKey () const
	BibTeX citation key for this article.
virtual std::string	bibTeX () const
	BibTeX citation entry for this article.
virtual std::string	status () const
	Whether this analysis is trusted (in any way!).
virtual std::vector< std::string >	todos () const
	Any work to be done on this analysis.
Run conditions
const ParticlePair &	beams () const
	Incoming beams for this run.
const PdgIdPair	beamIds () const
	Incoming beam IDs for this run.
double	sqrtS () const
	Centre of mass energy for this run.
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 (PdgId beam1, PdgId 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	_weightedTotalPartNum
double	_weightedTotalNumPiPlus
double	_weightedTotalNumKPlus
double	_weightedTotalNumP
double	_weightedTotalNumPhoton
double	_weightedTotalNumPi0
double	_weightedTotalNumEta
double	_weightedTotalNumEtaPrime
double	_weightedTotalNumK0
double	_weightedTotalNumLambda0
double	_weightedTotalNumXiMinus
double	_weightedTotalNumSigma1385Plus
double	_weightedTotalNumXi1530_0
double	_weightedTotalNumRho
double	_weightedTotalNumOmega782
double	_weightedTotalNumKStar892_0
double	_weightedTotalNumPhi
double	_weightedTotalNumKStar892Plus
double	_numChParticles
Histograms
AIDA::IHistogram1D *	_histSphericity
AIDA::IHistogram1D *	_histAplanarity
AIDA::IHistogram1D *	_hist1MinusT
AIDA::IHistogram1D *	_histTMinor
AIDA::IHistogram1D *	_histY3
AIDA::IHistogram1D *	_histHeavyJetMass
AIDA::IHistogram1D *	_histCParam
AIDA::IHistogram1D *	_histOblateness
AIDA::IHistogram1D *	_histScaledMom
AIDA::IHistogram1D *	_histRapidityT
AIDA::IHistogram1D *	_histPtSIn
AIDA::IHistogram1D *	_histPtSOut
AIDA::IHistogram1D *	_histJetRate2Durham
AIDA::IHistogram1D *	_histJetRate3Durham
AIDA::IHistogram1D *	_histJetRate4Durham
AIDA::IHistogram1D *	_histJetRate5Durham
AIDA::IHistogram1D *	_histLogScaledMom
AIDA::IHistogram1D *	_histChMult
AIDA::IHistogram1D *	_histMultiPiPlus
AIDA::IHistogram1D *	_histMultiKPlus
AIDA::IHistogram1D *	_histMultiP
AIDA::IHistogram1D *	_histMultiPhoton
AIDA::IHistogram1D *	_histMultiPi0
AIDA::IHistogram1D *	_histMultiEta
AIDA::IHistogram1D *	_histMultiEtaPrime
AIDA::IHistogram1D *	_histMultiK0
AIDA::IHistogram1D *	_histMultiLambda0
AIDA::IHistogram1D *	_histMultiXiMinus
AIDA::IHistogram1D *	_histMultiSigma1385Plus
AIDA::IHistogram1D *	_histMultiXi1530_0
AIDA::IHistogram1D *	_histMultiRho
AIDA::IHistogram1D *	_histMultiOmega782
AIDA::IHistogram1D *	_histMultiKStar892_0
AIDA::IHistogram1D *	_histMultiPhi
AIDA::IHistogram1D *	_histMultiKStar892Plus
AIDA::IHistogram1D *	_histMeanChMult
AIDA::IHistogram1D *	_histMeanChMultRapt05
AIDA::IHistogram1D *	_histMeanChMultRapt10
AIDA::IHistogram1D *	_histMeanChMultRapt15
AIDA::IHistogram1D *	_histMeanChMultRapt20
AIDA::IHistogram1D *	_histMeanMultiPi0
AIDA::IHistogram1D *	_histMeanMultiEta
AIDA::IHistogram1D *	_histMeanMultiEtaPrime
AIDA::IHistogram1D *	_histMeanMultiK0
AIDA::IHistogram1D *	_histMeanMultiRho
AIDA::IHistogram1D *	_histMeanMultiOmega782
AIDA::IHistogram1D *	_histMeanMultiPhi
AIDA::IHistogram1D *	_histMeanMultiKStar892Plus
AIDA::IHistogram1D *	_histMeanMultiKStar892_0
AIDA::IHistogram1D *	_histMeanMultiLambda0
AIDA::IHistogram1D *	_histMeanMultiSigma0
AIDA::IHistogram1D *	_histMeanMultiXiMinus
AIDA::IHistogram1D *	_histMeanMultiSigma1385Plus
AIDA::IHistogram1D *	_histMeanMultiXi1530_0
AIDA::IHistogram1D *	_histMeanMultiOmegaOmegaBar
Friends
class	Projectionhandler

Constructor & Destructor Documentation

ALEPH_1996_S3486095 ( ) [inline]

Constructor.

Definition at line 25 of file ALEPH_1996_S3486095.cc.

00026       : Analysis("ALEPH_1996_S3486095")
00027     {
00028       setBeams(ELECTRON, POSITRON);
00029       _numChParticles               = 0;
00030       _weightedTotalPartNum         = 0;
00031       _weightedTotalNumPiPlus       = 0;
00032       _weightedTotalNumKPlus        = 0;
00033       _weightedTotalNumP            = 0;
00034       _weightedTotalNumPhoton       = 0;
00035       _weightedTotalNumPi0          = 0;
00036       _weightedTotalNumEta          = 0;
00037       _weightedTotalNumEtaPrime     = 0;
00038       _weightedTotalNumK0           = 0;
00039       _weightedTotalNumLambda0      = 0;
00040       _weightedTotalNumXiMinus      = 0;
00041       _weightedTotalNumSigma1385Plus= 0;
00042       _weightedTotalNumXi1530_0     = 0;
00043       _weightedTotalNumRho          = 0;
00044       _weightedTotalNumOmega782     = 0;
00045       _weightedTotalNumKStar892_0   = 0;
00046       _weightedTotalNumPhi          = 0;
00047       _weightedTotalNumKStar892Plus = 0;
00048     }

Member Function Documentation

void init ( ) [inline, virtual]

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

Implements Analysis.

Definition at line 54 of file ALEPH_1996_S3486095.cc.

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

void analyze ( const Event & event ) [inline, virtual]

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

Implements Analysis.

Definition at line 135 of file ALEPH_1996_S3486095.cc.

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

void finalize ( ) [inline, virtual]

Finalize.

Implements Analysis.

Definition at line 359 of file ALEPH_1996_S3486095.cc.

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

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

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

Definition at line 126 of file Analysis.cc.

References Analysis::_info.

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

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

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

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

Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), Analysis::_makeHistoDir(), Analysis::crossSection(), Analysis::getLog(), Analysis::histoDir(), Analysis::normalize(), HistoHandler::registerAnalysisObject(), and Analysis::scale().

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

string spiresId ( ) const [virtual, inherited]

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

Definition at line 136 of file Analysis.cc.

References Analysis::_info.

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

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

References Analysis::_info.

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

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

References Analysis::_info.

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

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

References Analysis::_info.

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

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

References Analysis::_info.

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

string experiment ( ) const [virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 165 of file Analysis.cc.

References Analysis::_info.

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

string collider ( ) const [virtual, inherited]

Collider on which the experiment ran.

Definition at line 170 of file Analysis.cc.

References Analysis::_info.

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

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

Return the pair of incoming beams required by this analysis.

Definition at line 205 of file Analysis.cc.

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

Referenced by Analysis::isCompatible().

00205                                                          {
00206     return info().beams();
00207   }

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

Sets of valid beam energy pairs, in GeV.

Definition at line 161 of file Analysis.cc.

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

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

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

References Analysis::_info.

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

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

Journal, and preprint references.

Definition at line 180 of file Analysis.cc.

References Analysis::_info.

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

string bibKey ( ) const [virtual, inherited]

BibTeX citation key for this article.

Definition at line 185 of file Analysis.cc.

References Analysis::_info.

00185                                 {
00186     if (!_info) return "";
00187     return _info->bibKey();
00188   }

string bibTeX ( ) const [virtual, inherited]

BibTeX citation entry for this article.

Definition at line 190 of file Analysis.cc.

References Analysis::_info.

00190                                 {
00191     if (!_info) return "";
00192     return _info->bibTeX();
00193   }

string status ( ) const [virtual, inherited]

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

Definition at line 195 of file Analysis.cc.

References Analysis::_info.

00195                                 {
00196     if (!_info) return "UNVALIDATED";
00197     return _info->status();
00198   }

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

Any work to be done on this analysis.

Definition at line 200 of file Analysis.cc.

References Analysis::_info.

00200                                        {
00201     if (!_info) return vector<string>();
00202     return _info->todos();
00203   }

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_2004_S5765862::analyze(), ALEPH_1996_S3486095::analyze(), and Analysis::isCompatible().

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

const PdgIdPair 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(), ALEPH_2004_S5765862::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_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::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	(	PdgId	beam1,
		PdgId	beam2
	)			const `[virtual, inherited]`

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

Todo:: Deprecate?

Definition at line 220 of file Analysis.cc.

References Analysis::beams().

00220                                                             {
00221     PdgIdPair beams(beam1, beam2);
00222     return isCompatible(beams);
00223   }

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

Is this analysis able to run on the PdgIdPair beams ?

Todo:: Deprecate?

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

Definition at line 227 of file Analysis.cc.

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

00227                                                           {
00228     foreach (const PdgIdPair& bp, requiredBeams()) {
00229       if (compatible(beams, bp)) return true;
00230     }
00231     return false;
00232     /// @todo Need to also check internal consistency of the analysis'
00233     /// beam requirements with those of the projections it uses.
00234   }

AnalysisHandler & handler ( ) const [inherited]

Access the controlling AnalysisHandler object.

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

00269                                            {
00270     return *_analysishandler;
00271   }

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

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

Referenced by STAR_2008_S7869363::finalize(), OPAL_2001_S4553896::finalize(), OPAL_1998_S3780481::finalize(), H1_1994_S2919893::finalize(), ExampleAnalysis::finalize(), DELPHI_2002_069_CONF_603::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8671338::finalize(), D0_2009_S8320160::finalize(), D0_2008_S7554427::finalize(), D0_2008_S6879055::finalize(), D0_2007_S7075677::finalize(), D0_2004_S5992206::finalize(), D0_2001_S4674421::finalize(), D0_1996_S3324664::finalize(), D0_1996_S3214044::finalize(), CDF_2009_NOTE_9936::finalize(), CDF_2004_S5839831::finalize(), CDF_2002_S4796047::finalize(), CDF_2001_S4751469::finalize(), CDF_1997_S3541940::finalize(), CDF_1996_S3418421::finalize(), CDF_1996_S3349578::finalize(), CDF_1996_S3108457::finalize(), CDF_1993_S2742446::finalize(), ALEPH_2004_S5765862::finalize(), and ALEPH_1996_S3486095::finalize().

00512                                                                 {
00513     if (!histo) {
00514       getLog() << Log::ERROR << "Failed to normalise histo=NULL in analysis "
00515                << name() << " (norm=" << norm << ")" << endl;
00516       return;
00517     }
00518     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00519     getLog() << Log::TRACE << "Normalizing histo " << hpath << " to " << norm << endl;
00520 
00521     double oldintg = 0.0;
00522     int nBins = histo->axis().bins();
00523     for (int iBin = 0; iBin != nBins; ++iBin) {
00524       // Leaving out factor of binWidth because AIDA's "height" already includes a width factor.
00525       oldintg += histo->binHeight(iBin); // * histo->axis().binWidth(iBin);
00526     }
00527     if (oldintg == 0.0) {
00528       getLog() << Log::WARN << "Histo " << hpath << " has null integral during normalisation" << endl;
00529       return;
00530     }
00531 
00532     // Scale by the normalisation factor.
00533     scale(histo, norm/oldintg);
00534   }

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

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

Referenced by UA5_1989_S1926373::finalize(), UA5_1987_S1640666::finalize(), UA5_1986_S1583476::finalize(), UA5_1982_S875503::finalize(), UA1_1990_S2044935::finalize(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), SFM_1984_S1178091::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_2004_S6132243::finalize(), OPAL_1998_S3780481::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_TTBAR::finalize(), MC_PHOTONJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), JADE_OPAL_2000_S4300807::finalize(), H1_2000_S4129130::finalize(), H1_1994_S2919893::finalize(), E735_1998_S3905616::finalize(), D0_2010_S8570965::finalize(), D0_2009_S8349509::finalize(), D0_2009_S8202443::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7540469::finalize(), CDF_2007_S7057202::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), CDF_1990_S2089246::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2010_S8591806::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), ALEPH_1991_S2435284::finalize(), Analysis::normalize(), and BinnedHistogram::scale().

00537                                                              {
00538     if (!histo) {
00539       getLog() << Log::ERROR << "Failed to scale histo=NULL in analysis "
00540           << name() << " (scale=" << scale << ")" << endl;
00541       return;
00542     }
00543     const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo));
00544     getLog() << Log::TRACE << "Scaling histo " << hpath << endl;
00545 
00546     vector<double> x, y, ex, ey;
00547     for (size_t i = 0, N = histo->axis().bins(); i < N; ++i) {
00548       x.push_back(0.5 * (histo->axis().binLowerEdge(i) + histo->axis().binUpperEdge(i)));
00549       ex.push_back(histo->axis().binWidth(i)*0.5);
00550 
00551       // "Bin height" is a misnomer in the AIDA spec: width is neglected.
00552       // We'd like to do this: y.push_back(histo->binHeight(i) * scale);
00553       y.push_back(histo->binHeight(i)*scale/histo->axis().binWidth(i));
00554 
00555       // "Bin error" is a misnomer in the AIDA spec: width is neglected.
00556       // We'd like to do this: ey.push_back(histo->binError(i) * scale);
00557       ey.push_back(histo->binError(i)*scale/histo->axis().binWidth(i));
00558     }
00559 
00560     string title = histo->title();
00561     string xtitle = histo->xtitle();
00562     string ytitle = histo->ytitle();
00563 
00564     tree().mkdir("/tmpnormalize");
00565     tree().mv(hpath, "/tmpnormalize");
00566 
00567     AIDA::IDataPointSet* dps = datapointsetFactory().createXY(hpath, title, x, y, ex, ey);
00568     dps->setXTitle(xtitle);
00569     dps->setYTitle(ytitle);
00570 
00571     tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo)));
00572     tree().rmdir("/tmpnormalize");
00573 
00574     // Set histo pointer to null - it can no longer be used.
00575     histo = 0;
00576   }

Analysis & setCrossSection ( double xs ) [inherited]

Set the cross section from the generator.

Definition at line 237 of file Analysis.cc.

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

00237                                                {
00238     _crossSection = xs;
00239     _gotCrossSection = true;
00240     return *this;
00241   }

bool needsCrossSection ( ) const [inherited]

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

Todo:: Deprecate, eventually

Definition at line 244 of file Analysis.cc.

References Analysis::_needsCrossSection.

00244                                          {
00245     return _needsCrossSection;
00246   }

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

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

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

double crossSection ( ) const [protected, inherited]

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

Definition at line 254 of file Analysis.cc.

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

00254                                       {
00255     if (!_gotCrossSection || _crossSection < 0) {
00256       string errMsg = "You did not set the cross section for the analysis " + name();
00257       throw Error(errMsg);
00258     }
00259     return _crossSection;
00260   }

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

00262                                               {
00263     const double sumW = sumOfWeights();
00264     assert(sumW > 0);
00265     return _crossSection / sumW;
00266   }

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

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

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

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

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

Referenced by Analysis::crossSectionPerEvent(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_1998_S3780481::finalize(), OPAL_1993_S2692198::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_TTBAR::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), JADE_OPAL_2000_S4300807::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7541902::finalize(), CDF_2008_S7540469::finalize(), CDF_2006_S6653332::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2001_S4517016::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), and ALEPH_1991_S2435284::finalize().

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

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(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_GENERIC::finalize(), H1_1995_S3167097::finalize(), D0_2008_S7837160::finalize(), and D0_2008_S7719523::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.

Todo:: This doesn't change: calc and cache at Analysis construction!

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(), MC_WJETS::finalize(), D0_2008_S7837160::finalize(), and D0_2008_S7719523::finalize().

00088                                         {
00089     /// @todo This doesn't change: calc and cache at Analysis construction!
00090     string path = "/" + name();
00091     if (handler().runName().length() > 0) {
00092       path = "/" + handler().runName() + path;
00093     }
00094     while (find_first(path, "//")) {
00095       replace_all(path, "//", "/");
00096     }
00097     return path;
00098   }

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

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

Referenced by MC_JetAnalysis::finalize(), and MC_GENERIC::finalize().

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

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

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

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

Get bin edges with logarithmic widths.

Definition at line 317 of file Analysis.cc.

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

00317                                                                          {
00318     assert(lower>0.0);
00319     assert(upper>lower);
00320     double loglower=log10(lower);
00321     double logupper=log10(upper);
00322     vector<double> binedges;
00323     double stepwidth=(logupper-loglower)/double(nbins);
00324     for (size_t i=0; i<=nbins; ++i) {
00325       binedges.push_back(pow(10.0, loglower+double(i)*stepwidth));
00326     }
00327     return binedges;
00328   }

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_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), 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_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_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 ALEPH_2004_S5765862::finalize(), STAR_2006_S6860818::init(), OPAL_1993_S2692198::init(), MC_JetAnalysis::init(), JADE_OPAL_2000_S4300807::init(), D0_2001_S4674421::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_1996_S3418421::init(), CDF_1994_S2952106::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	(	PdgId	beam1,
		PdgId	beam2
	)			`[protected, inherited]`

Set the colliding beam pair.

Deprecated:: Use .info file and AnalysisInfo class instead

Todo:: Deprecate?

Definition at line 211 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_1993_S2742446::CDF_1993_S2742446(), CDF_1994_S2952106::CDF_1994_S2952106(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1996_S3349578::CDF_1996_S3349578(), CDF_1996_S3418421::CDF_1996_S3418421(), CDF_1997_S3541940::CDF_1997_S3541940(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2001_S4751469::CDF_2001_S4751469(), CDF_2002_S4796047::CDF_2002_S4796047(), CDF_2004_S5839831::CDF_2004_S5839831(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2005_S6217184::CDF_2005_S6217184(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_LEADINGJETS::CDF_2008_LEADINGJETS(), CDF_2008_NOTE_9351::CDF_2008_NOTE_9351(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7782535::CDF_2008_S7782535(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2008_S8095620::CDF_2008_S8095620(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8057893::CDF_2009_S8057893::CDF_2009_S8057893(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), 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(), D0_2010_S8671338::D0_2010_S8671338(), 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_1993_S2692198::OPAL_1993_S2692198(), OPAL_1998_S3780481::OPAL_1998_S3780481(), OPAL_2001_S4553896::OPAL_2001_S4553896(), PDG_HADRON_MULTIPLICITIES::PDG_HADRON_MULTIPLICITIES(), PDG_HADRON_MULTIPLICITIES_RATIOS::PDG_HADRON_MULTIPLICITIES_RATIOS(), SFM_1984_S1178091::SFM_1984_S1178091(), STAR_2006_S6500200::STAR_2006_S6500200(), STAR_2006_S6860818::STAR_2006_S6860818(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), STAR_2008_S7993412::STAR_2008_S7993412(), STAR_2009_UE_HELEN::STAR_2009_UE_HELEN(), UA1_1990_S2044935::UA1_1990_S2044935(), UA5_1986_S1583476::UA5_1986_S1583476(), UA5_1987_S1640666::UA5_1987_S1640666(), UA5_1988_S1867512::UA5_1988_S1867512(), UA5_1989_S1926373::UA5_1989_S1926373(), and ZEUS_2001_S4815815::ZEUS_2001_S4815815().

00211                                                        {
00212     assert(_info.get() != 0);
00213     _info->_beams.clear();
00214     _info->_beams += make_pair(beam1, beam2);
00215     return *this;
00216   }

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

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

Todo:: Deprecate, eventually

Definition at line 249 of file Analysis.cc.

References Analysis::_needsCrossSection.

Referenced by ATLAS_2010_S8591806::ATLAS_2010_S8591806(), CDF_1991_S2313472::CDF_1991_S2313472(), CDF_1993_S2742446::CDF_1993_S2742446(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_1998_S3711838::D0_1998_S3711838(), D0_2000_S4480767::D0_2000_S4480767(), D0_2001_S4674421::D0_2001_S4674421(), D0_2006_S6438750::D0_2006_S6438750(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7863608::D0_2008_S7863608(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), MC_DIPHOTON::MC_DIPHOTON(), MC_HJETS::MC_HJETS(), MC_JetAnalysis::MC_JetAnalysis(), MC_JETS::MC_JETS(), MC_PHOTONJETS::MC_PHOTONJETS(), MC_WJETS::MC_WJETS(), MC_WWJETS::MC_WWJETS(), MC_ZJETS::MC_ZJETS(), MC_ZZJETS::MC_ZZJETS(), OPAL_2001_S4553896::OPAL_2001_S4553896(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), and UA1_1990_S2044935::UA1_1990_S2044935().

00249                                                       {
00250     _needsCrossSection = needed;
00251     return *this;
00252   }

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

Get the contained projections, including recursion.

Definition at line 58 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

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

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

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

Definition at line 65 of file ProjectionApplier.hh.

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

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

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

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

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

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

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

Apply the supplied projection on event.

Definition at line 83 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

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

Apply the supplied projection on event.

Definition at line 90 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

Apply the named projection on event.

Definition at line 97 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

Get a reference to the ProjectionHandler for this thread.

Definition at line 110 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

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

00110                                               {
00111       assert(_projhandler);
00112       return *_projhandler;
00113     }

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 130 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_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), OPAL_1993_S2692198::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_LEADINGJETS::init(), MC_JETS::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), 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_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6217184::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_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(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< Rivet::ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), and VetoedFinalState::VetoedFinalState().

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

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   }

Friends And Related Function Documentation

friend class Projectionhandler [friend, inherited]

Definition at line 38 of file ProjectionApplier.hh.

Member Data Documentation

double _weightedTotalPartNum [private]

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

Definition at line 456 of file ALEPH_1996_S3486095.cc.

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

double _weightedTotalNumPiPlus [private]

Definition at line 457 of file ALEPH_1996_S3486095.cc.

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

double _weightedTotalNumKPlus [private]

Definition at line 458 of file ALEPH_1996_S3486095.cc.

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

double _weightedTotalNumP [private]

Definition at line 459 of file ALEPH_1996_S3486095.cc.

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

double _weightedTotalNumPhoton [private]