rivet is hosted by Hepforge, IPPP Durham
Public Member Functions | Protected Member Functions | Protected Attributes | Private Attributes | Friends
Hemispheres Class Reference

Calculate the hemisphere masses and broadenings. More...

#include <Hemispheres.hh>

Inheritance diagram for Hemispheres:
Inheritance graph
[legend]
Collaboration diagram for Hemispheres:
Collaboration graph
[legend]

List of all members.

Public Member Functions

 Hemispheres (const AxesDefinition &ax)
 Constructor.
virtual const Projection * clone () const
 Clone on the heap.
void clear ()
 Reset the projection.
void calc (const Vector3 &n, const std::vector< FourMomentum > &p4s)
 Use the projection manually (i.e. outside the projection mechanism) with raw 4-momentum inputs.
void calc (const Vector3 &n, const Particles &particles)
 Use the projection manually (i.e. outside the projection mechanism) with particle inputs.
void calc (const Vector3 &n, const Jets &jets)
 Use the projection manually (i.e. outside the projection mechanism) with jet inputs.
bool massMatchesBroadening ()
 Is the hemisphere with the max mass the same as the one with the max broadening?
bool before (const Projection &p) const
virtual const std::set< PdgIdPair > beamPairs () const
virtual std::string name () const
 Get the name of the projection.
Projection & addPdgIdPair (PdgId beam1, PdgId beam2)
 Add a colliding beam pair.
Log & getLog () 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.
Hemisphere masses (scaled by \f$ 1 / E^2_\mathrm{vis} \f$).
double E2vis () const
double Evis () const
double M2high () const
double Mhigh () const
double M2low () const
double Mlow () const
double M2diff () const
double Mdiff () const
double M2sum () const
double Msum () const
double scaledM2high () const
double scaledMhigh () const
double scaledM2low () const
double scaledMlow () const
double scaledM2diff () const
double scaledMdiff () const
Hemisphere broadenings.
double Bmax () const
double Bmin () const
double Bsum () const
double Bdiff () const
Projection "getting" functions
std::set< ConstProjectionPtr > getProjections () const
 Get the contained projections, including recursion.
template<typename PROJ >
const PROJ & getProjection (const std::string &name) const
 Get the named projection, specifying return type via a template argument.
const Projection & getProjection (const std::string &name) const
Projection applying functions
template<typename PROJ >
const PROJ & applyProjection (const Event &evt, const PROJ &proj) const
 Apply the supplied projection on event.
template<typename PROJ >
const PROJ & applyProjection (const Event &evt, const Projection &proj) const
 Apply the supplied projection on event.
template<typename PROJ >
const PROJ & applyProjection (const Event &evt, const std::string &name) const
 Apply the named projection on event.

Protected Member Functions

void project (const Event &e)
 Perform the projection on the Event.
int compare (const Projection &p) const
 Compare with other projections.
Cmp< Projection > mkNamedPCmp (const Projection &otherparent, const std::string &pname) const
Cmp< Projection > mkPCmp (const Projection &otherparent, const std::string &pname) const
ProjectionHandler & getProjHandler () 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

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

Private Attributes

double _E2vis
 Visible energy-squared, $ E^2_\mathrm{vis} $.
double _M2high
 Hemisphere mass variables.
double _M2low
double _Bmax
 Hemisphere broadening variables.
double _Bmin
bool _highMassEqMaxBroad
 Is the hemisphere with the max mass the same as the one with the max broadening?

Friends

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

Detailed Description

Calculate the hemisphere masses and broadenings.

Calculate the hemisphere masses and broadenings, with event hemispheres defined by the plane normal to the thrust vector, $ \vec{n}_\mathrm{T} $.

The "high" hemisphere mass, $ M^2_\mathrm{high} / E^2_\mathrm{vis} $, is defined as

\[ \frac{M^2_\mathrm{high}}{E^2_\mathrm{vis}} = \frac{1}{E^2_\mathrm{vis}} \max \left( \left| \sum_{\vec{p}_k \cdot \vec{n}_\mathrm{T} > 0} p_k \right|^2 , \left| \sum_{\vec{p}_k \cdot \vec{n}_\mathrm{T} < 0} p_k \right|^2 \right) \]

