WFinder Class Reference

Convenience finder of leptonically decaying Ws. More...

#include <WFinder.hh>

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

List of all members.

Public Types

typedef Particle entity_type
typedef ParticleVector collection_type

Public Member Functions

const FinalState & remainingFinalState () const
const FinalState & constituentsFinalState () const
 Access to the W constituent leptons and photons.
const FinalState & constituentLeptonsFinalState () const
 Access to the W constituent leptons.
const FinalState & clusteredPhotonsFinalState () const
 Access to the photons which have been clustered to the leptons.
void clear ()
 Clear the projection.
virtual const ParticleVector & particles () const
 Get the final-state particles.
template<typename F >
const ParticleVector & particles (F sorter) const
 Get the final-state particles, ordered by supplied sorting function object.
const ParticleVector & particlesByPt () const
 Get the final-state particles, ordered by decreasing $ p_T $.
const ParticleVector & particlesByE () const
 Get the final-state particles, ordered by decreasing $ E $.
const ParticleVector & particlesByEt () const
 Get the final-state particles, ordered by decreasing $ E_T $.
const ParticleVector & particlesByEta () const
 Get the final-state particles, ordered by increasing $ \eta $.
const ParticleVector & particlesByModEta () const
 Get the final-state particles, ordered by increasing $ |\eta| $.
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 size_t ptMin () const
 Minimum-$ p_\perp $ requirement.
const collection_type & entities () const
 Template-usable interface common to JetAlg.
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.
Constructors

 WFinder (const ChargedFinalState &fs_l, PdgId pid, double m2_min, double m2_max, double missingET, double dRmax)
 WFinder (double etaMin, double etaMax, double pTmin, PdgId pid, double m2_min, double m2_max, double missingET, double dRmax)
 WFinder (const std::vector< std::pair< double, double > > &etaRanges, double pTmin, PdgId pid, double m2_min, const double m2_max, double missingET, double dRmax)
virtual const Projection * clone () const
 Clone on the heap.
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)
 Apply the projection on the supplied event.
int compare (const Projection &p) const
 Compare projections.
bool accept (const Particle &p) const
 Decide if a particle is to be accepted or not.
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

vector< pair< double, double > > _etaRanges
 The ranges allowed for pseudorapidity.
double _ptmin
 The minimum allowed transverse momentum.
ParticleVector _theParticles
 The final-state particles.
bool _allowProjReg
 Flag to forbid projection registration in analyses until the init phase.

Private Member Functions

void _init (const std::vector< std::pair< double, double > > &etaRanges, double pTmin, PdgId pid, double m2_min, double m2_max, double missingET, double dRmax)
 Common implementation of constructor operation, taking FS params.
void _init (const ChargedFinalState &fs_l, PdgId pid, double m2_min, double m2_max, double missingET, double dRmax)
 Common implementation of constructor operation, taking FS.

Private Attributes

double _m2_min
double _m2_max
double _etMiss

Friends

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

Detailed Description

Convenience finder of leptonically decaying Ws.

Chain together different projections as convenience for finding W's from two leptons in the final state, including photon clustering.

Definition at line 19 of file WFinder.hh.

Member Typedef Documentation

typedef ParticleVector collection_type [inherited]

Definition at line 88 of file FinalState.hh.

typedef Particle entity_type [inherited]

Definition at line 87 of file FinalState.hh.

Constructor & Destructor Documentation

WFinder ( const ChargedFinalState &  fs_l,
PdgId  pid,
double  m2_min,
double  m2_max,
double  missingET,
double  dRmax 
)

Constructor taking a FinalState and type of the charged lepton, mass window, and maximum dR of photons around the charged lepton to take into account for W reconstruction.

Definition at line 15 of file WFinder.cc.

References WFinder::_init().

Referenced by WFinder::clone().

00019                                  {
00020     _init(fs_l, pid, m2_min, m2_max, missingET, dRmax);
00021   }

