CDF_1996_S3349578 Class Reference

Inheritance diagram for CDF_1996_S3349578:

[legend]

Collaboration diagram for CDF_1996_S3349578:

[legend]

List of all members.

Detailed Description

CDF properties of high-mass multi-jet events.

Definition at line 12 of file CDF_1996_S3349578.cc.


Public Member Functions
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.
Constructors etc.
	CDF_1996_S3349578 ()
	Constructor.
Analysis methods
void	init ()
	Book histograms and initialise projections before the run.
void	analyze (const Event &event)
void	_threeJetAnalysis (const Jets &jets, const double &weight)
void	_fourJetAnalysis (const Jets &jets, const double &weight)
void	_fiveJetAnalysis (const Jets &jets, const double &weight)
void	finalize ()
	Normalise histograms etc., after the run.
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 Member Functions
vector< FourMomentum >	_reduce (const vector< FourMomentum > &jets, FourMomentum &combined1, FourMomentum &combined2)
FourMomentum	_avg_beam_in_lab (const double &m, const double &y)
double	_psi (const FourMomentum &p1, const FourMomentum &p2, const FourMomentum &p3, const FourMomentum &p4)
double	_safeMass (const FourMomentum &p)
Private Attributes
Histograms
AIDA::IHistogram1D *	_h_3_mNJ
AIDA::IHistogram1D *	_h_3_X3
AIDA::IHistogram1D *	_h_3_X4
AIDA::IHistogram1D *	_h_3_costheta3
AIDA::IHistogram1D *	_h_3_psi3
AIDA::IHistogram1D *	_h_3_f3
AIDA::IHistogram1D *	_h_3_f4
AIDA::IHistogram1D *	_h_3_f5
AIDA::IHistogram1D *	_h_4_mNJ
AIDA::IHistogram1D *	_h_4_X3
AIDA::IHistogram1D *	_h_4_X4
AIDA::IHistogram1D *	_h_4_costheta3
AIDA::IHistogram1D *	_h_4_psi3
AIDA::IHistogram1D *	_h_4_f3
AIDA::IHistogram1D *	_h_4_f4
AIDA::IHistogram1D *	_h_4_f5
AIDA::IHistogram1D *	_h_4_XA
AIDA::IHistogram1D *	_h_4_psiAB
AIDA::IHistogram1D *	_h_4_fA
AIDA::IHistogram1D *	_h_4_fB
AIDA::IHistogram1D *	_h_5_mNJ
AIDA::IHistogram1D *	_h_5_X3
AIDA::IHistogram1D *	_h_5_X4
AIDA::IHistogram1D *	_h_5_costheta3
AIDA::IHistogram1D *	_h_5_psi3
AIDA::IHistogram1D *	_h_5_f3
AIDA::IHistogram1D *	_h_5_f4
AIDA::IHistogram1D *	_h_5_f5
AIDA::IHistogram1D *	_h_5_XA
AIDA::IHistogram1D *	_h_5_XC
AIDA::IHistogram1D *	_h_5_psiAB
AIDA::IHistogram1D *	_h_5_psiCD
AIDA::IHistogram1D *	_h_5_fA
AIDA::IHistogram1D *	_h_5_fB
AIDA::IHistogram1D *	_h_5_fC
AIDA::IHistogram1D *	_h_5_fD
Friends
class	Projectionhandler

Constructor & Destructor Documentation

CDF_1996_S3349578 ( ) [inline]

Constructor.

Definition at line 19 of file CDF_1996_S3349578.cc.

References Rivet::ANTIPROTON, Rivet::PROTON, and Analysis::setBeams().

00020       : Analysis("CDF_1996_S3349578")
00021     {
00022       setBeams(PROTON, ANTIPROTON);
00023     }

Member Function Documentation

void init ( ) [inline, virtual]

Book histograms and initialise projections before the run.

Initialise and register projections here

Book histograms here, e.g.:

Implements Analysis.

