rivet is hosted by Hepforge, IPPP Durham

Project out the last pre-decay b and c hadrons. More...

#include <HeavyHadrons.hh>

Inheritance diagram for HeavyHadrons:
Collaboration diagram for HeavyHadrons:

List of all members.

Public Member Functions

virtual int compare (const Projection &p) const
 Compare projections.
virtual bool accept (const Particle &p) const
 Decide if a particle is to be accepted or not.
bool before (const Projection &p) const
virtual const std::set< PdgIdPairbeamPairs () const
virtual std::string name () const
 Get the name of the projection.
ProjectionaddPdgIdPair (PdgId beam1, PdgId beam2)
 Add a colliding beam pair.
LoggetLog () const
 Get a Log object based on the getName() property of the calling projection object.
void setName (const std::string &name)
 Used by derived classes to set their name.
Constructors and destructors.
 HeavyHadrons (double mineta=-MAXDOUBLE, double maxeta=MAXDOUBLE, double minpt=0.0 *GeV)
virtual const Projectionclone () const
 Clone on the heap.
Particle accessors
const ParticlesbHadrons () const
 Get all weakly decaying b hadrons (return by reference)
Particles bHadrons (double pTmin) const
 Get weakly decaying b hadrons with a pTmin cut (return by value)
const ParticlescHadrons () const
 Get all weakly decaying c hadrons (return by reference)
const Particles cHadrons (double pTmin) const
 Get weakly decaying c hadrons with a pTmin cut (return by value)
Particle accessors
virtual size_t size () const
 Access the projected final-state particles.
virtual bool empty () const
 Is this final state empty?
virtual bool isEmpty () const
virtual const Particlesparticles () const
 Get the final-state particles in no particular order, with no cuts.
template<typename F >
Particles particles (Cut c, F sorter) const
 Get the final-state particles, with optional cuts.
Particles particlesByPt (Cut c=Cuts::open()) const
Particles particlesByP (Cut c=Cuts::open()) const
Particles particlesByE (Cut c=Cuts::open()) const
Particles particlesByEt (Cut c=Cuts::open()) const
Particles particlesByEta (Cut c=Cuts::open()) const
Particles particlesByModEta (Cut c=Cuts::open()) const
Particles particlesByRapidity (Cut c=Cuts::open()) const
Particles particlesByModRapidity (Cut c=Cuts::open()) 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

virtual void project (const Event &e)
 Apply the projection to the event.
Cmp< ProjectionmkNamedPCmp (const Projection &otherparent, const std::string &pname) const
Cmp< ProjectionmkPCmp (const Projection &otherparent, const std::string &pname) const
ProjectionHandlergetProjHandler () const
 Get a reference to the ProjectionHandler for this thread.
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

Particles _theBs
 b and c hadron containers
Particles _theCs
Cut _cuts
 The applicable cuts.
Particles _theParticles
 The final-state particles.
bool _allowProjReg
 Flag to forbid projection registration in analyses until the init phase.

Friends

class Event
 Event is a friend.
class Cmp< Projection >
 The Cmp specialization for Projection is a friend.

JetAlg compatibility

Minimum- $ p_\perp $ requirement.

Todo:
Replace with cuts() accessor virtual Cut cuts() const { return _cuts; }
const collection_typeentities () const
 Template-usable interface common to JetAlg.
typedef Particle entity_type
typedef Particles collection_type

Detailed Description

Project out the last pre-decay b and c hadrons.

This currently defines a c-hadron as one which contains a c quark and {not} a b quark.

Todo:
This assumes that the heavy hadrons are unstable... should we also look for stable ones in case the decays are disabled?

Definition at line 19 of file HeavyHadrons.hh.


Member Typedef Documentation

typedef Particles collection_type [inherited]

Definition at line 131 of file ParticleFinder.hh.

typedef Particle entity_type [inherited]

Definition at line 130 of file ParticleFinder.hh.


