TimeExtractionUtils.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost, Nils Braun                                 *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#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(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);
  };
}