Belle II Software light-2406-ragdoll
2_analysis_1110021010.py
1#!/usr/bin/env python3
2
3
10
11"""
12<header>
13 <input>../1110021010.dst.root</input>
14 <output>1110021010_Validation.root</output>
15 <contact>Frank Meier; frank.meier@duke.edu</contact>
16</header>
17"""
18# Runs a simple analysis on signal [B0 -> rho0 gamma] events. Saves histograms
19# in a format for the validation server. There are better example analysis
20# scripts in analysis/examples and analysis/examples/tutorials.
21
22import basf2
23from modularAnalysis import inputMdst, reconstructDecay, matchMCTruth
24from stdPhotons import stdPhotons
25from stdCharged import stdPi
26from validation_tools.metadata import create_validation_histograms
27
28DECAY_MODE_ID = 1110021010
29INPUT_FILENAME = f"../{DECAY_MODE_ID}.dst.root"
30OUTPUT_FILENAME = f"{DECAY_MODE_ID}_Validation.root"
31
32main = basf2.Path()
33inputMdst(INPUT_FILENAME, path=main)
34
35stdPi('loose', path=main)
36stdPhotons('loose', path=main)
37
38reconstructDecay('rho0:pipi -> pi+:loose pi-:loose', '0.35 < M < 1.2', path=main)
39matchMCTruth('rho0:pipi', path=main)
40
41reconstructDecay('B0 -> rho0:pipi gamma:loose', '4.8 < M < 5.8', path=main)
42matchMCTruth('B0', path=main)
43
44create_validation_histograms(
45 main, OUTPUT_FILENAME, "B0",
46 variables_1d=[
47 (
48 "Mbc", 180, 5.2, 5.29,
49 "Beam-constrained mass of #it{B^{0} #rightarrow #rho^{0}#gamma} candidates",
50 "Frank Meier <frank.meier@duke.edu>",
51 r"The beam-constrained mass distribution of $B^0\to\rho^0\gamma$ decays",
52 "Distribution should be peaking at the nominal $B^0$ mass. Tail towards low mass.",
53 "M_{bc} [GeV/c^{2}]", "Candidates"
54 ),
55 (
56 "useRestFrame(daughter(1, E))", 50, 1.5, 4.0,
57 "The photon energy distribution for #it{B^{0} #rightarrow #rho^{0}#gamma} decays",
58 "Frank Meier <frank.meier@duke.edu>; Torben Ferber <torben.ferber@desy.de>",
59 r"Photon energy distribution of $B^0\to\rho^0\gamma$ decays",
60 "Sharp distribution at 2.5 GeV. Look for differences in "
61 "the absolute scale. This could be an indication of a loss of photon efficiency.",
62 "E_{#gamma} [GeV]", "Candidates"
63 ),
64 ]
65)
66
67main.add_module('Progress')
68basf2.process(main)
69print(basf2.statistics)