Constructor & Destructor Documentation

HeavyHadrons ( double  mineta = -MAXDOUBLE,
double  maxeta = MAXDOUBLE,
double  minpt = 0.0*GeV 
) [inline]

Constructor with specification of the minimum and maximum pseudorapidity $ \eta $ and the min $ p_T $ (in GeV).

Definition at line 27 of file HeavyHadrons.hh.

References ProjectionApplier::addProjection(), and Projection::setName().

Referenced by HeavyHadrons::clone().

                                           {
      setName("HeavyHadrons");
      addProjection(UnstableFinalState(mineta, maxeta, minpt), "UFS");
    }

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;
  }
bool accept ( const Particle p) const [virtual, inherited]

Decide if a particle is to be accepted or not.

Decide if a particle is to be accepted or not.

Todo:
Rename to _accept or acceptFinal?

Definition at line 99 of file FinalState.cc.

References ParticleFinder::_cuts, and Particle::genParticle().

Referenced by InvMassFinalState::calc(), FinalState::project(), LeadingParticlesFinalState::project(), and DressedLeptons::project().

                                                 {
    // Not having s.c. == 1 should never happen!
    assert(p.genParticle() == NULL || p.genParticle()->status() == 1);

    return _cuts->accept(p);
  }
Projection& addPdgIdPair ( PdgId  beam1,
PdgId  beam2 
) [inline, inherited]

Add a colliding beam pair.

Definition at line 107 of file Projection.hh.

References Projection::_beamPairs.

