CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
l1t::HGCalClusterT< C > Class Template Reference

#include <HGCalClusterT.h>

Inheritance diagram for l1t::HGCalClusterT< C >:
l1t::L1Candidate reco::LeafCandidate reco::Candidate

Public Types

typedef std::unordered_map
< uint32_t, edm::Ptr< C >
>::const_iterator 
const_iterator
 
- 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

void addConstituent (const edm::Ptr< C > &c, bool updateCentre=true, float fraction=1.)
 
const GlobalPointcentre () const
 
const GlobalPointcentreProj () const
 
const std::unordered_map
< uint32_t, edm::Ptr< C > > & 
constituents () const
 
const_iterator constituents_begin () const
 
const_iterator constituents_end () const
 
int coreShowerLength () const
 
void coreShowerLength (int coreShowerLength)
 
uint32_t detId () const
 
double distance (const l1t::HGCalTriggerCell &tc) const
 
float eMax () const
 
void eMax (float eMax)
 
int firstLayer () const
 
void firstLayer (int firstLayer)
 
 HGCalClusterT ()
 
 HGCalClusterT (const LorentzVector p4, int pt=0, int eta=0, int phi=0)
 
 HGCalClusterT (const edm::Ptr< C > &c, float fraction=1.)
 
double hOverE () const
 
float layer10percent () const
 
void layer10percent (float layer10percent)
 
float layer50percent () const
 
void layer50percent (float layer50percent)
 
float layer90percent () const
 
void layer90percent (float layer90percent)
 
int maxLayer () const
 
void maxLayer (int maxLayer)
 
double mipPt () const
 
bool operator< (const HGCalClusterT< C > &cl) const
 
bool operator<= (const HGCalClusterT< C > &cl) const
 
bool operator> (const HGCalClusterT< C > &cl) const
 
bool operator>= (const HGCalClusterT< C > &cl) const
 
const GlobalPointposition () const
 
void removeConstituent (const edm::Ptr< C > &c, bool updateCentre=true)
 
double seedMipPt () const
 
void setDetId (uint32_t id)
 
void setPt (double pt)
 
void setValid (bool valid)
 
int showerLength () const
 
void showerLength (int showerLength)
 
float sigmaEtaEtaMax () const
 
void sigmaEtaEtaMax (float sigmaEtaEtaMax)
 
float sigmaEtaEtaTot () const
 
void sigmaEtaEtaTot (float sigmaEtaEtaTot)
 
float sigmaPhiPhiMax () const
 
void sigmaPhiPhiMax (float sigmaPhiPhiMax)
 
float sigmaPhiPhiTot () const
 
void sigmaPhiPhiTot (float sigmaPhiPhiTot)
 
float sigmaRRMax () const
 
void sigmaRRMax (float sigmaRRMax)
 
float sigmaRRMean () const
 
void sigmaRRMean (float sigmaRRMean)
 
float sigmaRRTot () const
 
void sigmaRRTot (float sigmaRRTot)
 
float sigmaZZ () const
 
void sigmaZZ (float sigmaZZ)
 
unsigned size () const
 
uint32_t subdetId () const
 
double sumPt () const
 
float triggerCells67percent () const
 
void triggerCells67percent (float triggerCells67percent)
 
float triggerCells90percent () const
 
void triggerCells90percent (float triggerCells90percent)
 
bool valid () const
 
float zBarycenter () const
 
void zBarycenter (float zBarycenter)
 
 ~HGCalClusterT () override
 
- Public Member Functions inherited from l1t::L1Candidate
int hwEta () const
 
int hwIso () const
 
int hwPhi () const
 
int hwPt () const
 
int hwQual () const
 
 L1Candidate ()
 
 L1Candidate (const LorentzVector &p4, int pt=0, int eta=0, int phi=0, int qual=0, int iso=0)
 
 L1Candidate (const PolarLorentzVector &p4, int pt=0, int eta=0, int phi=0, int qual=0, int iso=0)
 
