alexpearce

\begin{tikzpicture}
\pgfdeclareplotmark{cross} {
\pgfpathmoveto{\pgfpoint{-0.3\pgfplotmarksize}{\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+0.3\pgfplotmarksize}{\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+0.3\pgfplotmarksize}{0.3\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+1\pgfplotmarksize}{0.3\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+1\pgfplotmarksize}{-0.3\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+0.3\pgfplotmarksize}{-0.3\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{+0.3\pgfplotmarksize}{-1.\pgfplotmarksize}}
\pgfpathlineto{\pgfpoint{-0.3\pgfplotmarksize}{-1.\pgfplotmarksize}}

alexpearce

# Add this directory to PATH
import sys
import os
sys.path.append(os.getcwd())

from Configurables import (
    DaVinci,
    DecayTreeTuple,
    GaudiSequencer,
    ProcStatusCheck,

alexpearce

alexpearce

"""Test RooDataSet handling of TTree branch aliases.

The TTree::SetAlias method allows you to create "branch aliases", names that
point to formulas of other branches.
It seems RooFit will not load alias branches when importing a TTree in to a
RooDatSet, filling the variable with zeros.

This script demonstrates that, firstly by fitting the non-aliased branch 'x',
which works fine, and then fitting the alias to x, 'y', which fails because the
values are all zero.

alexpearce

#!/usr/bin/env python
from __future__ import print_function
import argparse
import os
import re
import ssl
import sys
import urllib
from xml.dom import minidom

alexpearce

# -*- coding: utf-8 -*-
from __future__ import print_function
from math import sqrt, sin
from array import array

import ROOT

# Proton mass in MeV
PROTON_MASS = 938.27

alexpearce

# Name of the delta mass variable in the workspace
dmass_var = 'Dst_delta_M'
# Species of the PDF we're about to build
species = 'sig'
workspace.factory('mu_dm_{0}[141, 146]'.format(species))
workspace.factory('sigma_one_dm_{0}[1, 0, 5]'.format(species))
workspace.factory('nsigma_two_dm_{0}[1.5, 1, 3]'.format(species))
workspace.factory('nsigma_three_dm_{0}[1.5, 1, 3]'.format(species))
workspace.factory((
    'expr::sigma_two_dm_{0}('

alexpearce

#!/usr/bin/env python
"""Inspect ROOT files."""
from __future__ import print_function

import argparse
import sys

parser = argparse.ArgumentParser(description=__doc__)
parser.add_argument('files', metavar='file', nargs='+',
                    help='ROOT file to inspect (ending in .root)')

alexpearce

from distutils.spawn import find_executable
import logging as log
import os
import shutil
from subprocess import call
import tempfile


class open_eos(object):
    """An EOS file wrapper that acts like the native `open`.

alexpearce