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 ParticleVectorbosons () const
 Access to the found bosons (currently either 0 or 1).
const vector< Particle > & constituentLeptons () const
const vector< Particle > & constituentNeutrinos () const
const FinalStateremainingFinalState () const
void clear ()
 Clear the projection.
virtual const ParticleVectorparticles () const
 Get the final-state particles.
template<typename F >
const ParticleVectorparticles (F sorter) const
 Get the final-state particles, ordered by supplied sorting function object.
const ParticleVectorparticlesByPt () const
 Get the final-state particles, ordered by decreasing $ p_T $.
const ParticleVectorparticlesByP () const
 Get the final-state particles, ordered by decreasing $ p $.
const ParticleVectorparticlesByE () const
 Get the final-state particles, ordered by decreasing $ E $.
const ParticleVectorparticlesByEt () const
 Get the final-state particles, ordered by decreasing $ E_T $.
const ParticleVectorparticlesByEta () const
 Get the final-state particles, ordered by increasing $ \eta $.
const ParticleVectorparticlesByModEta () 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 double ptMin () const
 Minimum-$ p_\perp $ requirement.
const collection_typeentities () const
 Template-usable interface common to JetAlg.
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

 WFinder (double etaMin, double etaMax, double pTmin, PdgId pid, double minmass, double maxmass, double missingET, double dRmax, bool clusterPhotons=true, bool trackPhotons=false, double masstarget=80.4, bool useTransverseMass=false)
 WFinder (const std::vector< std::pair< double, double > > &etaRanges, double pTmin, PdgId pid, double minmass, const double maxmass, double missingET, double dRmax, bool clusterPhotons=true, bool trackPhotons=false, double masstarget=80.4, bool useTransverseMass=false)
virtual const Projectionclone () const
 Clone on the heap.
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

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< 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

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 minmass, double maxmass, double missingET, double dRmax, bool clusterPhotons, bool trackPhotons, double masstarget, bool useTransverseMass)
 Common implementation of constructor operation, taking FS params.

Private Attributes

double _minmass
 Transverse mass cuts.
double _maxmass
double _masstarget
bool _useTransverseMass
double _etMiss
 Missing ET cut.
bool _trackPhotons
PdgId _pid
 Lepton flavour.
PdgId _nu_pid
 Neutrino flavour.
ParticleVector _bosons
 list of found bosons (currently either 0 or 1)
ParticleVector _constituentLeptons
 Constituent leptons (currently either 0 or 1).
ParticleVector _constituentNeutrinos
 Constituent neutrinos (currently either 0 or 1).

Friends

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

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 20 of file WFinder.hh.


Member Typedef Documentation

typedef ParticleVector collection_type [inherited]

Definition at line 93 of file FinalState.hh.

typedef Particle entity_type [inherited]

Definition at line 92 of file FinalState.hh.


Constructor & Destructor Documentation

WFinder ( double  etaMin,
double  etaMax,
double  pTmin,
PdgId  pid,
double  minmass,
double  maxmass,
double  missingET,
double  dRmax,
bool  clusterPhotons = true,
bool  trackPhotons = false,
double  masstarget = 80.4,
bool  useTransverseMass = false 
)

Constructor taking single eta/pT bounds and type of the charged lepton, transverse 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().

00022                                            {
00023     vector<pair<double, double> > etaRanges;
00024     etaRanges += std::make_pair(etaMin, etaMax);
00025     _init(etaRanges, pTmin, pid, minmass, maxmass, missingET,
00026           dRmax, clusterPhotons, trackPhotons, masstarget, useTransverseMass);
00027   }

WFinder ( const std::vector< std::pair< double, double > > &  etaRanges,
double  pTmin,
PdgId  pid,
double  minmass,
const double  maxmass,
double  missingET,
double  dRmax,
bool  clusterPhotons = true,
bool  trackPhotons = false,
double  masstarget = 80.4,
bool  useTransverseMass = false 
)

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