virtual bool operator!= (const l1t::L1Candidate &rhs) const
 
virtual bool operator== (const l1t::L1Candidate &rhs) const
 
void setHwEta (int eta)
 
void setHwIso (int iso)
 
void setHwPhi (int phi)
 
void setHwPt (int pt)
 
void setHwQual (int qual)
 
 ~L1Candidate () override
 
- Public Member Functions inherited from reco::LeafCandidate
Vector boostToCM () const final
 
int charge () const final
 electric charge More...
 
LeafCandidateclone () const override
 returns a clone of the Candidate object More...
 
void construct (int qx3, float pt, float eta, float phi, float mass, const Point &vtx, int pdgId, int status)
 
const Candidatedaughter (size_type) const override
 return daughter at a given position (throws an exception) More...
 
Candidatedaughter (size_type) override
 return daughter at a given position (throws an exception) More...
 
Candidatedaughter (const std::string &s) override
 return daughter with a specified role name More...
 
const Candidatedaughter (const std::string &s) const override
 return daughter with a specified role name More...
 
double energy () const final
 energy More...
 
double et () const final
 transverse energy More...
 
double et2 () const final
 transverse energy squared (use this for cut!) More...
 
double eta () const final
 momentum pseudorapidity More...
 
void fillVertexCovariance (CovarianceMatrix &v) const override
 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...
 
bool hasMasterClone () const override
 
bool hasMasterClonePtr () const override
 
bool isCaloMuon () const override
 
bool isConvertedPhoton () const override
 
bool isElectron () const override
 
bool isGlobalMuon () const override
 
bool isJet () const override
 
bool isMuon () const override
 
bool isPhoton () const override
 
bool isStandAloneMuon () const override
 
bool isTrackerMuon () const override
 
 LeafCandidate ()
 
 LeafCandidate (const Candidate &c)
 
template<typename... Args>
 LeafCandidate (Args &&...args)
 
 LeafCandidate (LeafCandidate &rh)
 
 LeafCandidate (LeafCandidate &&)=default
 
 LeafCandidate (LeafCandidate const &)=default
 
bool longLived () const final
 is long lived? More...
 
double mass () const final
 mass More...
 
bool massConstraint () const final
 do mass constraint? More...
 
double massSqr () const final
 mass squared More...
 
const CandidateBaseRefmasterClone () const override
 
const CandidatePtrmasterClonePtr () const override
 
template<typename Ref >
Ref masterRef () const
 cast master clone reference to a concrete type More...
 
Vector momentum () const final
 spatial momentum vector More...
 
const Candidatemother (size_type) const override
 return mother at a given position (throws an exception) More...
 
double mt () const final
 transverse mass More...
 
double mtSqr () const final
 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...
 
size_t numberOfDaughters () const override
 number of daughters More...
 
size_t numberOfMothers () const override
 number of mothers More...
 
size_t numberOfSourceCandidatePtrs () const override
 
LeafCandidateoperator= (LeafCandidate &&)=default
 
LeafCandidateoperator= (LeafCandidate const &)=default
 
double p () const final
 magnitude of momentum vector More...
 
const LorentzVectorp4 () const final
 four-momentum Lorentz vector More...
 
int pdgId () const final
 PDG identifier. More...
 
double phi () const final
 momentum azimuthal angle More...
 
const PolarLorentzVectorpolarP4 () const final
 four-momentum Lorentz vector More...
 
double pt () const final
 transverse momentum More...
 
double px () const final
 x coordinate of momentum vector More...
 
double py () const final
 y coordinate of momentum vector More...
 
double pz () const final
 z coordinate of momentum vector More...
 
double rapidity () const final
 rapidity More...
 
void setCharge (Charge q) final
 set electric charge More...
 
void setLongLived () final
 
void setMass (double m) final
 set particle mass More...
 
