示例#1
0
            'TTree_name': 'truth',
            'year': '2015+2016',
            'hard_cut': r'Muon $|#eta|$',
            'lepton': 'tau',
            'label': r'$W^-\rightarrow\tau\nu\rightarrow\mu\nu$',
        }
    }

    my_analysis = Analysis(datasets,
                           analysis_label='test_analysis',
                           force_rebuild=False,
                           log_level=10,
                           log_out='both',
                           timedatelog=False,
                           separate_loggers=False)
    my_analysis.print_latex_table(['wminmunu', 'wmintaunu'])

    # my_analysis.apply_cuts()
    # my_analysis.merge_datasets('wminmunu', 'wmintaunu', apply_cuts=r'Muon $|#eta|$')

    my_analysis.plot_hist(['wminmunu', 'wmintaunu'],
                          'MC_WZmu_el_eta_born',
                          bins=(30, -5, 5),
                          weight='truth_weight',
                          normalise='lumi',
                          lepton='muon',
                          yerr='rsumw2',
                          ratio_plot=True)
    my_analysis.plot_hist(['wminmunu', 'wmintaunu'],
                          'MC_WZ_dilep_m_born',
                          bins=(50, 120, 5000),
示例#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),