RecoNeutral.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 * External Contributor: Wouter Hulsbergen                                *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <analysis/VertexFitting/TreeFitter/RecoParticle.h>
 
#include <Eigen/Core>
 
namespace TreeFitter {

  class RecoNeutral : public RecoParticle {
 
  public:
    RecoNeutral(Belle2::Particle* bc, const ParticleBase* mother) ;

    virtual ~RecoNeutral() {};

    virtual ErrCode initParticleWithMother(FitParams& fitparams) override;

    virtual ErrCode initMotherlessParticle(FitParams& fitparams) override;

    ErrCode initCovariance(FitParams& fitparams) const override;

    ErrCode initParams();

    ErrCode projectRecoConstraint(const FitParams& fitparams, Projection& p) const override;

    virtual int dimM() const override { return dim(); }

    virtual bool hasEnergy() const override { return false; }

    virtual int dim() const override { return m_dim; }

    virtual int type() const override { return kRecoNeutral ; }

    virtual void addToConstraintList(constraintlist& alist, int depth) const override
    {
      alist.push_back(Constraint(this, Constraint::neutralHadron, depth, dimM())) ;
    }
 
  private:

    const int m_dim;

    bool m_init ;

    Eigen::Matrix<double, 1, 4> m_clusterPars;

    Eigen::Matrix<double, 4, 4> m_covariance;

    int m_i1;
    int m_i2;
    int m_i3;

    const float m_momentumScalingFactor;

    const double m_mass;

    const Belle2::Particle::EParticleSourceObject m_particleSource;
  };
 
}