Esempio n. 1
0
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)

Esempio n. 2
0
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