Definition at line 30 of file WFinder.cc.

References WFinder::_init().

00037                                            {
00038     _init(etaRanges, pTmin, pid, minmass, maxmass, missingET,
00039           dRmax, clusterPhotons, trackPhotons, masstarget, useTransverseMass);
00040   }


Member Function Documentation

const Projection & _addProjection ( const Projection proj,
const std::string &  name 
) [protected, inherited]

Untemplated function to do the work...

Definition at line 33 of file ProjectionApplier.cc.

References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().

Referenced by ProjectionApplier::addProjection().

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

void _init ( const std::vector< std::pair< double, double > > &  etaRanges,
double  pTmin,
PdgId  pid,
double  minmass,
double  maxmass,
double  missingET,
double  dRmax,
bool  clusterPhotons,
bool  trackPhotons,
double  masstarget,
bool  useTransverseMass 
) [private]

Common implementation of constructor operation, taking FS params.

Definition at line 43 of file WFinder.cc.

References WFinder::_etMiss, WFinder::_masstarget, WFinder::_maxmass, WFinder::_minmass, WFinder::_nu_pid, WFinder::_pid, WFinder::_trackPhotons, WFinder::_useTransverseMass, IdentifiedFinalState::acceptIdPair(), IdentifiedFinalState::acceptNeutrinos(), ProjectionApplier::addProjection(), VetoedFinalState::addVetoOnThisFinalState(), Rivet::ELECTRON, Rivet::MUON, Rivet::NU_E, Rivet::NU_MU, and Projection::setName().

Referenced by WFinder::WFinder().

00051   {
00052     setName("WFinder");
00053 
00054     _minmass = minmass;
00055     _maxmass = maxmass;
00056     _masstarget = masstarget;
00057     _pid = pid;
00058     _trackPhotons = trackPhotons;
00059     _useTransverseMass = useTransverseMass;
00060 
00061     // Check that the arguments are legal
00062     assert(abs(_pid) == ELECTRON || abs(_pid) == MUON);
00063     _nu_pid = abs(_pid) + 1;
00064     assert(abs(_nu_pid) == NU_E || abs(_nu_pid) == NU_MU);
00065 
00066     // Don't make pT or eta cuts on the neutrino
00067     IdentifiedFinalState neutrinos;
00068     neutrinos.acceptNeutrinos();
00069     addProjection(neutrinos, "Neutrinos");
00070 
00071     // Lepton clusters
00072     FinalState fs;
00073     IdentifiedFinalState bareleptons(fs);
00074     bareleptons.acceptIdPair(pid);
00075     LeptonClusters leptons(fs, bareleptons, dRmax,
00076                            clusterPhotons,
00077                            etaRanges, pTmin);
00078     addProjection(leptons, "LeptonClusters");
00079 
00080     // Add MissingMomentum proj to calc MET
00081     MissingMomentum vismom(fs);
00082     addProjection(vismom, "MissingET");
00083     // Set ETmiss
00084     _etMiss = missingET;
00085 
00086     VetoedFinalState remainingFS;
00087     remainingFS.addVetoOnThisFinalState(*this);
00088     addProjection(remainingFS, "RFS");
00089   }

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::hasGenParticle(), Particle::momentum(), FourMomentum::pT(), and Rivet::pT().

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

