DecayChain.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/ParticleBase.h>
#include <analysis/VertexFitting/TreeFitter/MergedConstraint.h>
#include <analysis/VertexFitting/TreeFitter/ConstraintConfiguration.h>
 
namespace TreeFitter {
 
  class FitParams ;
  class DecayChain {
 
  public:

    DecayChain() : m_dim(0), m_headOfChain(0), m_cand(0), m_isOwner(true) {} //Default constructor (actually never used)

    DecayChain(Belle2::Particle* bc,
               const ConstraintConfiguration& config,
               bool forceFitAll = false
              );

    ~DecayChain();

    ErrCode initialize(FitParams& par);

    ErrCode filter(FitParams& par);

    ErrCode filterWithReference(FitParams& par, const FitParams& ref);

    int dim() const { return m_dim;}

    void initConstraintList();

    void removeConstraintFromList();

    double chiSquare(const FitParams& par) const;

    ParticleBase* mother() { return m_headOfChain ; }

    const ParticleBase* cand() { return m_cand ; }

    const ParticleBase* mother() const { return m_headOfChain ; }

    const ParticleBase* locate(Belle2::Particle* bc) const ;

    int index(Belle2::Particle* bc) const ;

    int posIndex(Belle2::Particle* bc) const ;

    int momIndex(Belle2::Particle* bc) const ;

    int tauIndex(Belle2::Particle* bc) const ;

    int momIndex() const ;
 
  private:
 
    mutable int m_dim ; 
 
    mutable ParticleBase* m_headOfChain ; 
 
    const ParticleBase* m_cand ; 

    ParticleBase::constraintlist m_constraintlist ;

    typedef std::map<Belle2::Particle*, const ParticleBase*> ParticleMap ;

    ParticleMap m_particleMap ;

    const bool m_isOwner ;

    const ConstraintConfiguration m_config;
 
  };
 
}