WFinder ( double  etaMin,
double  etaMax,
double  pTmin,
PdgId  pid,
double  m2_min,
double  m2_max,
double  missingET,
double  dRmax 
)

Constructor taking single eta/pT bounds and type of the charged lepton, mass window, and maximum dR of photons around the charged lepton to take into account for W reconstruction.

Definition at line 24 of file WFinder.cc.

References WFinder::_init().

00029                                  {
00030     vector<pair<double, double> > etaRanges;
00031     etaRanges += std::make_pair(etaMin, etaMax);
00032     _init(etaRanges, pTmin, pid, m2_min, m2_max, missingET, dRmax);
00033   }

WFinder ( const std::vector< std::pair< double, double > > &  etaRanges,
double  pTmin,
PdgId  pid,
double  m2_min,
const double  m2_max,
double  missingET,
double  dRmax 
)

Constructor taking multiple eta/pT bounds and type of the charged lepton, mass window, and maximum dR of photons around the charged lepton to take into account for W reconstruction.

Definition at line 36 of file WFinder.cc.

References WFinder::_init().

00041                                  {
00042     _init(etaRanges, pTmin, pid, m2_min, m2_max, missingET, dRmax);
00043   }

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, Log::ERROR, ProjectionApplier::getLog(), ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().

Referenced by ProjectionApplier::addProjection().

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

void _init ( const ChargedFinalState &  fs_l,
PdgId  pid,
double  m2_min,
double  m2_max,
double  missingET,
double  dRmax 
) [private]

Common implementation of constructor operation, taking FS.

Definition at line 57 of file WFinder.cc.

References WFinder::_etMiss, WFinder::_m2_max, WFinder::_m2_min, IdentifiedFinalState::acceptNeutrinos(), ProjectionApplier::addProjection(), VetoedFinalState::addVetoOnThisFinalState(), Rivet::ELECTRON, FinalState::FinalState(), Rivet::MUON, Rivet::NU_E, Rivet::NU_MU, and Projection::setName().

00062   {
00063     setName("WFinder");
00064 
00065     // Check that the arguments are legal
00066     assert(abs(pid) == ELECTRON || abs(pid) == MUON);
00067     PdgId nu_pid = abs(pid) + 1;
00068     assert(abs(nu_pid) == NU_E || abs(nu_pid) == NU_MU);
00069 
00070     // Don't make pT or eta cuts on the neutrino
00071     IdentifiedFinalState fs_nu;
00072     fs_nu.acceptNeutrinos();
00073 
00074     // Make a merged final state projection for charged and neutral leptons
00075     MergedFinalState mergedFS(fs_l, fs_nu);
00076 
00077     // Mass range
00078     _m2_min = m2_min;
00079     _m2_max = m2_max;
00080 
00081     // Set ETmiss
00082     _etMiss = missingET;
00083 
00084     // Make and register an invariant mass final state for the W decay leptons
00085     vector<pair<PdgId, PdgId> > l_nu_ids;
00086     l_nu_ids += make_pair(abs(pid), -abs(nu_pid));
00087     l_nu_ids += make_pair(-abs(pid), abs(nu_pid));
00088     InvMassFinalState imfs(mergedFS, l_nu_ids, m2_min, m2_max);
00089     addProjection(imfs, "IMFS");
00090 
00091     // A projection for clustering photons on to the charged lepton
00092     ClusteredPhotons cphotons(FinalState(), imfs, dRmax);
00093     addProjection(cphotons, "CPhotons");
00094 
00095     // Projection for all signal constituents
00096     MergedFinalState signalFS(imfs, cphotons);
00097     addProjection(signalFS, "SignalParticles");
00098 
00099     // Add MissingMomentum proj to calc MET
00100     MissingMomentum vismom(signalFS);
00101     addProjection(vismom, "MissingET");
00102 
00103     // FS for non-signal bits of the event
00104     VetoedFinalState remainingFS;
00105     remainingFS.addVetoOnThisFinalState(signalFS);
00106     addProjection(remainingFS, "RFS");
00107   }