00099                                                  {
00100     // Not having s.c. == 1 should never happen!
00101     assert(!p.hasGenParticle() || 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 113 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

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

00113                                                                        {
00114       const Projection& reg = _addProjection(proj, name);
00115       const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
00116       return rtn;
00117     }

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

Apply the named projection on event.

Definition at line 81 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

00081                                                                                {
00082       return pcast<PROJ>(_applyProjection(evt, name));
00083     }

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

Apply the supplied projection on event.

Definition at line 74 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

00074                                                                                 {
00075       return pcast<PROJ>(_applyProjection(evt, proj));
00076     }

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

Apply the supplied projection on event.

Definition at line 67 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

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

00067                                                                           {
00068       return pcast<PROJ>(_applyProjection(evt, proj));
00069     }

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   }

const ParticleVector& bosons (  )  const [inline]

Access to the found bosons (currently either 0 or 1).

Definition at line 70 of file WFinder.hh.

References WFinder::_bosons.

Referenced by MC_WWJETS::analyze(), MC_WPOL::analyze(), MC_WJETS::analyze(), D0_2008_S7837160::analyze(), and D0_2000_S4480767::analyze().

00070 { return _bosons; }

void clear (  )  [inline]

Clear the projection.

Definition at line 97 of file WFinder.hh.

References WFinder::_bosons, WFinder::_constituentLeptons, WFinder::_constituentNeutrinos, and FinalState::_theParticles.

Referenced by WFinder::project().

00097                  {
00098       _theParticles.clear();
00099       _bosons.clear();
00100       _constituentLeptons.clear();
00101       _constituentNeutrinos.clear();
00102     }

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

Clone on the heap.

Reimplemented from FinalState.

Definition at line 63 of file WFinder.hh.

References WFinder::WFinder().

00063                                             {
00064       return new WFinder(*this);
00065     }

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

Compare projections.

Reimplemented from FinalState.

Definition at line 99 of file WFinder.cc.

References WFinder::_etMiss, WFinder::_maxmass, WFinder::_minmass, WFinder::_pid, WFinder::_trackPhotons, WFinder::_useTransverseMass, Rivet::cmp(), Rivet::EQUIVALENT, and Projection::mkNamedPCmp().

00099                                                 {
00100     PCmp LCcmp = mkNamedPCmp(p, "LeptonClusters");
00101     if (LCcmp != EQUIVALENT) return LCcmp;
00102 
00103     const WFinder& other = dynamic_cast<const WFinder&>(p);
00104     return (cmp(_minmass, other._minmass) || cmp(_maxmass, other._maxmass) ||
00105             cmp(_useTransverseMass, other._useTransverseMass) ||
00106             cmp(_etMiss, other._etMiss) ||
00107             cmp(_pid, other._pid) || cmp(_trackPhotons, other._trackPhotons));
00108 
00109   }

const vector<Particle>& constituentLeptons (  )  const [inline]

Access to the W constituent clustered leptons (currently either of size 0 if no boson was found or 1 if one boson was found)

Definition at line 74 of file WFinder.hh.

References WFinder::_constituentLeptons.

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

00074 { return _constituentLeptons; }

const vector<Particle>& constituentNeutrinos (  )  const [inline]

Access to the W constituent neutrinos (currently either of size 0 if no boson was found or 1 if one boson was found)

Definition at line 78 of file WFinder.hh.

References WFinder::_constituentNeutrinos.

Referenced by MC_WWJETS::analyze().

00078 { return _constituentNeutrinos; }

virtual bool empty (  )  const [inline, virtual, inherited]
const collection_type& entities (  )  const [inline, inherited]

Template-usable interface common to JetAlg.

Definition at line 96 of file FinalState.hh.

References FinalState::particles().

00096                                             {
00097       return particles();
00098     }

Log& getLog (  )  const [inline, inherited]
const Projection& getProjection ( const std::string &  name  )  const [inline, inherited]

Get the named projection (non-templated, so returns as a reference to a Projection base class).

Definition at line 57 of file ProjectionApplier.hh.

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

00057                                                                  {
00058       return getProjHandler().getProjection(*this, name);
00059     }

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

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

Definition at line 49 of file ProjectionApplier.hh.

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

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

00049                                                            {
00050       const Projection& p = getProjHandler().getProjection(*this, name);
00051       return pcast<PROJ>(p);
00052     }

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

Get the contained projections, including recursion.

Definition at line 42 of file ProjectionApplier.hh.

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

Referenced by Projection::beamPairs().

00042                                                       {
00043       return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
00044     }

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

Reimplemented in UnstableFinalState.

Definition at line 84 of file FinalState.hh.

References FinalState::_theParticles.

00084 { return _theParticles.empty(); }

Cmp< Projection > mkNamedPCmp ( const Projection otherparent,
const std::string &  pname 
) const [protected, inherited]
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]
const ParticleVector& particles ( 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.

Reimplemented in UnstableFinalState.

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(), TASSO_1990_S2148048::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_SUSY::analyze(), MC_PHOTONJETUE::analyze(), MC_PHOTONJETS::analyze(), MC_LEADINGJETS::analyze(), MC_GENERIC::analyze(), MC_DIJET::analyze(), JADE_1998_S3612880::analyze(), H1_2000_S4129130::analyze(), H1_1995_S3167097::analyze(), H1_1994_S2919893::analyze(), E735_1998_S3905616::analyze(), DELPHI_2003_WUD_03_11::analyze(), DELPHI_2002_069_CONF_603::analyze(), DELPHI_1996_S3430090::analyze(), DELPHI_1995_S3137023::analyze(), D0_2008_S7719523::analyze(), D0_2006_S6438750::analyze(), D0_2001_S4674421::analyze(), CMS_2011_S8884919::analyze(), CMS_2010_S8656010::analyze(), CMS_2010_S8547297::analyze(), CDF_2010_S8591881_QCD::analyze(), CDF_2010_S8591881_DY::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_1990_S2089246::analyze(), CDF_1988_S1865951::analyze(), ATLAS_2011_S9002537::analyze(), ATLAS_2010_S8919674::analyze(), ATLAS_2010_S8894728::analyze(), ATLAS_2010_S8591806::analyze(), ALICE_2010_S8706239::analyze(), ALICE_2010_S8625980::analyze(), ALEPH_2004_S5765862::analyze(), ALEPH_1996_S3486095::analyze(), Thrust::calc(), Spherocity::calc(), Sphericity::calc(), FParameter::calc(), BeamThrust::calc(), FinalState::entities(), ATLAS_2010_S8918562::fillPtEtaNch(), FinalState::particlesByE(), FinalState::particlesByEt(), FinalState::particlesByEta(), FinalState::particlesByModEta(), FinalState::particlesByP(), FinalState::particlesByPt(), VisibleFinalState::project(), VetoedFinalState::project(), TriggerUA5::project(), TriggerCDFRun2::project(), TriggerCDFRun0Run1::project(), TotalVisibleMomentum::project(), SVertex::project(), NeutralFinalState::project(), Multiplicity::project(), MissingMomentum::project(), MergedFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), LeptonClusters::project(), LeadingParticlesFinalState::project(), InvMassFinalState::project(), IdentifiedFinalState::project(), Hemispheres::project(), HadronicFinalState::project(), FoxWolframMoments::project(), FinalStateHCM::project(), FinalState::project(), DISLepton::project(), ClusteredPhotons::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 59 of file FinalState.hh.

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

