Ejemplo n.º 1
0
from JPsi.MuMu.common.roofit import *
from JPsi.MuMu.common.plotData import PlotData

_chains = esChains.getChains('v7')

_baseCuts = [
    'abs(1/kRatio - 1) < 0.5',
    'abs(mmgMass-91.2) < 4',
]

## ----------------------------------------------------------------------------
## Customize below
leftPlot = PlotData(
    name='EB_lowR9_data_new',
    title='Barrel, R9 < 0.94, data, new corrections',
    labels=['Barrel', 'R9 < 0.94', 'data', 'new corr.'],
    source=_chains['data'],
    expression='100 * (1/newCorrKRatio - 1)',
    cuts=_baseCuts + ['phoIsEB', 'phoR9 < 0.94'],
)

rightPlot = PlotData(
    name='EB_lowR9_data_old',
    title='Barrel, R9 < 0.94, data, default corrections',
    labels=['Barrel', 'R9 < 0.94', 'data', 'default corr.'],
    source=_chains['data'],
    expression='100 * (1/kRatio - 1)',
    cuts=_baseCuts + ['phoIsEB', 'phoR9 < 0.94'],
)

ws1 = RooWorkspace('ws1', 'mmg energy scale')