From 2af62e4979670556901490c3c6665fdd5eb447be Mon Sep 17 00:00:00 2001 From: rkansal47 Date: Sun, 15 Dec 2024 03:35:19 -0800 Subject: [PATCH] kt reclustering --- src/HHbbVV/postprocessing/TopAnalysis.ipynb | 1211 +---------------- src/HHbbVV/processors/corrections.py | 9 +- .../scale_factors/top_reweighting.ipynb | 82 +- .../submit_configs/ttsfs_24_12_15_lp.yaml | 26 + 4 files changed, 93 insertions(+), 1235 deletions(-) create mode 100644 src/condor/submit_configs/ttsfs_24_12_15_lp.yaml diff --git a/src/HHbbVV/postprocessing/TopAnalysis.ipynb b/src/HHbbVV/postprocessing/TopAnalysis.ipynb index c0169a85..89cbab91 100644 --- a/src/HHbbVV/postprocessing/TopAnalysis.ipynb +++ b/src/HHbbVV/postprocessing/TopAnalysis.ipynb @@ -367,36 +367,9 @@ }, { "cell_type": "code", - "execution_count": 18, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - " 0%| | 0/1 [00:00" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "plotting.plot_lund_plane_six(\n", " top_mc_sig_ratio,\n", @@ -480,20 +442,9 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/tmp/ipykernel_441066/1362833123.py:19: RuntimeWarning: divide by zero encountered in divide\n", - " top_mc_sig_ratio = np.nan_to_num(mc_density / sig_density, nan=1.0)\n", - "/tmp/ipykernel_441066/1362833123.py:19: RuntimeWarning: invalid value encountered in divide\n", - " top_mc_sig_ratio = np.nan_to_num(mc_density / sig_density, nan=1.0)\n" - ] - } - ], + "outputs": [], "source": [ "f = uproot.open(package_path / f\"corrections/lp_ratios/top_RW_2018_june17_ratio.root\")\n", "\n", @@ -518,1128 +469,6 @@ " top_mc_sig_ratio = np.clip(mc_sig_ratio, 0.05, 20.0)" ] }, - { - "cell_type": "code", - "execution_count": 13, - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([[[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]],\n", - "\n", - " [[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [2.09993494e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]],\n", - "\n", - " [[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [1.00000000e+000, 1.30434852e+000, 1.71632566e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [2.51013961e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.79769313e+308, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]],\n", - "\n", - " [[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [2.84667518e+000, 7.79620050e-004, 2.20348269e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.56798705e+000, 2.44714998e+000, 3.61465364e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.39415386e+001, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]],\n", - "\n", - " [[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [1.32434861e+000, 1.87885814e+000, 2.34885156e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.13724538e+000, 1.71413581e+000, 4.86292100e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [7.06610142e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]],\n", - "\n", - " [[1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [1.00000000e+000, 1.00000000e+000, 1.00000000e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " ...,\n", - " [1.00000000e+000, 2.01223679e+000, 1.30727249e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [6.63326847e-001, 3.97771103e+000, 1.56279968e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000],\n", - " [3.05064194e+000, 1.81662046e+000, 5.19524562e+000, ...,\n", - " 1.00000000e+000, 1.00000000e+000, 1.00000000e+000]]])" - ] - }, - "execution_count": 13, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "top_mc_sig_ratio" - ] - }, - { - "cell_type": "code", - "execution_count": 6, - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "['ratio_nom;3',\n", - " 'ratio_nom;2',\n", - " 'ratio_nom;1',\n", - " 'mc_nom;1',\n", - " 'bkg_nom;1',\n", - " 'data;1',\n", - " 'ratio_pdf_up;1',\n", - " 'ratio_pdf_down;1',\n", - " 'ratio_prefire_up;1',\n", - " 'ratio_prefire_down;1',\n", - " 'ratio_pileup_up;1',\n", - " 'ratio_pileup_down;1',\n", - " 'ratio_btag_up;1',\n", - " 'ratio_btag_down;1',\n", - " 'ratio_PS_ISR_up;1',\n", - " 'ratio_PS_ISR_down;1',\n", - " 'ratio_PS_FSR_up;1',\n", - " 'ratio_PS_FSR_down;1',\n", - " 'ratio_ttbar_F_up;1',\n", - " 'ratio_ttbar_F_down;1',\n", - " 'ratio_ttbar_R_up;1',\n", - " 'ratio_ttbar_R_down;1',\n", - " 'ratio_ttbar_RF_up;1',\n", - " 'ratio_ttbar_RF_down;1',\n", - " 'ratio_QCD_F_up;1',\n", - " 'ratio_QCD_F_down;1',\n", - " 'ratio_QCD_R_up;1',\n", - " 'ratio_QCD_R_down;1',\n", - " 'ratio_QCD_RF_up;1',\n", - " 'ratio_QCD_RF_down;1',\n", - " 'ratio_tW_F_up;1',\n", - " 'ratio_tW_F_down;1',\n", - " 'ratio_tW_R_up;1',\n", - " 'ratio_tW_R_down;1',\n", - " 'ratio_tW_RF_up;1',\n", - " 'ratio_tW_RF_down;1',\n", - " 'ratio_Single_F_up;1',\n", - " 'ratio_Single_F_down;1',\n", - " 'ratio_Single_R_up;1',\n", - " 'ratio_Single_R_down;1',\n", - " 'ratio_Single_RF_up;1',\n", - " 'ratio_Single_RF_down;1',\n", - " 'ratio_diboson_F_up;1',\n", - " 'ratio_diboson_F_down;1',\n", - " 'ratio_diboson_R_up;1',\n", - " 'ratio_diboson_R_down;1',\n", - " 'ratio_diboson_RF_up;1',\n", - " 'ratio_diboson_RF_down;1',\n", - " 'ratio_top_ptrw_up;1',\n", - " 'ratio_top_ptrw_down;1',\n", - " 'ratio_mu_trigger_up;1',\n", - " 'ratio_mu_trigger_down;1',\n", - " 'ratio_mu_id_up;1',\n", - " 'ratio_mu_id_down;1',\n", - " 'ratio_QCD_norm_up;1',\n", - " 'ratio_QCD_norm_down;1',\n", - " 'ratio_tW_norm_up;1',\n", - " 'ratio_tW_norm_down;1',\n", - " 'ratio_Diboson_norm_up;1',\n", - " 'ratio_Diboson_norm_down;1',\n", - " 'ratio_Single_norm_up;1',\n", - " 'ratio_Single_norm_down;1',\n", - " 'ratio_unmatched_norm_up;1',\n", - " 'ratio_unmatched_norm_down;1',\n", - " 'ratio_sys_tot_up;3',\n", - " 'ratio_sys_tot_up;2',\n", - " 'ratio_sys_tot_up;1',\n", - " 'ratio_sys_tot_down;3',\n", - " 'ratio_sys_tot_down;2',\n", - " 'ratio_sys_tot_down;1',\n", - " 'pt_extrap;1',\n", - " 'pt_extrap/func_covar_1_1;1',\n", - " 'pt_extrap/func_1_1;1',\n", - " 'pt_extrap/func_covar_1_2;1',\n", - " 'pt_extrap/func_1_2;1',\n", - " 'pt_extrap/func_covar_1_3;1',\n", - " 'pt_extrap/func_1_3;1',\n", - " 'pt_extrap/func_covar_1_4;1',\n", - " 'pt_extrap/func_1_4;1',\n", - " 'pt_extrap/func_covar_1_5;1',\n", - " 'pt_extrap/func_1_5;1',\n", - " 'pt_extrap/func_covar_1_6;1',\n", - " 'pt_extrap/func_1_6;1',\n", - " 'pt_extrap/func_covar_1_7;1',\n", - " 'pt_extrap/func_1_7;1',\n", - " 'pt_extrap/func_covar_1_8;1',\n", - " 'pt_extrap/func_1_8;1',\n", - " 'pt_extrap/func_covar_1_9;1',\n", - " 'pt_extrap/func_1_9;1',\n", - " 'pt_extrap/func_covar_1_10;1',\n", - " 'pt_extrap/func_1_10;1',\n", - " 'pt_extrap/func_covar_1_11;1',\n", - " 'pt_extrap/func_1_11;1',\n", - " 'pt_extrap/func_covar_1_12;1',\n", - " 'pt_extrap/func_1_12;1',\n", - " 'pt_extrap/func_covar_1_13;1',\n", - " 'pt_extrap/func_1_13;1',\n", - " 'pt_extrap/func_covar_1_14;1',\n", - " 'pt_extrap/func_1_14;1',\n", - " 'pt_extrap/func_covar_1_15;1',\n", - " 'pt_extrap/func_1_15;1',\n", - " 'pt_extrap/func_covar_1_16;1',\n", - " 'pt_extrap/func_1_16;1',\n", - " 'pt_extrap/func_covar_1_17;1',\n", - " 'pt_extrap/func_1_17;1',\n", - " 'pt_extrap/func_covar_1_18;1',\n", - " 'pt_extrap/func_1_18;1',\n", - " 'pt_extrap/func_covar_1_19;1',\n", - " 'pt_extrap/func_1_19;1',\n", - " 'pt_extrap/func_covar_1_20;1',\n", - " 'pt_extrap/func_1_20;1',\n", - " 'pt_extrap/func_covar_2_1;1',\n", - " 'pt_extrap/func_2_1;1',\n", - " 'pt_extrap/func_covar_2_2;1',\n", - " 'pt_extrap/func_2_2;1',\n", - " 'pt_extrap/func_covar_2_3;1',\n", - " 'pt_extrap/func_2_3;1',\n", - " 'pt_extrap/func_covar_2_4;1',\n", - " 'pt_extrap/func_2_4;1',\n", - " 'pt_extrap/func_covar_2_5;1',\n", - " 'pt_extrap/func_2_5;1',\n", - " 'pt_extrap/func_covar_2_6;1',\n", - " 'pt_extrap/func_2_6;1',\n", - " 'pt_extrap/func_covar_2_7;1',\n", - " 'pt_extrap/func_2_7;1',\n", - " 'pt_extrap/func_covar_2_8;1',\n", - " 'pt_extrap/func_2_8;1',\n", - " 'pt_extrap/func_covar_2_9;1',\n", - " 'pt_extrap/func_2_9;1',\n", - " 'pt_extrap/func_covar_2_10;1',\n", - " 'pt_extrap/func_2_10;1',\n", - " 'pt_extrap/func_covar_2_11;1',\n", - " 'pt_extrap/func_2_11;1',\n", - " 'pt_extrap/func_covar_2_12;1',\n", - " 'pt_extrap/func_2_12;1',\n", - " 'pt_extrap/func_covar_2_13;1',\n", - " 'pt_extrap/func_2_13;1',\n", - " 'pt_extrap/func_covar_2_14;1',\n", - " 'pt_extrap/func_2_14;1',\n", - " 'pt_extrap/func_covar_2_15;1',\n", - " 'pt_extrap/func_2_15;1',\n", - " 'pt_extrap/func_covar_2_16;1',\n", - " 'pt_extrap/func_2_16;1',\n", - " 'pt_extrap/func_covar_2_17;1',\n", - " 'pt_extrap/func_2_17;1',\n", - " 'pt_extrap/func_covar_2_18;1',\n", - " 'pt_extrap/func_2_18;1',\n", - " 'pt_extrap/func_covar_2_19;1',\n", - " 'pt_extrap/func_2_19;1',\n", - " 'pt_extrap/func_covar_2_20;1',\n", - " 'pt_extrap/func_2_20;1',\n", - " 'pt_extrap/func_covar_3_1;1',\n", - " 'pt_extrap/func_3_1;1',\n", - " 'pt_extrap/func_covar_3_2;1',\n", - " 'pt_extrap/func_3_2;1',\n", - " 'pt_extrap/func_covar_3_3;1',\n", - " 'pt_extrap/func_3_3;1',\n", - " 'pt_extrap/func_covar_3_4;1',\n", - " 'pt_extrap/func_3_4;1',\n", - " 'pt_extrap/func_covar_3_5;1',\n", - " 'pt_extrap/func_3_5;1',\n", - " 'pt_extrap/func_covar_3_6;1',\n", - " 'pt_extrap/func_3_6;1',\n", - " 'pt_extrap/func_covar_3_7;1',\n", - " 'pt_extrap/func_3_7;1',\n", - " 'pt_extrap/func_covar_3_8;1',\n", - " 'pt_extrap/func_3_8;1',\n", - " 'pt_extrap/func_covar_3_9;1',\n", - " 'pt_extrap/func_3_9;1',\n", - " 'pt_extrap/func_covar_3_10;1',\n", - " 'pt_extrap/func_3_10;1',\n", - " 'pt_extrap/func_covar_3_11;1',\n", - " 'pt_extrap/func_3_11;1',\n", - " 'pt_extrap/func_covar_3_12;1',\n", - " 'pt_extrap/func_3_12;1',\n", - " 'pt_extrap/func_covar_3_13;1',\n", - " 'pt_extrap/func_3_13;1',\n", - " 'pt_extrap/func_covar_3_14;1',\n", - " 'pt_extrap/func_3_14;1',\n", - " 'pt_extrap/func_covar_3_15;1',\n", - " 'pt_extrap/func_3_15;1',\n", - " 'pt_extrap/func_covar_3_16;1',\n", - " 'pt_extrap/func_3_16;1',\n", - " 'pt_extrap/func_covar_3_17;1',\n", - " 'pt_extrap/func_3_17;1',\n", - " 'pt_extrap/func_covar_3_18;1',\n", - " 'pt_extrap/func_3_18;1',\n", - " 'pt_extrap/func_covar_3_19;1',\n", - " 'pt_extrap/func_3_19;1',\n", - " 'pt_extrap/func_covar_3_20;1',\n", - " 'pt_extrap/func_3_20;1',\n", - " 'pt_extrap/func_covar_4_1;1',\n", - " 'pt_extrap/func_4_1;1',\n", - " 'pt_extrap/func_covar_4_2;1',\n", - " 'pt_extrap/func_4_2;1',\n", - " 'pt_extrap/func_covar_4_3;1',\n", - " 'pt_extrap/func_4_3;1',\n", - " 'pt_extrap/func_covar_4_4;1',\n", - " 'pt_extrap/func_4_4;1',\n", - " 'pt_extrap/func_covar_4_5;1',\n", - " 'pt_extrap/func_4_5;1',\n", - " 'pt_extrap/func_covar_4_6;1',\n", - " 'pt_extrap/func_4_6;1',\n", - " 'pt_extrap/func_covar_4_7;1',\n", - " 'pt_extrap/func_4_7;1',\n", - " 'pt_extrap/func_covar_4_8;1',\n", - " 'pt_extrap/func_4_8;1',\n", - " 'pt_extrap/func_covar_4_9;1',\n", - " 'pt_extrap/func_4_9;1',\n", - " 'pt_extrap/func_covar_4_10;1',\n", - " 'pt_extrap/func_4_10;1',\n", - " 'pt_extrap/func_covar_4_11;1',\n", - " 'pt_extrap/func_4_11;1',\n", - " 'pt_extrap/func_covar_4_12;1',\n", - " 'pt_extrap/func_4_12;1',\n", - " 'pt_extrap/func_covar_4_13;1',\n", - " 'pt_extrap/func_4_13;1',\n", - " 'pt_extrap/func_covar_4_14;1',\n", - " 'pt_extrap/func_4_14;1',\n", - " 'pt_extrap/func_covar_4_15;1',\n", - " 'pt_extrap/func_4_15;1',\n", - " 'pt_extrap/func_covar_4_16;1',\n", - " 'pt_extrap/func_4_16;1',\n", - " 'pt_extrap/func_covar_4_17;1',\n", - " 'pt_extrap/func_4_17;1',\n", - " 'pt_extrap/func_covar_4_18;1',\n", - " 'pt_extrap/func_4_18;1',\n", - " 'pt_extrap/func_covar_4_19;1',\n", - " 'pt_extrap/func_4_19;1',\n", - " 'pt_extrap/func_covar_4_20;1',\n", - " 'pt_extrap/func_4_20;1',\n", - " 'pt_extrap/func_covar_5_1;1',\n", - " 'pt_extrap/func_5_1;1',\n", - " 'pt_extrap/func_covar_5_2;1',\n", - " 'pt_extrap/func_5_2;1',\n", - " 'pt_extrap/func_covar_5_3;1',\n", - " 'pt_extrap/func_5_3;1',\n", - " 'pt_extrap/func_covar_5_4;1',\n", - " 'pt_extrap/func_5_4;1',\n", - " 'pt_extrap/func_covar_5_5;1',\n", - " 'pt_extrap/func_5_5;1',\n", - " 'pt_extrap/func_covar_5_6;1',\n", - " 'pt_extrap/func_5_6;1',\n", - " 'pt_extrap/func_covar_5_7;1',\n", - " 'pt_extrap/func_5_7;1',\n", - " 'pt_extrap/func_covar_5_8;1',\n", - " 'pt_extrap/func_5_8;1',\n", - " 'pt_extrap/func_covar_5_9;1',\n", - " 'pt_extrap/func_5_9;1',\n", - " 'pt_extrap/func_covar_5_10;1',\n", - " 'pt_extrap/func_5_10;1',\n", - " 'pt_extrap/func_covar_5_11;1',\n", - " 'pt_extrap/func_5_11;1',\n", - " 'pt_extrap/func_covar_5_12;1',\n", - " 'pt_extrap/func_5_12;1',\n", - " 'pt_extrap/func_covar_5_13;1',\n", - " 'pt_extrap/func_5_13;1',\n", - " 'pt_extrap/func_covar_5_14;1',\n", - " 'pt_extrap/func_5_14;1',\n", - " 'pt_extrap/func_covar_5_15;1',\n", - " 'pt_extrap/func_5_15;1',\n", - " 'pt_extrap/func_covar_5_16;1',\n", - " 'pt_extrap/func_5_16;1',\n", - " 'pt_extrap/func_covar_5_17;1',\n", - " 'pt_extrap/func_5_17;1',\n", - " 'pt_extrap/func_covar_5_18;1',\n", - " 'pt_extrap/func_5_18;1',\n", - " 'pt_extrap/func_covar_5_19;1',\n", - " 'pt_extrap/func_5_19;1',\n", - " 'pt_extrap/func_covar_5_20;1',\n", - " 'pt_extrap/func_5_20;1',\n", - " 'pt_extrap/func_covar_6_1;1',\n", - " 'pt_extrap/func_6_1;1',\n", - " 'pt_extrap/func_covar_6_2;1',\n", - " 'pt_extrap/func_6_2;1',\n", - " 'pt_extrap/func_covar_6_3;1',\n", - " 'pt_extrap/func_6_3;1',\n", - " 'pt_extrap/func_covar_6_4;1',\n", - " 'pt_extrap/func_6_4;1',\n", - " 'pt_extrap/func_covar_6_5;1',\n", - " 'pt_extrap/func_6_5;1',\n", - " 'pt_extrap/func_covar_6_6;1',\n", - " 'pt_extrap/func_6_6;1',\n", - " 'pt_extrap/func_covar_6_7;1',\n", - " 'pt_extrap/func_6_7;1',\n", - " 'pt_extrap/func_covar_6_8;1',\n", - " 'pt_extrap/func_6_8;1',\n", - " 'pt_extrap/func_covar_6_9;1',\n", - " 'pt_extrap/func_6_9;1',\n", - " 'pt_extrap/func_covar_6_10;1',\n", - " 'pt_extrap/func_6_10;1',\n", - " 'pt_extrap/func_covar_6_11;1',\n", - " 'pt_extrap/func_6_11;1',\n", - " 'pt_extrap/func_covar_6_12;1',\n", - " 'pt_extrap/func_6_12;1',\n", - " 'pt_extrap/func_covar_6_13;1',\n", - " 'pt_extrap/func_6_13;1',\n", - " 'pt_extrap/func_covar_6_14;1',\n", - " 'pt_extrap/func_6_14;1',\n", - " 'pt_extrap/func_covar_6_15;1',\n", - " 'pt_extrap/func_6_15;1',\n", - " 'pt_extrap/func_covar_6_16;1',\n", - " 'pt_extrap/func_6_16;1',\n", - " 'pt_extrap/func_covar_6_17;1',\n", - " 'pt_extrap/func_6_17;1',\n", - " 'pt_extrap/func_covar_6_18;1',\n", - " 'pt_extrap/func_6_18;1',\n", - " 'pt_extrap/func_covar_6_19;1',\n", - " 'pt_extrap/func_6_19;1',\n", - " 'pt_extrap/func_covar_6_20;1',\n", - " 'pt_extrap/func_6_20;1',\n", - " 'pt_extrap/func_covar_7_1;1',\n", - " 'pt_extrap/func_7_1;1',\n", - " 'pt_extrap/func_covar_7_2;1',\n", - " 'pt_extrap/func_7_2;1',\n", - " 'pt_extrap/func_covar_7_3;1',\n", - " 'pt_extrap/func_7_3;1',\n", - " 'pt_extrap/func_covar_7_4;1',\n", - " 'pt_extrap/func_7_4;1',\n", - " 'pt_extrap/func_covar_7_5;1',\n", - " 'pt_extrap/func_7_5;1',\n", - " 'pt_extrap/func_covar_7_6;1',\n", - " 'pt_extrap/func_7_6;1',\n", - " 'pt_extrap/func_covar_7_7;1',\n", - " 'pt_extrap/func_7_7;1',\n", - " 'pt_extrap/func_covar_7_8;1',\n", - " 'pt_extrap/func_7_8;1',\n", - " 'pt_extrap/func_covar_7_9;1',\n", - " 'pt_extrap/func_7_9;1',\n", - " 'pt_extrap/func_covar_7_10;1',\n", - " 'pt_extrap/func_7_10;1',\n", - " 'pt_extrap/func_covar_7_11;1',\n", - " 'pt_extrap/func_7_11;1',\n", - " 'pt_extrap/func_covar_7_12;1',\n", - " 'pt_extrap/func_7_12;1',\n", - " 'pt_extrap/func_covar_7_13;1',\n", - " 'pt_extrap/func_7_13;1',\n", - " 'pt_extrap/func_covar_7_14;1',\n", - " 'pt_extrap/func_7_14;1',\n", - " 'pt_extrap/func_covar_7_15;1',\n", - " 'pt_extrap/func_7_15;1',\n", - " 'pt_extrap/func_covar_7_16;1',\n", - " 'pt_extrap/func_7_16;1',\n", - " 'pt_extrap/func_covar_7_17;1',\n", - " 'pt_extrap/func_7_17;1',\n", - " 'pt_extrap/func_covar_7_18;1',\n", - " 'pt_extrap/func_7_18;1',\n", - " 'pt_extrap/func_covar_7_19;1',\n", - " 'pt_extrap/func_7_19;1',\n", - " 'pt_extrap/func_covar_7_20;1',\n", - " 'pt_extrap/func_7_20;1',\n", - " 'pt_extrap/func_covar_8_1;1',\n", - " 'pt_extrap/func_8_1;1',\n", - " 'pt_extrap/func_covar_8_2;1',\n", - " 'pt_extrap/func_8_2;1',\n", - " 'pt_extrap/func_covar_8_3;1',\n", - " 'pt_extrap/func_8_3;1',\n", - " 'pt_extrap/func_covar_8_4;1',\n", - " 'pt_extrap/func_8_4;1',\n", - " 'pt_extrap/func_covar_8_5;1',\n", - " 'pt_extrap/func_8_5;1',\n", - " 'pt_extrap/func_covar_8_6;1',\n", - " 'pt_extrap/func_8_6;1',\n", - " 'pt_extrap/func_covar_8_7;1',\n", - " 'pt_extrap/func_8_7;1',\n", - " 'pt_extrap/func_covar_8_8;1',\n", - " 'pt_extrap/func_8_8;1',\n", - " 'pt_extrap/func_covar_8_9;1',\n", - " 'pt_extrap/func_8_9;1',\n", - " 'pt_extrap/func_covar_8_10;1',\n", - " 'pt_extrap/func_8_10;1',\n", - " 'pt_extrap/func_covar_8_11;1',\n", - " 'pt_extrap/func_8_11;1',\n", - " 'pt_extrap/func_covar_8_12;1',\n", - " 'pt_extrap/func_8_12;1',\n", - " 'pt_extrap/func_covar_8_13;1',\n", - " 'pt_extrap/func_8_13;1',\n", - " 'pt_extrap/func_covar_8_14;1',\n", - " 'pt_extrap/func_8_14;1',\n", - " 'pt_extrap/func_covar_8_15;1',\n", - " 'pt_extrap/func_8_15;1',\n", - " 'pt_extrap/func_covar_8_16;1',\n", - " 'pt_extrap/func_8_16;1',\n", - " 'pt_extrap/func_covar_8_17;1',\n", - " 'pt_extrap/func_8_17;1',\n", - " 'pt_extrap/func_covar_8_18;1',\n", - " 'pt_extrap/func_8_18;1',\n", - " 'pt_extrap/func_covar_8_19;1',\n", - " 'pt_extrap/func_8_19;1',\n", - " 'pt_extrap/func_covar_8_20;1',\n", - " 'pt_extrap/func_8_20;1',\n", - " 'pt_extrap/func_covar_9_1;1',\n", - " 'pt_extrap/func_9_1;1',\n", - " 'pt_extrap/func_covar_9_2;1',\n", - " 'pt_extrap/func_9_2;1',\n", - " 'pt_extrap/func_covar_9_3;1',\n", - " 'pt_extrap/func_9_3;1',\n", - " 'pt_extrap/func_covar_9_4;1',\n", - " 'pt_extrap/func_9_4;1',\n", - " 'pt_extrap/func_covar_9_5;1',\n", - " 'pt_extrap/func_9_5;1',\n", - " 'pt_extrap/func_covar_9_6;1',\n", - " 'pt_extrap/func_9_6;1',\n", - " 'pt_extrap/func_covar_9_7;1',\n", - " 'pt_extrap/func_9_7;1',\n", - " 'pt_extrap/func_covar_9_8;1',\n", - " 'pt_extrap/func_9_8;1',\n", - " 'pt_extrap/func_covar_9_9;1',\n", - " 'pt_extrap/func_9_9;1',\n", - " 'pt_extrap/func_covar_9_10;1',\n", - " 'pt_extrap/func_9_10;1',\n", - " 'pt_extrap/func_covar_9_11;1',\n", - " 'pt_extrap/func_9_11;1',\n", - " 'pt_extrap/func_covar_9_12;1',\n", - " 'pt_extrap/func_9_12;1',\n", - " 'pt_extrap/func_covar_9_13;1',\n", - " 'pt_extrap/func_9_13;1',\n", - " 'pt_extrap/func_covar_9_14;1',\n", - " 'pt_extrap/func_9_14;1',\n", - " 'pt_extrap/func_covar_9_15;1',\n", - " 'pt_extrap/func_9_15;1',\n", - " 'pt_extrap/func_covar_9_16;1',\n", - " 'pt_extrap/func_9_16;1',\n", - " 'pt_extrap/func_covar_9_17;1',\n", - " 'pt_extrap/func_9_17;1',\n", - " 'pt_extrap/func_covar_9_18;1',\n", - " 'pt_extrap/func_9_18;1',\n", - " 'pt_extrap/func_covar_9_19;1',\n", - " 'pt_extrap/func_9_19;1',\n", - " 'pt_extrap/func_covar_9_20;1',\n", - " 'pt_extrap/func_9_20;1',\n", - " 'pt_extrap/func_covar_10_1;1',\n", - " 'pt_extrap/func_10_1;1',\n", - " 'pt_extrap/func_covar_10_2;1',\n", - " 'pt_extrap/func_10_2;1',\n", - " 'pt_extrap/func_covar_10_3;1',\n", - " 'pt_extrap/func_10_3;1',\n", - " 'pt_extrap/func_covar_10_4;1',\n", - " 'pt_extrap/func_10_4;1',\n", - " 'pt_extrap/func_covar_10_5;1',\n", - " 'pt_extrap/func_10_5;1',\n", - " 'pt_extrap/func_covar_10_6;1',\n", - " 'pt_extrap/func_10_6;1',\n", - " 'pt_extrap/func_covar_10_7;1',\n", - " 'pt_extrap/func_10_7;1',\n", - " 'pt_extrap/func_covar_10_8;1',\n", - " 'pt_extrap/func_10_8;1',\n", - " 'pt_extrap/func_covar_10_9;1',\n", - " 'pt_extrap/func_10_9;1',\n", - " 'pt_extrap/func_covar_10_10;1',\n", - " 'pt_extrap/func_10_10;1',\n", - " 'pt_extrap/func_covar_10_11;1',\n", - " 'pt_extrap/func_10_11;1',\n", - " 'pt_extrap/func_covar_10_12;1',\n", - " 'pt_extrap/func_10_12;1',\n", - " 'pt_extrap/func_covar_10_13;1',\n", - " 'pt_extrap/func_10_13;1',\n", - " 'pt_extrap/func_covar_10_14;1',\n", - " 'pt_extrap/func_10_14;1',\n", - " 'pt_extrap/func_covar_10_15;1',\n", - " 'pt_extrap/func_10_15;1',\n", - " 'pt_extrap/func_covar_10_16;1',\n", - " 'pt_extrap/func_10_16;1',\n", - " 'pt_extrap/func_covar_10_17;1',\n", - " 'pt_extrap/func_10_17;1',\n", - " 'pt_extrap/func_covar_10_18;1',\n", - " 'pt_extrap/func_10_18;1',\n", - " 'pt_extrap/func_covar_10_19;1',\n", - " 'pt_extrap/func_10_19;1',\n", - " 'pt_extrap/func_covar_10_20;1',\n", - " 'pt_extrap/func_10_20;1',\n", - " 'pt_extrap/func_covar_11_1;1',\n", - " 'pt_extrap/func_11_1;1',\n", - " 'pt_extrap/func_covar_11_2;1',\n", - " 'pt_extrap/func_11_2;1',\n", - " 'pt_extrap/func_covar_11_3;1',\n", - " 'pt_extrap/func_11_3;1',\n", - " 'pt_extrap/func_covar_11_4;1',\n", - " 'pt_extrap/func_11_4;1',\n", - " 'pt_extrap/func_covar_11_5;1',\n", - " 'pt_extrap/func_11_5;1',\n", - " 'pt_extrap/func_covar_11_6;1',\n", - " 'pt_extrap/func_11_6;1',\n", - " 'pt_extrap/func_covar_11_7;1',\n", - " 'pt_extrap/func_11_7;1',\n", - " 'pt_extrap/func_covar_11_8;1',\n", - " 'pt_extrap/func_11_8;1',\n", - " 'pt_extrap/func_covar_11_9;1',\n", - " 'pt_extrap/func_11_9;1',\n", - " 'pt_extrap/func_covar_11_10;1',\n", - " 'pt_extrap/func_11_10;1',\n", - " 'pt_extrap/func_covar_11_11;1',\n", - " 'pt_extrap/func_11_11;1',\n", - " 'pt_extrap/func_covar_11_12;1',\n", - " 'pt_extrap/func_11_12;1',\n", - " 'pt_extrap/func_covar_11_13;1',\n", - " 'pt_extrap/func_11_13;1',\n", - " 'pt_extrap/func_covar_11_14;1',\n", - " 'pt_extrap/func_11_14;1',\n", - " 'pt_extrap/func_covar_11_15;1',\n", - " 'pt_extrap/func_11_15;1',\n", - " 'pt_extrap/func_covar_11_16;1',\n", - " 'pt_extrap/func_11_16;1',\n", - " 'pt_extrap/func_covar_11_17;1',\n", - " 'pt_extrap/func_11_17;1',\n", - " 'pt_extrap/func_covar_11_18;1',\n", - " 'pt_extrap/func_11_18;1',\n", - " 'pt_extrap/func_covar_11_19;1',\n", - " 'pt_extrap/func_11_19;1',\n", - " 'pt_extrap/func_covar_11_20;1',\n", - " 'pt_extrap/func_11_20;1',\n", - " 'pt_extrap/func_covar_12_1;1',\n", - " 'pt_extrap/func_12_1;1',\n", - " 'pt_extrap/func_covar_12_2;1',\n", - " 'pt_extrap/func_12_2;1',\n", - " 'pt_extrap/func_covar_12_3;1',\n", - " 'pt_extrap/func_12_3;1',\n", - " 'pt_extrap/func_covar_12_4;1',\n", - " 'pt_extrap/func_12_4;1',\n", - " 'pt_extrap/func_covar_12_5;1',\n", - " 'pt_extrap/func_12_5;1',\n", - " 'pt_extrap/func_covar_12_6;1',\n", - " 'pt_extrap/func_12_6;1',\n", - " 'pt_extrap/func_covar_12_7;1',\n", - " 'pt_extrap/func_12_7;1',\n", - " 'pt_extrap/func_covar_12_8;1',\n", - " 'pt_extrap/func_12_8;1',\n", - " 'pt_extrap/func_covar_12_9;1',\n", - " 'pt_extrap/func_12_9;1',\n", - " 'pt_extrap/func_covar_12_10;1',\n", - " 'pt_extrap/func_12_10;1',\n", - " 'pt_extrap/func_covar_12_11;1',\n", - " 'pt_extrap/func_12_11;1',\n", - " 'pt_extrap/func_covar_12_12;1',\n", - " 'pt_extrap/func_12_12;1',\n", - " 'pt_extrap/func_covar_12_13;1',\n", - " 'pt_extrap/func_12_13;1',\n", - " 'pt_extrap/func_covar_12_14;1',\n", - " 'pt_extrap/func_12_14;1',\n", - " 'pt_extrap/func_covar_12_15;1',\n", - " 'pt_extrap/func_12_15;1',\n", - " 'pt_extrap/func_covar_12_16;1',\n", - " 'pt_extrap/func_12_16;1',\n", - " 'pt_extrap/func_covar_12_17;1',\n", - " 'pt_extrap/func_12_17;1',\n", - " 'pt_extrap/func_covar_12_18;1',\n", - " 'pt_extrap/func_12_18;1',\n", - " 'pt_extrap/func_covar_12_19;1',\n", - " 'pt_extrap/func_12_19;1',\n", - " 'pt_extrap/func_covar_12_20;1',\n", - " 'pt_extrap/func_12_20;1',\n", - " 'pt_extrap/func_covar_13_1;1',\n", - " 'pt_extrap/func_13_1;1',\n", - " 'pt_extrap/func_covar_13_2;1',\n", - " 'pt_extrap/func_13_2;1',\n", - " 'pt_extrap/func_covar_13_3;1',\n", - " 'pt_extrap/func_13_3;1',\n", - " 'pt_extrap/func_covar_13_4;1',\n", - " 'pt_extrap/func_13_4;1',\n", - " 'pt_extrap/func_covar_13_5;1',\n", - " 'pt_extrap/func_13_5;1',\n", - " 'pt_extrap/func_covar_13_6;1',\n", - " 'pt_extrap/func_13_6;1',\n", - " 'pt_extrap/func_covar_13_7;1',\n", - " 'pt_extrap/func_13_7;1',\n", - " 'pt_extrap/func_covar_13_8;1',\n", - " 'pt_extrap/func_13_8;1',\n", - " 'pt_extrap/func_covar_13_9;1',\n", - " 'pt_extrap/func_13_9;1',\n", - " 'pt_extrap/func_covar_13_10;1',\n", - " 'pt_extrap/func_13_10;1',\n", - " 'pt_extrap/func_covar_13_11;1',\n", - " 'pt_extrap/func_13_11;1',\n", - " 'pt_extrap/func_covar_13_12;1',\n", - " 'pt_extrap/func_13_12;1',\n", - " 'pt_extrap/func_covar_13_13;1',\n", - " 'pt_extrap/func_13_13;1',\n", - " 'pt_extrap/func_covar_13_14;1',\n", - " 'pt_extrap/func_13_14;1',\n", - " 'pt_extrap/func_covar_13_15;1',\n", - " 'pt_extrap/func_13_15;1',\n", - " 'pt_extrap/func_covar_13_16;1',\n", - " 'pt_extrap/func_13_16;1',\n", - " 'pt_extrap/func_covar_13_17;1',\n", - " 'pt_extrap/func_13_17;1',\n", - " 'pt_extrap/func_covar_13_18;1',\n", - " 'pt_extrap/func_13_18;1',\n", - " 'pt_extrap/func_covar_13_19;1',\n", - " 'pt_extrap/func_13_19;1',\n", - " 'pt_extrap/func_covar_13_20;1',\n", - " 'pt_extrap/func_13_20;1',\n", - " 'pt_extrap/func_covar_14_1;1',\n", - " 'pt_extrap/func_14_1;1',\n", - " 'pt_extrap/func_covar_14_2;1',\n", - " 'pt_extrap/func_14_2;1',\n", - " 'pt_extrap/func_covar_14_3;1',\n", - " 'pt_extrap/func_14_3;1',\n", - " 'pt_extrap/func_covar_14_4;1',\n", - " 'pt_extrap/func_14_4;1',\n", - " 'pt_extrap/func_covar_14_5;1',\n", - " 'pt_extrap/func_14_5;1',\n", - " 'pt_extrap/func_covar_14_6;1',\n", - " 'pt_extrap/func_14_6;1',\n", - " 'pt_extrap/func_covar_14_7;1',\n", - " 'pt_extrap/func_14_7;1',\n", - " 'pt_extrap/func_covar_14_8;1',\n", - " 'pt_extrap/func_14_8;1',\n", - " 'pt_extrap/func_covar_14_9;1',\n", - " 'pt_extrap/func_14_9;1',\n", - " 'pt_extrap/func_covar_14_10;1',\n", - " 'pt_extrap/func_14_10;1',\n", - " 'pt_extrap/func_covar_14_11;1',\n", - " 'pt_extrap/func_14_11;1',\n", - " 'pt_extrap/func_covar_14_12;1',\n", - " 'pt_extrap/func_14_12;1',\n", - " 'pt_extrap/func_covar_14_13;1',\n", - " 'pt_extrap/func_14_13;1',\n", - " 'pt_extrap/func_covar_14_14;1',\n", - " 'pt_extrap/func_14_14;1',\n", - " 'pt_extrap/func_covar_14_15;1',\n", - " 'pt_extrap/func_14_15;1',\n", - " 'pt_extrap/func_covar_14_16;1',\n", - " 'pt_extrap/func_14_16;1',\n", - " 'pt_extrap/func_covar_14_17;1',\n", - " 'pt_extrap/func_14_17;1',\n", - " 'pt_extrap/func_covar_14_18;1',\n", - " 'pt_extrap/func_14_18;1',\n", - " 'pt_extrap/func_covar_14_19;1',\n", - " 'pt_extrap/func_14_19;1',\n", - " 'pt_extrap/func_covar_14_20;1',\n", - " 'pt_extrap/func_14_20;1',\n", - " 'pt_extrap/func_covar_15_1;1',\n", - " 'pt_extrap/func_15_1;1',\n", - " 'pt_extrap/func_covar_15_2;1',\n", - " 'pt_extrap/func_15_2;1',\n", - " 'pt_extrap/func_covar_15_3;1',\n", - " 'pt_extrap/func_15_3;1',\n", - " 'pt_extrap/func_covar_15_4;1',\n", - " 'pt_extrap/func_15_4;1',\n", - " 'pt_extrap/func_covar_15_5;1',\n", - " 'pt_extrap/func_15_5;1',\n", - " 'pt_extrap/func_covar_15_6;1',\n", - " 'pt_extrap/func_15_6;1',\n", - " 'pt_extrap/func_covar_15_7;1',\n", - " 'pt_extrap/func_15_7;1',\n", - " 'pt_extrap/func_covar_15_8;1',\n", - " 'pt_extrap/func_15_8;1',\n", - " 'pt_extrap/func_covar_15_9;1',\n", - " 'pt_extrap/func_15_9;1',\n", - " 'pt_extrap/func_covar_15_10;1',\n", - " 'pt_extrap/func_15_10;1',\n", - " 'pt_extrap/func_covar_15_11;1',\n", - " 'pt_extrap/func_15_11;1',\n", - " 'pt_extrap/func_covar_15_12;1',\n", - " 'pt_extrap/func_15_12;1',\n", - " 'pt_extrap/func_covar_15_13;1',\n", - " 'pt_extrap/func_15_13;1',\n", - " 'pt_extrap/func_covar_15_14;1',\n", - " 'pt_extrap/func_15_14;1',\n", - " 'pt_extrap/func_covar_15_15;1',\n", - " 'pt_extrap/func_15_15;1',\n", - " 'pt_extrap/func_covar_15_16;1',\n", - " 'pt_extrap/func_15_16;1',\n", - " 'pt_extrap/func_covar_15_17;1',\n", - " 'pt_extrap/func_15_17;1',\n", - " 'pt_extrap/func_covar_15_18;1',\n", - " 'pt_extrap/func_15_18;1',\n", - " 'pt_extrap/func_covar_15_19;1',\n", - " 'pt_extrap/func_15_19;1',\n", - " 'pt_extrap/func_covar_15_20;1',\n", - " 'pt_extrap/func_15_20;1',\n", - " 'pt_extrap/func_covar_16_1;1',\n", - " 'pt_extrap/func_16_1;1',\n", - " 'pt_extrap/func_covar_16_2;1',\n", - " 'pt_extrap/func_16_2;1',\n", - " 'pt_extrap/func_covar_16_3;1',\n", - " 'pt_extrap/func_16_3;1',\n", - " 'pt_extrap/func_covar_16_4;1',\n", - " 'pt_extrap/func_16_4;1',\n", - " 'pt_extrap/func_covar_16_5;1',\n", - " 'pt_extrap/func_16_5;1',\n", - " 'pt_extrap/func_covar_16_6;1',\n", - " 'pt_extrap/func_16_6;1',\n", - " 'pt_extrap/func_covar_16_7;1',\n", - " 'pt_extrap/func_16_7;1',\n", - " 'pt_extrap/func_covar_16_8;1',\n", - " 'pt_extrap/func_16_8;1',\n", - " 'pt_extrap/func_covar_16_9;1',\n", - " 'pt_extrap/func_16_9;1',\n", - " 'pt_extrap/func_covar_16_10;1',\n", - " 'pt_extrap/func_16_10;1',\n", - " 'pt_extrap/func_covar_16_11;1',\n", - " 'pt_extrap/func_16_11;1',\n", - " 'pt_extrap/func_covar_16_12;1',\n", - " 'pt_extrap/func_16_12;1',\n", - " 'pt_extrap/func_covar_16_13;1',\n", - " 'pt_extrap/func_16_13;1',\n", - " 'pt_extrap/func_covar_16_14;1',\n", - " 'pt_extrap/func_16_14;1',\n", - " 'pt_extrap/func_covar_16_15;1',\n", - " 'pt_extrap/func_16_15;1',\n", - " 'pt_extrap/func_covar_16_16;1',\n", - " 'pt_extrap/func_16_16;1',\n", - " 'pt_extrap/func_covar_16_17;1',\n", - " 'pt_extrap/func_16_17;1',\n", - " 'pt_extrap/func_covar_16_18;1',\n", - " 'pt_extrap/func_16_18;1',\n", - " 'pt_extrap/func_covar_16_19;1',\n", - " 'pt_extrap/func_16_19;1',\n", - " 'pt_extrap/func_covar_16_20;1',\n", - " 'pt_extrap/func_16_20;1',\n", - " 'pt_extrap/func_covar_17_1;1',\n", - " 'pt_extrap/func_17_1;1',\n", - " 'pt_extrap/func_covar_17_2;1',\n", - " 'pt_extrap/func_17_2;1',\n", - " 'pt_extrap/func_covar_17_3;1',\n", - " 'pt_extrap/func_17_3;1',\n", - " 'pt_extrap/func_covar_17_4;1',\n", - " 'pt_extrap/func_17_4;1',\n", - " 'pt_extrap/func_covar_17_5;1',\n", - " 'pt_extrap/func_17_5;1',\n", - " 'pt_extrap/func_covar_17_6;1',\n", - " 'pt_extrap/func_17_6;1',\n", - " 'pt_extrap/func_covar_17_7;1',\n", - " 'pt_extrap/func_17_7;1',\n", - " 'pt_extrap/func_covar_17_8;1',\n", - " 'pt_extrap/func_17_8;1',\n", - " 'pt_extrap/func_covar_17_9;1',\n", - " 'pt_extrap/func_17_9;1',\n", - " 'pt_extrap/func_covar_17_10;1',\n", - " 'pt_extrap/func_17_10;1',\n", - " 'pt_extrap/func_covar_17_11;1',\n", - " 'pt_extrap/func_17_11;1',\n", - " 'pt_extrap/func_covar_17_12;1',\n", - " 'pt_extrap/func_17_12;1',\n", - " 'pt_extrap/func_covar_17_13;1',\n", - " 'pt_extrap/func_17_13;1',\n", - " 'pt_extrap/func_covar_17_14;1',\n", - " 'pt_extrap/func_17_14;1',\n", - " 'pt_extrap/func_covar_17_15;1',\n", - " 'pt_extrap/func_17_15;1',\n", - " 'pt_extrap/func_covar_17_16;1',\n", - " 'pt_extrap/func_17_16;1',\n", - " 'pt_extrap/func_covar_17_17;1',\n", - " 'pt_extrap/func_17_17;1',\n", - " 'pt_extrap/func_covar_17_18;1',\n", - " 'pt_extrap/func_17_18;1',\n", - " 'pt_extrap/func_covar_17_19;1',\n", - " 'pt_extrap/func_17_19;1',\n", - " 'pt_extrap/func_covar_17_20;1',\n", - " 'pt_extrap/func_17_20;1',\n", - " 'pt_extrap/func_covar_18_1;1',\n", - " 'pt_extrap/func_18_1;1',\n", - " 'pt_extrap/func_covar_18_2;1',\n", - " 'pt_extrap/func_18_2;1',\n", - " 'pt_extrap/func_covar_18_3;1',\n", - " 'pt_extrap/func_18_3;1',\n", - " 'pt_extrap/func_covar_18_4;1',\n", - " 'pt_extrap/func_18_4;1',\n", - " 'pt_extrap/func_covar_18_5;1',\n", - " 'pt_extrap/func_18_5;1',\n", - " 'pt_extrap/func_covar_18_6;1',\n", - " 'pt_extrap/func_18_6;1',\n", - " 'pt_extrap/func_covar_18_7;1',\n", - " 'pt_extrap/func_18_7;1',\n", - " 'pt_extrap/func_covar_18_8;1',\n", - " 'pt_extrap/func_18_8;1',\n", - " 'pt_extrap/func_covar_18_9;1',\n", - " 'pt_extrap/func_18_9;1',\n", - " 'pt_extrap/func_covar_18_10;1',\n", - " 'pt_extrap/func_18_10;1',\n", - " 'pt_extrap/func_covar_18_11;1',\n", - " 'pt_extrap/func_18_11;1',\n", - " 'pt_extrap/func_covar_18_12;1',\n", - " 'pt_extrap/func_18_12;1',\n", - " 'pt_extrap/func_covar_18_13;1',\n", - " 'pt_extrap/func_18_13;1',\n", - " 'pt_extrap/func_covar_18_14;1',\n", - " 'pt_extrap/func_18_14;1',\n", - " 'pt_extrap/func_covar_18_15;1',\n", - " 'pt_extrap/func_18_15;1',\n", - " 'pt_extrap/func_covar_18_16;1',\n", - " 'pt_extrap/func_18_16;1',\n", - " 'pt_extrap/func_covar_18_17;1',\n", - " 'pt_extrap/func_18_17;1',\n", - " 'pt_extrap/func_covar_18_18;1',\n", - " 'pt_extrap/func_18_18;1',\n", - " 'pt_extrap/func_covar_18_19;1',\n", - " 'pt_extrap/func_18_19;1',\n", - " 'pt_extrap/func_covar_18_20;1',\n", - " 'pt_extrap/func_18_20;1',\n", - " 'pt_extrap/func_covar_19_1;1',\n", - " 'pt_extrap/func_19_1;1',\n", - " 'pt_extrap/func_covar_19_2;1',\n", - " 'pt_extrap/func_19_2;1',\n", - " 'pt_extrap/func_covar_19_3;1',\n", - " 'pt_extrap/func_19_3;1',\n", - " 'pt_extrap/func_covar_19_4;1',\n", - " 'pt_extrap/func_19_4;1',\n", - " 'pt_extrap/func_covar_19_5;1',\n", - " 'pt_extrap/func_19_5;1',\n", - " 'pt_extrap/func_covar_19_6;1',\n", - " 'pt_extrap/func_19_6;1',\n", - " 'pt_extrap/func_covar_19_7;1',\n", - " 'pt_extrap/func_19_7;1',\n", - " 'pt_extrap/func_covar_19_8;1',\n", - " 'pt_extrap/func_19_8;1',\n", - " 'pt_extrap/func_covar_19_9;1',\n", - " 'pt_extrap/func_19_9;1',\n", - " 'pt_extrap/func_covar_19_10;1',\n", - " 'pt_extrap/func_19_10;1',\n", - " 'pt_extrap/func_covar_19_11;1',\n", - " 'pt_extrap/func_19_11;1',\n", - " 'pt_extrap/func_covar_19_12;1',\n", - " 'pt_extrap/func_19_12;1',\n", - " 'pt_extrap/func_covar_19_13;1',\n", - " 'pt_extrap/func_19_13;1',\n", - " 'pt_extrap/func_covar_19_14;1',\n", - " 'pt_extrap/func_19_14;1',\n", - " 'pt_extrap/func_covar_19_15;1',\n", - " 'pt_extrap/func_19_15;1',\n", - " 'pt_extrap/func_covar_19_16;1',\n", - " 'pt_extrap/func_19_16;1',\n", - " 'pt_extrap/func_covar_19_17;1',\n", - " 'pt_extrap/func_19_17;1',\n", - " 'pt_extrap/func_covar_19_18;1',\n", - " 'pt_extrap/func_19_18;1',\n", - " 'pt_extrap/func_covar_19_19;1',\n", - " 'pt_extrap/func_19_19;1',\n", - " 'pt_extrap/func_covar_19_20;1',\n", - " 'pt_extrap/func_19_20;1',\n", - " 'pt_extrap/func_covar_20_1;1',\n", - " 'pt_extrap/func_20_1;1',\n", - " 'pt_extrap/func_covar_20_2;1',\n", - " 'pt_extrap/func_20_2;1',\n", - " 'pt_extrap/func_covar_20_3;1',\n", - " 'pt_extrap/func_20_3;1',\n", - " 'pt_extrap/func_covar_20_4;1',\n", - " 'pt_extrap/func_20_4;1',\n", - " 'pt_extrap/func_covar_20_5;1',\n", - " 'pt_extrap/func_20_5;1',\n", - " 'pt_extrap/func_covar_20_6;1',\n", - " 'pt_extrap/func_20_6;1',\n", - " 'pt_extrap/func_covar_20_7;1',\n", - " 'pt_extrap/func_20_7;1',\n", - " 'pt_extrap/func_covar_20_8;1',\n", - " 'pt_extrap/func_20_8;1',\n", - " 'pt_extrap/func_covar_20_9;1',\n", - " 'pt_extrap/func_20_9;1',\n", - " 'pt_extrap/func_covar_20_10;1',\n", - " 'pt_extrap/func_20_10;1',\n", - " 'pt_extrap/func_covar_20_11;1',\n", - " 'pt_extrap/func_20_11;1',\n", - " 'pt_extrap/func_covar_20_12;1',\n", - " 'pt_extrap/func_20_12;1',\n", - " 'pt_extrap/func_covar_20_13;1',\n", - " 'pt_extrap/func_20_13;1',\n", - " 'pt_extrap/func_covar_20_14;1',\n", - " 'pt_extrap/func_20_14;1',\n", - " 'pt_extrap/func_covar_20_15;1',\n", - " 'pt_extrap/func_20_15;1',\n", - " 'pt_extrap/func_covar_20_16;1',\n", - " 'pt_extrap/func_20_16;1',\n", - " 'pt_extrap/func_covar_20_17;1',\n", - " 'pt_extrap/func_20_17;1',\n", - " 'pt_extrap/func_covar_20_18;1',\n", - " 'pt_extrap/func_20_18;1',\n", - " 'pt_extrap/func_covar_20_19;1',\n", - " 'pt_extrap/func_20_19;1',\n", - " 'pt_extrap/func_covar_20_20;1',\n", - " 'pt_extrap/func_20_20;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_1;1',\n", - " 'pt_extrap/func_sys_tot_up_1_1;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_2;1',\n", - " 'pt_extrap/func_sys_tot_up_1_2;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_3;1',\n", - " 'pt_extrap/func_sys_tot_up_1_3;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_4;1',\n", - " 'pt_extrap/func_sys_tot_up_1_4;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_5;1',\n", - " 'pt_extrap/func_sys_tot_up_1_5;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_6;1',\n", - " 'pt_extrap/func_sys_tot_up_1_6;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_7;1',\n", - " 'pt_extrap/func_sys_tot_up_1_7;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_8;1',\n", - " 'pt_extrap/func_sys_tot_up_1_8;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_9;1',\n", - " 'pt_extrap/func_sys_tot_up_1_9;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_10;1',\n", - " 'pt_extrap/func_sys_tot_up_1_10;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_11;1',\n", - " 'pt_extrap/func_sys_tot_up_1_11;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_12;1',\n", - " 'pt_extrap/func_sys_tot_up_1_12;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_13;1',\n", - " 'pt_extrap/func_sys_tot_up_1_13;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_14;1',\n", - " 'pt_extrap/func_sys_tot_up_1_14;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_15;1',\n", - " 'pt_extrap/func_sys_tot_up_1_15;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_16;1',\n", - " 'pt_extrap/func_sys_tot_up_1_16;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_17;1',\n", - " 'pt_extrap/func_sys_tot_up_1_17;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_18;1',\n", - " 'pt_extrap/func_sys_tot_up_1_18;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_19;1',\n", - " 'pt_extrap/func_sys_tot_up_1_19;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_1_20;1',\n", - " 'pt_extrap/func_sys_tot_up_1_20;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_1;1',\n", - " 'pt_extrap/func_sys_tot_up_2_1;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_2;1',\n", - " 'pt_extrap/func_sys_tot_up_2_2;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_3;1',\n", - " 'pt_extrap/func_sys_tot_up_2_3;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_4;1',\n", - " 'pt_extrap/func_sys_tot_up_2_4;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_5;1',\n", - " 'pt_extrap/func_sys_tot_up_2_5;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_6;1',\n", - " 'pt_extrap/func_sys_tot_up_2_6;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_7;1',\n", - " 'pt_extrap/func_sys_tot_up_2_7;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_8;1',\n", - " 'pt_extrap/func_sys_tot_up_2_8;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_9;1',\n", - " 'pt_extrap/func_sys_tot_up_2_9;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_10;1',\n", - " 'pt_extrap/func_sys_tot_up_2_10;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_11;1',\n", - " 'pt_extrap/func_sys_tot_up_2_11;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_12;1',\n", - " 'pt_extrap/func_sys_tot_up_2_12;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_13;1',\n", - " 'pt_extrap/func_sys_tot_up_2_13;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_14;1',\n", - " 'pt_extrap/func_sys_tot_up_2_14;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_15;1',\n", - " 'pt_extrap/func_sys_tot_up_2_15;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_16;1',\n", - " 'pt_extrap/func_sys_tot_up_2_16;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_17;1',\n", - " 'pt_extrap/func_sys_tot_up_2_17;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_18;1',\n", - " 'pt_extrap/func_sys_tot_up_2_18;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_19;1',\n", - " 'pt_extrap/func_sys_tot_up_2_19;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_2_20;1',\n", - " 'pt_extrap/func_sys_tot_up_2_20;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_1;1',\n", - " 'pt_extrap/func_sys_tot_up_3_1;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_2;1',\n", - " 'pt_extrap/func_sys_tot_up_3_2;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_3;1',\n", - " 'pt_extrap/func_sys_tot_up_3_3;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_4;1',\n", - " 'pt_extrap/func_sys_tot_up_3_4;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_5;1',\n", - " 'pt_extrap/func_sys_tot_up_3_5;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_6;1',\n", - " 'pt_extrap/func_sys_tot_up_3_6;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_7;1',\n", - " 'pt_extrap/func_sys_tot_up_3_7;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_8;1',\n", - " 'pt_extrap/func_sys_tot_up_3_8;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_9;1',\n", - " 'pt_extrap/func_sys_tot_up_3_9;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_10;1',\n", - " 'pt_extrap/func_sys_tot_up_3_10;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_11;1',\n", - " 'pt_extrap/func_sys_tot_up_3_11;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_12;1',\n", - " 'pt_extrap/func_sys_tot_up_3_12;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_13;1',\n", - " 'pt_extrap/func_sys_tot_up_3_13;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_14;1',\n", - " 'pt_extrap/func_sys_tot_up_3_14;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_15;1',\n", - " 'pt_extrap/func_sys_tot_up_3_15;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_16;1',\n", - " 'pt_extrap/func_sys_tot_up_3_16;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_17;1',\n", - " 'pt_extrap/func_sys_tot_up_3_17;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_18;1',\n", - " 'pt_extrap/func_sys_tot_up_3_18;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_19;1',\n", - " 'pt_extrap/func_sys_tot_up_3_19;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_3_20;1',\n", - " 'pt_extrap/func_sys_tot_up_3_20;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_4_1;1',\n", - " 'pt_extrap/func_sys_tot_up_4_1;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_4_2;1',\n", - " 'pt_extrap/func_sys_tot_up_4_2;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_4_3;1',\n", - " 'pt_extrap/func_sys_tot_up_4_3;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_4_4;1',\n", - " 'pt_extrap/func_sys_tot_up_4_4;1',\n", - " 'pt_extrap/func_sys_tot_up_covar_4_5;1',\n", - " ...]" - ] - }, - "execution_count": 6, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "f.keys()" - ] - }, { "cell_type": "code", "execution_count": null, @@ -1665,15 +494,6 @@ ")" ] }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [ - "h.variances()" - ] - }, { "cell_type": "code", "execution_count": null, @@ -1687,16 +507,6 @@ ")" ] }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [ - "print(np.mean(mc_sig_ratio_pt, axis=(1, 2)))\n", - "print(np.mean(mc_sig_old_ratio_pt, axis=(1, 2)))" - ] - }, { "cell_type": "code", "execution_count": null, @@ -2383,7 +1193,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3.10.4 ('python310')", + "display_name": "python39", "language": "python", "name": "python3" }, @@ -2399,12 +1209,7 @@ "pygments_lexer": "ipython3", "version": "3.9.15" }, - "orig_nbformat": 4, - "vscode": { - "interpreter": { - "hash": "31a7b1cb5f073f7a7d37b3db504c6954ce2b88e0f82e412b65ad0b5f2dd17394" - } - } + "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 diff --git a/src/HHbbVV/processors/corrections.py b/src/HHbbVV/processors/corrections.py index 99da0b5c..e036a46f 100644 --- a/src/HHbbVV/processors/corrections.py +++ b/src/HHbbVV/processors/corrections.py @@ -844,6 +844,7 @@ def _get_lund_arrays( fatjet_idx: int | ak.Array, num_prongs: int, min_pt: float = 1.0, + ca_recluster: bool = False, ): """ Gets the ``num_prongs`` subjet pTs and Delta and kT per primary LP splitting of fatjets at @@ -854,6 +855,8 @@ def _get_lund_arrays( jec_fatjets (FatJetArray): post-JEC fatjets, used to update subjet pTs. fatjet_idx (int | ak.Array): fatjet index num_prongs (int): number of prongs / subjets per jet to reweight + min_pt (float): minimum pT for pf candidates considered for LP declustering + ca_recluster (bool): whether to recluster subjets with CA after kT Returns: lds, kt_subjets_vec, kt_subjets_pt: lund declusterings, subjet 4-vectors, JEC-corrected subjet pTs @@ -870,6 +873,8 @@ def _get_lund_arrays( cadef = fastjet.JetDefinition(fastjet.cambridge_algorithm, dR) ktdef = fastjet.JetDefinition(fastjet.kt_algorithm, dR) + recluster_def = cadef if ca_recluster else ktdef + # get pfcands of the top-matched jets ak8_pfcands = events.FatJetPFCands ak8_pfcands = ak8_pfcands[ak8_pfcands.jetIdx == fatjet_idx] @@ -911,8 +916,8 @@ def _get_lund_arrays( # then re-cluster with CA # won't need to flatten once https://github.com/scikit-hep/fastjet/pull/145 is released - ca_clustering = fastjet.ClusterSequence(kt_subjet_consts, cadef) - lds = ca_clustering.exclusive_jets_lund_declusterings(1) + reclustering = fastjet.ClusterSequence(kt_subjet_consts, recluster_def) + lds = reclustering.exclusive_jets_lund_declusterings(1) return lds, kt_subjets_vec, kt_subjets_pt diff --git a/src/HHbbVV/scale_factors/top_reweighting.ipynb b/src/HHbbVV/scale_factors/top_reweighting.ipynb index 2be8baa0..7036142e 100644 --- a/src/HHbbVV/scale_factors/top_reweighting.ipynb +++ b/src/HHbbVV/scale_factors/top_reweighting.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -41,7 +41,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": {}, "outputs": [], "source": [ @@ -51,7 +51,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, "outputs": [], "source": [ @@ -132,7 +132,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, "outputs": [], "source": [ @@ -150,7 +150,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, "outputs": [], "source": [ @@ -196,7 +196,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, "outputs": [], "source": [ @@ -228,7 +228,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, "outputs": [], "source": [ @@ -252,7 +252,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], "source": [ @@ -498,7 +498,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, "outputs": [], "source": [ @@ -523,7 +523,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, "outputs": [], "source": [ @@ -538,7 +538,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, "outputs": [], "source": [ @@ -561,7 +561,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, "outputs": [], "source": [ @@ -604,7 +604,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, "outputs": [], "source": [ @@ -649,7 +649,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 18, "metadata": {}, "outputs": [], "source": [ @@ -658,7 +658,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "metadata": {}, "outputs": [], "source": [ @@ -677,7 +677,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 20, "metadata": {}, "outputs": [], "source": [ @@ -692,7 +692,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 21, "metadata": {}, "outputs": [], "source": [ @@ -725,7 +725,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 23, "metadata": {}, "outputs": [], "source": [ @@ -741,7 +741,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 24, "metadata": {}, "outputs": [], "source": [ @@ -763,7 +763,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 25, "metadata": {}, "outputs": [], "source": [ @@ -790,7 +790,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 26, "metadata": {}, "outputs": [], "source": [ @@ -805,7 +805,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 27, "metadata": {}, "outputs": [], "source": [ @@ -836,11 +836,38 @@ "# then re-cluster with CA\n", "# won't need to flatten once https://github.com/scikit-hep/fastjet/pull/145 is released\n", "# ca_clustering = fastjet.ClusterSequence(ak.flatten(kt_subjet_consts, axis=1), cadef)\n", - "ca_clustering = fastjet.ClusterSequence(filled_consts, cadef)\n", + "ca_clustering = fastjet.ClusterSequence(filled_consts, ktdef)\n", "lds = ca_clustering.exclusive_jets_lund_declusterings(1)\n", "lds_flat = ak.flatten(lds, axis=1)" ] }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "lds" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "merged_top_events" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "kt_subjet_consts" + ] + }, { "cell_type": "code", "execution_count": null, @@ -1432,7 +1459,7 @@ ], "metadata": { "kernelspec": { - "display_name": "python310", + "display_name": "python39", "language": "python", "name": "python3" }, @@ -1448,12 +1475,7 @@ "pygments_lexer": "ipython3", "version": "3.9.15" }, - "orig_nbformat": 4, - "vscode": { - "interpreter": { - "hash": "15adc7883e707560d0d9727709639b8fe3f3cff1f197d2d643742923ff23a29c" - } - } + "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 diff --git a/src/condor/submit_configs/ttsfs_24_12_15_lp.yaml b/src/condor/submit_configs/ttsfs_24_12_15_lp.yaml new file mode 100644 index 00000000..26b1143b --- /dev/null +++ b/src/condor/submit_configs/ttsfs_24_12_15_lp.yaml @@ -0,0 +1,26 @@ +{ + "TTbarScaleFactors": + { + "TTbar": + { + "files_per_job": 10, + "max_files": 80, + "chunksize": 80000, + "subsamples": ["TTToSemiLeptonic"], + }, + "SingleTop": + { + "files_per_job": 10, + "max_files": 80, + "chunksize": 80000, + "subsamples": + [ + "ST_tW_antitop_5f_NoFullyHadronicDecays", + "ST_tW_top_5f_NoFullyHadronicDecays", + "ST_s-channel_4f_leptonDecays", + "ST_t-channel_antitop_4f_InclusiveDecays", + "ST_t-channel_top_4f_InclusiveDecays", + ], + }, + }, +}