V0.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
#include <mdst/dataobjects/TrackFitResult.h>
#include <mdst/dataobjects/Track.h>
#include <framework/datastore/RelationsObject.h>
#include <framework/datastore/StoreArray.h>
#include <framework/gearbox/Const.h>
#include <utility>
 
namespace Belle2 {
  class V0 : public RelationsObject {
  public:
    V0() {};

    V0(const std::pair<const Belle2::Track*, const Belle2::TrackFitResult*>& trackPairPositive,
       const std::pair<const Belle2::Track*, const Belle2::TrackFitResult*>& trackPairNegative,
       double vertexX, double vertexY, double vertexZ);

    std::pair<Track*, Track*> getTracks() const
    {
      StoreArray<Track> tracks;
      return std::make_pair(tracks[m_trackIndexPositive], tracks[m_trackIndexNegative]);
    }

    std::pair<short, short> getTrackIndices() const
    {
      return std::make_pair(m_trackIndexPositive, m_trackIndexNegative);
    }

    std::pair<TrackFitResult*, TrackFitResult*> getTrackFitResults() const
    {
      StoreArray<TrackFitResult> trackFitResults;
      return std::make_pair(trackFitResults[m_trackFitResultIndexPositive], trackFitResults[m_trackFitResultIndexNegative]);
    }

    std::pair<short, short> getTrackFitResultIndices() const
    {
      return std::make_pair(m_trackFitResultIndexPositive, m_trackFitResultIndexNegative);
    }

    Const::ParticleType getV0Hypothesis() const;

    double getFittedVertexX() const { return m_fittedVertexX; }

    double getFittedVertexY() const { return m_fittedVertexY; }

    double getFittedVertexZ() const { return m_fittedVertexZ; }

    ROOT::Math::XYZVector getFittedVertexPosition() const { return ROOT::Math::XYZVector(m_fittedVertexX, m_fittedVertexY, m_fittedVertexZ); }
 
  private:
    Double32_t m_fittedVertexX = 0.0;

    Double32_t m_fittedVertexY = 0.0;

    Double32_t m_fittedVertexZ = 0.0;

    short m_trackIndexPositive = -1;

    short m_trackIndexNegative = -1;

    short m_trackFitResultIndexPositive = -1;

    short m_trackFitResultIndexNegative = -1;

    ClassDef(V0, 5);
 
    friend class FixMergedObjectsModule;
  };

}