Hola!

Lo lamento, pero me parece que no va a ser nada facil digerir lo que os
envio. El programa que creé, en un intento de ser suficientemente
flexible, resulta bastante complejo de manejar. De todas formas os lo
mando sin "censurar". Quiero decir que podría desactivar cosas para
hacer este "manual" más simple, pero que prefiero daros el todo.

El attachment contiene el archivo toCarlos.uu (para desempaquetarlo
teclead "csh toCarlos.uu", preferiblemente en un directorio vacio). A su
vez, este contiene los siguientes ficheros:

bsgexcl.F
BRbsg.in
bsgcore.F 
bsgprocs.f
bsgfuncs.f
Makefile

Os explico un poco como funciona todo y que hay que hacer para empezar a
producir numeros.

0. bsgexcl.F es la interfaz que produce las tablas de resultados, los
   plots, lo que haga falta. Lee los inputs del fichero BRbsg.in,
   escribe los resultados en xyBRbsg.dat, exclBRbsg.dat y
   bestofBRbsg.dat y reporta los errores/incidencias que encuentra (una
   combinacion de inputs que viola algun bound, por ejemplo) en
   BRbsg.err.

   0.1 xyBRbsg.dat: fichero de salida PRINCIPAL. Creo que podeis
     olvidar los otros dos (0.2, 0.3). Una primera linea muestra el
     nombre de cada columna, empezando por las variables X, Y (ved mas
     adelante) y siguiendo por el BR(b->s\gamma) al LO en el MSSM, el
     BR en el MSSM, el BR en el improved MSSM y TODOS los
     inputs que especifican este resultado. Los BR REALES se obtienen
     dividiendo por 10^4! Las lineas del fichero son bastante largas,
     tendreis que ensanchar mucho la ventana para verlas sin
     cortar.

   0.2 exclBRbsg.dat: fichero pensado para hacer contour plots de
     exclusion (zonas dentro y fuera de los limites
     experimentales). Tiene tres columnas: X, Y son dos variables
     genericas que se especifican al ejecutar el programa
     bsgexcl.exe (ved mas abajo). La tercera columna contiene el calculo
     del BR en el modelo (2HDM, MSSM, impMSSM) especificado al ejecutar
     bsgexcl.exe. El programa busca el punto mas proximo al valor
     medio experimental para ese modelo concreto. Si el par X,Y da un
     resultado fuera del rango experimental, la 3ra columna contiene
     el valor 10^5 (10^4 * BR siempre sera menor que 10^5). Este fichero
     solo tiene sentido si se utiliza Inofit=0 (ver mas abajo), es
     decir, cuando se esta haciendo un scan dentro de un rango de
     parametros. Yo creo que os podeis olvidar de el.

   0.3 bestofBRbsg.dat: esta relacionado con el anterior. Especifica
     los inputs empleados para el mejor punto del scan, y los
     resultados al LO y NLO, para el SM, el MSSM y el 2HDM para ese
     punto concreto.

   0.4 BRbsg.err: incidencias producidas durante el calculo. Por
     ejemplo, internamente se supone que los stops estan ordenados
     m_stop2<m_stop1. Si esto no es asi se escribe una linea en este
     fichero. Si un error se repite mucho se para deja de notificar para
     no ser demasiado pesado.

   Al compilar el fichero bsgexcl.F (ved como en el punto 5), se
   genera un ejecutable de nombre bsgexcl.exe. Al correr este
   programa, os pide algunos datos que no se encuentran en BRbsg.in:

   - El modelo utilizado en el fit: si tecleais 2 equivale al 2HDM, 3
   para MSSM, 4 para el improved MSSM. Si Inofit=0, el programa trata
   de minimizar la diferencia entre el modelo y el valor central
   experimental. Esto solo es importante para los ficheros
   exclBRbsg.dat y bestofBRbsg.dat. Los resultados en xyBRbsg.dat NO
   tendrian que depender de esto.

   - Las variables X e Y, dentro del conjunto tanb (teclead 1), mHc
     (teclead 2), McL (3), mstL (5), mhv (8). Teneis que entrar dos
     numeros, el primero es el eje x, el segundo el y. Si seleccionais
     "NPoints= 1 1", esto no tiene ningun efecto.