void _init ( const std::vector< std::pair< double, double > > &  etaRanges,
double  pTmin,
PdgId  pid,
double  m2_min,
double  m2_max,
double  missingET,
double  dRmax 
) [private]

Common implementation of constructor operation, taking FS params.

Definition at line 46 of file WFinder.cc.

Referenced by WFinder::WFinder().

00051                                     {
00052     ChargedFinalState fs_l(etaRanges, pTmin);
00053     _init(fs_l, pid, m2_min, m2_max, missingET, dRmax);
00054   }

bool accept ( const Particle &  p  )  const [protected, inherited]

Decide if a particle is to be accepted or not.

Definition at line 99 of file FinalState.cc.

References FinalState::_etaRanges, FinalState::_ptmin, FourVector::eta(), Rivet::eta(), Particle::genParticle(), Particle::momentum(), FourMomentum::pT(), and Rivet::pT().

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

00099                                                  {
00100     // Not having s.c. == 1 should never happen!
00101     assert(p.genParticle().status() == 1);
00102 
00103     // Check pT cut
00104     if (_ptmin > 0.0) {
00105       const double pT = p.momentum().pT();
00106       if (pT < _ptmin) return false;
00107     }
00108 
00109     // Check eta cuts
00110     if (!_etaRanges.empty()) {
00111       bool eta_pass = false;
00112       const double eta = p.momentum().eta();
00113       typedef pair<double,double> EtaPair;
00114       foreach (const EtaPair& etacuts, _etaRanges) {
00115         if (eta > etacuts.first && eta < etacuts.second) {
00116           eta_pass = true;
00117           break;
00118         }
00119       }
00120       if (!eta_pass) return false;
00121     }
00122  
00123     return true;
00124   }

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

00107                                                        {
00108       _beamPairs.insert(PdgIdPair(beam1, beam2));
00109       return *this;
00110     }

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 130 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

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

00130                                                                        {
00131       const Projection& reg = _addProjection(proj, name);
00132       const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
00133       return rtn;
00134     }

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

Apply the named projection on event.

Definition at line 97 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

00097                                                                                {
00098       return pcast<PROJ>(_applyProjection(evt, name));
00099     }

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

Apply the supplied projection on event.

Definition at line 90 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

00090                                                                                 {
00091       return pcast<PROJ>(_applyProjection(evt, proj));
00092     }

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

Apply the supplied projection on event.

Definition at line 83 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

00083                                                                           {
00084       return pcast<PROJ>(_applyProjection(evt, proj));
00085     }

const set< PdgIdPair > beamPairs (  )  const [virtual, inherited]

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

Definition at line 39 of file Projection.cc.

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

Referenced by Projection::beamPairs().

00039                                                    {
00040     set<PdgIdPair> ret = _beamPairs;
00041     set<ConstProjectionPtr> projs = getProjections();
00042     for (set<ConstProjectionPtr>::const_iterator ip = projs.begin(); ip != projs.end(); ++ip) {
00043       ConstProjectionPtr p = *ip;
00044       getLog() << Log::TRACE << "Proj addr = " << p << endl;
00045       if (p) ret = intersection(ret, p->beamPairs());
00046     }
00047     return ret;
00048   }

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 28 of file Projection.cc.

References Projection::compare().

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

00028                                                    {
00029     const std::type_info& thisid = typeid(*this);
00030     const std::type_info& otherid = typeid(p);
00031     if (thisid == otherid) {
00032       return compare(p) < 0;
00033     } else {
00034       return thisid.before(otherid);
00035     }
00036   }

void clear (  ) 

Clear the projection.

Definition at line 145 of file WFinder.cc.

References FinalState::_theParticles.

Referenced by WFinder::project().

00145                       {
00146     _theParticles.clear();
00147   }

