GLoBES Patrick Huber IPNAS, Virginia Tech International Neutrino - - PowerPoint PPT Presentation

GLoBES

Patrick Huber IPNAS, Virginia Tech

International Neutrino Summer School Fermilab, July 6-17, 2009

• P. Huber – p. 1

What? General Long Baseline Experiment Simulator

GLoBES is a software package designed for

• Simulation
• Analysis
• Comparison
• f neutrino oscillation experiments
• P. Huber – p. 2

Where?

It is developed and maintained by

• PH
• Joachim Kopp
• Manfred Lindner
• Walter Winter

URL – http://www.mpi-hd.mpg.de/lin/globes/ email – globes@mpi-hd.mpg.de

• P. Huber – p. 3

Why?

Wide band beam

Resolution δCP vs Sin22θ13

• 150
• 100
• 50

50 100 150 0.02 0.04 0.06 0.08 0.1 0.12

Sin22θ13 δCP (deg.)

νµ Running Only 90 % C.L. 3 cntrs: STAT only STAT+10% SYST STAT+20% SYST FNAL-HS 1290 km sin22θij (12,23,13) = 0.86/1.0/0.04, δCP=45o ∆mij2 (21,32) = 7.3e-5/2.5e-3 eV2 1 MW, 0.5 MT, 5e7 sec

103 102 101 True value of sin22Θ13 150 100 50 50 100 150 True value of ∆CP

• P. Huber – p. 4

Why?

APS study

105 104 103 102 101 sin22Θ13

Sensitivity reach in sin22Θ13

CP viol. sgnm31

2

sin22Θ13 CP viol. sgnm31

2

sin22Θ13 CP viol. sgnm31

2

sin22Θ13

7500 km 3000 km 7500 km 7500 km 3000 km 3000 km

No sensitivity

• P. Huber – p. 5

Why?

Fermilab’s Proton driver report

2005 2010 2015 2020 2025 2030 Year 105 104 103 102 101 100 sin22Θ13 discovery reach 3Σ CHOOZSolar excluded Branching point

• Conv. beams

SuperbeamsReactor exps Superbeam upgrades Νfactories MINOS CNGS DCHOOZ T2K NOA ReactorII NOAFPD 2ndGenPDExp NuFact

• P. Huber – p. 6

Why?

White paper on reactor neutrinos

10

2

10

3

10

4

10

5

luminosity [t GW y] 10

• 2

sin

22θ13 sensitivity limit

10

2

10

3

10

4

10

5

10

• 2

no BG, no bin-to-bin error 0.5% 5% 1% Reactor-I Reactor-II

. 1 %

BG in far detector bin-to-bin error

• P. Huber – p. 7

Why?

CERN strategy group

• P. Huber et al.

! 35

• P. Huber – p. 8

Why?

ISS

105 104 103 102 101 True value of sin22Θ13 0.2 0.4 0.6 0.8 1 Fraction of ∆CP SPL T2HK WBB NF BB

GLoBES 2006

• P. Huber – p. 9

Why?

Joint BNL-FNAL study group

NOΝA T2KK WBBWC WBBLAr

0.5 1 10 100 2 20 5 50 10.3. 10.2. exposure Mt MW 107s sin22Θ13 10.2. sin22Θ13 10.2. sin22Θ13 10.2. 10.1. sin22Θ13 10.3. 10.3. 10.3.

sin22Θ13 CPV CPV sgnm2

• P. Huber – p. 10

Reliability

• Re-use of code, the more a code has been used in

real world applications the less likely are severe bugs.

• Extensive testing
• Good documentation
• Intuitive API with error checking
• P. Huber – p. 11

Reproducibility

The information given in a publication or proposal is not sufficient to reproduce the sensitivity estimates.

• General data storage and exchange format for the

inputs ⇔ flexibility?

• All implicit assumptions and approximations

have to be documented, that includes the actual algorithms ⇔ accuracy of documentation?

• Version control and archiving
• P. Huber – p. 12

Documentation

Without good documentation, the best software is useless or will be after very short time (=memory decay constant of typical physicist). This is a general problem with legacy code! Document what you do – do what you document and make sure that the average user understands what is going on. Also documentation needs testing and debugging.

• P. Huber – p. 13

GLoBES history

• development started 2004 – PH, M. Lindner, W.

Winter

• major effort went into documentation
• first release August 2004 – version 2.0.0
• major bug fix release March 2005 version 2.0.11
• J. Kopp joined in July 2005
• January 2007 – version 3.0, addition of major

features

• 93 publications citing the GLoBES paper,

creating a total of 1514 citations

• P. Huber – p. 14

Design considerations

• GPL
• C-library – very portable, easy to interface,

numerically efficient

• Unix style separation of functionality – freedom

to design analysis and to use any graphics tools

• Experiments are defined using AEDL – relatively

complicated parser, transparent experiment definition

• Pull approach for systematics – flexible and

intuitive

• Local minimization instead of grids – much faster
• P. Huber – p. 15

Features

• Accurate treatment of systematical errors
• Arbitrary matter profile & uncertainties
• Arbitrary energy resolution function
• Single and multiple experiment simulation
• Simple χ2 calculation
• Inclusion of external input
• Projection of χ2 (minimization)
• User-defined systematics, oscillation probability

engine, priors

• Full support for lists in AEDL
• Interpolating functions in AEDL
• . . .
• P. Huber – p. 16

Summary

GLoBES

• is the only open source software of its kind
• has withstood the test of time (netx month, 5

years!)

• is at the core of most strategy documents
• completely in C
• flexibility to deal with complex many detector

setups and non-standard physics

• P. Huber – p. 17

Installation

If you have Linux

• Install GSL if you don’t have it already –

ftp://ftp.gnu.org/gnu/gsl/

• Go to

http://www.mpi-hd.mpg.de/lin/globes/downloads – download GLoBES

• ./configure

make sudo make install

• P. Huber – p. 18

Installation

If you have a Mac

• Install GSL if you don’t have it already –

ftp://ftp.gnu.org/gnu/gsl/

• Download

http://www.phys.vt.edu/ pahuber/globes- 3.1.3.tar.gz

• unpack it and change into the directory created by

this

./configure --disable-rpath --enable-no-binary=yes

make make install

• P. Huber – p. 19