Belle II Software light-2406-ragdoll
InternalParticle Class Reference

another unnecessary layer of abstraction More...

#include <InternalParticle.h>

Inheritance diagram for InternalParticle:
Collaboration diagram for InternalParticle:

Public Types

enum  TFParticleType {
  kInteractionPoint ,
  kOrigin ,
  kComposite ,
  kRecoResonance ,
  kInternalParticle ,
  kRecoTrack ,
  kResonance ,
  kRecoPhoton ,
  kRecoKlong ,
  kMissingParticle
}
 particle types
More...
 
typedef std::vector< Constraintconstraintlist
 alias
 
typedef std::vector< std::pair< const ParticleBase *, int > > indexmap
 alias
 

Public Member Functions

 InternalParticle (Belle2::Particle *particle, const ParticleBase *mother, const ConstraintConfiguration &config, bool forceFitAll)
 constructor
 
virtual ~InternalParticle ()
 destructor
 
virtual ErrCode initCovariance (FitParams &) const override
 init covariance
 
ErrCode projectKineConstraint (const FitParams &, Projection &) const
 project kinematical constraint
 
ErrCode projectBeamConstraint (const FitParams &, Projection &) const
 project beam four momentum constraint
 
virtual void forceP4Sum (FitParams &) const override
 enforce conservation of momentum sum
 
virtual ErrCode initParticleWithMother (FitParams &fitparams) override
 init particle in case it has a mother
 
virtual ErrCode initMotherlessParticle (FitParams &fitparams) override
 init particle in case it has no mother
 
ErrCode projectConstraint (const Constraint::Type type, const FitParams &fitparams, Projection &p) const override
 find out which constraint it is and project
 
virtual int dim () const override
 space reserved in fit params, if has mother then it has tau
 
virtual int type () const override
 type
 
virtual int posIndex () const override
 position index in fit params
 
virtual int tauIndex () const override
 tau index in fit params only if it has a mother
 
virtual int momIndex () const override
 momentum index in fit params depending on whether it has a mother

 
virtual bool hasEnergy () const override
 has energy in fitparams

 
virtual bool hasPosition () const override
 has position index

 
virtual std::string parname (int index) const override
 name

 
virtual void addToConstraintList (constraintlist &list, int depth) const override
 add to constraint list

 
void setMassConstraint (bool b)
 set mass constraint flag
 
virtual void updateIndex (int &offset)
 this sets the index for momentum, position, etc.
 
const ParticleBaselocate (Belle2::Particle *particle) const
 get particle base from basf2 particle
 
Belle2::Particleparticle () const
 get basf2 particle

 
int index () const
 get index

 
const ParticleBasemother () const
 getMother() / hasMother()
 
virtual ErrCode projectGeoConstraint (const FitParams &, Projection &) const
 project geometrical constraint
 
virtual ErrCode projectMassConstraintParticle (const FitParams &, Projection &) const
 project mass constraint using the particles parameters
 
virtual ErrCode projectMassConstraintDaughters (const FitParams &, Projection &) const
 project mass constraint using the parameters of the daughters
 
virtual ErrCode projectMassConstraint (const FitParams &, Projection &) const
 project mass constraint abstract
 
int eneIndex () const
 get energy index

 
virtual double chiSquare (const FitParams &) const
 get chi2
 
int charge () const
 get charge
 
virtual ParticleBaseaddDaughter (Belle2::Particle *, const ConstraintConfiguration &config, bool forceFitAll=false)
 add daughter

 
virtual void removeDaughter (const ParticleBase *pb)
 remove daughter
 
virtual void retrieveIndexMap (indexmap &anindexmap) const
 get index map

 
void setMother (const ParticleBase *m)
 set mother

 
void collectVertexDaughters (std::vector< ParticleBase * > &particles, int posindex)
 get vertex daughters
 
virtual int nFinalChargedCandidates () const
 number of charged candidates
 

Static Public Member Functions

static ParticleBasecreateParticle (Belle2::Particle *particle, const ParticleBase *mother, const ConstraintConfiguration &config, bool forceFitAll=false)
 create the according treeFitter particle obj for a basf2 particle type

 
static ParticleBasecreateOrigin (Belle2::Particle *daughter, const ConstraintConfiguration &config, bool forceFitAll)
 create a custom origin particle or a beamspot
 

Protected Types

typedef std::vector< ParticleBase * > ParticleContainer
 just an alias
 

Protected Member Functions

ErrCode initMomentum (FitParams &fitparams) const
 init momentum of *this and daughters
 
ErrCode initTau (FitParams &par) const
 initialises tau as a length

 
void setIndex (int i)
 set Index (in statevector)
 

Static Protected Member Functions

static bool isAResonance (Belle2::Particle *particle)
 controls if a particle is treated as a resonance(lifetime=0) or a particle that has a finite lifetime.
 

Protected Attributes

Belle2::Particlem_particle
 pointer to framework type

 
const ParticleBasem_mother
 motherparticle
 
std::vector< ParticleBase * > m_daughters
 daughter container

 
bool m_isStronglyDecayingResonance
 decay length less than 1 micron

 
const ConstraintConfigurationm_config
 has all the constraint config
 

Static Private Member Functions

static bool compTrkTransverseMomentum (const RecoTrack *lhs, const RecoTrack *rhs)
 compare transverse track momentum
 

Private Attributes

bool m_massconstraint
 has mass constraint
 
bool m_beamconstraint
 has beam constraint
 
bool m_shares_vertex_with_mother
 shares vertex with mother, that means decay vertex = productionvertex
 
bool m_geo_constraint
 use a geo metric constraint
 
bool m_lifetimeconstraint
 has lifetime constraint

 
bool m_isconversion
 is conversion

 
bool m_automatic_vertex_constraining
 automatically figure out if mother and particle vertex should be the same and also add geometric constraints
 
int m_index
 index
 
std::string m_name
 name

 

Detailed Description

another unnecessary layer of abstraction

Definition at line 17 of file InternalParticle.h.

Member Typedef Documentation

◆ constraintlist

typedef std::vector<Constraint> constraintlist
inherited

alias

Definition at line 52 of file ParticleBase.h.

◆ indexmap

typedef std::vector< std::pair<const ParticleBase*, int> > indexmap
inherited

alias

Definition at line 55 of file ParticleBase.h.

◆ ParticleContainer

typedef std::vector<ParticleBase*> ParticleContainer
protectedinherited

just an alias

Definition at line 178 of file ParticleBase.h.

Member Enumeration Documentation

◆ TFParticleType

enum TFParticleType
inherited

particle types

Definition at line 30 of file ParticleBase.h.