Referenced by Projection::Projection().

                                                       {
      _beamPairs.insert(PdgIdPair(beam1, beam2));
      return *this;
    }
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 TauFinder::_init(), 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(), ATLAS_2010_S8894728::init(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMS_2011_S8950903::init(), CMS_2012_PAS_QCD_11_010::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2010_CONF_2010_049::init(), ATLAS_2011_I894867::init(), LHCB_2013_I1208105::init(), CMS_2011_S8941262::init(), CMS_2011_S8968497::init(), CMS_2011_S8973270::init(), CMS_2011_S8978280::init(), ATLAS_2014_I1298811::init(), ALICE_2011_S8909580::init(), CMS_2012_I1090423::init(), CDF_2007_S7057202::init(), ALICE_2011_S8945144::init(), TOTEM_2012_I1115294::init(), ATLAS_2014_I1282441::init(), CMS_2011_S9086218::init(), CMS_2011_S9215166::init(), CMS_2012_I1087342::init(), CMS_2012_I1184941::init(), CMS_2012_I1193338::init(), ATLAS_2010_S8591806::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), LHCF_2012_I1115479::init(), ATLAS_2011_I930220::init(), CMS_2011_I954992::init(), TOTEM_2012_002::init(), CMS_2011_S8957746::init(), CMS_2011_S9088458::init(), ATLAS_2011_S9002537::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), CDF_2012_NOTE10874::init(), ATLAS_2012_I1124167::init(), STAR_2006_S6500200::init(), STAR_2008_S7993412::init(), UA5_1987_S1640666::init(), CMS_QCD_10_024::init(), CDF_1997_S3541940::init(), ATLAS_2011_S9128077::init(), BELLE_2008_I786560::init(), CDF_1993_S2742446::init(), ARGUS_1993_S2789213::init(), CDF_2000_S4155203::init(), CDF_2005_S6080774::init(), CDF_2006_S6450792::init(), SFM_1984_S1178091::init(), BABAR_2007_S7266081::init(), CMS_2013_I1261026::init(), D0_2008_S6879055::init(), ATLAS_2012_I1091481::init(), UA5_1982_S875503::init(), CDF_2008_S7782535::init(), CDF_2008_S8093652::init(), ATLAS_2012_I1204447::init(), STAR_2006_S6870392::init(), ARGUS_1993_S2669951::init(), UA5_1989_S1926373::init(), CMS_2013_I1272853::init(), D0_2000_I499943::init(), ALEPH_1991_S2435284::init(), E735_1998_S3905616::init(), ATLAS_2012_I1183818::init(), ATLAS_2013_I1243871::init(), ATLAS_2014_I1268975::init(), CMS_2012_I1102908::init(), ATLAS_2011_S9035664::init(), ALEPH_2002_S4823664::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), CDF_1988_S1865951::init(), CDF_1990_S2089246::init(), DELPHI_1999_S3960137::init(), ATLAS_2010_S8817804::init(), ATLAS_2011_S8924791::init(), ALEPH_1996_S3196992::init(), CDF_2005_S6217184::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(), UA5_1986_S1583476::init(), D0_2001_S4674421::init(), D0_2008_S7554427::init(), CDF_1994_S2952106::init(), H1_1995_S3167097::init(), ATLAS_2014_I1307756::init(), UA1_1990_S2044935::init(), ZEUS_2001_S4815815::init(), LHCB_2011_I919315::init(), ATLAS_2012_I1119557::init(), CDF_2008_S7540469::init(), CDF_2008_S7828950::init(), ATLAS_2012_I1188891::init(), JADE_1998_S3612880::init(), JADE_OPAL_2000_S4300807::init(), EXAMPLE::init(), OPAL_1997_S3396100::init(), ALEPH_2004_S5765862::init(), STAR_2006_S6860818::init(), D0_1996_S3214044::init(), STAR_2009_UE_HELEN::init(), ALICE_2010_S8624100::init(), CDF_2008_S8095620::init(), DELPHI_1995_S3137023::init(), ATLAS_2012_I1125575::init(), ALICE_2010_S8625980::init(), LHCB_2013_I1218996::init(), CDF_2009_NOTE_9936::init(), CMS_2013_I1218372::init(), ALICE_2010_S8706239::init(), D0_2009_S8320160::init(), CDF_2009_S8436959::init(), ATLAS_2013_I1263495::init(), CDF_1996_S3108457::init(), CDF_1998_S3618439::init(), CDF_2000_S4266730::init(), CDF_2001_S4563131::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), SLD_2002_S4869273::init(), D0_2010_S8566488::init(), D0_2006_S6438750::init(), CDF_1996_S3349578::init(), CDF_2001_S4517016::init(), DELPHI_2002_069_CONF_603::init(), ATLAS_2012_I1082936::init(), ALEPH_2001_S4656318::init(), LHCB_2012_I1119400::init(), ATLAS_2011_I921594::init(), CDF_1996_S3418421::init(), UA5_1988_S1867512::init(), LHCB_2011_I917009::init(), SLD_1996_S3398250::init(), CDF_2006_S6653332::init(), CDF_2008_S7541902::init(), D0_2008_S7662670::init(), ATLAS_2011_S9108483::init(), OPAL_2002_S5361494::init(), ATLAS_2012_I1082009::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I946427::init(), ATLAS_2012_I1199269::init(), DELPHI_2000_S4328825::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_DIJET::init(), CDF_2009_S8233977::init(), CDF_2010_S8591881_QCD::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2010_S8918562::init(), D0_2008_S7719523::init(), ATLAS_2011_S9120807::init(), CDF_2001_S4751469::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(), BELLE_2001_S4598261::init(), BABAR_2005_S6181155::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::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(), SLD_2004_S5693039::init(), BELLE_2006_S6265367::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(), 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 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));
    }
const set< PdgIdPair > beamPairs ( ) const [virtual, inherited]

Return the allowed beam pairs on which this projection can operate, not including recursion. Derived classes should ensure that all contained projections are registered in the _projections set for the beam constraint chaining to work.

Todo:
Remove the beam constraints system from projections.

Definition at line 33 of file Projection.cc.

References Projection::_beamPairs, Projection::beamPairs(), Projection::getLog(), ProjectionApplier::getProjections(), Rivet::intersection(), and Log::TRACE.