void setMassConstraint () final
 
void setP4 (const LorentzVector &p4) final
 set 4-momentum More...
 
void setP4 (const PolarLorentzVector &p4) final
 set 4-momentum More...
 
void setPdgId (int pdgId) final
 
void setPz (double pz) final
 
void setStatus (int status) final
 set status word More...
 
void setThreeCharge (Charge qx3) final
 set electric charge More...
 
void setVertex (const Point &vertex) override
 set vertex More...
 
CandidatePtr sourceCandidatePtr (size_type i) const override
 
int status () const final
 status word More...
 
double theta () const final
 momentum polar angle More...
 
int threeCharge () const final
 electric charge More...
 
const Pointvertex () const override
 vertex position (overwritten by PF...) More...
 
double vertexChi2 () const override
 chi-squares More...
 
double vertexCovariance (int i, int j) const override
 (i, j)-th element of error matrix, i, j = 0, ... 2 More...
 
CovarianceMatrix vertexCovariance () const final
 return SMatrix More...
 
double vertexNdof () const override
 
double vertexNormalizedChi2 () const override
 chi-squared divided by n.d.o.f. More...
 
double vx () const override
 x coordinate of vertex position More...
 
double vy () const override
 y coordinate of vertex position More...
 
double vz () const override
 z coordinate of vertex position More...
 
double y () const final
 rapidity More...
 
 ~LeafCandidate () override
 destructor More...
 
- Public Member Functions inherited from reco::Candidate
const_iterator begin () const
 first daughter const_iterator More...
 
iterator begin ()
 first daughter iterator More...
 
virtual const TrackbestTrack () const
 
 Candidate ()
 default constructor More...
 
virtual float dxyError () const
 uncertainty on dxy More...
 
virtual float dzError () const
 uncertainty on dz More...
 
const_iterator end () const
 last daughter const_iterator More...
 
iterator end ()
 last daughter iterator 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...
 
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...
 

Private Member Functions

void updateP4AndPosition (const edm::Ptr< C > &c, bool updateCentre=true, float fraction=1.)
 

Private Attributes

GlobalPoint centre_
 
GlobalPoint centreProj_
 
std::unordered_map< uint32_t,
edm::Ptr< C > > 
constituents_
 
std::unordered_map< uint32_t,
double > 
constituentsFraction_
 
int coreShowerLength_ = 0
 
DetId detId_
 
float eMax_ = 0.
 
int firstLayer_ = 0
 
float layer10percent_ = 0.
 
float layer50percent_ = 0.
 
float layer90percent_ = 0.
 
int maxLayer_ = 0
 
double mipPt_ = 0.
 
double seedMipPt_ = 0.
 
int showerLength_ = 0
 
float sigmaEtaEtaMax_ = 0.
 
float sigmaEtaEtaTot_ = 0.
 
float sigmaPhiPhiMax_ = 0.
 
float sigmaPhiPhiTot_ = 0.
 
float sigmaRRMax_ = 0.
 
float sigmaRRMean_ = 0.
 
float sigmaRRTot_ = 0.
 
float sigmaZZ_ = 0.
 
double sumPt_ = 0.
 
float triggerCells67percent_ = 0.
 
float triggerCells90percent_ = 0.
 
bool valid_ = false
 
float zBarycenter_ = 0.
 

Additional Inherited Members

Detailed Description

template<class C>
class l1t::HGCalClusterT< C >

Definition at line 20 of file HGCalClusterT.h.

Member Typedef Documentation

template<class C>
typedef std::unordered_map<uint32_t, edm::Ptr<C> >::const_iterator l1t::HGCalClusterT< C >::const_iterator

Definition at line 22 of file HGCalClusterT.h.

Constructor & Destructor Documentation

template<class C>
l1t::HGCalClusterT< C >::HGCalClusterT ( )
inline

Definition at line 25 of file HGCalClusterT.h.