virtual const Projection* clone (  )  const [inline, virtual]

Clone on the heap.

Reimplemented from FinalState.

Definition at line 58 of file WFinder.hh.

References WFinder::WFinder().

00058                                             {
00059       return new WFinder(*this);
00060     }

const FinalState & clusteredPhotonsFinalState (  )  const

Access to the photons which have been clustered to the leptons.

Definition at line 128 of file WFinder.cc.

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

00129   {
00130     return getProjection<FinalState>("CPhotons");
00131   }

int compare ( const Projection &  p  )  const [protected, virtual]

Compare projections.

Reimplemented from FinalState.

Definition at line 134 of file WFinder.cc.

References Rivet::cmp(), Rivet::EQUIVALENT, and Projection::mkNamedPCmp().

00134                                                 {
00135     PCmp cmp = mkNamedPCmp(p, "IMFS");
00136     if (cmp != EQUIVALENT) return cmp;
00137 
00138     cmp = mkNamedPCmp(p, "CPhotons");
00139     if (cmp != EQUIVALENT) return cmp;
00140 
00141     return EQUIVALENT;
00142   }

const FinalState & constituentLeptonsFinalState (  )  const

Access to the W constituent leptons.

Definition at line 123 of file WFinder.cc.

Referenced by MC_WJETS::analyze(), and D0_2008_S7837160::analyze().

00123                                                                 {
00124     return getProjection<FinalState>("IMFS");
00125   }

const FinalState & constituentsFinalState (  )  const

Access to the W constituent leptons and photons.

Definition at line 118 of file WFinder.cc.

Referenced by MC_WWJETS::analyze(), D0_1998_S3711838::analyze(), CDF_1991_S2313472::analyze(), and MC_WWJETS::init().

00118                                                           {
00119     return getProjection<FinalState>("SignalParticles");
00120   }

virtual bool empty (  )  const [inline, virtual, inherited]

Is this final state empty?

Definition at line 77 of file FinalState.hh.

References FinalState::_theParticles.

Referenced by STAR_2008_S7993412::analyze(), STAR_2006_S6870392::analyze(), MC_PHOTONJETS::analyze(), D0_2008_S7662670::analyze(), D0_2008_S6879055::analyze(), D0_2001_S4674421::analyze(), and CDF_2008_S7540469::analyze().

00077 { return _theParticles.empty(); }

const collection_type& entities (  )  const [inline, inherited]

Template-usable interface common to JetAlg.

Definition at line 91 of file FinalState.hh.

References FinalState::particles().

00091                                             {
00092       return particles();
00093     }

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 Sphericity::_calcSphericity(), Thrust::_calcThrust(), Projection::beamPairs(), FastJets::calc(), FastJets::FastJets(), FinalState::FinalState(), JetAlg::JetAlg(), ZFinder::project(), WFinder::project(), VetoedFinalState::project(), UnstableFinalState::project(), TriggerUA5::project(), TriggerCDFRun2::project(), TriggerCDFRun0Run1::project(), PVertex::project(), NeutralFinalState::project(), MergedFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), LeadingParticlesFinalState::project(), IsolationProjection< PROJ1, PROJ2, EST >::project(), InvMassFinalState::project(), InitialQuarks::project(), Hemispheres::project(), HadronicFinalState::project(), FoxWolframMoments::project(), ClusteredPhotons::project(), ChargedFinalState::project(), FastJets::splitJet(), and TotalVisibleMomentum::TotalVisibleMomentum().

00114                         {
00115       string logname = "Rivet.Projection." + name();
00116       return Log::getLog(logname);
00117     }

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 73 of file ProjectionApplier.hh.

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

00073                                                                  {
00074       return getProjHandler().getProjection(*this, name);
00075     }

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

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

Definition at line 65 of file ProjectionApplier.hh.

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

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

00065                                                            {
00066       const Projection& p = getProjHandler().getProjection(*this, name);
00067       return pcast<PROJ>(p);
00068     }

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