and the corresponding "low" hemisphere mass, $ M^2_\mathrm{low} / E^2_\mathrm{vis} $, is the sum of momentum vectors in the opposite hemisphere, i.e. $ \max \rightarrow \min $ in the formula above.

Finally, we define a hemisphere mass difference:

\[ \frac{M^2_\mathrm{diff} }{ E^2_\mathrm{vis}} = \frac{ M^2_\mathrm{high} - M^2_\mathrm{low} }{ E^2_\mathrm{vis}} . \]

Similarly to the masses, we also define hemisphere broadenings, using the momenta transverse to the thrust axis:

\[ B_\pm = \frac{ \sum{\pm \vec{p}_i \cdot \vec{n}_\mathrm{T} > 0} |\vec{p}_i \times \vec{n}_\mathrm{T} | }{ 2 \sum_i | \vec{p}_i | } \]

and then a set of the broadening maximum, minimum, sum and difference as follows:

\[ B_\mathrm{max} = \max(B_+, B_-) \]

\[ B_\mathrm{min} = \min(B_+, B_-) \]

\[ B_\mathrm{sum} = B_+ + B_- \]

\[ B_\mathrm{diff} = |B_+ - B_-| \]

Internally, this projection uses a Thrust or Sphericity projection to determine the hemisphere orientation.

Definition at line 56 of file Hemispheres.hh.

Constructor & Destructor Documentation

Hemispheres ( const AxesDefinition &  ax) [inline]

Constructor.

Definition at line 60 of file Hemispheres.hh.

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

Referenced by Hemispheres::clone().

                                          {
      setName("Hemispheres");
      addProjection(ax, "Axes");
      clear();
    }

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;
  }
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 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_S8941262::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), ATLAS_2010_CONF_2010_049::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(), 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(), D0_2010_S8570965::init(), ATLAS_2011_I925932::init(), MC_JETTAGS::init(), CDF_2012_NOTE10874::init(), CMS_2013_I1273574::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_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(), ALEPH_2002_S4823664::init(), DELPHI_1999_S3960137::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(), LHCB_2011_I919315::init(), ATLAS_2012_I1119557::init(), MC_GENERIC::init(), MC_LEADJETUE::init(), MC_TTBAR::init(), CDF_2008_S7828950::init(), ATLAS_2012_I1180197::init(), D0_2001_S4674421::init(), MC_PHOTONJETUE::init(), ATLAS_2014_I1300647::init(), ATLAS_2012_I1188891::init(), D0_2008_S7837160::init(), CDF_2008_S7540469::init(), ATLAS_2011_S9131140::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(), ATLAS_2012_CONF_2012_103::init(), JADE_1998_S3612880::init(), ATLAS_2012_CONF_2012_109::init(), ATLAS_2012_I1095236::init(), ATLAS_2011_S8983313::init(), ATLAS_2012_I1126136::init(), JADE_OPAL_2000_S4300807::init(), OPAL_1997_S3396100::init(), STAR_2006_S6860818::init(), CDF_2008_S8095620::init(), CDF_1996_S3108457::init(), CMS_2012_I941555::init(), ATLAS_2013_I1217863_Z::init(), ATLAS_2011_I926145::init(), D0_1996_S3214044::init(), ATLAS_2011_CONF_2011_098::init(), MC_QCD_PARTONS::init(), EXAMPLE::init(), ATLAS_2012_I1186556::init(), ALICE_2010_S8624100::init(), ATLAS_2013_I1217863_W::init(), CMS_2014_I1303894::init(), DELPHI_1995_S3137023::init(), STAR_2009_UE_HELEN::init(), ALEPH_2004_S5765862::init(), ATLAS_2011_S9212183::init(), ATLAS_2012_I1190891::init(), ATLAS_2012_CONF_2012_001::init(), CMS_2013_I1218372::init(), MC_HFJETS::init(), ATLAS_2012_I1125575::init(), ATLAS_2012_I1112263::init(), ATLAS_2012_I1125961::init(), LHCB_2013_I1218996::init(), ALICE_2010_S8706239::init(), ATLAS_2014_I1288706::init(), ALICE_2010_S8625980::init(), CDF_2009_NOTE_9936::init(), CDF_1996_S3349578::init(), ATLAS_2013_I1263495::init(), D0_2009_S8320160::init(), CDF_2000_S4266730::init(), CDF_2001_S4563131::init(), MC_SUSY::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), CDF_2009_S8436959::init(), D0_2009_S8349509::init(), CDF_1998_S3618439::init(), CDF_2001_S4517016::init(), SLD_2002_S4869273::init(), ALEPH_2001_S4656318::init(), ATLAS_2010_S8919674::init(), D0_2006_S6438750::init(), D0_2010_S8566488::init(), DELPHI_2002_069_CONF_603::init(), ATLAS_2012_I1082936::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(), MC_WPOL::init(), CDF_2008_S7541902::init(), CDF_2006_S6653332::init(), ATLAS_2011_S9108483::init(), ATLAS_2014_I1304688::init(), D0_2008_S7662670::init(), ATLAS_2014_I1306294::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_I1083318::init(), ATLAS_2012_CONF_2012_153::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_WJET::init(), CMS_2013_I1224539_DIJET::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 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));
    }
