SimpEquations.h

Go to the documentation of this file.

#include <iostream>
#include <TEfficiency.h>
#include <TMath.h>
#include <TH1F.h>
#include "json.hpp"
 
using json = nlohmann::json;
 
class SimpEquations {
 
    public:
            
 
        SimpEquations(int year);
 
        SimpEquations(int year, const std::string paramsConfigFile);
 
        void loadParametersConfig(const std::string paramsConfigFile);
 
        double rate_2pi(double m_Ap, double m_pi, double m_V, double alpha_dark);
 
        double rate_Vrho_pi(double m_Ap, double m_pi, double m_V,
                        double alpha_dark, double f_pi);
        double rate_Vphi_pi(double m_Ap, double m_pi, double m_V,
                        double alpha_dark, double f_pi);
        double rate_Vcharged_pi(double m_Ap, double m_pi, double m_V,
                        double alpha_dark, double f_pi);
        double br_2pi(double m_Ap, double m_pi, double m_V, 
                double alpha_dark, double f_pi);
        double br_Vrho_pi(double m_Ap, double m_pi, double m_V, double alpha_dark,
                double f_pi);
 
        double br_Vphi_pi(double m_Ap, double m_pi, double m_V, double alpha_dark,
                double f_pi);
        double br_Vcharged_pi(double m_Ap, double m_pi, double m_V, 
                double alpha_dark, double f_pi);
 
        double br_2V(double m_Ap,double m_pi,double m_V,double alpha_dark,double f_pi,double rho,double phi);
        
        double Tv(bool rho,bool phi);
        
        double Beta(double x,double y);
        
        double rate_2V(double m_Ap,double m_V,double alpha_dark);
        
        double f(double r);
        
        double rate_2l(double m_Ap,double m_pi,double m_V,double eps,double alpha_dark,double f_pi,double m_l,bool rho);
        
        double getCtau(double m_Ap,double m_pi,double m_V,double eps,double alpha_dark,double f_pi,double m_l,bool rho);
        
        double gamma(double m_V,double E_V);
 
        double radiativeFraction(double m_Ap);
 
        double radiativeAcceptance(double m_Ap);
 
        double massResolution(double m_V);
 
        double controlRegionBackgroundRate(double m_Ap);
 
        double expectedSignalCalculation(double m_V, double eps, bool rho,
                        double E_V, TEfficiency* effCalc_h, double target_pos, double zcut);
 
        double expectedSignalCalculation(double m_Ap, double m_pi, double m_V, double eps, double alpha_dark,
            double f_pi, double m_l, bool rho,double E_V, TEfficiency* effCalc_h, double target_pos, 
            double zcut);
 
        double expectedSignalCalculation(double m_V, double eps, bool rho, double E_V, 
                TEfficiency* effCalc_h, double dNdm, double radFrac, double radAcc, double target_pos, double zcut);
        
        double getAprimeMassFromVectorMass(double m_V){return m_V * mass_ratio_Ap_to_Vd_;};
 
    private:
 
        int year_ = 2016;
        json params_config_;
        double mass_ratio_Ap_to_Pid_ = 3.0;
        double mass_ratio_Ap_to_Vd_ = 3.0/1.8;
        double ratio_mPi_to_fPi_ = 4.*M_PI;
        double m_l_ = 0.511;
        double alpha_dark_ = 0.1;
 
};
class SimpEquations {…};