Get the contained projections, including recursion.

Definition at line 58 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

00058                                                       {
00059       return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
00060     }

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

Get a reference to the ProjectionHandler for this thread.

Definition at line 110 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

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

00110                                               {
00111       assert(_projhandler);
00112       return *_projhandler;
00113     }

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

Definition at line 79 of file FinalState.hh.

References FinalState::_theParticles.

00079 { return _theParticles.empty(); }

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 51 of file Projection.cc.

References Rivet::pcmp().

Referenced by ZFinder::compare(), WFinder::compare(), VetoedFinalState::compare(), TotalVisibleMomentum::compare(), Thrust::compare(), SVertex::compare(), Sphericity::compare(), ParisiTensor::compare(), NeutralFinalState::compare(), Multiplicity::compare(), MergedFinalState::compare(), LossyFinalState< ConstRandomFilter >::compare(), LeadingParticlesFinalState::compare(), JetShape::compare(), IsolationProjection< PROJ1, PROJ2, EST >::compare(), InvMassFinalState::compare(), IdentifiedFinalState::compare(), Hemispheres::compare(), FoxWolframMoments::compare(), FinalStateHCM::compare(), FastJets::compare(), DISLepton::compare(), DISKinematics::compare(), ClusteredPhotons::compare(), ClosestJetShape::compare(), ChargedLeptons::compare(), ChargedFinalState::compare(), and CentralEtHCM::compare().

00052                                                                      {
00053     return pcmp(*this, otherparent, pname);
00054   }

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 57 of file Projection.cc.

References Rivet::pcmp().

00058                                                                 {
00059     return pcmp(*this, otherparent, pname);
00060   }

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(), ZFinder::project(), ClusteredPhotons::project(), and ProjectionHandler::registerProjection().

00101                                    {
00102       return _name;
00103     }

const ParticleVector& particles ( F  sorter  )  const [inline, inherited]

Get the final-state particles, ordered by supplied sorting function object.

Definition at line 43 of file FinalState.hh.

References FinalState::_theParticles.

00043                                                     {
00044       std::sort(_theParticles.begin(), _theParticles.end(), sorter);
00045       return _theParticles;
00046     }

virtual const ParticleVector& particles (  )  const [inline, virtual, inherited]

Get the final-state particles.

Definition at line 39 of file FinalState.hh.

References FinalState::_theParticles.