00059                                                {
00060       return particles(cmpParticleByE);
00061     }

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

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

Definition at line 64 of file FinalState.hh.

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

00064                                                 {
00065       return particles(cmpParticleByEt);
00066     }

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

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

Definition at line 69 of file FinalState.hh.

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

00069                                                  {
00070       return particles(cmpParticleByAscPseudorapidity);
00071     }

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

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

Definition at line 74 of file FinalState.hh.

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

00074                                                     {
00075       return particles(cmpParticleByAscAbsPseudorapidity);
00076     }

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

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

Definition at line 54 of file FinalState.hh.

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

00054                                                {
00055       return particles(cmpParticleByP);
00056     }

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(), ATLAS_2011_S8994773::analyze(), and ATLAS_2010_S8894728::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 112 of file WFinder.cc.

References WFinder::_bosons, WFinder::_constituentLeptons, WFinder::_constituentNeutrinos, WFinder::_etMiss, WFinder::_masstarget, WFinder::_maxmass, WFinder::_minmass, WFinder::_nu_pid, WFinder::_pid, FinalState::_theParticles, WFinder::_trackPhotons, WFinder::_useTransverseMass, InvMassFinalState::calc(), WFinder::clear(), LeptonClusters::clusteredLeptons(), ClusteredLepton::constituentLepton(), ClusteredLepton::constituentPhotons(), Log::DEBUG, FinalState::FinalState(), Projection::getLog(), Rivet::GeV, Particle::momentum(), InvMassFinalState::particlePairs(), Particle::pdgId(), MissingMomentum::scalarET(), Rivet::PID::threeCharge(), InvMassFinalState::useTransverseMass(), Rivet::WMINUSBOSON, and Rivet::WPLUSBOSON.