Referenced by Projection::beamPairs().

                                                   {
    set<PdgIdPair> ret = _beamPairs;
    set<ConstProjectionPtr> projs = getProjections();
    for (set<ConstProjectionPtr>::const_iterator ip = projs.begin(); ip != projs.end(); ++ip) {
      ConstProjectionPtr p = *ip;
      getLog() << Log::TRACE << "Proj addr = " << p << endl;
      if (p) ret = intersection(ret, p->beamPairs());
    }
    return ret;
  }
bool before ( const Projection p) const [inherited]

Determine whether this object should be ordered before the object p given as argument. If p is of a different class than this, the before() function of the corresponding type_info objects is used. Otherwise, if the objects are of the same class, the virtual compare(const Projection &) will be returned.

Definition at line 22 of file Projection.cc.

References Projection::compare().

Referenced by less< const Rivet::Projection * >::operator()().

                                                   {
    const std::type_info& thisid = typeid(*this);
    const std::type_info& otherid = typeid(p);
    if (thisid == otherid) {
      return compare(p) < 0;
    } else {
      return thisid.before(otherid);
    }
  }
const Particles& bHadrons ( ) const [inline]

Get all weakly decaying b hadrons (return by reference)

Definition at line 46 of file HeavyHadrons.hh.

References HeavyHadrons::_theBs.

Referenced by HeavyHadrons::bHadrons().

                                      {
      return _theBs;
    }
Particles bHadrons ( double  pTmin) const [inline]

Get weakly decaying b hadrons with a pTmin cut (return by value)

Definition at line 51 of file HeavyHadrons.hh.

References HeavyHadrons::bHadrons(), and ParticleBase::pT().

                                           {
      Particles rtn;
      foreach (const Particle& p, bHadrons())
        if (p.pT() > pTmin) rtn += p;
      return rtn;
    }
const Particles& cHadrons ( ) const [inline]

Get all weakly decaying c hadrons (return by reference)

Definition at line 59 of file HeavyHadrons.hh.

References HeavyHadrons::_theCs.

Referenced by HeavyHadrons::cHadrons().

                                      {
      return _theCs;
    }
const Particles cHadrons ( double  pTmin) const [inline]

Get weakly decaying c hadrons with a pTmin cut (return by value)

Definition at line 64 of file HeavyHadrons.hh.

References HeavyHadrons::cHadrons(), and ParticleBase::pT().

                                                 {
      Particles rtn;
      foreach (const Particle& p, cHadrons())
        if (p.pT() > pTmin) rtn += p;
      return rtn;
    }
virtual const Projection* clone ( ) const [inline, virtual]

Clone on the heap.

Reimplemented from FinalState.

Definition at line 35 of file HeavyHadrons.hh.

References HeavyHadrons::HeavyHadrons().

                                            {
      return new HeavyHadrons(*this);
    }
int compare ( const Projection p) const [virtual, inherited]

Compare projections.

Todo:
HOW DO WE COMPARE CUTS OBJECTS?

Reimplemented from ParticleFinder.

Reimplemented in VetoedFinalState, IdentifiedFinalState, WFinder, ZFinder, DressedLeptons, InvMassFinalState, LossyFinalState< FILTER >, LossyFinalState< ConstRandomFilter >, LeadingParticlesFinalState, PromptFinalState, VisibleFinalState, DISFinalState, HadronicFinalState, NonHadronicFinalState, NeutralFinalState, MergedFinalState, ChargedLeptons, and ChargedFinalState.

Definition at line 47 of file FinalState.cc.