1. BRbsg.in es el fichero de entrada de datos para bsgexcl.F. Fijaos
   en el formato, HAY QUE RESPETARLO! Los inputs estan AGRUPADOS por
   SECCIONES, con nombres mas o menos acertados: SM, tHDM (vale por
   2HDM), MSSM, MSSMFlags y XY. Los grupos estan separados por una
   cabecera tipo &NOMBRE (e.g., &SM) y acaban en &END. No permuteis
   las secciones ni mezcleis los inputs! Los datos se especifican en
   la forma VARIABLE=VALOR, y van separados dentro de un mismo grupo
   por comas (excepto el ultimo).

   La mayoria de variables en BRbsg.in se pasan tal cual a la rutina
   principal, NLObsgMSSM, contenida en el fichero bsgcore.F (os la
   explico en el punto 2). Las pocas excepciones son:

   1.1 INLO: cuando INLO=1, los campos MSSM y impMSSM del fichero
     xyBRbsg.dat son calculados al NLO. Si INLO=0, el campo MSSM es
     igual al campo LO MSSM (calculo LO) i el impMSSM contiene el LO
     con el improvement del Yukawa del bottom.

   1.2 Inofit: especifica si se va a hacer un scan dentro de un rango
     de parametros, Inofit=0, o si se calcula para un unico punto
     Inofit=1. El caso del scan (quizas deberia decir fit) es
     complicado, y de momento no os doy mas detalles. DEJAD PUES
     Inofit=1! De hecho estaba pensando en desactivar toda esta parte
     antes de pasaros el programa...

   1.3 NPoints: numero de puntos en los ejes X e Y. Si poneis
     "NPoints= 1 1", como en la muestra que os mando, solo se calcula
     para la combinacion de parametros especificada en el resto de
     grupos de datos. SI EN CUALQUIER EJE PEDIS MAS DE 1 PUNTO, se
     utilizan unas variables X e Y que van tomando valores entre "X0 =
     x0 y0" y "X1 = x1 y1". Estas variables X, Y SUSTITUYEN (en ingles
     yo diria "override") las entradas de cualesquiera otras variables
     que hayais especificado al ejecutar el programa bsgexcl.exe.

   1.4 X0, X1: como acabais de leer, X0 contiene los valores iniciales
     de las variables genericas X e Y; X1 los valores finales.

   Nota bene: el formateado de los grupos depende del
     compilador/maquina. El formato que os paso me funciona con los
     compiladores g77, f77 y pgf77, pero puede que falle en otro
     compilador/maquina.

2. bsgcore.F  es el nucleo del programa, contiene la rutina 

      Subroutine NLObsgMSSM(IInit,dBS,IQED,IappDmb,IthAt,IMcmu !! INPUTS
     &           ,tb,mHc,McL,McH,IsM,IsC,mstL,mstH,thstAt,mhv  !! INPUTS
     &           ,LOBRg,BRg,IError,IMsg)                       !! OUTPUTS

que es la que realmente calcula todo y que es llamada desde
bsgexcl.F. Los inputs de la rutina son:

