Previous page


This page contain few instructions of how to generate varies MC samples

On this page:

Exclusive production (via proton or pomeron)


Madgraph is able to simulate both exclusive and inclusive γγ→XX processes (where either proton remain intact or dissociate). Here is an recommended settings for simulating MG exclusive events:

#common settings
set automatic_html_opening False

#set scale (new from 2.7.1)
set fixed_fac_scale True
set dsqrt_q2fact1 2
set dsqrt_q2fact2 2

# use elastic photons (will use EPA instead of PDF)
set lpp1 2  # proton 
set lpp2 2  # proton 

#if using inelastic photons, set up the PDF to LUXqed17_plus_PDF4LHC15_nnlo_100
#set lhaid 82200

Calculation of processes that include box diagrams at QED are available loop_qcd_qed_sm-full model, for example CEP of H500→tt can be calculated using:

import model loop_qcd_qed_sm-full
generate a a > h > t t~ [noborn=QED]
output cep_ttbar
launch cep_ttbar
set param_card mass 25 500.0
set param_card decay 25 68.0

Producing gridpack

Here is the short list of instructions on how to generate gridpacks with CMSSW. See QuickGuideMadGraph5aMCatNLO for more details about production with Madgraph.

Setup the code:

git clone
cd genproductions/bin/MadGraph5_aMCatNLO/
# e.g. ./ wplustest_4f_LO cards/examples/wplustest_4f_LO 

Madgraph generates collinear photons, and the protons are not stored in the LHE event info. The workaround is to add proton manually (using code) which will be later used by proton direct simulation, which is done .

The following steps are used to generate GEN event files:

Read LHE content

The LHE files can be explored via the LHEEventProduct objects:

Events->Scan("LHE.hepeup().IDUP : LHE.hepeup().ISTUP");

-- MichaelPitt - 2020-04-09

Edit | Attach | Watch | Print version | History: r7 < r6 < r5 < r4 < r3 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r7 - 2022-01-17 - MichaelPitt
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback