MC_PHOTONJETUE Class Reference
MC validation analysis for underlying event in jet + isolated photon events. More...
Public Member Functions | |
| MC_PHOTONJETUE () | |
| Constructor. | |
| AnalysisHandler & | handler () const |
| Access the controlling AnalysisHandler object. | |
| void | normalize (AIDA::IHistogram1D *&histo, double norm=1.0) |
| void | normalize (AIDA::IHistogram2D *&histo, double norm=1.0) |
| void | scale (AIDA::IHistogram1D *&histo, double scale) |
| void | scale (AIDA::IHistogram2D *&histo, double scale) |
| Analysis & | setCrossSection (double xs) |
| Set the cross section from the generator. | |
Metadata | |
Metadata is used for querying from the command line and also for building web pages and the analysis pages in the Rivet manual. | |
| const 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 std::string | year () const |
| When the original experimental analysis was published. | |
| virtual std::vector< std::string > | references () const |
| Journal, and preprint references. | |
| virtual std::string | bibKey () const |
| BibTeX citation key for this article. | |
| virtual std::string | bibTeX () const |
| BibTeX citation entry for this article. | |
| virtual std::string | status () const |
| Whether this analysis is trusted (in any way!). | |
| virtual std::vector< std::string > | todos () const |
| Any work to be done on this analysis. | |
| virtual const std::vector < PdgIdPair > & | requiredBeams () const |
| Return the allowed pairs of incoming beams required by this analysis. | |
| virtual Analysis & | setRequiredBeams (const std::vector< PdgIdPair > &requiredBeams) |
| Declare the allowed pairs of incoming beams required by this analysis. | |
| virtual const std::vector < std::pair< double, double > > & | requiredEnergies () const |
| Sets of valid beam energy pairs, in GeV. | |
| virtual Analysis & | setRequiredEnergies (const std::vector< std::pair< double, double > > &requiredEnergies) |
| Declare the list of valid beam energy pairs, in GeV. | |
| bool | needsCrossSection () const |
| Return true if this analysis needs to know the process cross-section. | |
| Analysis & | setNeedsCrossSection (bool needed=true) |
| Declare whether this analysis needs to know the process cross-section from the generator. | |
Internal metadata modifiying methods | |
| AnalysisInfo & | info () |
| Get the actual AnalysisInfo object in which all this metadata is stored (non-const). | |
| virtual Analysis & | setBeams (PdgId beam1, PdgId beam2) |
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. | |
Analysis / beam compatibility testing | |
| bool | isCompatible (const ParticlePair &beams) const |
| Check if analysis is compatible with the provided beam particle IDs and energies. | |
| bool | isCompatible (PdgId beam1, PdgId beam2, double e1, double e2) const |
| Check if analysis is compatible with the provided beam particle IDs and energies. | |
| bool | isCompatible (const PdgIdPair &beams, const std::pair< double, double > &energies) const |
| Check if analysis is compatible with the provided beam particle IDs and energies. | |
Projection "getting" functions | |
| std::set< 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 |
| 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="") |
| AIDA::IHistogram2D * | bookHistogram2D (const std::string &name, size_t nxbins, double xlower, double xupper, size_t nybins, double ylower, double yupper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="", const std::string &ztitle="") |
| AIDA::IHistogram2D * | bookHistogram2D (const std::string &name, const std::vector< double > &xbinedges, const std::vector< double > &ybinedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="", const std::string &ztitle="") |
Internal profile histogram booking (for use by Analysis sub-classes). | |
| AIDA::IProfile1D * | bookProfile1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
| AIDA::IProfile1D * | bookProfile1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
| AIDA::IProfile1D * | bookProfile1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
| AIDA::IProfile1D * | bookProfile1D (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
Internal data point set booking (for use by Analysis sub-classes). | |
| AIDA::IDataPointSet * | bookDataPointSet (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
| AIDA::IDataPointSet * | bookDataPointSet (const std::string &name, size_t npts, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
| AIDA::IDataPointSet * | bookDataPointSet (size_t datasetId, size_t xAxisId, size_t yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="") |
Projection registration functions | |
| template<typename PROJ > | |
| const PROJ & | addProjection (const PROJ &proj, const std::string &name) |
| const Projection & | _addProjection (const Projection &proj, const std::string &name) |
| Untemplated function to do the work... | |
Protected Attributes | |
| string | _defaultname |
| Name passed to constructor (used to find .info analysis data file, and as a fallback). | |
| shared_ptr< AnalysisInfo > | _info |
| Pointer to analysis metadata object. | |
| bool | _allowProjReg |
| Flag to forbid projection registration in analyses until the init phase. | |
Analysis methods | |
|
| |
| AIDA::IHistogram1D * | _hist_jetgamma_dR |
| AIDA::IHistogram1D * | _hist_jetgamma_dphi |
| AIDA::IProfile1D * | _hist_pnchg_jet |
| AIDA::IProfile1D * | _hist_pnchg_gamma |
| AIDA::IProfile1D * | _hist_pmaxnchg_jet |
| AIDA::IProfile1D * | _hist_pmaxnchg_gamma |
| AIDA::IProfile1D * | _hist_pminnchg_jet |
| AIDA::IProfile1D * | _hist_pminnchg_gamma |
| AIDA::IProfile1D * | _hist_pcptsum_jet |
| AIDA::IProfile1D * | _hist_pcptsum_gamma |
| AIDA::IProfile1D * | _hist_pmaxcptsum_jet |
| AIDA::IProfile1D * | _hist_pmaxcptsum_gamma |
| AIDA::IProfile1D * | _hist_pmincptsum_jet |
| AIDA::IProfile1D * | _hist_pmincptsum_gamma |
| AIDA::IProfile1D * | _hist_pcptave_jet |
| AIDA::IProfile1D * | _hist_pcptave_gamma |
| void | init () |
| void | analyze (const Event &evt) |
| void | finalize () |
Detailed Description
MC validation analysis for underlying event in jet + isolated photon events.
Definition at line 14 of file MC_PHOTONJETUE.cc.
Constructor & Destructor Documentation
| MC_PHOTONJETUE | ( | ) | [inline] |
Constructor.
Definition at line 18 of file MC_PHOTONJETUE.cc.
00019 : Analysis("MC_PHOTONJETUE") 00020 { }
Member Function Documentation
| const Projection & _addProjection | ( | const Projection & | proj, | |
| const std::string & | name | |||
| ) | [protected, inherited] |
Untemplated function to do the work...
Definition at line 33 of file ProjectionApplier.cc.
References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().
Referenced by ProjectionApplier::addProjection().
00034 { 00035 if (!_allowProjReg) { 00036 cerr << "Trying to register projection '" 00037 << proj.name() << "' before init phase in '" << this->name() << "'." << endl; 00038 exit(2); 00039 } 00040 const Projection& reg = getProjHandler().registerProjection(*this, proj, name); 00041 return reg; 00042 }
| const PROJ& addProjection | ( | const PROJ & | proj, | |
| const std::string & | name | |||
| ) | [inline, protected, inherited] |
Register a contained projection. The type of the argument is used to instantiate a new projection internally: this new object is applied to events rather than the argument object. Hence you are advised to only use locally-scoped Projection objects in your Projection and Analysis constructors, and to avoid polymorphism (e.g. handling ConcreteProjection via a pointer or reference to type Projection) since this will screw up the internal type management.
Definition at line 113 of file ProjectionApplier.hh.
References ProjectionApplier::_addProjection().
Referenced by ZFinder::_init(), WFinder::_init(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ClusteredPhotons::ClusteredPhotons(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), FinalStateHCM::FinalStateHCM(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), OPAL_1993_S2692198::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_LEADINGJETS::init(), MC_JETS::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), LHCB_2010_S8758301::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2003_WUD_03_11::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), DELPHI_1995_S3137023::init(), D0_2010_S8821313::init(), D0_2010_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CMS_2011_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::init(), CMS_2011_S8884919::init(), CMS_2010_S8656010::init(), CMS_2010_S8547297::init(), CDF_2010_S8591881_QCD::init(), CDF_2010_S8591881_DY::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6217184::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2011_S9120807::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S9002537::init(), ATLAS_2011_S8994773::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_S8971293::init(), ATLAS_2011_S8924791::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2010_S8919674::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8914702::init(), ATLAS_2010_S8894728::init(), ATLAS_2010_S8817804::init(), ATLAS_2010_S8591806::init(), ATLAS_2010_CONF_2010_049::init(), ALICE_2010_S8706239::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), ALEPH_1991_S2435284::init(), IsolationProjection< PROJ1, PROJ2, EST >::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LeptonClusters::LeptonClusters(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), Sphericity::Sphericity(), Spherocity::Spherocity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), and VisibleFinalState::VisibleFinalState().
00113 { 00114 const Projection& reg = _addProjection(proj, name); 00115 const PROJ& rtn = dynamic_cast<const PROJ&>(reg); 00116 return rtn; 00117 }
| IAnalysisFactory & analysisFactory | ( | ) | [protected, inherited] |
Access the AIDA analysis factory of the controlling AnalysisHandler object.
Definition at line 50 of file Analysis.cc.
References AnalysisHandler::analysisFactory(), and Analysis::handler().
00050 { 00051 return handler().analysisFactory(); 00052 }
| void analyze | ( | const Event & | event | ) | [inline, virtual] |
Analyze one event. A concrete class should here apply the necessary projections on the event and fill the relevant histograms. An overridden function must make sure it first calls the base class function.
Cut on back-to-backness of jet-photon
- Todo:
- Plot evolution of UE as a function of jet-photon angle
- Todo:
- Plot evolution of UE as a function of photon pT
Implements Analysis.
Definition at line 67 of file MC_PHOTONJETUE.cc.
References MC_PHOTONJETUE::_hist_jetgamma_dphi, MC_PHOTONJETUE::_hist_jetgamma_dR, MC_PHOTONJETUE::_hist_pcptave_gamma, MC_PHOTONJETUE::_hist_pcptave_jet, MC_PHOTONJETUE::_hist_pcptsum_gamma, MC_PHOTONJETUE::_hist_pcptsum_jet, MC_PHOTONJETUE::_hist_pmaxcptsum_gamma, MC_PHOTONJETUE::_hist_pmaxcptsum_jet, MC_PHOTONJETUE::_hist_pmaxnchg_gamma, MC_PHOTONJETUE::_hist_pmaxnchg_jet, MC_PHOTONJETUE::_hist_pmincptsum_gamma, MC_PHOTONJETUE::_hist_pmincptsum_jet, MC_PHOTONJETUE::_hist_pminnchg_gamma, MC_PHOTONJETUE::_hist_pminnchg_jet, MC_PHOTONJETUE::_hist_pnchg_gamma, MC_PHOTONJETUE::_hist_pnchg_jet, Log::DEBUG, Rivet::deltaPhi(), Rivet::deltaR(), FourVector::eta(), Analysis::getLog(), Rivet::GeV, Rivet::MINUS, Particle::momentum(), Jet::momentum(), FinalState::particles(), FinalState::particlesByPt(), FourVector::phi(), Rivet::PI, Rivet::pT(), FourMomentum::pT(), Rivet::sign(), FinalState::size(), and vetoEvent.
00067 { 00068 00069 // Get jets 00070 const Jets jets = applyProjection<FastJets>(evt, "Jets").jetsByPt(); 00071 getLog() << Log::DEBUG << "Jet multiplicity = " << jets.size() << endl; 00072 if (jets.size() < 1) { 00073 getLog() << Log::DEBUG << "No jets found" << endl; 00074 vetoEvent; 00075 } 00076 00077 // Get leading jet properties 00078 const FourMomentum pjet = jets.front().momentum(); 00079 const double jeteta = pjet.eta(); 00080 const double jetphi = pjet.phi(); 00081 const double jetpT = pjet.pT(); 00082 getLog() << Log::DEBUG << "Leading jet: pT = " << jetpT/GeV << " GeV" 00083 << ", eta = " << jeteta 00084 << ", phi = " << jetphi << endl; 00085 00086 // Require the leading jet to be within |eta| < 2 00087 if (fabs(jeteta) > 2) { 00088 getLog() << Log::DEBUG << "Failed leading jet eta cut" << endl; 00089 vetoEvent; 00090 } 00091 00092 // Get the leading photon 00093 const FinalState& photonfs = applyProjection<FinalState>(evt, "Photons"); 00094 if (photonfs.size() < 1) { 00095 getLog() << Log::DEBUG << "No hard photons found" << endl; 00096 vetoEvent; 00097 } 00098 const FourMomentum pgamma = photonfs.particlesByPt().front().momentum(); 00099 00100 // Check that leading photon is isolated from jets 00101 bool isolated = true; 00102 foreach (const Jet& j, jets) { 00103 if (deltaR(j.momentum(), pgamma) < 0.1) { 00104 isolated = false; 00105 break; 00106 } 00107 } 00108 if (!isolated) { 00109 getLog() << Log::DEBUG << "Leading photon is not isolated from jets" << endl; 00110 vetoEvent; 00111 } 00112 00113 // Get leading photon properties 00114 const double gammaeta = pgamma.eta(); 00115 const double gammaphi = pgamma.phi(); 00116 const double gammapT = pgamma.pT(); 00117 getLog() << Log::DEBUG << "Leading photon: pT = " << gammapT/GeV << " GeV" 00118 << ", eta = " << gammaeta 00119 << ", phi = " << gammaphi << endl; 00120 00121 // Get the event weight 00122 const double weight = evt.weight(); 00123 00124 // Fill jet1-photon separation histos 00125 const double dR_jetgamma = deltaR(pgamma, pjet); 00126 _hist_jetgamma_dR->fill(dR_jetgamma, weight); 00127 const double dPhi_jetgamma = deltaPhi(gammaphi, jetphi); 00128 _hist_jetgamma_dphi->fill(dPhi_jetgamma, weight); 00129 00130 /// Cut on back-to-backness of jet-photon 00131 if (dPhi_jetgamma < 0.5) vetoEvent; 00132 00133 /// @todo Plot evolution of UE as a function of jet-photon angle 00134 /// @todo Plot evolution of UE as a function of photon pT 00135 00136 // Get the final states to work with for filling the distributions 00137 const FinalState& cfs = applyProjection<ChargedFinalState>(evt, "Tracks"); 00138 00139 // Whole-event counters 00140 unsigned int numOverall(0); 00141 double ptSumOverall(0.0), ptMaxOverall(0.0); 00142 // Jet-oriented counters 00143 unsigned int numToward_jet(0), numTrans1_jet(0), numTrans2_jet(0), numAway_jet(0); 00144 double ptSumToward_jet(0.0), ptSumTrans1_jet(0.0), ptSumTrans2_jet(0.0), ptSumAway_jet(0.0); 00145 double ptMaxToward_jet(0.0), ptMaxTrans1_jet(0.0), ptMaxTrans2_jet(0.0), ptMaxAway_jet(0.0); 00146 // Photon-oriented counters 00147 unsigned int numToward_gamma(0), numTrans1_gamma(0), numTrans2_gamma(0), numAway_gamma(0); 00148 double ptSumToward_gamma(0.0), ptSumTrans1_gamma(0.0), ptSumTrans2_gamma(0.0), ptSumAway_gamma(0.0); 00149 double ptMaxToward_gamma(0.0), ptMaxTrans1_gamma(0.0), ptMaxTrans2_gamma(0.0), ptMaxAway_gamma(0.0); 00150 00151 // Calculate all the charged stuff 00152 foreach (const Particle& p, cfs.particles()) { 00153 ++numOverall; 00154 const double pT = p.momentum().pT(); 00155 ptSumOverall += pT; 00156 if (pT > ptMaxOverall) ptMaxOverall = pT; 00157 00158 // Increment jet-oriented variables 00159 const double dPhi_jet = jetphi - p.momentum().phi(); 00160 if (fabs(dPhi_jet) < PI/3.0) { 00161 ptSumToward_jet += pT; 00162 ++numToward_jet; 00163 if (pT > ptMaxToward_jet) ptMaxToward_jet = pT; 00164 } 00165 else if (fabs(dPhi_jet) < 2*PI/3.0) { 00166 if (sign(dPhi_jet) == MINUS) { 00167 ptSumTrans1_jet += pT; 00168 ++numTrans1_jet; 00169 if (pT > ptMaxTrans1_jet) ptMaxTrans1_jet = pT; 00170 } else { 00171 ptSumTrans2_jet += pT; 00172 ++numTrans2_jet; 00173 if (pT > ptMaxTrans2_jet) ptMaxTrans2_jet = pT; 00174 } 00175 } 00176 else { 00177 ptSumAway_jet += pT; 00178 ++numAway_jet; 00179 if (pT > ptMaxAway_jet) ptMaxAway_jet = pT; 00180 } 00181 00182 00183 // Increment photon-oriented variables 00184 const double dPhi_gamma = gammaphi - p.momentum().phi(); 00185 if (fabs(dPhi_gamma) < PI/3.0) { 00186 ptSumToward_gamma += pT; 00187 ++numToward_gamma; 00188 if (pT > ptMaxToward_gamma) ptMaxToward_gamma = pT; 00189 } 00190 else if (fabs(dPhi_gamma) < 2*PI/3.0) { 00191 if (sign(dPhi_gamma) == MINUS) { 00192 ptSumTrans1_gamma += pT; 00193 ++numTrans1_gamma; 00194 if (pT > ptMaxTrans1_gamma) ptMaxTrans1_gamma = pT; 00195 } else { 00196 ptSumTrans2_gamma += pT; 00197 ++numTrans2_gamma; 00198 if (pT > ptMaxTrans2_gamma) ptMaxTrans2_gamma = pT; 00199 } 00200 } 00201 else { 00202 ptSumAway_gamma += pT; 00203 ++numAway_gamma; 00204 if (pT > ptMaxAway_gamma) ptMaxAway_gamma = pT; 00205 } 00206 00207 00208 } 00209 00210 00211 // Fill the histograms 00212 _hist_pnchg_jet->fill(jetpT/GeV, (numTrans1_jet+numTrans2_jet)/(4*PI/3), weight); 00213 _hist_pmaxnchg_jet->fill(jetpT/GeV, (numTrans1_jet>numTrans2_jet ? numTrans1_jet : numTrans2_jet)/(2*PI/3), weight); 00214 _hist_pminnchg_jet->fill(jetpT/GeV, (numTrans1_jet<numTrans2_jet ? numTrans1_jet : numTrans2_jet)/(2*PI/3), weight); 00215 _hist_pcptsum_jet->fill(jetpT/GeV, (ptSumTrans1_jet+ptSumTrans2_jet)/GeV/(4*PI/3), weight); 00216 _hist_pmaxcptsum_jet->fill(jetpT/GeV, (ptSumTrans1_jet>ptSumTrans2_jet ? ptSumTrans1_jet : ptSumTrans2_jet)/GeV/(2*PI/3), weight); 00217 _hist_pmincptsum_jet->fill(jetpT/GeV, (ptSumTrans1_jet<ptSumTrans2_jet ? ptSumTrans1_jet : ptSumTrans2_jet)/GeV/(2*PI/3), weight); 00218 if ((numTrans1_jet+numTrans2_jet) > 0) { 00219 _hist_pcptave_jet->fill(jetpT/GeV, (ptSumTrans1_jet+ptSumTrans2_jet)/GeV/(numTrans1_jet+numTrans2_jet), weight); 00220 } 00221 // 00222 _hist_pnchg_gamma->fill(gammapT/GeV, (numTrans1_gamma+numTrans2_gamma)/(4*PI/3), weight); 00223 _hist_pmaxnchg_gamma->fill(gammapT/GeV, (numTrans1_gamma>numTrans2_gamma ? numTrans1_gamma : numTrans2_gamma)/(2*PI/3), weight); 00224 _hist_pminnchg_gamma->fill(gammapT/GeV, (numTrans1_gamma<numTrans2_gamma ? numTrans1_gamma : numTrans2_gamma)/(2*PI/3), weight); 00225 _hist_pcptsum_gamma->fill(gammapT/GeV, (ptSumTrans1_gamma+ptSumTrans2_gamma)/GeV/(4*PI/3), weight); 00226 _hist_pmaxcptsum_gamma->fill(gammapT/GeV, (ptSumTrans1_gamma>ptSumTrans2_gamma ? ptSumTrans1_gamma : ptSumTrans2_gamma)/GeV/(2*PI/3), weight); 00227 _hist_pmincptsum_gamma->fill(gammapT/GeV, (ptSumTrans1_gamma<ptSumTrans2_gamma ? ptSumTrans1_gamma : ptSumTrans2_gamma)/GeV/(2*PI/3), weight); 00228 if ((numTrans1_gamma+numTrans2_gamma) > 0) { 00229 _hist_pcptave_gamma->fill(gammapT/GeV, (ptSumTrans1_gamma+ptSumTrans2_gamma)/GeV/(numTrans1_gamma+numTrans2_gamma), weight); 00230 } 00231 00232 }
| const PROJ& applyProjection | ( | const Event & | evt, | |
| const std::string & | name | |||
| ) | const [inline, inherited] |
Apply the named projection on event.
Definition at line 81 of file ProjectionApplier.hh.
References ProjectionApplier::_applyProjection().
00081 { 00082 return pcast<PROJ>(_applyProjection(evt, name)); 00083 }
| const PROJ& applyProjection | ( | const Event & | evt, | |
| const Projection & | proj | |||
| ) | const [inline, inherited] |
Apply the supplied projection on event.
Definition at line 74 of file ProjectionApplier.hh.
References ProjectionApplier::_applyProjection().
00074 { 00075 return pcast<PROJ>(_applyProjection(evt, proj)); 00076 }
| const PROJ& applyProjection | ( | const Event & | evt, | |
| const PROJ & | proj | |||
| ) | const [inline, inherited] |
Apply the supplied projection on event.
Definition at line 67 of file ProjectionApplier.hh.
References ProjectionApplier::_applyProjection().
Referenced by HadronicFinalState::project(), and FinalStateHCM::project().
00067 { 00068 return pcast<PROJ>(_applyProjection(evt, proj)); 00069 }
| virtual std::vector<std::string> authors | ( | ) | const [inline, virtual, inherited] |
Names & emails of paper/analysis authors.
Names and email of authors in 'NAME <EMAIL>' format. The first name in the list should be the primary contact person.
Definition at line 125 of file Analysis.hh.
References AnalysisInfo::authors(), and Analysis::info().
00125 { 00126 return info().authors(); 00127 }
| const PdgIdPair beamIds | ( | ) | const [inherited] |
Incoming beam IDs for this run.
Definition at line 78 of file Analysis.cc.
References AnalysisHandler::beamIds(), and Analysis::handler().
Referenced by UA5_1982_S875503::finalize(), and UA5_1982_S875503::init().
00078 { 00079 return handler().beamIds(); 00080 }
| const ParticlePair & beams | ( | ) | const [inherited] |
Incoming beams for this run.
Definition at line 74 of file Analysis.cc.
References AnalysisHandler::beams(), and Analysis::handler().
Referenced by OPAL_1998_S3780481::analyze(), MC_WPOL::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1996_S3430090::analyze(), DELPHI_1995_S3137023::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_2004_S5765862::analyze(), ALEPH_1996_S3486095::analyze(), MC_WPOL::init(), and Analysis::isCompatible().
00074 { 00075 return handler().beams(); 00076 }
| virtual std::string bibKey | ( | ) | const [inline, virtual, inherited] |
BibTeX citation key for this article.
Definition at line 178 of file Analysis.hh.
References AnalysisInfo::bibKey(), and Analysis::info().
00178 { 00179 return info().bibKey(); 00180 }
| virtual std::string bibTeX | ( | ) | const [inline, virtual, inherited] |
BibTeX citation entry for this article.
Definition at line 183 of file Analysis.hh.
References AnalysisInfo::bibTeX(), and Analysis::info().
00183 { 00184 return info().bibTeX(); 00185 }
| const BinEdges & binEdges | ( | size_t | datasetId, | |
| size_t | xAxisId, | |||
| size_t | yAxisId | |||
| ) | const [protected, inherited] |
Get bin edges for a numbered histo (using ref AIDA caching).
Definition at line 223 of file Analysis.cc.
References Analysis::binEdges().
00223 { 00224 const string hname = makeAxisCode(datasetId, xAxisId, yAxisId); 00225 return binEdges(hname); 00226 }
| const BinEdges & binEdges | ( | const std::string & | hname | ) | const [protected, inherited] |
Get bin edges for a named histo (using ref AIDA caching).
Definition at line 208 of file Analysis.cc.
References Analysis::_cacheBinEdges(), Analysis::_histBinEdges, Analysis::getLog(), MSG_TRACE, Analysis::name(), and Log::TRACE.
Referenced by ATLAS_2011_S8994773::analyze(), ATLAS_2010_S8894728::analyze(), Analysis::binEdges(), Analysis::bookHistogram1D(), Analysis::bookProfile1D(), D0_2008_S7837160::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9002537::init(), and ATLAS_2010_S8894728::init().
00208 { 00209 _cacheBinEdges(); 00210 MSG_TRACE("Using histo bin edges for " << name() << ":" << hname); 00211 const BinEdges& edges = _histBinEdges.find(hname)->second; 00212 if (getLog().isActive(Log::TRACE)) { 00213 stringstream edges_ss; 00214 foreach (const double be, edges) { 00215 edges_ss << " " << be; 00216 } 00217 MSG_TRACE("Edges:" << edges_ss.str()); 00218 } 00219 return edges; 00220 }
| IDataPointSet * bookDataPointSet | ( | size_t | datasetId, | |
| size_t | xAxisId, | |||
| size_t | yAxisId, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 2-dimensional data point set based on the corresponding AIDA data file. The binnings (x-errors) will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property. Book a 2-dimensional data point set based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings (x-errors) will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.
Definition at line 418 of file Analysis.cc.
References Analysis::_cacheXAxisData(), Analysis::_dpsData, Analysis::bookDataPointSet(), MSG_TRACE, and Analysis::name().
00420 { 00421 // Get the bin edges (only read the AIDA file once) 00422 _cacheXAxisData(); 00423 // Build the axis code 00424 const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); 00425 //const map<string, vector<DPSXPoint> > xpoints = getDPSXValsErrs(papername); 00426 MSG_TRACE("Using DPS x-positions for " << name() << ":" << axisCode); 00427 IDataPointSet* dps = bookDataPointSet(axisCode, title, xtitle, ytitle); 00428 const vector<DPSXPoint> xpts = _dpsData.find(axisCode)->second; 00429 for (size_t pt = 0; pt < xpts.size(); ++pt) { 00430 dps->addPoint(); 00431 IMeasurement* meas = dps->point(pt)->coordinate(0); 00432 meas->setValue(xpts[pt].val); 00433 meas->setErrorPlus(xpts[pt].errplus); 00434 meas->setErrorMinus(xpts[pt].errminus); 00435 } 00436 MSG_TRACE("Made DPS " << axisCode << " for " << name()); 00437 return dps; 00438 }
| IDataPointSet * bookDataPointSet | ( | const std::string & | name, | |
| size_t | npts, | |||
| double | lower, | |||
| double | upper, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 2-dimensional data point set with equally spaced points in a range. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 400 of file Analysis.cc.
References Analysis::bookDataPointSet().
00403 { 00404 IDataPointSet* dps = bookDataPointSet(hname, title, xtitle, ytitle); 00405 for (size_t pt = 0; pt < npts; ++pt) { 00406 const double binwidth = (upper-lower)/npts; 00407 const double bincentre = lower + (pt + 0.5) * binwidth; 00408 dps->addPoint(); 00409 IMeasurement* meas = dps->point(pt)->coordinate(0); 00410 meas->setValue(bincentre); 00411 meas->setErrorPlus(binwidth/2.0); 00412 meas->setErrorMinus(binwidth/2.0); 00413 } 00414 return dps; 00415 }
| IDataPointSet * bookDataPointSet | ( | const std::string & | name, | |
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 2-dimensional data point set. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 388 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::datapointsetFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
Referenced by Analysis::bookDataPointSet(), ALEPH_2004_S5765862::finalize(), UA5_1988_S1867512::init(), STAR_2006_S6860818::init(), OPAL_1993_S2692198::init(), MC_XS::init(), MC_JetAnalysis::init(), JADE_OPAL_2000_S4300807::init(), D0_2001_S4674421::init(), CDF_2008_S7782535::init(), CDF_2008_S7541902::init(), CDF_2005_S6217184::init(), CDF_1996_S3418421::init(), CDF_1994_S2952106::init(), ATLAS_2011_S9002537::init(), ATLAS_2010_S8894728::init(), and ALEPH_2004_S5765862::init().
00389 { 00390 _makeHistoDir(); 00391 const string path = histoPath(hname); 00392 IDataPointSet* dps = datapointsetFactory().create(path, title, 2); 00393 MSG_TRACE("Made data point set " << hname << " for " << name()); 00394 dps->setXTitle(xtitle); 00395 dps->setYTitle(ytitle); 00396 return dps; 00397 }
| IHistogram1D * bookHistogram1D | ( | size_t | datasetId, | |
| size_t | xAxisId, | |||
| size_t | yAxisId, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D histogram based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.
Definition at line 242 of file Analysis.cc.
References Analysis::bookHistogram1D().
00245 { 00246 const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); 00247 return bookHistogram1D(axisCode, title, xtitle, ytitle); 00248 }
| IHistogram1D * bookHistogram1D | ( | const std::string & | name, | |
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D histogram based on the name in the corresponding AIDA file. The binnings will be obtained by reading the bundled AIDA data record file with the same filename as the analysis' name() property.
Definition at line 251 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00253 { 00254 // Get the bin edges (only read the AIDA file once) 00255 const BinEdges edges = binEdges(hname); 00256 _makeHistoDir(); 00257 const string path = histoPath(hname); 00258 IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, edges); 00259 MSG_TRACE("Made histogram " << hname << " for " << name()); 00260 hist->setXTitle(xtitle); 00261 hist->setYTitle(ytitle); 00262 return hist; 00263 }
| IHistogram1D * bookHistogram1D | ( | const std::string & | name, | |
| const std::vector< double > & | binedges, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D histogram with non-uniform bins defined by the vector of bin edges binedges . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 280 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00283 { 00284 _makeHistoDir(); 00285 const string path = histoPath(hname); 00286 IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, binedges); 00287 MSG_TRACE("Made histogram " << hname << " for " << name()); 00288 hist->setXTitle(xtitle); 00289 hist->setYTitle(ytitle); 00290 return hist; 00291 }
| IHistogram1D * bookHistogram1D | ( | const std::string & | name, | |
| size_t | nbins, | |||
| double | lower, | |||
| double | upper, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D histogram with nbins uniformly distributed across the range lower - upper . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 266 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
Referenced by Analysis::bookHistogram1D(), ZEUS_2001_S4815815::init(), UA5_1989_S1926373::init(), UA5_1987_S1640666::init(), UA5_1986_S1583476::init(), UA5_1982_S875503::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), STAR_2008_S7869363::init(), STAR_2006_S6870392::init(), STAR_2006_S6860818::init(), STAR_2006_S6500200::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), OPAL_2001_S4553896::init(), OPAL_1998_S3780481::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_TTBAR::init(), MC_SUSY::init(), MC_PHOTONJETUE::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), MC_GENERIC::init(), MC_DIPHOTON::init(), MC_DIJET::init(), LHCB_2010_S8758301::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), H1_2000_S4129130::init(), H1_1995_S3167097::init(), H1_1994_S2919893::init(), ExampleAnalysis::init(), E735_1998_S3905616::init(), DELPHI_2003_WUD_03_11::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), DELPHI_1995_S3137023::init(), D0_2010_S8821313::init(), D0_2010_S8671338::init(), D0_2010_S8570965::init(), D0_2010_S8566488::init(), D0_2009_S8349509::init(), D0_2009_S8320160::init(), D0_2009_S8202443::init(), D0_2008_S7863608::init(), D0_2008_S7837160::init(), D0_2008_S7719523::init(), D0_2008_S7662670::init(), D0_2008_S7554427::init(), D0_2008_S6879055::init(), D0_2007_S7075677::init(), D0_2006_S6438750::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), D0_2000_S4480767::init(), D0_1996_S3324664::init(), D0_1996_S3214044::init(), CMS_2011_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::init(), CMS_2011_S8884919::init(), CMS_2010_S8656010::init(), CMS_2010_S8547297::init(), CDF_2009_S8436959::init(), CDF_2009_S8383952::init(), CDF_2009_S8233977::init(), CDF_2009_NOTE_9936::init(), CDF_2008_S8095620::init(), CDF_2008_S8093652::init(), CDF_2008_S7828950::init(), CDF_2008_S7541902::init(), CDF_2008_S7540469::init(), CDF_2007_S7057202::init(), CDF_2006_S6653332::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), CDF_2000_S4266730::init(), CDF_2000_S4155203::init(), CDF_1998_S3618439::init(), CDF_1997_S3541940::init(), CDF_1996_S3418421::init(), CDF_1996_S3349578::init(), CDF_1996_S3108457::init(), CDF_1994_S2952106::init(), CDF_1993_S2742446::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), BELLE_2006_S6265367::init(), ATLAS_2011_S9120807::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S9002537::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_S8971293::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2010_S8919674::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8914702::init(), ATLAS_2010_S8817804::init(), ATLAS_2010_S8591806::init(), ATLAS_2010_CONF_2010_049::init(), ALICE_2010_S8706239::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), ALEPH_2004_S5765862::init(), ALEPH_1996_S3486095::init(), ALEPH_1996_S3196992::init(), and ALEPH_1991_S2435284::init().
00269 { 00270 _makeHistoDir(); 00271 const string path = histoPath(hname); 00272 IHistogram1D* hist = histogramFactory().createHistogram1D(path, title, nbins, lower, upper); 00273 MSG_TRACE("Made histogram " << hname << " for " << name()); 00274 hist->setXTitle(xtitle); 00275 hist->setYTitle(ytitle); 00276 return hist; 00277 }
| IHistogram2D * bookHistogram2D | ( | const std::string & | name, | |
| const std::vector< double > & | xbinedges, | |||
| const std::vector< double > & | ybinedges, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "", |
|||
| const std::string & | ztitle = "" | |||
| ) | [protected, inherited] |
Book a 2D histogram with non-uniform bins defined by the vectorx of bin edges xbinedges and ybinedges. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 313 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00317 { 00318 _makeHistoDir(); 00319 const string path = histoPath(hname); 00320 IHistogram2D* hist = 00321 histogramFactory().createHistogram2D(path, title, xbinedges, ybinedges); 00322 MSG_TRACE("Made 2D histogram " << hname << " for " << name()); 00323 hist->setXTitle(xtitle); 00324 hist->setYTitle(ytitle); 00325 hist->setZTitle(ztitle); 00326 return hist; 00327 }
| IHistogram2D * bookHistogram2D | ( | const std::string & | name, | |
| size_t | nxbins, | |||
| double | xlower, | |||
| double | xupper, | |||
| size_t | nybins, | |||
| double | ylower, | |||
| double | yupper, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "", |
|||
| const std::string & | ztitle = "" | |||
| ) | [protected, inherited] |
Book a 2D histogram with nxbins and nybins uniformly distributed across the ranges xlower - xupper and ylower - yupper respectively along the x- and y-axis. (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 294 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00298 { 00299 _makeHistoDir(); 00300 const string path = histoPath(hname); 00301 IHistogram2D* hist = 00302 histogramFactory().createHistogram2D(path, title, nxbins, xlower, xupper, 00303 nybins, ylower, yupper); 00304 MSG_TRACE("Made 2D histogram " << hname << " for " << name()); 00305 hist->setXTitle(xtitle); 00306 hist->setYTitle(ytitle); 00307 hist->setZTitle(ztitle); 00308 return hist; 00309 }
| IProfile1D * bookProfile1D | ( | size_t | datasetId, | |
| size_t | xAxisId, | |||
| size_t | yAxisId, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D profile histogram based on the paper, dataset and x/y-axis IDs in the corresponding HepData record. The binnings will be obtained by reading the bundled AIDA data record file of the same filename as the analysis' name() property.
Definition at line 333 of file Analysis.cc.
References Analysis::bookProfile1D().
00335 { 00336 const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId); 00337 return bookProfile1D(axisCode, title, xtitle, ytitle); 00338 }
| IProfile1D * bookProfile1D | ( | const std::string & | name, | |
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D profile histogram based on the name in the corresponding AIDA file. The binnings will be obtained by reading the bundled AIDA data record file with the same filename as the analysis' name() property.
Definition at line 341 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00343 { 00344 // Get the bin edges (only read the AIDA file once) 00345 const BinEdges edges = binEdges(hname); 00346 _makeHistoDir(); 00347 const string path = histoPath(hname); 00348 IProfile1D* prof = histogramFactory().createProfile1D(path, title, edges); 00349 MSG_TRACE("Made profile histogram " << hname << " for " << name()); 00350 prof->setXTitle(xtitle); 00351 prof->setYTitle(ytitle); 00352 return prof; 00353 }
| IProfile1D * bookProfile1D | ( | const std::string & | name, | |
| const std::vector< double > & | binedges, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D profile histogram with non-uniform bins defined by the vector of bin edges binedges . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 370 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
00373 { 00374 _makeHistoDir(); 00375 const string path = histoPath(hname); 00376 IProfile1D* prof = histogramFactory().createProfile1D(path, title, binedges); 00377 MSG_TRACE("Made profile histogram " << hname << " for " << name()); 00378 prof->setXTitle(xtitle); 00379 prof->setYTitle(ytitle); 00380 return prof; 00381 }
| IProfile1D * bookProfile1D | ( | const std::string & | name, | |
| size_t | nbins, | |||
| double | lower, | |||
| double | upper, | |||
| const std::string & | title = "", |
|||
| const std::string & | xtitle = "", |
|||
| const std::string & | ytitle = "" | |||
| ) | [protected, inherited] |
Book a 1D profile histogram with nbins uniformly distributed across the range lower - upper . (NB. this returns a pointer rather than a reference since it will have to be stored in the analysis class - there's no point in forcing users to explicitly get the pointer from a reference before they can use it!)
Definition at line 356 of file Analysis.cc.
References Analysis::_makeHistoDir(), Analysis::histogramFactory(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().
Referenced by Analysis::bookProfile1D(), UA1_1990_S2044935::init(), STAR_2009_UE_HELEN::init(), STAR_2008_S7993412::init(), STAR_2006_S6860818::init(), MC_WPOL::init(), MC_PHOTONJETUE::init(), MC_LEADINGJETS::init(), MC_GENERIC::init(), H1_2000_S4129130::init(), H1_1994_S2919893::init(), DELPHI_2002_069_CONF_603::init(), DELPHI_1996_S3430090::init(), D0_1996_S3324664::init(), CMS_2011_S8884919::init(), CDF_2010_S8591881_QCD::init(), CDF_2010_S8591881_DY::init(), CDF_2009_S8233977::init(), CDF_2008_S7782535::init(), CDF_2008_NOTE_9351::init(), CDF_2008_LEADINGJETS::init(), CDF_2005_S6217184::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_2001_S4751469::init(), ATLAS_2011_S8994773::init(), ATLAS_2011_S8924791::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8894728::init(), ATLAS_2010_S8591806::init(), and ALICE_2010_S8706239::init().
00359 { 00360 _makeHistoDir(); 00361 const string path = histoPath(hname); 00362 IProfile1D* prof = histogramFactory().createProfile1D(path, title, nbins, lower, upper); 00363 MSG_TRACE("Made profile histogram " << hname << " for " << name()); 00364 prof->setXTitle(xtitle); 00365 prof->setYTitle(ytitle); 00366 return prof; 00367 }
| virtual std::string collider | ( | ) | const [inline, virtual, inherited] |
Collider on which the experiment ran.
Definition at line 163 of file Analysis.hh.
References AnalysisInfo::collider(), and Analysis::info().
00163 { 00164 return info().collider(); 00165 }
| double crossSection | ( | ) | const [protected, inherited] |
Get the process cross-section in pb. Throws if this hasn't been set.
Definition at line 171 of file Analysis.cc.
References Analysis::_crossSection, Analysis::_gotCrossSection, and Analysis::name().
Referenced by STAR_2006_S6870392::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_XS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_DIPHOTON::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7541902::finalize(), CDF_2008_S7540469::finalize(), CDF_2006_S6653332::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2001_S4517016::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), and ATLAS_2010_CONF_2010_049::finalize().
00171 { 00172 if (!_gotCrossSection || std::isnan(_crossSection)) { 00173 string errMsg = "You did not set the cross section for the analysis " + name(); 00174 throw Error(errMsg); 00175 } 00176 return _crossSection; 00177 }
| double crossSectionPerEvent | ( | ) | const [protected, inherited] |
Get the process cross-section per generated event in pb. Throws if this hasn't been set.
Definition at line 179 of file Analysis.cc.
References Analysis::_crossSection, and Analysis::sumOfWeights().
Referenced by UA1_1990_S2044935::finalize(), MC_WPOL::finalize(), MC_PHOTONJETS::finalize(), D0_2010_S8671338::finalize(), D0_2009_S8349509::finalize(), D0_2001_S4674421::finalize(), CDF_2007_S7057202::finalize(), CDF_1988_S1865951::finalize(), and ATLAS_2010_S8817804::finalize().
00179 { 00180 const double sumW = sumOfWeights(); 00181 assert(sumW > 0); 00182 return _crossSection / sumW; 00183 }
| IDataPointSetFactory & datapointsetFactory | ( | ) | [protected, inherited] |
Access the AIDA histogram factory of the controlling AnalysisHandler object.
Definition at line 65 of file Analysis.cc.
References AnalysisHandler::datapointsetFactory(), and Analysis::handler().
Referenced by Analysis::bookDataPointSet(), and Analysis::scale().
00065 { 00066 return handler().datapointsetFactory(); 00067 }
| virtual std::string description | ( | ) | const [inline, virtual, inherited] |
Get a full description of the analysis.
Full textual description of this analysis, what it is useful for, what experimental techniques are applied, etc. Should be treated as a chunk of restructuredText (http://docutils.sourceforge.net/rst.html), with equations to be rendered as LaTeX with amsmath operators.
Definition at line 144 of file Analysis.hh.
References AnalysisInfo::description(), and Analysis::info().
00144 { 00145 return info().description(); 00146 }
| virtual std::string experiment | ( | ) | const [inline, virtual, inherited] |
Experiment which performed and published this analysis.
Definition at line 158 of file Analysis.hh.
References AnalysisInfo::experiment(), and Analysis::info().
00158 { 00159 return info().experiment(); 00160 }
| void finalize | ( | ) | [inline, virtual] |
Finalize 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 235 of file MC_PHOTONJETUE.cc.
| Log & getLog | ( | ) | const [protected, inherited] |
Get a Log object based on the name() property of the calling analysis object.
Reimplemented from ProjectionApplier.
Definition at line 102 of file Analysis.cc.
References Analysis::name().
Referenced by CDF_2004_S5839831::_calcTransCones(), CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), D0_1996_S3214044::_safeMass(), CDF_1997_S3541940::_safeMass(), CDF_1996_S3349578::_safeMass(), CDF_1996_S3349578::_threeJetAnalysis(), ZEUS_2001_S4815815::analyze(), UA5_1986_S1583476::analyze(), TASSO_1990_S2148048::analyze(), STAR_2009_UE_HELEN::analyze(), STAR_2008_S7993412::analyze(), STAR_2006_S6870392::analyze(), STAR_2006_S6860818::analyze(), STAR_2006_S6500200::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), OPAL_1998_S3780481::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_LEADINGJETS::analyze(), MC_JetAnalysis::analyze(), JADE_1998_S3612880::analyze(), H1_1994_S2919893::analyze(), DELPHI_2003_WUD_03_11::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1996_S3430090::analyze(), DELPHI_1995_S3137023::analyze(), D0_2009_S8349509::analyze(), D0_2009_S8202443::analyze(), D0_2008_S7863608::analyze(), D0_2008_S7837160::analyze(), D0_2008_S7719523::analyze(), D0_2008_S7662670::analyze(), D0_2008_S7554427::analyze(), D0_2007_S7075677::analyze(), D0_2004_S5992206::analyze(), CMS_2011_S8884919::analyze(), CDF_2010_S8591881_QCD::analyze(), CDF_2010_S8591881_DY::analyze(), CDF_2009_S8383952::analyze(), CDF_2008_S8095620::analyze(), CDF_2008_S7782535::analyze(), CDF_2008_S7540469::analyze(), CDF_2008_NOTE_9351::analyze(), CDF_2008_LEADINGJETS::analyze(), CDF_2006_S6653332::analyze(), CDF_2004_S5839831::analyze(), CDF_2001_S4751469::analyze(), BELLE_2006_S6265367::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_1991_S2435284::analyze(), Analysis::binEdges(), CDF_2009_S8057893::CDF_2009_S8057893::analyze(), UA5_1986_S1583476::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), DELPHI_2003_WUD_03_11::finalize(), D0_2001_S4674421::finalize(), CDF_2006_S6653332::finalize(), ZEUS_2001_S4815815::init(), and ALEPH_2004_S5765862::init().
00102 { 00103 string logname = "Rivet.Analysis." + name(); 00104 return Log::getLog(logname); 00105 }
| const Projection& getProjection | ( | const std::string & | name | ) | const [inline, inherited] |
Get the named projection (non-templated, so returns as a reference to a Projection base class).
Definition at line 57 of file ProjectionApplier.hh.
References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler().
00057 { 00058 return getProjHandler().getProjection(*this, name); 00059 }
| const PROJ& getProjection | ( | const std::string & | name | ) | const [inline, inherited] |
Get the named projection, specifying return type via a template argument.
Definition at line 49 of file ProjectionApplier.hh.
References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler().
Referenced by ProjectionApplier::_applyProjection(), Rivet::pcmp(), and Hemispheres::project().
00049 { 00050 const Projection& p = getProjHandler().getProjection(*this, name); 00051 return pcast<PROJ>(p); 00052 }
| std::set<ConstProjectionPtr> getProjections | ( | ) | const [inline, inherited] |
Get the contained projections, including recursion.
Definition at line 42 of file ProjectionApplier.hh.
References ProjectionHandler::DEEP, ProjectionHandler::getChildProjections(), and ProjectionApplier::getProjHandler().
Referenced by Projection::beamPairs().
00042 { 00043 return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP); 00044 }
| ProjectionHandler& getProjHandler | ( | ) | const [inline, protected, inherited] |
Get a reference to the ProjectionHandler for this thread.
Definition at line 94 of file ProjectionApplier.hh.
References ProjectionApplier::_projhandler.
Referenced by ProjectionApplier::_addProjection(), ProjectionApplier::getProjection(), ProjectionApplier::getProjections(), and ProjectionApplier::~ProjectionApplier().
00094 { 00095 return _projhandler; 00096 }
| AnalysisHandler& handler | ( | ) | const [inline, inherited] |
Access the controlling AnalysisHandler object.
Definition at line 285 of file Analysis.hh.
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().
00285 { return *_analysishandler; }
| const string histoDir | ( | ) | const [protected, inherited] |
Get the canonical histogram "directory" path for this analysis.
- Todo:
- This doesn't change: calc and cache at first use!
Definition at line 83 of file Analysis.cc.
References Analysis::handler(), Analysis::name(), and AnalysisHandler::runName().
Referenced by Analysis::_makeHistoDir(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), MC_WJETS::finalize(), D0_2008_S7837160::finalize(), D0_2008_S7719523::finalize(), CMS_2011_S8978280::finalize(), and Analysis::histoPath().
00083 { 00084 /// @todo This doesn't change: calc and cache at first use! 00085 string path = "/" + name(); 00086 if (handler().runName().length() > 0) { 00087 path = "/" + handler().runName() + path; 00088 } 00089 while (find_first(path, "//")) { 00090 replace_all(path, "//", "/"); 00091 } 00092 return path; 00093 }
| IHistogramFactory & histogramFactory | ( | ) | [protected, inherited] |
Access the AIDA histogram factory of the controlling AnalysisHandler object.
Definition at line 60 of file Analysis.cc.
References Analysis::handler(), and AnalysisHandler::histogramFactory().
Referenced by Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), 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(), D0_2008_S7719523::finalize(), CMS_2011_S8978280::finalize(), and ATLAS_2011_S9002537::finalize().
00060 { 00061 return handler().histogramFactory(); 00062 }
| const string histoPath | ( | const std::string & | hname | ) | const [protected, inherited] |
Get the canonical histogram path for the named histogram in this analysis.
Definition at line 96 of file Analysis.cc.
References Analysis::histoDir().
Referenced by Analysis::bookDataPointSet(), Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), MC_JetAnalysis::finalize(), and MC_GENERIC::finalize().
00096 { 00097 const string path = histoDir() + "/" + hname; 00098 return path; 00099 }
| AnalysisInfo& info | ( | ) | [inline, inherited] |
Get the actual AnalysisInfo object in which all this metadata is stored (non-const).
Definition at line 237 of file Analysis.hh.
References Analysis::_info.
| const AnalysisInfo& info | ( | ) | const [inline, inherited] |
Get the actual AnalysisInfo object in which all this metadata is stored.
Definition at line 102 of file Analysis.hh.
References Analysis::_info.
Referenced by Analysis::authors(), Analysis::bibKey(), Analysis::bibTeX(), Analysis::collider(), Analysis::description(), Analysis::experiment(), Analysis::name(), Analysis::needsCrossSection(), Analysis::references(), Analysis::requiredBeams(), Analysis::requiredEnergies(), Analysis::runInfo(), Analysis::setNeedsCrossSection(), Analysis::setRequiredBeams(), Analysis::setRequiredEnergies(), Analysis::spiresId(), Analysis::status(), Analysis::summary(), Analysis::todos(), and Analysis::year().
| void init | ( | ) | [inline, virtual] |
Initialize this analysis object. A concrete class should here book all necessary histograms. An overridden function must make sure it first calls the base class function.
Implements Analysis.
Definition at line 27 of file MC_PHOTONJETUE.cc.
References MC_PHOTONJETUE::_hist_jetgamma_dphi, MC_PHOTONJETUE::_hist_jetgamma_dR, MC_PHOTONJETUE::_hist_pcptave_gamma, MC_PHOTONJETUE::_hist_pcptave_jet, MC_PHOTONJETUE::_hist_pcptsum_gamma, MC_PHOTONJETUE::_hist_pcptsum_jet, MC_PHOTONJETUE::_hist_pmaxcptsum_gamma, MC_PHOTONJETUE::_hist_pmaxcptsum_jet, MC_PHOTONJETUE::_hist_pmaxnchg_gamma, MC_PHOTONJETUE::_hist_pmaxnchg_jet, MC_PHOTONJETUE::_hist_pmincptsum_gamma, MC_PHOTONJETUE::_hist_pmincptsum_jet, MC_PHOTONJETUE::_hist_pminnchg_gamma, MC_PHOTONJETUE::_hist_pminnchg_jet, MC_PHOTONJETUE::_hist_pnchg_gamma, MC_PHOTONJETUE::_hist_pnchg_jet, IdentifiedFinalState::acceptId(), ProjectionApplier::addProjection(), FastJets::ANTIKT, Analysis::bookHistogram1D(), Analysis::bookProfile1D(), Rivet::GeV, Rivet::PHOTON, and Rivet::PI.
00027 { 00028 // Final state for the jet finding 00029 const FinalState fsj(-4.0, 4.0, 0.1*GeV); 00030 addProjection(fsj, "FSJ"); 00031 addProjection(FastJets(fsj, FastJets::ANTIKT, 0.7), "Jets"); 00032 00033 // Charged final state for the distributions 00034 const ChargedFinalState cfs(-2.0, 2.0, 0.2*GeV); 00035 addProjection(cfs, "Tracks"); 00036 00037 // Photons (for isolation) 00038 const FinalState fsp(-2.0, 2.0, 20.0*GeV); 00039 IdentifiedFinalState photonfs(fsp); 00040 photonfs.acceptId(PHOTON); 00041 addProjection(photonfs, "Photons"); 00042 00043 // Histograms 00044 _hist_jetgamma_dR = bookHistogram1D("gammajet-dR", 52, 0.0, 5.2); 00045 _hist_jetgamma_dphi = bookHistogram1D("gammajet-dphi", 50, 0.0, PI); 00046 // 00047 const double MAXPT1 = 50.0; 00048 _hist_pnchg_jet = bookProfile1D("trans-nchg-jet", 50, 0.0, MAXPT1); 00049 _hist_pmaxnchg_jet = bookProfile1D("trans-maxnchg-jet", 50, 0.0, MAXPT1); 00050 _hist_pminnchg_jet = bookProfile1D("trans-minnchg-jet", 50, 0.0, MAXPT1); 00051 _hist_pcptsum_jet = bookProfile1D("trans-ptsum-jet", 50, 0.0, MAXPT1); 00052 _hist_pmaxcptsum_jet = bookProfile1D("trans-maxptsum-jet", 50, 0.0, MAXPT1); 00053 _hist_pmincptsum_jet = bookProfile1D("trans-minptsum-jet", 50, 0.0, MAXPT1); 00054 _hist_pcptave_jet = bookProfile1D("trans-ptavg-jet", 50, 0.0, MAXPT1); 00055 // 00056 _hist_pnchg_gamma = bookProfile1D("trans-nchg-gamma", 50, 0.0, MAXPT1); 00057 _hist_pmaxnchg_gamma = bookProfile1D("trans-maxnchg-gamma", 50, 0.0, MAXPT1); 00058 _hist_pminnchg_gamma = bookProfile1D("trans-minnchg-gamma", 50, 0.0, MAXPT1); 00059 _hist_pcptsum_gamma = bookProfile1D("trans-ptsum-gamma", 50, 0.0, MAXPT1); 00060 _hist_pmaxcptsum_gamma = bookProfile1D("trans-maxptsum-gamma", 50, 0.0, MAXPT1); 00061 _hist_pmincptsum_gamma = bookProfile1D("trans-minptsum-gamma", 50, 0.0, MAXPT1); 00062 _hist_pcptave_gamma = bookProfile1D("trans-ptavg-gamma", 50, 0.0, MAXPT1); 00063 }
| bool isCompatible | ( | const PdgIdPair & | beams, | |
| const std::pair< double, double > & | energies | |||
| ) | const [inherited] |
Check if analysis is compatible with the provided beam particle IDs and energies.
Check if analysis is compatible with the provided beam particle IDs and energies.
Definition at line 127 of file Analysis.cc.
References Analysis::beams(), and Analysis::isCompatible().
00127 { 00128 PdgIdPair beams(beam1, beam2); 00129 pair<double,double> energies(e1, e2); 00130 return isCompatible(beams, energies); 00131 }
| bool isCompatible | ( | const ParticlePair & | beams | ) | const [inherited] |
Check if analysis is compatible with the provided beam particle IDs and energies.
Definition at line 121 of file Analysis.cc.
Referenced by Analysis::isCompatible().
00121 { 00122 return isCompatible(beams.first.pdgId(), beams.second.pdgId(), 00123 beams.first.energy(), beams.second.energy()); 00124 }
| BinEdges logBinEdges | ( | size_t | nbins, | |
| double | lower, | |||
| double | upper | |||
| ) | [protected, inherited] |
Get bin edges with logarithmic widths.
Definition at line 229 of file Analysis.cc.
Referenced by MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), and MC_DIPHOTON::init().
00229 { 00230 assert(lower>0.0); 00231 assert(upper>lower); 00232 double loglower=log10(lower); 00233 double logupper=log10(upper); 00234 vector<double> binedges; 00235 double stepwidth=(logupper-loglower)/double(nbins); 00236 for (size_t i=0; i<=nbins; ++i) { 00237 binedges.push_back(pow(10.0, loglower+double(i)*stepwidth)); 00238 } 00239 return binedges; 00240 }
| virtual std::string name | ( | ) | const [inline, virtual, inherited] |
Get the name of the analysis.
By default this is computed by combining the results of the experiment, year and Spires ID metadata methods and you should only override it if there's a good reason why those won't work.
Implements ProjectionApplier.
Definition at line 112 of file Analysis.hh.
References Analysis::_defaultname, Analysis::info(), and AnalysisInfo::name().
Referenced by Analysis::_cacheBinEdges(), Analysis::_cacheXAxisData(), Analysis::_makeHistoDir(), Analysis::binEdges(), Analysis::bookDataPointSet(), Analysis::bookHistogram1D(), Analysis::bookHistogram2D(), Analysis::bookProfile1D(), Analysis::crossSection(), Analysis::getLog(), Analysis::histoDir(), Analysis::normalize(), HistoHandler::registerAnalysisObject(), and Analysis::scale().
00112 { 00113 return (info().name().empty()) ? _defaultname : info().name(); 00114 }
| bool needsCrossSection | ( | ) | const [inline, inherited] |
Return true if this analysis needs to know the process cross-section.
Definition at line 221 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::needsCrossSection().
00221 { 00222 return info().needsCrossSection(); 00223 }
| void normalize | ( | AIDA::IHistogram2D *& | histo, | |
| double | norm = 1.0 | |||
| ) | [inherited] |
Normalize the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path. It has the same effect as scale(histo, norm/sumOfWeights).
- Parameters:
-
histo The histogram to be normalised. norm The new area of the histogram.
- Warning:
- The old histogram will be deleted, and its pointer set to zero.
Definition at line 527 of file Analysis.cc.
References MSG_ERROR, MSG_TRACE, MSG_WARNING, Analysis::name(), Analysis::scale(), and Analysis::tree().
00527 { 00528 if (!histo) { 00529 MSG_ERROR("Failed to normalize histo=NULL in analysis " 00530 << name() << " (norm=" << norm << ")"); 00531 return; 00532 } 00533 const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo)); 00534 MSG_TRACE("Normalizing histo " << hpath << " to " << norm); 00535 00536 double oldintg = 0.0; 00537 int nxBins = histo->xAxis().bins(); 00538 int nyBins = histo->yAxis().bins(); 00539 for (int ixBin = 0; ixBin != nxBins; ++ixBin) 00540 for (int iyBin = 0; iyBin != nyBins; ++iyBin) { 00541 // Leaving out factor of binWidth because AIDA's "height" 00542 // already includes a width factor. 00543 oldintg += histo->binHeight(ixBin, iyBin); // * histo->axis().binWidth(iBin); 00544 } 00545 if (oldintg == 0.0) { 00546 MSG_WARNING("Histo " << hpath << " has null integral during normalization"); 00547 return; 00548 } 00549 00550 // Scale by the normalisation factor. 00551 scale(histo, norm/oldintg); 00552 }
| void normalize | ( | AIDA::IHistogram1D *& | histo, | |
| double | norm = 1.0 | |||
| ) | [inherited] |
Normalize the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path. It has the same effect as scale(histo, norm/sumOfWeights).
- Parameters:
-
histo The histogram to be normalised. norm The new area of the histogram.
- Warning:
- The old histogram will be deleted, and its pointer set to zero.
Definition at line 460 of file Analysis.cc.
References MSG_ERROR, MSG_TRACE, MSG_WARNING, Analysis::name(), Analysis::scale(), and Analysis::tree().
Referenced by TASSO_1990_S2148048::finalize(), STAR_2008_S7869363::finalize(), OPAL_2001_S4553896::finalize(), OPAL_1998_S3780481::finalize(), JADE_1998_S3612880::finalize(), H1_1994_S2919893::finalize(), ExampleAnalysis::finalize(), DELPHI_2003_WUD_03_11::finalize(), DELPHI_2002_069_CONF_603::finalize(), DELPHI_1996_S3430090::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8821313::finalize(), D0_2010_S8671338::finalize(), D0_2009_S8320160::finalize(), D0_2008_S7554427::finalize(), D0_2008_S6879055::finalize(), D0_2007_S7075677::finalize(), D0_2004_S5992206::finalize(), D0_2001_S4674421::finalize(), D0_1996_S3324664::finalize(), D0_1996_S3214044::finalize(), CMS_2011_S8968497::finalize(), CMS_2011_S8957746::finalize(), CMS_2011_S8884919::finalize(), CDF_2009_NOTE_9936::finalize(), CDF_2004_S5839831::finalize(), CDF_2002_S4796047::finalize(), CDF_2001_S4751469::finalize(), CDF_1997_S3541940::finalize(), CDF_1996_S3418421::finalize(), CDF_1996_S3349578::finalize(), CDF_1996_S3108457::finalize(), CDF_1994_S2952106::finalize(), CDF_1993_S2742446::finalize(), ATLAS_2011_S8971293::finalize(), ALICE_2010_S8625980::finalize(), ALICE_2010_S8624100::finalize(), ALEPH_2004_S5765862::finalize(), and ALEPH_1996_S3486095::finalize().
00460 { 00461 if (!histo) { 00462 MSG_ERROR("Failed to normalize histo=NULL in analysis " 00463 << name() << " (norm=" << norm << ")"); 00464 return; 00465 } 00466 const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo)); 00467 MSG_TRACE("Normalizing histo " << hpath << " to " << norm); 00468 00469 double oldintg = 0.0; 00470 int nBins = histo->axis().bins(); 00471 for (int iBin = 0; iBin != nBins; ++iBin) { 00472 // Leaving out factor of binWidth because AIDA's "height" already includes a width factor. 00473 oldintg += histo->binHeight(iBin); // * histo->axis().binWidth(iBin); 00474 } 00475 if (oldintg == 0.0) { 00476 MSG_WARNING("Histo " << hpath << " has null integral during normalization"); 00477 return; 00478 } 00479 00480 // Scale by the normalisation factor. 00481 scale(histo, norm/oldintg); 00482 }
| size_t numEvents | ( | ) | const [protected, inherited] |
Get the number of events seen (via the analysis handler). Use in the finalize phase only.
Definition at line 108 of file Analysis.cc.
References Analysis::handler(), and AnalysisHandler::numEvents().
00108 { 00109 return handler().numEvents(); 00110 }
| virtual std::vector<std::string> references | ( | ) | const [inline, virtual, inherited] |
Journal, and preprint references.
Definition at line 173 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::references().
00173 { 00174 return info().references(); 00175 }
| virtual const std::vector<PdgIdPair>& requiredBeams | ( | ) | const [inline, virtual, inherited] |
Return the allowed pairs of incoming beams required by this analysis.
Definition at line 199 of file Analysis.hh.
References AnalysisInfo::beams(), and Analysis::info().
00199 { 00200 return info().beams(); 00201 }
| virtual const std::vector<std::pair<double, double> >& requiredEnergies | ( | ) | const [inline, virtual, inherited] |
Sets of valid beam energy pairs, in GeV.
Definition at line 210 of file Analysis.hh.
References AnalysisInfo::energies(), and Analysis::info().
Referenced by Analysis::setRequiredEnergies().
00210 { 00211 return info().energies(); 00212 }
| virtual std::string runInfo | ( | ) | const [inline, virtual, inherited] |
Information about the events needed as input for this analysis.
Event types, energies, kinematic cuts, particles to be considered stable, etc. etc. Should be treated as a restructuredText bullet list (http://docutils.sourceforge.net/rst.html)
Definition at line 153 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::runInfo().
00153 { 00154 return info().runInfo(); 00155 }
| void scale | ( | AIDA::IHistogram2D *& | histo, | |
| double | scale | |||
| ) | [inherited] |
Multiplicatively scale the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path.
- Parameters:
-
histo The histogram to be scaled. scale The factor used to multiply the histogram bin heights.
- Warning:
- The old histogram will be deleted, and its pointer set to zero.
Definition at line 555 of file Analysis.cc.
References Analysis::datapointsetFactory(), MSG_ERROR, MSG_TRACE, Analysis::name(), and Analysis::tree().
00555 { 00556 if (!histo) { 00557 MSG_ERROR("Failed to scale histo=NULL in analysis " 00558 << name() << " (scale=" << scale << ")"); 00559 return; 00560 } 00561 const string hpath = 00562 tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo)); 00563 MSG_TRACE("Scaling histo " << hpath); 00564 00565 vector<double> x, y, z, ex, ey, ez; 00566 for (size_t ix = 0, Nx = histo->xAxis().bins(); ix < Nx; ++ix) 00567 for (size_t iy = 0, Ny = histo->yAxis().bins(); iy < Ny; ++iy) { 00568 x.push_back(0.5 * (histo->xAxis().binLowerEdge(ix) + 00569 histo->xAxis().binUpperEdge(ix))); 00570 ex.push_back(histo->xAxis().binWidth(ix)*0.5); 00571 y.push_back(0.5 * (histo->yAxis().binLowerEdge(iy) + 00572 histo->yAxis().binUpperEdge(iy))); 00573 ey.push_back(histo->yAxis().binWidth(iy)*0.5); 00574 00575 // "Bin height" is a misnomer in the AIDA spec: width is neglected. 00576 // We'd like to do this: y.push_back(histo->binHeight(i) * scale); 00577 z.push_back(histo->binHeight(ix, iy)*scale/ 00578 (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); 00579 // "Bin error" is a misnomer in the AIDA spec: width is neglected. 00580 // We'd like to do this: ey.push_back(histo->binError(i) * scale); 00581 ez.push_back(histo->binError(ix, iy)*scale/ 00582 (histo->xAxis().binWidth(ix)*histo->yAxis().binWidth(iy))); 00583 } 00584 00585 string title = histo->title(); 00586 string xtitle = histo->xtitle(); 00587 string ytitle = histo->ytitle(); 00588 string ztitle = histo->ztitle(); 00589 00590 tree().mkdir("/tmpnormalize"); 00591 tree().mv(hpath, "/tmpnormalize"); 00592 00593 AIDA::IDataPointSet* dps = 00594 datapointsetFactory().createXYZ(hpath, title, x, y, z, ex, ey, ez); 00595 dps->setXTitle(xtitle); 00596 dps->setYTitle(ytitle); 00597 dps->setZTitle(ztitle); 00598 00599 tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo))); 00600 tree().rmdir("/tmpnormalize"); 00601 00602 // Set histo pointer to null - it can no longer be used. 00603 histo = 0; 00604 }
| void scale | ( | AIDA::IHistogram1D *& | histo, | |
| double | scale | |||
| ) | [inherited] |
Multiplicatively scale the given histogram, histo. After this call the histogram will have been transformed to a DataPointSet with the same name and path.
- Parameters:
-
histo The histogram to be scaled. scale The factor used to multiply the histogram bin heights.
- Warning:
- The old histogram will be deleted, and its pointer set to zero.
Definition at line 485 of file Analysis.cc.
References Analysis::datapointsetFactory(), MSG_ERROR, MSG_TRACE, Analysis::name(), and Analysis::tree().
Referenced by UA5_1989_S1926373::finalize(), UA5_1987_S1640666::finalize(), UA5_1986_S1583476::finalize(), UA5_1982_S875503::finalize(), UA1_1990_S2044935::finalize(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), SFM_1984_S1178091::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_2004_S6132243::finalize(), OPAL_1998_S3780481::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_WWJETS::finalize(), MC_WPOL::finalize(), MC_WJETS::finalize(), MC_PHOTONJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), LHCB_2010_S8758301::finalize(), JADE_OPAL_2000_S4300807::finalize(), H1_2000_S4129130::finalize(), H1_1994_S2919893::finalize(), E735_1998_S3905616::finalize(), DELPHI_1996_S3430090::finalize(), D0_2010_S8671338::finalize(), D0_2010_S8570965::finalize(), D0_2009_S8349509::finalize(), D0_2009_S8202443::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S8978280::finalize(), CMS_2010_S8656010::finalize(), CMS_2010_S8547297::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7540469::finalize(), CDF_2007_S7057202::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), CDF_1990_S2089246::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2010_S8591806::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ALICE_2010_S8706239::finalize(), ALICE_2010_S8625980::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), ALEPH_1991_S2435284::finalize(), Analysis::normalize(), and BinnedHistogram< T >::scale().
00485 { 00486 if (!histo) { 00487 MSG_ERROR("Failed to scale histo=NULL in analysis " 00488 << name() << " (scale=" << scale << ")"); 00489 return; 00490 } 00491 const string hpath = tree().findPath(dynamic_cast<const AIDA::IManagedObject&>(*histo)); 00492 MSG_TRACE("Scaling histo " << hpath); 00493 00494 vector<double> x, y, ex, ey; 00495 for (size_t i = 0, N = histo->axis().bins(); i < N; ++i) { 00496 x.push_back(0.5 * (histo->axis().binLowerEdge(i) + histo->axis().binUpperEdge(i))); 00497 ex.push_back(histo->axis().binWidth(i)*0.5); 00498 00499 // "Bin height" is a misnomer in the AIDA spec: width is neglected. 00500 // We'd like to do this: y.push_back(histo->binHeight(i) * scale); 00501 y.push_back(histo->binHeight(i)*scale/histo->axis().binWidth(i)); 00502 00503 // "Bin error" is a misnomer in the AIDA spec: width is neglected. 00504 // We'd like to do this: ey.push_back(histo->binError(i) * scale); 00505 ey.push_back(histo->binError(i)*scale/histo->axis().binWidth(i)); 00506 } 00507 00508 string title = histo->title(); 00509 string xtitle = histo->xtitle(); 00510 string ytitle = histo->ytitle(); 00511 00512 tree().mkdir("/tmpnormalize"); 00513 tree().mv(hpath, "/tmpnormalize"); 00514 00515 AIDA::IDataPointSet* dps = datapointsetFactory().createXY(hpath, title, x, y, ex, ey); 00516 dps->setXTitle(xtitle); 00517 dps->setYTitle(ytitle); 00518 00519 tree().rm(tree().findPath(dynamic_cast<AIDA::IManagedObject&>(*histo))); 00520 tree().rmdir("/tmpnormalize"); 00521 00522 // Set histo pointer to null - it can no longer be used. 00523 histo = 0; 00524 }
Set the required beams
- Deprecated:
- To be removed in 2.0.0. Use .info file and AnalysisInfo class instead
- Todo:
- Print out a warning to use setRequiredBeams() instead (and really to use .info files)
Definition at line 244 of file Analysis.hh.
References Analysis::setRequiredBeams().
Referenced by ALEPH_1991_S2435284::ALEPH_1991_S2435284(), ALEPH_1996_S3196992::ALEPH_1996_S3196992(), ALEPH_1996_S3486095::ALEPH_1996_S3486095(), ALEPH_2004_S5765862::ALEPH_2004_S5765862(), ALICE_2010_S8624100::ALICE_2010_S8624100(), ALICE_2010_S8625980::ALICE_2010_S8625980(), ALICE_2010_S8706239::ALICE_2010_S8706239(), ATLAS_2010_S8817804::ATLAS_2010_S8817804(), ATLAS_2011_S8924791::ATLAS_2011_S8924791(), ATLAS_2011_S9120807::ATLAS_2011_S9120807(), BELLE_2006_S6265367::BELLE_2006_S6265367(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1990_S2089246::CDF_1990_S2089246(), CDF_1993_S2742446::CDF_1993_S2742446(), CDF_1994_S2952106::CDF_1994_S2952106(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1996_S3349578::CDF_1996_S3349578(), CDF_1996_S3418421::CDF_1996_S3418421(), CDF_1997_S3541940::CDF_1997_S3541940(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2001_S4751469::CDF_2001_S4751469(), CDF_2002_S4796047::CDF_2002_S4796047(), CDF_2004_S5839831::CDF_2004_S5839831(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2005_S6217184::CDF_2005_S6217184(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_LEADINGJETS::CDF_2008_LEADINGJETS(), CDF_2008_NOTE_9351::CDF_2008_NOTE_9351(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7782535::CDF_2008_S7782535(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2008_S8095620::CDF_2008_S8095620(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8057893::CDF_2009_S8057893::CDF_2009_S8057893(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), CDF_2010_S8591881_DY::CDF_2010_S8591881_DY(), CDF_2010_S8591881_QCD::CDF_2010_S8591881_QCD(), CMS_2011_S8957746::CMS_2011_S8957746(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_2001_S4674421::D0_2001_S4674421(), D0_2004_S5992206::D0_2004_S5992206(), D0_2006_S6438750::D0_2006_S6438750(), D0_2007_S7075677::D0_2007_S7075677(), D0_2008_S6879055::D0_2008_S6879055(), D0_2008_S7554427::D0_2008_S7554427(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7837160::D0_2008_S7837160(), D0_2008_S7863608::D0_2008_S7863608(), D0_2009_S8202443::D0_2009_S8202443(), D0_2009_S8320160::D0_2009_S8320160(), D0_2009_S8349509::D0_2009_S8349509(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), D0_2010_S8671338::D0_2010_S8671338(), DELPHI_1995_S3137023::DELPHI_1995_S3137023(), DELPHI_1996_S3430090::DELPHI_1996_S3430090(), DELPHI_2002_069_CONF_603::DELPHI_2002_069_CONF_603(), E735_1998_S3905616::E735_1998_S3905616(), H1_1994_S2919893::H1_1994_S2919893(), H1_1995_S3167097::H1_1995_S3167097(), H1_2000_S4129130::H1_2000_S4129130(), JADE_1998_S3612880::JADE_1998_S3612880(), JADE_OPAL_2000_S4300807::JADE_OPAL_2000_S4300807(), LHCB_2010_S8758301::LHCB_2010_S8758301(), OPAL_1993_S2692198::OPAL_1993_S2692198(), OPAL_1998_S3780481::OPAL_1998_S3780481(), OPAL_2001_S4553896::OPAL_2001_S4553896(), PDG_HADRON_MULTIPLICITIES::PDG_HADRON_MULTIPLICITIES(), PDG_HADRON_MULTIPLICITIES_RATIOS::PDG_HADRON_MULTIPLICITIES_RATIOS(), SFM_1984_S1178091::SFM_1984_S1178091(), STAR_2006_S6500200::STAR_2006_S6500200(), STAR_2006_S6860818::STAR_2006_S6860818(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), STAR_2008_S7993412::STAR_2008_S7993412(), STAR_2009_UE_HELEN::STAR_2009_UE_HELEN(), TASSO_1990_S2148048::TASSO_1990_S2148048(), UA1_1990_S2044935::UA1_1990_S2044935(), UA5_1986_S1583476::UA5_1986_S1583476(), UA5_1987_S1640666::UA5_1987_S1640666(), UA5_1988_S1867512::UA5_1988_S1867512(), UA5_1989_S1926373::UA5_1989_S1926373(), and ZEUS_2001_S4815815::ZEUS_2001_S4815815().
00244 { 00245 /// @todo Print out a warning to use setRequiredBeams() instead (and really to use .info files) 00246 return setRequiredBeams(std::vector<PdgIdPair>(1, make_pair(beam1, beam2))); 00247 }
| Analysis & setCrossSection | ( | double | xs | ) | [inherited] |
Set the cross section from the generator.
Definition at line 165 of file Analysis.cc.
References Analysis::_crossSection, and Analysis::_gotCrossSection.
00165 { 00166 _crossSection = xs; 00167 _gotCrossSection = true; 00168 return *this; 00169 }
| Analysis& setNeedsCrossSection | ( | bool | needed = true |
) | [inline, inherited] |
Declare whether this analysis needs to know the process cross-section from the generator.
Definition at line 225 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::setNeedsCrossSection().
Referenced by ALICE_2010_S8624100::ALICE_2010_S8624100(), ALICE_2010_S8625980::ALICE_2010_S8625980(), ALICE_2010_S8706239::ALICE_2010_S8706239(), ATLAS_2010_CONF_2010_049::ATLAS_2010_CONF_2010_049(), ATLAS_2010_S8591806::ATLAS_2010_S8591806(), ATLAS_2010_S8817804::ATLAS_2010_S8817804(), ATLAS_2010_S8914702::ATLAS_2010_S8914702(), ATLAS_2010_S8919674::ATLAS_2010_S8919674(), ATLAS_2011_CONF_2011_090::ATLAS_2011_CONF_2011_090(), ATLAS_2011_S8971293::ATLAS_2011_S8971293(), ATLAS_2011_S8983313::ATLAS_2011_S8983313(), ATLAS_2011_S9019561::ATLAS_2011_S9019561(), ATLAS_2011_S9120807::ATLAS_2011_S9120807(), CDF_1988_S1865951::CDF_1988_S1865951(), CDF_1993_S2742446::CDF_1993_S2742446(), CDF_1996_S3108457::CDF_1996_S3108457(), CDF_1998_S3618439::CDF_1998_S3618439(), CDF_2000_S4155203::CDF_2000_S4155203(), CDF_2000_S4266730::CDF_2000_S4266730(), CDF_2001_S4517016::CDF_2001_S4517016(), CDF_2001_S4563131::CDF_2001_S4563131(), CDF_2005_S6080774::CDF_2005_S6080774(), CDF_2006_S6450792::CDF_2006_S6450792(), CDF_2006_S6653332::CDF_2006_S6653332(), CDF_2007_S7057202::CDF_2007_S7057202(), CDF_2008_S7540469::CDF_2008_S7540469(), CDF_2008_S7541902::CDF_2008_S7541902(), CDF_2008_S7828950::CDF_2008_S7828950(), CDF_2008_S8093652::CDF_2008_S8093652(), CDF_2009_NOTE_9936::CDF_2009_NOTE_9936(), CDF_2009_S8233977::CDF_2009_S8233977(), CDF_2009_S8383952::CDF_2009_S8383952(), CDF_2009_S8436959::CDF_2009_S8436959(), CMS_2011_S8884919::CMS_2011_S8884919(), D0_1996_S3214044::D0_1996_S3214044(), D0_1996_S3324664::D0_1996_S3324664(), D0_2000_S4480767::D0_2000_S4480767(), D0_2001_S4674421::D0_2001_S4674421(), D0_2006_S6438750::D0_2006_S6438750(), D0_2008_S7662670::D0_2008_S7662670(), D0_2008_S7719523::D0_2008_S7719523(), D0_2008_S7863608::D0_2008_S7863608(), D0_2009_S8349509::D0_2009_S8349509(), D0_2010_S8566488::D0_2010_S8566488(), D0_2010_S8570965::D0_2010_S8570965(), D0_2010_S8671338::D0_2010_S8671338(), D0_2010_S8821313::D0_2010_S8821313(), LHCB_2010_S8758301::LHCB_2010_S8758301(), MC_DIPHOTON::MC_DIPHOTON(), MC_HJETS::MC_HJETS(), MC_JetAnalysis::MC_JetAnalysis(), MC_JETS::MC_JETS(), MC_PHOTONJETS::MC_PHOTONJETS(), MC_WJETS::MC_WJETS(), MC_WPOL::MC_WPOL(), MC_WWJETS::MC_WWJETS(), MC_XS::MC_XS(), MC_ZJETS::MC_ZJETS(), MC_ZZJETS::MC_ZZJETS(), OPAL_2001_S4553896::OPAL_2001_S4553896(), STAR_2006_S6870392::STAR_2006_S6870392(), STAR_2008_S7869363::STAR_2008_S7869363(), and UA1_1990_S2044935::UA1_1990_S2044935().
00225 { 00226 info().setNeedsCrossSection(needed); 00227 return *this; 00228 }
| virtual Analysis& setRequiredBeams | ( | const std::vector< PdgIdPair > & | requiredBeams | ) | [inline, virtual, inherited] |
Declare the allowed pairs of incoming beams required by this analysis.
Definition at line 203 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::setBeams().
Referenced by Analysis::setBeams().
00203 { 00204 info().setBeams(requiredBeams); 00205 return *this; 00206 }
| virtual Analysis& setRequiredEnergies | ( | const std::vector< std::pair< double, double > > & | requiredEnergies | ) | [inline, virtual, inherited] |
Declare the list of valid beam energy pairs, in GeV.
Definition at line 214 of file Analysis.hh.
References Analysis::info(), Analysis::requiredEnergies(), and AnalysisInfo::setEnergies().
00214 { 00215 info().setEnergies(requiredEnergies); 00216 return *this; 00217 }
| virtual std::string spiresId | ( | ) | const [inline, virtual, inherited] |
Get a the SPIRES/Inspire ID code for this analysis.
Definition at line 117 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::spiresId().
00117 { 00118 return info().spiresId(); 00119 }
| double sqrtS | ( | ) | const [inherited] |
Centre of mass energy for this run.
Definition at line 70 of file Analysis.cc.
References Analysis::handler(), and AnalysisHandler::sqrtS().
Referenced by UA1_1990_S2044935::analyze(), TASSO_1990_S2148048::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), JADE_1998_S3612880::analyze(), CDF_2004_S5839831::analyze(), ATLAS_2010_S8918562::analyze(), ATLAS_2010_S8894728::analyze(), ALICE_2010_S8625980::analyze(), ALEPH_2004_S5765862::analyze(), UA1_1990_S2044935::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), JADE_1998_S3612880::finalize(), CMS_2011_S8978280::finalize(), CDF_2004_S5839831::finalize(), CDF_2002_S4796047::finalize(), ATLAS_2010_S8918562::finalize(), ALICE_2010_S8625980::finalize(), ALEPH_2004_S5765862::finalize(), UA5_1989_S1926373::init(), UA5_1988_S1867512::init(), UA5_1986_S1583476::init(), UA1_1990_S2044935::init(), TASSO_1990_S2148048::init(), SFM_1984_S1178091::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), PDG_HADRON_MULTIPLICITIES::init(), OPAL_2004_S6132243::init(), MC_ZZJETS::init(), MC_ZJETS::init(), MC_WWJETS::init(), MC_WPOL::init(), MC_WJETS::init(), MC_PHOTONJETS::init(), MC_JetAnalysis::init(), MC_HJETS::init(), MC_DIPHOTON::init(), JADE_OPAL_2000_S4300807::init(), JADE_1998_S3612880::init(), CMS_2011_S8978280::init(), CMS_2011_S8884919::init(), CMS_2010_S8547297::init(), CDF_2004_S5839831::init(), CDF_2002_S4796047::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), ATLAS_2011_S8994773::init(), ATLAS_2010_S8918562::init(), ATLAS_2010_S8894728::init(), ALICE_2010_S8625980::init(), ALICE_2010_S8624100::init(), and ALEPH_2004_S5765862::init().
00070 { 00071 return handler().sqrtS(); 00072 }
| virtual std::string status | ( | ) | const [inline, virtual, inherited] |
Whether this analysis is trusted (in any way!).
Definition at line 188 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::status().
| virtual std::string summary | ( | ) | const [inline, virtual, inherited] |
Get a short description of the analysis.
Short (one sentence) description used as an index entry. Use description() to provide full descriptive paragraphs of analysis details.
Definition at line 134 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::summary().
00134 { 00135 return info().summary(); 00136 }
| double sumOfWeights | ( | ) | const [protected, inherited] |
Get the sum of event weights seen (via the analysis handler). Use in the finalize phase only.
Definition at line 113 of file Analysis.cc.
References Analysis::handler(), and AnalysisHandler::sumOfWeights().
Referenced by Analysis::crossSectionPerEvent(), STAR_2006_S6870392::finalize(), STAR_2006_S6860818::finalize(), STAR_2006_S6500200::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), OPAL_1998_S3780481::finalize(), OPAL_1993_S2692198::finalize(), MC_ZZJETS::finalize(), MC_ZJETS::finalize(), MC_WWJETS::finalize(), MC_WJETS::finalize(), MC_JetAnalysis::finalize(), MC_HJETS::finalize(), MC_GENERIC::finalize(), MC_DIPHOTON::finalize(), JADE_OPAL_2000_S4300807::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8570965::finalize(), D0_2010_S8566488::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7719523::finalize(), D0_2008_S7662670::finalize(), D0_2006_S6438750::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S8978280::finalize(), CMS_2010_S8656010::finalize(), CMS_2010_S8547297::finalize(), CDF_2009_S8436959::finalize(), CDF_2009_S8383952::finalize(), CDF_2009_S8233977::finalize(), CDF_2008_S8093652::finalize(), CDF_2008_S7828950::finalize(), CDF_2008_S7541902::finalize(), CDF_2008_S7540469::finalize(), CDF_2006_S6653332::finalize(), CDF_2006_S6450792::finalize(), CDF_2005_S6080774::finalize(), CDF_2001_S4563131::finalize(), CDF_2001_S4517016::finalize(), CDF_2000_S4266730::finalize(), CDF_2000_S4155203::finalize(), CDF_1998_S3618439::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ALEPH_2004_S5765862::finalize(), ALEPH_1996_S3486095::finalize(), ALEPH_1996_S3196992::finalize(), and ALEPH_1991_S2435284::finalize().
00113 { 00114 return handler().sumOfWeights(); 00115 }
| virtual std::vector<std::string> todos | ( | ) | const [inline, virtual, inherited] |
Any work to be done on this analysis.
Definition at line 193 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::todos().
00193 { 00194 return info().todos(); 00195 }
| ITree & tree | ( | ) | [protected, inherited] |
Access the AIDA tree of the controlling AnalysisHandler object.
Definition at line 55 of file Analysis.cc.
References Analysis::handler(), and AnalysisHandler::tree().
Referenced by Analysis::_makeHistoDir(), Analysis::normalize(), and Analysis::scale().
00055 { 00056 return handler().tree(); 00057 }
| virtual std::string year | ( | ) | const [inline, virtual, inherited] |
When the original experimental analysis was published.
Definition at line 168 of file Analysis.hh.
References Analysis::info(), and AnalysisInfo::year().
00168 { 00169 return info().year(); 00170 }
Member Data Documentation
bool _allowProjReg [protected, inherited] |
Flag to forbid projection registration in analyses until the init phase.
Definition at line 140 of file ProjectionApplier.hh.
Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().
string _defaultname [protected, inherited] |
Name passed to constructor (used to find .info analysis data file, and as a fallback).
Definition at line 533 of file Analysis.hh.
Referenced by Analysis::Analysis(), and Analysis::name().
AIDA::IHistogram1D* _hist_jetgamma_dphi [private] |
Definition at line 243 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IHistogram1D* _hist_jetgamma_dR [private] |
Definition at line 242 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pcptave_gamma [private] |
Definition at line 251 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pcptave_jet [private] |
Definition at line 251 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pcptsum_gamma [private] |
Definition at line 248 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pcptsum_jet [private] |
Definition at line 248 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pmaxcptsum_gamma [private] |
Definition at line 249 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pmaxcptsum_jet [private] |
Definition at line 249 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pmaxnchg_gamma [private] |
Definition at line 246 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pmaxnchg_jet [private] |
Definition at line 246 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pmincptsum_gamma [private] |
Definition at line 250 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pmincptsum_jet [private] |
Definition at line 250 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pminnchg_gamma [private] |
Definition at line 247 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pminnchg_jet [private] |
Definition at line 247 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D * _hist_pnchg_gamma [private] |
Definition at line 245 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
AIDA::IProfile1D* _hist_pnchg_jet [private] |
Definition at line 245 of file MC_PHOTONJETUE.cc.
Referenced by MC_PHOTONJETUE::analyze(), and MC_PHOTONJETUE::init().
shared_ptr<AnalysisInfo> _info [protected, inherited] |
Pointer to analysis metadata object.
Definition at line 536 of file Analysis.hh.
Referenced by Analysis::Analysis(), and Analysis::info().
The documentation for this class was generated from the following file: