Public Types |
Public Member Functions |
Protected Member Functions |
Protected Attributes |
Private Member Functions |
Private Attributes |
Friends
FastJets Class Reference Project out jets found using the FastJet package jet algorithms. More...
Inheritance diagram for FastJets:
Collaboration diagram for FastJets:
Detailed DescriptionProject out jets found using the FastJet package jet algorithms. Definition at line 47 of file FastJets.hh. Member Typedef Documentation
Member Enumeration Documentation
Wrapper enum for selected Fastjet jet algorithms.
Definition at line 51 of file FastJets.hh. Constructor & Destructor Documentation
"Wrapped" argument constructor using Rivet enums for most common jet alg choices (including some plugins). For the built-in algs, E-scheme recombination is used. For full control of FastJet built-in jet algs, use the native arg constructor. Definition at line 65 of file FastJets.hh. References FastJets::_init1(). Referenced by FastJets::clone().
Native argument constructor, using FastJet alg/scheme enums. Definition at line 70 of file FastJets.hh. References FastJets::_init2().
Explicitly pass in an externally-constructed plugin (must be heap-allocated, Rivet will delete) Definition at line 75 of file FastJets.hh. References FastJets::_init3().
Explicitly pass in an externally-constructed plugin (must be heap-allocated, Rivet will delete) Definition at line 78 of file FastJets.hh. References FastJets::_init3().
Same thing as above, but without an FS (for when we want to pass the particles directly to the calc method) Definition at line 83 of file FastJets.hh. References FastJets::_init1().
Same thing as above, but without an FS (for when we want to pass the particles directly to the calc method) Definition at line 86 of file FastJets.hh. References FastJets::_init2().
Same thing as above, but without an FS (for when we want to pass the particles directly to the calc method) Definition at line 89 of file FastJets.hh. References FastJets::_init3().
Same thing as above, but without an FS (for when we want to pass the particles directly to the calc method) Definition at line 92 of file FastJets.hh. References FastJets::_init3(). Member Function Documentation
Untemplated function to do the work... Definition at line 33 of file ProjectionApplier.cc. References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection(). Referenced by ProjectionApplier::addProjection(). {
if (!_allowProjReg) {
cerr << "Trying to register projection '"
<< proj.name() << "' before init phase in '" << this->name() << "'." << endl;
exit(2);
}
const Projection& reg = getProjHandler().registerProjection(*this, proj, name);
return reg;
}
Shared utility functions to implement constructor behaviour. Definition at line 11 of file FastJets.cc. References FastJets::_jdef, FastJets::_plugin, ProjectionApplier::addProjection(), FastJets::ANTIKT, FastJets::ATLASCONE, FastJets::CAM, FastJets::CDFJETCLU, FastJets::CDFMIDPOINT, FastJets::CMSCONE, FastJets::D0ILCONE, FastJets::DURHAM, FastJets::JADE, FastJets::KT, MSG_DEBUG, FastJets::PXCONE, Projection::setName(), FastJets::SISCONE, and FastJets::TRACKJET. Referenced by FastJets::FastJets(). {
setName("FastJets");
MSG_DEBUG("JetAlg = " << alg);
MSG_DEBUG("R parameter = " << rparameter);
MSG_DEBUG("Seed threshold = " << seed_threshold);
if (alg == KT) {
_jdef = fastjet::JetDefinition(fastjet::kt_algorithm, rparameter, fastjet::E_scheme);
} else if (alg == CAM) {
_jdef = fastjet::JetDefinition(fastjet::cambridge_algorithm, rparameter, fastjet::E_scheme);
} else if (alg == ANTIKT) {
_jdef = fastjet::JetDefinition(fastjet::antikt_algorithm, rparameter, fastjet::E_scheme);
} else if (alg == DURHAM) {
_jdef = fastjet::JetDefinition(fastjet::ee_kt_algorithm, fastjet::E_scheme);
} else {
// Plugins:
if (alg == SISCONE) {
const double OVERLAP_THRESHOLD = 0.75;
_plugin.reset(new fastjet::SISConePlugin(rparameter, OVERLAP_THRESHOLD));
} else if (alg == PXCONE) {
string msg = "PxCone currently not supported, since FastJet doesn't install it by default. ";
msg += "Please notify the Rivet authors if this behaviour should be changed.";
throw Error(msg);
//_plugin.reset(new fastjet::PxConePlugin(rparameter));
} else if (alg == ATLASCONE) {
const double OVERLAP_THRESHOLD = 0.5;
_plugin.reset(new fastjet::ATLASConePlugin(rparameter, seed_threshold, OVERLAP_THRESHOLD));
} else if (alg == CMSCONE) {
_plugin.reset(new fastjet::CMSIterativeConePlugin(rparameter, seed_threshold));
} else if (alg == CDFJETCLU) {
const double OVERLAP_THRESHOLD = 0.75;
_plugin.reset(new fastjet::CDFJetCluPlugin(rparameter, OVERLAP_THRESHOLD, seed_threshold));
} else if (alg == CDFMIDPOINT) {
const double OVERLAP_THRESHOLD = 0.5;
_plugin.reset(new fastjet::CDFMidPointPlugin(rparameter, OVERLAP_THRESHOLD, seed_threshold));
} else if (alg == D0ILCONE) {
const double min_jet_Et = 6.0;
_plugin.reset(new fastjet::D0RunIIConePlugin(rparameter, min_jet_Et));
} else if (alg == JADE) {
_plugin.reset(new fastjet::JadePlugin());
} else if (alg == TRACKJET) {
_plugin.reset(new fastjet::TrackJetPlugin(rparameter));
}
_jdef = fastjet::JetDefinition(_plugin.get());
}
addProjection(HeavyHadrons(), "HFHadrons");
addProjection(TauFinder(), "Taus");
}
Definition at line 59 of file FastJets.cc. References FastJets::_jdef, ProjectionApplier::addProjection(), and Projection::setName(). Referenced by FastJets::FastJets(). {
setName("FastJets");
_jdef = fastjet::JetDefinition(type, rparameter, recom);
addProjection(HeavyHadrons(), "HFHadrons");
addProjection(TauFinder(), "Taus");
}
Definition at line 67 of file FastJets.cc. References FastJets::_jdef, FastJets::_plugin, ProjectionApplier::addProjection(), and Projection::setName(). Referenced by FastJets::FastJets(). {
setName("FastJets");
/// @todo Should we be copying the plugin?
_plugin.reset(plugin);
_jdef = fastjet::JetDefinition(_plugin.get());
addProjection(HeavyHadrons(), "HFHadrons");
addProjection(TauFinder(), "Taus");
}
Get the jets (unordered).
Implements JetAlg. Definition at line 154 of file FastJets.cc. References FastJets::_particles, FastJets::clusterSeq(), and FastJets::pseudojets(). {
Jets rtn;
foreach (const fastjet::PseudoJet& pj, pseudojets()) {
assert(clusterSeq() != NULL);
const PseudoJets parts = clusterSeq()->constituents(pj);
vector<Particle> constituents, tags;
constituents.reserve(parts.size());
foreach (const fastjet::PseudoJet& p, parts) {
map<int, Particle>::const_iterator found = _particles.find(p.user_index());
assert(found != _particles.end());
assert(found->first != 0);
if (found->first > 0) constituents.push_back(found->second);
else if (found->first < 0) tags.push_back(found->second);
}
Jet j(pj, constituents, tags);
rtn.push_back(j);
}
/// @todo Cache?
return rtn;
}
Add a colliding beam pair. Definition at line 107 of file Projection.hh. References Projection::_beamPairs. Referenced by Projection::Projection(). {
_beamPairs.insert(PdgIdPair(beam1, beam2));
return *this;
}
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 Definition at line 116 of file ProjectionApplier.hh. References ProjectionApplier::_addProjection(). Referenced by TauFinder::_init(), FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), ATLAS_2010_S8894728::init(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMS_2011_S8950903::init(), CMS_2012_PAS_QCD_11_010::init(), CMS_2011_S8978280::init(), ATLAS_2014_I1298811::init(), ALICE_2011_S8909580::init(), CMS_2012_I1090423::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2010_CONF_2010_049::init(), ATLAS_2011_I894867::init(), LHCB_2013_I1208105::init(), CMS_2011_S8941262::init(), CMS_2011_S8968497::init(), CMS_2011_S8973270::init(), CDF_2007_S7057202::init(), ALICE_2011_S8945144::init(), TOTEM_2012_I1115294::init(), ATLAS_2014_I1282441::init(), CMS_2011_S9086218::init(), CMS_2011_S9215166::init(), CMS_2012_I1087342::init(), CMS_2012_I1184941::init(), CMS_2012_I1193338::init(), ATLAS_2010_S8591806::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), LHCF_2012_I1115479::init(), ATLAS_2011_I930220::init(), CMS_2011_I954992::init(), TOTEM_2012_002::init(), CMS_2011_S8957746::init(), CMS_2011_S9088458::init(), ATLAS_2011_S9002537::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), CDF_2012_NOTE10874::init(), STAR_2008_S7993412::init(), UA5_1987_S1640666::init(), CMS_QCD_10_024::init(), CDF_1997_S3541940::init(), ATLAS_2012_I1124167::init(), STAR_2006_S6500200::init(), ATLAS_2011_S9128077::init(), BELLE_2008_I786560::init(), CDF_1993_S2742446::init(), ARGUS_1993_S2789213::init(), CDF_2000_S4155203::init(), CDF_2005_S6080774::init(), CDF_2006_S6450792::init(), SFM_1984_S1178091::init(), BABAR_2007_S7266081::init(), CMS_2013_I1261026::init(), D0_2008_S6879055::init(), ATLAS_2012_I1091481::init(), UA5_1982_S875503::init(), CDF_2008_S7782535::init(), CDF_2008_S8093652::init(), ATLAS_2012_I1204447::init(), STAR_2006_S6870392::init(), ARGUS_1993_S2669951::init(), UA5_1989_S1926373::init(), CMS_2013_I1272853::init(), D0_2000_I499943::init(), ALEPH_1991_S2435284::init(), E735_1998_S3905616::init(), ATLAS_2012_I1183818::init(), ATLAS_2013_I1243871::init(), ATLAS_2014_I1268975::init(), CMS_2012_I1102908::init(), ATLAS_2011_S9035664::init(), ALEPH_2002_S4823664::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), CDF_1988_S1865951::init(), CDF_1990_S2089246::init(), DELPHI_1999_S3960137::init(), ATLAS_2010_S8817804::init(), ATLAS_2011_S8924791::init(), ALEPH_1996_S3196992::init(), CDF_2005_S6217184::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), UA5_1986_S1583476::init(), D0_2001_S4674421::init(), D0_2008_S7554427::init(), CDF_1994_S2952106::init(), H1_1995_S3167097::init(), CDF_2008_S7540469::init(), CDF_2008_S7828950::init(), ATLAS_2014_I1307756::init(), UA1_1990_S2044935::init(), ZEUS_2001_S4815815::init(), LHCB_2011_I919315::init(), ATLAS_2012_I1119557::init(), ATLAS_2012_I1188891::init(), CDF_2008_S8095620::init(), DELPHI_1995_S3137023::init(), STAR_2006_S6860818::init(), JADE_1998_S3612880::init(), JADE_OPAL_2000_S4300807::init(), ALICE_2010_S8624100::init(), OPAL_1997_S3396100::init(), ALEPH_2004_S5765862::init(), EXAMPLE::init(), D0_1996_S3214044::init(), STAR_2009_UE_HELEN::init(), ATLAS_2012_I1125575::init(), ALICE_2010_S8625980::init(), LHCB_2013_I1218996::init(), CMS_2013_I1218372::init(), ALICE_2010_S8706239::init(), CDF_2009_NOTE_9936::init(), D0_2009_S8320160::init(), CDF_1998_S3618439::init(), CDF_2000_S4266730::init(), CDF_2001_S4563131::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), CDF_2009_S8436959::init(), ATLAS_2013_I1263495::init(), CDF_1996_S3108457::init(), D0_2010_S8566488::init(), D0_2006_S6438750::init(), CDF_1996_S3349578::init(), DELPHI_2002_069_CONF_603::init(), ATLAS_2012_I1082936::init(), ALEPH_2001_S4656318::init(), CDF_2001_S4517016::init(), SLD_2002_S4869273::init(), LHCB_2012_I1119400::init(), ATLAS_2011_I921594::init(), CDF_1996_S3418421::init(), UA5_1988_S1867512::init(), LHCB_2011_I917009::init(), SLD_1996_S3398250::init(), CDF_2008_S7541902::init(), ATLAS_2011_S9108483::init(), CDF_2006_S6653332::init(), D0_2008_S7662670::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I1082009::init(), ATLAS_2012_I946427::init(), OPAL_2002_S5361494::init(), DELPHI_2000_S4328825::init(), ATLAS_2012_I1199269::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_I1084540::init(), CDF_2009_S8233977::init(), CDF_2010_S8591881_QCD::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), CDF_2001_S4751469::init(), D0_2008_S7719523::init(), ATLAS_2011_S9120807::init(), ATLAS_2010_S8918562::init(), ATLAS_2012_I1094568::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), ATLAS_2011_S9126244::init(), BABAR_2007_S6895344::init(), BELLE_2001_S4598261::init(), BABAR_2005_S6181155::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ATLAS_2012_I1094061::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), ATLAS_2012_I1093734::init(), OPAL_1993_S2692198::init(), ATLAS_2014_I1279489::init(), ATLAS_2012_I1094564::init(), SLD_2004_S5693039::init(), BELLE_2006_S6265367::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder(). {
const Projection& reg = _addProjection(proj, name);
const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
return rtn;
}
Apply the supplied projection on event. Definition at line 70 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). Referenced by DISFinalState::project(). {
return pcast<PROJ>(_applyProjection(evt, proj));
}
Apply the supplied projection on event. Definition at line 77 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). {
return pcast<PROJ>(_applyProjection(evt, proj));
}
Apply the named projection on event. Definition at line 84 of file ProjectionApplier.hh. References ProjectionApplier::_applyProjection(). {
return pcast<PROJ>(_applyProjection(evt, name));
}
Return the area definition.. May be null. Definition at line 181 of file FastJets.hh. References FastJets::_adef. Referenced by FastJets::clusterSeqArea(). {
return _adef;
}
Return the allowed beam pairs on which this projection can operate, not including recursion. Derived classes should ensure that all contained projections are registered in the _projections set for the beam constraint chaining to work.
Definition at line 33 of file Projection.cc. References Projection::_beamPairs, Projection::beamPairs(), Projection::getLog(), ProjectionApplier::getProjections(), Rivet::intersection(), and Log::TRACE. Referenced by Projection::beamPairs(). {
set<PdgIdPair> ret = _beamPairs;
set<ConstProjectionPtr> projs = getProjections();
for (set<ConstProjectionPtr>::const_iterator ip = projs.begin(); ip != projs.end(); ++ip) {
ConstProjectionPtr p = *ip;
getLog() << Log::TRACE << "Proj addr = " << p << endl;
if (p) ret = intersection(ret, p->beamPairs());
}
return ret;
}
Determine whether this object should be ordered before the object p given as argument. If p is of a different class than this, the before() function of the corresponding type_info objects is used. Otherwise, if the objects are of the same class, the virtual compare(const Projection &) will be returned. Definition at line 22 of file Projection.cc. References Projection::compare(). Referenced by less< const Rivet::Projection * >::operator()(). {
const std::type_info& thisid = typeid(*this);
const std::type_info& otherid = typeid(p);
if (thisid == otherid) {
return compare(p) < 0;
} else {
return thisid.before(otherid);
}
}
Do the calculation locally (no caching).
Implements JetAlg. Definition at line 104 of file FastJets.cc. References FastJets::_adef, FastJets::_cseq, FastJets::_jdef, FastJets::_particles, FourMomentum::E(), Particle::momentum(), MSG_DEBUG, FourMomentum::px(), FourMomentum::py(), and FourMomentum::pz(). Referenced by CDF_2008_S7540469::analyze(), and FastJets::project(). {
_particles.clear();
vector<fastjet::PseudoJet> pjs;
/// @todo Use FastJet3's UserInfo system
// Store 4 vector data about each particle into FastJet's PseudoJets
int counter = 1;
foreach (const Particle& p, fsparticles) {
const FourMomentum fv = p.momentum();
fastjet::PseudoJet pj(fv.px(), fv.py(), fv.pz(), fv.E()); //< @todo Eliminate?
pj.set_user_index(counter);
pjs.push_back(pj);
_particles[counter] = p;
counter += 1;
}
// And the same for ghost tagging particles (with negative user indices)
counter = 1;
foreach (const Particle& p, tagparticles) {
const FourMomentum fv = 1e-20 * p.momentum(); //< Ghostify the momentum
fastjet::PseudoJet pj(fv.px(), fv.py(), fv.pz(), fv.E()); //< @todo Eliminate?
pj.set_user_index(-counter);
pjs.push_back(pj);
_particles[-counter] = p;
counter += 1;
}
MSG_DEBUG("Running FastJet ClusterSequence construction");
// Choose cseq as basic or area-calculating
if (_adef == NULL) {
_cseq.reset(new fastjet::ClusterSequence(pjs, _jdef));
} else {
_cseq.reset(new fastjet::ClusterSequenceArea(pjs, _jdef, *_adef));
}
}
Clone on the heap. Implements JetAlg. Definition at line 97 of file FastJets.hh. References FastJets::FastJets(). {
return new FastJets(*this);
}
Return the cluster sequence. Definition at line 165 of file FastJets.hh. References FastJets::_cseq. Referenced by FastJets::_jets(), MC_JetSplittings::analyze(), ALEPH_1996_S3196992::analyze(), JADE_OPAL_2000_S4300807::analyze(), JADE_1998_S3612880::analyze(), OPAL_2001_S4553896::analyze(), OPAL_2004_S6132243::analyze(), DELPHI_2003_WUD_03_11::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), FastJets::filterJet(), FastJets::pseudoJets(), ATLAS_2012_I1094564::splitjet(), FastJets::splitJet(), and FastJets::ySubJet(). {
return _cseq.get();
}
Return the area-enabled cluster sequence (if an area defn exists, otherwise returns a null ptr). Definition at line 170 of file FastJets.hh. References FastJets::_cseq, and FastJets::areaDef().
Compare projections. Implements JetAlg. Definition at line 78 of file FastJets.cc. References FastJets::_adef, FastJets::_jdef, JetAlg::_useInvisibles, Rivet::cmp(), and Projection::mkNamedPCmp(). {
const FastJets& other = dynamic_cast<const FastJets&>(p);
return \
cmp(_useInvisibles, other._useInvisibles) ||
(_useInvisibles ? mkNamedPCmp(other, "FS") : mkNamedPCmp(other, "VFS")) ||
cmp(_jdef.jet_algorithm(), other._jdef.jet_algorithm()) ||
cmp(_jdef.recombination_scheme(), other._jdef.recombination_scheme()) ||
cmp(_jdef.plugin(), other._jdef.plugin()) ||
cmp(_jdef.R(), other._jdef.R()) ||
cmp(_adef, other._adef);
}
Determine if the jet collection is empty. Definition at line 125 of file JetAlg.hh. References JetAlg::size(). { return size() != 0; }
Template-usable interface common to FinalState. Definition at line 134 of file JetAlg.hh. References JetAlg::jets(). { return jets(); }
Filter a jet a la PRL100,242001(2008). Based on code from G.Salam, A.Davison. Definition at line 239 of file FastJets.cc. References FastJets::clusterSeq(). {
assert(clusterSeq());
if (jet.E() <= 0.0 || clusterSeq()->constituents(jet).size() == 0) {
return jet;
}
if (stingy_R == 0.0) {
stingy_R = def_R;
}
stingy_R = def_R < stingy_R ? def_R : stingy_R;
fastjet::JetDefinition stingy_jet_def(fastjet::cambridge_algorithm, stingy_R);
//FlavourRecombiner recom;
//stingy_jet_def.set_recombiner(&recom);
fastjet::ClusterSequence scs(clusterSeq()->constituents(jet), stingy_jet_def);
std::vector<fastjet::PseudoJet> stingy_jets = sorted_by_pt(scs.inclusive_jets());
fastjet::PseudoJet reconst_jet(0.0, 0.0, 0.0, 0.0);
for (unsigned isj = 0; isj < std::min(3U, (unsigned int) stingy_jets.size()); ++isj) {
reconst_jet += stingy_jets[isj];
}
return reconst_jet;
}
Get a Log object based on the getName() property of the calling projection object. Reimplemented from ProjectionApplier. Definition at line 114 of file Projection.hh. References Projection::name(). Referenced by Projection::beamPairs(), InvMassFinalState::calc(), ChargedFinalState::project(), InitialQuarks::project(), UnstableFinalState::project(), PromptFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), and VetoedFinalState::project(). {
string logname = "Rivet.Projection." + name();
return Log::getLog(logname);
}
Get the named projection, specifying return type via a template argument. Definition at line 52 of file ProjectionApplier.hh. References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler(). Referenced by ProjectionApplier::_applyProjection(), Rivet::pcmp(), and Hemispheres::project(). {
const Projection& p = getProjHandler().getProjection(*this, name);
return pcast<PROJ>(p);
}
Get the named projection (non-templated, so returns as a reference to a Projection base class). Definition at line 60 of file ProjectionApplier.hh. References ProjectionHandler::getProjection(), and ProjectionApplier::getProjHandler(). {
return getProjHandler().getProjection(*this, name);
}
Get the contained projections, including recursion. Definition at line 45 of file ProjectionApplier.hh. References ProjectionHandler::DEEP, ProjectionHandler::getChildProjections(), and ProjectionApplier::getProjHandler(). Referenced by Projection::beamPairs(). {
return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
}
Get a reference to the ProjectionHandler for this thread. Definition at line 97 of file ProjectionApplier.hh. References ProjectionApplier::_projhandler. Referenced by ProjectionApplier::_addProjection(), ProjectionApplier::getProjection(), ProjectionApplier::getProjections(), and ProjectionApplier::~ProjectionApplier(). {
return _projhandler;
}
Return the jet definition. Definition at line 176 of file FastJets.hh. References FastJets::_jdef. {
return _jdef;
}
Get jets in no guaranteed order, with optional cuts on
Definition at line 42 of file JetAlg.hh. References JetAlg::_jets(), FourVector::eta(), Rivet::GeV, Rivet::inRange(), Jet::momentum(), MSG_DEBUG, Rivet::PSEUDORAPIDITY, FourMomentum::pT(), Rivet::RAPIDITY, and FourMomentum::rapidity(). Referenced by CMS_2012_I1087342::analyze(), CMS_2011_S9086218::analyze(), CMS_2012_I1102908::analyze(), ZEUS_2001_S4815815::analyze(), CDF_2008_S7540469::analyze(), D0_2008_S7662670::analyze(), JetAlg::entities(), JetAlg::jets(), JetAlg::jetsByE(), JetAlg::jetsByEt(), JetAlg::jetsByP(), and JetAlg::jetsByPt(). {
const Jets rawjets = _jets(ptmin);
Jets rtn;
MSG_DEBUG("Raw jet size (with pTmin cut = " << ptmin/GeV << " GeV) = " << rawjets.size());
foreach (const Jet& j, rawjets) {
const FourMomentum pj = j.momentum();
if (!inRange(pj.pT(), ptmin, ptmax)) continue;
if (rapscheme == PSEUDORAPIDITY && !inRange(pj.eta(), rapmin, rapmax)) continue;
if (rapscheme == RAPIDITY && !inRange(pj.rapidity(), rapmin, rapmax)) continue;
rtn += j;
}
return rtn;
}
Get the jets, ordered by supplied sorting function object, with optional cuts on
Definition at line 62 of file JetAlg.hh. References JetAlg::jets().
Get the jets, ordered by
Definition at line 94 of file JetAlg.hh. References Rivet::cmpMomByE(), and JetAlg::jets().
Get the jets, ordered by
Definition at line 103 of file JetAlg.hh. References Rivet::cmpMomByEt(), and JetAlg::jets(). Referenced by CDF_2008_S7541902::analyze(). {
return jets(cmpMomByEt, ptmin, ptmax, rapmin, rapmax, rapscheme);
}
Get the jets, ordered by
Definition at line 85 of file JetAlg.hh. References Rivet::cmpMomByP(), and JetAlg::jets().
Get the jets, ordered by
Definition at line 76 of file JetAlg.hh. References Rivet::cmpMomByPt(), and JetAlg::jets(). Referenced by CDF_2008_S8093652::analyze(), STAR_2006_S6870392::analyze(), ATLAS_2010_CONF_2010_049::analyze(), CMS_2011_S9215166::analyze(), CMS_2013_I1272853::analyze(), MC_DIJET::analyze(), MC_LEADJETUE::analyze(), D0_2008_S7863608::analyze(), CMS_2013_I1209721::analyze(), D0_2004_S5992206::analyze(), CMS_2013_I1218372::analyze(), ATLAS_2012_I1183818::analyze(), CMS_2013_I1261026::analyze(), MC_TTBAR::analyze(), CDF_2001_S4751469::analyze(), ATLAS_2012_I1083318::analyze(), MC_SUSY::analyze(), CDF_2009_S8057893::CDF_2009_S8057893::analyze(), CMS_2013_I1261026::eventDecomp(), CMS_2013_I1258128::makePhotonCut(), CMS_2013_I1258128::makeZCut(), and ATLAS_2011_I945498::selectJets(). {
return jets(cmpMomByPt, ptmin, ptmax, rapmin, rapmax, rapscheme);
}
Shortcut to make a named Cmp<Projection> comparison with the Definition at line 45 of file Projection.cc. References Rivet::pcmp(). Referenced by ChargedFinalState::compare(), ChargedLeptons::compare(), BeamThrust::compare(), FParameter::compare(), CentralEtHCM::compare(), MergedFinalState::compare(), DISLepton::compare(), DISKinematics::compare(), NeutralFinalState::compare(), DISFinalState::compare(), VisibleFinalState::compare(), PromptFinalState::compare(), Spherocity::compare(), LeadingParticlesFinalState::compare(), ParisiTensor::compare(), FoxWolframMoments::compare(), LossyFinalState< ConstRandomFilter >::compare(), MissingMomentum::compare(), Thrust::compare(), InvMassFinalState::compare(), Sphericity::compare(), DressedLeptons::compare(), ZFinder::compare(), Hemispheres::compare(), WFinder::compare(), IdentifiedFinalState::compare(), VetoedFinalState::compare(), JetShape::compare(), and FastJets::compare(). {
return pcmp(*this, otherparent, pname);
}
Shortcut to make a named Cmp<Projection> comparison with the Definition at line 51 of file Projection.cc. References Rivet::pcmp(). {
return pcmp(*this, otherparent, pname);
}
Get the name of the projection. Implements ProjectionApplier. Definition at line 101 of file Projection.hh. References Projection::_name. Referenced by ProjectionApplier::_addProjection(), ProjectionHandler::_checkDuplicate(), ProjectionHandler::_clone(), ProjectionHandler::_getEquiv(), VetoedFinalState::addVetoOnThisFinalState(), Projection::getLog(), ProjectionHandler::registerProjection(), and Projection::setName(). {
return _name;
}
Number of jets above the Definition at line 148 of file FastJets.cc. References FastJets::_cseq. Referenced by FastJets::size(). Perform the projection on the Event. Implements JetAlg. Definition at line 92 of file FastJets.cc. References JetAlg::_useInvisibles, FastJets::calc(), and Rivet::particles(). {
// Final state particles
const string fskey = _useInvisibles ? "FS" : "VFS";
const Particles fsparticles = applyProjection<FinalState>(e, fskey).particles();
// Tagging particles
const Particles chadrons = applyProjection<HeavyHadrons>(e, "HFHadrons").cHadrons();
const Particles bhadrons = applyProjection<HeavyHadrons>(e, "HFHadrons").bHadrons();
const Particles taus = applyProjection<FinalState>(e, "Taus").particles();
calc(fsparticles, chadrons+bhadrons+taus);
}
Get the pseudo jets (unordered). Definition at line 176 of file FastJets.cc. References FastJets::clusterSeq(). Referenced by FastJets::pseudojets(), FastJets::pseudoJetsByE(), FastJets::pseudoJetsByPt(), and FastJets::pseudoJetsByRapidity(). {
return (clusterSeq() != NULL) ? clusterSeq()->inclusive_jets(ptmin) : PseudoJets();
}
Alias. Definition at line 135 of file FastJets.hh. References FastJets::pseudoJets(). Referenced by FastJets::_jets(). { return pseudoJets(ptmin); }
Get the pseudo jets, ordered by Definition at line 145 of file FastJets.hh. References FastJets::pseudoJets(). Referenced by FastJets::pseudojetsByE(). {
return sorted_by_E(pseudoJets(ptmin));
}
Alias. Definition at line 149 of file FastJets.hh. References FastJets::pseudoJetsByE(). { return pseudoJetsByE(ptmin); }
Get the pseudo jets, ordered by Definition at line 138 of file FastJets.hh. References FastJets::pseudoJets(). Referenced by CDF_2008_S8095620::analyze(), CDF_2006_S6653332::analyze(), ATLAS_2012_I1094564::analyze(), and FastJets::pseudojetsByPt(). {
return sorted_by_pt(pseudoJets(ptmin));
}
Alias. Definition at line 142 of file FastJets.hh. References FastJets::pseudoJetsByPt(). { return pseudoJetsByPt(ptmin); }
Get the pseudo jets, ordered by rapidity. Definition at line 152 of file FastJets.hh. References FastJets::pseudoJets(). Referenced by FastJets::pseudojetsByRapidity(). {
return sorted_by_rapidity(pseudoJets(ptmin));
}
Alias. Definition at line 156 of file FastJets.hh. References FastJets::pseudoJetsByRapidity(). { return pseudoJetsByRapidity(ptmin); }
Reset the projection. Jet def, etc. are unchanged.
Implements JetAlg. Definition at line 141 of file FastJets.cc. References FastJets::_particles, and FastJets::_yscales. {
_yscales.clear();
_particles.clear();
/// @todo _cseq = fastjet::ClusterSequence();
}
Used by derived classes to set their name. Definition at line 120 of file Projection.hh. References Projection::_name, and Projection::name(). Referenced by TauFinder::_init(), FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), Beam::Beam(), BeamThrust::BeamThrust(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ConstLossyFinalState::ConstLossyFinalState(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), InitialQuarks::InitialQuarks(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), UnstableFinalState::UnstableFinalState(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder().
Number of jets. Implements JetAlg. Definition at line 125 of file FastJets.hh. References FastJets::numJets(). Referenced by CDF_2008_S8095620::analyze(), and CDF_2006_S6653332::analyze(). {
return numJets();
}
Split a jet a la PRL100,242001(2008). Based on code from G.Salam, A.Davison. Definition at line 200 of file FastJets.cc. References FastJets::_cseq, FastJets::_jdef, FastJets::clusterSeq(), and MSG_DEBUG. {
// Sanity cuts
if (jet.E() <= 0 || _cseq->constituents(jet).size() <= 1) {
return jet;
}
// Build a new cluster sequence just using the consituents of this jet.
assert(clusterSeq());
fastjet::ClusterSequence cs(clusterSeq()->constituents(jet), _jdef);
// Get the jet back again
fastjet::PseudoJet remadeJet = cs.inclusive_jets()[0];
MSG_DEBUG("Jet2:" << remadeJet.m() << "," << remadeJet.e());
fastjet::PseudoJet parent1, parent2;
fastjet::PseudoJet split(0.0, 0.0, 0.0, 0.0);
while (cs.has_parents(remadeJet, parent1, parent2)) {
MSG_DEBUG("Parents:" << parent1.m() << "," << parent2.m());
if (parent1.m2() < parent2.m2()) {
fastjet::PseudoJet tmp;
tmp = parent1; parent1 = parent2; parent2 = tmp;
}
double ktdist = parent1.kt_distance(parent2);
double rtycut2 = 0.3*0.3;
if (parent1.m() < ((2.0*remadeJet.m())/3.0) && ktdist > rtycut2*remadeJet.m2()) {
break;
} else {
remadeJet = parent1;
}
}
last_R = 0.5 * sqrt(parent1.squared_distance(parent2));
split.reset(remadeJet.px(), remadeJet.py(), remadeJet.pz(), remadeJet.E());
return split;
}
Include invisible particles in jet construction. The default behaviour is that jets are only constructed from visible (i.e. charged under an SM gauge group) particles. Some jet studies, including those from ATLAS, use a definition in which neutrinos from hadron decays are included (via MC correction) in the experimental jet definition. Setting this flag to true avoids the automatic restriction to a VisibleFinalState. Definition at line 35 of file JetAlg.hh. References JetAlg::_useInvisibles. Referenced by CMS_2011_S8973270::init(), D0_2011_I895662::init(), ATLAS_2011_I930220::init(), ATLAS_2011_S9128077::init(), CDF_2008_S7782535::init(), D0_2000_I499943::init(), ATLAS_2014_I1268975::init(), ATLAS_2011_S8924791::init(), CDF_2005_S6217184::init(), ALEPH_2004_S5765862::init(), JADE_OPAL_2000_S4300807::init(), ATLAS_2011_I945498::init(), MC_HFJETS::init(), ATLAS_2012_I1082936::init(), ATLAS_2014_I1304688::init(), ATLAS_2013_I1230812::init(), ATLAS_2012_I1083318::init(), ATLAS_2013_I1217867::init(), and ATLAS_2012_I1093738::init(). {
_useInvisibles = useinvis;
}
Use provided jet area definition.
Provide an adef null pointer to re-disable jet area calculation Definition at line 114 of file FastJets.hh. References FastJets::_adef. Referenced by ATLAS_2014_I1307756::init(), ATLAS_2013_I1263495::init(), ATLAS_2011_I921594::init(), ATLAS_2012_I1199269::init(), ATLAS_2010_S8914702::init(), ATLAS_2011_S9120807::init(), and ATLAS_2012_I1093738::init(). {
_adef = adef;
}
Get the subjet splitting variables for the given jet. Definition at line 184 of file FastJets.cc. References FastJets::_jdef, FastJets::_yscales, and FastJets::clusterSeq(). {
assert(clusterSeq());
fastjet::ClusterSequence subjet_cseq(clusterSeq()->constituents(jet), _jdef);
vector<double> yMergeVals;
for (int i = 1; i < 4; ++i) {
// Multiply the dmerge value by R^2 so that it corresponds to a
// relative k_T (fastjet has 1/R^2 in the d_ij distance by default)
const double ktmerge = subjet_cseq.exclusive_dmerge(i) * _jdef.R()*_jdef.R();
yMergeVals.push_back(ktmerge/jet.perp2());
}
_yscales.insert(make_pair( jet.cluster_hist_index(), yMergeVals ));
return yMergeVals;
}
Friends And Related Function Documentation
The Cmp specialization for Projection is a friend. Definition at line 35 of file Projection.hh.
Event is a friend. Definition at line 32 of file Projection.hh. Member Data Documentation
Pointer to user-handled area definition. Definition at line 229 of file FastJets.hh. Referenced by FastJets::areaDef(), FastJets::calc(), FastJets::compare(), and FastJets::useJetArea().
Flag to forbid projection registration in analyses until the init phase. Definition at line 143 of file ProjectionApplier.hh. Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().
Cluster sequence. Definition at line 232 of file FastJets.hh. Referenced by FastJets::calc(), FastJets::clusterSeq(), FastJets::clusterSeqArea(), FastJets::numJets(), and FastJets::splitJet().
Jet definition. Definition at line 226 of file FastJets.hh. Referenced by FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), FastJets::calc(), FastJets::compare(), FastJets::jetDef(), FastJets::splitJet(), and FastJets::ySubJet().
set of particles sorted by their PT2 Definition at line 242 of file FastJets.hh. Referenced by FastJets::_jets(), FastJets::calc(), and FastJets::reset().
FastJet external plugin. Definition at line 235 of file FastJets.hh. Referenced by FastJets::_init1(), and FastJets::_init3().
Flag to determine whether or not the VFS wrapper is to be used. Definition at line 152 of file JetAlg.hh. Referenced by FastJets::compare(), FastJets::project(), and JetAlg::useInvisibles().
Map of vectors of y scales. This is mutable so we can use caching/lazy evaluation. Definition at line 238 of file FastJets.hh. Referenced by FastJets::reset(), and FastJets::ySubJet(). The documentation for this class was generated from the following files: Generated on Tue Sep 30 2014 19:46:00 for The Rivet MC analysis system by 
1.7.6.1
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
cut. 
. 
. 
, with optional cuts on 
, with optional cuts on