BoostVectorStandAlone.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/geometry/B2Vector3.h>
#include <vector>
#include <tuple>
#include <Eigen/Dense>
 
//If compiled within BASF2
#ifdef _PACKAGE_
#include <reconstruction/calibration/BeamSpotBoostInvMass/Splitter.h>
#else
#include <Splitter.h>
#endif
 
namespace Belle2::BoostVectorCalib {
 
  static const double realNaN = std::numeric_limits<double>::quiet_NaN();
  static const B2Vector3D vecNaN(realNaN, realNaN, realNaN);
 
  struct Track {
    B2Vector3D p = vecNaN; 
    double pid = realNaN; 
  };
 
  struct Event {
    int exp   = -1; 
    int run   = -1; 
    int evtNo = -1; 
 
 
    Track mu0; 
    Track mu1; 
 
    int nBootStrap = 1;     
    bool isSig = false;     
    double t = realNaN;     
 
  };
 
 
  std::vector<Event> getEvents(TTree* tr);
 
 
  double median(double* v, int n);
 
 
  std::tuple<std::vector<Eigen::VectorXd>, std::vector<Eigen::MatrixXd>, Eigen::MatrixXd>  runBoostVectorAnalysis(
    std::vector<Event> evts,
    const std::vector<double>& splitPoints);
 
}