GLoBES Patrick Huber Center for Neutrino Physics, Virginia Tech - - PowerPoint PPT Presentation

▶
globes
SMART_READER_LITE
LIVE PREVIEW

GLoBES Patrick Huber Center for Neutrino Physics, Virginia Tech - - PowerPoint PPT Presentation

Aug 05, 2023 120 likes •384 views

GLoBES Patrick Huber Center for Neutrino Physics, Virginia Tech What is ? Galileo Galilei Institute for Theoretical Physics, June 2012 P. Huber p. 1 What? General Long Baseline Experiment Simulator GLoBES is a software package

slide-1
SLIDE 1

GLoBES

Patrick Huber Center for Neutrino Physics, Virginia Tech

What is ν? Galileo Galilei Institute for Theoretical Physics, June 2012

  • P. Huber – p. 1
slide-2
SLIDE 2

What? General Long Baseline Experiment Simulator

GLoBES is a software package designed for

  • Simulation
  • Analysis
  • Comparison
  • f neutrino oscillation experiments
  • P. Huber – p. 2
slide-3
SLIDE 3

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
slide-4
SLIDE 4

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
slide-5
SLIDE 5

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
slide-6
SLIDE 6

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
slide-7
SLIDE 7

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
slide-8
SLIDE 8

Why?

CERN strategy group

  • P. Huber et al.

! 35

  • P. Huber – p. 8
slide-9
SLIDE 9

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
slide-10
SLIDE 10

Why?

LBNE

  • P. Huber – p. 10
slide-11
SLIDE 11

Why?

MINOS+

e µ

θ 2

2

sin

  • 4

10

  • 3

10

  • 2

10

  • 1

10 1

2

m ∆

  • 2

10

  • 1

10 1 10

2

10

LSND 90% CL LSND 99% CL BUGEY 90% CL* MiniBooNE 90% CL MiniBooNE 99% CL Combination 90% CL MINOS+ & Bugey

courtesy of P. Huber * GLoBES 2012 fit with new reactor fluxes,

  • P. Huber – p. 11
slide-12
SLIDE 12

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. 12
slide-13
SLIDE 13

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. 13
slide-14
SLIDE 14

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. 14
slide-15
SLIDE 15

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

  • >180 publications citing the GLoBES papers,

creating a total of >3000 citations

  • fall 2012 – GLoBES 4.0
  • P. Huber – p. 15
slide-16
SLIDE 16

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. 16
slide-17
SLIDE 17

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. 17
slide-18
SLIDE 18

GLoBES – overview

GLoBES

GLoBES User Interface

Application software to compute high−level sensitivities, precision etc.

AEDL

Abstract Experiment AEDL−file(s) and simulate experiment(s) provides functions to C−library which loads

AEDL− file(s)

Defines Experiments and modifies them Definition Language

  • P. Huber – p. 18
slide-19
SLIDE 19

GLoBES – AEDL

Cross Section

Flux

Energy− function Initial / final flavor, polarity Energy

Channel

Event rates

efficiencies dependent Resolution

  • P. Huber – p. 19
slide-20
SLIDE 20

GLoBES – AEDL

Signal Background

Channel 1 Channel 2

. . . . . .

Rule

Signal + Backgrounds with systematics

∆χ2

  • P. Huber – p. 20
slide-21
SLIDE 21

GLoBES – AEDL

Rule 2 Rule 1 Rule 3

Experiment

Σ ∆χ2

. . .

  • P. Huber – p. 21
slide-22
SLIDE 22

Summary

GLoBES

  • is the only open source software of its kind
  • has withstood the test of time
  • is at the core of most strategy documents
  • completely in C
  • flexibility to deal with complex many detector

setups and non-standard physics

  • v4.0 will greatly enhance the ability to deal with

realistic systematics

  • P. Huber – p. 22
slide-23
SLIDE 23

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/personalhomes/globes/download/globes- 3.1.11.tar.gz – download GLoBES

  • ./configure

make sudo make install

  • P. Huber – p. 23
slide-24
SLIDE 24

Installation

If you have a Mac

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

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

  • Download

http://www.mpi- hd.mpg.de/personalhomes/globes/download/globes- 3.1.11.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. 24
slide-25
SLIDE 25

Installation – Mac issues

On some of them, there are problems with the architecture (32 vs. 64 bit) since Mac OS tries (and fails) to allow mixing 32 and 64 bit code. The solution to this is to compiile both the GSL and GLoBES with the -m32 option enabled. Finally, people sometimes have trouble compiling the examples, which is due to some yet unresolved problem in our autoconf script. The solution is to modify the Makefile in the examples directory in such a way that globes-config --libs is not used. Instead, use whatever options globes-config

  • -libs proposes, minus any --rpath options.
  • P. Huber – p. 25