def CopyHists(infile, outfile, prefix = ''):
from load_data import TprimeData
_hists = TprimeData(1)
_hists.load_all_hists(infile, prefix)
_hists.SaveAllHists(outfile)
return
def ProjectOnHt(infile, outfile, prefix = ''):
from load_data import TprimeData
_hists = TprimeData(1)
_hists.load_all_hists(infile, projection_y = True,
rebin_ngroup = 4,
rebin_low = 200,
rebin_high = 2000,
prefix=prefix)
_hists.SaveAllHists(outfile)
return
def MakeSignalBoxDiagonal(infile, outfile, prefix = '',
mass = None,
signal_only = False):
from load_data import TprimeData
_hists = TprimeData(1)
_hists.load_all_hists(infile,
prefix=prefix)
_hists.ConvertToSignalDiagonal(mass,signal_only=signal_only)
_hists.SaveAllHists(outfile)
return
def MakeSignalBoxTemplates(infile, outfile, prefix = '',
xmin = None, xmax = None,
ymin = None, ymax = None,
signal_only = False):
from load_data import TprimeData
_hists = TprimeData(1)
_hists.load_all_hists(infile,
prefix=prefix)
_hists.ConvertToSignalBox(xmin, xmax, ymin, ymax, signal_only=signal_only)
_hists.SaveAllHists(outfile)
return
def MakeSignalMassDiagPerp(infile, outfile, prefix = '',
mass = None,
rebin_ngroup = None,
rebin_low = None,
rebin_high= None
):
from load_data import TprimeData
_hists = TprimeData(1)
_hists.load_all_hists(infile,
prefix=prefix)
_hists.ConvertToMassDiagPerp(mass,rebin_ngroup,rebin_low,rebin_high)
_hists.SaveAllHists(outfile)
return
def CreateTemplates(options):
#
# create two "standard" files with "plus" and "minus"
# shape templates for POWHEG-Nominal difference
# Since there is only one POWHEG template, the idea is
# to interpolate in some way between it and the nominal one
#
legend = '[tprime_powheg.CreateTemplates]:'
channels = ['e', 'mu']
input_nom = {}
input_pow = {}
input_nom[
'mu'] = 'data/mujets_4600/from_gueorgi/27jan2011v1/csv/v1/mujets_4601ipb_2D_v1_csv_jer_nom.root'
input_pow[
'mu'] = 'data/mujets_4600/from_gueorgi/27jan2011v1/csv/v1/mujets_4601ipb_2D_v1_csv_POWHEG.root'
input_nom[
'e'] = 'data/ejets_4700/ricardo_4683ipb_2d_csv_22jan2012v2/ht4jetsAfterFit_vs_fitMass.root'
input_pow[
'e'] = 'data/ejets_4700/ricardo_4683ipb_2d_csv_22jan2012v2/ht4jetsAfterFit_vs_fitMass_powheg.root'
hist_name = {}
hist_name['e'] = 'Top_ht4jetsAfterFit:fitMass'
hist_name['mu'] = 'TTjets_HtvsMfit'
for chan in channels:
data = TprimeData()
data_powheg = TprimeData()
_suffix = ''
data.load_all_hists(input_nom[chan],
merge_ewk=False,
merge_top=False,
strip_suffix=False,
suffix=_suffix)
data_powheg.load_all_hists(input_pow[chan],
merge_ewk=False,
merge_top=False,
strip_suffix=False,
suffix=_suffix)
hPlus = copy.deepcopy(data.hists[hist_name[chan]])
hMinus = copy.deepcopy(data_powheg.hists[hist_name[chan]])
hPlus.SetMarkerColor(ROOT.kRed)
#hMinus.Draw()
#hPlus.Draw('same')
#raw_input('press enter to continue')
hNom = copy.deepcopy(hMinus)
hNom.Add(hPlus)
hNom.Scale(0.5)
# save nominal
data.hists.pop(hist_name[chan])
iTop = 0
for name in data.hist_names:
if name == hist_name[chan]:
break
iTop += 1
data.hist_names.pop(iTop)
data.AddHist(hNom)
data.SaveAllHists(chan + '_powheg_nom.root')
# save minus
data.hists.pop(hist_name[chan])
iTop = 0
for name in data.hist_names:
if name == hist_name[chan]:
break
iTop += 1
data.hist_names.pop(iTop)
data.AddHist(hMinus)
data.SaveAllHists(chan + '_powheg_minus.root')
# save plus
data.hists.pop(hist_name[chan])
iTop = 0
for name in data.hist_names:
if name == hist_name[chan]:
break
iTop += 1
data.hist_names.pop(iTop)
data.AddHist(hPlus)
data.SaveAllHists(chan + '_powheg_plus.root')
projection_x=options.projection_x,
projection_y=options.projection_y,
rebin_ngroup=options.rebin_ngroup,
rebin_low=options.rebin_low,
rebin_high=options.rebin_high,
scale_bin_error=options.scale_bin_error)
i_suff += 1
# add pseudo data if requested
if options.add_pseudodata:
_map = {}
_map['ewk'] = 380
_map['top'] = 1766
ds.CreatePseudoData('DATA', _map)
ds.SaveAllHists(options.out_file)
if name == 'mujets_merge':
ds.load_all_hists('data/mujets/5jets/v2/mujets_1889ipb_2D_v2.root',
merge_ewk=merge_ewk)
ds.load_all_hists(
'data/mujets/5jets/v2/mujets_1889ipb_2D_v2_JES_minus.root',
"",
"_JESdown",
merge_ewk=merge_ewk)
ds.load_all_hists(
'data/mujets/5jets/v2/mujets_1889ipb_2D_v2_JES_plus.root',
"",
"_JESup",
merge_ewk=merge_ewk)
ds.SaveAllHists('data/mujets/tprime_mujets_2D_1889ipb_5j_v2.root')
