VetoedFinalState Class Reference

FS modifier to exclude classes of particles from the final state. More...

#include <VetoedFinalState.hh>

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

List of all members.

Public Types

typedef pair< double, double > BinaryCut
 Typedef for a pair of back-to-back cuts.
typedef map< long, BinaryCutVetoDetails
 Typedef for a vetoing entry.
typedef multimap< int, BinaryCutCompositeVeto
 Typedef for a veto on a composite particle mass.
typedef Particle entity_type
typedef ParticleVector collection_type

Public Member Functions

const VetoDetailsvetoDetails () const
 Get the list of particle IDs and $ p_T $ ranges to veto.
VetoedFinalStateaddVetoDetail (const long id, const double ptmin, const double ptmax)
VetoedFinalStateaddVetoPairDetail (const long id, const double ptmin, const double ptmax)
VetoedFinalStateaddVetoPairId (const long id)
VetoedFinalStateaddVetoId (const long id)
 Add a particle ID to veto (all $ p_T $ range will be vetoed).
VetoedFinalStatevetoNeutrinos ()
 Veto all neutrinos (convenience method).
VetoedFinalStateaddCompositeMassVeto (const double &mass, const double &width, int nProducts=2)
VetoedFinalStateaddDecayProductsVeto (const long id)
VetoedFinalStatesetVetoDetails (const VetoDetails &ids)
 Set the list of particle IDs and $ p_T $ ranges to veto.
VetoedFinalStatereset ()
 Clear the list of particle IDs and ranges to veto.
VetoedFinalStateaddVetoOnThisFinalState (const FinalState &fs)
 Veto particles from a supplied final state.
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

 VetoedFinalState ()
 Default constructor.
 VetoedFinalState (const FinalState &fsp)
 Constructor with specific FinalState.
 VetoedFinalState (const VetoDetails &vetocodes)
 VetoedFinalState (const FinalState &fsp, const VetoDetails &vetocodes)
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 Types

typedef set< long > ParentVetos

Private Attributes

VetoDetails _vetoCodes
 The final-state particles.
CompositeVeto _compositeVetoes
 Composite particle masses to veto.
set< int > _nCompositeDecays
ParentVetos _parentVetoes
 Set of decaying particle IDs to veto.
set< string > _vetofsnames
 Set of finalstate to be vetoed.

Friends

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

Detailed Description

FS modifier to exclude classes of particles from the final state.

Definition at line 16 of file VetoedFinalState.hh.


Member Typedef Documentation

typedef pair<double, double> BinaryCut

Typedef for a pair of back-to-back cuts.

Definition at line 21 of file VetoedFinalState.hh.

typedef ParticleVector collection_type [inherited]

Definition at line 93 of file FinalState.hh.

typedef multimap<int, BinaryCut> CompositeVeto

Typedef for a veto on a composite particle mass.

Definition at line 27 of file VetoedFinalState.hh.

typedef Particle entity_type [inherited]

Definition at line 92 of file FinalState.hh.

typedef set<long> ParentVetos [private]

Definition at line 179 of file VetoedFinalState.hh.

typedef map<long, BinaryCut> VetoDetails

Typedef for a vetoing entry.

Definition at line 24 of file VetoedFinalState.hh.


Constructor & Destructor Documentation

VetoedFinalState (  )  [inline]

Default constructor.

Definition at line 33 of file VetoedFinalState.hh.

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

Referenced by VetoedFinalState::clone().

00033                        {
00034       setName("VetoedFinalState");
00035       addProjection(FinalState(), "FS");
00036     }

VetoedFinalState ( const FinalState fsp  )  [inline]

Constructor with specific FinalState.

Definition at line 39 of file VetoedFinalState.hh.

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

00040     {
00041       setName("VetoedFinalState");
00042       addProjection(fsp, "FS");
00043     }

VetoedFinalState ( const VetoDetails vetocodes  )  [inline]

You can add a map of ID plus a pair containing $ p_{Tmin} $ and $ p_{Tmax} $ - these define the range of particles to be vetoed.

