from variables import *
from cuts import cuts_Bu, prntCuts, mctrue
from model import model_Bu
from data import mc_Pythia6, mc_Pythia8, mc_total
# Preparation
tBu = mc_total.data
cuts_Bu += mctrue
model_Bu.b.fix(0)
model_Bu.background.tau.fix(0)
for i in prntCuts(cuts_Bu, " CUTS B+ "):
logger.info(i)
# logger.info('Fill control B+ histogram (takes some time)')
# with timing():
# tBu.Project(h1.GetName(), 'DTFm_b', cuts_Bu)
# with rooSilent():
# logger.info('Fit Bc+ & B+ histogram (check the model)')
# r, f = model_Bu.fitHisto(h1)
from model import model_Bu
logger.info('Import data (for the first time it could take some time)')
from data import tSelection6 as tBu
cuts_s6 = " && MIPCHI2DV_k1 > 9 && MIPCHI2DV_k2 > 9 && MIPCHI2DV_k3 > 9"
cuts_Bu += cuts_s6
logger.info('DATA chain name is %s ' % (tBu.GetName()))
data_cuts_Bu = cuts_Bu
for i in prntCuts(data_cuts_Bu, " CUTS B+ "):
logger.info(i)
logger.info('Fill control B+ histogram (takes some time)')
with timing():
tBu.Project(h1.GetName(), 'DTFm_b', data_cuts_Bu)
with rooSilent():
logger.info('Fit Bc+ & B+ histogram (check the model)')
r, f = model_Bu.fitHisto(h1)
from PyPAW.Selectors import SelectorWithVars
from Variables import selector_variables