Definition at line 34 of file CDF_1996_S3349578.cc.

00034                 {
00035 
00036       /// Initialise and register projections here
00037       const FinalState fs(-4.2, 4.2);
00038       addProjection(FastJets(fs, FastJets::CDFJETCLU, 0.7), "Jets");
00039 
00040       /// Book histograms here, e.g.:
00041       _h_3_mNJ = bookHistogram1D(1, 1, 1);
00042       _h_3_X3 = bookHistogram1D(2, 1, 1);
00043       _h_3_X4 = bookHistogram1D(3, 1, 1);
00044       _h_3_costheta3 = bookHistogram1D(8, 1, 1);
00045       _h_3_psi3 = bookHistogram1D(9, 1, 1);
00046       _h_3_f3 = bookHistogram1D(14, 1, 1);
00047       _h_3_f4 = bookHistogram1D(14, 1, 2);
00048       _h_3_f5 = bookHistogram1D(14, 1, 3);
00049 
00050       _h_4_mNJ = bookHistogram1D(1, 1, 2);
00051       _h_4_X3 = bookHistogram1D(4, 1, 1);
00052       _h_4_X4 = bookHistogram1D(5, 1, 1);
00053       _h_4_costheta3 = bookHistogram1D(10, 1, 1);
00054       _h_4_psi3 = bookHistogram1D(11, 1, 1);
00055       _h_4_f3 = bookHistogram1D(15, 1, 1);
00056       _h_4_f4 = bookHistogram1D(15, 1, 2);
00057       _h_4_f5 = bookHistogram1D(15, 1, 3);
00058       _h_4_XA = bookHistogram1D(17, 1, 1);
00059       _h_4_psiAB = bookHistogram1D(19, 1, 1);
00060       _h_4_fA = bookHistogram1D(21, 1, 1);
00061       _h_4_fB = bookHistogram1D(21, 1, 2);
00062 
00063       _h_5_mNJ = bookHistogram1D(1, 1, 3);
00064       _h_5_X3 = bookHistogram1D(6, 1, 1);
00065       _h_5_X4 = bookHistogram1D(7, 1, 1);
00066       _h_5_costheta3 = bookHistogram1D(12, 1, 1);
00067       _h_5_psi3 = bookHistogram1D(13, 1, 1);
00068       _h_5_f3 = bookHistogram1D(16, 1, 1);
00069       _h_5_f4 = bookHistogram1D(16, 1, 2);
00070       _h_5_f5 = bookHistogram1D(16, 1, 3);
00071       _h_5_XA = bookHistogram1D(18, 1, 1);
00072       _h_5_XC = bookHistogram1D(18, 1, 2);
00073       _h_5_psiAB = bookHistogram1D(20, 1, 1);
00074       _h_5_psiCD = bookHistogram1D(20, 1, 2);
00075       _h_5_fA = bookHistogram1D(22, 1, 1);
00076       _h_5_fB = bookHistogram1D(23, 1, 1);
00077       _h_5_fC = bookHistogram1D(24, 1, 1);
00078       _h_5_fD = bookHistogram1D(25, 1, 1);
00079 
00080     }

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 83 of file CDF_1996_S3349578.cc.

References CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1996_S3349578::_threeJetAnalysis(), Rivet::deltaR(), Rivet::Et(), and Rivet::GeV.

00083                                      {
00084       const double weight = event.weight();
00085 
00086       Jets jets;
00087       FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0);
00088       foreach (const Jet& jet, applyProjection<FastJets>(event, "Jets").jetsByEt()) {
00089         double Et = jet.momentum().Et();
00090         if (Et > 20.0*GeV) {
00091           bool separated=true;
00092           foreach (const Jet& ref, jets) {
00093             if (deltaR(jet.momentum(), ref.momentum())<0.9) {
00094               separated=false;
00095               break;
00096             }
00097           }
00098           if (!separated) continue;
00099           jets.push_back(jet);
00100           jetsystem += jet.momentum();
00101         }
00102         if (jets.size()>=5) break;
00103       }
00104       /// @todo include gaussian jet energy resolution smearing?
00105 
00106       if (jets.size() > 4) {
00107         _fiveJetAnalysis(jets, weight);
00108         jets.resize(4);
00109       }
00110       if (jets.size() > 3) {
00111         _fourJetAnalysis(jets, weight);
00112         jets.resize(3);
00113       }
00114       if (jets.size() > 2) _threeJetAnalysis(jets, weight);
00115     }