References ParticleFinder::_cuts, Rivet::EQUIVALENT, and Rivet::UNDEFINED.

                                                   {
    const FinalState& other = dynamic_cast<const FinalState&>(p);

    //MSG_TRACE("FS::compare: " << 1 << " " << this << " " << &p);
    //    std::vector<std::pair<double, double> > eta1(_etaRanges);
    //std::vector<std::pair<double, double> > eta2(other._etaRanges);
    //std::sort(eta1.begin(), eta1.end());
    //std::sort(eta2.begin(), eta2.end());

    //MSG_TRACE("FS::compare: " << 2 << " " << this << " " << &p);
    //if (eta1 < eta2) return ORDERED;
    //else if (eta2 < eta1) return UNORDERED;

    //MSG_TRACE("FS::compare: " << 3 << " " << this << " " << &p);
    //return cmp(_ptmin, other._ptmin);
    return _cuts == other._cuts ? EQUIVALENT : UNDEFINED;
  }
const collection_type& entities ( ) const [inline, inherited]

Template-usable interface common to JetAlg.

Definition at line 134 of file ParticleFinder.hh.

References ParticleFinder::particles().

                                            {
      return particles();
    }
Log& getLog ( ) const [inline, inherited]

Get a Log object based on the getName() property of the calling projection object.

Reimplemented from ProjectionApplier.

Definition at line 114 of file Projection.hh.

References Projection::name().

Referenced by Projection::beamPairs(), InvMassFinalState::calc(), ChargedFinalState::project(), InitialQuarks::project(), UnstableFinalState::project(), PromptFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), and VetoedFinalState::project().

                        {
      string logname = "Rivet.Projection." + name();
      return Log::getLog(logname);
    }
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().

virtual bool isEmpty ( ) const [inline, virtual, inherited]
Deprecated:
Is this final state empty?

Definition at line 39 of file ParticleFinder.hh.

References ParticleFinder::_theParticles.

{ return _theParticles.empty(); }
Cmp< Projection > mkPCmp ( const Projection otherparent,
const std::string &  pname 
) const [protected, inherited]

Shortcut to make a named Cmp<Projection> comparison with the *this object automatically passed as one of the parent projections.

Definition at line 51 of file Projection.cc.

References Rivet::pcmp().

                                                                {
    return pcmp(*this, otherparent, pname);
  }
virtual const Particles& particles ( ) const [inline, virtual, inherited]

Get the final-state particles in no particular order, with no cuts.

Definition at line 42 of file ParticleFinder.hh.

References ParticleFinder::_theParticles.

