CDF_2004_S5839831 Class Reference

Inheritance diagram for CDF_2004_S5839831:

Inheritance graph
[legend]

Collaboration diagram for CDF_2004_S5839831:

Collaboration graph
[legend]

List of all members.

Detailed Description

"Acosta" CDF underlying event analysis

Author:
Andy Buckley

Definition at line 20 of file CDF_2004_S5839831.cc.


Public Member Functions

 CDF_2004_S5839831 ()
virtual bool isCompatible (const ParticleName &beam1, const ParticleName &beam2) const
 Is this analysis able to run on the supplied pair of beams?
virtual bool isCompatible (const BeamPair &beams) const
 Is this analysis able to run on the BeamPair beams ?
AnalysisHandler & handler () const
 Access the controlling AnalysisHandler object.
void normalize (AIDA::IHistogram1D *&histo, double norm=1.0)
void scale (AIDA::IHistogram1D *&histo, double scale)
Analysis & setCrossSection (double xs)
 Set the cross section from the generator.
bool needsCrossSection () const
 Return true if this analysis needs to know the process cross-section.
Metadata
Metadata is used for querying from the command line and also for building web pages and the analysis pages in the Rivet manual.

virtual const AnalysisInfo & info () const
 Get the actual AnalysisInfo object in which all this metadata is stored.
virtual std::string name () const
 Get the name of the analysis.
virtual std::string spiresId () const
 Get a the SPIRES/Inspire ID code for this analysis.
virtual std::vector< std::string > authors () const
 Names & emails of paper/analysis authors.
virtual std::string summary () const
 Get a short description of the analysis.
virtual std::string description () const
 Get a full description of the analysis.
virtual std::string runInfo () const
 Information about the events needed as input for this analysis.
virtual std::string experiment () const
 Experiment which performed and published this analysis.
virtual std::string collider () const
 Collider on which the experiment ran.
virtual const BeamPair requiredBeams () const
 Return the pair of incoming beams required by this analysis.
virtual const std::vector
< std::pair< double, double > > & 		energies () const
 Sets of valid beam energy pairs, in GeV.
virtual std::string year () const
 When the original experimental analysis was published.
virtual std::vector< std::string > references () const
 Journal, and preprint references.
virtual std::string status () const
 Whether this analysis is trusted (in any way!).
Run conditions
const ParticlePair & beams () const
 Incoming beams for this run.
const BeamPair beamIds () const
 Incoming beam IDs for this run.
double sqrtS () const
 Centre of mass energy for this run.
Projection "getting" functions
std::set< ConstProjectionPtr > getProjections () const
 Get the contained projections, including recursion.
template<typename PROJ>
const PROJ & getProjection (const std::string &name) const
 Get the named projection, specifying return type via a template argument.
const Projection & getProjection (const std::string &name) const
Projection applying functions
template<typename PROJ>
const PROJ & applyProjection (const Event &evt, const PROJ &proj) const
 Apply the supplied projection on event.
template<typename PROJ>
const PROJ & applyProjection (const Event &evt, const Projection &proj) const
 Apply the supplied projection on event.
template<typename PROJ>
const PROJ & applyProjection (const Event &evt, const std::string &name) const
 Apply the named projection on event.

Protected Member Functions

Log & getLog () const
 Get a Log object based on the name() property of the calling analysis object.
double crossSection () const
 Get the process cross-section in pb. Throws if this hasn't been set.
double crossSectionPerEvent () const
size_t numEvents () const
double sumOfWeights () const
Analysis & setBeams (const ParticleName &beam1, const ParticleName &beam2)
Analysis & setNeedsCrossSection (bool needed)
 Declare whether this analysis needs to know the process cross-section from the generator.
ProjectionHandler & getProjHandler () const
 Get a reference to the ProjectionHandler for this thread.
AIDA analysis infrastructure.
AIDA::IAnalysisFactory & analysisFactory ()
 Access the AIDA analysis factory of the controlling AnalysisHandler object.
AIDA::ITree & tree ()
 Access the AIDA tree of the controlling AnalysisHandler object.
AIDA::IHistogramFactory & histogramFactory ()
 Access the AIDA histogram factory of the controlling AnalysisHandler object.
AIDA::IDataPointSetFactory & datapointsetFactory ()
 Access the AIDA histogram factory of the controlling AnalysisHandler object.
const std::string histoDir () const
 Get the canonical histogram "directory" path for this analysis.
const std::string histoPath (const std::string &hname) const
 Get the canonical histogram path for the named histogram in this analysis.
Internal histogram booking (for use by Analysis sub-classes).
const BinEdges & binEdges (const std::string &hname) const
 Get bin edges for a named histo (using ref AIDA caching).
const BinEdges & binEdges (size_t datasetId, size_t xAxisId, size_t yAxisId) const
 Get bin edges for a numbered histo (using ref AIDA caching).
BinEdges logBinEdges (size_t nbins, double lower, double upper)
 Get bin edges with logarithmic widths.