Definition at line 47 of file VetoedFinalState.hh.

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

00048       : _vetoCodes(vetocodes)
00049     {
00050       setName("VetoedFinalState");
00051       addProjection(FinalState(), "FS");
00052     }

VetoedFinalState ( const FinalState fsp,
const VetoDetails vetocodes 
) [inline]

You can add a map of ID plus a pair containing $ p_{Tmin} $ and $ p_{Tmax} $ - these define the range of particles to be vetoed. This version also supplies a specifi FinalState to be used.

Definition at line 57 of file VetoedFinalState.hh.

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

00058       : _vetoCodes(vetocodes)
00059     {
00060       setName("VetoedFinalState");
00061       addProjection(fsp, "FS");
00062     }


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   }

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   }

VetoedFinalState& addCompositeMassVeto ( const double &  mass,
const double &  width,
int  nProducts = 2 
) [inline]

Add a veto on composite masses within a given width. The composite mass is composed of nProducts decay products @ todo might we want to specify a range of pdg ids for the decay products?

Definition at line 121 of file VetoedFinalState.hh.

References VetoedFinalState::_compositeVetoes, and VetoedFinalState::_nCompositeDecays.

00121                                                                                                     {
00122       double halfWidth = 0.5*width;
00123       BinaryCut massRange(mass - halfWidth, mass + halfWidth);
00124       _compositeVetoes.insert(make_pair(nProducts, massRange));
00125       _nCompositeDecays.insert(nProducts);
00126       return *this;
00127     }

VetoedFinalState& addDecayProductsVeto ( const long  id  )  [inline]

Veto the decay products of particle with pdg id

Todo:
Need HepMC to sort themselves out and keep vector bosons from the hard vtx in the event record before this will work reliably for all pdg ids

Definition at line 132 of file VetoedFinalState.hh.

References VetoedFinalState::_parentVetoes.

00132                                                          {
00133       _parentVetoes.insert(id);
00134       return *this;
00135     }

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     }

VetoedFinalState& addVetoDetail ( const long  id,
const double  ptmin,
const double  ptmax 
) [inline]

Add a particle ID and $ p_T $ range to veto. Particles with $ p_T $ IN the given range will be rejected.

Definition at line 81 of file VetoedFinalState.hh.

References VetoedFinalState::_vetoCodes.

Referenced by VetoedFinalState::addVetoPairDetail().

00081                                                                                            {
00082       BinaryCut ptrange(ptmin, ptmax);
00083       _vetoCodes.insert(make_pair(id, ptrange));
00084       return *this;
00085     }

VetoedFinalState& addVetoId ( const long  id  )  [inline]

Add a particle ID to veto (all $ p_T $ range will be vetoed).

Definition at line 104 of file VetoedFinalState.hh.

References VetoedFinalState::_vetoCodes.

Referenced by VetoedFinalState::addVetoPairId().

00104                                                {
00105       BinaryCut ptrange(0.0, numeric_limits<double>::max());
00106       _vetoCodes.insert(make_pair(id, ptrange));
00107       return *this;
00108     }

VetoedFinalState& addVetoOnThisFinalState ( const FinalState fs  )  [inline]

Veto particles from a supplied final state.

Definition at line 151 of file VetoedFinalState.hh.

References VetoedFinalState::_vetofsnames, ProjectionApplier::addProjection(), and Projection::name().

Referenced by ZFinder::_init(), WFinder::_init(), MC_ZZJETS::init(), MC_WWJETS::init(), MC_PHOTONJETS::init(), D0_2008_S7719523::init(), CDF_2008_S8095620::init(), CDF_2008_S7541902::init(), CDF_2006_S6653332::init(), CDF_1993_S2742446::init(), and ATLAS_2010_S8919674::init().

00151                                                                     {
00152       stringstream st_name;
00153       st_name << "FS_" << _vetofsnames.size();
00154       string name = st_name.str();
00155       addProjection(fs, name);
00156       _vetofsnames.insert(name);
00157       return *this;
00158     }

VetoedFinalState& addVetoPairDetail ( const long  id,
const double  ptmin,
const double  ptmax 
) [inline]

