TimeExtractionUtils.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 <framework/datastore/StoreObjPtr.h>
#include <framework/dataobjects/EventT0.h>
 
#include <vector>
#include <TVectorD.h>
#include <TMatrixDSym.h>
 
namespace Belle2 {
  class RecoTrack;

  class TimeExtractionUtils {
  public:
    static std::pair<double, double> getChi2WithFit(const std::vector<RecoTrack*>& recoTracks, bool setDirtyFlag);

    static std::pair<double, double> getExtractedTimeAndUncertaintyWithFit(const std::vector<RecoTrack*>& recoTracks,
        bool setDirtyFlag);

    static void addEventT0WithQuality(const std::vector<RecoTrack*>& recoTracks, StoreObjPtr<EventT0>& eventT0,
                                      std::vector<EventT0::EventT0Component>& eventT0WithQualityIndex);

    static double extractReducedChi2(const RecoTrack& recoTrack);

    static std::pair<double, double> getChi2Derivatives(const RecoTrack& recoTrack);
 
  private:
    static std::vector<int> getMeasurementDimensions(const RecoTrack& recoTrack);

    static bool buildFullCovarianceMatrix(const RecoTrack& recoTrack,
                                          TMatrixDSym& fullCovariance);

    static bool buildFullResidualCovarianceMatrix(const RecoTrack& recoTrack,
                                                  const std::vector<int>& vDimMeas,
                                                  const TMatrixDSym& fullCovariance,
                                                  TMatrixDSym& fullResidualCovariance,
                                                  TMatrixDSym& inverseFullMeasurementCovariance);

    static void buildResidualsAndTimeDerivative(const RecoTrack& recoTrack,
                                                const std::vector<int>& vDimMeas,
                                                TVectorD& residuals,
                                                TVectorD& residualTimeDerivative);
  };

}