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, FinalState inputfs=FinalState())
 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, FinalState inputfs=FinalState())
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, FinalState inputfs)
 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,
FinalState  inputfs = FinalState() 
)

Constructor taking single eta/pT bounds

Parameters:
pid type of the charged lepton
minmass,maxmass (transverse) mass window
missingET minimal amount of missing ET (neutrinos) required
dRmax maximum dR of photons around charged lepton to take into account for W reconstruction (only relevant if one of the following are true)
clusterPhotons whether such photons are supposed to be clustered to the lepton object and thus W mom
trackPhotons whether such photons should be added to _theParticles (cf. _trackPhotons)
useTransverseMass whether mass window should be applied using mT

Definition at line 15 of file WFinder.cc.

References WFinder::_init().

Referenced by WFinder::clone().

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

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,
FinalState  inputfs = FinalState() 
)

Constructor taking multiple eta/pT bounds

Parameters:
pid type of the charged lepton
minmass,maxmass (transverse) mass window
missingET minimal amount of missing ET (neutrinos) required
dRmax maximum dR of photons around charged lepton to take into account for W reconstruction (only relevant if one of the following are true)
clusterPhotons whether such photons are supposed to be clustered to the lepton object and thus W mom
trackPhotons whether such photons should be added to _theParticles (cf. _trackPhotons)
useTransverseMass whether mass window should be applied using mT

Definition at line 31 of file WFinder.cc.

References WFinder::_init().

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


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,
FinalState  inputfs 
) [private]

Common implementation of constructor operation, taking FS params.

Definition at line 45 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().

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

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

Decide if a particle is to be accepted or not.

Definition at line 95 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().

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

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(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), 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_VH2BB::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_S9120041::init(), CMS_2011_S9088458::init(), CMS_2011_S9086218::init(), CMS_2011_S8978280::init(), CMS_2011_S8968497::init(), CMS_2011_S8957746::init(), CMS_2011_S8950903::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_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_S9225137::init(), ATLAS_2011_S9212183::init(), ATLAS_2011_S9131140::init(), ATLAS_2011_S9128077::init(), ATLAS_2011_S9126244::init(), ATLAS_2011_S9120807::init(), ATLAS_2011_S9108483::init(), ATLAS_2011_S9041966::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_I925932::init(), ATLAS_2011_I919017::init(), ATLAS_2011_CONF_2011_098::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_2011_S8945144::init(), ALICE_2011_S8909580::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(), NonHadronicFinalState::NonHadronicFinalState(), 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 DISFinalState::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 78 of file WFinder.hh.

References WFinder::_bosons.

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

00078 { return _bosons; }

void clear (  )  [inline]

Clear the projection.

Definition at line 105 of file WFinder.hh.

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

Referenced by WFinder::project().

00105                  {
00106       _theParticles.clear();
00107       _bosons.clear();
00108       _constituentLeptons.clear();
00109       _constituentNeutrinos.clear();
00110     }

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

Clone on the heap.

Reimplemented from FinalState.

Definition at line 71 of file WFinder.hh.

References WFinder::WFinder().

00071                                             {
00072       return new WFinder(*this);
00073     }

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

Compare projections.

Reimplemented from FinalState.

Definition at line 101 of file WFinder.cc.

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

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

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

References WFinder::_constituentLeptons.

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

00082 { 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 86 of file WFinder.hh.

References WFinder::_constituentNeutrinos.

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

00086 { 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]

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(), VetoedFinalState::project(), UnstableFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), IsolationProjection< PROJ1, PROJ2, EST >::project(), InitialQuarks::project(), ChargedFinalState::project(), 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 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_VH2BB::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_1990_S2089246::analyze(), CDF_1988_S1865951::analyze(), ATLAS_2011_S9002537::analyze(), ATLAS_2011_I925932::analyze(), ATLAS_2010_S8919674::analyze(), ATLAS_2010_S8894728::analyze(), ATLAS_2010_S8591806::analyze(), ALICE_2011_S8945144::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(), NonHadronicFinalState::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(), 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 113 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(), Rivet::GeV, Particle::momentum(), MSG_DEBUG, InvMassFinalState::particlePairs(), Particle::pdgId(), Rivet::PID::threeCharge(), InvMassFinalState::useTransverseMass(), Rivet::WMINUSBOSON, and Rivet::WPLUSBOSON.

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

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 97 of file WFinder.cc.

Referenced by MC_WJETS::init().

00097                                                        {
00098     return getProjection<FinalState>("RFS");
00099   }

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(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), FinalState::FinalState(), 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(), NonHadronicFinalState::NonHadronicFinalState(), 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 150 of file WFinder.hh.

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

Constituent leptons (currently either 0 or 1).

Definition at line 153 of file WFinder.hh.

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

Constituent neutrinos (currently either 0 or 1).

Definition at line 156 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 133 of file WFinder.hh.

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

double _masstarget [private]

Definition at line 129 of file WFinder.hh.

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

double _maxmass [private]

Definition at line 129 of file WFinder.hh.

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

double _minmass [private]

Transverse mass cuts.

Definition at line 129 of file WFinder.hh.

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

PdgId _nu_pid [private]

Neutrino flavour.

Definition at line 147 of file WFinder.hh.

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

PdgId _pid [private]

Lepton flavour.

Definition at line 144 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 141 of file WFinder.hh.

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

bool _useTransverseMass [private]

Definition at line 130 of file WFinder.hh.

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


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