test2_feiSLB0.py

#!/usr/bin/env/python3
# -*-coding: utf-8-*-
 
"""
<header>
    <input>../feiSLB0.udst.root</input>
    <output>feiSLB0_Validation.root</output>
    <contact>philip.grace@adelaide.edu.au</contact>
</header>
"""
 
import basf2 as b2
import modularAnalysis as ma
from variables import variables as vm
from validation_tools.metadata import create_validation_histograms
 
path = b2.Path()
 
ma.inputMdst('default', '../feiSLB0.udst.root', path=path)
 
vm.addAlias('sigProb', 'extraInfo(SignalProbability)')
vm.addAlias('log10_sigProb', 'log10(extraInfo(SignalProbability))')
vm.addAlias('d0_massDiff', 'daughter(0,massDifference(0))')
vm.addAlias('d0_M', 'daughter(0,M)')
vm.addAlias('decayModeID', 'extraInfo(decayModeID)')
vm.addAlias('nDaug', 'countDaughters(1>0)')  # Dummy cut so all daughters are selected.
 
histogramFilename = 'feiSLB0_Validation.root'
myEmail = 'Phil Grace <philip.grace@adelaide.edu.au>'
 
create_validation_histograms(
    rootfile=histogramFilename,
    particlelist='B0:semileptonic',
    variables_1d=[
        ('sigProb', 100, 0.0, 1.0, 'Signal probability', myEmail,
         'Signal probability of the reconstructed tag B candidates', 'Most around zero, with a tail at non-zero values.',
         'Signal probability', 'Candidates', 'logy'),
        ('nDaug', 6, 0.0, 6, 'Number of daughters of tag B', myEmail,
         'Number of daughters of tag B', 'Some distribution of number of daughters', 'n_{daughters}', 'Candidates'),
        ('cosThetaBetweenParticleAndNominalB', 100, -6.0, 4.0, '#cos#theta_{BY}', myEmail,
         'Cosine of angle between the reconstructed B and the nominal B', 'Distribution peaking between -1 and 1',
         '#cos#theta_{BY}', 'Candidates'),
        ('d0_massDiff', 100, 0.0, 0.5, 'Mass difference of D* and D', myEmail,
         'Mass difference of $D^{*}$ and D', 'Peak at 0.14 GeV', 'm(D^{*})-m(D) [GeV]', 'Candidates', 'shifter'),
        ('d0_M', 100, 0.0, 3.0, 'Mass of zeroth daughter (D* or D)', myEmail,
         'Mass of zeroth daughter of tag B (either a $D^{*}$ or a D)', 'Peaks at 1.86 GeV and 2.00 GeV',
         'm(D^{(*)}) [GeV]', 'Candidates', 'shifter')],
    variables_2d=[('decayModeID', 8, 0, 8, 'log10_sigProb', 100, -3.0, 0.0,
                   'Signal probability for each decay mode ID', myEmail,
                   'Signal probability for each decay mode ID',
                   'Some distribtuion of candidates in the first few decay mode IDs',
                   'Decay mode ID', '#log_10(signal probability)', 'colz')],
    path=path)
 
b2.process(path)
print(b2.statistics)