Skip to content

Latest commit

 

History

History
40 lines (35 loc) · 1.4 KB

README.md

File metadata and controls

40 lines (35 loc) · 1.4 KB
Raw

tth_charge-flip_estimation

Installation

See https://github.com/HEP-KBFI/tth-htt/wiki/Setup#for-auxiliary-measurements

Running

  1. Create datacards for data and pseudodata:
create_pseudodata_datacard.py -i prepareDatacards_charge_flip_mass_ll.root \
  -o prepareDatacards_data_charge_flip_mass_ll.root \
  -O prepareDatacards_pseudodata_charge_flip_mass_ll.root \
  -t ele -s data_obs DY DY_fake Singletop Diboson TTbar
  1. Create text files from both datacards:
ChargeFlipDC -i prepareDatacards_data_charge_flip_mass_ll.root \
  -b DY_fake Singletop Diboson TTbar \
  -o fit_results_data -y 2016 -e -f

ChargeFlipDC -i prepareDatacards_pseudodata_charge_flip_mass_ll.root \
  -b DY_fake Singletop Diboson TTbar \
  -o fit_results_pseudodata -y 2016 -e -f
  1. Fit the 21 electron pair categories and make plots:
make_fits.py -i fit_results_data       -o fit_results_data/results_cat.txt
make_fits.py -i fit_results_pseudodata -o fit_results_pseudodata/results_cat.txt
  1. Create prefit and postfit plots:
make_fit_plots.py -i fit_results_data
make_fit_plots.py -i fit_results_pseudodata
  1. To obtain flip rates and create various pull plots:
plot_pulls.py -f fit_results_data -F fit_results_pseudodata -i hadd_stage2_Tight.root -o $PWD -l

Use fit_result_data_exclusions.root in the analysis to assign charge misidentification probabilities to the data events.