Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aefj_v1_0.tar.gz(3887 Kbytes)
Manuscript Title: Neutrino oscillation parameter sampling with MonteCUBES
Authors: Mattias Blennow, Enrique Fernandez-Martinez
Program title: MonteCUBES (Monte Carlo Utility Based Experiment Simulator).
Catalogue identifier: AEFJ_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 181(2010)227
Programming language: C.
Computer: MonteCUBES builds and installs on 32bit and 64bit Linux systems where GLoBES is installed.
Operating system: 32bit or 64bit Linux.
RAM: Typically a few MBs
Keywords: Neutrino oscillations, Long-baseline experiments, GLoBES, MonteCUBES.
PACS: 14.60.Pq.
Classification: 11.1.

External routines: GLoBES [1, 2] and routines/libraries used by GLoBES

Subprograms used:
Cat Id Title Reference
ADZI_v1_0 GLoBES CPC 177(2007)439

Nature of problem:
Since neutrino masses do not appear in the standard model of particle physics, many models of neutrino masses also induce other types of new physics, which could affect the outcome of neutrino oscillation experiments. In general, these new physics imply large dimensional parameter spaces that are difficult to explore using classical methods such as multi-dimensional projections and minimisations, such as those used in GLoBES [1, 2].

Solution method:
MonteCUBES is written as a plugin to the GLoBES software [1, 2] and provides the necessary methods to perform Markov Chain Monte Carlo sampling of the parameter space. This allows an efficient sampling of the parameter space and has a complexity which does not grow exponentially with the parameter space dimension. The integration of the MonteCUBES package with the GLoBES software makes sure that the experimental definitions already in use by the community can also be used with MonteCUBES, while also lowering the learning threshold for users who already know GLoBES.

Additional comments:
A Matlab GUI for interpretation of results is included in the distribution.

Running time:
The typical running time varies depending on the dimensionality of the parameter space, the complexity of the experiment, and how well the parameter space should be sampled. The running time for our simulations [3] with 15 free parameters at a Neutrino Factory with O(106) samples varied from a few hours to tens of hours.

[1] P. Huber, M. Lindner and W. Winter, Comput. Phys. Commun. 167 (2005) 195, hep-ph/0407333.
[2] P. Huber, J. Kopp, M. Lindner, M. Rolinec and W. Winter, Comput. Phys. Commun. 177 (2007) 432, hep-ph/0701187.
[3] S. Antusch, M. Blennow, E. Fernandez-Martinez, J. Lopez-Pavon, arXiv:0903.3986 [hep-ph].