Neutrino oscillation parameter sampling with MonteCUBES

We present MonteCUBES (“Monte Carlo Utility Based Experiment Simulator”), a software package designed to sample the neutrino oscillation parameter space through Markov Chain Monte Carlo algorithms. MonteCUBES makes use of the GLoBES software so that the existing experiment definitions for GLoBES, describing long baseline and reactor experiments, can be used with MonteCUBES. MonteCUBES consists of two main parts: The first is a C library, written as a plug-in for GLoBES, implementing the Markov Chain Monte Carlo algorithm to sample the parameter space. The second part is a user-friendly graphical Matlab interface to easily read, analyze, plot and export the results of the parameter space sampling.Program summaryProgram title: MonteCUBES (Monte Carlo Utility Based Experiment Simulator)Catalogue identifier: AEFJ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFJ_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public LicenceNo. of lines in distributed program, including test data, etc.: 69 634No. of bytes in distributed program, including test data, etc.: 3 980 776Distribution format: tar.gzProgramming language: CComputer: MonteCUBES builds and installs on 32 bit and 64 bit Linux systems where GLoBES is installedOperating system: 32 bit and 64 bit LinuxRAM: Typically a few MBsClassification: 11.1External routines: GLoBES [1,2] and routines/libraries used by GLoBESSubprograms used:Cat Id ADZI_v1_0, Title GLoBES, Reference CPC 177 (2007) 439Nature 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 high-dimensional parameter spaces that are difficult to explore using classical methods such as multi-dimensional projections and minimizations, such as those used in GLoBES [1,2].Solution method: MonteCUBES is written as a plug-in 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)O(106) samples varied from a few hours to tens of hours.References:[1]P. Huber, M. Lindner, W. Winter, Comput. Phys. Comm. 167 (2005) 195, hep-ph/0407333.[2]P. Huber, J. Kopp, M. Lindner, M. Rolinec, W. Winter, Comput. Phys. Comm. 177 (2007) 432, hep-ph/0701187.[3]S. Antusch, M. Blennow, E. Fernandez-Martinez, J. Lopez-Pavon, arXiv:0903.3986 [hep-ph].