30 {kInteractionPoint,
31 kOrigin,
32 kComposite,
33 kRecoResonance,
34 kInternalParticle,
35 kRecoTrack,
36 kResonance,
37 kRecoPhoton,
38 kRecoKlong,
39 kMissingParticle
40 };

Constructor & Destructor Documentation

◆ InternalParticle()

InternalParticle ( Belle2::Particle particle,
const ParticleBase mother,
const ConstraintConfiguration config,
bool  forceFitAll 
)

constructor

Definition at line 41 of file InternalParticle.cc.

45 :
46 ParticleBase(particle, mother, &config),// config pointer here to allow final states not to have it
47 m_massconstraint(false),
48 m_beamconstraint(false),
50 m_isconversion(false),
51 m_automatic_vertex_constraining(config.m_automatic_vertex_constraining)
52 {
53 if (particle) {
54 for (auto daughter : particle->getDaughters()) {
55 addDaughter(daughter, config, forceFitAll);
56 }
57 } else {
58 B2ERROR("Trying to create an InternalParticle from NULL. This should never happen.");
59 }
60
61 m_massconstraint = std::find(config.m_massConstraintListPDG.begin(), config.m_massConstraintListPDG.end(),
62 std::abs(m_particle->getPDGCode())) != config.m_massConstraintListPDG.end();
63
64 m_beamconstraint = (std::abs(m_particle->getPDGCode()) == config.m_beamConstraintPDG);
65
67 // if this is a hadronically decaying resonance it is useful to constraint the decay vertex to its mothers decay vertex.
68 //
69 m_shares_vertex_with_mother = std::find(config.m_fixedToMotherVertexListPDG.begin(),
70 config.m_fixedToMotherVertexListPDG.end(),
71 std::abs(m_particle->getPDGCode())) != config.m_fixedToMotherVertexListPDG.end() && this->mother();
72
73 // use geo constraint if this particle is in the list to constrain
74 m_geo_constraint = std::find(config.m_geoConstraintListPDG.begin(),
75 config.m_geoConstraintListPDG.end(),
76 std::abs(m_particle->getPDGCode())) != config.m_geoConstraintListPDG.end() && this->mother() && !m_shares_vertex_with_mother;
77 } else {
80 }
81 }
int getPDGCode(void) const
Returns PDG code.
Definition: Particle.h:454
std::vector< Belle2::Particle * > getDaughters() const
Returns a vector of pointers to daughter particles.
Definition: Particle.cc:637
bool m_isconversion
is conversion
bool m_beamconstraint
has beam constraint
bool m_massconstraint
has mass constraint
bool m_automatic_vertex_constraining
automatically figure out if mother and particle vertex should be the same and also add geometric cons...
bool m_lifetimeconstraint
has lifetime constraint
bool m_geo_constraint
use a geo metric constraint
bool m_shares_vertex_with_mother
shares vertex with mother, that means decay vertex = productionvertex
Belle2::Particle * particle() const
get basf2 particle
Definition: ParticleBase.h:92
virtual ParticleBase * addDaughter(Belle2::Particle *, const ConstraintConfiguration &config, bool forceFitAll=false)
add daughter
Definition: ParticleBase.cc:65
bool m_isStronglyDecayingResonance
decay length less than 1 micron
Definition: ParticleBase.h:201
Belle2::Particle * m_particle
pointer to framework type
Definition: ParticleBase.h:192
ParticleBase(Belle2::Particle *particle, const ParticleBase *mother, const ConstraintConfiguration *config=nullptr)
default constructor
Definition: ParticleBase.cc:30
const ParticleBase * mother() const
getMother() / hasMother()
Definition: ParticleBase.h:98

◆ ~InternalParticle()

virtual ~InternalParticle ( )
inlinevirtual

destructor

Definition at line 29 of file InternalParticle.h.

29{};

Member Function Documentation

◆ addDaughter()

ParticleBase * addDaughter ( Belle2::Particle cand,
const ConstraintConfiguration config,
bool  forceFitAll = false 
)
virtualinherited

add daughter

Definition at line 65 of file ParticleBase.cc.

66 {
67 auto newDaughter = ParticleBase::createParticle(cand, this, config, forceFitAll);
68 m_daughters.push_back(newDaughter);
69 return m_daughters.back();
70 }
static ParticleBase * createParticle(Belle2::Particle *particle, const ParticleBase *mother, const ConstraintConfiguration &config, bool forceFitAll=false)
create the according treeFitter particle obj for a basf2 particle type
std::vector< ParticleBase * > m_daughters
daughter container
Definition: ParticleBase.h:198

◆ addToConstraintList()

void addToConstraintList ( constraintlist list,
int  depth 
) const
overridevirtual

add to constraint list

Implements ParticleBase.

Definition at line 376 of file InternalParticle.cc.

378 {
379
380 for (auto daughter : m_daughters) {
381 daughter->addToConstraintList(list, depth - 1);
382 }
383 if (tauIndex() >= 0 && m_lifetimeconstraint) {
384 list.push_back(Constraint(this, Constraint::lifetime, depth, 1));
385 }
386 if (momIndex() >= 0) {
387 list.push_back(Constraint(this, Constraint::kinematic, depth, 4, 3));
388 }
389 if (m_geo_constraint) {
390 assert(m_config);
391 const int dim = m_config->m_originDimension == 2 && std::abs(m_particle->getPDGCode()) == m_config->m_headOfTreePDG ? 2 : 3;
392 list.push_back(Constraint(this, Constraint::geometric, depth, dim, 3));
393 }
394 if (m_massconstraint) {
395 list.push_back(Constraint(this, Constraint::mass, depth, 1, 3));
396 }
397 if (m_beamconstraint) {
398 assert(m_config);
399 list.push_back(Constraint(this, Constraint::beam, depth, 4, 3));
400 }
401
402 }
const int m_originDimension
dimension of the origin constraint and ALL geometric gcosntraints
int m_headOfTreePDG
PDG code of the head particle.
virtual int tauIndex() const override
tau index in fit params only if it has a mother
virtual int dim() const override
space reserved in fit params, if has mother then it has tau
virtual int momIndex() const override
momentum index in fit params depending on whether it has a mother
const ConstraintConfiguration * m_config
has all the constraint config
Definition: ParticleBase.h:204

◆ charge()

int charge ( ) const
inlineinherited

get charge

Definition at line 144 of file ParticleBase.h.

145 {
146 if (m_particle->getPDGCode()) {
147 double fltcharge = m_particle->getCharge();
148 return fltcharge < 0 ? int(fltcharge - 0.5) : int(fltcharge + 0.5);
149 } else {
150 return m_particle->getCharge() > 0 ? 1 : (m_particle->getCharge() < 0 ? -1 : 0);
151 }
152 }
double getCharge(void) const
Returns particle charge.
Definition: Particle.cc:622