void _threeJetAnalysis	(	const Jets &	jets,
		const double &	weight
	)			`[inline]`

Definition at line 120 of file CDF_1996_S3349578.cc.

Referenced by CDF_1996_S3349578::analyze().

00120                                                                    {
00121       getLog() << Log::DEBUG << "3 jet analysis" << std::endl;
00122 
00123       double sumEt=0.0;
00124       FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0);
00125       foreach (const Jet& jet, jets) {
00126         sumEt+=jet.momentum().Et();
00127         jetsystem+=jet.momentum();
00128       }
00129       if (sumEt < 420.0*GeV) return;
00130 
00131       const double m3J = _safeMass(jetsystem);
00132       if (m3J<600*GeV) {
00133         return;
00134       }
00135 
00136       LorentzTransform cms_boost(-jetsystem.boostVector());
00137       vector<FourMomentum> jets3;
00138       foreach (Jet jet, jets) {
00139         jets3.push_back(cms_boost.transform(jet.momentum()));
00140       }
00141       std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending());
00142       FourMomentum p3(jets3[0]);
00143       FourMomentum p4(jets3[1]);
00144       FourMomentum p5(jets3[2]);
00145 
00146       FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m3J, jetsystem.rapidity()));
00147       double costheta3=pAV.vector3().unit().dot(p3.vector3().unit());
00148       if (fabs(costheta3)>0.6) {
00149         return;
00150       }
00151 
00152       double X3 = 2.0*p3.E()/m3J;
00153       if (X3>0.9) {
00154         return;
00155       }
00156 
00157       const double X4 = 2.0*p4.E()/m3J;
00158       const double psi3 = _psi(p3, pAV, p4, p5);
00159       const double f3 = _safeMass(p3)/m3J;
00160       const double f4 = _safeMass(p4)/m3J;
00161       const double f5 = _safeMass(p5)/m3J;
00162 
00163       _h_3_mNJ->fill(m3J, weight);
00164       _h_3_X3->fill(X3, weight);
00165       _h_3_X4->fill(X4, weight);
00166       _h_3_costheta3->fill(costheta3, weight);
00167       _h_3_psi3->fill(psi3, weight);
00168       _h_3_f3->fill(f3, weight);
00169       _h_3_f4->fill(f4, weight);
00170       _h_3_f5->fill(f5, weight);
00171 
00172     }

void _fourJetAnalysis	(	const Jets &	jets,
		const double &	weight
	)			`[inline]`

Definition at line 177 of file CDF_1996_S3349578.cc.

Referenced by CDF_1996_S3349578::analyze().