Referenced by UA5_1986_S1583476::analyze(), UA5_1982_S875503::analyze(), UA1_1990_S2044935::analyze(), STAR_2009_UE_HELEN::analyze(), STAR_2008_S7993412::analyze(), STAR_2008_S7869363::analyze(), STAR_2006_S6500200::analyze(), SFM_1984_S1178091::analyze(), PDG_HADRON_MULTIPLICITIES_RATIOS::analyze(), PDG_HADRON_MULTIPLICITIES::analyze(), OPAL_1998_S3780481::analyze(), OPAL_1993_S2692198::analyze(), MC_ZZJETS::analyze(), MC_ZJETS::analyze(), MC_WWJETS::analyze(), MC_WJETS::analyze(), MC_TTBAR::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_PHOTONJETS::analyze(), MC_LEADINGJETS::analyze(), MC_HJETS::analyze(), MC_GENERIC::analyze(), MC_DIJET::analyze(), H1_2000_S4129130::analyze(), H1_1995_S3167097::analyze(), H1_1994_S2919893::analyze(), E735_1998_S3905616::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1995_S3137023::analyze(), D0_2010_S8671338::analyze(), D0_2009_S8349509::analyze(), D0_2009_S8202443::analyze(), D0_2008_S7863608::analyze(), D0_2008_S7837160::analyze(), D0_2008_S7719523::analyze(), D0_2008_S7554427::analyze(), D0_2008_S6879055::analyze(), D0_2007_S7075677::analyze(), D0_2006_S6438750::analyze(), D0_2001_S4674421::analyze(), D0_2000_S4480767::analyze(), D0_1998_S3711838::analyze(), CDF_2009_S8383952::analyze(), CDF_2009_S8233977::analyze(), CDF_2008_S8095620::analyze(), CDF_2008_S7541902::analyze(), CDF_2008_S7540469::analyze(), CDF_2008_NOTE_9351::analyze(), CDF_2008_LEADINGJETS::analyze(), CDF_2006_S6653332::analyze(), CDF_2002_S4796047::analyze(), CDF_2000_S4155203::analyze(), CDF_1991_S2313472::analyze(), CDF_1990_S2089246::analyze(), CDF_1988_S1865951::analyze(), ATLAS_2010_S8591806::analyze(), ATLAS_2010_CONF_2010_031::analyze(), ALICE_2010_S8706239::analyze(), ALICE_2010_S8625980::analyze(), ALEPH_2004_S5765862::analyze(), ALEPH_1996_S3486095::analyze(), Thrust::calc(), Sphericity::calc(), FinalState::entities(), FinalState::particlesByE(), FinalState::particlesByEt(), FinalState::particlesByEta(), FinalState::particlesByModEta(), FinalState::particlesByPt(), ZFinder::project(), WFinder::project(), VetoedFinalState::project(), TriggerUA5::project(), TriggerCDFRun2::project(), TriggerCDFRun0Run1::project(), TotalVisibleMomentum::project(), SVertex::project(), NeutralFinalState::project(), Multiplicity::project(), MergedFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), LeadingParticlesFinalState::project(), InvMassFinalState::project(), IdentifiedFinalState::project(), Hemispheres::project(), HadronicFinalState::project(), FoxWolframMoments::project(), FinalStateHCM::project(), FinalState::project(), DISLepton::project(), ClusteredPhotons::project(), ClosestJetShape::project(), ChargedLeptons::project(), ChargedFinalState::project(), and CentralEtHCM::project().

00039 { return _theParticles; }

const ParticleVector& particlesByE (  )  const [inline, inherited]

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

Definition at line 54 of file FinalState.hh.

References Rivet::cmpParticleByE(), and FinalState::particles().

00054                                                {
00055       return particles(cmpParticleByE);
00056     }

const ParticleVector& particlesByEt (  )  const [inline, inherited]

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

Definition at line 59 of file FinalState.hh.

References Rivet::cmpParticleByEt(), and FinalState::particles().

00059                                                 {
00060       return particles(cmpParticleByEt);
00061     }

const ParticleVector& particlesByEta (  )  const [inline, inherited]

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

Definition at line 64 of file FinalState.hh.

References Rivet::cmpParticleByAscPseudorapidity(), and FinalState::particles().

00064                                                  {
00065       return particles(cmpParticleByAscPseudorapidity);
00066     }

const ParticleVector& particlesByModEta (  )  const [inline, inherited]

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

Definition at line 69 of file FinalState.hh.

References Rivet::cmpParticleByAscAbsPseudorapidity(), and FinalState::particles().

00069                                                     {
00070       return particles(cmpParticleByAscAbsPseudorapidity);
00071     }

const ParticleVector& particlesByPt (  )  const [inline, inherited]

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

Definition at line 49 of file FinalState.hh.

References Rivet::cmpParticleByPt(), and FinalState::particles().

Referenced by MC_PHOTONJETUE::analyze(), and ATLAS_2010_CONF_2010_081::analyze().

00049                                                 {
00050       return particles(cmpParticleByPt);
00051     }

void project ( const Event &  e  )  [protected, virtual]

Apply the projection on the supplied event.

Todo:
Restrict missing momentum eta range? Use vectorET()?

Reimplemented from FinalState.

Definition at line 150 of file WFinder.cc.