◆ chiSquare()

double chiSquare ( const FitParams fitparams) const
virtualinherited

get chi2

Definition at line 231 of file ParticleBase.cc.

232 {
233 double rc = 0;
234 for (auto* daughter : m_daughters) {
235 rc += daughter->chiSquare(fitparams);
236 }
237 return rc;
238 }

◆ collectVertexDaughters()

void collectVertexDaughters ( std::vector< ParticleBase * > &  particles,
int  posindex 
)
inherited

get vertex daughters

Definition at line 156 of file ParticleBase.cc.

157 {
158 if (mother() && mother()->posIndex() == posindex) {
159 particles.push_back(this);
160 }
161
162 for (auto* daughter : m_daughters) {
163 daughter->collectVertexDaughters(particles, posindex);
164 }
165 }
virtual int posIndex() const
get vertex index (in statevector!)
Definition: ParticleBase.h:122

◆ compTrkTransverseMomentum()

bool compTrkTransverseMomentum ( const RecoTrack lhs,
const RecoTrack rhs 
)
staticprivate

compare transverse track momentum

Definition at line 83 of file InternalParticle.cc.

84 {
85
86 return lhs->particle()->getMomentum().Rho() > rhs->particle()->getMomentum().Rho();
87 }

◆ createOrigin()

ParticleBase * createOrigin ( Belle2::Particle daughter,
const ConstraintConfiguration config,
bool  forceFitAll 
)
staticinherited

create a custom origin particle or a beamspot

Definition at line 93 of file ParticleBase.cc.

98 {
99 return new Origin(daughter, config, forceFitAll);
100 }

◆ createParticle()

ParticleBase * createParticle ( Belle2::Particle particle,
const ParticleBase mother,
const ConstraintConfiguration config,
bool  forceFitAll = false 
)
staticinherited

create the according treeFitter particle obj for a basf2 particle type

Definition at line 102 of file ParticleBase.cc.