00112                                       {
00113     clear();
00114 
00115     const LeptonClusters& leptons = applyProjection<LeptonClusters>(e, "LeptonClusters");
00116     const FinalState& neutrinos = applyProjection<FinalState>(e, "Neutrinos");
00117 
00118     // Make and register an invariant mass final state for the W decay leptons
00119     vector<pair<PdgId, PdgId> > l_nu_ids;
00120     l_nu_ids += make_pair(abs(_pid), -abs(_nu_pid));
00121     l_nu_ids += make_pair(-abs(_pid), abs(_nu_pid));
00122     InvMassFinalState imfs(FinalState(), l_nu_ids, _minmass, _maxmass, _masstarget);
00123     imfs.useTransverseMass(_useTransverseMass);
00124     ParticleVector tmp;
00125     tmp.insert(tmp.end(), leptons.clusteredLeptons().begin(), leptons.clusteredLeptons().end());
00126     tmp.insert(tmp.end(), neutrinos.particles().begin(), neutrinos.particles().end());
00127     imfs.calc(tmp);
00128 
00129     if (imfs.particlePairs().size() < 1) return;
00130 
00131     ParticlePair Wconstituents(imfs.particlePairs()[0]);
00132     Particle p1(Wconstituents.first), p2(Wconstituents.second);
00133 
00134     if (PID::threeCharge(p1)==0) {
00135       _constituentLeptons += p2;
00136       _constituentNeutrinos += p1;
00137     }
00138     else {
00139       _constituentLeptons += p1;
00140       _constituentNeutrinos += p2;
00141     }
00142 
00143     FourMomentum pW = p1.momentum() + p2.momentum();
00144     const int w3charge = PID::threeCharge(p1) + PID::threeCharge(p2);
00145     assert(abs(w3charge) == 3);
00146     const int wcharge = w3charge/3;
00147 
00148     stringstream msg;
00149     string wsign = (wcharge == 1) ? "+" : "-";
00150     string wstr = "W" + wsign;
00151     msg << wstr << " reconstructed from: " << endl
00152         << "   " << p1.momentum() << " " << p1.pdgId() << endl
00153         << " + " << p2.momentum() << " " << p2.pdgId() << endl;
00154 
00155     // Check missing ET
00156     const MissingMomentum& vismom = applyProjection<MissingMomentum>(e, "MissingET");
00157     /// @todo Restrict missing momentum eta range? Use vectorET()?
00158     if (vismom.scalarET() < _etMiss) {
00159       getLog() << Log::DEBUG << "Not enough missing ET: " << vismom.scalarET()/GeV
00160                << " GeV vs. " << _etMiss/GeV << " GeV" << endl;
00161       return;
00162     }
00163 
00164     // Make W Particle and insert into particles list
00165     const PdgId wpid = (wcharge == 1) ? WPLUSBOSON : WMINUSBOSON;
00166     _bosons.push_back(Particle(wpid, pW));
00167 
00168     // find the LeptonClusters and neutrinos which survived the IMFS cut such that we can
00169     // extract their original particles
00170     foreach (const Particle& p, _constituentNeutrinos) {
00171       _theParticles.push_back(p);
00172     }
00173     foreach (const Particle& p, _constituentLeptons) {
00174       foreach (const ClusteredLepton& l, leptons.clusteredLeptons()) {
00175         if (p.pdgId()==l.pdgId() && p.momentum()==l.momentum()) {
00176           _theParticles.push_back(l.constituentLepton());
00177           if (_trackPhotons) {
00178             _theParticles.insert(_theParticles.end(),
00179                                  l.constituentPhotons().begin(),
00180                                  l.constituentPhotons().end());
00181           }
00182         }
00183       }
00184     }
00185   }

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