25 {}
template<class C>
l1t::HGCalClusterT< C >::HGCalClusterT ( const LorentzVector  p4,
int  pt = 0,
int  eta = 0,
int  phi = 0 
)
inline

Definition at line 26 of file HGCalClusterT.h.

27  : L1Candidate(p4, pt, eta, phi),
28  valid_(true),
29  detId_(0),
30  centre_(0, 0, 0),
31  centreProj_(0., 0., 0.),
32  mipPt_(0),
33  seedMipPt_(0),
34  sumPt_(0) {}
double pt() const final
transverse momentum
const LorentzVector & p4() const final
four-momentum Lorentz vector
GlobalPoint centreProj_
GlobalPoint centre_
double phi() const final
momentum azimuthal angle
double eta() const final
momentum pseudorapidity
template<class C>
l1t::HGCalClusterT< C >::HGCalClusterT ( const edm::Ptr< C > &  c,
float  fraction = 1. 
)
inline

Definition at line 36 of file HGCalClusterT.h.

37  : valid_(true),
38  detId_(c->detId()),
39  centre_(0., 0., 0.),
40  centreProj_(0., 0., 0.),
41  mipPt_(0.),
42  seedMipPt_(0.),
43  sumPt_(0.) {
44  addConstituent(c, true, fraction);
45  }
void addConstituent(const edm::Ptr< C > &c, bool updateCentre=true, float fraction=1.)
Definition: HGCalClusterT.h:54
GlobalPoint centreProj_
GlobalPoint centre_
template<class C>
l1t::HGCalClusterT< C >::~HGCalClusterT ( )
inlineoverride

Definition at line 47 of file HGCalClusterT.h.

47 {};

Member Function Documentation

template<class C>
void l1t::HGCalClusterT< C >::addConstituent ( const edm::Ptr< C > &  c,
bool  updateCentre = true,
float  fraction = 1. 
)
inline

Definition at line 54 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::HGCalClusterT(), HGCalClusteringImpl::mergeClusters(), and CaloTruthCellsProducer::produce().

54  {
55  double cMipt = c->mipPt() * fraction;
56 
57  if (constituents_.empty()) {
58  detId_ = DetId(c->detId());
59  seedMipPt_ = cMipt;
60  /* if the centre will not be dynamically calculated
61  the seed centre is considere as cluster centre */
62  if (!updateCentre) {
63  centre_ = c->position();
64  }
65  }
66  updateP4AndPosition(c, updateCentre, fraction);
67 
68  constituents_.emplace(c->detId(), c);
69  constituentsFraction_.emplace(c->detId(), fraction);
70  }
const edm::EventSetup & c
std::unordered_map< uint32_t, edm::Ptr< C > > constituents_
void updateP4AndPosition(const edm::Ptr< C > &c, bool updateCentre=true, float fraction=1.)
std::unordered_map< uint32_t, double > constituentsFraction_
GlobalPoint centre_
Definition: DetId.h:17
template<class C>
const GlobalPoint& l1t::HGCalClusterT< C >::centre ( ) const
inline

Definition at line 100 of file HGCalClusterT.h.

Referenced by CaloTruthCellsProducer::produce().

100 { return centre_; }
GlobalPoint centre_
template<class C>
const GlobalPoint& l1t::HGCalClusterT< C >::centreProj ( ) const
inline
template<class C>
const std::unordered_map<uint32_t, edm::Ptr<C> >& l1t::HGCalClusterT< C >::constituents ( ) const
inline
template<class C>
const_iterator l1t::HGCalClusterT< C >::constituents_begin ( ) const
inline

Definition at line 50 of file HGCalClusterT.h.

50 { return constituents_.begin(); }
std::unordered_map< uint32_t, edm::Ptr< C > > constituents_
template<class C>
const_iterator l1t::HGCalClusterT< C >::constituents_end ( ) const
inline

Definition at line 51 of file HGCalClusterT.h.

