Higgsino
The NUHM2 Higgsino study
- Add jupyter notebook
v55. New functionfor mll shape plot
- Add new function in
truth3_mll_study.py
compare_mll_shape()
v54. New scripts for making plots
- Add Judita's version cross-section file
NUHM2_Weak_Judita.txt
- Add new scripts:
plot_METOverHT_vs_mll.py
plot_eemumu_invariant_mass.py
plot_VRSS.py
- Modify the scripts for making plots in the dissertation
plot_distributions_in_SR.py
plot_truth3_kinematics.py
truth3_mll_study.py
upperlimit_xsec.py
reweight_mc_event_weight.C
reweight_truth3_mll.C
v53. New script for TMVA distributions
- Add
plot_TMVA_BDT_distributions.py
v52. Modify script
- Modify
reweight_truth3_mll.C
so it can make Higgsino and NUHM2 mll comparison plot as shown in Fig1 in paper.
v51. Align the files on lxplus
- Add new files:
ytHiggsinoAnalysis/misc/NUHM2_Strong.txt
ytHiggsinoAnalysis/misc/NUHM2_Strong_SUSYTools.txt
ytHiggsinoAnalysis/misc/NUHM2_Weak_SUSYTools.txt
- Add new folders:
ytHiggsinoAnalysis/misc/NUHM2/
ytHiggsinoAnalysis/misc/StrongSLHA/
ytHiggsinoAnalysis/misc/WeakSLHA/
- Modify script:
ytHiggsinoAnalysis/pythons/plot_truth3_kinematics.py
v50. Add new csv files
- Add new csv files:
NUHM2_all_points_upperlimit_upperlimit_only_MCprod_NoSys_doNewFFs.csv
NUHM2_all_points_upperlimit_upperlimit_only_MCprod_AllSys_doNewFFs.csv
upperlimit_xsec.py
: change input file
v49. Add new csv files
- Add new csv files:
NUHM2_all_points_upperlimit_upperlimit_only_Judita_MCprod_AllSys_doNewFFs.csv
NUHM2_all_points_upperlimit_upperlimit_only_reweighting_AllSys_doNewFFs.csv
NUHM2_all_points_upperlimit_upperlimit_only_reweighting_NoSys_doNewFFs.csv
upperlimit_xsec.py
:- re-draw y-axis to show ticks
- Instead of using
j = i + 1
, usem12
to find the correct element in list
v48. Script for upper limit plot
AtlasStyle.py
: Set fontNUHM2_Weak.txt
: New SUSYTools cross-section fileupperlimit_xsec.py
: New script to make cross-section upper limit plot.
v47. Modify the cutflow to fit the internal note
yt_cutflows.h
andyt_cutflows.cxx
: Add codes for cutflow table in the note.ytEventSelection.cxx
:- Add codes for cutflow table in the note.
- Add
bTagWeight
andFFWeight
in the weight calculation.
Run.cxx
andRun.sh
:- Use
/afs/cern.ch/work/j/jeliu/public/HiggsinoFitterTrees/v1_9/AllSignals_SusySkimHiggsino_v1.9b_SUSY16_tree_AllSys.root
.
- Use
v46. Modify script
truth3_mll_study.py
:- Fix the typo in efficiency
- Use the accurate value for BRs
plot_distributions_in_SR.py
:- Use thicker line size for NUHM2
- NUHM2 yields are multiply by 10 before making plot
v45. Modify script
plot_distributions_in_SR.py
: modify histogram xmax range
v44. Modify scripts
reweight_mc_event_weight.C
: change plot stylereweight_truth3_mll.C
:- rename
reweight_truth3.C
toreweight_truth3_mll.C
- change plot style
- rename
truth3_mll_study.py
: change plot styleplot_distributions_in_SR.py
: add SFOS results
v43. Modify scripts
Nevents.py
: change the output formatcompare_files.py
: Add commentsget_mass_difference.py
: Remove extra spaceget_yields.py
: Remove extra spaceplot_distributions_in_SR.py
: Add commentsplot_truth3_kinematics.py
: Add commentsread_Xsec.py
: Add commentstruth3_mll_study.py
: Add commentsreweight_mc_event_weight.C
: Change input pathreweight_truth3.C
: Change input path
v42. Move old scripts to obsolete directory
misc/obsolete/
: Create new directory and move the following files to this directory- from
python/
:BR_calculation.py
fit_common_ntuple.py
make_weight_plot.py
plots_making.py
reweight.py
show_AMI_info.py
test_fit.py
truth3_plots_comparison.py
- from
scripts/
:FuncMllDistr.C
file_path.C
fit_common_ntuple.C
mll_ratio.C
overlap_histograms.C
plot_reweight_common_ntuple_mll.C
reweight_common_ntuple.C
reweight_kinematic_variables.C
test_120_100.C
test_Higgsino_160_100.C
test_Higgsino_190_150_m12_600.C
test_fit.C
view_combine_hist.C
- from
v41. Add new cutflow results
make_cutflow.py
: modify scriptcutflows_20171027/
: renamecutflows/
tocutflows_20171027/
cutflows_20171028/
: new resultscutflows_20171029/
: new resultsplot_distributions_in_SR.py
: change outptu digitssignificance.py
: Use new results for nsig and nbkg
v40. modify python scripts
fit_common_ntuple.py
: Use new input filestest_fit.py
: Change the input pathAtlasStyle.py
: Change the sapce for input stringtruth3_mll_study.py
:- fix typo
- Change the input path, y title and draw option
v39. Add new cutflow results
ytHiggsinoAnalysis/misc/HF/
: New results for cutflows
v38. modify python scripts
plot_distributions_in_SR.py
: mergeplot_N_minus_one.py
into this scriptplot_N_minus_one.py
: delete this scripttruth3_plots_comparison.py
: Use new input files
v37. Add new script
plot_N_minus_one.py
: new script for N-1 plots
v36. modify python scripts
AtlasStyle.py
: change font sizeBR_calculation.py
: change input pathcompare_files.py
: Add new line at the end of fileget_mass_difference.py
: change the output formatNevents.py
: change the BR, Eff, weighted events, N_SR valuesytEventSelection.cxx
: Use b-jet pT>20add_weight_to_common_ntuple.C
: remove empty lineRun.sh
: remove empty line
v35. Add new python scripts:
get_truth3_information.py
: new script to get TRUTH3 informationplot_distributions_in_SR.py
: new script to make kinematic distribution in SRplot_truth3_kinematics.py
: new script to make kinematic distributions using official TRUTH3show_AMI_info.py
: new script to get filter efficiency
v34. modify python scripts
BR_calculation.py
: modify the commentscompare_files.py
: renamecheck_slha.py
and add new functionget_yields.py
: add SRee, SRmm, SRSF yields and weighted yieldsspectrum.py
: Bug fixedmake_weight_plot.py
: change the input pathplots_making.py
: change the input path and filesreweight.py
: change the input path and files
v33. Update to v1.9 and v1.9b ntuples
ytEventSelection.h
: Use v1.9b ntupleytEventSelection.cxx
: Use v1.9b ntuple and add weight to optimize SR histograms.Run.cxx
: Use v1.9 and v1.9b ntuplesRun.sh
: Use v1.9 and v1.9b ntuples
v32. Modify python scripts
BR_calculation.py
: modify script to run SLHA and SusyHitNevents.py
: Use official TRUTH3 SimpleAnalysis resultscheck_slha.py
: Add commons and fix typoget_mass_difference.py
: Add commonlist_trees.py
: Add commons and change to v1.9cutflow_v1.9.txt
: Add a new log file.file_path.C
: Use the official TRUTH3 SimpleAnalysis results
v31. New TRUTH3
BR_calculation.py
: Add code to read SusyHit output.RunSimpleAnalysis.sh
: Add new official TRUTH3 samples.EwkNUHM22016.cxx
: Use new SR, CR, and VR definitions
v30. Add Higgsino optimised SR
show_normalization_factor.C
: fix bugsplot_reweight_common_ntuple_mll.C
: change output file nameget_yields.py
: Change pathmake_weight_plot.py
: New scriptyt_regions.h
: Add new functionsyt_regions.cxx
: Implement new functionsytEventSelection.h
: Add new histogramsytEventSelection.cxx
:- Implement new histograms
- Add common cut, 2l channel cut, and Higgsino SR cuts
- B-jet uses
nBJet20_MV2c10
- Fill yields to new histograms
v29. add new branches to common ntuple
add_weight_to_common_ntuple.C
: New scriptplot_reweight_common_ntuple_mll.C
: New scriptget_mass_difference.py
: modify outputJoana_reweight.cxx
: Add DSID for C1C1, N2C1p, and N2C1mreweight_mc_event_weight.C
: Add weight when fill histogram
v28. reweight to luminosity
fit_common_ntuple.py
: move tomisc/
test_fit.py
: move tomisc/
Joana_reweight.cxx
: fix typoreweight_truth3.C
: removeTChain.h
show_normalization_factor.C
: fix typoreweight_kinematic_variables.C
: remove extra spaces.ytUtility.C
:get_cross_section()
: New functionget_filter_efficiency()
: New functionget_luminosity()
: New functionget_sum_of_event_weight()
: New functionget_mc_weight()
: New function
reweight.py
: New scriptreweight_mc_event_weight.C
: New script. This is based onreweight_kinematic_variables.C
v27. minor corrections
Xsec.pdf
: Change legend position and move tolog/
yields.txt
: Remove file, but add two files (yields_v1.6.txt and yields_v1.8.txt) inlog/
BR_calculation.py
: Add new functionslist_branching_ratio()
andget_slha_file()
Nevents.py
: Add new functiontotal_number_of_events_produced()
and change legend position.check_slha.py
: New end of lineread_Xsec.py
: Add new functioncross_section_ratio()
significance.py
: print informationtest_fit.py
: Use new directory20170817/
Run.sh
: Use v1.8b samplesfit_common_ntuple.C
: Change legend positionshow_normalization_factor.C
: Use new mappingRun.cxx
: Use v1.8b samplestest_Higgsino_170_150_m12_600.C
: Change name totest_Higgsino_190_150_m12_600.C
ytUtility.C
: New end of linetest_fit.C
: Remove unnecessary linestest_fit.py
: Remove unnecessary parts.
v26. v1.8 yields
yields.txt
: Add v1.8 yieldsget_yields.py
: Fix the v1.8 path
v25. Update to common ntuple v1.8
ytEventSelection.cxx
: Fix the bus error- All vector container pointers need to set to zero at the beginning
v24. New EwkNUHM22016.cxx
EwkNUHM22016.cxx
: Checkout the newest version from simple analysis.fit_common_ntuple.py
:- remove ";" in some lines
- Use Lorenzo's modified function
file_path.C
: Use new directory20170817/
mll_ratio.C
: Use new directory20170817/
overlap_histograms.C
: Use new directory20170817/
test_fit.C
: Use new directory20170817/
test_Higgsino_160_100.C
: Use new directory20170817/
view_combine_hist.C
: Use new directory20170817/
reweight_truth3.C
:- Use new directory
20170817/
- Add
preselection
and2LChannel
- Use new directory
show_normalization_factor.C
: Use Lorenzo's modified function
v23. Use 20170815 results and add preselection and 2L requirements
EwkNUHM22016.cxx
: New code modified fromEwkHiggsino2016
RunSimpleAnalysis.sh
: ChangeEwkHiggsino2016
toEwkNUHM22016
file_path.C
: Use new directory20170815/
get_combined_hist.C
: ChangeEwkHiggsino2016
toEwkNUHM22016
mll_ratio.C
: ChangeEwkHiggsino2016
toEwkNUHM22016
overlap_histograms.C
: Use new directory20170815/
view_combine_hist.C
: Use new directory20170815/
test_fit.C
:- Change function declaration.
- Use new directory
20170815/
- Include
ytUtility.C
so I can useget_dm_NUHM2()
- Add
xrange_max
forSetRangeUser()
- Fix the line color
- Change the legend coordinate
test_Higgsino_160_100.C
: Use new directory20170815/
reweight_truth3.C
:- Use new directory
20170815/
- Remove
get_func_ratio_parameters()
- Change
EwkHiggsino2016
toEwkNUHM22016
- Add
preselection
and2LChannel
- Use new directory
reweight_kinematic_variables.C
:- Use new directory
20170815/
- Add
y_scale_factor
and changey_max
- Change the legend coordinate and font size
- Fix the output file name.
- Change
EwkHiggsino2016
toEwkNUHM22016
- Add
preselection
and2LChannel
- Add
pTLep1
andnJet30
plots - Change some cosmetic of plots
- Use new directory
v22. reweight kinematic variables
AtlasStyle.py
: Comment ROOT colorBR_calculation.py
: Addsum_BR()
intomain()
but comment itNevents.py
: Add commentsget_yields.py
: Add horizontal lines to separate the function definition.reweight_truth3.C
:- Bug fix in
plot_Lorenzo()
- Use bin-by-bin weight in
plot_Mike()
- Remove unnecessary TH1 in
tree_reader()
- Bug fix in
get_combined_hist.C
: Add new function overloadget_histogram()
reweight_truth3.C
: Change the order of n2 and n1 in the input arguments.view_combine_hist.C
: Add new script to check the distributionsreweight_kinematic_variables.C
: Add new script to check the distributions of reweighted kinematic variables.
v21. add tree_reader()
FuncMllDistr_modified.C
: add commentfile_path.C
: add spaceget_combined_hist.C
:- remove unnecessary part
- add if condition when there are only histograms in the input.
mll_ratio.C
: remove unnecessary partoverlap_histograms.C
: usecombine_histogram()
reweight_common_ntuple.C
: remove unnecessary partreweight_truth3.C
:- remove unnecessary part
- Add
tree_reader()
and use it.
test_Higgsino_160_100.C
: remove unnecessary part
v20. Correct Lorenzo's approach
get_combined_hist.C
:- add new function
combine_histogram()
- modify
get_combined_hist()
- add new function
fit_common_ntuple.C
: Change to modified functionreweight_common_ntuple.C
:- Use modified function
- Use
ytUtility.C
- Use function to reweight
reweight_truth3.C
:- Correct Lorenzo's approach
- Add Mike's approach
FuncMllDistr_modified.C
: add new scriptfile_path.C
:: add new scriptmll_ratio.C
:: add new scripttest_Higgsino_160_100.C
: add new scriptytUtility.C
: add new script
v19. Rewrite the reweighting script
show_normalization_factor.C
: add new scripttest_120_100.C
: add new scripttest_Higgsino_170_150_m12_600.C
: add new scriptFuncMllDistr.C
: Add commentRunSimpleAnalysis.sh
: Add Higgsino 190_150 samplesconvert_EVNT_to_TRUTH.sh
: Add Higgsino 190_150 samplesoverlap_histograms.C
: Change input pathfit_common_ntuple.C
:- Add Higgsino 190_150
- New fit plot method
reweight_truth3.C
: Complete rewrite this script.test_fit.py
: rewrite this script
v18. Add scripts for reweighting
test_fit.C
: Add new script. This works for FuncMllDistr.get_combined_hist.C
: Add new script.fit_common_ntuple.C
: Add new script. This works for FuncMllDistr.FuncMllDistr.C
: Add new script.GetWeightFunction.C
: Add Lorenzo's script as an example.test_fit.py
:- Add new functions
get_histogram()
,funcMllDistr()
,reweight_method_1()
, andreweight_method_2()
- Add new functions
overlap_histograms.C
: Add new script.reweight_common_ntuple.C
: Add new script.reweight_truth3.C
: Add new script.
v17. Add new Higgsino samples and scripts for these new samples
RunSimpleAnalysis.sh
: Add new Higgsino samplescheck_slha.py
: Add new scriptget_mass_difference.py
: Add new scriptXsec.pdf
: Cross section for N2N1, C1C1, N2C1p, N2C1mconvert_EVNT_to_TRUTH.sh
: modified for new Higgsino samplestest_fit.py
: Add new script
v16. Add new branches
NUHM2.h
: New TSelector for v1.8ytEventSelection.h
: Add new branchesytEventSelection.cxx
: Add new branches- Joey removed METOverHTLep12 on Jun/20. Fix
METOverHTLep12
toMETOverHT
v15. Add new scripts
NUHM2_DSID.tar
: Judita's JOBR_calculation.py
: new scriptNevents.py
: new scriptsum_BR.py
: remove this script (move the function toBR_calculation.py
)list_trees.py
: change to v1.8Run.sh
: Add v1.8 new samplesRun.cxx
: Add v1.8 new samplestruth3_plots_comparison.py
:- Read the cross-section, branching ratio, and filter efficiency from the files.
- Use 20170716 results to compare the plots
- Remove
met>50
cut because it already applied in MET50 filter. - Add new function
weight()
v14. Use new SR definitions
RunSimpleAnalysis.sh
: Add new samples for the truth studyconvert_EVNT_to_TRUTH.sh
: Add new scriptMGPy8EG_A14N23LO_Higgsino_2L.txt
: Add Higgsino simplified model fileread_Xsec.py
: Add new python script.EwkHiggsino2016.cxx
: Use new SR definitions
v13. Add new function
RunSimpleAnalysis.sh
: Add new samples for the truth studytruth3_mll_study.py
: Add new functioncompare_n2_decay()
v12. Add new script
sum_BR.py
: new script to calculate the summation of BR.spectrum.py
:Xsec_plots_weak()
function call with different arguments.truth3_mll_study.py
: new functioncompare_pythia_madgraph_madspin()
.
v11. Finish spectrum.py
spectrum.py
:- Modify
read_cross_section()
- Rename
Xsec_plots()
toXsec_plots_all_FS()
- Add new functions
Xsec_plots_weak()
,sum_of_Xsec()
,read_slha()
,mass_spectra()
, andmass_difference_spectra()
- Modify
v10. Variables are nitialized to -99999
truth3_plots_comparison.py
:- Change the samples for making plots.
- Add met > 50 cut for NUHM2
- Remove most of other cuts and only keep HT30, HTIncl, and mT2 because they have 0
EwkHiggsino2016.cxx
:- Comment
addHistogram()
- Sort baseline and signal objects after we get them.
- Remove met cut.
- Add number of baseline electrons and baseline muons in the ntuple
- Change the initializatin value from 0 to -99999.
- Comment
v09. Add stack plots function and mll study python script
RunSimpleAnalysis.sh
: Add new samples for the truth studyEwkHiggsino2016.cxx
:- Add
baselineTaus
andnTaus
- Set
MTauTau
initial to -99999 - Comment met cut
- Move
MTauTau
calculation toif
condition
- Add
truth3_plots_comparison.py
:- Add more samples
- Add new function
compare_stack_curves
and make stack plots - Bugs fix
truth3_mll_study.py
: New python script
EwkHiggsino2016.cxx
v08. Bugs fixed in RunSimpleAnalysis.sh
: Add new samples for the truth studytruth3_plots_comparison.py
:- Add new samples
- Set log Y for some plots and modify the Y maximum of the plots
- Change the X range for mll plot
- Add cuts for some plots when project
TTree
to histogram
EwkHiggsino2016.cxx
:- Add some histograms which Jesse used in his code
- Copy
calc_massive_MT2
andcalc_MTauTau
from Jesse's code - Comment the radius for OR
- Add OR(jet, mu), and OR(el, mu) but comment them
- Add MET < 50 cuts and it only works on NUHM2 samples
- Fix the bug of
HTLep12
andMETOverHTLep12
calculations - Use Jesse's function to calculate
MTauTau
v07. Add codes for SimpleAnalysis
- Add new codes
EwkHiggsino2016.cxx
for SimpleAnalysis - Add new script
RunSimpleAnalysis.sh
for running SimpleAnalysis - Add new plotting python script
truth3_plots_comparison.py
- Add
spectrum.py
v06. New script for significance
AtlasStyle.py
: New scriptsignificance.py
: New scriptplots_making.py
:- Add version for the input file
- Add ATLASLabel
- Add plot cosmetics
get_yields.py
: change version of the input file to v1.5.gitignore
: ignores python.pyc
fileytEventSelection.cxx
:- Use emul triggers:
HLT_mu4_j125_xe90_mht_emul
andHLT_2mu4_j85_xe50_mht_emul
- Use emul triggers:
v05. Update to common ntuple v1.5
Run.sh
: Update to common ntuple v1.5Run.cxx
: Update to common ntuple v1.5ytEventSelection.h
: Comment some histograms related toMETOverHT
andMETOverHTLep12
ytEventSelection.cxx
:- Comment some histograms related to
METOverHT
andMETOverHTLep12
- Add FS protection to reject FS=86 in NUHM2 m12=600
- Remove
bTagWeight
from the weight calculation - Move cutflow related code into
if
condition - Modify same flavor, opposit sign, and number of baseline and signal lepton requirement
- Use all SR cuts including METOverHT cut
- Comment some histograms related to
get_yields.py
: Recoding this scriptplots_making.py
: Add some new histograms
v04. Finish implement the cutflow and add all signal samples
yt_cutflows.h
:- Add numbers for enumeration
- Remove argument of
print()
method - Remove
float weight
from sequential cut methods
ytEventSelection.h
:- Add
isMC
andsample
- Add
set_isMC()
andset_sample()
- Add
yt_cutflows.cxx
:- Add numbers for
cut_name
- Modify
print()
- Modify all
pass_xxx()
methods
- Add numbers for
ytEventSelection.cxx
:- Select sample depends on DSID for cutflow study
- Calculate weight for MC only
- Add sequential cuts
- print out cutflow results
Run.cxx
- Add
isMC
and related codes - Fix typo for
submitDir
for cutflow - Change the data
inputFilePath
- Add all signal samples for Higgsino and Slepton
- Add
- Add new script
Run.sh
v03. Add two new classes
yt_regions
:- Add this new class files:
yt_regions.h
andyt_regions.cxx
- Add this new class into
LinkDef.h
- Add this new class files:
yt_cutflows
:- Add this new class files:
yt_cutflows.h
andyt_cutflows.cxx
- Add this new class into
LinkDef.h
- Add this new class files:
ytEventSelection.h
:- Add
yt_regions
andyt_cutflows
objects and related codes - Add
isCutflow
,sample_type
, andluminosity
- Add histograms for the distributions and yields
- Add
set_isCutflow()
,set_sample_type()
, andset_luminosity()
methods - Add
debug_print()
method
- Add
ytEventSelection.cxx
:- Add
yt_regions
andyt_cutflows
objects - Add histograms for the distributions and yields
- Add weight calculation
- Add cutflows
- Implement
debug_print()
- Add
Run.cxx
:- Add descriptions for the usage
- Add
isCutflow
and the related things - Change the input files path and file name
- Specify the tree name
- Add the luminosity
- Add python script
plots_making.py
to make distribution plots setupHiggsino.sh
: Add path for the scripts- Add
.gitignore
- Add
get_yields.py
to show the yields and weighted yields - Add
list_trees.py
to list the name of trees in root files
v02. Add new directories and files
- LinkDef.h: Add
ytEventSelection
- Add new directories and files
- ytHiggsinoAnalysis/Root/ytEventSelection.cxx
- ytHiggsinoAnalysis/misc/NUHM2.C
- ytHiggsinoAnalysis/misc/NUHM2.h
- ytHiggsinoAnalysis/util/Run.cxx
- ytHiggsinoAnalysis/ytHiggsinoAnalysis/ytEventSelection.h