Validate_B2Vlnu.py

#!/usr/bin/env python3
 
 
 
"""
<header>
  <contact>P. Urquijo phillip.urquijo@unimelb.edu.au</contact>
  <output>Validate_B2Vlnu.root</output>
</header>
"""
 
import basf2
from generators import add_evtgen_generator
from modularAnalysis import findMCDecay
from variables import variables
from validation_tools.metadata import create_validation_histograms
 
basf2.set_random_seed('B+ => [anti-D*0 -> anti-D0 pi0] e+ nu_e')
 
variables.addAlias('q2', 'formula(daughterInvM(1, 2)**2)')
 
path = basf2.Path()
path.add_module('EventInfoSetter', evtNumList=[10000])
add_evtgen_generator(path, 'signal', basf2.find_file('decfiles/dec/1263210000.dec'))
findMCDecay('B+:sig', 'B+ -> [anti-D*0 -> anti-D0 pi0] e+ nu_e', path=path)
create_validation_histograms(
    path, 'Validate_B2Vlnu.root', 'B+:sig',
    [
        ('q2', 50, -0.5, 12.5, '', 'P. Urquijo <phillip.urquijo@unimelb.edu.au>',
         r'B2Vlnu $q^2$ of the $e^+\nu_{e}$ system in $B^+ \to \bar{D^{*0}} e^+ \nu_e$ (truth values)',
         'should follow the reference', '#q^2'),
    ],
    description='B2Vlnu Validation to check generator level distributions',
)
path.add_module('Progress')
basf2.process(path)