51 { return constituents_.end(); }
std::unordered_map< uint32_t, edm::Ptr< C > > constituents_
template<class C>
int l1t::HGCalClusterT< C >::coreShowerLength ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::coreShowerLength ( int  coreShowerLength)
inline
template<class C>
uint32_t l1t::HGCalClusterT< C >::detId ( ) const
inline
template<class C>
double l1t::HGCalClusterT< C >::distance ( const l1t::HGCalTriggerCell tc) const
inline

Definition at line 97 of file HGCalClusterT.h.

Referenced by HGCalClusteringImpl::isPertinent().

97 { return (tc.position() - centre_).mag(); }
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
const GlobalPoint & position() const
GlobalPoint centre_
template<class C>
float l1t::HGCalClusterT< C >::eMax ( ) const
inline

Definition at line 146 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

146 { return eMax_; }
template<class C>
void l1t::HGCalClusterT< C >::eMax ( float  eMax)
inline

Definition at line 166 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::eMax().

166 { eMax_ = eMax; }
float eMax() const
template<class C>
int l1t::HGCalClusterT< C >::firstLayer ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::firstLayer ( int  firstLayer)
inline
template<class C>
double l1t::HGCalClusterT< C >::hOverE ( ) const
inline

Definition at line 103 of file HGCalClusterT.h.

Referenced by l1t::HGCalMulticluster::hOverE(), and l1t::HGCalMulticluster::saveHOverE().

