representation of the Klong constraint More...

#include <RecoKlong.h>

Inheritance diagram for RecoKlong:

Collaboration diagram for RecoKlong:

[legend]

Public Types
enum	TFParticleType { kInteractionPoint , kOrigin , kComposite , kRecoResonance , kInternalParticle , kRecoTrack , kResonance , kRecoPhoton , kRecoKlong , kMissingParticle }
	particle types

typedef std::vector< Constraint >	constraintlist
	alias

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

Public Member Functions
	RecoKlong (Belle2::Particle bc, const ParticleBase mother)
	constructor

virtual	~RecoKlong ()
	destructor

virtual ErrCode	initParticleWithMother (FitParams &fitparams) override
	init particle with mother

virtual ErrCode	initMotherlessParticle (FitParams &fitparams) override
	init particle without mother

ErrCode	initCovariance (FitParams &fitparams) const override
	init covariance

ErrCode	initParams ()
	update or init params More...

ErrCode	projectRecoConstraint (const FitParams &fitparams, Projection &p) const override
	project klong constraint

virtual int	dimM () const override
	sets the size of the corresponding residual projection

virtual bool	hasEnergy () const override
	how should the energy be calculated ? from momentum or from E ?

virtual int	dim () const override
	set the size of the particle in the statevector

virtual int	type () const override
	type

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

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

virtual int	momIndex () const override
	get momentum index

virtual ErrCode	projectConstraint (Constraint::Type, const FitParams &, Projection &) const override
	abstract abstract projection

virtual void	updateIndex (int &offset)
	this sets the index for momentum, position, etc. More...

const ParticleBase *	locate (Belle2::Particle *particle) const
	get particle base from basf2 particle

Belle2::Particle *	particle () const
	get basf2 particle

int	index () const
	get index

const ParticleBase *	mother () const
	getMother() / hasMother()

virtual ErrCode	projectGeoConstraint (const FitParams &, Projection &) const
	project geometrical constraint More...

virtual ErrCode	projectMassConstraintParticle (const FitParams &, Projection &) const
	project mass constraint using the particles parameters More...

virtual ErrCode	projectMassConstraintDaughters (const FitParams &, Projection &) const
	project mass constraint using the parameters of the daughters More...

virtual ErrCode	projectMassConstraint (const FitParams &, Projection &) const
	project mass constraint abstract

virtual void	forceP4Sum (FitParams &) const
	force p4 sum conservation all along the tree

virtual int	posIndex () const
	get vertex index (in statevector!)

virtual int	tauIndex () const
	get tau index

virtual bool	hasPosition () const
	get false

int	eneIndex () const
	get energy index

virtual double	chiSquare (const FitParams &) const
	get chi2

int	charge () const
	get charge

virtual ParticleBase *	addDaughter (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 bool	useEnergy (Belle2::Particle &cand)
	has energy in fit params?

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

static ParticleBase *	createOrigin (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	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. More...

Protected Attributes
Belle2::Particle *	m_particle
	pointer to framework type

const ParticleBase *	m_mother
	motherparticle

std::vector< ParticleBase * >	m_daughters
	daughter container

bool	m_isStronglyDecayingResonance
	decay length less than 1 micron

const ConstraintConfiguration *	m_config
	has all the constraint config

Private Attributes
const int	m_dim
	dimension of residuals and 'width' of H

bool	m_init
	was initialized*

bool	m_useEnergy
	has energy ins statevector

Eigen::Matrix< double, 1, 4 >	m_clusterPars
	constrains measured params (x_c, y_c, z_c, E_c)

Eigen::Matrix< double, 4, 4 >	m_covariance
	covariance (x_c,y_c,z_c,E_c) of measured pars

int	m_i1
	index with the highest momentum. More...

int	m_i2
	random index

int	m_i3
	another random index

int	m_index
	index

std::string	m_name
	name

Detailed Description

representation of the Klong constraint

Definition at line 16 of file RecoKlong.h.

Member Function Documentation

◆ initParams()

ErrCode initParams ( )

update or init params

currently the energy in KLM is calculated as n2dHits in cluster times 0.214 GeV at time of writing - 8.3.18 - the KLMCluster returns 0 for the E covariance

Definition at line 79 of file RecoKlong.cc.

   {
     const Belle2::KLMCluster* cluster = particle()->getKLMCluster();
  
     const Belle2::B2Vector3D cluster_pos = cluster->getClusterPosition();
  
     m_init = true;
     m_covariance =  Eigen::Matrix<double, 4, 4>::Zero(4, 4);
  
     TMatrixDSym cov7 = cluster->getError7x7();
  
     for (int row = 0; row < 3; ++row) {
       for (int col = 0; col < 3; ++col) {
         m_covariance(row, col) = cov7[row + 4][col + 4] ;
       }
     }
  
     if (0 == m_covariance(3, 3)) {m_covariance(3, 3) = .214;}
  
     m_clusterPars(0) = cluster_pos.X();
     m_clusterPars(1) = cluster_pos.Y();
     m_clusterPars(2) = cluster_pos.Z();
     m_clusterPars(3) = sqrt(particle()->getPDGMass() * particle()->getPDGMass() +
                             cluster->getMomentumMag() * cluster->getMomentumMag());
  
     auto p_vec = particle()->getMomentum();
     // find highest momentum, eliminate dim with highest mom
     if ((std::abs(p_vec.X()) >= std::abs(p_vec.Y())) && (std::abs(p_vec.X()) >= std::abs(p_vec.Z()))) {
       m_i1 = 0; m_i2 = 1; m_i3 = 2;
     } else if ((std::abs(p_vec.Y()) >= std::abs(p_vec.X())) && (std::abs(p_vec.Y()) >= std::abs(p_vec.Z()))) {
       m_i1 = 1; m_i2 = 0; m_i3 = 2;
     } else if ((std::abs(p_vec.Z()) >= std::abs(p_vec.Y())) && (std::abs(p_vec.Z()) >= std::abs(p_vec.X()))) {
       m_i1 = 2; m_i2 = 1; m_i3 = 0;
     } else {
       B2ERROR("Could not estimate highest momentum for Klong constraint. Aborting this fit.\n px: "
               << p_vec.X() << " py: " << p_vec.Y() << " pz: " << p_vec.Z() << " calculated from Ec: " << m_clusterPars(3));
       return ErrCode(ErrCode::Status::photondimerror);
     }
  
     return ErrCode(ErrCode::Status::success);
   }

◆ 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.

◆ 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.

◆ 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.

◆ 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.

◆ 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.

Member Data Documentation

◆ m_i1

int m_i1

private

index with the highest momentum.

We have to make sure this does not change during the fit.

Definition at line 79 of file RecoKlong.h.

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

analysis/VertexFitting/TreeFitter/include/RecoKlong.h
analysis/VertexFitting/TreeFitter/src/RecoKlong.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Function Documentation

◆ initParams()

◆ isAResonance()

◆ projectGeoConstraint()

◆ projectMassConstraintDaughters()

◆ projectMassConstraintParticle()

◆ updateIndex()

Member Data Documentation

◆ m_i1