AIDA::IHistogram1D * bookHistogram1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D * bookHistogram1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D * bookHistogram1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IHistogram1D * bookHistogram1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
Internal profile histogram booking (for use by Analysis sub-classes).
AIDA::IProfile1D * bookProfile1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D * bookProfile1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D * bookProfile1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IProfile1D * bookProfile1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
Internal data point set booking (for use by Analysis sub-classes).
AIDA::IDataPointSet * bookDataPointSet (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IDataPointSet * bookDataPointSet (const std::string &name, size_t npts, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
AIDA::IDataPointSet * bookDataPointSet (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
Projection registration functions
template<typename PROJ>
const PROJ & addProjection (const PROJ &proj, const std::string &name)
const Projection & _addProjection (const Projection &proj, const std::string &name)
 Untemplated function to do the work...

Protected Attributes

string _defaultname
 Name passed to constructor (used to find .info analysis data file, and as a fallback).
shared_ptr< AnalysisInfo > _info
 Pointer to analysis metadata object.
bool _allowProjReg
 Flag to forbid projection registration in analyses until the init phase.

Private Member Functions

ConesInfo _calcTransCones (const double etaLead, const double phiLead, const ParticleVector &tracks)
ConesInfo _calcTransCones (const FourMomentum &leadvec, const ParticleVector &tracks)
Analysis methods
void init ()
void analyze (const Event &event)
 Do the analysis.
void finalize ()

Private Attributes

Histogram collections
AIDA::IProfile1D * _pt90MaxAvg1800
AIDA::IProfile1D * _pt90MinAvg1800
AIDA::IProfile1D * _pt90Max1800
AIDA::IProfile1D * _pt90Min1800
AIDA::IProfile1D * _pt90Diff1800
AIDA::IProfile1D * _pt90Max630
AIDA::IProfile1D * _pt90Min630
AIDA::IProfile1D * _pt90Diff630
AIDA::IProfile1D * _num90Max1800
AIDA::IProfile1D * _num90Min1800
AIDA::IProfile1D * _pTSum1800_2Jet
AIDA::IProfile1D * _pTSum1800_3Jet
AIDA::IProfile1D * _pTSum630_2Jet
AIDA::IProfile1D * _pTSum630_3Jet
AIDA::IHistogram1D * _pt90Dbn1800Et40
AIDA::IHistogram1D * _pt90Dbn1800Et80
AIDA::IHistogram1D * _pt90Dbn1800Et120
AIDA::IHistogram1D * _pt90Dbn1800Et160
AIDA::IHistogram1D * _pt90Dbn1800Et200
AIDA::IHistogram1D * _numTracksDbn1800MB
AIDA::IHistogram1D * _ptDbn1800MB
AIDA::IHistogram1D * _numTracksDbn630MB
AIDA::IHistogram1D * _ptDbn630MB

Friends

class Projectionhandler

Classes

struct  ConesInfo

Constructor & Destructor Documentation

CDF_2004_S5839831 (  )  [inline]

Constructor: cuts on charged final state are $ -1 < \eta < 1 $ and $ p_T > 0.4 $ GeV.

Definition at line 25 of file CDF_2004_S5839831.cc.

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

00026       : Analysis("CDF_2004_S5839831")
00027     {
00028       setBeams(PROTON, ANTIPROTON);
00029     }

Member Function Documentation

ConesInfo _calcTransCones ( const double  etaLead,
const double  phiLead,
const ParticleVector &  tracks 
) [inline, private]

Definition at line 41 of file CDF_2004_S5839831.cc.

References Log::DEBUG, Rivet::deltaR(), Analysis::getLog(), Rivet::GeV, Rivet::mapAngle0To2Pi(), Particle::momentum(), CDF_2004_S5839831::ConesInfo::numMax, CDF_2004_S5839831::ConesInfo::numMin, Rivet::PI, FourMomentum::pT(), CDF_2004_S5839831::ConesInfo::ptDiff, CDF_2004_S5839831::ConesInfo::ptMax, and CDF_2004_S5839831::ConesInfo::ptMin.

Referenced by CDF_2004_S5839831::_calcTransCones(), and CDF_2004_S5839831::analyze().

00042                                                             {
00043       const double phiTransPlus = mapAngle0To2Pi(phiLead + PI/2.0);
00044       const double phiTransMinus = mapAngle0To2Pi(phiLead - PI/2.0);
00045       getLog() << Log::DEBUG << "phi_lead = " << phiLead
00046                << " -> trans = (" << phiTransPlus
00047                << ", " << phiTransMinus << ")" << endl;
00048    
00049       unsigned int numPlus(0), numMinus(0);
00050       double ptPlus(0), ptMinus(0);
00051       // Run over all charged tracks
00052       foreach (const Particle& t, tracks) {
00053         FourMomentum trackMom = t.momentum();
00054         const double pt = trackMom.pT();
00055      
00056         // Find if track mom is in either transverse cone
00057         if (deltaR(trackMom, etaLead, phiTransPlus) < 0.7) {
00058           ptPlus += pt;
00059           numPlus += 1;
00060         } else if (deltaR(trackMom, etaLead, phiTransMinus) < 0.7) {
00061           ptMinus += pt;
00062           numMinus += 1;
00063         }
00064       }
00065    
00066       ConesInfo rtn;
00067       // Assign N_{min,max} from N_{plus,minus}
00068       rtn.numMax = (ptPlus >= ptMinus) ? numPlus : numMinus;
00069       rtn.numMin = (ptPlus >= ptMinus) ? numMinus : numPlus;
00070       // Assign pT_{min,max} from pT_{plus,minus}
00071       rtn.ptMax = (ptPlus >= ptMinus) ? ptPlus : ptMinus;
00072       rtn.ptMin = (ptPlus >= ptMinus) ? ptMinus : ptPlus;
00073       rtn.ptDiff = fabs(rtn.ptMax - rtn.ptMin);
00074    
00075       getLog() << Log::DEBUG << "Min cone has " << rtn.numMin << " tracks -> "
00076                << "pT_min = " << rtn.ptMin/GeV << " GeV" << endl;
00077       getLog() << Log::DEBUG << "Max cone has " << rtn.numMax << " tracks -> "
00078                << "pT_max = " << rtn.ptMax/GeV << " GeV" << endl;
00079    
00080       return rtn;
00081     }

ConesInfo _calcTransCones ( const FourMomentum &  leadvec,
const ParticleVector &  tracks 
) [inline, private]

Definition at line 84 of file CDF_2004_S5839831.cc.

References CDF_2004_S5839831::_calcTransCones(), FourVector::azimuthalAngle(), and FourVector::pseudorapidity().

00085                                                             {
00086       const double etaLead = leadvec.pseudorapidity();
00087       const double phiLead = leadvec.azimuthalAngle();
00088       return _calcTransCones(etaLead, phiLead, tracks);
00089     }

void init (  )  [inline, private, 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 95 of file CDF_2004_S5839831.cc.

References CDF_2004_S5839831::_num90Max1800, CDF_2004_S5839831::_num90Min1800, CDF_2004_S5839831::_numTracksDbn1800MB, CDF_2004_S5839831::_numTracksDbn630MB, CDF_2004_S5839831::_pt90Dbn1800Et120, CDF_2004_S5839831::_pt90Dbn1800Et160, CDF_2004_S5839831::_pt90Dbn1800Et200, CDF_2004_S5839831::_pt90Dbn1800Et40, CDF_2004_S5839831::_pt90Dbn1800Et80, CDF_2004_S5839831::_pt90Diff1800, CDF_2004_S5839831::_pt90Diff630, CDF_2004_S5839831::_pt90Max1800, CDF_2004_S5839831::_pt90Max630, CDF_2004_S5839831::_pt90MaxAvg1800, CDF_2004_S5839831::_pt90Min1800, CDF_2004_S5839831::_pt90Min630, CDF_2004_S5839831::_pt90MinAvg1800, CDF_2004_S5839831::_ptDbn1800MB, CDF_2004_S5839831::_ptDbn630MB, CDF_2004_S5839831::_pTSum1800_2Jet, CDF_2004_S5839831::_pTSum1800_3Jet, CDF_2004_S5839831::_pTSum630_2Jet, CDF_2004_S5839831::_pTSum630_3Jet, ProjectionApplier::addProjection(), Analysis::bookHistogram1D(), Analysis::bookProfile1D(), FastJets::CDFJETCLU, Rivet::fuzzyEquals(), Rivet::GeV, and Analysis::sqrtS().

00095                 {
00096       // Set up projections
00097       addProjection(TriggerCDFRun0Run1(), "Trigger");
00098       addProjection(Beam(), "Beam");
00099       const FinalState calofs(-1.2, 1.2);
00100       addProjection(calofs, "CaloFS");
00101       addProjection(FastJets(calofs, FastJets::CDFJETCLU, 0.7), "Jets");
00102       const ChargedFinalState trackfs(-1.2, 1.2, 0.4*GeV);
00103       addProjection(trackfs, "TrackFS");
00104       // Restrict tracks to |eta| < 0.7 for the min bias part.
00105       const ChargedFinalState mbfs(-0.7, 0.7, 0.4*GeV);
00106       addProjection(mbfs, "MBFS");
00107       // Restrict tracks to |eta| < 1 for the Swiss-Cheese part.
00108       const ChargedFinalState cheesefs(-1.0, 1.0, 0.4*GeV);
00109       addProjection(cheesefs, "CheeseFS");
00110       addProjection(FastJets(cheesefs, FastJets::CDFJETCLU, 0.7), "CheeseJets");
00111 
00112       // Book histograms
00113       if (fuzzyEquals(sqrtS()/GeV, 1800, 1E-3)) {
00114         _pt90MaxAvg1800 = bookProfile1D(1, 1, 1);
00115         _pt90MinAvg1800 = bookProfile1D(1, 1, 2);
00116         _pt90Max1800 = bookProfile1D(2, 1, 1);
00117         _pt90Min1800 = bookProfile1D(2, 1, 2);
00118         _pt90Diff1800 = bookProfile1D(2, 1, 3);
00119         _num90Max1800 = bookProfile1D(4, 1, 1);
00120         _num90Min1800 = bookProfile1D(4, 1, 2);
00121         _pTSum1800_2Jet = bookProfile1D(7, 1, 1);
00122         _pTSum1800_3Jet = bookProfile1D(7, 1, 2);
00123 
00124         _pt90Dbn1800Et40 = bookHistogram1D(3, 1, 1);
00125         _pt90Dbn1800Et80 = bookHistogram1D(3, 1, 2);
00126         _pt90Dbn1800Et120 = bookHistogram1D(3, 1, 3);
00127         _pt90Dbn1800Et160 = bookHistogram1D(3, 1, 4);
00128         _pt90Dbn1800Et200 = bookHistogram1D(3, 1, 5);
00129         _numTracksDbn1800MB = bookHistogram1D(5, 1, 1);
00130         _ptDbn1800MB = bookHistogram1D(6, 1, 1);
00131       } else if (fuzzyEquals(sqrtS()/GeV, 630, 1E-3)) {
00132         _pt90Max630 = bookProfile1D(8, 1, 1);
00133         _pt90Min630 = bookProfile1D(8, 1, 2);
00134         _pt90Diff630 = bookProfile1D(8, 1, 3);
00135         _pTSum630_2Jet = bookProfile1D(9, 1, 1);
00136         _pTSum630_3Jet = bookProfile1D(9, 1, 2);
00137 
00138         _numTracksDbn630MB = bookHistogram1D(10, 1, 1);
00139         _ptDbn630MB = bookHistogram1D(11, 1, 1);
00140       }
00141     }

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

Do the analysis.

Implements Analysis.

Definition at line 145 of file CDF_2004_S5839831.cc.

References CDF_2004_S5839831::_calcTransCones(), CDF_2004_S5839831::_num90Max1800, CDF_2004_S5839831::_num90Min1800, CDF_2004_S5839831::_numTracksDbn1800MB, CDF_2004_S5839831::_numTracksDbn630MB, CDF_2004_S5839831::_pt90Dbn1800Et120, CDF_2004_S5839831::_pt90Dbn1800Et160, CDF_2004_S5839831::_pt90Dbn1800Et200, CDF_2004_S5839831::_pt90Dbn1800Et40, CDF_2004_S5839831::_pt90Dbn1800Et80, CDF_2004_S5839831::_pt90Diff1800, CDF_2004_S5839831::_pt90Diff630, CDF_2004_S5839831::_pt90Max1800, CDF_2004_S5839831::_pt90Max630, CDF_2004_S5839831::_pt90MaxAvg1800, CDF_2004_S5839831::_pt90Min1800, CDF_2004_S5839831::_pt90Min630, CDF_2004_S5839831::_pt90MinAvg1800, CDF_2004_S5839831::_ptDbn1800MB, CDF_2004_S5839831::_ptDbn630MB, CDF_2004_S5839831::_pTSum1800_2Jet, CDF_2004_S5839831::_pTSum1800_3Jet, CDF_2004_S5839831::_pTSum630_2Jet, CDF_2004_S5839831::_pTSum630_3Jet, FourVector::azimuthalAngle(), Log::DEBUG, Rivet::deltaR(), Rivet::fuzzyEquals(), Analysis::getLog(), Rivet::GeV, Rivet::inRange(), Rivet::isZero(), Rivet::momentum(), Particle::momentum(), CDF_2004_S5839831::ConesInfo::numMax, CDF_2004_S5839831::ConesInfo::numMin, Rivet::particles(), FourVector::pseudorapidity(), Rivet::pseudorapidity(), FourMomentum::pT(), CDF_2004_S5839831::ConesInfo::ptDiff, CDF_2004_S5839831::ConesInfo::ptMax, CDF_2004_S5839831::ConesInfo::ptMin, Analysis::sqrtS(), Log::TRACE, and vetoEvent.

00145                                      {
00146       // Trigger
00147       const bool trigger = applyProjection<TriggerCDFRun0Run1>(event, "Trigger").minBiasDecision();
00148       if (!trigger) vetoEvent;
00149 
00150       // Get sqrt(s) and event weight
00151       const double sqrtS = applyProjection<Beam>(event, "Beam").sqrtS();
00152       const double weight = event.weight();
00153    
00154       {
00155         getLog() << Log::DEBUG << "Running max/min analysis" << endl;
00156         vector<Jet> jets = applyProjection<JetAlg>(event, "Jets").jetsByE();
00157         if (!jets.empty()) {
00158           // Leading jet must be in central |eta| < 0.5 region
00159           const Jet leadingjet = jets.front();
00160           const double etaLead = leadingjet.momentum().eta();
00161           // Get Et of the leading jet: used to bin histograms
00162           const double ETlead = leadingjet.EtSum();
00163           getLog() << Log::DEBUG << "Leading Et = " << ETlead/GeV << " GeV" << endl;
00164           if (fabs(etaLead) > 0.5 && ETlead < 15*GeV) {
00165             getLog() << Log::DEBUG << "Leading jet eta = " << etaLead
00166                      << " not in |eta| < 0.5 & pT > 15 GeV" << endl;
00167           } else {
00168             // Multiplicity & pT distributions for sqrt(s) = 630 GeV, 1800 GeV
00169             const ParticleVector tracks = applyProjection<FinalState>(event, "TrackFS").particles();
00170             const ConesInfo cones = _calcTransCones(leadingjet.momentum(), tracks);
00171             if (fuzzyEquals(sqrtS/GeV, 630)) {
00172               _pt90Max630->fill(ETlead/GeV, cones.ptMax/GeV, weight);
00173               _pt90Min630->fill(ETlead/GeV, cones.ptMin/GeV, weight);
00174               _pt90Diff630->fill(ETlead/GeV, cones.ptDiff/GeV, weight);
00175             } else if (fuzzyEquals(sqrtS/GeV, 1800)) {
00176               _num90Max1800->fill(ETlead/GeV, cones.numMax, weight);
00177               _num90Min1800->fill(ETlead/GeV, cones.numMin, weight);
00178               _pt90Max1800->fill(ETlead/GeV, cones.ptMax/GeV, weight);
00179               _pt90Min1800->fill(ETlead/GeV, cones.ptMin/GeV, weight);
00180               _pt90Diff1800->fill(ETlead/GeV, cones.ptDiff/GeV, weight);
00181               _pt90MaxAvg1800->fill(ETlead/GeV, cones.ptMax/GeV, weight); // /numMax
00182               _pt90MinAvg1800->fill(ETlead/GeV, cones.ptMin/GeV, weight); // /numMin
00183               //
00184               const double ptTransTotal = cones.ptMax + cones.ptMin;
00185               if (inRange(ETlead/GeV, 40., 80.)) {
00186                 _pt90Dbn1800Et40->fill(ptTransTotal/GeV, weight);
00187               } else if (inRange(ETlead/GeV, 80., 120.)) {
00188                 _pt90Dbn1800Et80->fill(ptTransTotal/GeV, weight);
00189               } else if (inRange(ETlead/GeV, 120., 160.)) {
00190                 _pt90Dbn1800Et120->fill(ptTransTotal/GeV, weight);
00191               } else if (inRange(ETlead/GeV, 160., 200.)) {
00192                 _pt90Dbn1800Et160->fill(ptTransTotal/GeV, weight);
00193               } else if (inRange(ETlead/GeV, 200., 270.)) {
00194                 _pt90Dbn1800Et200->fill(ptTransTotal/GeV, weight);
00195               }
00196             }
00197          
00198           }
00199         }
00200       }
00201    
00202    
00203       // Fill min bias total track multiplicity histos
00204       {
00205         getLog() << Log::DEBUG << "Running min bias multiplicity analysis" << endl;
00206         const ParticleVector mbtracks = applyProjection<FinalState>(event, "MBFS").particles();
00207         if (fuzzyEquals(sqrtS/GeV, 1800)) {
00208           _numTracksDbn1800MB->fill(mbtracks.size(), weight);
00209         } else if (fuzzyEquals(sqrtS/GeV, 630)) {
00210           _numTracksDbn630MB->fill(mbtracks.size(), weight);
00211         }
00212         // Run over all charged tracks
00213         foreach (const Particle& t, mbtracks) {
00214           FourMomentum trackMom = t.momentum();
00215           const double pt = trackMom.pT();
00216           // Plot total pT distribution for min bias
00217           if (fuzzyEquals(sqrtS/GeV, 1800)) {
00218             _ptDbn1800MB->fill(pt/GeV, weight);
00219           } else if (fuzzyEquals(sqrtS/GeV, 630)) {
00220             _ptDbn630MB->fill(pt/GeV, weight);
00221           }
00222         }
00223       }
00224    
00225    
00226    
00227       // Construct "Swiss Cheese" pT distributions, with pT contributions from
00228       // tracks within R = 0.7 of the 1st, 2nd (and 3rd) jets being ignored. A
00229       // different set of charged tracks, with |eta| < 1.0, is used here, and all
00230       // the removed jets must have Et > 5 GeV.
00231       {
00232         getLog() << Log::DEBUG << "Running Swiss Cheese analysis" << endl;
00233         const ParticleVector cheesetracks = applyProjection<FinalState>(event, "CheeseFS").particles();
00234         vector<Jet> cheesejets = applyProjection<JetAlg>(event, "Jets").jetsByE();
00235         if (cheesejets.empty()) {
00236           getLog() << Log::DEBUG << "No 'cheese' jets found in event" << endl;
00237           return;
00238         }
00239         if (cheesejets.size() > 1 &&
00240             fabs(cheesejets[0].momentum().pseudorapidity()) <= 0.5 &&
00241             cheesejets[0].momentum().Et()/GeV > 5.0 &&
00242             cheesejets[1].momentum().Et()/GeV > 5.0) {
00243        
00244           const double cheeseETlead = cheesejets[0].momentum().Et();
00245        
00246           const double eta1 = cheesejets[0].momentum().pseudorapidity();
00247           const double phi1 = cheesejets[0].momentum().azimuthalAngle();
00248           const double eta2 = cheesejets[1].momentum().pseudorapidity();
00249           const double phi2 = cheesejets[1].momentum().azimuthalAngle();
00250        
00251           double ptSumSub2(0), ptSumSub3(0);
00252           foreach (const Particle& t, cheesetracks) {
00253             FourMomentum trackMom = t.momentum();
00254             const double pt = trackMom.pT();
00255          
00256             // Subtracting 2 leading jets
00257             const double deltaR1 = deltaR(trackMom, eta1, phi1);
00258             const double deltaR2 = deltaR(trackMom, eta2, phi2);
00259             getLog() << Log::TRACE << "Track vs jet(1): "
00260                      << "|(" << trackMom.pseudorapidity() << ", " << trackMom.azimuthalAngle() << ") - "
00261                      << "|(" << eta1 << ", " << phi1 << ")| = " << deltaR1 << endl;
00262             getLog() << Log::TRACE << "Track vs jet(2): "
00263                      << "|(" << trackMom.pseudorapidity() << ", " << trackMom.azimuthalAngle() << ") - "
00264                      << "|(" << eta2 << ", " << phi2 << ")| = " << deltaR2 << endl;
00265             if (deltaR1 > 0.7 && deltaR2 > 0.7) {
00266               ptSumSub2 += pt;
00267            
00268               // Subtracting 3rd leading jet
00269               if (cheesejets.size() > 2 &&
00270                   cheesejets[2].momentum().Et()/GeV > 5.0) {
00271                 const double eta3 = cheesejets[2].momentum().pseudorapidity();
00272                 const double phi3 = cheesejets[2].momentum().azimuthalAngle();
00273                 const double deltaR3 = deltaR(trackMom, eta3, phi3);
00274                 getLog() << Log::TRACE << "Track vs jet(3): "
00275                          << "|(" << trackMom.pseudorapidity() << ", " << trackMom.azimuthalAngle() << ") - "
00276                          << "|(" << eta3 << ", " << phi3 << ")| = " << deltaR3 << endl;
00277                 if (deltaR3 > 0.7) {
00278                   ptSumSub3 += pt;
00279                 }
00280               }
00281             }
00282           }
00283        
00284           // Swiss Cheese sub 2,3 jets distributions for sqrt(s) = 630 GeV, 1800 GeV
00285           if (fuzzyEquals(sqrtS/GeV, 630)) {
00286             if (!isZero(ptSumSub2)) _pTSum630_2Jet->fill(cheeseETlead/GeV, ptSumSub2/GeV, weight);
00287             if (!isZero(ptSumSub3))_pTSum630_3Jet->fill(cheeseETlead/GeV, ptSumSub3/GeV, weight);
00288           } else if (fuzzyEquals(sqrtS/GeV, 1800)) {
00289             if (!isZero(ptSumSub2))_pTSum1800_2Jet->fill(cheeseETlead/GeV, ptSumSub2/GeV, weight);
00290             if (!isZero(ptSumSub3))_pTSum1800_3Jet->fill(cheeseETlead/GeV, ptSumSub3/GeV, weight);
00291           }
00292        
00293         }
00294       }
00295    
00296     }

void finalize (  )  [inline, private, virtual]

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

Implements Analysis.

Definition at line 299 of file CDF_2004_S5839831.cc.

References CDF_2004_S5839831::_numTracksDbn1800MB, CDF_2004_S5839831::_numTracksDbn630MB, CDF_2004_S5839831::_pt90Dbn1800Et120, CDF_2004_S5839831::_pt90Dbn1800Et160, CDF_2004_S5839831::_pt90Dbn1800Et200, CDF_2004_S5839831::_pt90Dbn1800Et40, CDF_2004_S5839831::_pt90Dbn1800Et80, CDF_2004_S5839831::_ptDbn1800MB, CDF_2004_S5839831::_ptDbn630MB, Rivet::fuzzyEquals(), Rivet::GeV, Analysis::normalize(), and Analysis::sqrtS().

00299                     {
00300       /// @todo Take these normalisations from the data histo (it can't come from just the MC)
00301 
00302       if (fuzzyEquals(sqrtS()/GeV, 1800, 1E-3)) {
00303         // Normalize to actual number of entries in pT dbn histos...
00304         normalize(_pt90Dbn1800Et40,  1656.75); // norm OK
00305         normalize(_pt90Dbn1800Et80,  4657.5); // norm OK
00306         normalize(_pt90Dbn1800Et120, 5395.5); // norm OK
00307         normalize(_pt90Dbn1800Et160, 7248.75); // norm OK
00308         normalize(_pt90Dbn1800Et200, 2442.0); // norm OK
00309       }
00310    
00311       // ...and for min bias distributions:
00312       if (fuzzyEquals(sqrtS()/GeV, 1800, 1E-3)) {
00313         normalize(_numTracksDbn1800MB, 309718.25); // norm OK
00314         normalize(_ptDbn1800MB, 33600.0); // norm OK
00315       } else if (fuzzyEquals(sqrtS()/GeV, 630, 1E-3)) {
00316         normalize(_numTracksDbn630MB, 1101024.0); // norm OK
00317         normalize(_ptDbn630MB, 105088.0); // norm OK
00318       }
00319     }

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

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

Definition at line 125 of file Analysis.cc.

References Analysis::_info.

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

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

string name (  )  const [virtual, inherited]

Get the name of the analysis.

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

Implements ProjectionApplier.

Definition at line 130 of file Analysis.cc.

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

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

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

string spiresId (  )  const [virtual, inherited]

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

Definition at line 135 of file Analysis.cc.

References Analysis::_info.

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

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

Names & emails of paper/analysis authors.

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

Definition at line 140 of file Analysis.cc.

References Analysis::_info.

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

string summary (  )  const [virtual, inherited]

Get a short description of the analysis.

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

Definition at line 145 of file Analysis.cc.

References Analysis::_info.

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

string description (  )  const [virtual, inherited]

Get a full description of the analysis.

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

Definition at line 150 of file Analysis.cc.

References Analysis::_info.

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

string runInfo (  )  const [virtual, inherited]

Information about the events needed as input for this analysis.

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

Definition at line 155 of file Analysis.cc.

References Analysis::_info.

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

string experiment (  )  const [virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 164 of file Analysis.cc.

References Analysis::_info.

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

string collider (  )  const [virtual, inherited]

Collider on which the experiment ran.

Definition at line 169 of file Analysis.cc.

References Analysis::_info.

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

const BeamPair requiredBeams (  )  const [virtual, inherited]

Return the pair of incoming beams required by this analysis.

Definition at line 189 of file Analysis.cc.

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

Referenced by Analysis::isCompatible().

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

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

Sets of valid beam energy pairs, in GeV.

Definition at line 160 of file Analysis.cc.

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

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

string year (  )  const [virtual, inherited]

When the original experimental analysis was published.

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

Definition at line 174 of file Analysis.cc.

References Analysis::_info.

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

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

Journal, and preprint references.

Definition at line 179 of file Analysis.cc.

References Analysis::_info.

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

string status (  )  const [virtual, inherited]

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

Definition at line 184 of file Analysis.cc.

References Analysis::_info.

Referenced by AnalysisHandler::init().

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

const ParticlePair & beams (  )  const [inherited]

Incoming beams for this run.

Definition at line 79 of file Analysis.cc.

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

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

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

const BeamPair beamIds (  )  const [inherited]

Incoming beam IDs for this run.

Definition at line 83 of file Analysis.cc.

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

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

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

double sqrtS (  )  const [inherited]

Centre of mass energy for this run.

Definition at line 75 of file Analysis.cc.

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

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

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

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

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

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

Definition at line 200 of file Analysis.cc.

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

Referenced by AnalysisHandler::removeIncompatibleAnalyses().

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

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

Is this analysis able to run on the BeamPair beams ?

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

Definition at line 207 of file Analysis.cc.

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

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

AnalysisHandler & handler (  )  const [inherited]

Access the controlling AnalysisHandler object.

Definition at line 244 of file Analysis.cc.

References Analysis::_analysishandler.

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

00244                                            {
00245     return *_analysishandler;
00246   }

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

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

Parameters:
histo The histogram to be normalised.
norm The new area of the histogram.
Warning:
The old histogram will be deleted, and its pointer set to zero.

Definition at line 487 of file Analysis.cc.

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

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

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

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

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

Parameters:
histo The histogram to be scaled.
scale The factor used to multiply the histogram bin heights.
Warning:
The old histogram will be deleted, and its pointer set to zero.

Definition at line 512 of file Analysis.cc.

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

Referenced by UA5_1989_S1926373::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_ZJETS::finalize(), MC_WJETS::finalize(), MC_TTBAR::finalize(), MC_PHOTONJETS::finalize(), MC_JetAnalysis::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_2010_S8566488::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(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7540469::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(), 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(), and Analysis::normalize().

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

Analysis & setCrossSection ( double  xs  )  [inherited]

Set the cross section from the generator.

Definition at line 214 of file Analysis.cc.

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

Referenced by AnalysisHandler::setCrossSection().

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

bool needsCrossSection (  )  const [inherited]

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

Definition at line 220 of file Analysis.cc.

References Analysis::_needsCrossSection.

Referenced by AnalysisHandler::needCrossSection().

00220                                          {
00221     return _needsCrossSection;
00222   }

Log & getLog (  )  const [protected, inherited]

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

Reimplemented from ProjectionApplier.

Definition at line 106 of file Analysis.cc.

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

Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), CDF_2004_S5839831::_calcTransCones(), CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), D0_1996_S3214044::_safeMass(), CDF_1997_S3541940::_safeMass(), CDF_1996_S3349578::_safeMass(), CDF_1996_S3349578::_threeJetAnalysis(), ZEUS_2001_S4815815::analyze(), UA5_1986_S1583476::analyze(), UA1_1990_S2044935::analyze(), STAR_2009_UE_HELEN::analyze(), STAR_2008_S7993412::analyze(), STAR_2006_S6870392::analyze(), STAR_2006_S6860818::analyze(), STAR_2006_S6500200::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), OPAL_1998_S3780481::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_LEADINGJETS::analyze(), JADE_OPAL_2000_S4300807::analyze(), H1_1994_S2919893::analyze(), ExampleAnalysis::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1995_S3137023::analyze(), D0_2009_S8349509::analyze(), D0_2009_S8202443::analyze(), D0_2008_S7863608::analyze(), D0_2008_S7837160::analyze(), D0_2008_S7719523::analyze(), D0_2008_S7662670::analyze(), D0_2008_S7554427::analyze(), D0_2008_S6879055::analyze(), D0_2007_S7075677::analyze(), D0_2006_S6438750::analyze(), D0_2004_S5992206::analyze(), D0_2001_S4674421::analyze(), D0_1998_S3711838::analyze(), CDF_2009_S8383952::analyze(), CDF_2009_S8233977::analyze(), CDF_2008_S8095620::analyze(), CDF_2008_S7782535::analyze(), CDF_2008_S7540469::analyze(), CDF_2008_NOTE_9351::analyze(), CDF_2008_LEADINGJETS::analyze(), CDF_2006_S6653332::analyze(), CDF_2005_S6217184::analyze(), CDF_2004_S5839831::analyze(), CDF_2001_S4751469::analyze(), CDF_2000_S4155203::analyze(), CDF_1994_S2952106::analyze(), CDF_1991_S2313472::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_1991_S2435284::analyze(), CDF_2009_S8057893::CDF_2009_S8057893::analyze(), UA5_1986_S1583476::finalize(), UA1_1990_S2044935::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), D0_2001_S4674421::finalize(), CDF_2009_S8233977::finalize(), CDF_2007_S7057202::finalize(), CDF_2006_S6653332::finalize(), ZEUS_2001_S4815815::init(), JADE_OPAL_2000_S4300807::init(), ALEPH_2004_S5765862::init(), Analysis::normalize(), and Analysis::scale().

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

double crossSection (  )  const [protected, inherited]

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

Definition at line 229 of file Analysis.cc.

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

Referenced by STAR_2006_S6870392::finalize(), MC_ZJETS::finalize(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_DIPHOTON::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), 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_2007_S7057202::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(), and CDF_1998_S3618439::finalize().

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

double crossSectionPerEvent (  )  const [protected, inherited]

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

Definition at line 237 of file Analysis.cc.

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

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

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

size_t numEvents (  )  const [protected, inherited]

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

Definition at line 112 of file Analysis.cc.

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

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

double sumOfWeights (  )  const [protected, inherited]

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

Definition at line 117 of file Analysis.cc.

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

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

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

IAnalysisFactory & analysisFactory (  )  [protected, inherited]

Access the AIDA analysis factory of the controlling AnalysisHandler object.

Definition at line 55 of file Analysis.cc.

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

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

ITree & tree (  )  [protected, inherited]

Access the AIDA tree of the controlling AnalysisHandler object.

Definition at line 60 of file Analysis.cc.

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

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

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

IHistogramFactory & histogramFactory (  )  [protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 65 of file Analysis.cc.

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

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

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

IDataPointSetFactory & datapointsetFactory (  )  [protected, inherited]

Access the AIDA histogram factory of the controlling AnalysisHandler object.

Definition at line 70 of file Analysis.cc.

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

Referenced by Analysis::scale().

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

const string histoDir (  )  const [protected, inherited]

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

Definition at line 88 of file Analysis.cc.

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

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

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

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

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

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

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

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

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

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

Definition at line 286 of file Analysis.cc.

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

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

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

Get bin edges with logarithmic widths.

Definition at line 292 of file Analysis.cc.

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

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

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

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

Referenced by ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_1998_S3780481::init(), MC_ZJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_DIPHOTON::init(), MC_DIJET::init(), JADE_OPAL_2000_S4300807::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1995_S3137023::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1991_S2313472::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2010_S8591806::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), and ALEPH_1991_S2435284::init().

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

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

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

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

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

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

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

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

Referenced by UA5_1988_S1867512::init(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2006_S6860818::init(), MC_PHOTONJETUE::init(), MC_LEADINGJETS::init(), H1_2000_S4129130::init(), H1_1994_S2919893::init(), DELPHI_2002_069_CONF_603::init(), D0_1996_S3324664::init(), CDF_2009_S8233977::init(), CDF_2008_S7782535::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2005_S6217184::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), and ATLAS_2010_S8591806::init().

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Set the colliding beam pair.

Deprecated:
Use .info file and AnalysisInfo class instead

Definition at line 193 of file Analysis.cc.

References Analysis::_info.

Referenced by ALEPH_1991_S2435284::ALEPH_1991_S2435284(), ALEPH_1996_S3196992::ALEPH_1996_S3196992(), ALEPH_1996_S3486095::ALEPH_1996_S3486095(), ALEPH_2004_S5765862::ALEPH_2004_S5765862(), BELLE_2006_S6265367::BELLE_2006_S6265367(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1990_S2089246::CDF_1990_S2089246(), CDF_1991_S2313472::CDF_1991_S2313472(), CDF_1994_S2952106::CDF_1994_S2952106(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1996_S3349578::CDF_1996_S3349578(), CDF_1996_S3418421::CDF_1996_S3418421(), CDF_1997_S3541940::CDF_1997_S3541940(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2001_S4751469::CDF_2001_S4751469(), CDF_2002_S4796047::CDF_2002_S4796047(), CDF_2004_S5839831::CDF_2004_S5839831(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2005_S6217184::CDF_2005_S6217184(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_LEADINGJETS::CDF_2008_LEADINGJETS(), CDF_2008_NOTE_9351::CDF_2008_NOTE_9351(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7782535::CDF_2008_S7782535(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2008_S8095620::CDF_2008_S8095620(), CDF_2009_S8057893::CDF_2009_S8057893::CDF_2009_S8057893(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_1998_S3711838::D0_1998_S3711838(), D0_2001_S4674421::D0_2001_S4674421(), D0_2004_S5992206::D0_2004_S5992206(), D0_2006_S6438750::D0_2006_S6438750(), D0_2007_S7075677::D0_2007_S7075677(), D0_2008_S6879055::D0_2008_S6879055(), D0_2008_S7554427::D0_2008_S7554427(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7837160::D0_2008_S7837160(), D0_2008_S7863608::D0_2008_S7863608(), D0_2009_S8202443::D0_2009_S8202443(), D0_2009_S8320160::D0_2009_S8320160(), D0_2009_S8349509::D0_2009_S8349509(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), DELPHI_1995_S3137023::DELPHI_1995_S3137023(), DELPHI_2002_069_CONF_603::DELPHI_2002_069_CONF_603(), E735_1998_S3905616::E735_1998_S3905616(), H1_1994_S2919893::H1_1994_S2919893(), H1_1995_S3167097::H1_1995_S3167097(), H1_2000_S4129130::H1_2000_S4129130(), JADE_OPAL_2000_S4300807::JADE_OPAL_2000_S4300807(), OPAL_1998_S3780481::OPAL_1998_S3780481(), PDG_HADRON_MULTIPLICITIES::PDG_HADRON_MULTIPLICITIES(), PDG_HADRON_MULTIPLICITIES_RATIOS::PDG_HADRON_MULTIPLICITIES_RATIOS(), SFM_1984_S1178091::SFM_1984_S1178091(), STAR_2006_S6500200::STAR_2006_S6500200(), STAR_2006_S6860818::STAR_2006_S6860818(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7993412::STAR_2008_S7993412(), STAR_2009_UE_HELEN::STAR_2009_UE_HELEN(), UA1_1990_S2044935::UA1_1990_S2044935(), UA5_1986_S1583476::UA5_1986_S1583476(), UA5_1988_S1867512::UA5_1988_S1867512(), UA5_1989_S1926373::UA5_1989_S1926373(), and ZEUS_2001_S4815815::ZEUS_2001_S4815815().

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

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

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

Definition at line 224 of file Analysis.cc.

References Analysis::_needsCrossSection.

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

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

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

Get the contained projections, including recursion.

Definition at line 43 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

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

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

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

Definition at line 50 of file ProjectionApplier.hh.

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

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

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

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

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

Definition at line 58 of file ProjectionApplier.hh.

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

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

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

Apply the supplied projection on event.

Definition at line 68 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

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

Apply the supplied projection on event.

Definition at line 75 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

Apply the named projection on event.

Definition at line 82 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

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

Get a reference to the ProjectionHandler for this thread.

Definition at line 95 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

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

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

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

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

Definition at line 115 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

Referenced by ZFinder::_init(), WFinder::_init(), VetoedFinalState::addVetoOnThisFinalState(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ClusteredPhotons::ClusteredPhotons(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), FinalStateHCM::FinalStateHCM(), FoxWolframMoments::FoxWolframMoments(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_1998_S3780481::init(), MC_ZJETS::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_LEADINGJETS::init(), MC_JETS::init(), MC_DIPHOTON::init(), MC_DIJET::init(), JADE_OPAL_2000_S4300807::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1995_S3137023::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_1998_S3711838::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6217184::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1991_S2313472::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2010_S8591806::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), ALEPH_1991_S2435284::init(), IsolationProjection::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), KtJets::KtJets(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState::LossyFinalState(), MergedFinalState::MergedFinalState(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerUA5::TriggerUA5(), and VetoedFinalState::VetoedFinalState().

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

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

Untemplated function to do the work...

Definition at line 33 of file ProjectionApplier.cc.

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

Referenced by ProjectionApplier::addProjection().

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

Friends And Related Function Documentation

friend class Projectionhandler [friend, inherited]

Definition at line 23 of file ProjectionApplier.hh.

Member Data Documentation

AIDA::IProfile1D* _pt90MaxAvg1800 [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the average $ p_T $ in the toward, transverse and away regions at $ \sqrt{s} = 1800 \text{GeV} $. Corresponds to Table 1, and HepData table 1.

Definition at line 332 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pt90MinAvg1800 [private]

Definition at line 332 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D* _pt90Max1800 [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the $ p_T $ sum in the toward, transverse and away regions at $ \sqrt{s} = 1800 \text{GeV} $. Corresponds to figure 2/3, and HepData table 2.

Definition at line 338 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pt90Min1800 [private]

Definition at line 338 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pt90Diff1800 [private]

Definition at line 338 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D* _pt90Max630 [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the $ p_T $ sum in the toward, transverse and away regions at at $ \sqrt{s} = 630 \text{GeV} $. Corresponds to figure 8, and HepData table 8.

Definition at line 344 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pt90Min630 [private]

Definition at line 344 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pt90Diff630 [private]

Definition at line 344 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D* _num90Max1800 [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the cone track multiplicity at $ \sqrt{s} = 1800 \text{GeV} $. Corresponds to figure 5, and HepData table 4.

Definition at line 349 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _num90Min1800 [private]

Definition at line 349 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D* _pTSum1800_2Jet [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the $ p_T $ sum at $ \sqrt{s} = 1800 \text{GeV} $. Corresponds to figure 7, and HepData table 7.

Definition at line 354 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pTSum1800_3Jet [private]

Definition at line 354 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D* _pTSum630_2Jet [private]

Profile histograms, binned in the $ E_T $ of the leading jet, for the $ p_T $ sum at $ \sqrt{s} = 630 \text{GeV} $. Corresponds to figure 9, and HepData table 9.

Definition at line 359 of file CDF_2004_S5839831.cc.

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

AIDA::IProfile1D * _pTSum630_3Jet [private]

Definition at line 359 of file CDF_2004_S5839831.cc.

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

AIDA::IHistogram1D* _pt90Dbn1800Et40 [private]

Histogram of $ p_{T\text{sum}} $ distribution for 5 different $ E_{T1} $ bins. Corresponds to figure 4, and HepData table 3.

Definition at line 364 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _pt90Dbn1800Et80 [private]

Definition at line 364 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _pt90Dbn1800Et120 [private]

Definition at line 364 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _pt90Dbn1800Et160 [private]

Definition at line 364 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _pt90Dbn1800Et200 [private]

Definition at line 364 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D* _numTracksDbn1800MB [private]

Histograms of track multiplicity and $ p_T $ distributions for minimum bias events. Figure 6, and HepData tables 5 & 6. Figure 10, and HepData tables 10 & 11.

Definition at line 371 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _ptDbn1800MB [private]

Definition at line 371 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D* _numTracksDbn630MB [private]

Definition at line 372 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

AIDA::IHistogram1D * _ptDbn630MB [private]

Definition at line 372 of file CDF_2004_S5839831.cc.

Referenced by CDF_2004_S5839831::analyze(), CDF_2004_S5839831::finalize(), and CDF_2004_S5839831::init().

string _defaultname [protected, inherited]

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

Definition at line 412 of file Analysis.hh.

Referenced by Analysis::name().

shared_ptr<AnalysisInfo> _info [protected, inherited]

Pointer to analysis metadata object.

Definition at line 415 of file Analysis.hh.

Referenced by Analysis::authors(), Analysis::collider(), Analysis::description(), Analysis::experiment(), Analysis::info(), Analysis::name(), Analysis::references(), Analysis::runInfo(), Analysis::setBeams(), Analysis::spiresId(), Analysis::status(), Analysis::summary(), and Analysis::year().

bool _allowProjReg [protected, inherited]

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

Definition at line 141 of file ProjectionApplier.hh.

Referenced by ProjectionApplier::_addProjection(), and AnalysisHandler::init().

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