Skip to content

Instantly share code, notes, and snippets.

View MarkusPrim's full-sized avatar

Markus Prim MarkusPrim

  • University of Bonn
  • Bonn
View GitHub Profile
@MarkusPrim
MarkusPrim / VxbOverTime.py
Created August 13, 2019 07:07
Vub and Vcb development from inclusive and exclusive measurements over time, with additional points from leptonic B decays and CKM unitarity.
import matplotlib.pyplot as plt
import numpy
import uncertainties
from uncertainties import ufloat
from wg1template.histogram_plots import create_solo_figure, add_descriptions_to_plot
from wg1template.plot_style import TangoColors
from wg1template.plot_utilities import export
from wg1template.point_plots import DataPoints, DataVariable, DataPointsPlot
"""The values are taken from from the PDG reviews of the corresponding year."""
@MarkusPrim
MarkusPrim / FlavorSensitivity.py
Created August 13, 2019 07:05
Flavor Sensitivity to new Physics
import matplotlib.pyplot as plt
import numpy
from wg1template.histogram_plots import create_solo_figure, add_descriptions_to_plot
from wg1template.plot_style import TangoColors
from wg1template.plot_utilities import export
from wg1template.point_plots import DataVariable, DataPointsPlot
def cpv_and_mixing(k, h_np=1.):
# deltaF = 2
@MarkusPrim
MarkusPrim / cramer-von-mises.py
Last active July 27, 2017 13:00
Cramer-von-Mises distance for two weighted samples
import numpy
def cvm_2samp_weigted(data1, data2, wei1, wei2):
"""Calculate Cramer-vonMises test statistics with weighted samples."""
# Indices for a sorted array
ix1 = numpy.argsort(data1)
ix2 = numpy.argsort(data2)
# Sorted data
data1 = data1[ix1]
data2 = data2[ix2]