kopecna/EWP-BplusToKstMuMu-AngAna
kopecna/EWP-BplusToKstMuMu-AngAna
2
1
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity
Page: CompareVariables.cpp
Pages
AddVariable.cpp B mass model B2Kstmumu B2Kstmumu.pro BDTSelection BDTcutScanner.cpp BackgroundPdf Class BackgroundType Class CompareVariables.cpp Comparison tool Design.cpp Efficiency Class Efficiency.cpp Event generator FCNC Fitter Ganga GlobalFunctions.hh Home InclusiveCheck.cpp Input MVA variables MCtruthmatching.cpp MVA Class MassFit.cpp ParamValues Class Paths.hpp PlotTMVA.cpp RemoveMultipleCandidates.py ReorganizeTGraph.py Rescale pi0 momentum Resolution Plotting Running the FCNC fitter Selection code Signal Study SignalPdf Class SignalType Class Standalone Scripts TMVA Class application TMVA TMsource Utils.hpp Utils3.py flav.io getPathForPython nTrackWeights.cpp
3 CompareVariables.cpp
Renata Kopecná edited this page 2022-02-10 11:29:09 +01:00
Table of Contents

Older function created before my time, so you will find new lines without brackets here and other lovely features.

Table of Contents

Global variables:

  • const bool CutInQ2 = false //If true, select only events between minQ2 and maxQ2
  • double minQ2 = 15.0
  • double maxQ2 = 100.0
  • const bool CutInMuPt = false //If true, select only events between minPT and maxPT
  • double minPT = 1100.0 //2500
  • double maxPT = 2500.0 //10000
  • const bool separateCharge = false //Check for B+ and B- separately
  • bool plusCharge = false //If separateCharge, do you want to plot B+ or B-?

Functions

int evalOperation()

Before I discovered that ROOT can take formulas, I had to figure it out on my own. So this takes x and y and performs the desired operation. It is useful when comparing eg momentum-momentum or 1-log(momentum).

  • Parameters
    • bool useTwoVars = false
    • int x = 0
    • int y = 0
    • std::string operation = "plus"
    • bool bLogVar = false
    • bool bLog1minusVar = false
  • Return
    • Operation(x,y)

double evalOperation()

Before I discovered that ROOT can take formulas, I had to figure it out on my own. So this takes x and y and performs the desired operation.

  • Parameters
    • bool useTwoVars = false
    • double x = 0.0
    • double y = 0.0
    • std::string operation = "plus"
    • bool bLogVar = false
    • bool bLog1minusVar = false
  • Return
    • Operation(x,y)

int compareVariables()

Compute and plot a desired variable for given year and subdecay channel. Choose weither it is sWeighted, not weighted, weighted by nTracks or 2 * 1D reweighted by first and second weighting variables defined in GlobalFunctions.hh. The data sample can be sweighted, so the weighting only concerns the MC sample. Choose whether to print the plots to pdf and root files. Variables can be plotted in log(var) or log(1 - var), using the according booleans.

  • Parameters
    • int year = 2011
    • std::string variable = "B_plus_PT"
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2Dweighted = true
    • bool bLogVar = false
    • bool bLog1minusVar = false
    • bool bPrint = true
    • bool KshortDecayInVelo = true
    • std::string weightBranch = "CB"
    • bool useTwoVars = false
    • string operation = ""
    • string variable2 = ""
    • string xAxis = ""
  • Return
    • 0 if something fails, 1 otherwise

int compareVariablesKplusKshort()

Seems to do the same thing as compareVariables(). Probably obsolete.

  • Parameters
    • int year = 2011
    • std::string variable = "B_plus_PT"
    • bool sWeightKplus = false
    • bool sWeightKshort = true
    • bool bWeighted = true
    • bool b2Dweighted = true
    • bool bLogVar = false
    • bool bLog1minusVar = false
    • bool bPrint = true
    • bool KshortDecayInVelo = true
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKshorts()

Calls compareVariables() for all variables needed by the KS channel.

  • Parameters
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2DWeighted = true
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKshortsWeightedAndUnweighted()

Calls compareAllKshorts(), Run 1, with no weights, 1D weights and 2D weights.

  • Return
    • Always 1

