rivet is hosted by Hepforge, IPPP Durham
Inheritance diagram for ATLAS_2012_I1204447:
Collaboration diagram for ATLAS_2012_I1204447:

List of all members.

Public Member Functions

 ATLAS_2012_I1204447 ()
 Constructor.
void init ()
 Book histograms and initialise projections before the run.
void analyze (const Event &event)
 Perform the per-event analysis.
void finalize ()
 Normalise histograms etc., after the run.
AnalysissetCrossSection (double xs)
 Set the cross section from the generator.
AnalysisHandlerhandler () const
 Access the controlling AnalysisHandler object.
const vector< AnalysisObjectPtr > & analysisObjects () const
 List of registered analysis data objects.
vector< string > getSignalRegions ()
 Helper functions.
vector< int > getCutsPerSignalRegion (const string &signal_region, int onZ=0)
 Function giving all cut vales per kinematic variable (taking onZ for MET into account)
void fillEventCountsPerSR (const string &basic_signal_region, int onZ, double HTlep, double eTmiss, double HTjets, double meff, double weight)
FourMomentum get_tau_neutrino_mom (const Particle &p)
void get_prong_number (const GenParticle *p, unsigned int &nprong, bool &lep_decaying_tau)
double apply_reco_eff (int flavor, const Particle &p)
 Function giving fiducial lepton efficiency.
double getUpperLimit (const string &signal_region, bool observed)
 Function giving observed upper limit (visible cross-section)
int isonZ (const Particles &particles)
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 AnalysisInfoinfo () 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 inspireId () const
 Get the Inspire ID code for this analysis.