00177                                                                   {
00178       getLog() << Log::DEBUG << "4 jet analysis" << std::endl;
00179 
00180       double sumEt=0.0;
00181       FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0);
00182       foreach (const Jet& jet, jets) {
00183         sumEt+=jet.momentum().Et();
00184         jetsystem+=jet.momentum();
00185       }
00186       if (sumEt < 420.0*GeV) return;
00187 
00188       const double m4J = _safeMass(jetsystem);
00189       if (m4J < 650*GeV) return;
00190 
00191       LorentzTransform cms_boost(-jetsystem.boostVector());
00192       vector<FourMomentum> jets4;
00193       foreach (Jet jet, jets) {
00194         jets4.push_back(cms_boost.transform(jet.momentum()));
00195       }
00196       std::sort(jets4.begin(), jets4.end(), FourMomentum::byEDescending());
00197 
00198       FourMomentum pA, pB;
00199       vector<FourMomentum> jets3(_reduce(jets4, pA, pB));
00200       std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending());
00201       FourMomentum p3(jets3[0]);
00202       FourMomentum p4(jets3[1]);
00203       FourMomentum p5(jets3[2]);
00204 
00205       FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m4J, jetsystem.rapidity()));
00206       double costheta3=pAV.vector3().unit().dot(p3.vector3().unit());
00207       if (fabs(costheta3)>0.8) {
00208         return;
00209       }
00210 
00211       const double X3 = 2.0*p3.E()/m4J;
00212       if (X3>0.9) {
00213         return;
00214       }
00215 
00216       // fill histograms
00217       const double X4 = 2.0*p4.E()/m4J;
00218       const double psi3 = _psi(p3, pAV, p4, p5);
00219       const double f3 = _safeMass(p3)/m4J;
00220       const double f4 = _safeMass(p4)/m4J;
00221       const double f5 = _safeMass(p5)/m4J;
00222       const double fA = _safeMass(pA)/m4J;
00223       const double fB = _safeMass(pB)/m4J;
00224       const double XA = pA.E()/(pA.E()+pB.E());
00225       const double psiAB = _psi(pA, pB, pA+pB, pAV);
00226 
00227       _h_4_mNJ->fill(m4J, weight);
00228       _h_4_X3->fill(X3, weight);
00229       _h_4_X4->fill(X4, weight);
00230       _h_4_costheta3->fill(costheta3, weight);
00231       _h_4_psi3->fill(psi3, weight);
00232       _h_4_f3->fill(f3, weight);
00233       _h_4_f4->fill(f4, weight);
00234       _h_4_f5->fill(f5, weight);
00235       _h_4_XA->fill(XA, weight);
00236       _h_4_psiAB->fill(psiAB, weight);
00237       _h_4_fA->fill(fA, weight);
00238       _h_4_fB->fill(fB, weight);
00239     }

void _fiveJetAnalysis	(	const Jets &	jets,
		const double &	weight
	)			`[inline]`

Definition at line 244 of file CDF_1996_S3349578.cc.

Referenced by CDF_1996_S3349578::analyze().

