kopecna/EWP-BplusToKstMuMu-AngAna
kopecna/EWP-BplusToKstMuMu-AngAna
2
1
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Update 'FCNC Fitter'

Renata Kopecná 2022-02-04 14:20:13 +01:00
parent a45612d06e
commit 06cecf4a15
1 changed files with 35 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
FCNC-Fitter.md

@ -110,4 +110,38 @@ Currently, the file structure is as follows. All the `.cc` source files also hav
## Running ## Running
python script For the detailed explanation of possible parameters, [[click here|Running-the-FCNC-fitter]].
For quickly running everything, run the following commands
```
python run.py -convert -all -Run 1
python run.py -dontCompile -convert -all -Run 2
python run.py -dontCompile -MC -angCorr -Run 1 -scan
python run.py -dontCompile -MC -angCorr -Run 2 -scan
python run.py -dontCompile -MC -angCorr -Run 1
python run.py -dontCompile -MC -angCorr -Run 2
python run.py -dontCompile -MC -angRes -Run 1
python run.py -dontCompile -MC -angRes -Run 2
python run.py -dontCompile -MC -Run 12 -fit -Ref -nBins 1
python run.py -dontCompile -MC -Run 12 -fit -nBins 5
python run.py -dontCompile -Run 12 -v 2 -fit -Ref -nBins 1 -genMC -Ref
python run.py -dontCompile -Run 12 -v 2 -fit -nBins 5 -genMC
python run.py -dontCompile -Data -Run 12 -fit -Ref -nBins 1 -massDim
python run.py -dontCompile -Data -Run 12 -fit -nBins 5 -massDim
python run.py -dontCompile -Data -Run 12 -fit -Ref -nBins 1 -bkgOnly -upper
python run.py -dontCompile -Data -Run 12 -fit -nBins 5 -onlyBkg -upper
python run.py -dontCompile -Run 12 -v 2 -fit -Ref -nBins 1
```
**TODO** Create PHSP weights
* First, it compiles the code and prepares the tuples with Run 1.
* Then, since we already compiled the code, we use the -dontCompile option for running the next commands.
* We prepare the tuples for Run2. Then, it scans the order of polynomial used for the PHSP angular description.
* Next step is to create the angular acceptance weights using the optimal polynomial.
* Then, the signal Monte Carlo sample is fitted: first fit the reference channel, then the signal channel.
* After this step is done, it fits the generator level Mont sample in the J/psi Q2 bin and rare channel in 5 Q2 bins.
* The next step is to fit the data in the B+ and K*+ masses: first, fit the J/psi channel, then the rare channel in the 5 bins.
* Then, fit the angular background distributions (upper B+ mass sideband) using a chebyschev polynomial.
* The last step is to fit the data, J/psi channel in 4 dimensions (B+ mass, cos(thetal), cos(thetak), phi)