virtual std::string spiresId () const
 Get the SPIRES ID code for this analysis (~deprecated).
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 AnalysissetRequiredBeams (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 AnalysissetRequiredEnergies (const std::vector< std::pair< double, double > > &requiredEnergies)
 Declare the list of valid beam energy pairs, in GeV.
bool needsCrossSection () const
AnalysissetNeedsCrossSection (bool needed=true)
Internal metadata modifying methods
AnalysisInfoinfo ()
 Get the actual AnalysisInfo object in which all this metadata is stored (non-const).
Run conditions
const ParticlePairbeams () 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.
Histogram manipulation
Todo:
What follows should really be protected: only public to keep BinnedHistogram happy for now...
void normalize (Histo1DPtr histo, double norm=1.0, bool includeoverflows=true)
void normalize (Histo2DPtr histo, double norm=1.0, bool includeoverflows=true)
void scale (Histo1DPtr histo, double scale)
void scale (Histo2DPtr histo, double scale)
void divide (Histo1DPtr h1, Histo1DPtr h2, Scatter2DPtr s) const
void divide (const YODA::Histo1D &h1, const YODA::Histo1D &h2, Scatter2DPtr s) const
void divide (Profile1DPtr p1, Profile1DPtr p2, Scatter2DPtr s) const
void divide (const YODA::Profile1D &p1, const YODA::Profile1D &p2, Scatter2DPtr s) const
void divide (Histo2DPtr h1, Histo2DPtr h2, Scatter3DPtr s) const
void divide (const YODA::Histo2D &h1, const YODA::Histo2D &h2, Scatter3DPtr s) const
void divide (Profile2DPtr p1, Profile2DPtr p2, Scatter3DPtr s) const
void divide (const YODA::Profile2D &p1, const YODA::Profile2D &p2, Scatter3DPtr s) const
void efficiency (Histo1DPtr h1, Histo1DPtr h2, Scatter2DPtr s) const
void efficiency (const YODA::Histo1D &h1, const YODA::Histo1D &h2, Scatter2DPtr s) const
void asymm (Histo1DPtr h1, Histo1DPtr h2, Scatter2DPtr s) const
void asymm (const YODA::Histo1D &h1, const YODA::Histo1D &h2, Scatter2DPtr s) const
void integrate (Histo1DPtr h, Scatter2DPtr s) const
void integrate (const Histo1D &h, Scatter2DPtr s) const
Projection "getting" functions
std::set< ConstProjectionPtrgetProjections () 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 ProjectiongetProjection (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

LoggetLog () 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
ProjectionHandlergetProjHandler () const
 Get a reference to the ProjectionHandler for this thread.
Histogram paths
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.
const std::string histoPath (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get the canonical histogram path for the numbered histogram in this analysis.
const std::string makeAxisCode (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get the internal histogram name for given d, x and y (cf. HepData)
Histogram reference data
const YODA::Scatter2D & refData (const string &hname) const
 Get reference data for a named histo.
const YODA::Scatter2D & refData (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get reference data for a numbered histo.
1D histogram booking
Todo:
Provide 3D versions as well? (How to distinguish the signatures? Template magic or explicit name?)
Histo1DPtr bookHisto1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D histogram with nbins uniformly distributed across the range lower - upper .
Histo1DPtr bookHisto1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D histogram with non-uniform bins defined by the vector of bin edges binedges .
Histo1DPtr bookHisto1D (const std::string &name, const Scatter2D &refscatter, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D histogram with binning from a reference scatter.
Histo1DPtr bookHisto1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D histogram, using the binnings in the reference data histogram.
Histo1DPtr bookHisto1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
2D histogram booking
Histo2DPtr bookHisto2D (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="")
Histo2DPtr bookHisto2D (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="")
1D profile histogram booking
Profile1DPtr bookProfile1D (const std::string &name, size_t nbins, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D profile histogram with nbins uniformly distributed across the range lower - upper .
Profile1DPtr bookProfile1D (const std::string &name, const std::vector< double > &binedges, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D profile histogram with non-uniform bins defined by the vector of bin edges binedges .
Profile1DPtr bookProfile1D (const std::string &name, const Scatter2D &refscatter, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D profile histogram with binning from a reference scatter.
Profile1DPtr bookProfile1D (const std::string &name, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 1D profile histogram, using the binnings in the reference data histogram.
Profile1DPtr bookProfile1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
2D profile histogram booking
Profile2DPtr bookProfile2D (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="")
Profile2DPtr bookProfile2D (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="")
2D scatter booking

Book a 2D profile histogram with binning from a reference scatter.

Scatter2DPtr bookScatter2D (const std::string &name, bool copy_pts=false, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 2-dimensional data point set with the given name.
Scatter2DPtr bookScatter2D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId, bool copy_pts=false, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 2-dimensional data point set, using the binnings in the reference data histogram.
Scatter2DPtr bookScatter2D (const std::string &name, size_t npts, double lower, double upper, const std::string &title="", const std::string &xtitle="", const std::string &ytitle="")
 Book a 2-dimensional data point set with equally spaced x-points in a range.
Scatter2DPtr bookScatter2D (const std::string &hname, const std::vector< double > &binedges, const std::string &title, const std::string &xtitle, const std::string &ytitle)
 Book a 2-dimensional data point set based on provided contiguous "bin edges".
Data object registration, retrieval, and removal
void addAnalysisObject (AnalysisObjectPtr ao)
 Register a data object in the histogram system.
template<typename AO >
const shared_ptr< AO > getAnalysisObject (const std::string &name) const
template<typename AO >
shared_ptr< AO > getAnalysisObject (const std::string &name)
void removeAnalysisObject (const std::string &path)
 Unregister a data object from the histogram system (by name)
void removeAnalysisObject (AnalysisObjectPtr ao)
 Unregister a data object from the histogram system (by pointer)
const Histo1DPtr getHisto1D (const std::string &name) const
 Get a named Histo1D object from the histogram system.
Histo1DPtr getHisto1D (const std::string &name)
 Get a named Histo1D object from the histogram system (non-const)
const Histo1DPtr getHisto1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get a Histo1D object from the histogram system by axis ID codes (non-const)
Histo1DPtr getHisto1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId)
 Get a Histo1D object from the histogram system by axis ID codes (non-const)
const Profile1DPtr getProfile1D (const std::string &name) const
 Get a named Profile1D object from the histogram system.
Profile1DPtr getProfile1D (const std::string &name)
 Get a named Profile1D object from the histogram system (non-const)
const Profile1DPtr getProfile1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get a Profile1D object from the histogram system by axis ID codes (non-const)
Profile1DPtr getProfile1D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId)
 Get a Profile1D object from the histogram system by axis ID codes (non-const)
const Scatter2DPtr getScatter2D (const std::string &name) const
 Get a named Scatter2D object from the histogram system.
Scatter2DPtr getScatter2D (const std::string &name)
 Get a named Scatter2D object from the histogram system (non-const)
const Scatter2DPtr getScatter2D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId) const
 Get a Scatter2D object from the histogram system by axis ID codes (non-const)
Scatter2DPtr getScatter2D (unsigned int datasetId, unsigned int xAxisId, unsigned int yAxisId)
 Get a Scatter2D object from the histogram system by axis ID codes (non-const)
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

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

Private Attributes

bool _use_fiducial_lepton_efficiency
 Fiducial efficiencies to model the effects of the ATLAS detector.
vector< string > _signal_regions
 List of signal regions and event counts per signal region.
map< string, double > _eventCountsPerSR
Histo1DPtr _h_HTlep_all
 Histograms.
Histo1DPtr _h_HTjets_all
Histo1DPtr _h_MET_all
Histo1DPtr _h_Meff_all
Histo1DPtr _h_pt_1_3l
Histo1DPtr _h_pt_2_3l
Histo1DPtr _h_pt_3_3l
Histo1DPtr _h_pt_1_2ltau
Histo1DPtr _h_pt_2_2ltau
Histo1DPtr _h_pt_3_2ltau
Histo1DPtr _h_e_n
Histo1DPtr _h_mu_n
Histo1DPtr _h_tau_n
Histo1DPtr _h_excluded

Detailed Description

Definition at line 14 of file ATLAS_2012_I1204447.cc.


Constructor & Destructor Documentation

ATLAS_2012_I1204447 ( ) [inline]

Constructor.

Definition at line 18 of file ATLAS_2012_I1204447.cc.

      : Analysis("ATLAS_2012_I1204447") {  }

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

                                                                             {
    if (!_allowProjReg) {
      cerr << "Trying to register projection '"
           << proj.name() << "' before init phase in '" << this->name() << "'." << endl;
      exit(2);
    }
    const Projection& reg = getProjHandler().registerProjection(*this, proj, name);
    return reg;
  }
void addAnalysisObject ( AnalysisObjectPtr  ao) [protected, inherited]

Register a data object in the histogram system.

Todo:
2D versions of integrate... defined how, exactly?!?

Definition at line 721 of file Analysis.cc.

References Analysis::_analysisobjects.

Referenced by Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), and H1_1995_S3167097::finalize().

                                                       {
    _analysisobjects.push_back(ao);
  }
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 116 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

Referenced by FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), ATLAS_2010_S8894728::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMS_2011_S8950903::init(), CMS_2012_PAS_QCD_11_010::init(), CMS_2011_S8968497::init(), CMS_2012_I1193338::init(), CMS_2012_I1184941::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2014_I1298811::init(), ATLAS_2011_I894867::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), CMS_2011_S8941262::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), ATLAS_2010_CONF_2010_049::init(), CDF_2007_S7057202::init(), CMS_2011_S9086218::init(), TOTEM_2012_I1115294::init(), CMS_2011_S9215166::init(), CMS_2012_I1087342::init(), ATLAS_2014_I1282441::init(), ALICE_2011_S8945144::init(), ATLAS_2012_I1091481::init(), MC_MUONS::init(), MC_TAUS::init(), ATLAS_2010_S8591806::init(), ATLAS_2011_S9002537::init(), CMS_2011_I954992::init(), CMS_2011_S8957746::init(), CMS_2011_S8978280::init(), TOTEM_2012_002::init(), CMS_2011_S9088458::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), ATLAS_2011_I930220::init(), LHCF_2012_I1115479::init(), MC_ELECTRONS::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), ATLAS_2011_I925932::init(), MC_JETTAGS::init(), CDF_2012_NOTE10874::init(), ATLAS_2012_I1124167::init(), STAR_2006_S6500200::init(), STAR_2008_S7993412::init(), CMS_2013_I1256943::init(), CMS_QCD_10_024::init(), UA5_1987_S1640666::init(), CDF_1997_S3541940::init(), MC_HINC::init(), MC_JETS::init(), MC_KTSPLITTINGS::init(), CDF_2005_S6080774::init(), CDF_2006_S6450792::init(), SFM_1984_S1178091::init(), ARGUS_1993_S2789213::init(), ATLAS_2011_S9128077::init(), CMS_2013_I1261026::init(), CMS_2013_I1258128::init(), CMS_2014_I1298810::init(), D0_2000_S4480767::init(), BABAR_2007_S7266081::init(), BELLE_2008_I786560::init(), CDF_1993_S2742446::init(), CDF_2000_S4155203::init(), MC_WWINC::init(), CDF_2008_S7782535::init(), MC_ZJETS::init(), CDF_2008_S8093652::init(), MC_ZZINC::init(), STAR_2006_S6870392::init(), UA5_1982_S875503::init(), UA5_1989_S1926373::init(), ATLAS_2013_I1190187::init(), CMS_2013_I1209721::init(), D0_2008_S6879055::init(), ATLAS_2012_I1204447::init(), ARGUS_1993_S2669951::init(), LHCB_2012_I1208102::init(), ATLAS_2011_I944826::init(), MC_ZINC::init(), ALEPH_1991_S2435284::init(), MC_WINC::init(), MC_WJETS::init(), MC_WWKTSPLITTINGS::init(), MC_ZZJETS::init(), ATLAS_2012_I1183818::init(), ATLAS_2011_S9035664::init(), ATLAS_2012_I1204784::init(), CMS_2012_I1102908::init(), ATLAS_2013_I1219109::init(), CMS_2013_I1272853::init(), ATLAS_2013_I1243871::init(), ATLAS_2014_I1268975::init(), D0_2000_I499943::init(), ATLAS_2014_I1325553::init(), E735_1998_S3905616::init(), MC_ZZKTSPLITTINGS::init(), MC_DIJET::init(), MC_DIPHOTON::init(), ALEPH_2002_S4823664::init(), DELPHI_1999_S3960137::init(), ALEPH_1996_S3196992::init(), UA5_1986_S1583476::init(), EXAMPLE_CUTS::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), MC_ZKTSPLITTINGS::init(), MC_HJETS::init(), MC_HKTSPLITTINGS::init(), CDF_2005_S6217184::init(), MC_PHOTONINC::init(), MC_WKTSPLITTINGS::init(), ATLAS_2011_S8924791::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), ATLAS_2010_S8817804::init(), PDG_TAUS::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), D0_2009_S8202443::init(), D0_2010_S8821313::init(), ATLAS_2011_I945498::init(), MC_PHOTONKTSPLITTINGS::init(), MC_PHOTONJETS::init(), MC_WWJETS::init(), CMS_2012_I1107658::init(), D0_2008_S7554427::init(), D0_2008_S7863608::init(), D0_2010_S8671338::init(), MC_IDENTIFIED::init(), H1_1995_S3167097::init(), CDF_1994_S2952106::init(), ATLAS_2011_I954993::init(), ATLAS_2011_S9212353::init(), ATLAS_2014_I1306615::init(), ATLAS_2014_I1307756::init(), ATLAS_2012_CONF_2012_104::init(), ATLAS_2012_CONF_2012_105::init(), UA1_1990_S2044935::init(), ZEUS_2001_S4815815::init(), MC_LEADJETUE::init(), LHCB_2011_I919315::init(), MC_TTBAR::init(), ATLAS_2012_I1119557::init(), CDF_2008_S7828950::init(), MC_GENERIC::init(), ATLAS_2012_I1180197::init(), MC_PHOTONJETUE::init(), ATLAS_2012_I1188891::init(), CDF_2008_S7540469::init(), D0_2001_S4674421::init(), ATLAS_2014_I1300647::init(), D0_2008_S7837160::init(), ATLAS_2011_S9131140::init(), CDF_2008_S8095620::init(), DELPHI_1995_S3137023::init(), ATLAS_2012_CONF_2012_109::init(), ATLAS_2012_I1126136::init(), ATLAS_2011_S8983313::init(), ATLAS_2013_I1217863_Z::init(), STAR_2009_UE_HELEN::init(), OPAL_1997_S3396100::init(), STAR_2006_S6860818::init(), ATLAS_2011_CONF_2011_098::init(), JADE_1998_S3612880::init(), CMS_2012_I941555::init(), ATLAS_2012_I1095236::init(), JADE_OPAL_2000_S4300807::init(), D0_1996_S3214044::init(), ATLAS_2012_CONF_2012_103::init(), CDF_1996_S3108457::init(), EXAMPLE::init(), ATLAS_2011_I926145::init(), ALICE_2010_S8624100::init(), MC_QCD_PARTONS::init(), ATLAS_2012_I1186556::init(), ATLAS_2013_I1217863_W::init(), CMS_2014_I1303894::init(), ALEPH_2004_S5765862::init(), ALICE_2010_S8706239::init(), MC_HFJETS::init(), ATLAS_2014_I1288706::init(), ATLAS_2012_I1112263::init(), ATLAS_2011_S9212183::init(), ATLAS_2012_I1190891::init(), ATLAS_2012_CONF_2012_001::init(), CMS_2013_I1218372::init(), ATLAS_2012_I1125575::init(), ATLAS_2012_I1125961::init(), LHCB_2013_I1218996::init(), ALICE_2010_S8625980::init(), CDF_2009_NOTE_9936::init(), CDF_1996_S3349578::init(), MC_SUSY::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), CDF_2009_S8436959::init(), CDF_1998_S3618439::init(), ATLAS_2013_I1263495::init(), D0_2009_S8320160::init(), CDF_2000_S4266730::init(), D0_2009_S8349509::init(), CDF_2001_S4563131::init(), CDF_2001_S4517016::init(), D0_2006_S6438750::init(), ALEPH_2001_S4656318::init(), SLD_2002_S4869273::init(), ATLAS_2010_S8919674::init(), ATLAS_2012_I1082936::init(), D0_2010_S8566488::init(), DELPHI_2002_069_CONF_603::init(), LHCB_2014_I1281685::init(), LHCB_2011_I917009::init(), UA5_1988_S1867512::init(), SLD_1996_S3398250::init(), MC_PHOTONS::init(), ATLAS_2011_I921594::init(), CDF_1996_S3418421::init(), LHCB_2012_I1119400::init(), ATLAS_2011_S9108483::init(), MC_WPOL::init(), CDF_2006_S6653332::init(), CDF_2008_S7541902::init(), ATLAS_2014_I1306294::init(), ATLAS_2014_I1304688::init(), D0_2008_S7662670::init(), ATLAS_2013_I1230812::init(), ATLAS_2012_I1199269::init(), ATLAS_2014_I1307243::init(), DELPHI_2000_S4328825::init(), ATLAS_2012_I943401::init(), OPAL_2002_S5361494::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I946427::init(), ATLAS_2012_I1117704::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2011_S9225137::init(), ATLAS_2012_I1082009::init(), ATLAS_2012_CONF_2012_153::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_WJET::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_I1083318::init(), CDF_2009_S8233977::init(), ATLAS_2013_I1217867::init(), CMS_2013_I1224539_ZJET::init(), CDF_2010_S8591881_DY::init(), CDF_2010_S8591881_QCD::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2014_I1312627::init(), ATLAS_2010_S8918562::init(), ATLAS_2011_S9120807::init(), CDF_2001_S4751469::init(), D0_2008_S7719523::init(), ATLAS_2011_S9041966::init(), ATLAS_2012_I1094568::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), ATLAS_2011_S9126244::init(), BABAR_2007_S6895344::init(), MC_VH2BB::init(), BABAR_2005_S6181155::init(), BELLE_2001_S4598261::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), ATLAS_2012_I1203852::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ATLAS_2012_I1094061::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), ATLAS_2012_I1093734::init(), OPAL_1993_S2692198::init(), ATLAS_2014_I1279489::init(), ATLAS_2012_I1094564::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), TauFinder::TauFinder(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder().

                                                                       {
      const Projection& reg = _addProjection(proj, name);
      const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
      return rtn;
    }
const vector<AnalysisObjectPtr>& analysisObjects ( ) const [inline, inherited]

List of registered analysis data objects.

Definition at line 694 of file Analysis.hh.

References Analysis::_analysisobjects.

Referenced by Analysis::getAnalysisObject().

                                                             {
      return _analysisobjects;
    }
void analyze ( const Event event) [inline, virtual]

Perform the per-event analysis.

Todo:
This could benefit from a tau finder projection
Todo:
Huh? Unstable taus with no decay vtx? Can use Particle.isStable()? But why in this situation?

Implements Analysis.

Definition at line 88 of file ATLAS_2012_I1204447.cc.

References ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_use_fiducial_lepton_efficiency, ParticleBase::abseta(), FourVector::abseta(), Particle::abspid(), ATLAS_2012_I1204447::apply_reco_eff(), ParticleBase::azimuthalAngle(), Rivet::deltaR(), Rivet::Cuts::eta, ATLAS_2012_I1204447::fillEventCountsPerSR(), Particle::genParticle(), ATLAS_2012_I1204447::get_prong_number(), ATLAS_2012_I1204447::get_tau_neutrino_mom(), Rivet::GeV, Particle::hasAncestor(), Rivet::inRange(), ATLAS_2012_I1204447::isonZ(), Rivet::MINUSPI_PLUSPI, Particle::momentum(), Jet::momentum(), MSG_DEBUG, Rivet::PID::MUON, Rivet::particles(), ParticleBase::perp(), Rivet::Cuts::phi, Rivet::Cuts::pT, ParticleBase::pT(), FourMomentum::pT(), ParticleBase::rapidity(), Rivet::sortByPt(), Rivet::PID::TAU, and vetoEvent.

                                     {

      // Muons
      Particles muon_candidates;
      const Particles charged_tracks    = applyProjection<ChargedFinalState>(event, "CFS").particles();
      const Particles visible_particles = applyProjection<VisibleFinalState>(event, "VFS").particles();
      foreach (const Particle& mu, applyProjection<IdentifiedFinalState>(event, "muons").particlesByPt()) {
        // Calculate pTCone30 variable (pT of all tracks within dR<0.3 - pT of muon itself)
        double pTinCone = -mu.pT();
        foreach (const Particle& track, charged_tracks) {
          if (deltaR(mu.momentum(), track.momentum()) < 0.3)
            pTinCone += track.pT();
        }

        // Calculate eTCone30 variable (pT of all visible particles within dR<0.3)
        double eTinCone = 0.;
        foreach (const Particle& visible_particle, visible_particles) {
          if (visible_particle.abspid() != PID::MUON && inRange(deltaR(mu.momentum(), visible_particle.momentum()), 0.1, 0.3))
            eTinCone += visible_particle.pT();
        }

        // Apply reconstruction efficiency and simulate reco
        int muon_id = 13;
        if ( mu.hasAncestor(15) || mu.hasAncestor(-15)) muon_id = 14;
        const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(muon_id, mu) : 1.0;
        const bool keep_muon = rand()/static_cast<double>(RAND_MAX) <= eff;

        // Keep muon if pTCone30/pT < 0.15 and eTCone30/pT < 0.2 and reconstructed
        if (keep_muon && pTinCone/mu.pT() <= 0.15 && eTinCone/mu.pT() < 0.2)
          muon_candidates.push_back(mu);
      }


      // Electrons
      Particles electron_candidates;
      foreach (const Particle& e, applyProjection<IdentifiedFinalState>(event, "elecs").particlesByPt()) {
        // Neglect electrons in crack regions
        if (inRange(e.abseta(), 1.37, 1.52)) continue;

        // Calculate pTCone30 variable (pT of all tracks within dR<0.3 - pT of electron itself)
        double pTinCone = -e.pT();
        foreach (const Particle& track, charged_tracks) {
          if (deltaR(e.momentum(), track.momentum()) < 0.3) pTinCone += track.pT();
        }

        // Calculate eTCone30 variable (pT of all visible particles (except muons) within dR<0.3)
        double eTinCone = 0.;
        foreach (const Particle& visible_particle, visible_particles) {
          if (visible_particle.abspid() != PID::MUON && inRange(deltaR(e.momentum(), visible_particle.momentum()), 0.1, 0.3))
            eTinCone += visible_particle.pT();
        }

        // Apply reconstruction efficiency and simulate reco
        int elec_id = 11;
        if (e.hasAncestor(15) || e.hasAncestor(-15)) elec_id = 12;
        const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(elec_id, e) : 1.0;
        const bool keep_elec = rand()/static_cast<double>(RAND_MAX) <= eff;

        // Keep electron if pTCone30/pT < 0.13 and eTCone30/pT < 0.2 and reconstructed
        if (keep_elec && pTinCone/e.pT() <= 0.13 && eTinCone/e.pT() < 0.2)
          electron_candidates.push_back(e);
      }


      // Taus
      /// @todo This could benefit from a tau finder projection
      Particles tau_candidates;
      foreach (const Particle& tau, applyProjection<UnstableFinalState>(event, "UFS").particlesByPt()) {
        // Only pick taus out of all unstable particles
        if (tau.abspid() != PID::TAU) continue;

        // Check that tau has decayed into daughter particles
        /// @todo Huh? Unstable taus with no decay vtx? Can use Particle.isStable()? But why in this situation?
        if (tau.genParticle()->end_vertex() == 0) continue;

        // Calculate visible tau pT from pT of tau neutrino in tau decay for pT and |eta| cuts
        const FourMomentum tau_neutrino_mom = get_tau_neutrino_mom(tau);
        const FourMomentum visible_tau_mom = tau.momentum() - tau_neutrino_mom;
        if (visible_tau_mom.pT() < 15*GeV || visible_tau_mom.abseta() > 2.5) continue;

        // Get prong number (number of tracks) in tau decay and check if tau decays leptonically
        unsigned int nprong = 0;
        bool lep_decaying_tau = false;
        get_prong_number(tau.genParticle(), nprong, lep_decaying_tau);

        // Apply reconstruction efficiency
        int tau_id = 15;
        if (nprong == 1) tau_id = 15;
        else if (nprong == 3) tau_id = 16;

        // Get fiducial lepton efficiency simulate reco efficiency
        const double eff = (_use_fiducial_lepton_efficiency) ? apply_reco_eff(tau_id, tau) : 1.0;
        const bool keep_tau = rand()/static_cast<double>(RAND_MAX) <= eff;

        // Keep tau if nprong = 1, it decays hadronically, and it's reconstructed by the detector
        if ( !lep_decaying_tau && nprong == 1 && keep_tau) tau_candidates.push_back(tau);
      }


      // Jets (all anti-kt R=0.4 jets with pT > 25 GeV and eta < 4.9)
      Jets jet_candidates;
      foreach (const Jet& jet, applyProjection<FastJets>(event, "AntiKtJets04").jetsByPt(25*GeV)) {
        if (jet.abseta() < 4.9) jet_candidates.push_back(jet);
      }


      // ETmiss
      Particles vfs_particles = applyProjection<VisibleFinalState>(event, "VFS").particles();
      FourMomentum pTmiss;
      foreach (const Particle& p, vfs_particles) pTmiss -= p.momentum();
      double eTmiss = pTmiss.pT()/GeV;


      //------------------
      // Overlap removal

      // electron - electron
      Particles electron_candidates_2;
      for (size_t ie = 0; ie < electron_candidates.size(); ++ie) {
        const Particle & e = electron_candidates[ie];
        bool away = true;
        // If electron pair within dR < 0.1: remove electron with lower pT
        for (size_t ie2=0; ie2 < electron_candidates_2.size(); ++ie2) {
          if ( deltaR( e.momentum(), electron_candidates_2[ie2].momentum()) < 0.1 ) {
            away = false;
            break;
          }
        }
        // If isolated keep it
        if ( away )
          electron_candidates_2.push_back( e );
      }
      // jet - electron
      Jets recon_jets;
      foreach (const Jet& jet, jet_candidates) {
        bool away = true;
        // if jet within dR < 0.2 of electron: remove jet
        foreach (const Particle& e, electron_candidates_2) {
          if (deltaR(e.momentum(), jet.momentum()) < 0.2) {
            away = false;
            break;
          }
        }
        // jet - tau
        if (away)  {
          // If jet within dR < 0.2 of tau: remove jet
          foreach (const Particle& tau, tau_candidates) {
            if (deltaR(tau.momentum(), jet.momentum()) < 0.2) {
              away = false;
              break;
            }
          }
        }
        // If isolated keep it
        if ( away )
          recon_jets.push_back( jet );
      }


      // electron - jet
      Particles recon_leptons, recon_e;
      for (size_t ie = 0; ie < electron_candidates_2.size(); ++ie) {
        const Particle& e = electron_candidates_2[ie];
        // If electron within 0.2 < dR < 0.4 from any jets: remove electron
        bool away = true;
        foreach (const Jet& jet, recon_jets) {
          if (deltaR(e.momentum(), jet.momentum()) < 0.4) {
            away = false;
            break;
          }
        }
        // electron - muon
        // if electron within dR < 0.1 of a muon: remove electron
        if (away) {
          foreach (const Particle& mu, muon_candidates) {
            if (deltaR(mu.momentum(), e.momentum()) < 0.1) {
              away = false;
              break;
            }
          }
        }
        // If isolated keep it
        if (away)  {
          recon_e += e;
          recon_leptons += e;
        }
      }


      // tau - electron
      Particles recon_tau;
      foreach ( const Particle& tau, tau_candidates ) {
        bool away = true;
        // If tau within dR < 0.2 of an electron: remove tau
        foreach ( const Particle& e, recon_e ) {
          if (deltaR( tau.momentum(), e.momentum()) < 0.2) {
            away = false;
            break;
          }
        }
        // tau - muon
        // If tau within dR < 0.2 of a muon: remove tau
        if (away)  {
          foreach (const Particle& mu, muon_candidates) {
            if (deltaR(tau.momentum(), mu.momentum()) < 0.2) {
              away = false;
              break;
            }
          }
        }
        // If isolated keep it
        if (away) recon_tau.push_back( tau );
      }

      // Muon - jet isolation
      Particles recon_mu, trigger_mu;
      // If muon within dR < 0.4 of a jet, remove muon
      foreach (const Particle& mu, muon_candidates) {
        bool away = true;
        foreach (const Jet& jet, recon_jets) {
          if ( deltaR( mu.momentum(), jet.momentum()) < 0.4 ) {
            away = false;
            break;
          }
        }
        if (away) {
          recon_mu.push_back( mu );
          recon_leptons.push_back( mu );
          if (mu.abseta() < 2.4) trigger_mu.push_back( mu );
        }
      }

      // End overlap removal
      //------------------


      // Jet cleaning
      if (rand()/static_cast<double>(RAND_MAX) <= 0.42) {
        foreach (const Jet& jet, recon_jets) {
          const double eta = jet.rapidity();
          const double phi = jet.azimuthalAngle(MINUSPI_PLUSPI);
          if (jet.pT() > 25*GeV && inRange(eta, -0.1, 1.5) && inRange(phi, -0.9, -0.5)) vetoEvent;
        }
      }


      // Post-isolation event cuts
      // Require at least 3 charged tracks in event
      if (charged_tracks.size() < 3) vetoEvent;

      // And at least one e/mu passing trigger
      if (!( !recon_e   .empty() && recon_e[0]   .pT() > 25*GeV)  &&
          !( !trigger_mu.empty() && trigger_mu[0].pT() > 25*GeV) ) {
        MSG_DEBUG("Hardest lepton fails trigger");
        vetoEvent;
      }

      // And only accept events with at least 2 electrons and muons and at least 3 leptons in total
      if (recon_mu.size() + recon_e.size() + recon_tau.size() < 3 || recon_leptons.size() < 2) vetoEvent;

      // Now it's worth getting the event weight
      const double weight = event.weight();

      // Sort leptons by decreasing pT
      sortByPt(recon_leptons);
      sortByPt(recon_tau);

      // Calculate HTlep, fill lepton pT histograms & store chosen combination of 3 leptons
      double HTlep = 0.;
      Particles chosen_leptons;
      if ( recon_leptons.size() > 2 ) {
        _h_pt_1_3l->fill(recon_leptons[0].perp()/GeV, weight);
        _h_pt_2_3l->fill(recon_leptons[1].perp()/GeV, weight);
        _h_pt_3_3l->fill(recon_leptons[2].perp()/GeV, weight);
        HTlep = (recon_leptons[0].pT() + recon_leptons[1].pT() + recon_leptons[2].pT())/GeV;
        chosen_leptons.push_back( recon_leptons[0] );
        chosen_leptons.push_back( recon_leptons[1] );
        chosen_leptons.push_back( recon_leptons[2] );
      }
      else {
        _h_pt_1_2ltau->fill(recon_leptons[0].perp()/GeV, weight);
        _h_pt_2_2ltau->fill(recon_leptons[1].perp()/GeV, weight);
        _h_pt_3_2ltau->fill(recon_tau[0].perp()/GeV,     weight);
        HTlep = (recon_leptons[0].pT() + recon_leptons[1].pT() + recon_tau[0].pT())/GeV ;
        chosen_leptons.push_back( recon_leptons[0] );
        chosen_leptons.push_back( recon_leptons[1] );
        chosen_leptons.push_back( recon_tau[0] );
      }

      // Number of prompt e/mu and had taus
      _h_e_n  ->fill(recon_e.size()  , weight);
      _h_mu_n ->fill(recon_mu.size() , weight);
      _h_tau_n->fill(recon_tau.size(), weight);

      // Calculate HTjets
      double HTjets = 0.;
      foreach ( const Jet & jet, recon_jets )
        HTjets += jet.perp()/GeV;

      // Calculate meff
      double meff = eTmiss + HTjets;
      Particles all_leptons;
      foreach ( const Particle & e , recon_e  )  {
        meff += e.perp()/GeV;
        all_leptons.push_back( e );
      }
      foreach ( const Particle & mu, recon_mu )  {
        meff += mu.perp()/GeV;
        all_leptons.push_back( mu );
      }
      foreach ( const Particle & tau, recon_tau )  {
        meff += tau.perp()/GeV;
        all_leptons.push_back( tau );
      }

      // Fill histogram of kinematic variables
      _h_HTlep_all ->fill(HTlep , weight);
      _h_HTjets_all->fill(HTjets, weight);
      _h_MET_all   ->fill(eTmiss, weight);
      _h_Meff_all  ->fill(meff  , weight);

      // Determine signal region (3l/2ltau, onZ/offZ)
      string basic_signal_region;
      if ( recon_mu.size() + recon_e.size() > 2 )
        basic_signal_region += "3l_";
      else if ( (recon_mu.size() + recon_e.size() == 2) && (recon_tau.size() > 0))
        basic_signal_region += "2ltau_";
      // Is there an OSSF pair or a three lepton combination with an invariant mass close to the Z mass
      int onZ = isonZ(chosen_leptons);
      if      (onZ == 1)   basic_signal_region += "onZ";
      else if (onZ == 0)   basic_signal_region += "offZ";
      // Check in which signal regions this event falls and adjust event counters
      fillEventCountsPerSR(basic_signal_region, onZ, HTlep, eTmiss, HTjets, meff, weight);
    }
double apply_reco_eff ( int  flavor,
const Particle p 
) [inline]

Function giving fiducial lepton efficiency.

Definition at line 659 of file ATLAS_2012_I1204447.cc.

References Rivet::Cuts::eta, ParticleBase::eta(), Rivet::GeV, Rivet::Cuts::pt, and ParticleBase::pT().

Referenced by ATLAS_2012_I1204447::analyze().

                                                         {
      float pt = p.pT()/GeV;
      float eta = p.eta();

      double eff = 0.;
      //double err = 0.;

      if (flavor == 11) { // weight prompt electron -- now including data/MC ID SF in eff.
        //float rho = 0.820;
        float p0 = 7.34;  float p1 = 0.8977;
        //float ep0= 0.5 ;  float ep1= 0.0087;
        eff = p1 - p0/pt;

        //double err0 = ep0/pt; // d(eff)/dp0
        //double err1 = ep1;    // d(eff)/dp1
        //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1);

        double avgrate = 0.6867;
        float wz_ele_eta[] = {0.588717,0.603674,0.666135,0.747493,0.762202,0.675051,0.751606,0.745569,0.665333,0.610432,0.592693,};
        //float ewz_ele_eta[] ={0.00292902,0.002476,0.00241209,0.00182319,0.00194339,0.00299785,0.00197339,0.00182004,0.00241793,0.00245997,0.00290394,};
        int ibin = 3;

        if (eta >= -2.5 && eta < -2.0) ibin = 0;
        if (eta >= -2.0 && eta < -1.5) ibin = 1;
        if (eta >= -1.5 && eta < -1.0) ibin = 2;
        if (eta >= -1.0 && eta < -0.5) ibin = 3;
        if (eta >= -0.5 && eta < -0.1) ibin = 4;
        if (eta >= -0.1 && eta <  0.1) ibin = 5;
        if (eta >=  0.1 && eta <  0.5) ibin = 6;
        if (eta >=  0.5 && eta <  1.0) ibin = 7;
        if (eta >=  1.0 && eta <  1.5) ibin = 8;
        if (eta >=  1.5 && eta <  2.0) ibin = 9;
        if (eta >=  2.0 && eta <  2.5) ibin = 10;

        double eff_eta = wz_ele_eta[ibin];
        //double err_eta = ewz_ele_eta[ibin];

        eff = (eff*eff_eta)/avgrate;
      }

      if (flavor == 12)  { // weight electron from tau
        //float rho = 0.884;
        float p0 = 6.799;  float p1 = 0.842;
        //float ep0= 0.664;  float ep1= 0.016;
        eff = p1 - p0/pt;

        //double err0 = ep0/pt; // d(eff)/dp0
        //double err1 = ep1;    // d(eff)/dp1
        //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1);

        double avgrate = 0.5319;
        float wz_elet_eta[] = {0.468945,0.465953,0.489545,0.58709,0.59669,0.515829,0.59284,0.575828,0.498181,0.463536,0.481738,};
        //float ewz_elet_eta[] ={0.00933795,0.00780868,0.00792679,0.00642083,0.00692652,0.0101568,0.00698452,0.00643524,0.0080002,0.00776238,0.0094699,};
        int ibin = 3;

        if (eta >= -2.5 && eta < -2.0) ibin = 0;
        if (eta >= -2.0 && eta < -1.5) ibin = 1;
        if (eta >= -1.5 && eta < -1.0) ibin = 2;
        if (eta >= -1.0 && eta < -0.5) ibin = 3;
        if (eta >= -0.5 && eta < -0.1) ibin = 4;
        if (eta >= -0.1 && eta <  0.1) ibin = 5;
        if (eta >=  0.1 && eta <  0.5) ibin = 6;
        if (eta >=  0.5 && eta <  1.0) ibin = 7;
        if (eta >=  1.0 && eta <  1.5) ibin = 8;
        if (eta >=  1.5 && eta <  2.0) ibin = 9;
        if (eta >=  2.0 && eta <  2.5) ibin = 10;

        double eff_eta = wz_elet_eta[ibin];
        //double err_eta = ewz_elet_eta[ibin];

        eff = (eff*eff_eta)/avgrate;

      }

      if (flavor == 13)  {// weight prompt muon

        //if eta>0.1
        float p0 = -18.21;  float p1 = 14.83;  float p2 = 0.9312;
        //float ep0= 5.06;    float ep1= 1.9;    float ep2=0.00069;

        if ( fabs(eta) < 0.1)  {
          p0  = 7.459; p1 = 2.615; p2  = 0.5138;
          //ep0 = 10.4; ep1 = 4.934; ep2 = 0.0034;
        }

        double arg = ( pt-p0 )/( 2.*p1 ) ;
        eff = 0.5 * p2 * (1.+erf(arg));
        //err = 0.1*eff;
      }

      if (flavor == 14)  {// weight muon from tau

        if (fabs(eta) < 0.1) {
          float p0 = -1.756;  float p1 = 12.38;  float p2 = 0.4441;
          //float ep0= 10.39;   float ep1= 7.9;  float ep2=0.022;
          double arg = ( pt-p0 )/( 2.*p1 ) ;
          eff = 0.5 * p2 * (1.+erf(arg));
          //err = 0.1*eff;
        }
        else {
          float p0 = 2.102;  float p1 = 0.8293;
          //float ep0= 0.271;  float ep1= 0.0083;
          eff = p1 - p0/pt;
          //double err0 = ep0/pt; // d(eff)/dp0
          //double err1 = ep1;    // d(eff)/dp1
          //err = sqrt(err0*err0 + err1*err1 - 2*rho*err0*err1);
        }
      }

      if (flavor == 15)  {// weight hadronic tau 1p

        float wz_tau1p[] = {0.0249278,0.146978,0.225049,0.229212,0.21519,0.206152,0.201559,0.197917,0.209249,0.228336,0.193548,};
        //float ewz_tau1p[] ={0.00178577,0.00425252,0.00535052,0.00592126,0.00484684,0.00612941,0.00792099,0.0083006,0.0138307,0.015568,0.0501751,};
        int ibin = 0;
        if (pt > 15)  ibin = 1;
        if (pt > 20)  ibin = 2;
        if (pt > 25)  ibin = 3;
        if (pt > 30)  ibin = 4;
        if (pt > 40)  ibin = 5;
        if (pt > 50)  ibin = 6;
        if (pt > 60)  ibin = 7;
        if (pt > 80)  ibin = 8;
        if (pt > 100) ibin = 9;
        if (pt > 200) ibin = 10;

        eff = wz_tau1p[ibin];
        //err = ewz_tau1p[ibin];


        double avgrate = 0.1718;
        float wz_tau1p_eta[] = {0.162132,0.176393,0.139619,0.178813,0.185144,0.210027,0.203937,0.178688,0.137034,0.164216,0.163713,};
        //float ewz_tau1p_eta[] ={0.00706705,0.00617989,0.00506798,0.00525172,0.00581865,0.00865675,0.00599245,0.00529877,0.00506368,0.00617025,0.00726219,};

        ibin = 3;
        if (eta >= -2.5 && eta < -2.0) ibin = 0;
        if (eta >= -2.0 && eta < -1.5) ibin = 1;
        if (eta >= -1.5 && eta < -1.0) ibin = 2;
        if (eta >= -1.0 && eta < -0.5) ibin = 3;
        if (eta >= -0.5 && eta < -0.1) ibin = 4;
        if (eta >= -0.1 && eta <  0.1) ibin = 5;
        if (eta >=  0.1 && eta <  0.5) ibin = 6;
        if (eta >=  0.5 && eta <  1.0) ibin = 7;
        if (eta >=  1.0 && eta <  1.5) ibin = 8;
        if (eta >=  1.5 && eta <  2.0) ibin = 9;
        if (eta >=  2.0 && eta <  2.5) ibin = 10;

        double eff_eta = wz_tau1p_eta[ibin];
        //double err_eta = ewz_tau1p_eta[ibin];

        eff = (eff*eff_eta)/avgrate;
      }

      if (flavor == 16)  { //weight hadronic tau 3p

        float wz_tau3p[] = {0.000587199,0.00247181,0.0013031,0.00280112,};
        //float ewz_tau3p[] ={0.000415091,0.000617187,0.000582385,0.00197792,};

        int ibin = 0;
        if (pt > 15) ibin = 1;
        if (pt > 20) ibin = 2;
        if (pt > 40) ibin = 3;
        if (pt > 80) ibin = 4;

        eff = wz_tau3p[ibin];
        //err = ewz_tau3p[ibin];
      }

      return eff;
    }
const PROJ& applyProjection ( const Event evt,
const PROJ &  proj 
) const [inline, inherited]

Apply the supplied projection on event.

Definition at line 70 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

Referenced by DISFinalState::project().

                                                                          {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }
const PROJ& applyProjection ( const Event evt,
const Projection proj 
) const [inline, inherited]

Apply the supplied projection on event.

Definition at line 77 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                                {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }
const PROJ& applyProjection ( const Event evt,
const std::string &  name 
) const [inline, inherited]

Apply the named projection on event.

Definition at line 84 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                               {
      return pcast<PROJ>(_applyProjection(evt, name));
    }
void asymm ( Histo1DPtr  h1,
Histo1DPtr  h2,
Scatter2DPtr  s 
) const [inherited]

Helper for histogram asymmetry calculation.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.

Definition at line 639 of file Analysis.cc.

                                                                         {
    const string path = s->path();
    *s = YODA::asymm(*h1, *h2);
    s->setPath(path);
  }
void asymm ( const YODA::Histo1D &  h1,
const YODA::Histo1D &  h2,
Scatter2DPtr  s 
) const [inherited]

Helper for histogram asymmetry calculation.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
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 132 of file Analysis.hh.

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

                                                 {
      return info().authors();
    }
const PdgIdPair beamIds ( ) const [inherited]

Incoming beam IDs for this run.

Definition at line 35 of file Analysis.cc.

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

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

                                          {
    return handler().beamIds();
  }
virtual std::string bibKey ( ) const [inline, virtual, inherited]

BibTeX citation key for this article.

Definition at line 185 of file Analysis.hh.

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

                                     {
      return info().bibKey();
    }
virtual std::string bibTeX ( ) const [inline, virtual, inherited]

BibTeX citation entry for this article.

Definition at line 190 of file Analysis.hh.

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

                                     {
      return info().bibTeX();
    }
Histo1DPtr bookHisto1D ( 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 .

Definition at line 200 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by ATLAS_2014_I1304688::_histogram(), Analysis::bookHisto1D(), ATLAS_2012_I1118269::init(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), LHCB_2010_I867355::init(), ALICE_2012_I1181770::init(), CMS_2011_S8950903::init(), ATLAS_2010_CONF_2010_049::init(), CMS_2011_S8941262::init(), CMS_2011_S8968497::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), CMS_2012_I1184941::init(), CMS_2012_I1193338::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2014_I1298811::init(), ATLAS_2011_I894867::init(), ALICE_2011_S8909580::init(), TOTEM_2012_I1115294::init(), CMS_2011_S9086218::init(), CMS_2012_I1087342::init(), CMS_2011_S9215166::init(), ATLAS_2014_I1282441::init(), ALICE_2011_S8945144::init(), ATLAS_2012_I1091481::init(), CDF_2007_S7057202::init(), ATLAS_2010_S8591806::init(), CMS_2011_I954992::init(), CMS_2011_S8957746::init(), CMS_2011_S8978280::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), LHCF_2012_I1115479::init(), ATLAS_2011_I930220::init(), TOTEM_2012_002::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), MC_JETTAGS::init(), ATLAS_2011_I925932::init(), STAR_2006_S6500200::init(), ATLAS_2012_I1124167::init(), UA5_1987_S1640666::init(), CMS_2013_I1256943::init(), CMS_QCD_10_024::init(), CDF_1997_S3541940::init(), SFM_1984_S1178091::init(), ARGUS_1993_S2789213::init(), CMS_2013_I1261026::init(), CMS_2013_I1258128::init(), CMS_2014_I1298810::init(), D0_2000_S4480767::init(), ATLAS_2011_S9128077::init(), BABAR_2007_S7266081::init(), BELLE_2008_I786560::init(), MC_HINC::init(), CDF_1993_S2742446::init(), CDF_2000_S4155203::init(), CDF_2006_S6450792::init(), CDF_2005_S6080774::init(), LHCB_2012_I1208102::init(), MC_JetSplittings::init(), MC_PDFS::init(), MC_WWINC::init(), MC_ZJETS::init(), MC_ZZINC::init(), CDF_2008_S8093652::init(), STAR_2006_S6870392::init(), UA5_1982_S875503::init(), UA5_1989_S1926373::init(), ARGUS_1993_S2669951::init(), ATLAS_2012_I1204447::init(), CMS_2013_I1209721::init(), ATLAS_2013_I1190187::init(), D0_2008_S6879055::init(), ATLAS_2011_I944826::init(), ALEPH_1991_S2435284::init(), ATLAS_2012_I1183818::init(), CMS_2012_I1102908::init(), ATLAS_2011_S9035664::init(), ATLAS_2012_I1204784::init(), CMS_2013_I1272853::init(), ATLAS_2013_I1219109::init(), D0_2000_I499943::init(), ATLAS_2014_I1268975::init(), ATLAS_2014_I1325553::init(), E735_1998_S3905616::init(), MC_DIJET::init(), MC_DIPHOTON::init(), MC_WINC::init(), MC_WJETS::init(), MC_ZINC::init(), MC_ZZJETS::init(), MC_ParticleAnalysis::init(), PDG_TAUS::init(), UA5_1986_S1583476::init(), ALEPH_1996_S3196992::init(), ATLAS_2010_S8817804::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), D0_2009_S8202443::init(), D0_2010_S8821313::init(), DELPHI_1999_S3960137::init(), EXAMPLE_CUTS::init(), ALEPH_2002_S4823664::init(), MC_HJETS::init(), CDF_1988_S1865951::init(), MC_JetAnalysis::init(), CDF_1990_S2089246::init(), MC_PHOTONINC::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), ATLAS_2011_I945498::init(), ATLAS_2011_I954993::init(), CMS_2012_I1107658::init(), D0_2008_S7554427::init(), D0_2008_S7863608::init(), D0_2010_S8671338::init(), H1_1995_S3167097::init(), MC_IDENTIFIED::init(), CDF_1994_S2952106::init(), MC_PHOTONJETS::init(), MC_WWJETS::init(), CDF_2008_S7828950::init(), ATLAS_2012_I1119557::init(), ZEUS_2001_S4815815::init(), ATLAS_2012_I1180197::init(), ATLAS_2012_I1188891::init(), MC_PHOTONJETUE::init(), D0_2001_S4674421::init(), ATLAS_2011_S9131140::init(), D0_2008_S7837160::init(), ATLAS_2011_S9212353::init(), ATLAS_2014_I1306615::init(), CDF_2008_S7540469::init(), ATLAS_2012_CONF_2012_105::init(), ATLAS_2012_CONF_2012_104::init(), LHCB_2011_I919315::init(), ATLAS_2014_I1300647::init(), UA1_1990_S2044935::init(), MC_GENERIC::init(), MC_TTBAR::init(), OPAL_1997_S3396100::init(), ATLAS_2012_I1095236::init(), ALICE_2010_S8624100::init(), ATLAS_2012_I1126136::init(), ATLAS_2012_CONF_2012_103::init(), JADE_1998_S3612880::init(), CDF_1996_S3108457::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_CONF_2011_098::init(), ATLAS_2012_I1186556::init(), JADE_OPAL_2000_S4300807::init(), STAR_2006_S6860818::init(), D0_1996_S3214044::init(), CMS_2012_I941555::init(), CDF_2008_S8095620::init(), EXAMPLE::init(), ATLAS_2013_I1217863_W::init(), ALEPH_2004_S5765862::init(), ATLAS_2011_I926145::init(), CMS_2014_I1303894::init(), ATLAS_2012_CONF_2012_109::init(), MC_XS::init(), ATLAS_2013_I1217863_Z::init(), DELPHI_1995_S3137023::init(), ATLAS_2012_I1190891::init(), ATLAS_2012_I1112263::init(), ATLAS_2012_I1125961::init(), ALICE_2010_S8625980::init(), ATLAS_2012_CONF_2012_001::init(), CMS_2013_I1218372::init(), ALICE_2010_S8706239::init(), ATLAS_2014_I1288706::init(), ATLAS_2011_S9212183::init(), LHCB_2013_I1218996::init(), MC_HFJETS::init(), CDF_1996_S3349578::init(), CDF_2009_NOTE_9936::init(), ATLAS_2013_I1263495::init(), CDF_2009_S8436959::init(), CDF_2000_S4266730::init(), MC_SUSY::init(), D0_2009_S8320160::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), D0_2009_S8349509::init(), CDF_2001_S4563131::init(), CDF_1998_S3618439::init(), CDF_2001_S4517016::init(), ATLAS_2012_I1082936::init(), ALEPH_2001_S4656318::init(), DELPHI_2002_069_CONF_603::init(), D0_2006_S6438750::init(), D0_2010_S8566488::init(), ATLAS_2010_S8919674::init(), SLD_2002_S4869273::init(), SLD_1996_S3398250::init(), MC_PHOTONS::init(), CDF_1996_S3418421::init(), ATLAS_2011_I921594::init(), LHCB_2012_I1119400::init(), LHCB_2014_I1281685::init(), CDF_2008_S7541902::init(), CDF_2006_S6653332::init(), ATLAS_2011_S9108483::init(), D0_2008_S7662670::init(), MC_WPOL::init(), ATLAS_2014_I1306294::init(), ATLAS_2013_I1230812::init(), ATLAS_2014_I1307243::init(), ATLAS_2012_I943401::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I946427::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2011_S9225137::init(), ATLAS_2012_I1117704::init(), ATLAS_2012_I1199269::init(), ATLAS_2011_S9019561::init(), ATLAS_2012_I1082009::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_CONF_2012_153::init(), ATLAS_2012_I1083318::init(), CMS_2013_I1224539_WJET::init(), CDF_2009_S8233977::init(), ATLAS_2013_I1217867::init(), CMS_2013_I1224539_ZJET::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2014_I1312627::init(), CDF_2001_S4751469::init(), ATLAS_2010_S8918562::init(), ATLAS_2011_S9120807::init(), D0_2008_S7719523::init(), ATLAS_2011_S9041966::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), BABAR_2007_S6895344::init(), MC_VH2BB::init(), BABAR_2005_S6181155::init(), BELLE_2001_S4598261::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), ATLAS_2012_I1203852::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), OPAL_1993_S2692198::init(), ATLAS_2012_I1094564::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), ATLAS_2014_I1279489::initialisePlots(), ATLAS_2012_I1125575::initializeHistograms(), ATLAS_2012_I1094568::initializePlots(), and ATLAS_2011_S9126244::initializePlots().

                                                         {
    const string path = histoPath(hname);
    Histo1DPtr hist( new Histo1D(nbins, lower, upper, path, title) );
    addAnalysisObject(hist);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setAnnotation("XLabel", xtitle);
    hist->setAnnotation("YLabel", ytitle);
    return hist;
  }
