InvariantMassBhadStandAlone.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::InvariantMassBhadCalib {
 
  static const double realNaN = std::numeric_limits<double>::quiet_NaN();
  static const B2Vector3D vecNaN(realNaN, realNaN, realNaN);
 
 
  struct Candidate {
    double pBcms = realNaN;   
    double mB    = realNaN;   
    int pdg = 0;              
    int mode = -1;            
    double Kpid = realNaN;    
    double R2 = realNaN;      
    double mD = realNaN;      
    double dmDstar = realNaN; 
 
    bool isSig = false;       
  };
 
 
 
  struct Event {
 
    int exp   = -1;       
    int run   = -1;       
    int evtNo = -1;       
    double t = realNaN;   
 
 
    std::vector<Candidate>  cand = {}; 
 
 
    int nBootStrap = 1;  
  };
 
 
 
 
 
  std::vector<Event> getEvents(TTree* tr);
 
 
  std::vector<std::vector<double>> doBhadFit(const std::vector<Event>& evts, std::vector<std::pair<double, double>> limits,
                                             std::vector<std::pair<double, double>> mumuVals,  const std::vector<double>& startPars);
 
  std::vector<double> doBhadOnlyFit(const std::vector<Event>& evts, const std::vector<std::pair<double, double>>& limits);
 
 
  std::tuple<std::vector<Eigen::VectorXd>, std::vector<Eigen::MatrixXd>, Eigen::MatrixXd>  runInvariantMassAnalysis(
    std::vector<Event> evts,
    const std::vector<double>& splitPoints);
 
}