|
|
|
Base class for all types of Jets. More...
#include <Jet.h>
Classes | |
| class | EtaPhiMoments |
| record to store eta-phi first and second moments More... | |
Public Types | |
| typedef edm::Ptr< Candidate > | Constituent |
| typedef std::vector< Constituent > | Constituents |
Public Member Functions | |
| float | constituentEtaPhiSpread () const |
| float | constituentPtDistribution () const |
| float | etaetaMoment () const |
| eta-eta second moment, ET weighted | |
| float | etaphiMoment () const |
| eta-phi second moment, ET weighted | |
| EtaPhiMoments | etaPhiStatistics () const |
| eta-phi statistics, ET weighted | |
| float | etInAnnulus (float fRmin, float fRmax) const |
| ET in annulus between rmin and rmax around jet direction. | |
| virtual Constituents | getJetConstituents () const |
| list of constituents | |
| virtual std::vector< const reco::Candidate * > | getJetConstituentsQuick () const |
| quick list of constituents | |
| bool | isJet () const |
| Jet () | |
| Default constructor. | |
| Jet (const LorentzVector &fP4, const Point &fVertex) | |
| Initiator. | |
| Jet (const LorentzVector &fP4, const Point &fVertex, const Constituents &fConstituents) | |
| virtual float | jetArea () const |
| get jet area | |
| float | maxDistance () const |
| maximum distance from jet to constituent | |
| int | nCarrying (float fFraction) const |
| return # of constituent carrying fraction of energy | |
| virtual int | nConstituents () const |
| # of constituents | |
| virtual int | nPasses () const |
| number of passes taken by algorithm | |
| float | phiphiMoment () const |
| phi-phi second moment, ET weighted | |
| virtual float | pileup () const |
| pileup energy contribution as calculated by algorithm | |
| virtual std::string | print () const |
| Print object. | |
| virtual void | scaleEnergy (double fScale) |
| scale energy of the jet | |
| virtual void | setJetArea (float fArea) |
| set jet area | |
| virtual void | setNPasses (int fPasses) |
| Set number of passes taken by algorithm. | |
| virtual void | setPileup (float fEnergy) |
| Set pileup energy contribution as calculated by algorithm. | |
| virtual | ~Jet () |
| Destructor. | |
Static Public Member Functions | |
| static float | detectorEta (float fZVertex, float fPhysicsEta) |
| static function to convert physics eta to detector eta | |
| static Candidate::LorentzVector | detectorP4 (const Candidate::Point &vertex, const Candidate &inParticle) |
| static float | physicsEta (float fZVertex, float fDetectorEta) |
| static function to convert detector eta to physics eta | |
| static Candidate::LorentzVector | physicsP4 (const Candidate::Point &newVertex, const Candidate &inParticle, const Candidate::Point &oldVertex=Candidate::Point(0, 0, 0)) |
Private Attributes | |
| float | mJetArea |
| int | mPassNumber |
| float | mPileupEnergy |
Base class for all types of Jets.
GenericJet describes jets made from arbitrary constituents, No direct constituents reference is stored for now
Jet describes properties common for all kinds of jets, essentially kinematics. Base class for all types of Jets.
| typedef edm::Ptr<Candidate> reco::Jet::Constituent |
| typedef std::vector<Constituent> reco::Jet::Constituents |
| reco::Jet::Jet | ( | ) | [inline] |
Default constructor.
Reimplemented in pat::Jet.
Definition at line 37 of file Jet.h.
: mJetArea (0), mPileupEnergy (0), mPassNumber (0) {}
| Jet::Jet | ( | const LorentzVector & | fP4, |
| const Point & | fVertex | ||
| ) |
Initiator.
Definition at line 175 of file Jet.cc.
: CompositePtrCandidate (0, fP4, fVertex), mJetArea (0), mPileupEnergy (0), mPassNumber (0) {}
| Jet::Jet | ( | const LorentzVector & | fP4, |
| const Point & | fVertex, | ||
| const Constituents & | fConstituents | ||
| ) |
Definition at line 162 of file Jet.cc.
References reco::CompositePtrCandidate::addDaughter(), and i.
: CompositePtrCandidate (0, fP4, fVertex), mJetArea (0), mPileupEnergy (0), mPassNumber (0) { for (unsigned i = 0; i < fConstituents.size (); i++) { addDaughter (fConstituents [i]); } }
| virtual reco::Jet::~Jet | ( | ) | [inline, virtual] |
| float Jet::constituentEtaPhiSpread | ( | ) | const |
Definition at line 395 of file Jet.cc.
References reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, getJetConstituents(), and reco::LeafCandidate::pt().
{
Jet::Constituents constituents = this->getJetConstituents();
float sum_pt2 = 0.;
float sum_pt2deltaR2 = 0.;
for( unsigned iConst=0; iConst<constituents.size(); ++iConst ) {
LorentzVector thisConstituent = constituents[iConst]->p4();
float pt = thisConstituent.Pt();
float pt2 = pt*pt;
double dR = deltaR (*this, *(constituents[iConst]));
float pt2deltaR2 = pt*pt*dR*dR;
sum_pt2 += pt2;
sum_pt2deltaR2 += pt2deltaR2;
} //for constituents
float rmsCand_value = (sum_pt2>0.) ? sum_pt2deltaR2/sum_pt2 : 0.;
return rmsCand_value;
} //constituentEtaPhiSpread
| float Jet::constituentPtDistribution | ( | ) | const |
Definition at line 370 of file Jet.cc.
References getJetConstituents(), reco::LeafCandidate::pt(), and mathSSE::sqrt().
{
Jet::Constituents constituents = this->getJetConstituents();
float sum_pt2 = 0.;
float sum_pt = 0.;
for( unsigned iConst=0; iConst<constituents.size(); ++iConst ) {
float pt = constituents[iConst]->p4().Pt();
float pt2 = pt*pt;
sum_pt += pt;
sum_pt2 += pt2;
} //for constituents
float ptD_value = (sum_pt>0.) ? sqrt( sum_pt2 / (sum_pt*sum_pt) ) : 0.;
return ptD_value;
} //constituentPtDistribution
| float Jet::detectorEta | ( | float | fZVertex, |
| float | fPhysicsEta | ||
| ) | [static] |
| Candidate::LorentzVector Jet::detectorP4 | ( | const Candidate::Point & | vertex, |
| const Candidate & | inParticle | ||
| ) | [static] |
Definition at line 340 of file Jet.cc.
References reco::Candidate::energy(), reco::Candidate::momentum(), np, reco::LeafCandidate::p(), and p3.
{
CaloPoint3D<Point> caloPoint(vertex,inParticle.momentum()); // Jet position in Calo.
static const Point np(0,0,0);
Vector detectorDir = caloPoint.caloPoint() - np;
double p = inParticle.momentum().r();
Vector p3 = p * detectorDir.unit();
LorentzVector returnVector(p3.x(), p3.y(), p3.z(), inParticle.energy());
return returnVector;
}
| float Jet::etaetaMoment | ( | ) | const |
eta-eta second moment, ET weighted
eta-eta second moment
Definition at line 224 of file Jet.cc.
References getJetConstituentsQuick(), i, relativeConstraints::value, and CommonMethods::weight().
Referenced by fireworks::makeEveJetCone(), and JetIDSelectionFunctor::operator()().
{
std::vector<const Candidate*> towers = getJetConstituentsQuick ();
double sumw = 0;
double sum = 0;
double sum2 = 0;
int i = towers.size();
while (--i >= 0) {
double value = towers[i]->eta();
double weight = towers[i]->et();
sumw += weight;
sum += value * weight;
sum2 += value * value * weight;
}
return sumw > 0 ? (sum2 - sum*sum/sumw ) / sumw : 0;
}
| float Jet::etaphiMoment | ( | ) | const |
eta-phi second moment, ET weighted
eta-phi second moment
Definition at line 259 of file Jet.cc.
References reco::deltaPhi(), getJetConstituentsQuick(), i, reco::LeafCandidate::phi(), and CommonMethods::weight().
{
double phiRef = phi();
std::vector<const Candidate*> towers = getJetConstituentsQuick ();
double sumw = 0;
double sumA = 0;
double sumB = 0;
double sumAB = 0;
int i = towers.size();
while (--i >= 0) {
double valueA = towers[i]->eta();
double valueB = deltaPhi (towers[i]->phi(), phiRef);
double weight = towers[i]->et();
sumw += weight;
sumA += valueA * weight;
sumB += valueB * weight;
sumAB += valueA * valueB * weight;
}
return sumw > 0 ? (sumAB - sumA*sumB/sumw ) / sumw : 0;
}
| Jet::EtaPhiMoments Jet::etaPhiStatistics | ( | ) | const |
eta-phi statistics, ET weighted
eta-phi statistics
Definition at line 184 of file Jet.cc.
References reco::deltaPhi(), reco::LeafCandidate::eta(), reco::Jet::EtaPhiMoments::etaEtaMoment, reco::Jet::EtaPhiMoments::etaMean, reco::Jet::EtaPhiMoments::etaPhiMoment, getJetConstituentsQuick(), i, reco::LeafCandidate::phi(), reco::Jet::EtaPhiMoments::phiMean, reco::Jet::EtaPhiMoments::phiPhiMoment, query::result, and CommonMethods::weight().
{
std::vector<const Candidate*> towers = getJetConstituentsQuick ();
double phiRef = phi();
double sumw = 0;
double sumEta = 0;
double sumEta2 = 0;
double sumPhi = 0;
double sumPhi2 = 0;
double sumEtaPhi = 0;
int i = towers.size();
while (--i >= 0) {
double eta = towers[i]->eta();
double phi = deltaPhi (towers[i]->phi(), phiRef);
double weight = towers[i]->et();
sumw += weight;
sumEta += eta * weight;
sumEta2 += eta * eta * weight;
sumPhi += phi * weight;
sumPhi2 += phi * phi * weight;
sumEtaPhi += eta * phi * weight;
}
Jet::EtaPhiMoments result;
if (sumw > 0) {
result.etaMean = sumEta / sumw;
result.phiMean = deltaPhi (phiRef + sumPhi, 0.);
result.etaEtaMoment = (sumEta2 - sumEta * sumEta / sumw) / sumw;
result.phiPhiMoment = (sumPhi2 - sumPhi * sumPhi / sumw) / sumw;
result.etaPhiMoment = (sumEtaPhi - sumEta * sumPhi / sumw) / sumw;
}
else {
result.etaMean = 0;
result.phiMean = 0;
result.etaEtaMoment = 0;
result.phiPhiMoment = 0;
result.etaPhiMoment = 0;
}
return result;
}
| float Jet::etInAnnulus | ( | float | fRmin, |
| float | fRmax | ||
| ) | const |
ET in annulus between rmin and rmax around jet direction.
et in annulus between rmin and rmax around jet direction
Definition at line 280 of file Jet.cc.
References reco::deltaR(), getJetConstituentsQuick(), i, and query::result.
| Jet::Constituents Jet::getJetConstituents | ( | ) | const [virtual] |
list of constituents
Definition at line 351 of file Jet.cc.
References reco::CompositePtrCandidate::daughterPtr(), i, reco::CompositePtrCandidate::numberOfDaughters(), and query::result.
Referenced by reco::helper::CastorJetIDHelper::calculate(), constituentEtaPhiSpread(), constituentPtDistribution(), fireworks::makeEveJetCone(), reco::tau::helpers::nCharged(), reco::tau::helpers::nGammas(), SmearedJetProducer_namespace::JetResolutionExtractorT< pat::Jet >::operator()(), CATopJetHelper::operator()(), JetPlusTrackProducerAA::produce(), JetPlusTrackProducer::produce(), and CastorClusterProducer::produce().
{
Jet::Constituents result;
for (unsigned i = 0; i < numberOfDaughters(); i++) {
result.push_back (daughterPtr (i));
}
return result;
}
| std::vector< const Candidate * > Jet::getJetConstituentsQuick | ( | ) | const [virtual] |
quick list of constituents
Definition at line 359 of file Jet.cc.
References reco::CompositePtrCandidate::daughter(), i, reco::CompositePtrCandidate::numberOfDaughters(), and query::result.
Referenced by FWPFTauProxyBuilder::buildViewType(), etaetaMoment(), etaphiMoment(), etaPhiStatistics(), etInAnnulus(), maxDistance(), nCarrying(), and phiphiMoment().
| bool Jet::isJet | ( | ) | const [virtual] |
| virtual float reco::Jet::jetArea | ( | ) | const [inline, virtual] |
get jet area
Definition at line 106 of file Jet.h.
References mJetArea.
Referenced by L1FastjetCorrector::correction(), and HiL1Subtractor::produce().
{return mJetArea;}
| float Jet::maxDistance | ( | ) | const |
maximum distance from jet to constituent
Definition at line 307 of file Jet.cc.
References reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, getJetConstituentsQuick(), i, and query::result.
| int Jet::nCarrying | ( | float | fFraction | ) | const |
return # of constituent carrying fraction of energy
return # of constituent carring fraction of energy. Assume ordered towers
Definition at line 292 of file Jet.cc.
References reco::LeafCandidate::et(), getJetConstituentsQuick(), i, and query::result.
Referenced by pat::Jet::n60(), reco::CaloJet::n60(), pat::Jet::n90(), and reco::CaloJet::n90().
{
std::vector<const Candidate*> towers = getJetConstituentsQuick ();
if (fFraction >= 1) return towers.size();
double totalEt = 0;
for (unsigned i = 0; i < towers.size(); ++i) totalEt += towers[i]->et();
double fractionEnergy = totalEt * fFraction;
unsigned result = 0;
for (; result < towers.size(); ++result) {
fractionEnergy -= towers[result]->et();
if (fractionEnergy <= 0) return result+1;
}
return 0;
}
| virtual int reco::Jet::nConstituents | ( | ) | const [inline, virtual] |
# of constituents
Definition at line 66 of file Jet.h.
References reco::CompositePtrCandidate::numberOfDaughters().
Referenced by print(), and JetIdSelector< T >::produce().
{return numberOfDaughters();}
| virtual int reco::Jet::nPasses | ( | ) | const [inline, virtual] |
number of passes taken by algorithm
Definition at line 116 of file Jet.h.
References mPassNumber.
{return mPassNumber;}
| float Jet::phiphiMoment | ( | ) | const |
phi-phi second moment, ET weighted
phi-phi second moment
Definition at line 241 of file Jet.cc.
References reco::deltaPhi(), getJetConstituentsQuick(), i, reco::LeafCandidate::phi(), relativeConstraints::value, and CommonMethods::weight().
Referenced by fireworks::makeEveJetCone(), and JetIDSelectionFunctor::operator()().
{
double phiRef = phi();
std::vector<const Candidate*> towers = getJetConstituentsQuick ();
double sumw = 0;
double sum = 0;
double sum2 = 0;
int i = towers.size();
while (--i >= 0) {
double value = deltaPhi (towers[i]->phi(), phiRef);
double weight = towers[i]->et();
sumw += weight;
sum += value * weight;
sum2 += value * value * weight;
}
return sumw > 0 ? (sum2 - sum*sum/sumw ) / sumw : 0;
}
| float Jet::physicsEta | ( | float | fZVertex, |
| float | fDetectorEta | ||
| ) | [static] |
static function to convert detector eta to physics eta
Definition at line 319 of file Jet.cc.
Referenced by JetTracksAssociationDRVertexAssigned::produce().
{
CaloPointZ refPoint (0., fDetectorEta);
return refPoint.etaReference (fZVertex);
}
| Candidate::LorentzVector Jet::physicsP4 | ( | const Candidate::Point & | newVertex, |
| const Candidate & | inParticle, | ||
| const Candidate::Point & | oldVertex = Candidate::Point(0,0,0) |
||
| ) | [static] |
Definition at line 331 of file Jet.cc.
References reco::Candidate::energy(), reco::Candidate::momentum(), reco::LeafCandidate::p(), and p3.
{
CaloPoint3D<Point> caloPoint(oldVertex,inParticle.momentum()); // Jet position in Calo.
Vector physicsDir = caloPoint.caloPoint() - newVertex;
double p = inParticle.momentum().r();
Vector p3 = p * physicsDir.unit();
LorentzVector returnVector(p3.x(), p3.y(), p3.z(), inParticle.energy());
return returnVector;
}
| virtual float reco::Jet::pileup | ( | ) | const [inline, virtual] |
pileup energy contribution as calculated by algorithm
Definition at line 111 of file Jet.h.
References mPileupEnergy.
{return mPileupEnergy;}
| std::string Jet::print | ( | void | ) | const [virtual] |
Print object.
Reimplemented in reco::BasicJet, reco::CaloJet, reco::GenJet, reco::JPTJet, reco::PFClusterJet, reco::PFJet, and reco::TrackJet.
Definition at line 425 of file Jet.cc.
References reco::CompositePtrCandidate::daughterPtr(), reco::LeafCandidate::eta(), edm::Ptr< T >::id(), edm::Ptr< T >::isNonnull(), edm::Ptr< T >::key(), nConstituents(), reco::CompositePtrCandidate::numberOfDaughters(), dbtoconf::out, reco::LeafCandidate::p(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::pz().
Referenced by reco::TrackJet::print(), and reco::PFClusterJet::print().
{
std::ostringstream out;
out << "Jet p/px/py/pz/pt: " << p() << '/' << px () << '/' << py() << '/' << pz() << '/' << pt() << std::endl
<< " eta/phi: " << eta () << '/' << phi () << std::endl
<< " # of constituents: " << nConstituents () << std::endl;
out << " Constituents:" << std::endl;
for (unsigned index = 0; index < numberOfDaughters(); index++) {
Constituent constituent = daughterPtr (index); // deref
if (constituent.isNonnull()) {
out << " #" << index << " p/pt/eta/phi: "
<< constituent->p() << '/' << constituent->pt() << '/' << constituent->eta() << '/' << constituent->phi()
<< " productId/index: " << constituent.id() << '/' << constituent.key() << std::endl;
}
else {
out << " #" << index << " constituent is not available in the event" << std::endl;
}
}
return out.str ();
}
| void Jet::scaleEnergy | ( | double | fScale | ) | [virtual] |
scale energy of the jet
Definition at line 445 of file Jet.cc.
References reco::LeafCandidate::p4(), and reco::LeafCandidate::setP4().
Referenced by JPTJetTester::analyze(), CMSDAS11DijetAnalyzer::analyze(), PFJetTester::analyze(), CaloJetTester::analyze(), TopSingleLepton::MonitorEnsemble::fill(), TopHLTSingleLepton::MonitorEnsemble::fill(), TopDiLeptonOffline::MonitorEnsemble::fill(), hitfit::JetTranslatorBase< AJet >::operator()(), CorrectJet::operator()(), JetPlusTrackProducerAA::produce(), JetPlusTrackProducer::produce(), and HiL1Subtractor::produce().
| virtual void reco::Jet::setJetArea | ( | float | fArea | ) | [inline, virtual] |
set jet area
Definition at line 104 of file Jet.h.
References mJetArea.
Referenced by FastjetJetProducer::produceTrackJets(), cms::SubEventGenJetProducer::runAlgorithm(), cms::CompoundJetProducer::writeCompoundJets(), and VirtualJetProducer::writeCompoundJets().
{mJetArea = fArea;}
| virtual void reco::Jet::setNPasses | ( | int | fPasses | ) | [inline, virtual] |
Set number of passes taken by algorithm.
Definition at line 114 of file Jet.h.
References mPassNumber.
{mPassNumber = fPasses;}
| virtual void reco::Jet::setPileup | ( | float | fEnergy | ) | [inline, virtual] |
Set pileup energy contribution as calculated by algorithm.
Definition at line 109 of file Jet.h.
References mPileupEnergy.
Referenced by HiL1Subtractor::produce(), FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().
{mPileupEnergy = fEnergy;}
float reco::Jet::mJetArea [private] |
Definition at line 121 of file Jet.h.
Referenced by jetArea(), and setJetArea().
int reco::Jet::mPassNumber [private] |
Definition at line 123 of file Jet.h.
Referenced by nPasses(), and setNPasses().
float reco::Jet::mPileupEnergy [private] |
Definition at line 122 of file Jet.h.
Referenced by pileup(), and setPileup().
1.7.3