Add a particle/antiparticle pair to veto in a given $ p_T $ range. Given a single ID, both the particle and its conjugate antiparticle will be rejected if their $ p_T $ is IN the given range.

Definition at line 89 of file VetoedFinalState.hh.

References VetoedFinalState::addVetoDetail().

Referenced by D0_2004_S5992206::init(), and ATLAS_2011_S9041966::init().

00089                                                                                                {
00090       addVetoDetail(id,  ptmin, ptmax);
00091       addVetoDetail(-id, ptmin, ptmax);
00092       return *this;
00093     }

VetoedFinalState& addVetoPairId ( const long  id  )  [inline]

Add a particle/antiparticle pair to veto. Given a single ID, both the particle and its corresponding antiparticle (for all $ p_T $ values) will be vetoed.

Definition at line 97 of file VetoedFinalState.hh.

References VetoedFinalState::addVetoId().

Referenced by STAR_2009_UE_HELEN::init(), ATLAS_2011_S9225137::init(), ATLAS_2011_S9108483::init(), ATLAS_2011_S9019561::init(), ATLAS_2011_S8983313::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2010_S8919674::init(), and VetoedFinalState::vetoNeutrinos().

00097                                                    {
00098       addVetoId(id);
00099       addVetoId(-id);
00100       return *this;
00101     }

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   }

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

Clone on the heap.

Reimplemented from FinalState.

Definition at line 66 of file VetoedFinalState.hh.

References VetoedFinalState::VetoedFinalState().

00066                                             {
00067       return new VetoedFinalState(*this);
00068     }

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

Compare projections.

Reimplemented from FinalState.

Definition at line 11 of file VetoedFinalState.cc.

References VetoedFinalState::_compositeVetoes, VetoedFinalState::_parentVetoes, VetoedFinalState::_vetoCodes, VetoedFinalState::_vetofsnames, Rivet::cmp(), Rivet::EQUIVALENT, Projection::mkNamedPCmp(), and Rivet::UNDEFINED.

00011                                                          {
00012     const PCmp fscmp = mkNamedPCmp(p, "FS");
00013     if (fscmp != EQUIVALENT) return fscmp;
00014     if (_vetofsnames.size() != 0) return UNDEFINED;
00015     const VetoedFinalState& other = dynamic_cast<const VetoedFinalState&>(p);
00016     return \
00017       cmp(_vetoCodes, other._vetoCodes) ||
00018       cmp(_compositeVetoes, other._compositeVetoes) ||
00019       cmp(_parentVetoes, other._parentVetoes);
00020   }

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:
Improve!

Reimplemented from FinalState.

Definition at line 23 of file VetoedFinalState.cc.

References VetoedFinalState::_compositeVetoes, VetoedFinalState::_nCompositeDecays, VetoedFinalState::_parentVetoes, FinalState::_theParticles, VetoedFinalState::_vetoCodes, VetoedFinalState::_vetofsnames, Projection::getLog(), Log::isActive(), std::join(), Rivet::mass(), Rivet::mass2(), Particle::momentum(), MSG_TRACE, FinalState::particles(), Particle::pdgId(), FourMomentum::pT(), and Log::TRACE.