Minimum-$ p_\perp $ requirement.

Definition at line 87 of file FinalState.hh.

References FinalState::_ptmin.

00087 { 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 95 of file WFinder.cc.

Referenced by MC_WJETS::init().

00095                                                        {
00096     return getProjection<FinalState>("RFS");
00097   }

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(), FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), Beam::Beam(), BeamThrust::BeamThrust(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ClusteredPhotons::ClusteredPhotons(), ConstLossyFinalState::ConstLossyFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), FinalStateHCM::FinalStateHCM(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), InitialQuarks::InitialQuarks(), IsolationProjection< PROJ1, PROJ2, EST >::IsolationProjection(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LeptonClusters::LeptonClusters(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), Multiplicity::Multiplicity(), NeutralFinalState::NeutralFinalState(), ParisiTensor::ParisiTensor(), PVertex::PVertex(), Sphericity::Sphericity(), Spherocity::Spherocity(), SVertex::SVertex(), Thrust::Thrust(), TotalVisibleMomentum::TotalVisibleMomentum(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), UnstableFinalState::UnstableFinalState(), VetoedFinalState::VetoedFinalState(), and VisibleFinalState::VisibleFinalState().

00120                                         {
00121       _name = name;
00122     }

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

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.


Member Data Documentation

bool _allowProjReg [protected, inherited]

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

Definition at line 140 of file ProjectionApplier.hh.

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

list of found bosons (currently either 0 or 1)

Definition at line 141 of file WFinder.hh.

Referenced by WFinder::bosons(), WFinder::clear(), and WFinder::project().

Constituent leptons (currently either 0 or 1).

Definition at line 144 of file WFinder.hh.

Referenced by WFinder::clear(), WFinder::constituentLeptons(), and WFinder::project().

Constituent neutrinos (currently either 0 or 1).

Definition at line 147 of file WFinder.hh.

Referenced by WFinder::clear(), WFinder::constituentNeutrinos(), and WFinder::project().

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

The ranges allowed for pseudorapidity.

Definition at line 116 of file FinalState.hh.

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

double _etMiss [private]

Missing ET cut.

Definition at line 124 of file WFinder.hh.

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

double _masstarget [private]

Definition at line 120 of file WFinder.hh.

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

double _maxmass [private]

Definition at line 120 of file WFinder.hh.

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

double _minmass [private]

Transverse mass cuts.

Definition at line 120 of file WFinder.hh.

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

PdgId _nu_pid [private]

Neutrino flavour.

Definition at line 138 of file WFinder.hh.

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

PdgId _pid [private]

Lepton flavour.

Definition at line 135 of file WFinder.hh.

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

double _ptmin [protected, inherited]

The minimum allowed transverse momentum.

Reimplemented in UnstableFinalState.

Definition at line 119 of file FinalState.hh.

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

ParticleVector _theParticles [mutable, protected, inherited]
bool _trackPhotons [private]

Switch for tracking of photons (whether to add them to _theParticles) This is relevant when the ZFinder::_theParticles are to be excluded from e.g. the input to a jet finder, to specify whether the clustered photons are to be excluded as well. (Yes, some experiments make a difference between clusterPhotons and trackPhotons!)

Definition at line 132 of file WFinder.hh.

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

bool _useTransverseMass [private]

Definition at line 121 of file WFinder.hh.

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


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