示例#1
0

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 )
示例#2
0
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')