00023                                                {
00024     const FinalState& fs = applyProjection<FinalState>(e, "FS");
00025     _theParticles.clear();
00026     _theParticles.reserve(fs.particles().size());
00027     foreach (const Particle& p, fs.particles()) {
00028       if (getLog().isActive(Log::TRACE)) {
00029         vector<long> codes;
00030         for (VetoDetails::const_iterator code = _vetoCodes.begin(); code != _vetoCodes.end(); ++code) {
00031           codes.push_back(code->first);
00032         }
00033         const string codestr = "{ " + join(codes) + " }";
00034         MSG_TRACE(p.pdgId() << " vs. veto codes = " << codestr << " (" << codes.size() << ")");
00035       }
00036       const long pdgid = p.pdgId();
00037       const double pt = p.momentum().pT();
00038       VetoDetails::iterator iter = _vetoCodes.find(pdgid);
00039       if (iter == _vetoCodes.end()) {
00040         MSG_TRACE("Storing with PDG code = " << pdgid << ", pT = " << pt);
00041         _theParticles.push_back(p);
00042       } else {
00043         // This particle code is listed as a possible veto... check pT.
00044         // Make sure that the pT range is sensible:
00045         BinaryCut ptrange = iter->second;
00046         assert(ptrange.first <= ptrange.second);
00047         stringstream rangess;
00048         if (ptrange.first < numeric_limits<double>::max()) rangess << ptrange.second;
00049         rangess << " - ";
00050         if (ptrange.second < numeric_limits<double>::max()) rangess << ptrange.second;
00051         MSG_TRACE("ID = " << pdgid << ", pT range = " << rangess.str());
00052         stringstream debugline;
00053         debugline << "with PDG code = " << pdgid << " pT = " << p.momentum().pT();
00054         if (pt < ptrange.first || pt > ptrange.second) {
00055           MSG_TRACE("Storing " << debugline.str());
00056           _theParticles.push_back(p);
00057         } else {
00058           MSG_TRACE("Vetoing " << debugline.str());
00059         }
00060       }
00061     }
00062 
00063     set<ParticleVector::iterator> toErase;
00064     for (set<int>::iterator nIt = _nCompositeDecays.begin();
00065          nIt != _nCompositeDecays.end() && !_theParticles.empty(); ++nIt) {
00066       map<set<ParticleVector::iterator>, FourMomentum> oldMasses;
00067       map<set<ParticleVector::iterator>, FourMomentum> newMasses;
00068       set<ParticleVector::iterator> start;
00069       start.insert(_theParticles.begin());
00070       oldMasses.insert(pair<set<ParticleVector::iterator>, FourMomentum>
00071                        (start, _theParticles.begin()->momentum()));
00072 
00073       for (int nParts = 1; nParts != *nIt; ++nParts) {
00074         for (map<set<ParticleVector::iterator>, FourMomentum>::iterator mIt = oldMasses.begin();
00075              mIt != oldMasses.end(); ++mIt) {
00076           ParticleVector::iterator pStart = *(mIt->first.rbegin());
00077           for (ParticleVector::iterator pIt = pStart + 1; pIt != _theParticles.end(); ++pIt) {
00078             FourMomentum cMom = mIt->second + pIt->momentum();
00079             set<ParticleVector::iterator> pList(mIt->first);
00080             pList.insert(pIt);
00081             newMasses[pList] = cMom;
00082           }
00083         }
00084         oldMasses = newMasses;
00085         newMasses.clear();
00086       }
00087       for (map<set<ParticleVector::iterator>, FourMomentum>::iterator mIt = oldMasses.begin();
00088            mIt != oldMasses.end(); ++mIt) {
00089         double mass2 = mIt->second.mass2();
00090         if (mass2 >= 0.0) {
00091           double mass = sqrt(mass2);
00092           for (CompositeVeto::iterator cIt = _compositeVetoes.lower_bound(*nIt);
00093                cIt != _compositeVetoes.upper_bound(*nIt); ++cIt) {
00094             BinaryCut massRange = cIt->second;
00095             if (mass < massRange.second && mass > massRange.first) {
00096               for (set<ParticleVector::iterator>::iterator lIt = mIt->first.begin();
00097                    lIt != mIt->first.end(); ++lIt) {
00098                 toErase.insert(*lIt);
00099               }
00100             }
00101           }
00102         }
00103       }
00104     }
00105 
00106     for (set<ParticleVector::iterator>::reverse_iterator p = toErase.rbegin(); p != toErase.rend(); ++p) {
00107       _theParticles.erase(*p);
00108     }
00109 
00110     /// @todo Improve!
00111     for (ParentVetos::const_iterator vIt = _parentVetoes.begin(); vIt != _parentVetoes.end(); ++vIt) {
00112       for (ParticleVector::iterator p = _theParticles.begin(); p != _theParticles.end(); ++p) {
00113         GenVertex* startVtx = ((*p).genParticle()).production_vertex();
00114         bool veto = false;
00115         if (startVtx!=0) {
00116           for (GenVertex::particle_iterator pIt = startVtx->particles_begin(HepMC::ancestors);
00117                    pIt != startVtx->particles_end(HepMC::ancestors) && !veto; ++pIt) {
00118             if (*vIt == (*pIt)->pdg_id()) {
00119               veto = true;
00120               p = _theParticles.erase(p);
00121               --p;
00122             }
00123           }
00124         }
00125       }
00126     }
00127 
00128     // Now veto on the FS
00129     foreach (const string& ifs, _vetofsnames) {
00130       const FinalState& vfs = applyProjection<FinalState>(e, ifs);
00131       const ParticleVector& vfsp = vfs.particles();
00132       for (ParticleVector::iterator icheck = _theParticles.begin(); icheck != _theParticles.end(); ++icheck) {
00133         if (!icheck->hasGenParticle()) continue;
00134         bool found = false;
00135         for (ParticleVector::const_iterator ipart = vfsp.begin(); ipart != vfsp.end(); ++ipart){
00136           if (!ipart->hasGenParticle()) continue;
00137           MSG_TRACE("Comparing barcode " << icheck->genParticle().barcode()
00138                    << " with veto particle " << ipart->genParticle().barcode());
00139           if (ipart->genParticle().barcode() == icheck->genParticle().barcode()){
00140             found = true;
00141             break;
00142           }
00143         }
00144         if (found) {
00145           _theParticles.erase(icheck);
00146           --icheck;
00147         }
00148       }
00149     }
00150   }

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

