CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PackedGenParticle.cc
Go to the documentation of this file.
4 
5 
6 
7 void pat::PackedGenParticle::pack(bool unpackAfterwards) {
9  packedY_ = int16_t(p4_.Rapidity()/6.0f*std::numeric_limits<int16_t>::max());
10  packedPhi_ = int16_t(p4_.Phi()/3.2f*std::numeric_limits<int16_t>::max());
12  if (unpackAfterwards) unpack(); // force the values to match with the packed ones
13 }
14 
15 
17  float y = int16_t(packedY_)*6.0f/std::numeric_limits<int16_t>::max();
18  float pt=MiniFloatConverter::float16to32(packedPt_);
19  float m=MiniFloatConverter::float16to32(packedM_);
20  float pz = std::tanh(y)*std::sqrt((m*m+pt*pt)/(1.-std::tanh(y)*std::tanh(y)));
21  float eta = std::asinh(pz/pt);
22  double shift = (pt<1. ? 0.1*pt : 0.1/pt); // shift particle phi to break degeneracies in angular separations
23  double sign = ( ( int(pt*10) % 2 == 0 ) ? 1 : -1 ); // introduce a pseudo-random sign of the shift
24  double phi = int16_t(packedPhi_)*3.2f/std::numeric_limits<int16_t>::max() + sign*shift*3.2/std::numeric_limits<int16_t>::max();
25  p4_ = PolarLorentzVector(pt,eta,phi,m);
26  p4c_ = p4_;
27  unpacked_ = true;
28 }
29 
31 
32 
33 float pat::PackedGenParticle::dxy(const Point &p) const {
34  unpack();
35  return -(vertex_.X()-p.X()) * std::sin(float(p4_.Phi())) + (vertex_.Y()-p.Y()) * std::cos(float(p4_.Phi()));
36 }
37 float pat::PackedGenParticle::dz(const Point &p) const {
38  unpack();
39  return (vertex_.Z()-p.X()) - ((vertex_.X()-p.X()) * std::cos(float(p4_.Phi())) + (vertex_.Y()-p.Y()) * std::sin(float(p4_.Phi()))) * p4_.Pz()/p4_.Pt();
40 }
41 
42 
44 
46  throw cms::Exception("Invalid Reference")
47  << "this Candidate has no master clone reference."
48  << "Can't call masterClone() method.\n";
49 }
50 
52  return false;
53 }
54 
56  return false;
57 }
58 
59 
61  throw cms::Exception("Invalid Reference")
62  << "this Candidate has no master clone ptr."
63  << "Can't call masterClonePtr() method.\n";
64 }
65 
67  return 0;
68 }
69 
71  if(mother_.isNonnull()) return 1;
72  return 0;
73 }
74 
76  return p4() == o.p4() && vertex() == o.vertex() && charge() == o.charge();
77 // return p4() == o.p4() && charge() == o.charge();
78 }
79 
81  return 0;
82 }
83 
85  return mother_.get();
86 }
87 
90  << "This Candidate type does not implement daughter(std::string). "
91  << "Please use CompositeCandidate or NamedCompositeCandidate.\n";
92 }
93 
96  << "This Candidate type does not implement daughter(std::string). "
97  << "Please use CompositeCandidate or NamedCompositeCandidate.\n";
98 }
99 
100 
101 
103  return 0;
104 }
105 
107  return 0;
108 }
109 
111  return 0;
112 }
113 
115  return 0;
116 }
117 
120  << "reco::ConcreteCandidate does not implement vertex covariant matrix.\n";
121 }
122 
125  << "reco::ConcreteCandidate does not implement vertex covariant matrix.\n";
126 }
127 
128 bool pat::PackedGenParticle::isElectron() const { return false; }
129 
130 bool pat::PackedGenParticle::isMuon() const { return false; }
131 
132 bool pat::PackedGenParticle::isGlobalMuon() const { return false; }
133 
134 bool pat::PackedGenParticle::isStandAloneMuon() const { return false; }
135 
136 bool pat::PackedGenParticle::isTrackerMuon() const { return false; }
137 
138 bool pat::PackedGenParticle::isCaloMuon() const { return false; }
139 
140 bool pat::PackedGenParticle::isPhoton() const { return false; }
141 
142 bool pat::PackedGenParticle::isConvertedPhoton() const { return false; }
143 
144 bool pat::PackedGenParticle::isJet() const { return false; }
145 
146 bool pat::PackedGenParticle::longLived() const {return false;}
147 
148 bool pat::PackedGenParticle::massConstraint() const {return false;}
149 
150 
151 
152 
int i
Definition: DBlmapReader.cc:9
virtual void fillVertexCovariance(CovarianceMatrix &v) const
fill SMatrix
CovarianceMatrix vertexCovariance() const
return SMatrix
virtual size_t numberOfMothers() const
number of mothers
size_t size_type
Definition: Candidate.h:30
virtual double vertexNdof() const
double sign(double x)
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
virtual bool isTrackerMuon() const
virtual bool isJet() const
virtual ~PackedGenParticle()
destructor
virtual const reco::Candidate * daughter(size_type) const
return daughter at a given position (throws an exception)
T eta() const
virtual bool isElectron() const
get a component
virtual bool isPhoton() const
static float float16to32(uint16_t h)
Definition: libminifloat.h:10
PolarLorentzVector p4_
the four vector
void pack(bool unpackAfterwards=true)
virtual double vertexNormalizedChi2() const
chi-squared divided by n.d.o.f.
virtual bool isMuon() const
virtual bool longLived() const
is long lived?
T sqrt(T t)
Definition: SSEVec.h:48
double p4[4]
Definition: TauolaWrapper.h:92
static uint16_t float32to16(float x)
Definition: libminifloat.h:15
virtual const reco::CandidateBaseRef & masterClone() const
virtual double vertexChi2() const
chi-squares
virtual size_t numberOfDaughters() const
number of daughters
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
virtual bool isStandAloneMuon() const
int j
Definition: DBlmapReader.cc:9
virtual const reco::CandidatePtr & masterClonePtr() const
virtual const Point & vertex() const =0
vertex position
virtual int charge() const =0
electric charge
virtual bool hasMasterClonePtr() const
virtual bool overlap(const reco::Candidate &) const
check overlap with another Candidate
virtual bool hasMasterClone() const
virtual bool isCaloMuon() const
virtual bool isGlobalMuon() const
virtual const reco::Candidate * mother(size_type) const
return mother at a given position (throws an exception)
virtual bool isConvertedPhoton() const
T get() const
get a component
Definition: Candidate.h:217
static unsigned int const shift
math::XYZPoint Point
point in the space
Definition: Candidate.h:41
virtual bool massConstraint() const
do mass constraint?
virtual float dz() const
dz with respect to the PV ref
virtual float dxy() const
dxy with respect to the PV ref
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
Definition: DDAxes.h:10
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
Definition: Candidate.h:39