BeamSpotStandAlone.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2020 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Radek Zlebcik                                            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <vector>
#include <tuple>
#include <TTree.h>
#include <Eigen/Dense>
 
//If compiled within BASF2
#ifdef _PACKAGE_
#include <tracking/calibration/Splitter.h>
#else
#include <Splitter.h>
#endif
 
namespace Belle2 {
  namespace BeamSpotCalib {
 
 
    struct Track {
      double d0;   
      double z0;   
      double phi0; 
      double tanlambda; 
    };
 
    struct Event {
      int exp;   
      int run;   
      int evtNo; 
 
 
      Track mu0; 
      Track mu1; 
 
      int nBootStrap; 
      bool isSig;     
      double t;       
 
      void toMicroM()
      {
        mu0.d0 *= 1e4;
        mu1.d0 *= 1e4;
        mu0.z0 *= 1e4;
        mu1.z0 *= 1e4;
      }
    };
 
 
 
    std::vector<Event> getEvents(TTree* tr);
 
    std::tuple<std::vector<Eigen::VectorXd>, std::vector<Eigen::MatrixXd>, Eigen::MatrixXd>  runBeamSpotAnalysis(
      std::vector<Event> evts,
      const std::vector<double>& splitPoints);
 
 
  }
}