Referenced by BABAR_2003_I593379::analyze(), BELLE_2001_S4598261::analyze(), BABAR_2005_S6181155::analyze(), CLEO_2004_S5809304::analyze(), BELLE_2006_S6265367::analyze(), BABAR_2007_S6895344::analyze(), ARGUS_1993_S2653028::analyze(), BABAR_2013_I1238276::analyze(), BELLE_2013_I1216515::analyze(), OPAL_1994_S2927284::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), OPAL_1993_S2692198::analyze(), TOTEM_2012_002::analyze(), H1_2000_S4129130::analyze(), BELLE_2008_I786560::analyze(), SLD_2004_S5693039::analyze(), TOTEM_2012_I1115294::analyze(), CMS_2010_S8656010::analyze(), OPAL_1998_S3780481::analyze(), ATLAS_2010_S8591806::analyze(), CMS_2012_PAS_QCD_11_010::analyze(), STAR_2008_S7993412::analyze(), LHCF_2012_I1115479::analyze(), ALICE_2011_S8909580::analyze(), SLD_1999_S3743934::analyze(), ALICE_2011_S8945144::analyze(), OPAL_2000_S4418603::analyze(), CMSTOTEM_2014_I1294140::analyze(), ARGUS_1993_S2669951::analyze(), OPAL_1995_S3198391::analyze(), OPAL_1997_S3608263::analyze(), ALEPH_2002_S4823664::analyze(), DELPHI_1999_S3960137::analyze(), OPAL_1996_S3257789::analyze(), OPAL_1998_S3702294::analyze(), E735_1998_S3905616::analyze(), H1_1994_S2919893::analyze(), ATLAS_2011_S9002537::analyze(), CMS_QCD_10_024::analyze(), BABAR_2007_S7266081::analyze(), DELPHI_1995_S3137023::analyze(), UA5_1982_S875503::analyze(), CMS_2010_S8547297::analyze(), SLD_2002_S4869273::analyze(), CDF_1990_S2089246::analyze(), CDF_1988_S1865951::analyze(), ALEPH_2001_S4656318::analyze(), ATLAS_2011_S9035664::analyze(), DELPHI_2002_069_CONF_603::analyze(), ALICE_2010_S8706239::analyze(), ATLAS_2011_I925932::analyze(), SLD_1996_S3398250::analyze(), DELPHI_2000_S4328825::analyze(), CDF_2008_S7540469::analyze(), CMS_2011_S8973270::analyze(), OPAL_2002_S5361494::analyze(), CMS_2011_S9215166::analyze(), SFM_1984_S1178091::analyze(), LHCB_2011_I919315::analyze(), OPAL_1998_S3749908::analyze(), CMS_2013_I1272853::analyze(), ALICE_2010_S8625980::analyze(), D0_2006_S6438750::analyze(), MC_PHOTONINC::analyze(), OPAL_1997_S3396100::analyze(), MC_PHOTONKTSPLITTINGS::analyze(), MC_DIJET::analyze(), MC_LEADJETUE::analyze(), STAR_2006_S6500200::analyze(), UA5_1986_S1583476::analyze(), CDF_2009_S8233977::analyze(), MC_IDENTIFIED::analyze(), CMS_2011_S8978280::analyze(), MC_PHOTONJETS::analyze(), STAR_2006_S6860818::analyze(), CMS_2011_S8884919::analyze(), ATLAS_2012_I1082009::analyze(), CDF_2008_S8095620::analyze(), LHCB_2010_S8758301::analyze(), LHCB_2011_I917009::analyze(), CMS_2013_I1256943::analyze(), STAR_2009_UE_HELEN::analyze(), CDF_2006_S6653332::analyze(), LHCB_2013_I1208105::analyze(), ARGUS_1993_S2789213::analyze(), ALEPH_1999_S4193598::analyze(), UA1_1990_S2044935::analyze(), D0_2001_S4674421::analyze(), H1_1995_S3167097::analyze(), CMS_2013_I1218372::analyze(), CDF_2008_S7541902::analyze(), LHCB_2013_I1218996::analyze(), ATLAS_2012_I1183818::analyze(), JADE_1998_S3612880::analyze(), MC_PHOTONJETUE::analyze(), CDF_2010_S8591881_QCD::analyze(), CDF_2010_S8591881_DY::analyze(), TASSO_1990_S2148048::analyze(), MC_GENERIC::analyze(), LHCB_2012_I1119400::analyze(), ATLAS_2012_I1084540::analyze(), STAR_2008_S7869363::analyze(), D0_2008_S7719523::analyze(), ATLAS_2011_I926145::analyze(), ATLAS_2010_S8894728::analyze(), DELPHI_2003_WUD_03_11::analyze(), MC_SUSY::analyze(), ATLAS_2012_I1093738::analyze(), ALEPH_1996_S3486095::analyze(), ATLAS_2011_I944826::analyze(), ALEPH_2004_S5765862::analyze(), ATLAS_2012_I1094061::analyze(), DELPHI_1996_S3430090::analyze(), ATLAS_2012_I1093734::analyze(), BeamThrust::calc(), FParameter::calc(), Spherocity::calc(), Thrust::calc(), Sphericity::calc(), ParticleFinder::entities(), CMS_2013_I1261026::eventDecomp(), ATLAS_2010_S8918562::fillPtEtaNch(), CMS_2013_I1258128::makePhotonCut(), ParticleFinder::particlesByE(), ParticleFinder::particlesByEt(), ParticleFinder::particlesByEta(), ParticleFinder::particlesByModEta(), ParticleFinder::particlesByModRapidity(), ParticleFinder::particlesByP(), ParticleFinder::particlesByPt(), ParticleFinder::particlesByRapidity(), ChargedFinalState::project(), ChargedLeptons::project(), CentralEtHCM::project(), MergedFinalState::project(), FinalState::project(), DISLepton::project(), TriggerCDFRun0Run1::project(), TriggerCDFRun2::project(), NeutralFinalState::project(), HadronicFinalState::project(), PrimaryHadrons::project(), NonHadronicFinalState::project(), TauFinder::project(), VisibleFinalState::project(), PromptFinalState::project(), TriggerUA5::project(), LossyFinalState< ConstRandomFilter >::project(), LeadingParticlesFinalState::project(), FoxWolframMoments::project(), MissingMomentum::project(), InvMassFinalState::project(), DressedLeptons::project(), HeavyHadrons::project(), Hemispheres::project(), IdentifiedFinalState::project(), and VetoedFinalState::project().

