EWP-BplusToKstMuMu-AngAna

kopecna/EWP-BplusToKstMuMu-AngAna

Fork 0

Table of Contents

Table of Contents
Introduction
Backup of the tuples and code
Running the code
Q2 binning

In the selection
In the fit

Available Monte Carlo samples

Signal Monte Carlo
Reference Monte Carlo
Phase-space Monte Carlo
Background Monte Carlo

Feynman graphs

Introduction
Backup of the tuples and code
Running the code
Q2 binning
- In the selection
- In the fit
Available Monte Carlo samples
Feynman graphs

The Tables of Content in this wiki were created by gh-md-toc

The detialed markdown of functions and classes was created by Extract_functions_with_parameters.py. It takes either C++ or python files as input and prints out the functions with their parameters. It also adds link to the function definitions as comments (handy when you need to reference them quickly). The script is definintely not perfect and relies on the user using 4-space tab (also in the C++ code!). However, when one needs to create a markdown, it is handy.

Introduction

This wiki descripes the repository containing the code that was used for the analysis of B+->K*+mumu, K*+->K+pi0. The analysis started as a joint analysis with the decay B+->K*+mumu, K*+->Kshortpi+. A lot of tweaks and changes were made on the way with little time to make the code pretty: be patient.

The FCNC Fitter used for the angular analysis have been updated, cleaned up and it got a new library for logging. It was created by Ch. Langenbruch ages ago, then tweaked significantly by David and then by me. It is included in case some cross-check is needed.

The details are described in my thesis

Backup of the tuples and code

In order to be able to perform cross-checks and just to have a general backup, the tuples are saved in Heidelberg at /auto/data/kopecna/B2KstMuMu/data/. The code is stored here, however it was developed here: https://gitlab.cern.ch/LHCb-RD/ewp-Bplus2Kstmumu-AngAna/-/tree/work_renata

The tuples and some other files used in the KshortPiplus analysis are saved at /auto/data/dgerick/B2Kstmumu/. This was created with the code on the master branch in the ewp-Bplus2Kstmumu-AngAna project on gitlab: https://gitlab.cern.ch/LHCb-RD/ewp-Bplus2Kstmumu-AngAna/-/tree/master

All the PDFs with meeting contributions are saved in

/auto/data/kopecna/B2KstMuMu/meetings

Some files related to the toys I used are saved in

/auto/data/kopecna/B2KstMuMu/data/Toys/

Running the code

Q2 binning

In the selection

The range of the Q2 is defined in Utils.hpp for the Jpsi region and the rare regions.

The Jpsi Q2 region is 8.68-10.09 GeV^2

The rare Q2 region is 0-0.98, 1.1-8.0, 11.0 - 12.5, 15.0-20 Gev^2

The rare region is dived so to remove the contribution of the phi, Jpsi and psi2s resonances.

The binning in Q2 is defined in TMefficiency class for the selection part, the name tag is q2_binned_fit. The name of the class was chosen when testing how different variables behave under truth matching. Since the class was generally useful, it became the default.

They are set to be divided into five rare bins: * 0.1 - 4.0 GeV^2 with the exclusion of 0.98-1.1 GeV^2 (phi resonance) * 4.0 - 8.0 GeV^2 * 11.0 - 12.5 GeV^2 * 15.0 - 20.0 GeV^2

In the fit

The ranges are defined in constant.cc. Different binning options are defined there in the functions get_TheQ2binsmin(int nBins, bool Reference) and get_TheQ2binsmax(int nBins, bool Reference). The typically used bins are

Jpsi
- 8.68-10.09 GeV^2
Rare
- 0.25 - 4.0 GeV^2 with the exclusion of 0.98-1.1 GeV^2 (phi resonance)
- 4.0 - 8.0 GeV^2
- 11.0 - 12.5 GeV^2
- 15.0 - 18.0 GeV^2
- 1.1 - 6.0 GeV^2 (extra wide bin very much below J/psi, useful for theorists)

When running the fit, the default number of bins is -1, hence the number of desired bins has to be set by using the parameter -nBins <int>!

Available Monte Carlo samples

Signal Monte Carlo

Decay number: 12113100

Available years: 2011, 2012, 2016, 2017, 2018

The mising 2015 is caused by missing Calo tables. The problem is likely in the MC production where the CaloHypo2MC table was not stored as it should. Same problem in 2016, but somehow the rawEvent is available so CaloHypo2MC is recovered and we can use 2016.