double Bdiff ( ) const [inline]

Definition at line 148 of file Hemispheres.hh.

References Hemispheres::_Bmax, and Hemispheres::_Bmin.

Referenced by DELPHI_1996_S3430090::analyze().

{ return fabs(_Bmax - _Bmin); } // <- fabs(), just in case...
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);
    }
  }
double Bmax ( ) const [inline]

Definition at line 145 of file Hemispheres.hh.

References Hemispheres::_Bmax.

Referenced by JADE_1998_S3612880::analyze(), OPAL_2004_S6132243::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _Bmax; }
double Bmin ( ) const [inline]

Definition at line 146 of file Hemispheres.hh.

References Hemispheres::_Bmin.

Referenced by OPAL_2004_S6132243::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _Bmin; }
double Bsum ( ) const [inline]

Definition at line 147 of file Hemispheres.hh.

References Hemispheres::_Bmax, and Hemispheres::_Bmin.

Referenced by JADE_1998_S3612880::analyze(), OPAL_2004_S6132243::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _Bmax + _Bmin; }
void calc ( const Vector3 &  n,
const std::vector< FourMomentum > &  p4s 
)

Use the projection manually (i.e. outside the projection mechanism) with raw 4-momentum inputs.

Definition at line 34 of file Hemispheres.cc.

References Hemispheres::_Bmax, Hemispheres::_Bmin, Hemispheres::_E2vis, Hemispheres::_highMassEqMaxBroad, Hemispheres::_M2high, Hemispheres::_M2low, Rivet::dot(), FourMomentum::E(), Hemispheres::Evis(), Vector< N >::mod(), MSG_DEBUG, MSG_WARNING, Rivet::PID::n, Rivet::sqr(), and FourVector::vector3().