References WFinder::_etMiss, WFinder::_m2_max, WFinder::_m2_min, FinalState::_theParticles, WFinder::clear(), Log::DEBUG, Projection::getLog(), Rivet::GeV, Rivet::inRange(), FourMomentum::mass(), Particle::momentum(), InvMassFinalState::particlePairs(), FinalState::particles(), Particle::pdgId(), Rivet::PID::threeCharge(), Rivet::WMINUSBOSON, and Rivet::WPLUSBOSON.

00150                                       {
00151     clear();
00152     const InvMassFinalState& imfs = applyProjection<InvMassFinalState>(e, "IMFS");
00153 
00154     Particle p1,p2;
00155     if(imfs.particles().size() == 2) {
00156       p1 = imfs.particles()[0];
00157       p2 = imfs.particles()[1];
00158     }
00159     else {
00160       // no candiate, return
00161       if(imfs.particles().empty()) {
00162     getLog() << Log::DEBUG << "No W+- candidates found" << " "
00163          << imfs.particles().size() << " " << endl;
00164     return;
00165       }
00166       // more than one, pick the one nearer the W mass
00167       double deltaM = 1e30;
00168       for(unsigned int ix=0;ix<imfs.particlePairs().size();++ix) {
00169     FourMomentum pW = imfs.particlePairs()[ix].first .momentum()+
00170       imfs.particlePairs()[ix].second.momentum();
00171     double mW = pW.mass();
00172     if(abs(mW-80.403)<deltaM) {
00173       deltaM = abs(mW-80.4);
00174       p1 = imfs.particlePairs()[ix].first ;
00175       p2 = imfs.particlePairs()[ix].second;
00176     }
00177       }
00178     }
00179 
00180     FourMomentum pW = p1.momentum() + p2.momentum();
00181     const int w3charge = PID::threeCharge(p1) + PID::threeCharge(p2);
00182     assert(abs(w3charge) == 3);
00183     const int wcharge = w3charge/3;
00184 
00185     stringstream msg;
00186     string wsign = (wcharge == 1) ? "+" : "-";
00187     string wstr = "W" + wsign;
00188     msg << wstr << " reconstructed from: " << endl
00189         << "   " << p1.momentum() << " " << p1.pdgId() << endl
00190         << " + " << p2.momentum() << " " << p2.pdgId() << endl;
00191 
00192     // Add in clustered photons
00193     const FinalState& photons = applyProjection<FinalState>(e, "CPhotons");
00194     foreach (const Particle& photon, photons.particles()) {
00195       msg << " + " << photon.momentum() << " " << photon.pdgId() << endl;
00196       pW += photon.momentum();
00197     }
00198     msg << " = " << pW;
00199 
00200     // Check missing ET
00201     const MissingMomentum& vismom = applyProjection<MissingMomentum>(e, "MissingET");
00202     /// @todo Restrict missing momentum eta range? Use vectorET()?
00203     if (vismom.scalarET() < _etMiss) {
00204       getLog() << Log::DEBUG << "Not enough missing ET: " << vismom.scalarET()/GeV
00205                << " GeV vs. " << _etMiss/GeV << " GeV" << endl;
00206       return;
00207     }
00208 
00209     // Check mass range again
00210     if (!inRange(pW.mass()/GeV, _m2_min, _m2_max)) return;
00211     getLog() << Log::DEBUG << msg.str() << endl;
00212 
00213     // Make W Particle and insert into particles list
00214     const PdgId wpid = (wcharge == 1) ? WPLUSBOSON : WMINUSBOSON;
00215     _theParticles.push_back(Particle(wpid, pW));
00216   }

virtual size_t ptMin (  )  const [inline, virtual, inherited]

Minimum-$ p_\perp $ requirement.

Definition at line 82 of file FinalState.hh.

References FinalState::_ptmin.

00082 { return _ptmin; }

const FinalState & remainingFinalState (  )  const

Access to the remaining particles, after the W and clustered photons have been removed from the full final state (e.g. for running a jet finder on it)

Definition at line 113 of file WFinder.cc.

