|
| | rng = np.random.default_rng() |
| |
| list | correction_labels = ["metphicorr_pfmet_mc", "metphicorr_puppimet_mc", "metphicorr_pfmet_data", "metphicorr_puppimet_data"] |
| |
| list | eras = ["2018_UL", "2017_UL", "2016postVFP_UL", "2016preVFP_UL"] |
| |
| list | run_ranges = [[315252, 325274], [297020, 306463], [278769, 284045], [272007, 278771]] |
| |
| string | infile = "met.json.gz" |
| |
| | is_data = None |
| |
| | ceval = correctionlib.CorrectionSet.from_file("../POG/JME/{}/{}".format(era,infile)) |
| |
| tuple | pts = (rng.pareto(1.5,size=1000000))*100 |
| |
| | phis = rng.uniform(low=-3.14,high=3.14,size=1000000) |
| |
| | npvs = rng.integers(low=0,high=200,size=1000000) |
| |
| | runs = None |
| |
| | corrected_pts = ceval["pt_{}".format(correction_label)].evaluate(pts,phis,npvs,runs) |
| |
| | corrected_phis = ceval["phi_{}".format(correction_label)].evaluate(pts,phis,npvs,runs) |
| |
| | fig |
| |
| | axs |
| |
| | sharey |
| |
| | True |
| |
| | tight_layout |
| |
| | xlabel |
| |
| | bins |
| |
| | range |
| |
| | ylabel |
| |