VetoedFinalState& reset (  )  [inline]

Clear the list of particle IDs and ranges to veto.

Definition at line 144 of file VetoedFinalState.hh.

References VetoedFinalState::_vetoCodes.

00144                               {
00145       _vetoCodes.clear();
00146       return *this;
00147     }

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     }

VetoedFinalState& setVetoDetails ( const VetoDetails ids  )  [inline]

Set the list of particle IDs and $ p_T $ ranges to veto.

Definition at line 138 of file VetoedFinalState.hh.

References VetoedFinalState::_vetoCodes.

00138                                                              {
00139       _vetoCodes = ids;
00140       return *this;
00141     }

virtual size_t size (  )  const [inline, virtual, inherited]
const VetoDetails& vetoDetails (  )  const [inline]

Get the list of particle IDs and $ p_T $ ranges to veto.

Definition at line 75 of file VetoedFinalState.hh.

References VetoedFinalState::_vetoCodes.

00075                                            {
00076       return _vetoCodes;
00077     }

VetoedFinalState& vetoNeutrinos (  )  [inline]

Veto all neutrinos (convenience method).

Definition at line 111 of file VetoedFinalState.hh.

References VetoedFinalState::addVetoPairId(), Rivet::NU_E, Rivet::NU_MU, and Rivet::NU_TAU.

Referenced by STAR_2009_UE_HELEN::init(), D0_2004_S5992206::init(), D0_2001_S4674421::init(), and ATLAS_2010_S8919674::init().

00111                                       {
00112       addVetoPairId(NU_E);
00113       addVetoPairId(NU_MU);
00114       addVetoPairId(NU_TAU);
00115       return *this;
00116     }


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

Composite particle masses to veto.

Definition at line 176 of file VetoedFinalState.hh.

Referenced by VetoedFinalState::addCompositeMassVeto(), VetoedFinalState::compare(), and VetoedFinalState::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().

set<int> _nCompositeDecays [private]

Set of decaying particle IDs to veto.

Definition at line 182 of file VetoedFinalState.hh.

Referenced by VetoedFinalState::addDecayProductsVeto(), VetoedFinalState::compare(), and VetoedFinalState::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]
set<string> _vetofsnames [private]

Set of finalstate to be vetoed.

Definition at line 185 of file VetoedFinalState.hh.

Referenced by VetoedFinalState::addVetoOnThisFinalState(), VetoedFinalState::compare(), and VetoedFinalState::project().


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