Referenced by MC_WJETS::init().

00113                                                        {
00114     return getProjection<FinalState>("RFS");
00115   }

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.

Referenced by ZFinder::_init(), WFinder::_init(), Beam::Beam(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ClosestJetShape::ClosestJetShape(), ClusteredPhotons::ClusteredPhotons(), ConstLossyFinalState::ConstLossyFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FastJets::FastJets(), FinalState::FinalState(), FinalStateHCM::FinalStateHCM(), FoxWolframMoments::FoxWolframMoments(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), InitialQuarks::InitialQuarks(), IsolationProjection< PROJ1, PROJ2, EST >::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), PVertex::PVertex(), Sphericity::Sphericity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), UnstableFinalState::UnstableFinalState(), and VetoedFinalState::VetoedFinalState().

00120                                         {
00121       _name = name;
00122     }

virtual size_t size (  )  const [inline, virtual, inherited]

Access the projected final-state particles.

Definition at line 74 of file FinalState.hh.

References FinalState::_theParticles.

Referenced by UA5_1986_S1583476::analyze(), UA5_1982_S875503::analyze(), UA1_1990_S2044935::analyze(), STAR_2006_S6860818::analyze(), STAR_2006_S6500200::analyze(), OPAL_2004_S6132243::analyze(), MC_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_GENERIC::analyze(), D0_2008_S7837160::analyze(), D0_2000_S4480767::analyze(), D0_1998_S3711838::analyze(), CDF_2009_S8233977::analyze(), CDF_2009_NOTE_9936::analyze(), CDF_1991_S2313472::analyze(), ATLAS_2010_S8591806::analyze(), ATLAS_2010_CONF_2010_081::analyze(), ATLAS_2010_CONF_2010_031::analyze(), ALICE_2010_S8625980::analyze(), and ALICE_2010_S8624100::analyze().

00074 { return _theParticles.size(); }

Friends And Related Function Documentation

friend class Cmp< Projection > [friend, inherited]

The Cmp specialization for Projection is a friend.

Definition at line 36 of file Projection.hh.

friend class Event [friend, inherited]

Event is a friend.

Definition at line 33 of file Projection.hh.

friend class Projectionhandler [friend, inherited]

Definition at line 38 of file ProjectionApplier.hh.

Member Data Documentation

bool _allowProjReg [protected, inherited]

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

Definition at line 157 of file ProjectionApplier.hh.

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

vector<pair<double,double> > _etaRanges [protected, inherited]

The ranges allowed for pseudorapidity.

Definition at line 111 of file FinalState.hh.

Referenced by FinalState::accept(), FinalState::compare(), FinalState::FinalState(), and FinalState::project().

double _etMiss [private]

Definition at line 117 of file WFinder.hh.

Referenced by WFinder::_init(), and WFinder::project().

double _m2_max [private]

Definition at line 114 of file WFinder.hh.

Referenced by WFinder::_init(), and WFinder::project().

double _m2_min [private]

Definition at line 114 of file WFinder.hh.

Referenced by WFinder::_init(), and WFinder::project().

double _ptmin [protected, inherited]

The minimum allowed transverse momentum.

Definition at line 114 of file FinalState.hh.

Referenced by FinalState::accept(), FinalState::compare(), FinalState::project(), and FinalState::ptMin().

ParticleVector _theParticles [mutable, protected, inherited]

The final-state particles.

Definition at line 117 of file FinalState.hh.

Referenced by WFinder::clear(), FinalState::empty(), FinalState::isEmpty(), FinalState::particles(), ZFinder::project(), WFinder::project(), VetoedFinalState::project(), NeutralFinalState::project(), MergedFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), LeadingParticlesFinalState::project(), InvMassFinalState::project(), IdentifiedFinalState::project(), HadronicFinalState::project(), FinalStateHCM::project(), FinalState::project(), ClusteredPhotons::project(), ChargedFinalState::project(), and FinalState::size().

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