RecoComposite.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2018 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributor: Wouter Hulsbergen, Francesco Tenchini, Jo-Frederik Krohn  *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <analysis/VertexFitting/TreeFitter/ParticleBase.h>
 
namespace TreeFitter {
 
  class RecoComposite : public ParticleBase {
 
  public:
    RecoComposite(Belle2::Particle* bc, const ParticleBase* mother) ;
 
    virtual ~RecoComposite() {};
 
    virtual ErrCode initParticleWithMother(FitParams& fitparams) override;
 
    virtual ErrCode initMotherlessParticle(FitParams& fitparams) override;
 
    void updateParams();
 
    ErrCode projectRecoComposite(const FitParams& fitparams, Projection& p) const;
 
    int dimMeas() const        { return m_hasEnergy ? 7 : 6 ; }
 
    virtual int dim() const override { return m_hasEnergy ? 8 : 7 ; }// (x,y,z,t,px,py,pz,(E))
 
    int dimM() const        { return m_hasEnergy ? 7 : 6 ; }
 
    virtual ErrCode projectConstraint(Constraint::Type, const FitParams&, Projection&) const override ;
 
    virtual int type() const override { return kRecoComposite ; }
 
    virtual int posIndex() const override { return index()   ; }
 
    virtual int tauIndex() const override { return index() + 3 ; }
 
    virtual int momIndex() const override { return index() + 4 ; }
 
    virtual bool hasEnergy() const override { return m_hasEnergy ; }
 
    virtual bool hasPosition() const override { return true ; }
 
    virtual void addToConstraintList(constraintlist& alist, int depth) const override
    {
      alist.push_back(Constraint(this, Constraint::composite, depth, dimM())) ;
      alist.push_back(Constraint(this, Constraint::geometric, depth, 3)) ;
    }
 
  protected:
 
    Eigen::Matrix<double, 7, 1> m_params;
 
    Eigen::Matrix < double, -1, -1, 0, 7, 7 >  m_covariance;
 
    bool m_hasEnergy;
  };
}