InvariantMassStandAlone.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2021 - 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 <utility>
#include <tuple>
#include <TVector3.h>
#include <Eigen/Dense>
 
//If compiled within BASF2
#ifdef _PACKAGE_
#include <tracking/calibration/Splitter.h>
#else
#include <Splitter.h>
#endif
 
namespace Belle2 {
  namespace InvariantMassCalib {
 
    static const double realNaN = std::numeric_limits<double>::quiet_NaN();
    static const int intNaN     = std::numeric_limits<int>::quiet_NaN();
    static const TVector3 vecNaN(realNaN, realNaN, realNaN);
 
    struct Track {
      TVector3 p = vecNaN; 
      double pid = realNaN; 
    };
 
    struct Event {
      int exp   = intNaN;  
      int run   = intNaN;  
      int evtNo = intNaN;  
 
 
      Track mu0; 
      Track mu1; 
 
      int nBootStrap = intNaN; 
      bool isSig = false;      
      double t   = realNaN;    
 
    };
 
 
    std::vector<Event> getEvents(TTree* tr);
 
    std::tuple<std::vector<Eigen::VectorXd>, std::vector<Eigen::MatrixXd>, Eigen::MatrixXd>  runInvariantMassAnalysis(
      std::vector<Event> evts,
      const std::vector<double>& splitPoints);
 
  }
}