Referenced by Hemispheres::calc(), and Hemispheres::project().

                                                                          {
    MSG_DEBUG("Hemisphere axis = " << n);
    MSG_DEBUG("Number of constituents = " << p4s.size());

    FourMomentum p4With, p4Against;
    double Evis(0), broadWith(0), broadAgainst(0), broadDenom(0);
    foreach (const FourMomentum& p4, p4s) {
      const Vector3 p3 = p4.vector3();
      const double p3Para = dot(p3, n);
      const double p3Trans = (p3 - p3Para * n).mod();

      // Update normalisations
      Evis += p4.E();
      broadDenom += 2.0 * p3.mod();

      // Update the mass and broadening variables
      if (p3Para > 0) {
        p4With += p4;
        broadWith += p3Trans;
      } else if (p3Para < 0) {
        p4Against += p4;
        broadAgainst += p3Trans;
      } else {
        // In the incredibly unlikely event that a particle goes exactly along the
        // thrust plane, add half to each hemisphere.
        MSG_WARNING("Particle split between hemispheres");
        p4With += 0.5 * p4;
        p4Against += 0.5 * p4;
        broadWith += 0.5 * p3Trans;
        broadAgainst += 0.5 * p3Trans;
      }
    }

    // Visible energy squared.
    _E2vis = sqr(Evis);

    // Calculate masses.
    const double mass2With = p4With.mass2();
    const double mass2Against = p4Against.mass2();
    _M2high = max(mass2With, mass2Against);
    _M2low = min(mass2With, mass2Against);

    // Calculate broadenings.
    broadWith /= broadDenom;
    broadAgainst /= broadDenom;
    _Bmax = max(broadWith, broadAgainst);
    _Bmin = min(broadWith, broadAgainst);

    // Calculate high-max correlation flag.
    const int maxMassID = (mass2With >= mass2Against);
    const int maxBroadID = (broadWith >= broadAgainst);
    _highMassEqMaxBroad = (maxMassID == maxBroadID);
  }
void calc ( const Vector3 &  n,
const Particles &  particles 
)

Use the projection manually (i.e. outside the projection mechanism) with particle inputs.

Definition at line 20 of file Hemispheres.cc.

References Hemispheres::calc(), and ParticleBase::mom().

                                                                     {
    vector<FourMomentum> p4s; p4s.reserve(particles.size());
    foreach (const Particle& p, particles) p4s.push_back(p.mom());
    calc(n, p4s);
  }
void calc ( const Vector3 &  n,
const Jets &  jets 
)

Use the projection manually (i.e. outside the projection mechanism) with jet inputs.

Definition at line 27 of file Hemispheres.cc.

References Hemispheres::calc(), and ParticleBase::mom().

                                                           {
    vector<FourMomentum> p4s; p4s.reserve(jets.size());
    foreach (const Jet& j, jets) p4s.push_back(j.mom());
    calc(n, p4s);
  }
void clear ( ) [inline]

Reset the projection.

Definition at line 72 of file Hemispheres.hh.

References Hemispheres::_Bmax, Hemispheres::_Bmin, Hemispheres::_E2vis, Hemispheres::_highMassEqMaxBroad, Hemispheres::_M2high, and Hemispheres::_M2low.

Referenced by Hemispheres::Hemispheres(), and Hemispheres::project().

                 {
      _E2vis = -1;
      _M2high = -1;
      _M2low = -1;
      _Bmax = -1;
      _Bmin = -1;
      _highMassEqMaxBroad = true;
    }
virtual const Projection* clone ( ) const [inline, virtual]

Clone on the heap.

Implements Projection.

Definition at line 67 of file Hemispheres.hh.

References Hemispheres::Hemispheres().

                                            {
      return new Hemispheres(*this);
    }
int compare ( const Projection &  p) const [inline, protected, virtual]

Compare with other projections.

Implements Projection.

Definition at line 95 of file Hemispheres.hh.

References Projection::mkNamedPCmp().

                                           {
      return mkNamedPCmp(p, "Axes");
    }
double E2vis ( ) const [inline]

Definition at line 105 of file Hemispheres.hh.

References Hemispheres::_E2vis.

{ return _E2vis; }
double Evis ( ) const [inline]

Definition at line 106 of file Hemispheres.hh.

References Hemispheres::_E2vis.

Referenced by Hemispheres::calc().

{ return sqrt(_E2vis); }
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(), PromptFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), UnstableFinalState::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().

                                                      {
      return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
    }
ProjectionHandler& getProjHandler ( ) const [inline, protected, inherited]

Get a reference to the ProjectionHandler for this thread.

Definition at line 97 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

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

                                              {
      return _projhandler;
    }
double M2diff ( ) const [inline]

Definition at line 114 of file Hemispheres.hh.

References Hemispheres::_M2high, and Hemispheres::_M2low.

Referenced by Hemispheres::Mdiff(), and Hemispheres::scaledM2diff().