00244                                                                   {
00245       getLog() << Log::DEBUG << "5 jet analysis" << std::endl;
00246 
00247       double sumEt=0.0;
00248       FourMomentum jetsystem(0.0, 0.0, 0.0, 0.0);
00249       foreach (const Jet& jet, jets) {
00250         sumEt+=jet.momentum().Et();
00251         jetsystem+=jet.momentum();
00252       }
00253       if (sumEt < 420.0*GeV) return;
00254 
00255       const double m5J = _safeMass(jetsystem);
00256       if (m5J < 750*GeV) return;
00257 
00258       LorentzTransform cms_boost(-jetsystem.boostVector());
00259       vector<FourMomentum> jets5;
00260       foreach (Jet jet, jets) {
00261         jets5.push_back(cms_boost.transform(jet.momentum()));
00262       }
00263       std::sort(jets5.begin(), jets5.end(), FourMomentum::byEDescending());
00264 
00265       FourMomentum pC, pD;
00266       vector<FourMomentum> jets4(_reduce(jets5, pC, pD));
00267       std::sort(jets4.begin(), jets4.end(), FourMomentum::byEDescending());
00268 
00269       FourMomentum pA, pB;
00270       vector<FourMomentum> jets3(_reduce(jets4, pA, pB));
00271       std::sort(jets3.begin(), jets3.end(), FourMomentum::byEDescending());
00272       FourMomentum p3(jets3[0]);
00273       FourMomentum p4(jets3[1]);
00274       FourMomentum p5(jets3[2]);
00275 
00276       // fill histograms
00277       FourMomentum pAV = cms_boost.transform(_avg_beam_in_lab(m5J, jetsystem.rapidity()));
00278       const double costheta3 = pAV.vector3().unit().dot(p3.vector3().unit());
00279       const double X3 = 2.0*p3.E()/m5J;
00280       const double X4 = 2.0*p4.E()/m5J;
00281       const double psi3 = _psi(p3, pAV, p4, p5);
00282       const double f3 = _safeMass(p3)/m5J;
00283       const double f4 = _safeMass(p4)/m5J;
00284       const double f5 = _safeMass(p5)/m5J;
00285       const double fA = _safeMass(pA)/m5J;
00286       const double fB = _safeMass(pB)/m5J;
00287       const double XA = pA.E()/(pA.E()+pB.E());
00288       const double psiAB = _psi(pA, pB, pA+pB, pAV);
00289       const double fC = _safeMass(pC)/m5J;
00290       const double fD = _safeMass(pD)/m5J;
00291       const double XC = pC.E()/(pC.E()+pD.E());
00292       const double psiCD = _psi(pC, pD, pC+pD, pAV);
00293 
00294       _h_5_mNJ->fill(m5J, weight);
00295       _h_5_X3->fill(X3, weight);
00296       _h_5_X4->fill(X4, weight);
00297       _h_5_costheta3->fill(costheta3, weight);
00298       _h_5_psi3->fill(psi3, weight);
00299       _h_5_f3->fill(f3, weight);
00300       _h_5_f4->fill(f4, weight);
00301       _h_5_f5->fill(f5, weight);
00302       _h_5_XA->fill(XA, weight);
00303       _h_5_psiAB->fill(psiAB, weight);
00304       _h_5_fA->fill(fA, weight);
00305       _h_5_fB->fill(fB, weight);
00306       _h_5_XC->fill(XC, weight);
00307       _h_5_psiCD->fill(psiCD, weight);
00308       _h_5_fC->fill(fC, weight);
00309       _h_5_fD->fill(fD, weight);
00310     }

void finalize ( ) [inline, virtual]

Normalise histograms etc., after the run.

Normalise, scale and otherwise manipulate histograms here

Implements Analysis.

Definition at line 314 of file CDF_1996_S3349578.cc.

00314                     {
00315 
00316       /// Normalise, scale and otherwise manipulate histograms here
00317       normalize(_h_3_mNJ, 1.0);
00318       normalize(_h_3_X3, 1.0);
00319       normalize(_h_3_X4, 1.0);
00320       normalize(_h_3_costheta3, 1.0);
00321       normalize(_h_3_psi3, 1.0);
00322       normalize(_h_3_f3, 1.0);
00323       normalize(_h_3_f4, 1.0);
00324       normalize(_h_3_f5, 1.0);
00325 
00326       normalize(_h_4_mNJ, 1.0);
00327       normalize(_h_4_X3, 1.0);
00328       normalize(_h_4_X4, 1.0);
00329       normalize(_h_4_costheta3, 1.0);
00330       normalize(_h_4_psi3, 1.0);
00331       normalize(_h_4_f3, 1.0);
00332       normalize(_h_4_f4, 1.0);
00333       normalize(_h_4_f5, 1.0);
00334       normalize(_h_4_XA, 1.0);
00335       normalize(_h_4_psiAB, 1.0);
00336       normalize(_h_4_fA, 1.0);
00337       normalize(_h_4_fB, 1.0);
00338 
00339       normalize(_h_5_mNJ, 1.0);
00340       normalize(_h_5_X3, 1.0);
00341       normalize(_h_5_X4, 1.0);
00342       normalize(_h_5_costheta3, 1.0);
00343       normalize(_h_5_psi3, 1.0);
00344       normalize(_h_5_f3, 1.0);
00345       normalize(_h_5_f4, 1.0);
00346       normalize(_h_5_f5, 1.0);
00347       normalize(_h_5_XA, 1.0);
00348       normalize(_h_5_XC, 1.0);
00349       normalize(_h_5_psiAB, 1.0);
00350       normalize(_h_5_psiCD, 1.0);
00351       normalize(_h_5_fA, 1.0);
00352       normalize(_h_5_fB, 1.0);
00353       normalize(_h_5_fC, 1.0);
00354       normalize(_h_5_fD, 1.0);
00355 
00356     }

