From 415ca089d5a02cef840b6cfa9f46391e5ac41e39 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Renata=20Kopecn=C3=A1?= <kopecna@physi.uni-heidelberg.de>
Date: Thu, 27 Jan 2022 19:05:21 +0100
Subject: [PATCH] Update 'GlobalFunctions.hh'

---
 GlobalFunctions.cpp.md |   1 -
 GlobalFunctions.hh.md  | 372 +++++++++++++++++++++++++++++++++++++++++
 2 files changed, 372 insertions(+), 1 deletion(-)
 delete mode 100644 GlobalFunctions.cpp.md
 create mode 100644 GlobalFunctions.hh.md

diff --git a/GlobalFunctions.cpp.md b/GlobalFunctions.cpp.md
deleted file mode 100644
index 5d08b7b..0000000
--- a/GlobalFunctions.cpp.md
+++ /dev/null
@@ -1 +0,0 @@
-Welcome to the Wiki.
\ No newline at end of file
diff --git a/GlobalFunctions.hh.md b/GlobalFunctions.hh.md
new file mode 100644
index 0000000..01c30df
--- /dev/null
+++ b/GlobalFunctions.hh.md
@@ -0,0 +1,372 @@
+
+
+## Global variables:
+
+* const bool SplitDDandLL = KshortChannel ? true : false //for pi0 SplitDDandLL is always false!
+* const bool Kst2Kpluspi0Resolved = !KshortChannel //keep the old variables
+* but set global!
+* const bool Kst2Kpluspi0Merged = false
+* const bool Kst2Kspiplus = KshortChannel
+* const bool smallSample = false
+
+### Paths
+* const std::string path_to_data = "/auto/data/dgerick/B2Kstmumu"
+* const std::string path_to_output_KshortPiplus = "/auto/data/dgerick/B2Kstmumu"
+* const std::string path_to_output_KplusPizero = "/home/lhcb/kopecna/B2KstarMuMu/data"
+* const std::string thePath = Kst2Kspiplus ? path_to_output_KshortPiplus : path_to_output_KplusPizero
+### TODO
+* const bool UseDTF = KshortChannel ? true : true //use mass and momenta of particles obtained by DTF?
+* const bool UsePDGmIn4Dvectors = KshortChannel ? false : true
+* const bool ReweightByRefChannel = KshortChannel ? true : true
+* const bool AlwaysUseRefChannelData = KshortChannel ? true : true
+* const bool Use2017 = KshortChannel ? true : true
+* const bool Use2018 = KshortChannel ? true : true
+* const bool UsePIDgen = KshortChannel ? true : false
+* const UInt_t firstnBins = 75
+* const UInt_t secondnBins = 50
+* const double firstMCrange[2] = {-0.5
+* 149.5}
+* const std::string firstMClatex = KshortChannel ? "N^{Long}_{Tracks}" : "N^{Long}_{Tracks}"
+* const bool UseBKGCAT = KshortChannel ? true : false
+* const bool pi0TM = true
+* const bool gammaTMdefault = false
+* const bool SplitInQ2 = KshortChannel ? false : false
+* TString TMVAmethod =  KshortChannel ? "BDTG" : "MLP" //used TMVA method
+* int verboseLevel = 2
+
+### Cuts
+* const double cut_DTF_status = 0.5
+* const double cut_DTF_chi2   = 200.0
+* const double K_star_plus_mass_difference = 100.0 //MeV  //difference between K* mass and PDG K* mass
+* const double cut_B_plus_M_low_basic = 4900.  //MeV
+* to account for cut on non-DTF mass in stripping
+* const double cut_B_plus_M_low =  KshortChannel ? 5150. : 5150.  //MeV
+* const double cut_B_plus_M_high = KshortChannel ? 6000. : 5800.  //MeV
+* const double B_plus_M_signal_window = KshortChannel ? 50. : 100 	//MeV -- as PDGmass +/- the window
+* const double cut_kin_Theta_low = 0.0
+* const double cut_kin_Theta_up = 0.4
+* const double cut_K_star_plus_pT = 1350.0 //MeV
+* const double cut_B_plus_pT = 2000.0 //MeV
+* const double cut_pi_zero_PT = 800.0 //MeV
+* const double cut_B_plus_DIRA = 0.99996
+* const double cut_muon_angle = 0.001
+* const double cut_B_plus_FDCHI2 = 121.0
+* const double cut_K_star_plus_FDCHI2 = 9.0
+* const double cut_B_plus_ConePTasym = -0.5
+* const double cut_DiMuon_M = 7100.0
+* const double cut_muon_IPCHI2_OWNPV_low  = 9.0
+* const double cut_IPCHI2_OWNPV_low  = 0.0
+* const double cut_IPCHI2_OWNPV_high = 12.0
+* const double cut_muon_ProbNNmu = 0.25
+* const double cut_K_plus_ProbNNk = 0.25
+* const double cut_gamma_CL = 0.15
+
+
+# Classes
+
+* struct SignalFitParameters{ //TODO: possibly add the K1 to ParamValues at some point
+* struct FixedParameters{//TODO: possibly add the K1 to ParamValues at some point
+* struct stat buffer
+
+# Functions and their parameters:
+### bool GetKSfromExePath()
+
+* **Parameters**
+    * 
+* **Return**
+
+### if ()
+
+* **Parameters**
+    * getcwd
+* **Return**
+
+### const bool KshortChannel = GetKSfromExePath()
+
+* **Parameters**
+    * ); // false := pi0 channel
+* **Return**
+
+### std::string TheDecay = ()
+
+* **Parameters**
+    * KshortChannel ? "KshortPiplus" : "KplusPi0Resolved");
+* **Return**
+
+### const std::string firstMCweight = ()
+
+* **Parameters**
+    * KshortChannel ? "nLongTracks" : "nLongTracks");
+* **Return**
+
+### const std::string seconMCweight = KshortChannel ? ()
+
+* **Parameters**
+    * UseDTF ? "B_plus_PT_DTF" : "B_plus_PT") : 
+* **Return**
+
+### const double seconMCrange[2] = {TMath::Log()
+
+* **Parameters**
+    * 1500.)
+    * TMath::Log
+* **Return**
+
+### const std::string seconMClatex = KshortChannel ? "p_{T}()
+
+* **Parameters**
+    * B^{+}) 
+* **Return**
+
+### const std::string TMtag = UseBKGCAT ? "" : "_TM" + std::string()
+
+* **Parameters**
+    * pi0TM ? "" : "_noPi0TM");
+* **Return**
+
+### const std::string TMbranch = "TMed" + std::string()
+
+* **Parameters**
+    * UseBKGCAT ? "BKGCAT" : "") + std::string
+* **Return**
+
+### const std::string gammaTMbranch = "TM_gammas" + std::string()
+
+* **Parameters**
+    * pi0TM ? "" : "_noPi0");
+* **Return**
+
+### bool isTM()
+
+* **Parameters**
+    * std::string customTMbranch
+    * bool TM
+    * bool gammaTM
+    * int gTM
+* **Return**
+
+### inline bool exists_test ()
+
+* **Parameters**
+    * const std::string& name
+* **Return**
+
+### return ()
+
+* **Parameters**
+    * stat 
+* **Return**
+
+### //verboseLevel 4: cout Errors ()
+
+* **Parameters**
+    * always cout errors)
+* **Return**
+
+### void setVerboseLevel()
+
+* **Parameters**
+    * int level
+* **Return**
+
+### bool set_gErrorIgnoreLevel()
+
+* **Parameters**
+    * 
+* **Return**
+
+### bool get_gErrorIgnoreLevel = set_gErrorIgnoreLevel()
+
+* **Parameters**
+    * );
+* **Return**
+
+### void coutTest()
+
+* **Parameters**
+    * string message
+* **Return**
+
+### void coutDebug()
+
+* **Parameters**
+    * string message
+* **Return**
+
+### void coutInfo()
+
+* **Parameters**
+    * string message
+* **Return**
+
+### void coutWarning()
+
+* **Parameters**
+    * string message
+* **Return**
+
+### void coutERROR()
+
+* **Parameters**
+    * string message
+* **Return**
+
+### double get_cut_B_plus_M_low()
+
+* **Parameters**
+    * int Run = 2 //setting lower B mass cut for different run/years
+* **Return**
+
+### double get_cut_B_plus_M_low()
+
+* **Parameters**
+    * std::string year = "2015"
+* **Return**
+
+### bool IsDST()
+
+* **Parameters**
+    * std::string year
+    * bool MC
+    * bool ReferenceChannel
+    * bool PHSP
+* **Return**
+
+### void addLHCbtag()
+
+* **Parameters**
+    * Float_t x = 0.6
+    * Float_t y = 0.85
+    * std::string suffix = ""
+    * Int_t color = 1
+    * Float_t Scaling = 1.0
+* **Return**
+
+### std::string Float2Comma()
+
+* **Parameters**
+    * Float_t f
+    * Int_t d
+* **Return**
+
+### if()
+
+* **Parameters**
+    * d > 5
+* **Return**
+
+### if()
+
+* **Parameters**
+    * d <= 0)return out;
+* **Return**
+
+### for()
+
+* **Parameters**
+    * Int_t n = 0; n <= d; n++
+* **Return**
+
+### N[n] = static_cast<Int_t>()
+
+* **Parameters**
+    * 
+* **Return**
+
+### for()
+
+* **Parameters**
+    * Int_t m = 0; m < n; m++)N[n] -= N[m] \* TMath::Power
+* **Return**
+
+### if()
+
+* **Parameters**
+    * n != 0)out.append
+* **Return**
+
+### bool replace()
+
+* **Parameters**
+    * std::string& str
+    * const std::string& from
+    * const std::string& to
+* **Return**
+
+### void printVector ()
+
+* **Parameters**
+    * std::vector<int> vector
+* **Return**
+
+### void printVector ()
+
+* **Parameters**
+    * std::vector<double> vector
+* **Return**
+
+### void printVector ()
+
+* **Parameters**
+    * std::vector<float> vector
+* **Return**
+
+### void printVector ()
+
+* **Parameters**
+    * std::vector<string> vector
+* **Return**
+
+### std::vector<string> yearsMC ()
+
+* **Parameters**
+    * bool Reference = false
+    * bool PHSP = false
+    * int Run = 1
+* **Return**
+
+### std::vector<string> yearsData ()
+
+* **Parameters**
+    * int Run = 1
+* **Return**
+
+### std::vector<string> yearsVector ()
+
+* **Parameters**
+    * bool MC = false
+    * bool Reference = false
+    * bool PHSP = false
+    * int Run = 1
+* **Return**
+
+### std::vector<int> yearsVectorInt()
+
+* **Parameters**
+    * bool MC = false
+    * bool Reference = false
+    * bool PHSP = false
+    * int Run = 1
+* **Return**
+
+### std::vector<string> yearsInc()
+
+* **Parameters**
+    * int Run
+* **Return**
+
+### std::vector<string> yearsBkgMC()
+
+* **Parameters**
+    * bool Reference
+    * bool B0
+    * bool K1
+    * bool Inc
+    * int Run
+* **Return**
+
+### int getRunID()
+
+* **Parameters**
+    * string year
+* **Return**