{ return _M2high -_M2low; }
double M2high ( ) const [inline]

Definition at line 108 of file Hemispheres.hh.

References Hemispheres::_M2high.

Referenced by Hemispheres::Mhigh().

{ return _M2high; }
double M2low ( ) const [inline]

Definition at line 111 of file Hemispheres.hh.

References Hemispheres::_M2low.

Referenced by Hemispheres::Mlow().

{ return _M2low; }
double M2sum ( ) const [inline]

Definition at line 117 of file Hemispheres.hh.

References Hemispheres::_M2high, and Hemispheres::_M2low.

Referenced by Hemispheres::Msum().

{ return _M2high +_M2low; }
bool massMatchesBroadening ( ) [inline]

Is the hemisphere with the max mass the same as the one with the max broadening?

Definition at line 153 of file Hemispheres.hh.

References Hemispheres::_highMassEqMaxBroad.

                                 {
      return _highMassEqMaxBroad;
    }
double Mdiff ( ) const [inline]

Definition at line 115 of file Hemispheres.hh.

References Hemispheres::M2diff().

{ return sqrt(M2diff()); }
double Mhigh ( ) const [inline]

Definition at line 109 of file Hemispheres.hh.

References Hemispheres::M2high().

{ return sqrt(M2high()); }
Cmp< Projection > mkNamedPCmp ( 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 45 of file Projection.cc.

References Rivet::pcmp().

Referenced by BeamThrust::compare(), ChargedLeptons::compare(), FParameter::compare(), CentralEtHCM::compare(), MergedFinalState::compare(), ChargedFinalState::compare(), DISLepton::compare(), DISKinematics::compare(), DISFinalState::compare(), VisibleFinalState::compare(), TauFinder::compare(), PromptFinalState::compare(), NeutralFinalState::compare(), Spherocity::compare(), ParisiTensor::compare(), LeadingParticlesFinalState::compare(), FoxWolframMoments::compare(), MissingMomentum::compare(), LossyFinalState< ConstRandomFilter >::compare(), Thrust::compare(), InvMassFinalState::compare(), Sphericity::compare(), HeavyHadrons::compare(), ZFinder::compare(), Hemispheres::compare(), DressedLeptons::compare(), WFinder::compare(), IdentifiedFinalState::compare(), VetoedFinalState::compare(), JetShape::compare(), and FastJets::compare().

                                                                     {
    return pcmp(*this, otherparent, pname);
  }
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);
  }
double Mlow ( ) const [inline]

Definition at line 112 of file Hemispheres.hh.

References Hemispheres::M2low().

{ return sqrt(M2low()); }
double Msum ( ) const [inline]

Definition at line 118 of file Hemispheres.hh.

References Hemispheres::M2sum().

{ return sqrt(M2sum()); }
virtual std::string name ( ) const [inline, virtual, inherited]

Get the name of the projection.

Implements ProjectionApplier.

Definition at line 101 of file Projection.hh.

References Projection::_name.

Referenced by ProjectionApplier::_addProjection(), ProjectionHandler::_checkDuplicate(), ProjectionHandler::_clone(), ProjectionHandler::_getEquiv(), VetoedFinalState::addVetoOnThisFinalState(), Projection::getLog(), ProjectionHandler::registerProjection(), and Projection::setName().

                                   {
      return _name;
    }
void project ( const Event &  e) [protected, virtual]

Perform the projection on the Event.

Implements Projection.

Definition at line 8 of file Hemispheres.cc.

References AxesDefinition::axis1(), Hemispheres::calc(), Hemispheres::clear(), ProjectionApplier::getProjection(), Rivet::PID::n, Rivet::particles(), and ParticleFinder::particles().

                                          {
    clear();

    // Get thrust axes.
    const AxesDefinition& ax = applyProjection<AxesDefinition>(e, "Axes");
    const Vector3 n = ax.axis1();
    const FinalState& fs = applyProjection<FinalState>(e, ax.getProjection("FS"));
    const Particles& particles = fs.particles();
    calc(n, particles);
  }
