CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Types | Public Member Functions | Static Public Member Functions | Private Attributes
reco::Jet Class Reference

Base class for all types of Jets. More...

#include <Jet.h>

Inheritance diagram for reco::Jet:
reco::CompositePtrCandidate reco::LeafCandidate reco::Candidate pat::PATObject< reco::Jet > reco::BasicJet reco::CaloJet reco::GenJet reco::JPTJet reco::PFClusterJet reco::PFJet reco::TrackJet

Classes

class  EtaPhiMoments
 record to store eta-phi first and second moments More...
 

Public Types

typedef edm::Ptr< CandidateConstituent
 
typedef std::vector< ConstituentConstituents
 
- Public Types inherited from reco::CompositePtrCandidate
typedef std::vector< CandidatePtrdaughters
 collection of references to daughters More...
 
typedef std::vector< CandidatePtrmothers
 collection of references to daughters More...
 
- Public Types inherited from reco::LeafCandidate
typedef int Charge
 electric charge type More...
 
typedef CandidateCollection daughters
 collection of daughter candidates More...
 
typedef unsigned int index
 
typedef math::XYZTLorentzVector LorentzVector
 Lorentz vector. More...
 
typedef math::XYZPoint Point
 point in the space More...
 
typedef
math::PtEtaPhiMLorentzVector 
PolarLorentzVector
 Lorentz vector. More...
 
typedef math::XYZVector Vector
 point in the space More...
 
- Public Types inherited from reco::Candidate
enum  { dimension = 3 }
 
enum  { size = dimension * (dimension + 1)/2 }
 matix size More...
 
typedef int Charge
 electric charge type More...
 
typedef candidate::const_iterator const_iterator
 
typedef math::Error< dimension >
::type 
CovarianceMatrix
 covariance error matrix (3x3) More...
 
typedef unsigned int index
 index type More...
 
typedef candidate::iterator iterator
 
typedef math::XYZTLorentzVector LorentzVector
 Lorentz vector. More...
 
typedef math::XYZPoint Point
 point in the space More...
 
typedef
math::PtEtaPhiMLorentzVector 
PolarLorentzVector
 Lorentz vector. More...
 
typedef size_t size_type
 
typedef math::XYZVector Vector
 point in the space More...
 

Public Member Functions

float constituentEtaPhiSpread () const
 
float constituentPtDistribution () const
 
float etaetaMoment () const
 eta-eta second moment, ET weighted More...
 
float etaphiMoment () const
 eta-phi second moment, ET weighted More...
 
EtaPhiMoments etaPhiStatistics () const
 eta-phi statistics, ET weighted More...
 
float etInAnnulus (float fRmin, float fRmax) const
 ET in annulus between rmin and rmax around jet direction. More...
 
virtual Constituents getJetConstituents () const
 list of constituents More...
 
virtual std::vector< const
reco::Candidate * > 
getJetConstituentsQuick () const
 quick list of constituents More...
 
bool isJet () const
 
 Jet ()
 Default constructor. More...
 
 Jet (const LorentzVector &fP4, const Point &fVertex)
 Initiator. More...
 
 Jet (const LorentzVector &fP4, const Point &fVertex, const Constituents &fConstituents)
 
virtual float jetArea () const
 get jet area More...
 
float maxDistance () const
 maximum distance from jet to constituent More...
 
int nCarrying (float fFraction) const
 return # of constituent carrying fraction of energy More...
 
virtual int nConstituents () const
 

of constituents

More...
 
virtual int nPasses () const
 number of passes taken by algorithm More...
 
float phiphiMoment () const
 phi-phi second moment, ET weighted More...
 
virtual float pileup () const
 pileup energy contribution as calculated by algorithm More...
 
virtual std::string print () const
 Print object. More...
 
virtual void scaleEnergy (double fScale)
 scale energy of the jet More...
 
virtual void setJetArea (float fArea)
 set jet area More...
 
virtual void setNPasses (int fPasses)
 Set number of passes taken by algorithm. More...
 
virtual void setPileup (float fEnergy)
 Set pileup energy contribution as calculated by algorithm. More...
 
virtual ~Jet ()
 Destructor. More...
 
- Public Member Functions inherited from reco::CompositePtrCandidate
void addDaughter (const CandidatePtr &)
 add a daughter via a reference More...
 