vector<FourMomentum> _reduce	(	const vector< FourMomentum > &	jets,
		FourMomentum &	combined1,
		FourMomentum &	combined2
	)			`[inline, private]`

Definition at line 362 of file CDF_1996_S3349578.cc.

References Rivet::mass2().

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), and CDF_1996_S3349578::_fourJetAnalysis().

00364                                                           {
00365       double minMass2 = 1e9;
00366       size_t idx1(jets.size()), idx2(jets.size());
00367       for (size_t i=0; i<jets.size(); ++i) {
00368         for (size_t j=i+1; j<jets.size(); ++j) {
00369           double mass2 = FourMomentum(jets[i]+jets[j]).mass2();
00370           if (mass2<minMass2) {
00371             idx1=i;
00372             idx2=j;
00373           }
00374         }
00375       }
00376       vector<FourMomentum> newjets;
00377       for (size_t i=0; i<jets.size(); ++i) {
00378         if (i!=idx1 && i!=idx2) newjets.push_back(jets[i]);
00379       }
00380       newjets.push_back(jets[idx1]+jets[idx2]);
00381       combined1 = jets[idx1];
00382       combined2 = jets[idx2];
00383       return newjets;
00384     }

FourMomentum _avg_beam_in_lab	(	const double &	m,
		const double &	y
	)			`[inline, private]`

Definition at line 386 of file CDF_1996_S3349578.cc.

References FourMomentum::boostVector(), and FourMomentum::E().

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), and CDF_1996_S3349578::_threeJetAnalysis().

00386                                                                     {
00387       const double mt = m/2.0;
00388       FourMomentum beam1(mt, 0, 0, mt);
00389       FourMomentum beam2(mt, 0, 0, -mt);
00390       if (fabs(y)>1e-3) {
00391         FourMomentum boostvec(cosh(y), 0.0, 0.0, sinh(y));
00392         LorentzTransform cms_boost(-boostvec.boostVector());
00393         cms_boost = cms_boost.inverse();
00394         beam1=cms_boost.transform(beam1);
00395         beam2=cms_boost.transform(beam2);
00396       }
00397       if (beam1.E()>beam2.E()) {
00398         return beam1-beam2;
00399       }
00400       else {
00401         return beam2-beam1;
00402       }
00403     }

double _psi	(	const FourMomentum &	p1,
		const FourMomentum &	p2,
		const FourMomentum &	p3,
		const FourMomentum &	p4
	)			`[inline, private]`

Definition at line 405 of file CDF_1996_S3349578.cc.

References Vector3::cross(), Vector3::dot(), Rivet::mapAngle0ToPi(), Vector3::unit(), and FourVector::vector3().

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), and CDF_1996_S3349578::_threeJetAnalysis().

00406                                                                 {
00407       Vector3 p1xp2 = p1.vector3().cross(p2.vector3());
00408       Vector3 p3xp4 = p3.vector3().cross(p4.vector3());
00409       return mapAngle0ToPi(acos(p1xp2.unit().dot(p3xp4.unit())));
00410     }

double _safeMass ( const FourMomentum & p ) [inline, private]

Definition at line 412 of file CDF_1996_S3349578.cc.

References Analysis::getLog(), Rivet::m2, FourMomentum::mass2(), Rivet::mass2(), and Log::WARNING.

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), and CDF_1996_S3349578::_threeJetAnalysis().

00412                                             {
00413       double mass2=p.mass2();
00414       if (mass2>0.0) return sqrt(mass2);
00415       else if (mass2<-1.0e-5) {
00416         getLog() << Log::WARNING << "m2 = " << m2 << ". Assuming m2=0." << endl;
00417         return 0.0;
00418       }
00419       else return 0.0;
00420     }

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   }