See https://github.com/HEP-KBFI/tth-htt/wiki/Setup#for-auxiliary-measurements
- 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- 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- 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- Create prefit and postfit plots:
make_fit_plots.py -i fit_results_data
make_fit_plots.py -i fit_results_pseudodata- 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 -lUse fit_result_data_exclusions.root in the analysis to assign charge misidentification probabilities to the data events.