QualityEstimatorMC.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Jonas Wagner                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <tracking/trackFindingVXD/trackQualityEstimators/QualityEstimatorBase.h>
#include <framework/datastore/StoreArray.h>
#include <tracking/dataobjects/RecoTrack.h>
#include <tracking/spacePointCreation/SpacePoint.h>
 
namespace Belle2 {
  class QualityEstimatorMC : public QualityEstimatorBase {
  public:
    // some typedefs to increase readability
    typedef int MCRecoTrackIndex; 
    typedef unsigned int NMatches; 
    typedef std::pair<MCRecoTrackIndex, NMatches> MatchInfo; 
    QualityEstimatorMC(const std::string& mcRecoTracksStoreArrayName = "MCRecoTracks",
                       bool strictQualityIndicator = true, bool mva_target = false):
      QualityEstimatorBase(), m_strictQualityIndicator(strictQualityIndicator), m_mva_target(mva_target),
      m_mcRecoTracksStoreArrayName(mcRecoTracksStoreArrayName)
    {
      m_mcRecoTracks.isRequired(m_mcRecoTracksStoreArrayName);
    };
 
    virtual double estimateQuality(std::vector<SpacePoint const*> const& measurements) final;
 
    virtual QualityEstimationResults estimateQualityAndProperties(std::vector<SpacePoint const*> const& measurements) override final;
 
    void setClustersNames(const std::string& svdClustersName, const std::string& pxdClustersName)
    { m_svdClustersName = svdClustersName; m_pxdClustersName = pxdClustersName; };
 
  protected:
    MatchInfo getBestMatchToMCClusters(std::vector<SpacePoint const*> const& measurements);
 
    double calculateQualityIndicator(unsigned int nClusters, MatchInfo& match);
 
    // parameters
    bool m_strictQualityIndicator;
 
    bool m_mva_target;
 
    // module members
    std::string m_mcRecoTracksStoreArrayName; 
    std::string m_svdClustersName = ""; 
    std::string m_pxdClustersName = ""; 
    MatchInfo m_match;
 
    StoreArray<RecoTrack> m_mcRecoTracks;
  };
}