103  {
104  double pt_em = 0.;
105  double pt_had = 0.;
106  double hOe = 0.;
107 
108  for (const auto& id_constituent : constituents()) {
109  DetId id(id_constituent.first);
110  auto id_fraction = constituentsFraction_.find(id_constituent.first);
111  double fraction = (id_fraction != constituentsFraction_.end() ? id_fraction->second : 1.);
112  if ((id.det() == DetId::HGCalEE) ||
113  (id.det() == DetId::HGCalTrigger &&
115  (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose && HFNoseDetId(id).isEE()) ||
116  (id.det() == DetId::HGCalTrigger &&
118  HFNoseTriggerDetId(id).isEE())) {
119  pt_em += id_constituent.second->pt() * fraction;
120  } else if ((id.det() == DetId::HGCalHSi) || (id.det() == DetId::HGCalHSc) ||
121  (id.det() == DetId::HGCalTrigger &&
123  (id.det() == DetId::Forward && id.subdetId() == ForwardSubdetector::HFNose &&
124  HFNoseDetId(id).isHE()) ||
125  (id.det() == DetId::HGCalTrigger &&
127  HFNoseTriggerDetId(id).isHSilicon())) {
128  pt_had += id_constituent.second->pt() * fraction;
129  }
130  }
131  if (pt_em > 0)
132  hOe = pt_had / pt_em;
133  else
134  hOe = -1.;
135  return hOe;
136  }
HGCalTriggerSubdetector subdet() const
get the subdetector
uint16_t *__restrict__ id
bool isEE() const
consistency check : no bits left =&gt; no overhead
Definition: HFNoseDetId.h:104
uint32_t subdetId() const
bool isEE() const
consistency check : no bits left =&gt; no overhead
bool isHE(int etabin, int depth)
std::unordered_map< uint32_t, double > constituentsFraction_
Definition: DetId.h:17
const std::unordered_map< uint32_t, edm::Ptr< C > > & constituents() const
Definition: HGCalClusterT.h:49
template<class C>
float l1t::HGCalClusterT< C >::layer10percent ( ) const
inline

Definition at line 156 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

156 { return layer10percent_; }
template<class C>
void l1t::HGCalClusterT< C >::layer10percent ( float  layer10percent)
inline
template<class C>
float l1t::HGCalClusterT< C >::layer50percent ( ) const
inline

Definition at line 157 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

157 { return layer50percent_; }
template<class C>
void l1t::HGCalClusterT< C >::layer50percent ( float  layer50percent)
inline
template<class C>
float l1t::HGCalClusterT< C >::layer90percent ( ) const
inline

Definition at line 158 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

158 { return layer90percent_; }
template<class C>
void l1t::HGCalClusterT< C >::layer90percent ( float  layer90percent)
inline
template<class C>
int l1t::HGCalClusterT< C >::maxLayer ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::maxLayer ( int  maxLayer)
inline

Definition at line 165 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::maxLayer().

165 { maxLayer_ = maxLayer; }
int maxLayer() const
template<class C>
double l1t::HGCalClusterT< C >::mipPt ( ) const
inline
template<class C>
bool l1t::HGCalClusterT< C >::operator< ( const HGCalClusterT< C > &  cl) const
inline

Definition at line 183 of file HGCalClusterT.h.

183 { return mipPt() < cl.mipPt(); }
tuple cl
Definition: haddnano.py:49
double mipPt() const
Definition: HGCalClusterT.h:90
template<class C>
bool l1t::HGCalClusterT< C >::operator<= ( const HGCalClusterT< C > &  cl) const
inline

Definition at line 185 of file HGCalClusterT.h.

185 { return !(cl > *this); }
tuple cl
Definition: haddnano.py:49
template<class C>
bool l1t::HGCalClusterT< C >::operator> ( const HGCalClusterT< C > &  cl) const
inline

Definition at line 184 of file HGCalClusterT.h.

184 { return cl < *this; }
tuple cl
Definition: haddnano.py:49
template<class C>
bool l1t::HGCalClusterT< C >::operator>= ( const HGCalClusterT< C > &  cl) const
inline

Definition at line 186 of file HGCalClusterT.h.

186 { return !(cl < *this); }
tuple cl
Definition: haddnano.py:49
template<class C>
const GlobalPoint& l1t::HGCalClusterT< C >::position ( ) const
inline

Definition at line 99 of file HGCalClusterT.h.

99 { return centre_; }
GlobalPoint centre_
template<class C>
void l1t::HGCalClusterT< C >::removeConstituent ( const edm::Ptr< C > &  c,
bool  updateCentre = true 
)
inline

Definition at line 72 of file HGCalClusterT.h.

Referenced by HGCalClusteringImpl::removeUnconnectedTCinCluster().

72  {
73  /* remove the pointer to c from the std::vector */
74  double fraction = 0;
75  const auto& constituent_itr = constituents_.find(c->detId());
76  const auto& fraction_itr = constituentsFraction_.find(c->detId());
77  if (constituent_itr != constituents_.end()) {
78  // remove constituent and get its fraction in the cluster
79  fraction = fraction_itr->second;
80  constituents_.erase(constituent_itr);
81  constituentsFraction_.erase(fraction_itr);
82 
83  updateP4AndPosition(c, updateCentre, -fraction);
84  }
85  }
std::unordered_map< uint32_t, edm::Ptr< C > > constituents_
void updateP4AndPosition(const edm::Ptr< C > &c, bool updateCentre=true, float fraction=1.)
std::unordered_map< uint32_t, double > constituentsFraction_
template<class C>
double l1t::HGCalClusterT< C >::seedMipPt ( ) const
inline

Definition at line 91 of file HGCalClusterT.h.

91 { return seedMipPt_; }
template<class C>
void l1t::HGCalClusterT< C >::setDetId ( uint32_t  id)
inline

Definition at line 93 of file HGCalClusterT.h.

Referenced by CaloTruthCellsProducer::produce().

93 { detId_ = id; }
uint16_t *__restrict__ id
template<class C>
void l1t::HGCalClusterT< C >::setPt ( double  pt)
inline

Definition at line 94 of file HGCalClusterT.h.

Referenced by HGCalClusteringDummyImpl::calibratePt().

double pt() const final
transverse momentum
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
double mass() const final
mass
double phi() const final
momentum azimuthal angle
void setP4(const LorentzVector &p4) final
set 4-momentum
double eta() const final
momentum pseudorapidity
template<class C>
void l1t::HGCalClusterT< C >::setValid ( bool  valid)
inline

Definition at line 88 of file HGCalClusterT.h.

88 { valid_ = valid; }
bool valid() const
Definition: HGCalClusterT.h:87
template<class C>
int l1t::HGCalClusterT< C >::showerLength ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::showerLength ( int  showerLength)
inline
template<class C>
float l1t::HGCalClusterT< C >::sigmaEtaEtaMax ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaEtaEtaMax ( float  sigmaEtaEtaMax)
inline
template<class C>
float l1t::HGCalClusterT< C >::sigmaEtaEtaTot ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaEtaEtaTot ( float  sigmaEtaEtaTot)
inline
template<class C>
float l1t::HGCalClusterT< C >::sigmaPhiPhiMax ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaPhiPhiMax ( float  sigmaPhiPhiMax)
inline
template<class C>
float l1t::HGCalClusterT< C >::sigmaPhiPhiTot ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaPhiPhiTot ( float  sigmaPhiPhiTot)
inline
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRMax ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaRRMax ( float  sigmaRRMax)
inline

Definition at line 171 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::sigmaRRMax().

171 { sigmaRRMax_ = sigmaRRMax; }
float sigmaRRMax() const
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRMean ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaRRMean ( float  sigmaRRMean)
inline

Definition at line 173 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::sigmaRRMean().

float sigmaRRMean() const
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRTot ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaRRTot ( float  sigmaRRTot)
inline

Definition at line 172 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::sigmaRRTot().

172 { sigmaRRTot_ = sigmaRRTot; }
float sigmaRRTot() const
template<class C>
float l1t::HGCalClusterT< C >::sigmaZZ ( ) const
inline
template<class C>
void l1t::HGCalClusterT< C >::sigmaZZ ( float  sigmaZZ)
inline

Definition at line 174 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::sigmaZZ().

174 { sigmaZZ_ = sigmaZZ; }
float sigmaZZ() const
template<class C>
unsigned l1t::HGCalClusterT< C >::size ( void  ) const
inline

Definition at line 52 of file HGCalClusterT.h.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

52 { return constituents_.size(); }
std::unordered_map< uint32_t, edm::Ptr< C > > constituents_
template<class C>
uint32_t l1t::HGCalClusterT< C >::subdetId ( ) const
inline

Definition at line 138 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::hOverE().

138 { return detId_.subdetId(); }
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
template<class C>
double l1t::HGCalClusterT< C >::sumPt ( ) const
inline
template<class C>
float l1t::HGCalClusterT< C >::triggerCells67percent ( ) const
inline

Definition at line 159 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

159 { return triggerCells67percent_; }
template<class C>
void l1t::HGCalClusterT< C >::triggerCells67percent ( float  triggerCells67percent)
inline
template<class C>
float l1t::HGCalClusterT< C >::triggerCells90percent ( ) const
inline

Definition at line 160 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

160 { return triggerCells90percent_; }
template<class C>
void l1t::HGCalClusterT< C >::triggerCells90percent ( float  triggerCells90percent)
inline
template<class C>
void l1t::HGCalClusterT< C >::updateP4AndPosition ( const edm::Ptr< C > &  c,
bool  updateCentre = true,
float  fraction = 1. 
)
inlineprivate

Definition at line 224 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::addConstituent(), and l1t::HGCalClusterT< l1t::HGCalCluster >::removeConstituent().

224  {
225  double cMipt = c->mipPt() * fraction;
226  double cPt = c->pt() * fraction;
227  /* update cluster positions (IF requested) */
228  if (updateCentre) {
229  Basic3DVector<float> constituentCentre(c->position());
230  Basic3DVector<float> clusterCentre(centre_);
231 
232  clusterCentre = clusterCentre * mipPt_ + constituentCentre * cMipt;
233  if ((mipPt_ + cMipt) > 0) {
234  clusterCentre /= (mipPt_ + cMipt);
235  }
236  centre_ = GlobalPoint(clusterCentre);
237 
238  if (clusterCentre.z() != 0) {
239  centreProj_ = GlobalPoint(clusterCentre / std::abs(clusterCentre.z()));
240  }
241  }
242 
243  /* update cluster energies */
244  mipPt_ += cMipt;
245  sumPt_ += cPt;
246  int updatedPt = hwPt() + (int)(c->hwPt() * fraction);
247  setHwPt(updatedPt);
248 
249  math::PtEtaPhiMLorentzVector updatedP4(p4());
250  updatedP4 += (c->p4() * fraction);
251  setP4(updatedP4);
252  }
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
const LorentzVector & p4() const final
four-momentum Lorentz vector
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
GlobalPoint centreProj_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
GlobalPoint centre_
int hwPt() const
Definition: L1Candidate.h:35
void setHwPt(int pt)
Definition: L1Candidate.h:28
void setP4(const LorentzVector &p4) final
set 4-momentum
template<class C>
bool l1t::HGCalClusterT< C >::valid ( ) const
inline

Definition at line 87 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::setValid().

87 { return valid_; }
template<class C>
float l1t::HGCalClusterT< C >::zBarycenter ( ) const
inline

Definition at line 155 of file HGCalClusterT.h.

Referenced by HGCalShowerShape::fillShapes().

155 { return zBarycenter_; }
template<class C>
void l1t::HGCalClusterT< C >::zBarycenter ( float  zBarycenter)
inline

Definition at line 175 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::zBarycenter().

float zBarycenter() const

Member Data Documentation

template<class C>
GlobalPoint l1t::HGCalClusterT< C >::centre_
private
template<class C>
GlobalPoint l1t::HGCalClusterT< C >::centreProj_
private
template<class C>
std::unordered_map<uint32_t, edm::Ptr<C> > l1t::HGCalClusterT< C >::constituents_
private
template<class C>
std::unordered_map<uint32_t, double> l1t::HGCalClusterT< C >::constituentsFraction_
private
template<class C>
int l1t::HGCalClusterT< C >::coreShowerLength_ = 0
private
template<class C>
DetId l1t::HGCalClusterT< C >::detId_
private
template<class C>
float l1t::HGCalClusterT< C >::eMax_ = 0.
private

Definition at line 208 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::eMax().

template<class C>
int l1t::HGCalClusterT< C >::firstLayer_ = 0
private
template<class C>
float l1t::HGCalClusterT< C >::layer10percent_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::layer50percent_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::layer90percent_ = 0.
private
template<class C>
int l1t::HGCalClusterT< C >::maxLayer_ = 0
private

Definition at line 207 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::maxLayer().

template<class C>
double l1t::HGCalClusterT< C >::mipPt_ = 0.
private
template<class C>
double l1t::HGCalClusterT< C >::seedMipPt_ = 0.
private
template<class C>
int l1t::HGCalClusterT< C >::showerLength_ = 0
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaEtaEtaMax_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaEtaEtaTot_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaPhiPhiMax_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaPhiPhiTot_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRMax_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRMean_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaRRTot_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::sigmaZZ_ = 0.
private

Definition at line 216 of file HGCalClusterT.h.

Referenced by l1t::HGCalClusterT< l1t::HGCalCluster >::sigmaZZ().

template<class C>
double l1t::HGCalClusterT< C >::sumPt_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::triggerCells67percent_ = 0.
private
template<class C>
float l1t::HGCalClusterT< C >::triggerCells90percent_ = 0.
private
template<class C>
bool l1t::HGCalClusterT< C >::valid_ = false
private
template<class C>
float l1t::HGCalClusterT< C >::zBarycenter_ = 0.
private