virtual const_iterator begin () const
 first daughter const_iterator More...
 
virtual iterator begin ()
 first daughter iterator More...
 
void clearDaughters ()
 clear daughter references More...
 
virtual CompositePtrCandidateclone () const
 returns a clone of the candidate More...
 
 CompositePtrCandidate ()
 default constructor More...
 
 CompositePtrCandidate (Charge q, const LorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 CompositePtrCandidate (Charge q, const PolarLorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 CompositePtrCandidate (const Candidate &p)
 constructor from a Candidate More...
 
virtual const Candidatedaughter (size_type) const
 return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) More...
 
virtual Candidatedaughter (size_type)
 return daughter at a given position, i = 0, ... numberOfDaughters() - 1 More...
 
CandidatePtr daughterPtr (size_type i) const
 reference to daughter at given position More...
 
const daughtersdaughterPtrVector () const
 references to daughtes More...
 
virtual const_iterator end () const
 last daughter const_iterator More...
 
virtual iterator end ()
 last daughter iterator More...
 
virtual const Candidatemother (size_t i=0) const
 return pointer to mother More...
 
virtual size_t numberOfDaughters () const
 number of daughters More...
 
virtual size_t numberOfMothers () const
 number of mothers More...
 
virtual size_type numberOfSourceCandidatePtrs () const
 
virtual CandidatePtr sourceCandidatePtr (size_type i) const
 
virtual ~CompositePtrCandidate ()
 destructor More...
 
- Public Member Functions inherited from reco::LeafCandidate
template<typename S >
daughter_iterator< S >::type beginFilter (const S &s) const
 
virtual Vector boostToCM () const
 
virtual int charge () const
 electric charge More...
 
virtual Candidatedaughter (const std::string &s)
 return daughter with a specified role name More...
 
virtual const Candidatedaughter (const std::string &s) const
 return daughter with a specified role name More...
 
template<typename S >
daughter_iterator< S >::type endFilter (const S &s) const
 
virtual double energy () const
 energy More...
 
virtual double et () const
 transverse energy More...
 
virtual float eta () const
 momentum pseudorapidity More...
 
virtual void fillVertexCovariance (CovarianceMatrix &v) const
 fill SMatrix More...
 
template<typename T >
T get () const
 get a component More...
 
template<typename T , typename Tag >
T get () const
 get a component More...
 
template<typename T >
T get (size_type i) const
 get a component More...
 
template<typename T , typename Tag >
T get (size_type i) const
 get a component More...
 
virtual bool hasMasterClone () const
 
virtual bool hasMasterClonePtr () const
 
virtual bool isCaloMuon () const
 
virtual bool isConvertedPhoton () const
 
virtual bool isElectron () const
 
virtual bool isGlobalMuon () const
 
virtual bool isMuon () const
 
virtual bool isPhoton () const
 
virtual bool isStandAloneMuon () const
 
virtual bool isTrackerMuon () const
 
 LeafCandidate ()
 default constructor More...
 
 LeafCandidate (const Candidate &c)
 
template<typename P4 >
 LeafCandidate (Charge q, const P4 &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from Any values More...
 
 LeafCandidate (Charge q, const PtEtaPhiMass &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 LeafCandidate (Charge q, const LorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 LeafCandidate (Charge q, const PolarLorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 LeafCandidate (Charge q, const GlobalVector &p3, float iEnergy, bool massless, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
 LeafCandidate (Charge q, const GlobalVector &p3, float iEnergy, float imass, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
 constructor from values More...
 
virtual bool longLived () const
 is long lived? More...
 
virtual float mass () const
 mass More...
 
virtual bool massConstraint () const
 do mass constraint? More...
 
virtual float massSqr () const
 mass squared More...
 
virtual const CandidateBaseRefmasterClone () const
 
virtual const CandidatePtrmasterClonePtr () const
 
template<typename Ref >
Ref masterRef () const
 cast master clone reference to a concrete type More...
 
virtual Vector momentum () const
 spatial momentum vector More...
 
virtual double mt () const
 transverse mass More...
 
virtual double mtSqr () const
 transverse mass squared More...
 
template<typename T >
size_type numberOf () const
 number of components More...
 
template<typename T , typename Tag >
size_type numberOf () const
 number of components More...
 
virtual double p () const
 magnitude of momentum vector More...
 
virtual const LorentzVectorp4 () const
 four-momentum Lorentz vector More...
 
virtual int pdgId () const
 PDG identifier. More...
 
virtual float phi () const
 momentum azimuthal angle More...
 
virtual const PolarLorentzVectorpolarP4 () const
 four-momentum Lorentz vector More...
 
virtual float pt () const
 transverse momentum More...
 
virtual double px () const
 x coordinate of momentum vector More...
 
virtual double py () const
 y coordinate of momentum vector More...
 
virtual double pz () const
 z coordinate of momentum vector More...
 
virtual double rapidity () const
 rapidity More...
 
virtual void setCharge (Charge q)
 set electric charge More...
 
virtual void setLongLived ()
 set long lived flag More...
 
virtual void setMass (double m)
 set particle mass More...
 
virtual void setMassConstraint ()
 set mass constraint flag More...
 
virtual void setP4 (const LorentzVector &p4)
 set 4-momentum More...
 
virtual void setP4 (const PolarLorentzVector &p4)
 set 4-momentum More...
 
virtual void setPdgId (int pdgId)
 
virtual void setPz (double pz)
 
virtual void setStatus (int status)
 set status word More...
 
virtual void setThreeCharge (Charge qx3)
 set electric charge More...
 
virtual void setVertex (const Point &vertex)
 set vertex More...
 
virtual int status () const
 status word More...
 
virtual double theta () const
 momentum polar angle More...
 
virtual int threeCharge () const
 electric charge More...
 
virtual const Pointvertex () const
 vertex position (overwritten by PF...) More...
 
virtual double vertexChi2 () const
 chi-squares More...
 
virtual double vertexCovariance (int i, int j) const
 (i, j)-th element of error matrix, i, j = 0, ... 2 More...
 
CovarianceMatrix vertexCovariance () const
 return SMatrix More...
 
virtual double vertexNdof () const
 
virtual double vertexNormalizedChi2 () const
 chi-squared divided by n.d.o.f. More...
 
virtual double vx () const
 x coordinate of vertex position More...
 
virtual double vy () const
 y coordinate of vertex position More...
 
virtual double vz () const
 z coordinate of vertex position More...
 
virtual double y () const
 rapidity More...
 
virtual ~LeafCandidate ()
 destructor More...
 
- Public Member Functions inherited from reco::Candidate
template<typename S >
daughter_iterator< S >::type beginFilter (const S &s) const
 
 Candidate ()
 default constructor More...
 
template<typename S >
daughter_iterator< S >::type endFilter (const S &s) const
 
template<typename T >
T get () const
 get a component More...
 
template<typename T , typename Tag >
T get () const
 get a component More...
 
template<typename T >
T get (size_type i) const
 get a component More...
 
template<typename T , typename Tag >
T get (size_type i) const
 get a component More...
 
template<typename Ref >
Ref masterRef () const
 cast master clone reference to a concrete type More...
 
template<typename T >
size_type numberOf () const
 number of components More...
 
template<typename T , typename Tag >
size_type numberOf () const
 number of components More...
 
virtual void setSourceCandidatePtr (const CandidatePtr &ptr)
 Set the ptr to the source Candidate. More...
 
virtual ~Candidate ()
 destructor More...
 

Static Public Member Functions

static float detectorEta (float fZVertex, float fPhysicsEta)
 static function to convert physics eta to detector eta More...
 
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 More...
 
static Candidate::LorentzVector physicsP4 (const Candidate::Point &newVertex, const Candidate &inParticle, const Candidate::Point &oldVertex=Candidate::Point(0, 0, 0))
 
- Static Public Member Functions inherited from reco::LeafCandidate
static double dmass (GlobalVector v, double e)
 
static double magd (GlobalVector v)
 

Private Attributes

float mJetArea
 
int mPassNumber
 
float mPileupEnergy
 

Additional Inherited Members

- Static Public Attributes inherited from reco::LeafCandidate
static const unsigned int longLivedTag = 65536
 long lived flag More...
 
static const unsigned int massConstraintTag = 131072
 do mass constraint flag More...
 
- Protected Member Functions inherited from reco::LeafCandidate
void cacheCartesian () const
 set internal cache More...
 
void cachePolar () const
 set internal cache More...
 
void clearCache () const
 clear internal cache More...
 
- Protected Attributes inherited from reco::LeafCandidate
bool cacheCartesianFixed_
 
bool cachePolarFixed_
 has cache been set? More...
 
float eta_
 
float mass_
 
LorentzVector p4Cartesian_
 internal cache for p4 More...
 
PolarLorentzVector p4Polar_
 internal cache for p4 More...
 
int pdgId_
 PDG identifier. More...
 
float phi_
 
float pt_
 four-momentum Lorentz vector More...
 
Charge qx3_
 electric charge More...
 
int status_
 status word More...
 
Point vertex_
 vertex position More...
 

Detailed Description

Base class for all types of Jets.

GenericJet describes jets made from arbitrary constituents, No direct constituents reference is stored for now

Author
Fedor Ratnikov, UMd
Version
Mar 23, 2007 by F.R.

Jet describes properties common for all kinds of jets, essentially kinematics. Base class for all types of Jets.

Author
Fedor Ratnikov, UMd
Version
Original: April 22, 2005 by Fernando Varela Rodriguez.
May 23, 2006 by F.R.

Definition at line 20 of file Jet.h.

Member Typedef Documentation

Definition at line 22 of file Jet.h.

typedef std::vector<Constituent> reco::Jet::Constituents

Definition at line 23 of file Jet.h.

Constructor & Destructor Documentation

reco::Jet::Jet ( )
inline

Default constructor.

Definition at line 36 of file Jet.h.

36 : mJetArea (0), mPileupEnergy (0), mPassNumber (0) {}
int mPassNumber
Definition: Jet.h:122
float mJetArea
Definition: Jet.h:120
float mPileupEnergy
Definition: Jet.h:121
Jet::Jet ( const LorentzVector fP4,
const Point fVertex 
)

Initiator.

Definition at line 174 of file Jet.cc.

176  : CompositePtrCandidate (0, fP4, fVertex),
177  mJetArea (0),
178  mPileupEnergy (0),
179  mPassNumber (0)
180 {}
int mPassNumber
Definition: Jet.h:122
CompositePtrCandidate()
default constructor
float mJetArea
Definition: Jet.h:120
float mPileupEnergy
Definition: Jet.h:121
Jet::Jet ( const LorentzVector fP4,
const Point fVertex,
const Constituents fConstituents 
)

Definition at line 161 of file Jet.cc.

References reco::CompositePtrCandidate::addDaughter(), and i.

164  : CompositePtrCandidate (0, fP4, fVertex),
165  mJetArea (0),
166  mPileupEnergy (0),
167  mPassNumber (0)
168 {
169  for (unsigned i = 0; i < fConstituents.size (); i++) {
170  addDaughter (fConstituents [i]);
171  }
172 }
int i
Definition: DBlmapReader.cc:9
int mPassNumber
Definition: Jet.h:122
CompositePtrCandidate()
default constructor
float mJetArea
Definition: Jet.h:120
void addDaughter(const CandidatePtr &)
add a daughter via a reference
float mPileupEnergy
Definition: Jet.h:121
virtual reco::Jet::~Jet ( )
inlinevirtual

Destructor.

Reimplemented in pat::Jet.

Definition at line 41 of file Jet.h.

41 {}

Member Function Documentation

float Jet::constituentEtaPhiSpread ( ) const

Definition at line 394 of file Jet.cc.

References reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, getJetConstituents(), and reco::LeafCandidate::pt().

Referenced by FFTGenericScaleCalculator::mapFFTJet().

394  {
395 
396  Jet::Constituents constituents = this->getJetConstituents();
397 
398 
399  float sum_pt2 = 0.;
400  float sum_pt2deltaR2 = 0.;
401 
402  for( unsigned iConst=0; iConst<constituents.size(); ++iConst ) {
403 
404  LorentzVector thisConstituent = constituents[iConst]->p4();
405 
406  float pt = thisConstituent.Pt();
407  float pt2 = pt*pt;
408  double dR = deltaR (*this, *(constituents[iConst]));
409  float pt2deltaR2 = pt*pt*dR*dR;
410 
411  sum_pt2 += pt2;
412  sum_pt2deltaR2 += pt2deltaR2;
413 
414  } //for constituents
415 
416  float rmsCand_value = (sum_pt2>0.) ? sum_pt2deltaR2/sum_pt2 : 0.;
417 
418  return rmsCand_value;
419 
420 } //constituentEtaPhiSpread
virtual float pt() const
transverse momentum
std::vector< Constituent > Constituents
Definition: Jet.h:23
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
virtual Constituents getJetConstituents() const
list of constituents
Definition: Jet.cc:350
float Jet::constituentPtDistribution ( ) const

Definition at line 369 of file Jet.cc.

References getJetConstituents(), reco::LeafCandidate::pt(), and mathSSE::sqrt().

Referenced by FFTGenericScaleCalculator::mapFFTJet().

369  {
370 
371  Jet::Constituents constituents = this->getJetConstituents();
372 
373  float sum_pt2 = 0.;
374  float sum_pt = 0.;
375 
376  for( unsigned iConst=0; iConst<constituents.size(); ++iConst ) {
377 
378  float pt = constituents[iConst]->p4().Pt();
379  float pt2 = pt*pt;
380 
381  sum_pt += pt;
382  sum_pt2 += pt2;
383 
384  } //for constituents
385 
386  float ptD_value = (sum_pt>0.) ? sqrt( sum_pt2 / (sum_pt*sum_pt) ) : 0.;
387 
388  return ptD_value;
389 
390 } //constituentPtDistribution
virtual float pt() const
transverse momentum
std::vector< Constituent > Constituents
Definition: Jet.h:23
T sqrt(T t)
Definition: SSEVec.h:48
virtual Constituents getJetConstituents() const
list of constituents
Definition: Jet.cc:350
float Jet::detectorEta ( float  fZVertex,
float  fPhysicsEta 
)
static

static function to convert physics eta to detector eta

Definition at line 325 of file Jet.cc.

Referenced by reco::GenJet::detectorEta().

325  {
326  CaloPointZ refPoint (fZVertex, fPhysicsEta);
327  return refPoint.etaReference (0.);
328 }
Candidate::LorentzVector Jet::detectorP4 ( const Candidate::Point vertex,
const Candidate inParticle 
)
static

Definition at line 339 of file Jet.cc.

References reco::Candidate::energy(), reco::Candidate::momentum(), np, reco::LeafCandidate::p(), and p3.

Referenced by reco::CaloJet::detectorP4().

339  {
340  CaloPoint3D<Point> caloPoint(vertex,inParticle.momentum()); // Jet position in Calo.
341  static const Point np(0,0,0);
342  Vector detectorDir = caloPoint.caloPoint() - np;
343  double p = inParticle.momentum().r();
344  Vector p3 = p * detectorDir.unit();
345  LorentzVector returnVector(p3.x(), p3.y(), p3.z(), inParticle.energy());
346  return returnVector;
347 }
virtual double p() const
magnitude of momentum vector
virtual double energy() const =0
energy
virtual const Point & vertex() const
vertex position (overwritten by PF...)
virtual Vector momentum() const =0
spatial momentum vector
int np
Definition: AMPTWrapper.h:33
fixed size vector
Definition: Vector.h:31
math::XYZPoint Point
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
double p3[4]
Definition: TauolaWrapper.h:91
float Jet::etaetaMoment ( ) const

eta-eta second moment, ET weighted

eta-eta second moment

Definition at line 223 of file Jet.cc.

References getJetConstituentsQuick(), i, relativeConstraints::value, and histoStyle::weight.

Referenced by fireworks::makeEveJetCone(), and JetIDSelectionFunctor::operator()().

223  {
224  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
225  double sumw = 0;
226  double sum = 0;
227  double sum2 = 0;
228  int i = towers.size();
229  while (--i >= 0) {
230  double value = towers[i]->eta();
231  double weight = towers[i]->et();
232  sumw += weight;
233  sum += value * weight;
234  sum2 += value * value * weight;
235  }
236  return sumw > 0 ? (sum2 - sum*sum/sumw ) / sumw : 0;
237 }
int i
Definition: DBlmapReader.cc:9
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
int weight
Definition: histoStyle.py:50
float Jet::etaphiMoment ( ) const

eta-phi second moment, ET weighted

eta-phi second moment

Definition at line 258 of file Jet.cc.

References reco::deltaPhi(), getJetConstituentsQuick(), i, reco::LeafCandidate::phi(), and histoStyle::weight.

258  {
259  double phiRef = phi();
260  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
261  double sumw = 0;
262  double sumA = 0;
263  double sumB = 0;
264  double sumAB = 0;
265  int i = towers.size();
266  while (--i >= 0) {
267  double valueA = towers[i]->eta();
268  double valueB = deltaPhi (towers[i]->phi(), phiRef);
269  double weight = towers[i]->et();
270  sumw += weight;
271  sumA += valueA * weight;
272  sumB += valueB * weight;
273  sumAB += valueA * valueB * weight;
274  }
275  return sumw > 0 ? (sumAB - sumA*sumB/sumw ) / sumw : 0;
276 }
int i
Definition: DBlmapReader.cc:9
virtual float phi() const
momentum azimuthal angle
double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:12
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
int weight
Definition: histoStyle.py:50
Jet::EtaPhiMoments Jet::etaPhiStatistics ( ) const

eta-phi statistics, ET weighted

eta-phi statistics

Definition at line 183 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 histoStyle::weight.

183  {
184  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
185  double phiRef = phi();
186  double sumw = 0;
187  double sumEta = 0;
188  double sumEta2 = 0;
189  double sumPhi = 0;
190  double sumPhi2 = 0;
191  double sumEtaPhi = 0;
192  int i = towers.size();
193  while (--i >= 0) {
194  double eta = towers[i]->eta();
195  double phi = deltaPhi (towers[i]->phi(), phiRef);
196  double weight = towers[i]->et();
197  sumw += weight;
198  sumEta += eta * weight;
199  sumEta2 += eta * eta * weight;
200  sumPhi += phi * weight;
201  sumPhi2 += phi * phi * weight;
202  sumEtaPhi += eta * phi * weight;
203  }
205  if (sumw > 0) {
206  result.etaMean = sumEta / sumw;
207  result.phiMean = deltaPhi (phiRef + sumPhi, 0.);
208  result.etaEtaMoment = (sumEta2 - sumEta * sumEta / sumw) / sumw;
209  result.phiPhiMoment = (sumPhi2 - sumPhi * sumPhi / sumw) / sumw;
210  result.etaPhiMoment = (sumEtaPhi - sumEta * sumPhi / sumw) / sumw;
211  }
212  else {
213  result.etaMean = 0;
214  result.phiMean = 0;
215  result.etaEtaMoment = 0;
216  result.phiPhiMoment = 0;
217  result.etaPhiMoment = 0;
218  }
219  return result;
220 }
int i
Definition: DBlmapReader.cc:9
virtual float phi() const
momentum azimuthal angle
record to store eta-phi first and second moments
Definition: Jet.h:26
virtual float eta() const
momentum pseudorapidity
tuple result
Definition: query.py:137
double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:12
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
int weight
Definition: histoStyle.py:50
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 279 of file Jet.cc.

References reco::deltaR(), getJetConstituentsQuick(), i, and query::result.

279  {
280  float result = 0;
281  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
282  int i = towers.size ();
283  while (--i >= 0) {
284  double r = deltaR (*this, *(towers[i]));
285  if (r >= fRmin && r < fRmax) result += towers[i]->et ();
286  }
287  return result;
288 }
int i
Definition: DBlmapReader.cc:9
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
tuple result
Definition: query.py:137
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
Jet::Constituents Jet::getJetConstituents ( ) const
virtual
std::vector< const Candidate * > Jet::getJetConstituentsQuick ( ) const
virtual

quick list of constituents

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

358  {
359  std::vector<const Candidate*> result;
360  int nDaughters = numberOfDaughters();
361  for (int i = 0; i < nDaughters; ++i) {
362  result.push_back (daughter (i));
363  }
364  return result;
365 }
int i
Definition: DBlmapReader.cc:9
virtual size_t numberOfDaughters() const
number of daughters
tuple result
Definition: query.py:137
virtual const Candidate * daughter(size_type) const
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) ...
bool Jet::isJet ( ) const
virtual

Reimplemented from reco::LeafCandidate.

Definition at line 448 of file Jet.cc.

448  {
449  return true;
450 }
virtual float reco::Jet::jetArea ( ) const
inlinevirtual
float Jet::maxDistance ( ) const

maximum distance from jet to constituent

Definition at line 306 of file Jet.cc.

References reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, getJetConstituentsQuick(), i, and query::result.

306  {
307  float result = 0;
308  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
309  for (unsigned i = 0; i < towers.size(); ++i) {
310  float dR = deltaR (*(towers[i]), *this);
311  if (dR > result) result = dR;
312  }
313  return result;
314 }
int i
Definition: DBlmapReader.cc:9
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
tuple result
Definition: query.py:137
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
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 291 of file Jet.cc.

References reco::LeafCandidate::et(), getJetConstituentsQuick(), i, and query::result.

Referenced by reco::CaloJet::n60(), pat::Jet::n60(), reco::CaloJet::n90(), and pat::Jet::n90().

291  {
292  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
293  if (fFraction >= 1) return towers.size();
294  double totalEt = 0;
295  for (unsigned i = 0; i < towers.size(); ++i) totalEt += towers[i]->et();
296  double fractionEnergy = totalEt * fFraction;
297  unsigned result = 0;
298  for (; result < towers.size(); ++result) {
299  fractionEnergy -= towers[result]->et();
300  if (fractionEnergy <= 0) return result+1;
301  }
302  return 0;
303 }
int i
Definition: DBlmapReader.cc:9
virtual double et() const
transverse energy
tuple result
Definition: query.py:137
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
virtual int reco::Jet::nConstituents ( ) const
inlinevirtual
virtual int reco::Jet::nPasses ( ) const
inlinevirtual

number of passes taken by algorithm

Definition at line 115 of file Jet.h.

References mPassNumber.

115 {return mPassNumber;}
int mPassNumber
Definition: Jet.h:122
float Jet::phiphiMoment ( ) const

phi-phi second moment, ET weighted

phi-phi second moment

Definition at line 240 of file Jet.cc.

References reco::deltaPhi(), getJetConstituentsQuick(), i, reco::LeafCandidate::phi(), relativeConstraints::value, and histoStyle::weight.

Referenced by fireworks::makeEveJetCone(), and JetIDSelectionFunctor::operator()().

240  {
241  double phiRef = phi();
242  std::vector<const Candidate*> towers = getJetConstituentsQuick ();
243  double sumw = 0;
244  double sum = 0;
245  double sum2 = 0;
246  int i = towers.size();
247  while (--i >= 0) {
248  double value = deltaPhi (towers[i]->phi(), phiRef);
249  double weight = towers[i]->et();
250  sumw += weight;
251  sum += value * weight;
252  sum2 += value * value * weight;
253  }
254  return sumw > 0 ? (sum2 - sum*sum/sumw ) / sumw : 0;
255 }
int i
Definition: DBlmapReader.cc:9
virtual float phi() const
momentum azimuthal angle
double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:12
virtual std::vector< const reco::Candidate * > getJetConstituentsQuick() const
quick list of constituents
Definition: Jet.cc:358
int weight
Definition: histoStyle.py:50
float Jet::physicsEta ( float  fZVertex,
float  fDetectorEta 
)
static

static function to convert detector eta to physics eta

Definition at line 318 of file Jet.cc.

Referenced by JetTracksAssociationDRVertexAssigned::produce().

318  {
319  CaloPointZ refPoint (0., fDetectorEta);
320  return refPoint.etaReference (fZVertex);
321 }
Candidate::LorentzVector Jet::physicsP4 ( const Candidate::Point newVertex,
const Candidate inParticle,
const Candidate::Point oldVertex = Candidate::Point(0,0,0) 
)
static

Definition at line 330 of file Jet.cc.

References reco::Candidate::energy(), reco::Candidate::momentum(), reco::LeafCandidate::p(), and p3.

Referenced by reco::CaloJet::physicsP4().

330  {
331  CaloPoint3D<Point> caloPoint(oldVertex,inParticle.momentum()); // Jet position in Calo.
332  Vector physicsDir = caloPoint.caloPoint() - newVertex;
333  double p = inParticle.momentum().r();
334  Vector p3 = p * physicsDir.unit();
335  LorentzVector returnVector(p3.x(), p3.y(), p3.z(), inParticle.energy());
336  return returnVector;
337 }
virtual double p() const
magnitude of momentum vector
virtual double energy() const =0
energy
virtual Vector momentum() const =0
spatial momentum vector
fixed size vector
Definition: Vector.h:31
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
double p3[4]
Definition: TauolaWrapper.h:91
virtual float reco::Jet::pileup ( ) const
inlinevirtual

pileup energy contribution as calculated by algorithm

Definition at line 110 of file Jet.h.

References mPileupEnergy.

110 {return mPileupEnergy;}
float mPileupEnergy
Definition: Jet.h:121
std::string Jet::print ( void  ) const
virtual

Print object.

Reimplemented in reco::PFJet, reco::CaloJet, reco::JPTJet, reco::GenJet, reco::TrackJet, reco::PFClusterJet, and reco::BasicJet.

Definition at line 424 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::BasicJet::print(), reco::PFClusterJet::print(), reco::TrackJet::print(), reco::GenJet::print(), reco::JPTJet::print(), reco::CaloJet::print(), and reco::PFJet::print().

424  {
425  std::ostringstream out;
426  out << "Jet p/px/py/pz/pt: " << p() << '/' << px () << '/' << py() << '/' << pz() << '/' << pt() << std::endl
427  << " eta/phi: " << eta () << '/' << phi () << std::endl
428  << " # of constituents: " << nConstituents () << std::endl;
429  out << " Constituents:" << std::endl;
430  for (unsigned index = 0; index < numberOfDaughters(); index++) {
431  Constituent constituent = daughterPtr (index); // deref
432  if (constituent.isNonnull()) {
433  out << " #" << index << " p/pt/eta/phi: "
434  << constituent->p() << '/' << constituent->pt() << '/' << constituent->eta() << '/' << constituent->phi()
435  << " productId/index: " << constituent.id() << '/' << constituent.key() << std::endl;
436  }
437  else {
438  out << " #" << index << " constituent is not available in the event" << std::endl;
439  }
440  }
441  return out.str ();
442 }
virtual double p() const
magnitude of momentum vector
virtual float pt() const
transverse momentum
CandidatePtr daughterPtr(size_type i) const
reference to daughter at given position
virtual float phi() const
momentum azimuthal angle
virtual size_t numberOfDaughters() const
number of daughters
virtual float eta() const
momentum pseudorapidity
unsigned int index
Definition: LeafCandidate.h:36
tuple out
Definition: dbtoconf.py:99
edm::Ptr< Candidate > Constituent
Definition: Jet.h:22
virtual double px() const
x coordinate of momentum vector
virtual double pz() const
z coordinate of momentum vector
virtual int nConstituents() const
of constituents
Definition: Jet.h:65
virtual double py() const
y coordinate of momentum vector
void Jet::scaleEnergy ( double  fScale)
virtual
virtual void reco::Jet::setJetArea ( float  fArea)
inlinevirtual
virtual void reco::Jet::setNPasses ( int  fPasses)
inlinevirtual

Set number of passes taken by algorithm.

Definition at line 113 of file Jet.h.

References mPassNumber.

Referenced by pat::PATJetSlimmer::produce().

113 {mPassNumber = fPasses;}
int mPassNumber
Definition: Jet.h:122
virtual void reco::Jet::setPileup ( float  fEnergy)
inlinevirtual

Set pileup energy contribution as calculated by algorithm.

Definition at line 108 of file Jet.h.

References mPileupEnergy.

Referenced by HiL1Subtractor::produce(), FastjetJetProducer::produceTrackJets(), and cms::SubEventGenJetProducer::runAlgorithm().

108 {mPileupEnergy = fEnergy;}
float mPileupEnergy
Definition: Jet.h:121

Member Data Documentation

float reco::Jet::mJetArea
private

Definition at line 120 of file Jet.h.

Referenced by jetArea(), and setJetArea().

int reco::Jet::mPassNumber
private

Definition at line 122 of file Jet.h.

Referenced by nPasses(), and setNPasses().

float reco::Jet::mPileupEnergy
private

Definition at line 121 of file Jet.h.

Referenced by pileup(), and setPileup().