double scaledM2diff ( ) const [inline]

Definition at line 134 of file Hemispheres.hh.

References Hemispheres::_E2vis, and Hemispheres::M2diff().

Referenced by ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), and Hemispheres::scaledMdiff().

                                {
      if (M2diff() == 0.0) return 0.0;
      if (_E2vis != 0.0) return M2diff()/_E2vis;
      else return std::numeric_limits<double>::max();
    }
double scaledM2high ( ) const [inline]

Definition at line 120 of file Hemispheres.hh.

References Hemispheres::_E2vis, Hemispheres::_M2high, and Rivet::isZero().

Referenced by JADE_1998_S3612880::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), and Hemispheres::scaledMhigh().

                                {
      if (isZero(_M2high)) return 0.0;
      if (!isZero(_E2vis)) return _M2high/_E2vis;
      else return std::numeric_limits<double>::max();
    }
double scaledM2low ( ) const [inline]

Definition at line 127 of file Hemispheres.hh.

References Hemispheres::_E2vis, Hemispheres::_M2low, and Rivet::isZero().

Referenced by DELPHI_1996_S3430090::analyze(), and Hemispheres::scaledMlow().

                               {
      if (isZero(_M2low)) return 0.0;
      if (!isZero(_E2vis)) return _M2low/_E2vis;
      else return std::numeric_limits<double>::max();
    }
double scaledMdiff ( ) const [inline]

Definition at line 139 of file Hemispheres.hh.

References Hemispheres::scaledM2diff().

{ return sqrt(scaledM2diff()); }
double scaledMhigh ( ) const [inline]

Definition at line 125 of file Hemispheres.hh.

References Hemispheres::scaledM2high().

Referenced by OPAL_2004_S6132243::analyze().

{ return sqrt(scaledM2high()); }
double scaledMlow ( ) const [inline]

Definition at line 132 of file Hemispheres.hh.

References Hemispheres::scaledM2low().

Referenced by OPAL_2004_S6132243::analyze().

{ return sqrt(scaledM2low()); }
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 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(), TauFinder::TauFinder(), 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().

double _Bmax [private]

Hemisphere broadening variables.

Definition at line 167 of file Hemispheres.hh.

Referenced by Hemispheres::Bdiff(), Hemispheres::Bmax(), Hemispheres::Bsum(), Hemispheres::calc(), and Hemispheres::clear().

double _Bmin [private]

Definition at line 167 of file Hemispheres.hh.

Referenced by Hemispheres::Bdiff(), Hemispheres::Bmin(), Hemispheres::Bsum(), Hemispheres::calc(), and Hemispheres::clear().

double _E2vis [private]

Visible energy-squared, $ E^2_\mathrm{vis} $.

Definition at line 161 of file Hemispheres.hh.

Referenced by Hemispheres::calc(), Hemispheres::clear(), Hemispheres::E2vis(), Hemispheres::Evis(), Hemispheres::scaledM2diff(), Hemispheres::scaledM2high(), and Hemispheres::scaledM2low().

bool _highMassEqMaxBroad [private]

Is the hemisphere with the max mass the same as the one with the max broadening?

Definition at line 170 of file Hemispheres.hh.

Referenced by Hemispheres::calc(), Hemispheres::clear(), and Hemispheres::massMatchesBroadening().

double _M2high [private]

Hemisphere mass variables.

Definition at line 164 of file Hemispheres.hh.

Referenced by Hemispheres::calc(), Hemispheres::clear(), Hemispheres::M2diff(), Hemispheres::M2high(), Hemispheres::M2sum(), and Hemispheres::scaledM2high().

double _M2low [private]

Definition at line 164 of file Hemispheres.hh.

Referenced by Hemispheres::calc(), Hemispheres::clear(), Hemispheres::M2diff(), Hemispheres::M2low(), Hemispheres::M2sum(), and Hemispheres::scaledM2low().

The documentation for this class was generated from the following files:
Generated on Tue Mar 24 2015 17:35:44 for The Rivet MC analysis system by   doxygen 1.7.6.1