logger.info('DATA chain name is %s ' % (tBu.GetName()))
data = []
for cutval in xrange(1,9):
data_cuts_Bu = cuts_Bu + "&& c2ip_b < %d " % cutval
# 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)
model_Bu.signal.mean .release()
model_Bu.signal.sigma.release()
r, f = model_Bu.fitHisto(h1)
r, f = model_Bu.fitHisto(h1, draw=True)
from Selectors import SBT
sel_Bu = SBT(m_Bu, cuts_Bu)
# logger.info(
# 'Build RooFit dataset for B+ , it could take as long as 3-5 minutes')
with timing():
## logger.warning('almost skip B+')
## with rooSilent () : tBu.process ( sel_Bu , 100000 )
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
sel_Bu = SelectorWithVars(
variables=selector_variables,
selection=cuts_Bu
)
logger.info(
'Build RooFit dataset for B+ , it could take as long as 3-5 minutes')