int compareKshortPreSelectionAndStrippingVariables()

Calls compareVariables() for all pi+ PIDK and mu+- ProbNNmu.

  • Parameters
    • int Run = 1
  • Return

int compareKshortMVAonly()

Calls compareVariables() for all variables used in MVA training defined in [[TMVA_variables_TheDecayDL.txt|Input-MVA-variables]].

  • Parameters
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2DWeighted = true
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareMuIPKplusResolved()

Calls compareAllKshorts() for mu+- IP chi2.

  • Parameters
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKplusResolved()

Calls compareAllKshorts() for all useful variables in the K+pi0 channel.

  • Parameters
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2DWeighted = false
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKplusMerged()

Obsolete, was checking the merged pi0 subchannel.

  • Parameters
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2DWeighted = false
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareKplusResolvedMVAonly()

Calls compareVariables() for all variables used in MVA training defined in [[TMVA_variables_TheDecayDL.txt|Input-MVA-variables]].* used by the K+pi0 channel.

  • Parameters
    • bool sWeight = true
    • bool bWeighted = true
    • bool b2DWeighted = true
    • int Run = 1
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKplus()

Calls compareAllKplusResolved() for both runs, with different weights and specified weightBranch.

  • Parameters
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAll()

Calls either compareKplusResolvedMVAonly() or compareKshortMVAonly()

  • Parameters
    • int Run = 1
  • Return
    • 0 if something fails, 1 otherwise

int compareTrackMultiplicities()

Comparison of track multiplicities between Run I and Run II.

  • Parameters
    • bool KshortDecaysInVelo = false
    • bool MC = false
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAnglesInPHSP()

Comparison of the angles between (2D-)weighted and non-weighted PHSP MC.

  • Parameters
    • int year = 2011
    • bool sWeight = true
    • bool b2Dweighted = true
    • bool KshortDecayInVelo = true
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllAnglesInPHSP()

Calls compareAnglesInPHSP() for all years and LL/DD.

  • Parameters
    • int Run = 1
  • Return
    • 0 if something fails, 1 otherwise

int compareAnglesInMC()

Cmparison of the ngles between (2D-)weighted and non-weighted signal MC.

  • Parameters
    • int year = 2011
    • bool sWeight = true
    • bool b2Dweighted = true
    • bool KshortDecayInVelo = true
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllAnglesInMC()

Calls compareAllAnglesInMC() for all years and LL/DD.

  • Parameters
    • int Run = 1
  • Return
    • 0 if something fails, 1 otherwise

int compareAnglesInPHSPbeforeAndAfterBDTcut()

Cmparison of the angles in PHSP MC before and after the MVA cut.

  • Parameters
    • int year = 2011
    • bool bweighted = true
    • bool b2Dweighted = true
    • bool KshortDecayInVelo = true
    • std::string weightBranch = firstMCweight
  • Return
    • 0 if something fails, 1 otherwise

int compareAllPHSPAnglesBeforeAndAfterBDT()

Calls compareAnglesInPHSPbeforeAndAfterBDTcut() for all years and LL/DD.

  • Parameters
    • int Run = 1
  • Return
    • 0 if something fails, 1 otherwise

int compareNTracksAndBplusPT()

Comparison of nTracks and B_plus_PT between signal and reference channel MC.

  • Parameters
    • int year = 2011
    • bool KshortDecayInVelo = true
  • Return
    • 0 if something fails, 1 otherwise

int compareAllNtracksAndBplusPT()

Calls compareNTracksAndBplusPT() for all years and LL/DD.

  • Parameters
    • int Run = 1
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKplusKshort()

Comparison of Kplus and Kshort.

  • Parameters
    • bool sWeightKplus = false
    • bool sWeightKshort = true
    • bool bWeighted = true
    • bool b2DWeighted = false
    • int Run = 1
    • bool KshortDecayInVelo = true
    • std::string weightBranch = "nTracks"
  • Return
    • 0 if something fails, 1 otherwise

int compareAllKplusAllRunsAndWeights()

Calls compareAllKplusKshort() for both runs and not/using weight sin K+ and/or KS.

  • Parameters
    • bool bWeighted = true
    • bool b2DWeighted = false
  • Return
    • 0 if something fails, 1 otherwise