Histo1DPtr bookHisto1D ( 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 .

Definition at line 215 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                         {
    const string path = histoPath(hname);
    Histo1DPtr hist( new Histo1D(binedges, path, title) );
    addAnalysisObject(hist);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setAnnotation("XLabel", xtitle);
    hist->setAnnotation("YLabel", ytitle);
    return hist;
  }
Histo1DPtr bookHisto1D ( const std::string &  name,
const Scatter2D &  refscatter,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D histogram with binning from a reference scatter.

Definition at line 230 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                         {
    const string path = histoPath(hname);
    Histo1DPtr hist( new Histo1D(refscatter, path) );
    addAnalysisObject(hist);
    MSG_TRACE("Made histogram " << hname <<  " for " << name());
    hist->setTitle(title);
    hist->setAnnotation("XLabel", xtitle);
    hist->setAnnotation("YLabel", ytitle);
    return hist;
  }
Histo1DPtr bookHisto1D ( const std::string &  name,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D histogram, using the binnings in the reference data histogram.

Definition at line 246 of file Analysis.cc.

References Analysis::bookHisto1D(), and Analysis::refData().

                                                         {
    const Scatter2D& refdata = refData(hname);
    return bookHisto1D(hname, refdata, title, xtitle, ytitle);
  }
Histo1DPtr bookHisto1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D histogram, using the binnings in the reference data histogram.

The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way.

Definition at line 255 of file Analysis.cc.

References Analysis::bookHisto1D(), and Analysis::makeAxisCode().

                                                         {
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    return bookHisto1D(axisCode, title, xtitle, ytitle);
  }
Histo2DPtr bookHisto2D ( 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.

Todo:
Add booking methods which take a path, titles and *a reference Scatter from which to book*

Definition at line 270 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by MC_PDFS::init().

  {
    const string path = histoPath(hname);
    Histo2DPtr hist( new Histo2D(nxbins, xlower, xupper, nybins, ylower, yupper, path, title) );
    addAnalysisObject(hist);
    MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
    hist->setAnnotation("XLabel", xtitle);
    hist->setAnnotation("YLabel", ytitle);
    hist->setAnnotation("ZLabel", ztitle);
    return hist;
  }
Histo2DPtr bookHisto2D ( 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.

Definition at line 289 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

  {
    const string path = histoPath(hname);
    Histo2DPtr hist( new Histo2D(xbinedges, ybinedges, path, title) );
    addAnalysisObject(hist);
    MSG_TRACE("Made 2D histogram " << hname <<  " for " << name());
    hist->setAnnotation("XLabel", xtitle);
    hist->setAnnotation("YLabel", ytitle);
    hist->setAnnotation("ZLabel", ztitle);
    return hist;
  }
Profile1DPtr 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 .

Definition at line 349 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

Referenced by Analysis::bookProfile1D(), ATLAS_2010_S8894728::init(), ATLAS_2011_S8994773::init(), CMS_2012_PAS_QCD_11_010::init(), ATLAS_2014_I1298811::init(), ALICE_2011_S8945144::init(), ATLAS_2010_S8591806::init(), CDF_2012_NOTE10874::init(), STAR_2008_S7993412::init(), ATLAS_2012_I1124167::init(), CMS_2013_I1261026::init(), CDF_2008_S7782535::init(), ATLAS_2013_I1243871::init(), D0_1996_S3324664::init(), ATLAS_2011_S8924791::init(), CDF_2005_S6217184::init(), CMS_2012_I1107658::init(), MC_GENERIC::init(), MC_LEADJETUE::init(), UA1_1990_S2044935::init(), MC_PHOTONJETUE::init(), STAR_2009_UE_HELEN::init(), STAR_2006_S6860818::init(), ALICE_2010_S8706239::init(), ALEPH_2001_S4656318::init(), DELPHI_2002_069_CONF_603::init(), MC_PHOTONS::init(), MC_WPOL::init(), ATLAS_2014_I1307243::init(), CDF_2009_S8233977::init(), CDF_2010_S8591881_DY::init(), CDF_2010_S8591881_QCD::init(), CDF_2001_S4751469::init(), ATLAS_2010_S8918562::init(), DELPHI_1996_S3430090::init(), CDF_2004_S5839831::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), ATLAS_2014_I1279489::initialisePlots(), ATLAS_2011_S9126244::initializePlots(), and ATLAS_2012_I1125575::initializeProfiles().

                                                             {
    const string path = histoPath(hname);
    Profile1DPtr prof( new Profile1D(nbins, lower, upper, path, title) );
    addAnalysisObject(prof);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setAnnotation("XLabel", xtitle);
    prof->setAnnotation("YLabel", ytitle);
    return prof;
  }
Profile1DPtr 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 .

Definition at line 364 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                             {
    const string path = histoPath(hname);
    Profile1DPtr prof( new Profile1D(binedges, path, title) );
    addAnalysisObject(prof);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setAnnotation("XLabel", xtitle);
    prof->setAnnotation("YLabel", ytitle);
    return prof;
  }
Profile1DPtr bookProfile1D ( const std::string &  name,
const Scatter2D &  refscatter,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D profile histogram with binning from a reference scatter.

Definition at line 379 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

                                                             {
    const string path = histoPath(hname);
    Profile1DPtr prof( new Profile1D(refscatter, path) );
    addAnalysisObject(prof);
    MSG_TRACE("Made profile histogram " << hname <<  " for " << name());
    prof->setTitle(title);
    prof->setAnnotation("XLabel", xtitle);
    prof->setAnnotation("YLabel", ytitle);
    return prof;
  }
Profile1DPtr bookProfile1D ( const std::string &  name,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D profile histogram, using the binnings in the reference data histogram.

Definition at line 395 of file Analysis.cc.

References Analysis::bookProfile1D(), and Analysis::refData().

                                                             {
    const Scatter2D& refdata = refData(hname);
    return bookProfile1D(hname, refdata, title, xtitle, ytitle);
  }
Profile1DPtr bookProfile1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 1D profile histogram, using the binnings in the reference data histogram.

The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way.

Definition at line 404 of file Analysis.cc.

References Analysis::bookProfile1D(), and Analysis::makeAxisCode().

                                                             {
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    return bookProfile1D(axisCode, title, xtitle, ytitle);
  }
Profile2DPtr bookProfile2D ( 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 profile histogram with nxbins and nybins uniformly distributed across the ranges xlower - xupper and ylower - yupper respectively along the x- and y-axis.

Definition at line 417 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

  {
    const string path = histoPath(hname);
    Profile2DPtr prof( new Profile2D(nxbins, xlower, xupper, nybins, ylower, yupper, path, title) );
    addAnalysisObject(prof);
    MSG_TRACE("Made 2D profile histogram " << hname <<  " for " << name());
    prof->setAnnotation("XLabel", xtitle);
    prof->setAnnotation("YLabel", ytitle);
    prof->setAnnotation("ZLabel", ztitle);
    return prof;
  }
Profile2DPtr bookProfile2D ( 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 profile histogram with non-uniform bins defined by the vectorx of bin edges xbinedges and ybinedges.

Definition at line 436 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, and Analysis::name().

  {
    const string path = histoPath(hname);
    Profile2DPtr prof( new Profile2D(xbinedges, ybinedges, path, title) );
    addAnalysisObject(prof);
    MSG_TRACE("Made 2D profile histogram " << hname <<  " for " << name());
    prof->setAnnotation("XLabel", xtitle);
    prof->setAnnotation("YLabel", ytitle);
    prof->setAnnotation("ZLabel", ztitle);
    return prof;
  }
Scatter2DPtr bookScatter2D ( const std::string &  name,
bool  copy_pts = false,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 2-dimensional data point set with the given name.

Note:
Unlike histogram booking, scatter booking by default makes no attempt to use reference data to pre-fill the data object. If you want this, which is sometimes useful e.g. when the x-position is not really meaningful and can't be extracted from the data, then set the copy_pts parameter to true. This creates points to match the reference data's x values and errors, but with the y values and errors zeroed... assuming that there is a reference histo with the same name: if there isn't, an exception will be thrown.

Definition at line 506 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Analysis::refData(), and Rivet::s.

Referenced by Analysis::bookScatter2D(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), DELPHI_2000_S4328825::finalize(), OPAL_2002_S5361494::finalize(), BABAR_2007_S7266081::finalize(), LHCB_2011_I917009::finalize(), ARGUS_1993_S2669951::finalize(), ATLAS_2014_I1307243::finalize(), ALEPH_2004_S5765862::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2011_S9126244::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2014_I1279489::finalizeEfficiencies(), ATLAS_2014_I1312627::hInit(), ATLAS_2010_S8894728::init(), ALICE_2012_I1181770::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), ATLAS_2011_S9002537::init(), CMS_2011_S8978280::init(), CMS_2011_S9088458::init(), STAR_2006_S6500200::init(), CMS_2014_I1298810::init(), ATLAS_2011_S9128077::init(), MC_JetSplittings::init(), CDF_2008_S7782535::init(), MC_WINC::init(), ATLAS_2011_I944826::init(), CMS_2012_I1102908::init(), MC_JetAnalysis::init(), CDF_2005_S6217184::init(), MC_ParticleAnalysis::init(), ATLAS_2011_I945498::init(), CDF_1994_S2952106::init(), ATLAS_2012_I1188891::init(), D0_2001_S4674421::init(), D0_2008_S7837160::init(), MC_GENERIC::init(), STAR_2006_S6860818::init(), MC_XS::init(), CMS_2013_I1218372::init(), LHCB_2012_I1119400::init(), CDF_1996_S3418421::init(), UA5_1988_S1867512::init(), CDF_2008_S7541902::init(), ATLAS_2013_I1230812::init(), ATLAS_2014_I1307243::init(), ATLAS_2012_I1083318::init(), D0_2008_S7719523::init(), ATLAS_2012_I1093734::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), ATLAS_2012_I1094568::initializePlots(), ATLAS_2011_S9126244::initializePlots(), and SLD_1996_S3398250::multiplicity_subtract().

                                                             {
    Scatter2DPtr s;
    const string path = histoPath(hname);
    if (copy_pts) {
      const Scatter2D& refdata = refData(hname);
      s.reset( new Scatter2D(refdata, path) );
      foreach (Point2D& p, s->points()) p.setY(0, 0);
    } else {
      s.reset( new Scatter2D(path) );
    }
    addAnalysisObject(s);
    MSG_TRACE("Made scatter " << hname <<  " for " << name());
    s->setTitle(title);
    s->setAnnotation("XLabel", xtitle);
    s->setAnnotation("YLabel", ytitle);
    return s;
  }
Scatter2DPtr bookScatter2D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId,
bool  copy_pts = false,
const std::string &  title = "",
const std::string &  xtitle = "",
const std::string &  ytitle = "" 
) [protected, inherited]

Book a 2-dimensional data point set, using the binnings in the reference data histogram.

The paper, dataset and x/y-axis IDs will be used to build the histo name in the HepData standard way.

Note:
Unlike histogram booking, scatter booking by default makes no attempt to use reference data to pre-fill the data object. If you want this, which is sometimes useful e.g. when the x-position is not really meaningful and can't be extracted from the data, then set the copy_pts parameter to true. This creates points to match the reference data's x values and errors, but with the y values and errors zeroed.

Definition at line 496 of file Analysis.cc.

References Analysis::bookScatter2D(), and Analysis::makeAxisCode().

                                                             {
    const string axisCode = makeAxisCode(datasetId, xAxisId, yAxisId);
    return bookScatter2D(axisCode, copy_pts, title, xtitle, ytitle);
  }
Scatter2DPtr bookScatter2D ( 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 x-points in a range.

The y values and errors will be set to 0.

Definition at line 529 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Rivet::Cuts::pt, and Rivet::s.

                                                             {
    const string path = histoPath(hname);
    Scatter2DPtr s( new Scatter2D(path) );
    const double binwidth = (upper-lower)/npts;
    for (size_t pt = 0; pt < npts; ++pt) {
      const double bincentre = lower + (pt + 0.5) * binwidth;
      s->addPoint(bincentre, 0, binwidth/2.0, 0);
    }
    addAnalysisObject(s);
    MSG_TRACE("Made scatter " << hname <<  " for " << name());
    s->setTitle(title);
    s->setAnnotation("XLabel", xtitle);
    s->setAnnotation("YLabel", ytitle);
    return s;
  }
Scatter2DPtr bookScatter2D ( const std::string &  hname,
const std::vector< double > &  binedges,
const std::string &  title,
const std::string &  xtitle,
const std::string &  ytitle 
) [protected, inherited]

Book a 2-dimensional data point set based on provided contiguous "bin edges".

The y values and errors will be set to 0.

Definition at line 550 of file Analysis.cc.

References Analysis::addAnalysisObject(), Analysis::histoPath(), MSG_TRACE, Analysis::name(), Rivet::Cuts::pt, and Rivet::s.

                                                             {
    const string path = histoPath(hname);
    Scatter2DPtr s( new Scatter2D(path) );
    for (size_t pt = 0; pt < binedges.size()-1; ++pt) {
      const double bincentre = (binedges[pt] + binedges[pt+1]) / 2.0;
      const double binwidth = binedges[pt+1] - binedges[pt];
      s->addPoint(bincentre, 0, binwidth/2.0, 0);
    }
    addAnalysisObject(s);
    MSG_TRACE("Made scatter " << hname <<  " for " << name());
    s->setTitle(title);
    s->setAnnotation("XLabel", xtitle);
    s->setAnnotation("YLabel", ytitle);
    return s;
  }
virtual std::string collider ( ) const [inline, virtual, inherited]

Collider on which the experiment ran.

Definition at line 170 of file Analysis.hh.

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

                                       {
      return info().collider();
    }
double crossSection ( ) const [protected, inherited]

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

Definition at line 155 of file Analysis.cc.

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

Referenced by MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), D0_2000_S4480767::finalize(), TOTEM_2012_002::finalize(), CMS_2011_S9086218::finalize(), CDF_2001_S4563131::finalize(), CDF_2008_S7828950::finalize(), D0_2010_S8671338::finalize(), MC_XS::finalize(), CDF_2008_S8093652::finalize(), BELLE_2013_I1216515::finalize(), CDF_2000_S4155203::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), CMS_2012_I1184941::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1087342::finalize(), LHCB_2012_I1208102::finalize(), MC_HINC::finalize(), MC_ZJETS::finalize(), CMS_2011_I954992::finalize(), MC_HJETS::finalize(), MC_WJETS::finalize(), ATLAS_2012_I1118269::finalize(), ATLAS_2011_I894867::finalize(), BABAR_2005_S6181155::finalize(), CDF_2009_S8383952::finalize(), CDF_2001_S4517016::finalize(), D0_2010_S8566488::finalize(), STAR_2006_S6870392::finalize(), MC_ZINC::finalize(), LHCB_2010_I867355::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), D0_2011_I895662::finalize(), CDF_2005_S6080774::finalize(), ATLAS_2014_I1288706::finalize(), ATLAS_2014_I1325553::finalize(), ATLAS_2014_I1282441::finalize(), D0_2006_S6438750::finalize(), CMS_2013_I1273574::finalize(), MC_DIPHOTON::finalize(), D0_2000_I499943::finalize(), ALICE_2012_I1181770::finalize(), CLEO_2004_S5809304::finalize(), D0_2008_S7662670::finalize(), D0_2008_S7863608::finalize(), MC_ZZJETS::finalize(), MC_WINC::finalize(), ATLAS_2011_S9035664::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), LHCB_2013_I1218996::finalize(), ATLAS_2013_I1219109::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), CDF_2008_S7541902::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_ZZINC::finalize(), MC_WWINC::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ATLAS_2013_I1217863_Z::finalize(), D0_2008_S7719523::finalize(), ATLAS_2013_I1217863_W::finalize(), ATLAS_2014_I1307756::finalize(), ATLAS_2011_I926145::finalize(), CMS_2013_I1256943::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2014_I1307243::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_CONF_2012_105::finalize(), CMS_2014_I1303894::finalize(), ATLAS_2012_CONF_2012_104::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2012_I1117704::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ATLAS_2012_CONF_2012_153::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), and ATLAS_2011_I919017::finalize().

                                      {
    if (!_gotCrossSection || std::isnan(_crossSection)) {
      string errMsg = "You did not set the cross section for the analysis " + name();
      throw Error(errMsg);
    }
    return _crossSection;
  }
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 151 of file Analysis.hh.

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

                                          {
      return info().description();
    }
void divide ( const YODA::Histo1D &  h1,
const YODA::Histo1D &  h2,
Scatter2DPtr  s 
) const [inherited]

Helper for histogram division with raw YODA objects.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( Profile1DPtr  p1,
Profile1DPtr  p2,
Scatter2DPtr  s 
) const [inherited]

Helper for profile histogram division.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( const YODA::Profile1D &  p1,
const YODA::Profile1D &  p2,
Scatter2DPtr  s 
) const [inherited]

Helper for profile histogram division with raw YODA objects.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( Histo2DPtr  h1,
Histo2DPtr  h2,
Scatter3DPtr  s 
) const [inherited]

Helper for 2D histogram division.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( const YODA::Histo2D &  h1,
const YODA::Histo2D &  h2,
Scatter3DPtr  s 
) const [inherited]

Helper for 2D histogram division with raw YODA objects.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( Profile2DPtr  p1,
Profile2DPtr  p2,
Scatter3DPtr  s 
) const [inherited]

Helper for 2D profile histogram division.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void divide ( const YODA::Profile2D &  p1,
const YODA::Profile2D &  p2,
Scatter3DPtr  s 
) const [inherited]

Helper for 2D profile histogram division with raw YODA objects

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
void efficiency ( Histo1DPtr  h1,
Histo1DPtr  h2,
Scatter2DPtr  s 
) const [inherited]

Helper for histogram efficiency calculation.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.

Definition at line 626 of file Analysis.cc.

Referenced by CMS_2012_I1102908::finalize(), and ATLAS_2014_I1307243::finalize().

                                                                              {
    const string path = s->path();
    *s = YODA::efficiency(*h1, *h2);
    s->setPath(path);
  }
void efficiency ( const YODA::Histo1D &  h1,
const YODA::Histo1D &  h2,
Scatter2DPtr  s 
) const [inherited]

Helper for histogram efficiency calculation.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
virtual std::string experiment ( ) const [inline, virtual, inherited]

Experiment which performed and published this analysis.

Definition at line 165 of file Analysis.hh.

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

                                         {
      return info().experiment();
    }
void fillEventCountsPerSR ( const string &  basic_signal_region,
int  onZ,
double  HTlep,
double  eTmiss,
double  HTjets,
double  meff,
double  weight 
) [inline]

function fills map EventCountsPerSR by looping over all signal regions and looking if the event falls into this signal region

Definition at line 579 of file ATLAS_2012_I1204447.cc.

References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::getCutsPerSignalRegion(), and Rivet::toString().

Referenced by ATLAS_2012_I1204447::analyze().

                                              {

      // Get cut values for HTlep, loop over them and add event if cut is passed
      vector<int> cut_values = getCutsPerSignalRegion("HTlep", onZ);
      for (size_t i = 0; i < cut_values.size(); i++)  {
        if (HTlep > cut_values[i])
          _eventCountsPerSR[("HTlep_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight;
      }

      // Get cut values for METStrong, loop over them and add event if cut is passed
      cut_values = getCutsPerSignalRegion("METStrong", onZ);
      for (size_t i = 0; i < cut_values.size(); i++)  {
        if (eTmiss > cut_values[i] && HTjets > 100.)
          _eventCountsPerSR[("METStrong_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight;
      }

      // Get cut values for METWeak, loop over them and add event if cut is passed
      cut_values = getCutsPerSignalRegion("METWeak", onZ);
      for (size_t i = 0; i < cut_values.size(); i++)  {
        if (eTmiss > cut_values[i] && HTjets <= 100.)
          _eventCountsPerSR[("METWeak_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight;
      }

      // Get cut values for Meff, loop over them and add event if cut is passed
      cut_values = getCutsPerSignalRegion("Meff", onZ);
      for (size_t i = 0; i < cut_values.size(); i++)  {
        if (meff > cut_values[i])
          _eventCountsPerSR[("Meff_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight;
      }

      // Get cut values for MeffStrong, loop over them and add event if cut is passed
      cut_values = getCutsPerSignalRegion("MeffStrong", onZ);
      for (size_t i = 0; i < cut_values.size(); i++)  {
        if (meff > cut_values[i] && eTmiss > 75.)
          _eventCountsPerSR[("MeffStrong_" + basic_signal_region + "_cut_" + toString(cut_values[i]))] += weight;
      }
    }
void finalize ( ) [inline, virtual]

Normalise histograms etc., after the run.

Reimplemented from Analysis.

Definition at line 425 of file ATLAS_2012_I1204447.cc.

References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_excluded, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_signal_regions, Analysis::crossSection(), Rivet::femtobarn, ATLAS_2012_I1204447::getUpperLimit(), Analysis::scale(), and Analysis::sumOfWeights().

                    {

      // Normalize to an integrated luminosity of 1 fb-1
      double norm = crossSection()/femtobarn/sumOfWeights();
      string best_signal_region = "";
      double ratio_best_SR = 0.;

      // Loop over all signal regions and find signal region with best sensitivity (ratio signal events/visible cross-section)
      for (size_t i = 0; i < _signal_regions.size(); i++)  {
        double signal_events = _eventCountsPerSR[_signal_regions[i]] * norm;
        // Use expected upper limits to find best signal region
        double UL95  = getUpperLimit(_signal_regions[i], false);
        double ratio = signal_events / UL95;
        if (ratio > ratio_best_SR)  {
          best_signal_region = _signal_regions[i];
          ratio_best_SR = ratio;
        }
      }

      double signal_events_best_SR = _eventCountsPerSR[best_signal_region] * norm;
      double exp_UL_best_SR = getUpperLimit(best_signal_region, false);
      double obs_UL_best_SR = getUpperLimit(best_signal_region, true);

      // Print out result
      cout << "----------------------------------------------------------------------------------------" << endl;
      cout << "Best signal region: " << best_signal_region << endl;
      cout << "Normalized number of signal events in this best signal region (per fb-1): " << signal_events_best_SR << endl;
      cout << "Efficiency*Acceptance: " <<  _eventCountsPerSR[best_signal_region]/sumOfWeights() << endl;
      cout << "Cross-section [fb]: " << crossSection()/femtobarn << endl;
      cout << "Expected visible cross-section (per fb-1): " << exp_UL_best_SR << endl;
      cout << "Ratio (signal events / expected visible cross-section): " << ratio_best_SR << endl;
      cout << "Observed visible cross-section (per fb-1): " << obs_UL_best_SR << endl;
      cout << "Ratio (signal events / observed visible cross-section): " <<  signal_events_best_SR/obs_UL_best_SR << endl;
      cout << "----------------------------------------------------------------------------------------" << endl;

      cout << "Using the EXPECTED limits (visible cross-section) of the analysis: " << endl;
      if (signal_events_best_SR > exp_UL_best_SR)  {
        cout << "Since the number of signal events > the visible cross-section, this model/grid point is EXCLUDED with 95% CL." << endl;
        _h_excluded->fill(1);
      }
      else  {
        cout << "Since the number of signal events < the visible cross-section, this model/grid point is NOT EXCLUDED." << endl;
        _h_excluded->fill(0);
      }
      cout << "----------------------------------------------------------------------------------------" << endl;

      cout << "Using the OBSERVED limits (visible cross-section) of the analysis: " << endl;
      if (signal_events_best_SR > obs_UL_best_SR)  {
        cout << "Since the number of signal events > the visible cross-section, this model/grid point is EXCLUDED with 95% CL." << endl;
        _h_excluded->fill(1);
      }
      else  {
        cout << "Since the number of signal events < the visible cross-section, this model/grid point is NOT EXCLUDED." << endl;
        _h_excluded->fill(0);
      }
      cout << "----------------------------------------------------------------------------------------" << endl;


      // Normalize to cross section
      if (norm != 0)  {
        scale(_h_HTlep_all,  norm);
        scale(_h_HTjets_all, norm);
        scale(_h_MET_all,    norm);
        scale(_h_Meff_all,   norm);

        scale(_h_pt_1_3l,    norm);
        scale(_h_pt_2_3l,    norm);
        scale(_h_pt_3_3l,    norm);
        scale(_h_pt_1_2ltau, norm);
        scale(_h_pt_2_2ltau, norm);
        scale(_h_pt_3_2ltau, norm);

        scale(_h_e_n,        norm);
        scale(_h_mu_n,       norm);
        scale(_h_tau_n,      norm);

        scale(_h_excluded,   signal_events_best_SR);
      }
    }
void get_prong_number ( const GenParticle *  p,
unsigned int &  nprong,
bool &  lep_decaying_tau 
) [inline]

Function calculating the prong number of taus

Todo:
Move to TauFinder and make less HepMC-ish

Definition at line 636 of file ATLAS_2012_I1204447.cc.

References Rivet::PID::charge(), Rivet::PID::ELECTRON, Rivet::PID::MUON, and Rivet::PID::TAU.

Referenced by ATLAS_2012_I1204447::analyze().

                                                                                              {
      assert(p != NULL);
      //const int tau_barcode = p->barcode();
      const GenVertex* dv = p->end_vertex();
      assert(dv != NULL);
      for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin(); pp != dv->particles_out_const_end(); ++pp) {
        // If they have status 1 and are charged they will produce a track and the prong number is +1
        if ((*pp)->status() == 1 )  {
          const int id = (*pp)->pdg_id();
          if (Rivet::PID::charge(id) != 0 ) ++nprong;
          // Check if tau decays leptonically
          // @todo Can a tau decay include a tau in its decay daughters?!
          if ((abs(id) == PID::ELECTRON || abs(id) == PID::MUON || abs(id) == PID::TAU) && abs(p->pdg_id()) == PID::TAU) lep_decaying_tau = true;
        }
        // If the status of the daughter particle is 2 it is unstable and the further decays are checked
        else if ((*pp)->status() == 2 )  {
          get_prong_number(*pp, nprong, lep_decaying_tau);
        }
      }
    }
FourMomentum get_tau_neutrino_mom ( const Particle p) [inline]

Function returning 4-vector of daughter-particle if it is a tau neutrino

Todo:
Move to TauFinder and make less HepMC-ish

Definition at line 623 of file ATLAS_2012_I1204447.cc.

References Particle::abspid(), Particle::genParticle(), Rivet::PID::NU_TAU, and Rivet::PID::TAU.

Referenced by ATLAS_2012_I1204447::analyze().

                                                          {
      assert(p.abspid() == PID::TAU);
      const GenVertex* dv = p.genParticle()->end_vertex();
      assert(dv != NULL);
      for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin(); pp != dv->particles_out_const_end(); ++pp) {
        if ((*pp)->pdg_id() == PID::NU_TAU) return FourMomentum((*pp)->momentum());
      }
      return FourMomentum();
    }
const shared_ptr<AO> getAnalysisObject ( const std::string &  name) const [inline, protected, inherited]

Get a data object from the histogram system

Todo:
Use this default function template arg in C++11

Definition at line 711 of file Analysis.hh.

References Analysis::analysisObjects(), and Analysis::histoPath().

                                                                        {
      foreach (const AnalysisObjectPtr& ao, analysisObjects()) {
        if (ao->path() == histoPath(name)) return dynamic_pointer_cast<AO>(ao);
      }
      throw Exception("Data object " + histoPath(name) + " not found");
    }
shared_ptr<AO> getAnalysisObject ( const std::string &  name) [inline, protected, inherited]

Get a data object from the histogram system (non-const)

Todo:
Use this default function template arg in C++11

Definition at line 722 of file Analysis.hh.

References Analysis::analysisObjects(), and Analysis::histoPath().

                                                            {
      foreach (const AnalysisObjectPtr& ao, analysisObjects()) {
        if (ao->path() == histoPath(name)) return dynamic_pointer_cast<AO>(ao);
      }
      throw Exception("Data object " + histoPath(name) + " not found");
    }
vector<int> getCutsPerSignalRegion ( const string &  signal_region,
int  onZ = 0 
) [inline]

Function giving all cut vales per kinematic variable (taking onZ for MET into account)

Definition at line 547 of file ATLAS_2012_I1204447.cc.

Referenced by ATLAS_2012_I1204447::fillEventCountsPerSR(), and ATLAS_2012_I1204447::getSignalRegions().

                                                                                {
      vector<int> cutValues;

      // Cut values for HTlep
      if (signal_region.compare("HTlep") == 0)  {
        cutValues.push_back(0);
        cutValues.push_back(100);
        cutValues.push_back(150);
        cutValues.push_back(200);
        cutValues.push_back(300);
      }
      // Cut values for METStrong (HTjets > 100 GeV) and METWeak (HTjets < 100 GeV)
      else if (signal_region.compare("METStrong") == 0 || signal_region.compare("METWeak") == 0)  {
        if      (onZ == 0) cutValues.push_back(0);
        else if (onZ == 1) cutValues.push_back(20);
        cutValues.push_back(50);
        cutValues.push_back(75);
      }
      // Cut values for Meff and MeffStrong (MET > 75 GeV)
      if (signal_region.compare("Meff") == 0 || signal_region.compare("MeffStrong") == 0)  {
        cutValues.push_back(0);
        cutValues.push_back(150);
        cutValues.push_back(300);
        cutValues.push_back(500);
      }

      return cutValues;
    }
const Histo1DPtr getHisto1D ( const std::string &  name) const [inline, protected, inherited]

Get a named Histo1D object from the histogram system.

Definition at line 737 of file Analysis.hh.

References Analysis::name().

                                                             {
      return getAnalysisObject<Histo1D>(name);
    }
Histo1DPtr getHisto1D ( const std::string &  name) [inline, protected, inherited]

Get a named Histo1D object from the histogram system (non-const)

Definition at line 742 of file Analysis.hh.

References Analysis::name().

                                                 {
      return getAnalysisObject<Histo1D>(name);
    }
const Histo1DPtr getHisto1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [inline, protected, inherited]

Get a Histo1D object from the histogram system by axis ID codes (non-const)

Definition at line 747 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                                          {
      return getAnalysisObject<Histo1D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
Histo1DPtr getHisto1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) [inline, protected, inherited]

Get a Histo1D object from the histogram system by axis ID codes (non-const)

Definition at line 752 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                              {
      return getAnalysisObject<Histo1D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
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 82 of file Analysis.cc.

References Analysis::name().

                              {
    string logname = "Rivet.Analysis." + name();
    return Log::getLog(logname);
  }
const Profile1DPtr getProfile1D ( const std::string &  name) const [inline, protected, inherited]

Get a named Profile1D object from the histogram system.

Definition at line 779 of file Analysis.hh.

References Analysis::name().

                                                                 {
      return getAnalysisObject<Profile1D>(name);
    }
Profile1DPtr getProfile1D ( const std::string &  name) [inline, protected, inherited]

Get a named Profile1D object from the histogram system (non-const)

Definition at line 784 of file Analysis.hh.

References Analysis::name().

                                                     {
      return getAnalysisObject<Profile1D>(name);
    }
const Profile1DPtr getProfile1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [inline, protected, inherited]

Get a Profile1D object from the histogram system by axis ID codes (non-const)

Definition at line 789 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                                              {
      return getAnalysisObject<Profile1D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
Profile1DPtr getProfile1D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) [inline, protected, inherited]

Get a Profile1D object from the histogram system by axis ID codes (non-const)

Definition at line 794 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                                  {
      return getAnalysisObject<Profile1D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
const PROJ& getProjection ( const std::string &  name) const [inline, inherited]

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

Definition at line 52 of file ProjectionApplier.hh.

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

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

                                                           {
      const Projection& p = getProjHandler().getProjection(*this, name);
      return pcast<PROJ>(p);
    }
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 60 of file ProjectionApplier.hh.

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

                                                                 {
      return getProjHandler().getProjection(*this, name);
    }
std::set<ConstProjectionPtr> getProjections ( ) const [inline, inherited]

Get the contained projections, including recursion.

Definition at line 45 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

const Scatter2DPtr getScatter2D ( const std::string &  name) const [inline, protected, inherited]

Get a named Scatter2D object from the histogram system.

Definition at line 821 of file Analysis.hh.

References Analysis::name().

Referenced by ATLAS_2014_I1279489::finalizeEfficiencies().

                                                                 {
      return getAnalysisObject<Scatter2D>(name);
    }
Scatter2DPtr getScatter2D ( const std::string &  name) [inline, protected, inherited]

Get a named Scatter2D object from the histogram system (non-const)

Definition at line 826 of file Analysis.hh.

References Analysis::name().

                                                     {
      return getAnalysisObject<Scatter2D>(name);
    }
const Scatter2DPtr getScatter2D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [inline, protected, inherited]

Get a Scatter2D object from the histogram system by axis ID codes (non-const)

Definition at line 831 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                                              {
      return getAnalysisObject<Scatter2D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
Scatter2DPtr getScatter2D ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) [inline, protected, inherited]

Get a Scatter2D object from the histogram system by axis ID codes (non-const)

Definition at line 836 of file Analysis.hh.

References Analysis::makeAxisCode().

                                                                                                  {
      return getAnalysisObject<Scatter2D>(makeAxisCode(datasetId, xAxisId, yAxisId));
    }
vector<string> getSignalRegions ( ) [inline]

Helper functions.

Function giving a list of all signal regions

Definition at line 510 of file ATLAS_2012_I1204447.cc.

References ATLAS_2012_I1204447::getCutsPerSignalRegion(), and Rivet::toString().

Referenced by ATLAS_2012_I1204447::init().

                                       {

      // List of basic signal regions
      vector<string> basic_signal_regions;
      basic_signal_regions.push_back("3l_offZ");
      basic_signal_regions.push_back("3l_onZ");
      basic_signal_regions.push_back("2ltau_offZ");
      basic_signal_regions.push_back("2ltau_onZ");

      // List of kinematic variables
      vector<string> kinematic_variables;
      kinematic_variables.push_back("HTlep");
      kinematic_variables.push_back("METStrong");
      kinematic_variables.push_back("METWeak");
      kinematic_variables.push_back("Meff");
      kinematic_variables.push_back("MeffStrong");

      vector<string> signal_regions;
      // Loop over all kinematic variables and basic signal regions
      for (size_t i0 = 0; i0 < kinematic_variables.size(); i0++)  {
        for (size_t i1 = 0; i1 < basic_signal_regions.size(); i1++)  {
          // Is signal region onZ?
          int onZ = (basic_signal_regions[i1].find("onZ") != string::npos) ? 1 : 0;
          // Get cut values for this kinematic variable
          vector<int> cut_values = getCutsPerSignalRegion(kinematic_variables[i0], onZ);
          // Loop over all cut values
          for (size_t i2 = 0; i2 < cut_values.size(); i2++)  {
            // push signal region into vector
            signal_regions.push_back( (kinematic_variables[i0] + "_" + basic_signal_regions[i1] + "_cut_" + toString(i2)) );
          }
        }
      }
      return signal_regions;
    }
double getUpperLimit ( const string &  signal_region,
bool  observed 
) [inline]

Function giving observed upper limit (visible cross-section)

Definition at line 831 of file ATLAS_2012_I1204447.cc.

Referenced by ATLAS_2012_I1204447::finalize().

                                                                      {
      map<string,double> upperLimitsObserved;
      upperLimitsObserved["HTlep_3l_offZ_cut_0"] = 11.;
      upperLimitsObserved["HTlep_3l_offZ_cut_100"] = 8.7;
      upperLimitsObserved["HTlep_3l_offZ_cut_150"] = 4.0;
      upperLimitsObserved["HTlep_3l_offZ_cut_200"] = 4.4;
      upperLimitsObserved["HTlep_3l_offZ_cut_300"] = 1.6;
      upperLimitsObserved["HTlep_2ltau_offZ_cut_0"] = 25.;
      upperLimitsObserved["HTlep_2ltau_offZ_cut_100"] = 14.;
      upperLimitsObserved["HTlep_2ltau_offZ_cut_150"] = 6.1;
      upperLimitsObserved["HTlep_2ltau_offZ_cut_200"] = 3.3;
      upperLimitsObserved["HTlep_2ltau_offZ_cut_300"] = 1.2;
      upperLimitsObserved["HTlep_3l_onZ_cut_0"] = 48.;
      upperLimitsObserved["HTlep_3l_onZ_cut_100"] = 38.;
      upperLimitsObserved["HTlep_3l_onZ_cut_150"] = 14.;
      upperLimitsObserved["HTlep_3l_onZ_cut_200"] = 7.2;
      upperLimitsObserved["HTlep_3l_onZ_cut_300"] = 4.5;
      upperLimitsObserved["HTlep_2ltau_onZ_cut_0"] = 85.;
      upperLimitsObserved["HTlep_2ltau_onZ_cut_100"] = 53.;
      upperLimitsObserved["HTlep_2ltau_onZ_cut_150"] = 11.0;
      upperLimitsObserved["HTlep_2ltau_onZ_cut_200"] = 5.2;
      upperLimitsObserved["HTlep_2ltau_onZ_cut_300"] = 3.0;
      upperLimitsObserved["METStrong_3l_offZ_cut_0"] = 2.6;
      upperLimitsObserved["METStrong_3l_offZ_cut_50"] = 2.1;
      upperLimitsObserved["METStrong_3l_offZ_cut_75"] = 2.1;
      upperLimitsObserved["METStrong_2ltau_offZ_cut_0"] = 4.2;
      upperLimitsObserved["METStrong_2ltau_offZ_cut_50"] = 3.1;
      upperLimitsObserved["METStrong_2ltau_offZ_cut_75"] = 2.6;
      upperLimitsObserved["METStrong_3l_onZ_cut_20"] = 11.0;
      upperLimitsObserved["METStrong_3l_onZ_cut_50"] = 6.4;
      upperLimitsObserved["METStrong_3l_onZ_cut_75"] = 5.1;
      upperLimitsObserved["METStrong_2ltau_onZ_cut_20"] = 5.9;
      upperLimitsObserved["METStrong_2ltau_onZ_cut_50"] = 3.4;
      upperLimitsObserved["METStrong_2ltau_onZ_cut_75"] = 1.2;
      upperLimitsObserved["METWeak_3l_offZ_cut_0"] = 11.;
      upperLimitsObserved["METWeak_3l_offZ_cut_50"] = 5.3;
      upperLimitsObserved["METWeak_3l_offZ_cut_75"] = 3.1;
      upperLimitsObserved["METWeak_2ltau_offZ_cut_0"] = 23.;
      upperLimitsObserved["METWeak_2ltau_offZ_cut_50"] = 4.3;
      upperLimitsObserved["METWeak_2ltau_offZ_cut_75"] = 3.1;
      upperLimitsObserved["METWeak_3l_onZ_cut_20"] = 41.;
      upperLimitsObserved["METWeak_3l_onZ_cut_50"] = 16.;
      upperLimitsObserved["METWeak_3l_onZ_cut_75"] = 8.0;
      upperLimitsObserved["METWeak_2ltau_onZ_cut_20"] = 80.;
      upperLimitsObserved["METWeak_2ltau_onZ_cut_50"] = 4.4;
      upperLimitsObserved["METWeak_2ltau_onZ_cut_75"] = 1.8;
      upperLimitsObserved["Meff_3l_offZ_cut_0"] = 11.;
      upperLimitsObserved["Meff_3l_offZ_cut_150"] = 8.1;
      upperLimitsObserved["Meff_3l_offZ_cut_300"] = 3.1;
      upperLimitsObserved["Meff_3l_offZ_cut_500"] = 2.1;
      upperLimitsObserved["Meff_2ltau_offZ_cut_0"] = 25.;
      upperLimitsObserved["Meff_2ltau_offZ_cut_150"] = 12.;
      upperLimitsObserved["Meff_2ltau_offZ_cut_300"] = 3.9;
      upperLimitsObserved["Meff_2ltau_offZ_cut_500"] = 2.2;
      upperLimitsObserved["Meff_3l_onZ_cut_0"] = 48.;
      upperLimitsObserved["Meff_3l_onZ_cut_150"] = 37.;
      upperLimitsObserved["Meff_3l_onZ_cut_300"] = 11.;
      upperLimitsObserved["Meff_3l_onZ_cut_500"] = 4.8;
      upperLimitsObserved["Meff_2ltau_onZ_cut_0"] = 85.;
      upperLimitsObserved["Meff_2ltau_onZ_cut_150"] = 28.;
      upperLimitsObserved["Meff_2ltau_onZ_cut_300"] = 5.9;
      upperLimitsObserved["Meff_2ltau_onZ_cut_500"] = 1.9;
      upperLimitsObserved["MeffStrong_3l_offZ_cut_0"] = 3.8;
      upperLimitsObserved["MeffStrong_3l_offZ_cut_150"] = 3.8;
      upperLimitsObserved["MeffStrong_3l_offZ_cut_300"] = 2.8;
      upperLimitsObserved["MeffStrong_3l_offZ_cut_500"] = 2.1;
      upperLimitsObserved["MeffStrong_2ltau_offZ_cut_0"] = 3.9;
      upperLimitsObserved["MeffStrong_2ltau_offZ_cut_150"] = 4.0;
      upperLimitsObserved["MeffStrong_2ltau_offZ_cut_300"] = 2.9;
      upperLimitsObserved["MeffStrong_2ltau_offZ_cut_500"] = 1.5;
      upperLimitsObserved["MeffStrong_3l_onZ_cut_0"] = 10.0;
      upperLimitsObserved["MeffStrong_3l_onZ_cut_150"] = 10.0;
      upperLimitsObserved["MeffStrong_3l_onZ_cut_300"] = 6.8;
      upperLimitsObserved["MeffStrong_3l_onZ_cut_500"] = 3.9;
      upperLimitsObserved["MeffStrong_2ltau_onZ_cut_0"] = 1.6;
      upperLimitsObserved["MeffStrong_2ltau_onZ_cut_150"] = 1.4;
      upperLimitsObserved["MeffStrong_2ltau_onZ_cut_300"] = 1.5;
      upperLimitsObserved["MeffStrong_2ltau_onZ_cut_500"] = 0.9;

      // Expected upper limits are also given but not used in this analysis
      map<string,double> upperLimitsExpected;
      upperLimitsExpected["HTlep_3l_offZ_cut_0"] = 11.;
      upperLimitsExpected["HTlep_3l_offZ_cut_100"] = 8.5;
      upperLimitsExpected["HTlep_3l_offZ_cut_150"] = 4.6;
      upperLimitsExpected["HTlep_3l_offZ_cut_200"] = 3.6;
      upperLimitsExpected["HTlep_3l_offZ_cut_300"] = 1.9;
      upperLimitsExpected["HTlep_2ltau_offZ_cut_0"] = 23.;
      upperLimitsExpected["HTlep_2ltau_offZ_cut_100"] = 14.;
      upperLimitsExpected["HTlep_2ltau_offZ_cut_150"] = 6.4;
      upperLimitsExpected["HTlep_2ltau_offZ_cut_200"] = 3.6;
      upperLimitsExpected["HTlep_2ltau_offZ_cut_300"] = 1.5;
      upperLimitsExpected["HTlep_3l_onZ_cut_0"] = 33.;
      upperLimitsExpected["HTlep_3l_onZ_cut_100"] = 25.;
      upperLimitsExpected["HTlep_3l_onZ_cut_150"] = 12.;
      upperLimitsExpected["HTlep_3l_onZ_cut_200"] = 6.5;
      upperLimitsExpected["HTlep_3l_onZ_cut_300"] = 3.1;
      upperLimitsExpected["HTlep_2ltau_onZ_cut_0"] = 94.;
      upperLimitsExpected["HTlep_2ltau_onZ_cut_100"] = 61.;
      upperLimitsExpected["HTlep_2ltau_onZ_cut_150"] = 9.9;
      upperLimitsExpected["HTlep_2ltau_onZ_cut_200"] = 4.5;
      upperLimitsExpected["HTlep_2ltau_onZ_cut_300"] = 1.9;
      upperLimitsExpected["METStrong_3l_offZ_cut_0"] = 3.1;
      upperLimitsExpected["METStrong_3l_offZ_cut_50"] = 2.4;
      upperLimitsExpected["METStrong_3l_offZ_cut_75"] = 2.3;
      upperLimitsExpected["METStrong_2ltau_offZ_cut_0"] = 4.8;
      upperLimitsExpected["METStrong_2ltau_offZ_cut_50"] = 3.3;
      upperLimitsExpected["METStrong_2ltau_offZ_cut_75"] = 2.1;
      upperLimitsExpected["METStrong_3l_onZ_cut_20"] = 8.7;
      upperLimitsExpected["METStrong_3l_onZ_cut_50"] = 4.9;
      upperLimitsExpected["METStrong_3l_onZ_cut_75"] = 3.8;
      upperLimitsExpected["METStrong_2ltau_onZ_cut_20"] = 7.3;
      upperLimitsExpected["METStrong_2ltau_onZ_cut_50"] = 2.8;
      upperLimitsExpected["METStrong_2ltau_onZ_cut_75"] = 1.5;
      upperLimitsExpected["METWeak_3l_offZ_cut_0"] = 10.;
      upperLimitsExpected["METWeak_3l_offZ_cut_50"] = 4.7;
      upperLimitsExpected["METWeak_3l_offZ_cut_75"] = 3.0;
      upperLimitsExpected["METWeak_2ltau_offZ_cut_0"] = 21.;
      upperLimitsExpected["METWeak_2ltau_offZ_cut_50"] = 4.0;
      upperLimitsExpected["METWeak_2ltau_offZ_cut_75"] = 2.6;
      upperLimitsExpected["METWeak_3l_onZ_cut_20"] = 30.;
      upperLimitsExpected["METWeak_3l_onZ_cut_50"] = 10.;
      upperLimitsExpected["METWeak_3l_onZ_cut_75"] = 5.4;
      upperLimitsExpected["METWeak_2ltau_onZ_cut_20"] = 88.;
      upperLimitsExpected["METWeak_2ltau_onZ_cut_50"] = 5.5;
      upperLimitsExpected["METWeak_2ltau_onZ_cut_75"] = 2.2;
      upperLimitsExpected["Meff_3l_offZ_cut_0"] = 11.;
      upperLimitsExpected["Meff_3l_offZ_cut_150"] = 8.8;
      upperLimitsExpected["Meff_3l_offZ_cut_300"] = 3.7;
      upperLimitsExpected["Meff_3l_offZ_cut_500"] = 2.1;
      upperLimitsExpected["Meff_2ltau_offZ_cut_0"] = 23.;
      upperLimitsExpected["Meff_2ltau_offZ_cut_150"] = 13.;
      upperLimitsExpected["Meff_2ltau_offZ_cut_300"] = 4.9;
      upperLimitsExpected["Meff_2ltau_offZ_cut_500"] = 2.4;
      upperLimitsExpected["Meff_3l_onZ_cut_0"] = 33.;
      upperLimitsExpected["Meff_3l_onZ_cut_150"] = 25.;
      upperLimitsExpected["Meff_3l_onZ_cut_300"] = 9.;
      upperLimitsExpected["Meff_3l_onZ_cut_500"] = 3.9;
      upperLimitsExpected["Meff_2ltau_onZ_cut_0"] = 94.;
      upperLimitsExpected["Meff_2ltau_onZ_cut_150"] = 35.;
      upperLimitsExpected["Meff_2ltau_onZ_cut_300"] = 6.8;
      upperLimitsExpected["Meff_2ltau_onZ_cut_500"] = 2.5;
      upperLimitsExpected["MeffStrong_3l_offZ_cut_0"] = 3.9;
      upperLimitsExpected["MeffStrong_3l_offZ_cut_150"] = 3.9;
      upperLimitsExpected["MeffStrong_3l_offZ_cut_300"] = 3.0;
      upperLimitsExpected["MeffStrong_3l_offZ_cut_500"] = 2.0;
      upperLimitsExpected["MeffStrong_2ltau_offZ_cut_0"] = 3.8;
      upperLimitsExpected["MeffStrong_2ltau_offZ_cut_150"] = 3.9;
      upperLimitsExpected["MeffStrong_2ltau_offZ_cut_300"] = 3.1;
      upperLimitsExpected["MeffStrong_2ltau_offZ_cut_500"] = 1.6;
      upperLimitsExpected["MeffStrong_3l_onZ_cut_0"] = 6.9;
      upperLimitsExpected["MeffStrong_3l_onZ_cut_150"] = 7.1;
      upperLimitsExpected["MeffStrong_3l_onZ_cut_300"] = 4.9;
      upperLimitsExpected["MeffStrong_3l_onZ_cut_500"] = 3.0;
      upperLimitsExpected["MeffStrong_2ltau_onZ_cut_0"] = 2.4;
      upperLimitsExpected["MeffStrong_2ltau_onZ_cut_150"] = 2.5;
      upperLimitsExpected["MeffStrong_2ltau_onZ_cut_300"] = 2.0;
      upperLimitsExpected["MeffStrong_2ltau_onZ_cut_500"] = 1.1;

      if (observed) return upperLimitsObserved[signal_region];
      else          return upperLimitsExpected[signal_region];
    }
AnalysisHandler& handler ( ) const [inline, inherited]
const string histoDir ( ) const [protected, inherited]

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

Todo:
Cache in a member variable

Definition at line 40 of file Analysis.cc.

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

Referenced by Analysis::histoPath().

                                        {
    /// @todo Cache in a member variable
    string _histoDir;
    if (_histoDir.empty()) {
      _histoDir = "/" + name();
      if (handler().runName().length() > 0) {
        _histoDir = "/" + handler().runName() + _histoDir;
      }
      while (find_first(_histoDir, "//")) {
        replace_all(_histoDir, "//", "/");
      }
    }
    return _histoDir;
  }
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 56 of file Analysis.cc.

References Analysis::histoDir().

Referenced by Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), ATLAS_2012_I1094061::HistoPair::finalize(), H1_1995_S3167097::finalize(), and Analysis::getAnalysisObject().

                                                            {
    const string path = histoDir() + "/" + hname;
    return path;
  }
const string histoPath ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [protected, inherited]

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

Definition at line 62 of file Analysis.cc.

References Analysis::histoDir(), and Analysis::makeAxisCode().

                                                                                                           {
    return histoDir() + "/" + makeAxisCode(datasetId, xAxisId, yAxisId);
  }
AnalysisInfo& info ( ) [inline, inherited]

Get the actual AnalysisInfo object in which all this metadata is stored (non-const).

Definition at line 246 of file Analysis.hh.

References Analysis::_info.

                         {
      assert(_info.get() != 0 && "No AnalysisInfo object :O");
      return *_info;
    }
void init ( ) [inline, virtual]

Book histograms and initialise projections before the run.

Reimplemented from Analysis.

Definition at line 23 of file ATLAS_2012_I1204447.cc.

References ATLAS_2012_I1204447::_eventCountsPerSR, ATLAS_2012_I1204447::_h_e_n, ATLAS_2012_I1204447::_h_excluded, ATLAS_2012_I1204447::_h_HTjets_all, ATLAS_2012_I1204447::_h_HTlep_all, ATLAS_2012_I1204447::_h_Meff_all, ATLAS_2012_I1204447::_h_MET_all, ATLAS_2012_I1204447::_h_mu_n, ATLAS_2012_I1204447::_h_pt_1_2ltau, ATLAS_2012_I1204447::_h_pt_1_3l, ATLAS_2012_I1204447::_h_pt_2_2ltau, ATLAS_2012_I1204447::_h_pt_2_3l, ATLAS_2012_I1204447::_h_pt_3_2ltau, ATLAS_2012_I1204447::_h_pt_3_3l, ATLAS_2012_I1204447::_h_tau_n, ATLAS_2012_I1204447::_signal_regions, ATLAS_2012_I1204447::_use_fiducial_lepton_efficiency, Rivet::Cuts::abseta, IdentifiedFinalState::acceptIdPair(), ProjectionApplier::addProjection(), VetoedFinalState::addVetoPairId(), FastJets::ANTIKT, Analysis::bookHisto1D(), Rivet::PID::ELECTRON, ATLAS_2012_I1204447::getSignalRegions(), Rivet::GeV, and Rivet::PID::MUON.

                {

      // To calculate the acceptance without having the fiducial lepton efficiencies included, this part can be turned off
      _use_fiducial_lepton_efficiency = true;

      // Random numbers for simulation of ATLAS detector reconstruction efficiency
      srand(160385);

      // Read in all signal regions
      _signal_regions = getSignalRegions();

      // Set number of events per signal region to 0
      for (size_t i = 0; i < _signal_regions.size(); i++)
        _eventCountsPerSR[_signal_regions[i]] = 0.0;

      // Final state including all charged and neutral particles
      const FinalState fs(-5.0, 5.0, 1*GeV);
      addProjection(fs, "FS");

      // Final state including all charged particles
      addProjection(ChargedFinalState(Cuts::abseta < 2.5 && Cuts::pT > 1*GeV), "CFS");

      // Final state including all visible particles (to calculate MET, Jets etc.)
      addProjection(VisibleFinalState(Cuts::abseta < 5.0), "VFS");

      // Final state including all AntiKt 04 Jets
      VetoedFinalState vfs;
      vfs.addVetoPairId(PID::MUON);
      addProjection(FastJets(vfs, FastJets::ANTIKT, 0.4), "AntiKtJets04");

      // Final state including all unstable particles (including taus)
      addProjection(UnstableFinalState(Cuts::abseta < 5.0 && Cuts::pT > 5*GeV), "UFS");

      // Final state including all electrons
      IdentifiedFinalState elecs(Cuts::abseta < 2.47 && Cuts::pT > 10*GeV);
      elecs.acceptIdPair(PID::ELECTRON);
      addProjection(elecs, "elecs");

      // Final state including all muons
      IdentifiedFinalState muons(Cuts::abseta < 2.5 && Cuts::pT > 10*GeV);
      muons.acceptIdPair(PID::MUON);
      addProjection(muons, "muons");

      // Book histograms
      _h_HTlep_all  = bookHisto1D("HTlep_all" , 30, 0, 1500);
      _h_HTjets_all = bookHisto1D("HTjets_all", 30, 0, 1500);
      _h_MET_all    = bookHisto1D("MET_all"   , 20, 0, 1000);
      _h_Meff_all   = bookHisto1D("Meff_all"  , 30, 0, 3000);

      _h_e_n        = bookHisto1D("e_n"  , 10, -0.5, 9.5);
      _h_mu_n       = bookHisto1D("mu_n" , 10, -0.5, 9.5);
      _h_tau_n      = bookHisto1D("tau_n", 10, -0.5, 9.5);

      _h_pt_1_3l    = bookHisto1D("pt_1_3l", 100, 0, 2000);
      _h_pt_2_3l    = bookHisto1D("pt_2_3l", 100, 0, 2000);
      _h_pt_3_3l    = bookHisto1D("pt_3_3l", 100, 0, 2000);
      _h_pt_1_2ltau = bookHisto1D("pt_1_2ltau", 100, 0, 2000);
      _h_pt_2_2ltau = bookHisto1D("pt_2_2ltau", 100, 0, 2000);
      _h_pt_3_2ltau = bookHisto1D("pt_3_2ltau", 100, 0, 2000);

      _h_excluded   = bookHisto1D("excluded", 2, -0.5, 1.5);
    }
virtual std::string inspireId ( ) const [inline, virtual, inherited]

Get the Inspire ID code for this analysis.

Definition at line 119 of file Analysis.hh.

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

                                        {
      return info().inspireId();
    }
void integrate ( Histo1DPtr  h,
Scatter2DPtr  s 
) const [inherited]

Helper for converting a differential histo to an integral one.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.
Todo:
2D versions of scale and normalize...

Definition at line 700 of file Analysis.cc.

                                                             {
    // preserve the path info
    const string path = s->path();
    *s = toIntegralHisto(*h);
    s->setPath(path);
  }
void integrate ( const Histo1D &  h,
Scatter2DPtr  s 
) const [inherited]

Helper for converting a differential histo to an integral one.

Note:
Assigns to the (already registered) output scatter, s. Preserves the path information of the target.

Definition at line 707 of file Analysis.cc.

                                                                 {
    // preserve the path info
    const string path = s->path();
    *s = toIntegralHisto(h);
    s->setPath(path);
  }
bool isCompatible ( const ParticlePair beams) const [inherited]

Check if analysis is compatible with the provided beam particle IDs and energies.

Definition at line 101 of file Analysis.cc.

Referenced by Analysis::isCompatible().

                                                             {
    return isCompatible(beams.first.pid(),  beams.second.pid(),
                        beams.first.energy(), beams.second.energy());
  }
bool isCompatible ( PdgId  beam1,
PdgId  beam2,
double  e1,
double  e2 
) const [inherited]

Check if analysis is compatible with the provided beam particle IDs and energies.

Definition at line 107 of file Analysis.cc.

References Analysis::beams(), and Analysis::isCompatible().

                                                                                  {
    PdgIdPair beams(beam1, beam2);
    pair<double,double> energies(e1, e2);
    return isCompatible(beams, energies);
  }
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.

int isonZ ( const Particles particles) [inline]

Function checking if there is an OSSF lepton pair or a combination of 3 leptons with an invariant mass close to the Z mass

Todo:
Should the reference Z mass be 91.2?

Definition at line 996 of file ATLAS_2012_I1204447.cc.

References Rivet::GeV, Rivet::Cuts::mass, Particle::momentum(), and Particle::pid().

Referenced by ATLAS_2012_I1204447::analyze().

                                            {
      int onZ = 0;
      double best_mass_2 = 999.;
      double best_mass_3 = 999.;

      // Loop over all 2 particle combinations to find invariant mass of OSSF pair closest to Z mass
      foreach ( const Particle& p1, particles  )  {
        foreach ( const Particle& p2, particles  )  {
          double mass_difference_2_old = fabs(91.0 - best_mass_2);
          double mass_difference_2_new = fabs(91.0 - (p1.momentum() + p2.momentum()).mass()/GeV);

          // If particle combination is OSSF pair calculate mass difference to Z mass
          if ( (p1.pid()*p2.pid() == -121 || p1.pid()*p2.pid() == -169) )  {

            // Get invariant mass closest to Z mass
            if (mass_difference_2_new < mass_difference_2_old)
              best_mass_2 = (p1.momentum() + p2.momentum()).mass()/GeV;

            // In case there is an OSSF pair take also 3rd lepton into account (e.g. from FSR and photon to electron conversion)
            foreach ( const Particle & p3 , particles  )  {
              double mass_difference_3_old = fabs(91.0 - best_mass_3);
              double mass_difference_3_new = fabs(91.0 - (p1.momentum() + p2.momentum() + p3.momentum()).mass()/GeV);
              if (mass_difference_3_new < mass_difference_3_old)
                best_mass_3 = (p1.momentum() + p2.momentum() + p3.momentum()).mass()/GeV;
            }
          }
        }
      }

      // Pick the minimum invariant mass of the best OSSF pair combination and the best 3 lepton combination
      // If this mass is in a 20 GeV window around the Z mass, the event is classified as onZ
      double best_mass = min(best_mass_2, best_mass_3);
      if (fabs(91.0 - best_mass) < 20) onZ = 1;
      return onZ;
    }
const string makeAxisCode ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [protected, inherited]

Get the internal histogram name for given d, x and y (cf. HepData)

Definition at line 67 of file Analysis.cc.

Referenced by Analysis::bookHisto1D(), Analysis::bookProfile1D(), Analysis::bookScatter2D(), Analysis::getHisto1D(), Analysis::getProfile1D(), Analysis::getScatter2D(), Analysis::histoPath(), and Analysis::refData().

                                                                                                              {
    stringstream axisCode;
    axisCode << "d";
    if (datasetId < 10) axisCode << 0;
    axisCode << datasetId;
    axisCode << "-x";
    if (xAxisId < 10) axisCode << 0;
    axisCode << xAxisId;
    axisCode << "-y";
    if (yAxisId < 10) axisCode << 0;
    axisCode << yAxisId;
    return axisCode.str();
  }
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 114 of file Analysis.hh.

References Analysis::_defaultname, Analysis::info(), and AnalysisInfo::name().

Referenced by Analysis::_cacheRefData(), Analysis::Analysis(), Analysis::bookHisto1D(), Analysis::bookHisto2D(), Analysis::bookProfile1D(), Analysis::bookProfile2D(), Analysis::bookScatter2D(), Analysis::crossSection(), Analysis::getHisto1D(), Analysis::getLog(), Analysis::getProfile1D(), Analysis::getScatter2D(), Analysis::histoDir(), Analysis::normalize(), Analysis::refData(), and Analysis::scale().

                                   {
      return (info().name().empty()) ? _defaultname : info().name();
    }
bool needsCrossSection ( ) const [inline, inherited]

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

Todo:
Remove this and require HepMC >= 2.06

Definition at line 229 of file Analysis.hh.

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

                                   {
      return info().needsCrossSection();
    }
void normalize ( Histo1DPtr  histo,
double  norm = 1.0,
bool  includeoverflows = true 
) [inherited]

Normalize the given histogram, histo, to area = norm.

Note:
The histogram is no longer invalidated by this procedure.

Definition at line 652 of file Analysis.cc.

References MSG_ERROR, MSG_TRACE, MSG_WARNING, and Analysis::name().

Referenced by CMS_2011_S8950903::finalize(), TOTEM_2012_002::finalize(), CMS_2011_S8968497::finalize(), CMS_2012_I1090423::finalize(), D0_2010_S8671338::finalize(), D0_2007_S7075677::finalize(), MC_JETTAGS::finalize(), CDF_2009_NOTE_9936::finalize(), D0_2008_S7554427::finalize(), MC_HINC::finalize(), MC_ZJETS::finalize(), EXAMPLE_CUTS::finalize(), MC_HJETS::finalize(), BABAR_2005_S6181155::finalize(), CMS_2013_I1265659::finalize(), MC_ZINC::finalize(), ALICE_2010_S8624100::finalize(), ATLAS_2011_S8971293::finalize(), ALICE_2010_S8625980::finalize(), CDF_1996_S3418421::finalize(), D0_2009_S8320160::finalize(), CDF_1993_S2742446::finalize(), CMS_2011_S8957746::finalize(), D0_1996_S3324664::finalize(), ATLAS_2012_I1119557::finalize(), CMS_2013_I1273574::finalize(), DELPHI_1995_S3137023::finalize(), D0_2010_S8821313::finalize(), SLD_2002_S4869273::finalize(), MC_ZZJETS::finalize(), CMS_2013_I1272853::finalize(), ATLAS_2011_I925932::finalize(), D0_1996_S3214044::finalize(), D0_2008_S6879055::finalize(), CDF_1996_S3108457::finalize(), ATLAS_2014_I1300647::finalize(), ATLAS_2012_I1204784::finalize(), EXAMPLE::finalize(), CDF_1994_S2952106::finalize(), CMS_2013_I1224539_DIJET::finalize(), JADE_1998_S3612880::finalize(), ALEPH_2001_S4656318::finalize(), DELPHI_2002_069_CONF_603::finalize(), D0_2004_S5992206::finalize(), OPAL_1997_S3396100::finalize(), MC_PHOTONS::finalize(), MC_HFJETS::finalize(), MC_GENERIC::finalize(), D0_2001_S4674421::finalize(), MC_ZZINC::finalize(), OPAL_2001_S4553896::finalize(), ATLAS_2013_I1217867::finalize(), CMS_2012_I1107658::finalize(), TASSO_1990_S2148048::finalize(), CMS_2013_I1209721::finalize(), ATLAS_2013_I1217863_Z::finalize(), STAR_2008_S7869363::finalize(), ATLAS_2013_I1217863_W::finalize(), OPAL_1998_S3780481::finalize(), CDF_1997_S3541940::finalize(), CMS_2013_I1224539_WJET::finalize(), CMS_2013_I1256943::finalize(), CMS_2013_I1224539_ZJET::finalize(), DELPHI_2003_WUD_03_11::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1298811::finalize(), ATLAS_2012_I1124167::finalize(), CMS_2014_I1303894::finalize(), H1_1994_S2919893::finalize(), CDF_2001_S4751469::finalize(), MC_TTBAR::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2013_I1190187::finalize(), CDF_2004_S5839831::finalize(), CDF_1996_S3349578::finalize(), ATLAS_2012_I1094564::finalize(), ATLAS_2012_I1203852::finalize(), ALEPH_1996_S3486095::finalize(), DELPHI_1996_S3430090::finalize(), ATLAS_2014_I1279489::finalizePlots(), CMS_2013_I1258128::normalizeByContents(), and ATLAS_2011_I919017::safeinvscale().

                                                                               {
    if (!histo) {
      MSG_ERROR("Failed to normalize histo=NULL in analysis " << name() << " (norm=" << norm << ")");
      return;
    }
    MSG_TRACE("Normalizing histo " << histo->path() << " to " << norm);
    try {
      histo->normalize(norm, includeoverflows);
    } catch (YODA::Exception& we) {
      MSG_WARNING("Could not normalize histo " << histo->path());
      return;
    }
  }
void normalize ( Histo2DPtr  histo,
double  norm = 1.0,
bool  includeoverflows = true 
) [inherited]

Normalize the given histogram, histo, to area = norm.

Note:
The histogram is no longer invalidated by this procedure.
const Scatter2D & refData ( const string &  hname) const [protected, inherited]
const Scatter2D & refData ( unsigned int  datasetId,
unsigned int  xAxisId,
unsigned int  yAxisId 
) const [protected, inherited]

Get reference data for a numbered histo.

Definition at line 194 of file Analysis.cc.

References Analysis::makeAxisCode(), and Analysis::refData().

                                                                                                             {
    const string hname = makeAxisCode(datasetId, xAxisId, yAxisId);
    return refData(hname);
  }
virtual std::vector<std::string> references ( ) const [inline, virtual, inherited]

Journal, and preprint references.

Definition at line 180 of file Analysis.hh.

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

                                                    {
      return info().references();
    }
void removeAnalysisObject ( const std::string &  path) [protected, inherited]

Unregister a data object from the histogram system (by name)

Definition at line 725 of file Analysis.cc.

References Analysis::_analysisobjects.

Referenced by ATLAS_2014_I1279489::finalizeEfficiencies().

                                                        {
    for (vector<AnalysisObjectPtr>::iterator it = _analysisobjects.begin();  it != _analysisobjects.end(); ++it) {
      if ((*it)->path() == path) {
        _analysisobjects.erase(it);
        break;
      }
    }
  }
void removeAnalysisObject ( AnalysisObjectPtr  ao) [protected, inherited]

Unregister a data object from the histogram system (by pointer)

Definition at line 734 of file Analysis.cc.

References Analysis::_analysisobjects.

                                                          {
    for (vector<AnalysisObjectPtr>::iterator it = _analysisobjects.begin();  it != _analysisobjects.end(); ++it) {
      if (*it == ao) {
        _analysisobjects.erase(it);
        break;
      }
    }
 }
virtual const std::vector<PdgIdPair>& requiredBeams ( ) const [inline, virtual, inherited]

Return the allowed pairs of incoming beams required by this analysis.

Definition at line 206 of file Analysis.hh.

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

                                                              {
      return info().beams();
    }
virtual const std::vector<std::pair<double, double> >& requiredEnergies ( ) const [inline, virtual, inherited]

Sets of valid beam energy pairs, in GeV.

Definition at line 217 of file Analysis.hh.

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

Referenced by Analysis::setRequiredEnergies().

                                                                                {
      return info().energies();
    }
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 160 of file Analysis.hh.

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

                                      {
      return info().runInfo();
    }
void scale ( Histo1DPtr  histo,
double  scale 
) [inherited]

Multiplicatively scale the given histogram, histo, by factor scale.

Note:
The histogram is no longer invalidated by this procedure.

Definition at line 667 of file Analysis.cc.

References MSG_ERROR, MSG_TRACE, MSG_WARNING, and Analysis::name().

Referenced by MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), ALEPH_1991_S2435284::finalize(), D0_2000_S4480767::finalize(), TOTEM_2012_I1115294::finalize(), E735_1998_S3905616::finalize(), D0_2010_S8671338::finalize(), ATLAS_2010_S8591806::finalize(), CDF_2001_S4563131::finalize(), MC_XS::finalize(), UA5_1987_S1640666::finalize(), CDF_2007_S7057202::finalize(), CDF_2000_S4155203::finalize(), CDF_2008_S8093652::finalize(), BELLE_2013_I1216515::finalize(), BABAR_2007_S6895344::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), BELLE_2008_I786560::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1184941::finalize(), LHCB_2012_I1208102::finalize(), CDF_1990_S2089246::finalize(), CMS_2012_I1087342::finalize(), CMS_2010_S8656010::finalize(), CMS_2011_I954992::finalize(), CDF_1988_S1865951::finalize(), ATLAS_2012_I1118269::finalize(), MC_WJETS::finalize(), UA5_1982_S875503::finalize(), BABAR_2005_S6181155::finalize(), ATLAS_2011_I894867::finalize(), CMSTOTEM_2014_I1294140::finalize(), CDF_2009_S8383952::finalize(), STAR_2006_S6870392::finalize(), CMS_QCD_10_024::finalize(), LHCB_2010_I867355::finalize(), BELLE_2001_S4598261::finalize(), OPAL_1995_S3198391::finalize(), CMS_2010_S8547297::finalize(), OPAL_1994_S2927284::finalize(), OPAL_1997_S3608263::finalize(), OPAL_2000_S4418603::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), LHCB_2011_I919315::finalize(), ALICE_2010_S8625980::finalize(), D0_2011_I895662::finalize(), MC_PDFS::finalize(), ALEPH_2002_S4823664::finalize(), ALICE_2010_S8706239::finalize(), OPAL_1996_S3257789::finalize(), UA5_1989_S1926373::finalize(), ATLAS_2014_I1288706::finalize(), CDF_2005_S6080774::finalize(), DELPHI_1999_S3960137::finalize(), ATLAS_2014_I1282441::finalize(), ALICE_2011_S8945144::finalize(), MC_IDENTIFIED::finalize(), OPAL_1998_S3702294::finalize(), BABAR_2013_I1238276::finalize(), D0_2006_S6438750::finalize(), MC_DIPHOTON::finalize(), CMS_2013_I1273574::finalize(), ALICE_2011_S8909580::finalize(), STAR_2006_S6500200::finalize(), UA5_1986_S1583476::finalize(), MC_PHOTONINC::finalize(), ATLAS_2011_S9131140::finalize(), CMS_2011_S8978280::finalize(), D0_2000_I499943::finalize(), H1_1995_S3167097::finalize(), BABAR_2003_I593379::finalize(), PDG_TAUS::finalize(), ALICE_2012_I1181770::finalize(), D0_2008_S7863608::finalize(), D0_2008_S7662670::finalize(), MC_WINC::finalize(), ARGUS_1993_S2653028::finalize(), CLEO_2004_S5809304::finalize(), ATLAS_2011_S9035664::finalize(), D0_2009_S8202443::finalize(), D0_2008_S6879055::finalize(), CMS_2011_S9215166::finalize(), MC_PHOTONJETS::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), ALEPH_1999_S4193598::finalize(), BABAR_2007_S7266081::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), SLD_1996_S3398250::finalize(), CMS_2013_I1261026::finalize(), ATLAS_2011_I930220::finalize(), ATLAS_2013_I1219109::finalize(), SFM_1984_S1178091::finalize(), CMS_2012_I941555::finalize(), OPAL_1998_S3749908::finalize(), D0_2009_S8349509::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), ATLAS_2012_I1082009::finalize(), CDF_2008_S7541902::finalize(), ARGUS_1993_S2669951::finalize(), ATLAS_2010_S8919674::finalize(), UA1_1990_S2044935::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_WPOL::finalize(), MC_WWINC::finalize(), CMS_2012_I1107658::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ATLAS_2013_I1217863_Z::finalize(), ALEPH_1996_S3196992::finalize(), OPAL_1993_S2692198::finalize(), D0_2008_S7719523::finalize(), ATLAS_2013_I1217863_W::finalize(), CDF_2008_S8095620::finalize(), OPAL_1998_S3780481::finalize(), STAR_2006_S6860818::finalize(), CMS_2013_I1218372::finalize(), ATLAS_2012_I1091481::finalize(), ATLAS_2011_S9128077::finalize(), ATLAS_2011_I926145::finalize(), LHCB_2014_I1281685::finalize(), ATLAS_2010_S8918562::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_I1183818::finalize(), ATLAS_2012_CONF_2012_105::finalize(), OPAL_2004_S6132243::finalize(), ARGUS_1993_S2789213::finalize(), ATLAS_2012_CONF_2012_104::finalize(), H1_1994_S2919893::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2011_I944826::finalize(), H1_2000_S4129130::finalize(), ATLAS_2013_I1230812::finalize(), ATLAS_2012_I1117704::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2011_I945498::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2014_I1306615::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ALEPH_1996_S3486095::finalize(), ATLAS_2012_CONF_2012_153::finalize(), DELPHI_1996_S3430090::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), ATLAS_2012_I1125575::finalizeHistograms(), ATLAS_2011_I919017::safeinvscale(), and BinnedHistogram< T >::scale().

                                                     {
    if (!histo) {
      MSG_ERROR("Failed to scale histo=NULL in analysis " << name() << " (scale=" << scale << ")");
      return;
    }
    if (std::isnan(scale) || std::isinf(scale)) {
      MSG_ERROR("Failed to scale histo=" << histo->path() << " in analysis: " << name() << " (invalid scale factor = " << scale << ")");
      scale = 0;
    }
    MSG_TRACE("Scaling histo " << histo->path() << " by factor " << scale);
    try {
      histo->scaleW(scale);
    } catch (YODA::Exception& we) {
      MSG_WARNING("Could not scale histo " << histo->path());
      return;
    }
    // // Transforming the histo into a scatter after scaling
    // vector<double> x, y, ex, ey;
    // for (size_t i = 0, N = histo->numBins(); i < N; ++i) {
    //   x.push_back( histo->bin(i).midpoint() );
    //   ex.push_back(histo->bin(i).width()*0.5);
    //   y.push_back(histo->bin(i).height()*scale);
    //   ey.push_back(histo->bin(i).heightErr()*scale);
    // }
    // string title = histo->title();
    // Scatter2DPtr dps( new Scatter2D(x, y, ex, ey, hpath, title) );
    // addAnalysisObject(dps);
  }
void scale ( Histo2DPtr  histo,
double  scale 
) [inherited]

Multiplicatively scale the given histogram, histo, by factor scale.

Note:
The histogram is no longer invalidated by this procedure.
Analysis & setCrossSection ( double  xs) [inherited]

Set the cross section from the generator.

Definition at line 149 of file Analysis.cc.

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

                                               {
    _crossSection = xs;
    _gotCrossSection = true;
    return *this;
  }
Analysis& setNeedsCrossSection ( bool  needed = true) [inline, inherited]

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

Todo:
Remove this and require HepMC >= 2.06

Definition at line 234 of file Analysis.hh.

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

Referenced by ATLAS_2011_I954993::ATLAS_2011_I954993(), ATLAS_2014_I1306294::ATLAS_2014_I1306294(), ATLAS_2014_I1312627::ATLAS_2014_I1312627(), MC_JetAnalysis::MC_JetAnalysis(), MC_JetSplittings::MC_JetSplittings(), and MC_ParticleAnalysis::MC_ParticleAnalysis().

                                                     {
      info().setNeedsCrossSection(needed);
      return *this;
    }
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 210 of file Analysis.hh.

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

                                                                                  {
      info().setBeams(requiredBeams);
      return *this;
    }
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 221 of file Analysis.hh.

References Analysis::info(), Analysis::requiredEnergies(), and AnalysisInfo::setEnergies().

                                                                                                       {
      info().setEnergies(requiredEnergies);
      return *this;
    }
virtual std::string spiresId ( ) const [inline, virtual, inherited]

Get the SPIRES ID code for this analysis (~deprecated).

Definition at line 124 of file Analysis.hh.

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

                                       {
      return info().spiresId();
    }
double sqrtS ( ) const [inherited]

Centre of mass energy for this run.

Definition at line 27 of file Analysis.cc.

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

Referenced by ATLAS_2011_I894867::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), CMS_2012_I1193338::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), TOTEM_2012_002::analyze(), CMS_2012_I1184941::analyze(), SLD_2004_S5693039::analyze(), CMS_2011_I954992::analyze(), ALICE_2012_I1181770::analyze(), CMS_2011_S9215166::analyze(), SFM_1984_S1178091::analyze(), ALICE_2010_S8625980::analyze(), UA1_1990_S2044935::analyze(), CMS_2013_I1218372::analyze(), JADE_1998_S3612880::analyze(), TASSO_1990_S2148048::analyze(), ATLAS_2013_I1190187::analyze(), MC_TTBAR::analyze(), ATLAS_2010_S8894728::analyze(), ATLAS_2010_S8918562::analyze(), ALEPH_2004_S5765862::analyze(), CDF_2004_S5839831::analyze(), ATLAS_2012_I1203852::analyze(), ALICE_2010_S8625980::finalize(), DELPHI_2000_S4328825::finalize(), OPAL_2002_S5361494::finalize(), JADE_1998_S3612880::finalize(), UA1_1990_S2044935::finalize(), ATLAS_2010_S8918562::finalize(), ATLAS_2012_I1093734::finalize(), ALEPH_2004_S5765862::finalize(), CDF_2004_S5839831::finalize(), PDG_HADRON_MULTIPLICITIES_RATIOS::finalize(), PDG_HADRON_MULTIPLICITIES::finalize(), ATLAS_2010_S8894728::init(), CMS_2010_S8547297::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMSTOTEM_2014_I1294140::init(), CMS_2011_S9215166::init(), ATLAS_2012_I1091481::init(), CMS_2011_S8978280::init(), CDF_2012_NOTE10874::init(), CMS_QCD_10_024::init(), SFM_1984_S1178091::init(), MC_HINC::init(), UA5_1989_S1926373::init(), MC_JetSplittings::init(), MC_WWINC::init(), MC_ZZINC::init(), ATLAS_2011_I944826::init(), MC_WINC::init(), MC_DIPHOTON::init(), MC_ZINC::init(), MC_ZZJETS::init(), MC_PHOTONINC::init(), CDF_1988_S1865951::init(), UA5_1986_S1583476::init(), MC_ParticleAnalysis::init(), CDF_1990_S2089246::init(), MC_JetAnalysis::init(), MC_WWJETS::init(), UA1_1990_S2044935::init(), ALEPH_2004_S5765862::init(), JADE_1998_S3612880::init(), JADE_OPAL_2000_S4300807::init(), ALICE_2010_S8624100::init(), ALICE_2010_S8625980::init(), CMS_2013_I1218372::init(), LHCB_2012_I1119400::init(), UA5_1988_S1867512::init(), LHCB_2011_I917009::init(), MC_PHOTONS::init(), MC_WPOL::init(), TASSO_1990_S2148048::init(), ATLAS_2010_S8918562::init(), OPAL_2004_S6132243::init(), CDF_2004_S5839831::init(), ATLAS_2012_I1094061::init(), ATLAS_2012_I1093734::init(), PDG_HADRON_MULTIPLICITIES::init(), and PDG_HADRON_MULTIPLICITIES_RATIOS::init().

                               {
    return handler().sqrtS();
  }