{ return _theParticles; }
Particles particles ( Cut  c,
sorter 
) const [inline, inherited]

Get the final-state particles, with optional cuts.

Note:
Returns a copy rather than a reference, due to cuts Get the final-state particles, ordered by supplied sorting function object.
Returns a copy rather than a reference, due to cuts and sorting
Todo:
Update to use Cuts

Definition at line 52 of file ParticleFinder.hh.

References ParticleFinder::_theParticles, Rivet::Cuts::open(), and ParticleFinder::size().

                                               {
      Particles result;
      result.reserve(size());
      if ( c == Cuts::open() )
        result.assign(_theParticles.begin(), _theParticles.end());
      else 
    foreach (const Particle& p, _theParticles)
          if ( c->accept(p) ) result.push_back(p);
      std::sort(result.begin(), result.end(), sorter);
      return result;
    }
Particles particlesByE ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by decreasing $ E $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 80 of file ParticleFinder.hh.

References Rivet::cmpMomByE(), and ParticleFinder::particles().

                                                     {
      return particles(c, cmpMomByE);
    }
Particles particlesByEt ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by decreasing $ E_T $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 87 of file ParticleFinder.hh.

References Rivet::cmpMomByEt(), and ParticleFinder::particles().

                                                      {
      return particles(c, cmpMomByEt);
    }
Particles particlesByEta ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by increasing $ \eta $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 94 of file ParticleFinder.hh.

References Rivet::cmpMomByAscPseudorapidity(), and ParticleFinder::particles().

Referenced by ATLAS_2011_I894867::analyze(), CMS_2012_I1193338::analyze(), CMS_2012_I1184941::analyze(), ALICE_2012_I1181770::analyze(), ATLAS_2012_I1091481::analyze(), and ATLAS_2012_I1084540::fillMap().

Particles particlesByModEta ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by increasing $ |\eta| $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 101 of file ParticleFinder.hh.

References Rivet::cmpMomByAscAbsPseudorapidity(), and ParticleFinder::particles().

Particles particlesByModRapidity ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by increasing $ |y| $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 115 of file ParticleFinder.hh.

References Rivet::cmpMomByAscAbsRapidity(), and ParticleFinder::particles().

Particles particlesByP ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by decreasing $ p $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 73 of file ParticleFinder.hh.

References Rivet::cmpMomByP(), and ParticleFinder::particles().

                                                     {
      return particles(c, cmpMomByP);
    }
Particles particlesByPt ( Cut  c = Cuts::open()) const [inline, inherited]
Particles particlesByRapidity ( Cut  c = Cuts::open()) const [inline, inherited]

Get the final-state particles, ordered by increasing $ y $.

Todo:

Update to use Cuts

Remove, since there is the templated method or sortByX methods available for these unusual cases?

Definition at line 108 of file ParticleFinder.hh.

References Rivet::cmpMomByAscRapidity(), and ParticleFinder::particles().