------ INPUTS ------

 2.1 IInit: hay que ejecutar NLObsgMSSM una primera vez con IInit=1
   para inicializar una serie de constantes. A partir de ahi se debe
   utilizar siempre IInit=0. Esto esta ya implementado automaticamente
   en bsgexcl.F

 2.2 dBS: valor del parametro \delta relacionado con el calculo de los
   efectos de bremsstrahlung. Aunque es una variable de entrada, NO
   GARANTIZO el resultado si no utilizais dBS=0.9. Tendria que
   revisarlo con calma para asegurarme que se pueden utilizar otros
   valores. Fijaos que \delta=0.99 en el 1er paper de Ciuchini et al!

 2.3 IQED: si IQED=1 se incluyen las correcciones QED, si IQED=0 no se
   incluyen

 2.4 IappDmb: tiene que ser =1. Si IappDmb=0 se toma Dmb=0.

 2.5 IthAt: si IthAt=0, la variable thstAt representa el angulo de
   mixing de los stops, thst; si IthAt=1, thstAt representa At, el
   valor del trilineal. Un poco extravagante, quizas, pero permite 
   escoger At o thst como input sin necesidad de recompilar cada vez.

 2.6 IMcmu: mas de lo mismo. Si IMcmu=0, las variables McL, McH
   representan las masas del chargino Light y Heavy,
   respectivamente. Si IMcmu=1, entonces |McH| es M2,
   |McL|*Sign(McH*McL) es \mu. Fijaos que el signo de \mu esta en el
   PRODUCTO de McH*McL.

 2.7 tb, mHc: tan\beta, masa del H+

 2.8 McL, McH: ver punto 2.6

 2.9 IsM, IsC: si se decide trabajar con las masas de los charginos
   como inputs (IMcmu=0) aparece una indeterminacion quadruple al
   deducir los valores de \mu, M2 correspondientes. En este caso, IsM,
   IsC contienen los signos (IsM,IsC=1 o -1) adecuados. Si IMcmu=1
   estos valores se ignoran. Si se trabaja con IMcmu=0 y se hace un
   fit (Inofit=1), el interfaz bsgexcl.F se encarga de PROBAR las
   cuatro cobinaciones antes de aceptar(el resultado esta dentro de la
   horquilla experimental)/rechazar(el resultado esta excluido
   experimentalmente) un punto del espacio de parametros.

 2.10 mstL, mstH: massa del stop Light y Heavy. 

 2.11 thstAt: ver 2.5

 2.12 mhv: el calculo de Ciuchini et al. parte de un lagrangiano
   efectivo con una serie de particulas pesadas (el gluino, los
   sbottoms, un stop, un chargino) y otras ligeras (un stop y un
   chargino, el W...). A todas aquellas particulas pesadas para las
   cuales NO SE ESPECIFICA SU MASA EXPLICITAMENTE, se les asigna mhv
   (mheavy) como masa.

 WARNING: las variables mstH, McH y thstAt se interpretan, cuando
   Inofit=0, como los LIMITES SUPERIORES del scan de parametros (para
   thst, At, se toma el valor absoluto). De todas maneras, los valores
   "reales" de mstH y McH aparecen escritos para cada punto en el
   fichero xyBRbsg.dat.

------ OUTPUTS ------

 2.13 LOBRg: array conteniendo el calculo al LO. Toma valores desde
   LOBRg(1) hasta LOBRg(4), correspondientes a los modelos: SM, 2HDM,
   MSSM, improved MSSM.

 2.14 BRg: array BRg(1->4) que contiene el calculo al NLO.

 2.15 IError: IError=0 si no pasa nada especial, IError!=0 si ha
   habido algun problema. El valor indica el tipo de problema, si
   utilizais bsgexcl.F como interfaz, podeis ver qué ha ido pasando en
   el fichero BRbsg.err .


3. bsgprocs.f: fichero que contiene subrutinas utilizadas en bsgcore.F

4. bsgfuncs.f: lo mismo pero para funciones

5. Makefile: fichero de datos que se utiliza para facilitar el proceso
   de compilacion. Y sí, lo facilita, aunque no lo creais. Para compilar
   bsgexcl.F, teneis que poneros en el directorio donde este el
   Makefile (y el resto de ficheros .f y .F) y teclear

   make

   Si algo falla, enviadme el output y os dire de que se puede
   tratar. Lo mas probable es que haga falta editar las lineas que
   especifican el compilador (buscad "FC" en Makefile) o las opciones
   de compilacion ("FFLAGS" en Makefile).

   Si todo va bien, tendrian que aparecer varios ficheros .o (codigo
   objeto, los podeis olvidar) y un ejecutable: el programa

   bsgexcl.exe 

   Solo hace falta que pongais los inputs que mas os gusten en
   BRbsg.in (ver puntos 1 y 2) y ejecuteis el programa. bsgexcl.exe
   pide unos datos adicionales por pantalla, ved punto 0. Si algo falla
   aqui, puede ser problema del formato del fichero BRbsg.in.



Y bueno, si habeis llegado hasta aqui, enhorabuena! 


Saludos,


David
_________________________________________________________
 David Garcia Muchart       dgarcia@mail.cern.ch
 Theory Division, CERN      Phone(Off.): +41 22 767 2466
 1211 Geneva 23             Fax:         +41 22 767 3850
                            Phone(Home): +33  4 50408939 
_________________________________________________________