virtual std::string status ( ) const [inline, virtual, inherited]

Whether this analysis is trusted (in any way!)

Definition at line 195 of file Analysis.hh.

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

                                     {
      return (info().status().empty()) ? "UNVALIDATED" : info().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 141 of file Analysis.hh.

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

                                      {
      return info().summary();
    }
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 93 of file Analysis.cc.

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

Referenced by Analysis::crossSectionPerEvent(), MC_JetSplittings::finalize(), MC_JetAnalysis::finalize(), MC_ParticleAnalysis::finalize(), CDF_2006_S6450792::finalize(), ALEPH_1991_S2435284::finalize(), D0_2000_S4480767::finalize(), CMS_2011_S9086218::finalize(), CDF_2001_S4563131::finalize(), CDF_2008_S7828950::finalize(), D0_2010_S8671338::finalize(), LHCF_2012_I1115479::finalize(), MC_XS::finalize(), CDF_2008_S8093652::finalize(), CDF_2000_S4155203::finalize(), BELLE_2013_I1216515::finalize(), BABAR_2007_S6895344::finalize(), CDF_2000_S4266730::finalize(), CMS_2011_S8941262::finalize(), CDF_1998_S3618439::finalize(), CMS_2012_I1184941::finalize(), CMS_2012_I1087342::finalize(), LHCB_2012_I1208102::finalize(), CMS_2010_S8656010::finalize(), CMS_2011_I954992::finalize(), ATLAS_2012_I1118269::finalize(), MC_WJETS::finalize(), BABAR_2005_S6181155::finalize(), ATLAS_2011_I894867::finalize(), CDF_2001_S4517016::finalize(), CDF_2009_S8383952::finalize(), STAR_2006_S6870392::finalize(), D0_2010_S8566488::finalize(), LHCB_2010_I867355::finalize(), OPAL_1995_S3198391::finalize(), CMS_2010_S8547297::finalize(), OPAL_1994_S2927284::finalize(), OPAL_1997_S3608263::finalize(), OPAL_2000_S4418603::finalize(), CDF_2009_S8436959::finalize(), CMS_2012_I1193338::finalize(), D0_2011_I895662::finalize(), MC_PDFS::finalize(), ALEPH_2002_S4823664::finalize(), OPAL_1996_S3257789::finalize(), CDF_2005_S6080774::finalize(), ATLAS_2014_I1288706::finalize(), ATLAS_2014_I1325553::finalize(), BABAR_2013_I1238276::finalize(), OPAL_1998_S3702294::finalize(), DELPHI_1999_S3960137::finalize(), ATLAS_2014_I1282441::finalize(), ALICE_2011_S8945144::finalize(), MC_IDENTIFIED::finalize(), D0_2006_S6438750::finalize(), MC_DIPHOTON::finalize(), CMS_2013_I1273574::finalize(), DELPHI_1995_S3137023::finalize(), ALICE_2011_S8909580::finalize(), CMS_2011_S8978280::finalize(), D0_2000_I499943::finalize(), ALICE_2012_I1181770::finalize(), CLEO_2004_S5809304::finalize(), D0_2008_S7662670::finalize(), D0_2008_S7863608::finalize(), MC_WINC::finalize(), ATLAS_2011_S9035664::finalize(), ATLAS_2010_CONF_2010_049::finalize(), ATLAS_2011_I954993::finalize(), MC_WWJETS::finalize(), CDF_2009_S8233977::finalize(), CMS_2014_I1298810::finalize(), LHCB_2013_I1218996::finalize(), ATLAS_2013_I1219109::finalize(), OPAL_1998_S3749908::finalize(), OPAL_1997_S3396100::finalize(), D0_2010_S8570965::finalize(), CMS_2011_S8973270::finalize(), ATLAS_2012_I946427::finalize(), CDF_2008_S7541902::finalize(), ATLAS_2010_S8919674::finalize(), ATLAS_2011_I921594::finalize(), ATLAS_2013_I1263495::finalize(), MC_WWINC::finalize(), ATLAS_2014_I1312627::finalize(), CDF_2006_S6653332::finalize(), ALEPH_1996_S3196992::finalize(), ATLAS_2013_I1217863_Z::finalize(), OPAL_1993_S2692198::finalize(), ATLAS_2013_I1217863_W::finalize(), D0_2008_S7719523::finalize(), OPAL_1998_S3780481::finalize(), STAR_2006_S6860818::finalize(), ATLAS_2014_I1307756::finalize(), ATLAS_2011_I926145::finalize(), CMS_2013_I1256943::finalize(), JADE_OPAL_2000_S4300807::finalize(), ATLAS_2014_I1306294::finalize(), ATLAS_2011_S9108483::finalize(), CDF_2008_S7540469::finalize(), ATLAS_2012_I1199269::finalize(), ATLAS_2010_S8914702::finalize(), ATLAS_2014_I1307243::finalize(), ATLAS_2012_I1083318::finalize(), ATLAS_2014_I1304688::finalize(), ATLAS_2012_CONF_2012_105::finalize(), CMS_2014_I1303894::finalize(), ATLAS_2012_CONF_2012_104::finalize(), ATLAS_2011_S9120807::finalize(), ATLAS_2012_I1186556::finalize(), ATLAS_2012_I1084540::finalize(), MC_VH2BB::finalize(), ATLAS_2012_CONF_2012_103::finalize(), ATLAS_2012_I1117704::finalize(), ALEPH_2004_S5765862::finalize(), ATLAS_2013_I1190187::finalize(), ATLAS_2012_I1093738::finalize(), SLD_2004_S5693039::finalize(), ATLAS_2012_I1125961::finalize(), ATLAS_2011_CONF_2011_098::finalize(), ATLAS_2011_S8983313::finalize(), ATLAS_2012_I1126136::finalize(), ATLAS_2012_I1095236::finalize(), ATLAS_2012_I1190891::finalize(), ATLAS_2011_S9019561::finalize(), ATLAS_2012_I1112263::finalize(), ATLAS_2011_S9212183::finalize(), ATLAS_2012_CONF_2012_109::finalize(), ATLAS_2011_CONF_2011_090::finalize(), ATLAS_2011_S9225137::finalize(), ATLAS_2012_CONF_2012_001::finalize(), ATLAS_2012_I1203852::finalize(), ATLAS_2012_I943401::finalize(), ATLAS_2012_I1180197::finalize(), ALEPH_1996_S3486095::finalize(), ATLAS_2012_CONF_2012_153::finalize(), SLD_1999_S3743934::finalize(), ATLAS_2012_I1204447::finalize(), ATLAS_2011_S9212353::finalize(), ATLAS_2011_S9041966::finalize(), and PDG_HADRON_MULTIPLICITIES::finalize().

                                      {
    return handler().sumOfWeights();
  }
virtual std::vector<std::string> todos ( ) const [inline, virtual, inherited]

Any work to be done on this analysis.

Definition at line 200 of file Analysis.hh.

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

                                               {
      return info().todos();
    }
virtual std::string year ( ) const [inline, virtual, inherited]

When the original experimental analysis was published.

Definition at line 175 of file Analysis.hh.

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

                                   {
      return info().year();
    }

Member Data Documentation

bool _allowProjReg [protected, inherited]

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

Definition at line 143 of file ProjectionApplier.hh.

Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().

vector<string> _signal_regions [private]

List of signal regions and event counts per signal region.

Definition at line 1049 of file ATLAS_2012_I1204447.cc.

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

Fiducial efficiencies to model the effects of the ATLAS detector.

Definition at line 1046 of file ATLAS_2012_I1204447.cc.

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


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