MadGraph
previous 4 of 6 next




Simulación Hadrónica

Como se dijo anteriormente, la simulación hadrónica, se hace luego de que MadEvent(ME) genera los eventos a nivel partónico. En éste tutorial se utilizará el paquete Pythia8 previamente instalado. Recordando que Pythia se puede ejecutar directamente desde MG, creando una tarjeta de configuración pythya8.dat en la carpeta Cards y ejecutando ./bin/mg5_aMC, o ejecutando en la carpeta de salida ./bin/madevent opciones 0 y 0. A continuación se muestran las tres opciones suponiendo una carpeta de salida (output) llamada "InclusiveTTbar1J".

Directamente desde MG:
Primero se ejecuta MG,
./PATH_TO_MADGRAPH/bin/mg5_aMC

Luego se introducen cuales procesos se quieren simular 
define p = u u~ d d~ c c~ s s~ b b~ g
define j = p
generate p p > t t~ @0
add process p p > t t~ j @1
output InclusiveTTbar1J
launch

Nota: En la primera linea se define el elemento j (jet) como una multi-partícula igual al proton (se puede definir como se quiera).
En la segunda linea se pide que el proceso que se genere sea el de dos protones yendo a un par quarks top anti-top, donde "@0" quiere decir que es no hay elementos extras.
Por último en la tercera linea se le pide genere eventos de dos protones yendo a un par quarks top anti-top, donde "@1" quiere decir que el jet es un elemento extra.

A continuación, MadGraph mostrará un set de opciones, así: 
************************************************************
*                                                          *
*                      W E L C O M E to                    *
*             M A D G R A P H 5 _ a M C @ N L O            *
*                      M A D E V E N T                     *
*                                                          *
*                 *                       *                *
*                   *        * *        *                  *
*                     * * * * 5 * * * *                    *
*                   *        * *        *                  *
*                 *                       *                *
*                                                          *
*         VERSION 2.5.2                 2016-12-10         *
*                                                          *
*    The MadGraph5_aMC@NLO Development Team - Find us at   *
*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
*                                                          *
*               Type 'help' for in-line help.              *
*                                                          *
************************************************************
INFO: load configuration from /scratch/cms/Datos/Trash/ejemplo/Cards/me5_configuration.txt  
INFO: load configuration from /home/camilo/HEPTools/MadGraph/MG5_aMC_v2_5_2-root6/input/mg5_configuration.txt  
INFO: load configuration from /scratch/cms/Datos/Trash/ejemplo/Cards/me5_configuration.txt  
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/--------------------------------------------------------------------------------------------------\
|  1. Choose the shower/hadronization program:                           shower = OFF              |
|  2. Choose the detector simulation program:                          detector = OFF              |
|  3. Run an analysis package on the events generated:                 analysis = No analysis tool interfaced to MG5aMC.|
|  4. Decay particles with the MadSpin module:                          madspin = OFF              |
|  5. Add weights to events for different model hypothesis:            reweight = OFF              |
\--------------------------------------------------------------------------------------------------/
  Either type the switch number (1 to 5) to change its setting,
  Set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
  Type 'help' for the list of all valid option
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, auto, done, delphes, delphes=ON, delphes=OFF, shower=OFF, shower=PYTHIA8, shower, detector=OFF, ... ][60s to answer] 
>


donde se pone 1 para activar la opción de shower, así:
The following switches determine which programs are run:
/--------------------------------------------------------------------------------------------------\
|  1. Choose the shower/hadronization program:                           shower = PYTHIA8          |
|  2. Choose the detector simulation program:                          detector = OFF              |
|  3. Run an analysis package on the events generated:                 analysis = No analysis tool interfaced to MG5aMC.|
|  4. Decay particles with the MadSpin module:                          madspin = OFF              |
|  5. Add weights to events for different model hypothesis:            reweight = OFF              |
\--------------------------------------------------------------------------------------------------/
  Either type the switch number (1 to 5) to change its setting,
  Set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
  Type 'help' for the list of all valid option
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, auto, done, delphes, delphes=ON, delphes=OFF, shower=OFF, shower=PYTHIA8, shower, detector=OFF, ... ][60s to answer] 
>

finalizando con la cantidad de enters necesarios para que el programa se ejecute. Pythia entrega como resultado un archivo comprimido llamado (por defecto) "tag_1_pythia8_events.hepmc.gz"

Usando una targeta de configuración:
Dentro de la carpeta "InclusiveTTbar1J/Cards/" está la tarjeta por defecto de pythia8 llamada "pythia8_card_default.dat", si se copia esta tarjeta con el nombre de "pythia8_card.dat" (con los cambios que se necesiten) en este mismo directorio y se ejecuta Madgraph, la opción de usar Pythia8 estará por defecto. (muy útil cuando se hace scripting)


Usando MadEvent?
Luego de haber generado la carpeta de salida, se procede a dicha carpeta y se llama a madevent
./bin/madevent

luego que se abra la linea de comandos del madevent, se pone lo siguiente:
ExclusiveTTbar>pythia8 run_##
>0
>0

donde ## es el número del run que se quiere hadronizar con Pythia.
previous 4 of 6 next