DecayChain.h

/**************************************************************************
 *
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2013 - 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>
#include <analysis/VertexFitting/TreeFitter/MergedConstraint.h>
#include <analysis/VertexFitting/TreeFitter/ConstraintConfiguration.h>
 
namespace TreeFitter {
 
  class FitParams ;
  class ParticleBase ;
  class ConstraintConfiguration;
  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;
 
  };
 
}