'cutfile_path': '../options/jesal_cutflow/DY_peak.txt',
            'lepton': 'tau',
            'label': r'$W^+\rightarrow\tau\nu\rightarrow\mu\nu$',
        }
    }

    analysis = Analysis(datasets,
                        'mutau_compare',
                        log_level=10,
                        log_out='both',
                        timedatelog=False,
                        year='2015+2016',
                        force_rebuild=False,
                        TTree_name='nominal_Loose')

    analysis.merge_datasets("wminmunu", "wminmunu_hm", verify=True)
    analysis.merge_datasets("wmintaunu", "wmintaunu_hm", verify=True)
    analysis.merge_datasets("wplusmunu", "wplusmunu_hm", verify=True)
    analysis.merge_datasets("wplustaunu", "wplustaunu_hm", verify=True)

    # normalised
    analysis.plot_hist(['wminmunu', 'wmintaunu'],
                       'met_met',
                       weight='reco_weight',
                       title='reco 139fb$^{-1}$',
                       bins=(30, 150, 5000),
                       logbins=True,
                       logx=True,
                       lepton='muon',
                       normalise='lumi',
                       yerr='rsumw2',
Exemplo n.º 2
0
    analysis = Analysis(
        datasets,
        'mutau_compare_full',
        data_dir='/data/keanu/framework_outputs/mutau_compare_full/',
        log_level=10,
        log_out='both',
        timedatelog=True,
        year='2015+2016',
        force_rebuild=False,
        TTree_name='truth',
        hard_cut='M_W')

    analysis.print_latex_table()

    analysis.merge_datasets("wminmunu", "wminmunu_hm")
    analysis.merge_datasets("wmintaunu", "wmintaunu_hm")
    analysis.merge_datasets("wplusmunu", "wplusmunu_hm")
    analysis.merge_datasets("wplustaunu", "wplustaunu_hm")

    # =========================
    # ===== TRUTH - UNCUT =====
    # =========================
    # normalised
    analysis.plot_hist(['wminmunu', 'wmintaunu'],
                       'MC_WZ_dilep_m_born',
                       weight='truth_weight',
                       title='truth - 139fb$^{-1}$',
                       bins=(30, 50, 5000),
                       logbins=True,
                       logx=True,