SuperIso: A program for calculating the isospin asymmetry of B→K∗γB→K∗γ in the MSSM

We present a program for calculating the isospin symmetry breaking in B→K∗γB→K∗γ decays in the MSSM with minimal flavor violation. This program calculates the NLO supersymmetric contributions to the isospin asymmetry, using the effective Hamiltonian approach and within the QCD factorization method. We show that isospin symmetry breaking proves to be a very restrictive observable, in particular in the mSUGRA parameter space. The program also calculates the inclusive branching ratio associated to b→sγb→sγ transition, as a comparison reference.Program summaryProgram title:SuperIsoCatalogue identifier:AEAN_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAN_v1_0.htmlProgram obtainable from:CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions:Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.:3537No. of bytes in distributed program, including test data, etc.:23 595Distribution format:tar.gzProgramming language:C (C99 Standard compliant)Computer:32- or 64-bit PC, MacOperating system:Linux, MacOSRAM:less than 1 MBClassification:11.6External routines:ISASUGRA/ISAJET (http://www.hep.fsu.edu/~isajet/) and/or SOFTSUSY (http://projects.hepforge.org/softsusy)Nature of problem:  Calculation of isospin asymmetry in B→K∗γB→K∗γ decays in the Minimal Supersymmetric extension of the Standard Model with minimal flavor violation. This observable brings severe constraints on the supersymmetric parameter space.Solution method:  SuperIso uses a SUSY Les Houches Accord file, which can be either generated automatically by calling SOFTSUSY or ISAJET, or provided by the user. This file contains the masses and couplings of the supersymmetric particles. SuperIso then computes Wilson coefficients, which are used to calculate the isospin symmetry breaking in B→K∗γB→K∗γ decays, as well as the inclusive branching ratio of B→XsγB→Xsγ. SuperIso is able to perform the calculation in different supersymmetry breaking scenarios, such as mSUGRA, AMSB and GMSB.Restrictions:This calculation is valid only in the MSSM with minimal flavor violation.Unusual features:  The code is very modular, and new routines for calculating new physical observables such as the branching ratio of Bs→μ+μ−Bs→μ+μ−, can be easily added.Additional comments:The module which reads the SUSY Les Houches Accord file can be used independently in other programs.Running time:0.1 sec