104 {
105 ParticleBase* rc = nullptr;
106
107 if (!mother) { // If there is no mother, this is the 'head of tree' particle (is never a resonance)
108 rc = new InternalParticle(particle, nullptr, config, forceFitAll);
109 } else if (particle->hasExtraInfo("bremsCorrected") // Has Bremsstrahlungs-recovery
110 && particle->getExtraInfo("bremsCorrected") != 0) { // and gammas are attached
111 rc = new Composite(particle, mother, config, true);
112 // if no gamma is attached, it is treated as a RecoTrack
113 } else if (particle->hasExtraInfo("treeFitterTreatMeAsInvisible")
114 && particle->getExtraInfo("treeFitterTreatMeAsInvisible") == 1) { // dummy particles with invisible flag
115 rc = new RecoResonance(particle, mother, config);
116 } else if (particle->getTrack()) { // external reconstructed track
117 rc = new RecoTrack(particle, mother);
118 } else if (particle->getECLCluster()) { // external reconstructed photon
119 rc = new RecoPhoton(particle, mother);
120 } else if (particle->getKLMCluster()) { // external reconstructed klong
121 rc = new RecoKlong(particle, mother);
122 } else if (particle->getMdstArrayIndex()) { // external composite e.g. V0
123 rc = new InternalParticle(particle, mother, config, forceFitAll);
124 } else { // 'internal' particles
125 if (isAResonance(particle)) {
126 rc = new Resonance(particle, mother, config, forceFitAll);
127 } else {
128 rc = new InternalParticle(particle, mother, config, forceFitAll);
129 }
130 }
131 return rc;
132 }
const KLMCluster * getKLMCluster() const
Returns the pointer to the KLMCluster object that was used to create this Particle (ParticleType == c...
Definition: Particle.cc:926
const Track * getTrack() const
Returns the pointer to the Track object that was used to create this Particle (ParticleType == c_Trac...
Definition: Particle.cc:845
const ECLCluster * getECLCluster() const
Returns the pointer to the ECLCluster object that was used to create this Particle (if ParticleType =...
Definition: Particle.cc:891
bool hasExtraInfo(const std::string &name) const
Return whether the extra info with the given name is set.
Definition: Particle.cc:1266
unsigned getMdstArrayIndex(void) const
Returns 0-based index of MDST store array object (0 for composite particles)
Definition: Particle.h:487
double getExtraInfo(const std::string &name) const
Return given value if set.
Definition: Particle.cc:1289
static bool isAResonance(Belle2::Particle *particle)
controls if a particle is treated as a resonance(lifetime=0) or a particle that has a finite lifetime...

◆ dim()

int dim ( ) const
overridevirtual

space reserved in fit params, if has mother then it has tau

Implements ParticleBase.

Reimplemented in Resonance.

Definition at line 336 of file InternalParticle.cc.

337 {
338 // { x, y, z, tau, px, py, pz, E }
339 // the last 4 always exists for composite particles
340 // tau index only exist with a vertex and geo constraint
341 //
342 if (m_shares_vertex_with_mother) { return 4; }
343 else if (!m_geo_constraint) { return 7; }
344 else { return 8; }
345 }

◆ eneIndex()

int eneIndex ( ) const
inlineinherited

get energy index

Definition at line 138 of file ParticleBase.h.

138{ return hasEnergy() ? momIndex() + 3 : -1 ; }
virtual int momIndex() const
get momentum index
Definition: ParticleBase.h:128
virtual bool hasEnergy() const
get momentum dimension
Definition: ParticleBase.h:132

◆ forceP4Sum()

void forceP4Sum ( FitParams fitparams) const
overridevirtual

enforce conservation of momentum sum

Reimplemented from ParticleBase.

Definition at line 411 of file InternalParticle.cc.

412 {
413 for (const auto daughter : m_daughters) {
414 daughter->forceP4Sum(fitparams);
415 }
416 const int momindex = momIndex();
417 if (momindex > 0) {
418 const int dim = hasEnergy() ? 4 : 3;
419 Projection p(fitparams.getDimensionOfState(), dim);
420 projectKineConstraint(fitparams, p);
421 fitparams.getStateVector().segment(momindex, dim) -= p.getResiduals().segment(0, dim);
422 }
423 }
ErrCode projectKineConstraint(const FitParams &, Projection &) const
project kinematical constraint
virtual bool hasEnergy() const override
has energy in fitparams

◆ hasEnergy()

virtual bool hasEnergy ( ) const
inlineoverridevirtual

has energy in fitparams

Reimplemented from ParticleBase.

Definition at line 68 of file InternalParticle.h.

68{ return true ; }

◆ hasPosition()

bool hasPosition ( ) const
overridevirtual

has position index

Reimplemented from ParticleBase.

Reimplemented in Resonance.

Definition at line 369 of file InternalParticle.cc.

370 {
371 // does this particle have a position index (decayvertex)
372 // true in any case
373 return true;
374 }

◆ index()

int index ( ) const
inlineinherited

get index

Definition at line 95 of file ParticleBase.h.

95{ return m_index ; }

◆ initCovariance()

ErrCode initCovariance ( FitParams fitparams) const
overridevirtual

init covariance

Reimplemented from ParticleBase.

Definition at line 226 of file InternalParticle.cc.

227 {
228 ErrCode status;
230 for (auto daughter : m_daughters) {
231 status |= daughter->initCovariance(fitparams);
232 }
233 return status;
234 }
virtual ErrCode initCovariance(FitParams &) const
init covariance matrix

◆ initMomentum()

ErrCode initMomentum ( FitParams fitparams) const
protected

init momentum of *this and daughters

Definition at line 206 of file InternalParticle.cc.

207 {
208 int momindex = momIndex();
209 fitparams.getStateVector().segment(momindex, 4) = Eigen::Matrix<double, 4, 1>::Zero(4);
210
211 for (auto daughter : m_daughters) {
212 int daumomindex = daughter->momIndex();
213 int maxrow = daughter->hasEnergy() ? 4 : 3;
214
215 double e2 = fitparams.getStateVector().segment(daumomindex, maxrow).squaredNorm();
216 fitparams.getStateVector().segment(momindex, maxrow) += fitparams.getStateVector().segment(daumomindex, maxrow);
217
218 if (maxrow == 3) {
219 double mass = daughter->particle()->getPDGMass();
220 fitparams.getStateVector()(momindex + 3) += std::sqrt(e2 + mass * mass);
221 }
222 }
223 return ErrCode(ErrCode::Status::success);
224 }

◆ initMotherlessParticle()

ErrCode initMotherlessParticle ( FitParams fitparams)
overridevirtual

init particle in case it has no mother

(0,0,0) is the best guess in any other case

Implements ParticleBase.

Reimplemented in Resonance.

Definition at line 89 of file InternalParticle.cc.

90 {
91 ErrCode status ;
92 for (auto daughter : m_daughters) {
93 status |= daughter->initMotherlessParticle(fitparams);
94 }
95
96 int posindex = posIndex();
97 // FIXME update this and check if this position is already initialised
98 // lower stream vertices might be better
99 if (hasPosition()) {
100 fitparams.getStateVector().segment(posindex, 3) = Eigen::Matrix<double, 3, 1>::Zero(3);
101
102 std::vector<ParticleBase*> alldaughters;
103 ParticleBase::collectVertexDaughters(alldaughters, posindex);
104
105 std::vector<ParticleBase*> vtxdaughters;
106
107 vector<RecoTrack*> trkdaughters;
108 for (auto daughter : alldaughters) {
109 if (daughter->type() == ParticleBase::TFParticleType::kRecoTrack) {
110 trkdaughters.push_back(static_cast<RecoTrack*>(daughter));
111 } else if (daughter->hasPosition()
112 && fitparams.getStateVector()(daughter->posIndex()) != 0) {
113 vtxdaughters.push_back(daughter);
114 }
115 }
116
117 if (trkdaughters.size() >= 2) {
118
119 auto v0 = particle()->getV0();
120 auto dummy_vertex = ROOT::Math::XYZVector(0, 0, 0);
121
122 bool initWithV0 = false;
123 if (v0 && v0->getFittedVertexPosition() != dummy_vertex) {
124 auto part_dau1 = particle()->getDaughter(0);
125 auto part_dau2 = particle()->getDaughter(1);
126
127 auto recotrack_dau1 = std::find_if(trkdaughters.begin(), trkdaughters.end(),
128 [&part_dau1](RecoTrack * rt) { return rt->particle()->getMdstArrayIndex() == part_dau1->getMdstArrayIndex(); });
129 auto recotrack_dau2 = std::find_if(trkdaughters.begin(), trkdaughters.end(),
130 [&part_dau2](RecoTrack * rt) { return rt->particle()->getMdstArrayIndex() == part_dau2->getMdstArrayIndex(); });
131
132 if (recotrack_dau1 == trkdaughters.end() || recotrack_dau2 == trkdaughters.end()) {
133 B2WARNING("V0 daughter particles do not match with RecoTracks.");
134 } else {
135 double X_V0(v0->getFittedVertexX()), Y_V0(v0->getFittedVertexY()), Z_V0(v0->getFittedVertexZ());
136 fitparams.getStateVector()(posindex) = X_V0;
137 fitparams.getStateVector()(posindex + 1) = Y_V0;
138 fitparams.getStateVector()(posindex + 2) = Z_V0;
139
140 Belle2::Helix helix1 = v0->getTrackFitResults().first->getHelix();
141 Belle2::Helix helix2 = v0->getTrackFitResults().second->getHelix();
142
143 (*recotrack_dau1)->setFlightLength(helix1.getArcLength2DAtXY(X_V0, Y_V0));
144 (*recotrack_dau2)->setFlightLength(helix2.getArcLength2DAtXY(X_V0, Y_V0));
145
146 initWithV0 = true;
147 }
148 }
149
150 if (!initWithV0) {
151 std::sort(trkdaughters.begin(), trkdaughters.end(), compTrkTransverseMomentum);
152
153 RecoTrack* dau1 = trkdaughters[0];
154 RecoTrack* dau2 = trkdaughters[1];
155
156 Belle2::Helix helix1 = dau1->particle()->getTrackFitResult()->getHelix();
157 Belle2::Helix helix2 = dau2->particle()->getTrackFitResult()->getHelix();
158
159 double flt1(0), flt2(0);
161 HelixUtils::helixPoca(helix1, helix2, flt1, flt2, v, m_isconversion);
162
163 fitparams.getStateVector()(posindex) = v.X();
164 fitparams.getStateVector()(posindex + 1) = v.Y();
165 fitparams.getStateVector()(posindex + 2) = v.Z();
166
167 dau1->setFlightLength(flt1);
168 dau2->setFlightLength(flt2);
169 }
170
171 } else if (false && trkdaughters.size() + vtxdaughters.size() >= 2) {
172 // TODO switched off waiting for refactoring of init1 and init2 functions (does not affect performance)
173 } else if (mother() && mother()->posIndex() >= 0) {
174 const int posindexmother = mother()->posIndex();
175 const int dim = m_config->m_originDimension; //TODO access mother
176 fitparams.getStateVector().segment(posindex, dim) = fitparams.getStateVector().segment(posindexmother, dim);
177 } else {
179 fitparams.getStateVector().segment(posindex, 3) = Eigen::Matrix<double, 1, 3>::Zero(3);
180 }
181 }
182
183 for (auto daughter : m_daughters) {
184 daughter->initParticleWithMother(fitparams);
185 }
186
187 initMomentum(fitparams);
188 return status;
189 }
DataType Z() const
access variable Z (= .at(2) without boundary check)
Definition: B2Vector3.h:435
DataType X() const
access variable X (= .at(0) without boundary check)
Definition: B2Vector3.h:431
DataType Y() const
access variable Y (= .at(1) without boundary check)
Definition: B2Vector3.h:433
This class represents an ideal helix in perigee parameterization.
Definition: Helix.h:59
double getArcLength2DAtXY(const double &x, const double &y) const
Calculates the two dimensional arc length at which the circle in the xy projection is closest to the ...
Definition: Helix.cc:141
const V0 * getV0() const
Returns the pointer to the V0 object that was used to create this Particle (if ParticleType == c_V0).
Definition: Particle.cc:880
const Particle * getDaughter(unsigned i) const
Returns a pointer to the i-th daughter particle.
Definition: Particle.cc:631
static double helixPoca(const Belle2::Helix &helix1, const Belle2::Helix &helix2, double &flt1, double &flt2, Belle2::B2Vector3D &vertex, bool parallel=false)
POCA between two tracks.
Definition: HelixUtils.cc:214
virtual bool hasPosition() const override
has position index
ErrCode initMomentum(FitParams &fitparams) const
init momentum of *this and daughters
static bool compTrkTransverseMomentum(const RecoTrack *lhs, const RecoTrack *rhs)
compare transverse track momentum
virtual int posIndex() const override
position index in fit params
void collectVertexDaughters(std::vector< ParticleBase * > &particles, int posindex)
get vertex daughters

◆ initParticleWithMother()

ErrCode initParticleWithMother ( FitParams fitparams)
overridevirtual

init particle in case it has a mother

Implements ParticleBase.

Reimplemented in Resonance.

Definition at line 191 of file InternalParticle.cc.

192 {
193 int posindex = posIndex();
194 if (hasPosition() &&
195 mother() &&
196 fitparams.getStateVector()(posindex) == 0 &&
197 fitparams.getStateVector()(posindex + 1) == 0 && \
198 fitparams.getStateVector()(posindex + 2) == 0) {
199 const int posindexmom = mother()->posIndex();
200 const int dim = m_config->m_originDimension; //TODO access mother?
201 fitparams.getStateVector().segment(posindex, dim) = fitparams.getStateVector().segment(posindexmom, dim);
202 }
203 return initTau(fitparams);
204 }
ErrCode initTau(FitParams &par) const
initialises tau as a length

◆ initTau()

ErrCode initTau ( FitParams par) const
protectedinherited

initialises tau as a length

Definition at line 426 of file ParticleBase.cc.

427 {
428 const int tauindex = tauIndex();
429 if (tauindex >= 0 && hasPosition()) {
430
431 const int posindex = posIndex();
432 const int mother_ps_index = mother()->posIndex();
433 const int dim = m_config->m_originDimension; // TODO can we configure this to be particle specific?
434
435 // tau has different meaning depending on the dimension of the constraint
436 // 2-> use x-y projection
437 const Eigen::Matrix < double, 1, -1, 1, 1, 3 > vertex_dist =
438 fitparams.getStateVector().segment(posindex, dim) - fitparams.getStateVector().segment(mother_ps_index, dim);
439 const Eigen::Matrix < double, 1, -1, 1, 1, 3 >
440 mom = fitparams.getStateVector().segment(posindex, dim) - fitparams.getStateVector().segment(mother_ps_index, dim);
441
442 // if an intermediate vertex is not well defined by a track or so it will be initialised with 0
443 // same for the momentum of for example B0, it might be initialised with 0
444 // in those cases use pdg value
445 const double mom_norm = mom.norm();
446 const double dot = std::abs(vertex_dist.dot(mom));
447 const double tau = dot / mom_norm;
448 if (0 == mom_norm || 0 == dot) {
449 const double mass = m_particle->getPDGMass();
450 if (mass > 0)
451 fitparams.getStateVector()(tauindex) = m_particle->getPDGLifetime() * 1e9 * Belle2::Const::speedOfLight / mass;
452 else
453 fitparams.getStateVector()(tauindex) = 0;
454 } else {
455 fitparams.getStateVector()(tauindex) = tau;
456 }
457 }
458
459 return ErrCode(ErrCode::Status::success);
460 }
static const double speedOfLight
[cm/ns]
Definition: Const.h:695
double getPDGMass(void) const
Returns uncertainty on the invariant mass (requires valid momentum error matrix)
Definition: Particle.cc:604
double getPDGLifetime() const
Returns particle nominal lifetime.
Definition: Particle.cc:613
virtual int dim() const =0
get dimension of constraint
virtual int tauIndex() const
get tau index
Definition: ParticleBase.h:125
virtual bool hasPosition() const
get false
Definition: ParticleBase.h:135

◆ isAResonance()

bool isAResonance ( Belle2::Particle particle)
staticprotectedinherited

controls if a particle is treated as a resonance(lifetime=0) or a particle that has a finite lifetime.

A finite life time means it will register a geo constraint for this particle

Definition at line 134 of file ParticleBase.cc.

135 {
136 bool rc = false ;
137 const int pdgcode = std::abs(particle->getPDGCode());
138
139 if (pdgcode && !(particle->getMdstArrayIndex())) {
140 switch (pdgcode) {
141 case 22: //photon conversion
142 rc = false;
143 break ;
144
145 case -11: //bremsstrahlung
146 case 11:
147 rc = true ;
148 break ;
149 default: //everything with boosted flight length less than 1 micrometer
150 rc = (pdgcode && particle->getPDGLifetime() < 1e-14);
151 }
152 }
153 return rc ;
154 }

◆ locate()

const ParticleBase * locate ( Belle2::Particle particle) const
inherited

get particle base from basf2 particle

Definition at line 211 of file ParticleBase.cc.

212 {
213 const ParticleBase* rc = (m_particle == particle) ? this : nullptr;
214 if (!rc) {
215 for (auto* daughter : m_daughters) {
216 rc = daughter->locate(particle);
217 if (rc) {break;}
218 }
219 }
220 return rc;
221 }

◆ momIndex()

int momIndex ( ) const
overridevirtual

momentum index in fit params depending on whether it has a mother

indexing in { x, y, z, tau, px, py, pz, E } but tau is not existing for all InternalParticles

Reimplemented from ParticleBase.

Reimplemented in Resonance.

Definition at line 353 of file InternalParticle.cc.

354 {
359 if (m_geo_constraint && !m_shares_vertex_with_mother) { return this->index() + 4; }
360
361 if (m_shares_vertex_with_mother) { return this->index(); }
362
363 if (!m_geo_constraint) {return index() + 3 ;}
364
365 // this will crash the initialisation
366 return -1;
367 }
int index() const
get index
Definition: ParticleBase.h:95

◆ mother()

const ParticleBase * mother ( ) const
inlineinherited

getMother() / hasMother()

Definition at line 98 of file ParticleBase.h.

98{ return m_mother; };
const ParticleBase * m_mother
motherparticle
Definition: ParticleBase.h:195

◆ nFinalChargedCandidates()

int nFinalChargedCandidates ( ) const
virtualinherited

number of charged candidates

Reimplemented in RecoTrack.

Definition at line 240 of file ParticleBase.cc.

241 {
242 int rc = 0;
243 for (auto* daughter : m_daughters) {
244 rc += daughter->nFinalChargedCandidates();
245 }
246 return rc;
247 }

◆ parname()

std::string parname ( int  index) const
overridevirtual

name

Reimplemented from ParticleBase.

Reimplemented in Resonance.

Definition at line 404 of file InternalParticle.cc.

405 {
406 int id = thisindex;
407 if (!mother() && id >= 3) {++id;}
408 return ParticleBase::parname(id);
409 }
virtual std::string parname(int index) const
get name of parameter i

◆ particle()

Belle2::Particle * particle ( ) const
inlineinherited

get basf2 particle

Definition at line 92 of file ParticleBase.h.

92{ return m_particle ; }

◆ posIndex()

int posIndex ( ) const
overridevirtual

position index in fit params

Reimplemented from ParticleBase.

Reimplemented in Resonance.

Definition at line 330 of file InternalParticle.cc.

331 {
332 // for example B0 and D* can share the same vertex
333 return m_shares_vertex_with_mother ? this->mother()->posIndex() : index();
334 }

◆ projectBeamConstraint()

ErrCode projectBeamConstraint ( const FitParams fitparams,
Projection p 
) const

project beam four momentum constraint

Definition at line 285 of file InternalParticle.cc.

287 {
288
289 const int momindex = momIndex() ;
290
291 const Eigen::Matrix<double, 4, 1> fitMomE = fitparams.getStateVector().segment(momindex, 4);
292
293 p.getResiduals() = m_config->m_beamMomE - fitMomE;
294
295 for (int row = 0; row < 4; ++row) {
296 p.getH()(row, momindex + row) = -1;
297 }
298
299 p.getV() = m_config->m_beamCovariance;
300
301 return ErrCode(ErrCode::Status::success) ;
302 }
Eigen::Matrix< double, 4, 1 > m_beamMomE
Beam four-momentum.
Eigen::Matrix< double, 4, 4 > m_beamCovariance
Beam Covariance.

◆ projectConstraint()

ErrCode projectConstraint ( const Constraint::Type  type,
const FitParams fitparams,
Projection p 
) const
overridevirtual

find out which constraint it is and project

Reimplemented from ParticleBase.

Definition at line 305 of file InternalParticle.cc.

308 {
309 ErrCode status;
310 switch (type) {
311 case Constraint::mass:
312 status |= projectMassConstraint(fitparams, p);
313 break;
314 case Constraint::geometric:
315 status |= projectGeoConstraint(fitparams, p);
316 break;
317 case Constraint::kinematic:
318 status |= projectKineConstraint(fitparams, p);
319 break;
320 case Constraint::beam:
321 status |= projectBeamConstraint(fitparams, p);
322 break;
323 default:
324 status |= ParticleBase::projectConstraint(type, fitparams, p);
325 }
326
327 return status;
328 }
ErrCode projectBeamConstraint(const FitParams &, Projection &) const
project beam four momentum constraint
virtual int type() const override
type
virtual ErrCode projectMassConstraint(const FitParams &, Projection &) const
project mass constraint abstract
virtual ErrCode projectConstraint(Constraint::Type, const FitParams &, Projection &) const
project constraint.
virtual ErrCode projectGeoConstraint(const FitParams &, Projection &) const
project geometrical constraint

◆ projectGeoConstraint()

ErrCode projectGeoConstraint ( const FitParams fitparams,
Projection p 
) const
virtualinherited

project geometrical constraint

the direction of the momentum is very well known from the kinematic constraints that is why we do not extract the distance as a vector here

Definition at line 249 of file ParticleBase.cc.

250 {
251 assert(m_config);
252 // only allow 2d for head of tree particles that are beam constrained
253 const int dim = m_config->m_originDimension == 2 && std::abs(m_particle->getPDGCode()) == m_config->m_headOfTreePDG ? 2 : 3;
254 const int posindexmother = mother()->posIndex();
255 const int posindex = posIndex();
256 const int tauindex = tauIndex();
257 const int momindex = momIndex();
258
259 const double tau = fitparams.getStateVector()(tauindex);
260 Eigen::Matrix < double, 1, -1, 1, 1, 3 > x_vec = fitparams.getStateVector().segment(posindex, dim);
261 Eigen::Matrix < double, 1, -1, 1, 1, 3 > x_m = fitparams.getStateVector().segment(posindexmother, dim);
262 Eigen::Matrix < double, 1, -1, 1, 1, 3 > p_vec = fitparams.getStateVector().segment(momindex, dim);
263 const double mom = p_vec.norm();
264 const double mom3 = mom * mom * mom;
265
266 if (3 == dim) {
267 // we can already set these
268 //diagonal momentum
269 p.getH()(0, momindex) = tau * (p_vec(1) * p_vec(1) + p_vec(2) * p_vec(2)) / mom3 ;
270 p.getH()(1, momindex + 1) = tau * (p_vec(0) * p_vec(0) + p_vec(2) * p_vec(2)) / mom3 ;
271 p.getH()(2, momindex + 2) = tau * (p_vec(0) * p_vec(0) + p_vec(1) * p_vec(1)) / mom3 ;
272
273 //offdiagonal momentum
274 p.getH()(0, momindex + 1) = - tau * p_vec(0) * p_vec(1) / mom3 ;
275 p.getH()(0, momindex + 2) = - tau * p_vec(0) * p_vec(2) / mom3 ;
276
277 p.getH()(1, momindex + 0) = - tau * p_vec(1) * p_vec(0) / mom3 ;
278 p.getH()(1, momindex + 2) = - tau * p_vec(1) * p_vec(2) / mom3 ;
279
280 p.getH()(2, momindex + 0) = - tau * p_vec(2) * p_vec(0) / mom3 ;
281 p.getH()(2, momindex + 1) = - tau * p_vec(2) * p_vec(1) / mom3 ;
282
283 } else if (2 == dim) {
284
285 // NOTE THAT THESE ARE DIFFERENT IN 2d
286 p.getH()(0, momindex) = tau * (p_vec(1) * p_vec(1)) / mom3 ;
287 p.getH()(1, momindex + 1) = tau * (p_vec(0) * p_vec(0)) / mom3 ;
288
289 //offdiagonal momentum
290 p.getH()(0, momindex + 1) = - tau * p_vec(0) * p_vec(1) / mom3 ;
291 p.getH()(1, momindex + 0) = - tau * p_vec(1) * p_vec(0) / mom3 ;
292 } else {
293 B2FATAL("Dimension of Geometric constraint is not 2 or 3. This will crash many things. You should feel bad.");
294 }
295
296 for (int row = 0; row < dim; ++row) {
297
298 double posxmother = x_m(row);
299 double posx = x_vec(row);
300 double momx = p_vec(row);
301
305 p.getResiduals()(row) = posxmother + tau * momx / mom - posx ;
306 p.getH()(row, posindexmother + row) = 1;
307 p.getH()(row, posindex + row) = -1;
308 p.getH()(row, tauindex) = momx / mom;
309 }
310
311 return ErrCode(ErrCode::Status::success);
312 }

◆ projectKineConstraint()

ErrCode projectKineConstraint ( const FitParams fitparams,
Projection p 
) const

project kinematical constraint

Definition at line 237 of file InternalParticle.cc.

239 {
240 const int momindex = momIndex();
241
242 // `this` always has an energy row
243 p.getResiduals().segment(0, 4) = fitparams.getStateVector().segment(momindex, 4);
244
245 for (int imom = 0; imom < 4; ++imom) {
246 p.getH()(imom, momindex + imom) = 1;
247 }
248
249 for (const auto daughter : m_daughters) {
250 const int daumomindex = daughter->momIndex();
251 const Eigen::Matrix<double, 1, 3> p3_vec = fitparams.getStateVector().segment(daumomindex, 3);
252
253 // three momentum is easy just subtract the vectors
254 p.getResiduals().segment(0, 3) -= p3_vec;
255
256 // energy depends on the parametrisation!
257 if (daughter->hasEnergy()) {
258 p.getResiduals()(3) -= fitparams.getStateVector()(daumomindex + 3);
259 p.getH()(3, daumomindex + 3) = -1; // d/dE -E
260 } else {
261 // m^2 + p^2 = E^2
262 // so
263 // E = sqrt(m^2 + p^2)
264 const double mass = daughter->particle()->getPDGMass();
265 const double p2 = p3_vec.squaredNorm();
266 const double energy = std::sqrt(mass * mass + p2);
267 p.getResiduals()(3) -= energy;
268
269 for (unsigned i = 0; i < 3; ++i) {
270 // d/dpx_i sqrt(m^2 + p^2)
271 p.getH()(3, daumomindex + i) = -1 * p3_vec(i) / energy;
272 }
273 }
274
275 // this has to be in any case
276 // d/dp_i p_i
277 for (unsigned i = 0; i < 3; ++i) {
278 p.getH()(i, daumomindex + i) = -1;
279 }
280 }
281 return ErrCode(ErrCode::Status::success);
282 }

◆ projectMassConstraint()

ErrCode projectMassConstraint ( const FitParams fitparams,
Projection p 
) const
virtualinherited

project mass constraint abstract

Definition at line 405 of file ParticleBase.cc.

407 {
408 assert(m_config);
409 if (m_config->m_massConstraintType == 0) {
410 return projectMassConstraintParticle(fitparams, p);
411 } else {
412 return projectMassConstraintDaughters(fitparams, p);
413 }
414 }
const bool m_massConstraintType
const flag for the type of the mass constraint
virtual ErrCode projectMassConstraintParticle(const FitParams &, Projection &) const
project mass constraint using the particles parameters
virtual ErrCode projectMassConstraintDaughters(const FitParams &, Projection &) const
project mass constraint using the parameters of the daughters

◆ projectMassConstraintDaughters()

ErrCode projectMassConstraintDaughters ( const FitParams fitparams,
Projection p 
) const
virtualinherited

project mass constraint using the parameters of the daughters

be aware that the signs here are important E-|p|-m extracts a negative mass and messes with the momentum !

Definition at line 314 of file ParticleBase.cc.

316 {
317 const double mass = particle()->getPDGMass();
318 const double mass2 = mass * mass;
319 double px = 0;
320 double py = 0;
321 double pz = 0;
322 double E = 0;
323
324 // the parameters of the daughters must be used otherwise the mass constraint does not have an effect on the extracted daughter momenta
325 for (const auto* daughter : m_daughters) {
326 const int momindex = daughter->momIndex();
327 // in most cases the daughters will be final states so we cache the value to use it in the energy column
328 const double px_daughter = fitparams.getStateVector()(momindex);
329 const double py_daughter = fitparams.getStateVector()(momindex + 1);
330 const double pz_daughter = fitparams.getStateVector()(momindex + 2);
331
332 px += px_daughter;
333 py += py_daughter;
334 pz += pz_daughter;
335 if (daughter->hasEnergy()) {
336 E += fitparams.getStateVector()(momindex + 3);
337 } else {
338 // final states dont have an energy index
339 const double m = daughter->particle()->getPDGMass();
340 E += std::sqrt(m * m + px_daughter * px_daughter + py_daughter * py_daughter + pz_daughter * pz_daughter);
341 }
342 }
343
347 p.getResiduals()(0) = mass2 - E * E + px * px + py * py + pz * pz;
348
349 for (const auto* daughter : m_daughters) {
350 //dr/dx = d/dx m2-{E1+E2+...}^2+{p1+p2+...}^2 = 2*x (x= E or p)
351 const int momindex = daughter->momIndex();
352 p.getH()(0, momindex) = 2.0 * px;
353 p.getH()(0, momindex + 1) = 2.0 * py;
354 p.getH()(0, momindex + 2) = 2.0 * pz;
355
356 if (daughter->hasEnergy()) {
357 p.getH()(0, momindex + 3) = -2.0 * E;
358 } else {
359 const double px_daughter = fitparams.getStateVector()(momindex);
360 const double py_daughter = fitparams.getStateVector()(momindex + 1);
361 const double pz_daughter = fitparams.getStateVector()(momindex + 2);
362 const double m = daughter->particle()->getPDGMass();
363
364 const double E_daughter = std::sqrt(m * m + px_daughter * px_daughter + py_daughter * py_daughter + pz_daughter * pz_daughter);
365 const double E_by_E_daughter = E / E_daughter;
366 p.getH()(0, momindex) -= 2.0 * E_by_E_daughter * px_daughter;
367 p.getH()(0, momindex + 1) -= 2.0 * E_by_E_daughter * py_daughter;
368 p.getH()(0, momindex + 2) -= 2.0 * E_by_E_daughter * pz_daughter;
369 }
370
371 }
372 return ErrCode(ErrCode::Status::success);
373 }

◆ projectMassConstraintParticle()

ErrCode projectMassConstraintParticle ( const FitParams fitparams,
Projection p 
) const
virtualinherited

project mass constraint using the particles parameters

be aware that the signs here are important E-|p|-m extracts a negative mass and messes with the momentum !

Definition at line 375 of file ParticleBase.cc.

377 {
378 const double mass = particle()->getPDGMass();
379 const double mass2 = mass * mass;
380 const int momindex = momIndex();
381 const double px = fitparams.getStateVector()(momindex);
382 const double py = fitparams.getStateVector()(momindex + 1);
383 const double pz = fitparams.getStateVector()(momindex + 2);
384 const double E = fitparams.getStateVector()(momindex + 3);
385
389 p.getResiduals()(0) = mass2 - E * E + px * px + py * py + pz * pz;
390
391 p.getH()(0, momindex) = 2.0 * px;
392 p.getH()(0, momindex + 1) = 2.0 * py;
393 p.getH()(0, momindex + 2) = 2.0 * pz;
394 p.getH()(0, momindex + 3) = -2.0 * E;
395
396 // TODO 0 in most cases -> needs special treatment if width=0 to not crash chi2 calculation
397 // const double width = TDatabasePDG::Instance()->GetParticle(particle()->getPDGCode())->Width();
398 // transport measurement uncertainty into residual system
399 // f' = sigma_x^2 * (df/dx)^2
400 // p.getV()(0) = width * width * 4 * mass2;
401
402 return ErrCode(ErrCode::Status::success);
403 }

◆ removeDaughter()

void removeDaughter ( const ParticleBase pb)
virtualinherited

remove daughter

Definition at line 73 of file ParticleBase.cc.

74 {
75 auto iter = std::find(m_daughters.begin(), m_daughters.end(), pb);
76 if (iter != m_daughters.end()) {
77 delete *iter;
78 m_daughters.erase(iter);
79 } else {
80 B2ERROR("Cannot remove particle, because not found ...");
81 }
82 }

◆ retrieveIndexMap()

void retrieveIndexMap ( indexmap anindexmap) const
virtualinherited

get index map

Definition at line 223 of file ParticleBase.cc.

224 {
225 map.push_back(std::pair<const ParticleBase*, int>(this, index()));
226 for (auto* daughter : m_daughters) {
227 daughter->retrieveIndexMap(map);
228 }
229 }

◆ setIndex()

void setIndex ( int  i)
inlineprotectedinherited

set Index (in statevector)

Definition at line 189 of file ParticleBase.h.

189{ m_index = i ; }

◆ setMassConstraint()

void setMassConstraint ( bool  b)
inline

set mass constraint flag

Definition at line 80 of file InternalParticle.h.

80{ m_massconstraint = b ; }

◆ setMother()

void setMother ( const ParticleBase m)
inlineinherited

set mother

Definition at line 164 of file ParticleBase.h.

164{ m_mother = m ; }

◆ tauIndex()

int tauIndex ( ) const
overridevirtual

tau index in fit params only if it has a mother

only exists if particle is geo constraint and has a mother

Reimplemented from ParticleBase.

Reimplemented in Resonance.

Definition at line 347 of file InternalParticle.cc.

348 {
350 return m_geo_constraint ? index() + 3 : -1;
351 }

◆ type()

virtual int type ( ) const
inlineoverridevirtual

type

Implements ParticleBase.

Reimplemented in Resonance.

Definition at line 56 of file InternalParticle.h.

56{ return kInternalParticle ; }

◆ updateIndex()

void updateIndex ( int &  offset)
virtualinherited

this sets the index for momentum, position, etc.

in the statevector

Definition at line 84 of file ParticleBase.cc.

85 {
86 for (auto* daughter : m_daughters) {
87 daughter->updateIndex(offset);
88 }
89 m_index = offset;
90 offset += dim();
91 }

Member Data Documentation

◆ m_automatic_vertex_constraining

bool m_automatic_vertex_constraining
private

automatically figure out if mother and particle vertex should be the same and also add geometric constraints

Definition at line 112 of file InternalParticle.h.

◆ m_beamconstraint

bool m_beamconstraint
private

has beam constraint

Definition at line 96 of file InternalParticle.h.

◆ m_config

const ConstraintConfiguration* m_config
protectedinherited

has all the constraint config

Definition at line 204 of file ParticleBase.h.

◆ m_daughters

std::vector<ParticleBase*> m_daughters
protectedinherited

daughter container

Definition at line 198 of file ParticleBase.h.

◆ m_geo_constraint

bool m_geo_constraint
private

use a geo metric constraint

Definition at line 102 of file InternalParticle.h.

◆ m_index

int m_index
privateinherited

index

Definition at line 208 of file ParticleBase.h.

◆ m_isconversion

bool m_isconversion
private

is conversion

Definition at line 108 of file InternalParticle.h.

◆ m_isStronglyDecayingResonance

bool m_isStronglyDecayingResonance
protectedinherited

decay length less than 1 micron

Definition at line 201 of file ParticleBase.h.

◆ m_lifetimeconstraint

bool m_lifetimeconstraint
private

has lifetime constraint

Definition at line 105 of file InternalParticle.h.

◆ m_massconstraint

bool m_massconstraint
private

has mass constraint

Definition at line 93 of file InternalParticle.h.

◆ m_mother

const ParticleBase* m_mother
protectedinherited

motherparticle

Definition at line 195 of file ParticleBase.h.

◆ m_name

std::string m_name
privateinherited

name

Definition at line 211 of file ParticleBase.h.

◆ m_particle

Belle2::Particle* m_particle
protectedinherited

pointer to framework type

Definition at line 192 of file ParticleBase.h.

◆ m_shares_vertex_with_mother

bool m_shares_vertex_with_mother
private

shares vertex with mother, that means decay vertex = productionvertex

Definition at line 99 of file InternalParticle.h.


The documentation for this class was generated from the following files: