InvariantMassMuMuStandAlone.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 <vector>
#include <tuple>
#include <functional>
#include <framework/geometry/B2Vector3.h>
#include <TTree.h>
#include <Eigen/Dense>
 
//If compiled within BASF2
#ifdef _PACKAGE_
#include <reconstruction/calibration/BeamSpotBoostInvMass/Splitter.h>
#else
#include <Splitter.h>
#endif
 
namespace Belle2::InvariantMassMuMuCalib {
 
  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; 
    double m; 
 
    int nBootStrap = 1;      
    bool isSig = false;      
    double t   = realNaN;    
    bool is4S = false;       
 
  };
 
 
  std::vector<Event> getEvents(TTree* tr, bool is4S);
 
  std::tuple<std::vector<Eigen::VectorXd>, std::vector<Eigen::MatrixXd>, Eigen::MatrixXd>  runMuMuInvariantMassAnalysis(
    std::vector<Event> evts,
    const std::vector<double>& splitPoints);
 
 
  double gausExpConv(double mean, double sigma, double bMean, double bDelta, double tauL, double tauR, double sigmaK, double x);
 
 
  double integrate(std::function<double(double)> f, double a, double b);
 
}