def get_electron_normalisation(met_bin, b_tag):
global electron_data_file
input_file = File(electron_data_file)
histogram_for_estimation = 'TTbarPlusMetAnalysis/EPlusJets/QCD e+jets PFRelIso/BinnedMETAnalysis/Electron_patType1CorrectedPFMet_bin_%s/electron_pfIsolation_03_%s' % (met_bin, b_tag)
input_histogram = input_file.Get(histogram_for_estimation)
result = estimate_with_fit_to_relative_isolation(input_histogram)
value, error = result['value'], result['error']
return value, error
def get_electron_normalisation(met_bin, b_tag):
global electron_data_file
input_file = File(electron_data_file)
histogram_for_estimation = 'TTbarPlusMetAnalysis/EPlusJets/QCD e+jets PFRelIso/BinnedMETAnalysis/Electron_patType1CorrectedPFMet_bin_%s/electron_pfIsolation_03_%s' % (
met_bin, b_tag)
input_histogram = input_file.Get(histogram_for_estimation)
result = estimate_with_fit_to_relative_isolation(input_histogram)
value, error = result['value'], result['error']
return value, error
1528. / 0.48,
'QCDStudy/PFIsolation_3jets_WithMETCutAndAsymJetCuts_1orMoreBtag':
491. / 0.67,
'QCDStudy/PFIsolation_WithMETCutAndAsymJetCuts_DR03_0btag':
1398. / 0.61,
'QCDStudy/PFIsolation_WithMETCutAndAsymJetCuts_DR03_1btag':
504. / 0.79,
'QCDStudy/PFIsolation_WithMETCutAndAsymJetCuts_DR03_2orMoreBtags':
210 / 0.73
}
QCD_rate_estimation.relative_isolation_bias = 0.0
results = {}
for histogram_for_estimation in histograms_for_estimation:
input_histogram = input_file.Get(histogram_for_estimation)
result = estimate_with_fit_to_relative_isolation(input_histogram)
results[histogram_for_estimation] = result
#QCD normalised to lumi
input_histogram_qcd = File(histogram_files['QCD']).Get(
histogram_for_estimation.replace('_DR03_', '_'))
input_histogram_qcd.Scale((1.0 * lumi) / 1959.75)
input_histogram_qcd.Rebin(10)
nQCD = input_histogram_qcd.GetBinContent(1)
# print nQCD
# nQCD = qcd_expected[histogram_for_estimation]
fQCD = '---'
if not nQCD == 0:
fQCD = (result['value']) / nQCD
print histogram_for_estimation, ':', result['value'], '+-', result[
'error'], 'QCD: ', nQCD, ',f_QCD:', fQCD
''' Created on 15 Jan 2013 @author: kreczko This module will provide estimates for the number of QCD events after event selection. It will be able to do the closure test (option --do-closure-test) as well get the calibration curve for the estimates (option --make-calibration-curves). In addition to above it should provide an easy way to present the binned estimates (MET, b-tag, other bins) ''' from tools.QCD_rate_estimation import estimate_with_fit_to_relative_isolation from rootpy.io import File path_to_files = '/storage/TopQuarkGroup/results/histogramfiles/AN-12-241_V4/' electron_data_file = path_to_files + 'central/ElectronHad_5050pb_PFElectron_PFMuon_PF2PATJets_PFMET.root' input_file = File(electron_data_file) histogram_for_estimation = 'TTbarPlusMetAnalysis/EPlusJets/QCD e+jets PFRelIso/BinnedMETAnalysis/Electron_patType1CorrectedPFMet_bin_0-25/electron_pfIsolation_03_0btag' input_histogram = input_file.Get(histogram_for_estimation) result = estimate_with_fit_to_relative_isolation(input_histogram) print result