Referenced by ALICE_2012_I1181770::analyze(), and CMS_2013_I1218372::analyze().

                                                            {
      return particles(c, cmpMomByAscRapidity);
    }
void project ( const Event e) [protected, virtual]

Apply the projection to the event.

Todo:
Allow user to choose whether primary or final HF hadrons are to be returned
Todo:
What about charm hadrons coming from bottom hadron decays?

Reimplemented from FinalState.

Definition at line 7 of file HeavyHadrons.cc.

References HeavyHadrons::_theBs, HeavyHadrons::_theCs, ParticleFinder::_theParticles, Particle::genParticle(), Rivet::PID::hasBottom(), Rivet::PID::hasCharm(), Rivet::PID::isHadron(), MSG_DEBUG, ParticleFinder::particles(), Rivet::particles_out(), and Particle::pid().

                                           {
    _theParticles.clear();
    _theBs.clear();
    _theCs.clear();

    /// @todo Allow user to choose whether primary or final HF hadrons are to be returned

    const Particles& unstables = applyProjection<FinalState>(e, "UFS").particles();
    foreach (const Particle& p, unstables) {
      // Exclude non-b/c-hadrons
      if (!isHadron(p)) continue;
      if (!hasCharm(p) && !hasBottom(p)) continue;
      MSG_DEBUG("Found a heavy (b or c) unstable hadron: " << p.pid());

      // An unbound, or undecayed status 2 hadron: this is weird, but I guess is allowed...
      if (!p.genParticle() || !p.genParticle()->end_vertex()) {
        MSG_DEBUG("Heavy hadron " << p.pid() << " with no GenParticle or decay found");
        _theParticles.push_back(p);
        if (hasBottom(p)) _theBs.push_back(p); else _theCs.push_back(p);
        continue;
      }
      // There are descendants -- check them for b or c content
      /// @todo What about charm hadrons coming from bottom hadron decays?
      const vector<GenParticle*> children = particles_out(p.genParticle(), HepMC::children);
      if (hasBottom(p)) {
        bool has_b_child = false;
        foreach (const GenParticle* p2, children) {
          if (PID::hasBottom(p2->pdg_id())) {
            has_b_child = true;
            break;
          }
        }
        if (!has_b_child) {
          _theParticles.push_back(p);
          _theBs.push_back(p);
        }
      } else if (hasCharm(p)) {
        bool has_c_child = false;
        foreach (const GenParticle* p2, children) {
          if (PID::hasCharm(p2->pdg_id())) {
            has_c_child = true;
            break;
          }
        }
        if (!has_c_child) {
          _theParticles.push_back(p);
          _theCs.push_back(p);
        }
      }
    }

    MSG_DEBUG("Num b hadrons = " << _theBs.size() <<
              ", num c hadrons = " << _theCs.size() <<
              ", total = " << _theParticles.size());
  }
void setName ( const std::string &  name) [inline, inherited]

Used by derived classes to set their name.

Definition at line 120 of file Projection.hh.

References Projection::_name, and Projection::name().

Referenced by TauFinder::_init(), FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), Beam::Beam(), BeamThrust::BeamThrust(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ConstLossyFinalState::ConstLossyFinalState(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), InitialQuarks::InitialQuarks(), 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(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), UnstableFinalState::UnstableFinalState(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder().

                                        {
      _name = name;
    }

Friends And Related Function Documentation

friend class Cmp< Projection > [friend, inherited]

The Cmp specialization for Projection is a friend.

Definition at line 35 of file Projection.hh.

friend class Event [friend, inherited]

Event is a friend.

Definition at line 32 of file Projection.hh.


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

Particles _theBs [protected]

b and c hadron containers

Definition at line 80 of file HeavyHadrons.hh.

Referenced by HeavyHadrons::bHadrons(), and HeavyHadrons::project().

Particles _theCs [protected]

Definition at line 80 of file HeavyHadrons.hh.

Referenced by HeavyHadrons::cHadrons(), and HeavyHadrons::project().


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