| ... | ... | @@ -12,7 +12,7 @@ MadGraph needs some set of cards to run. |
|
|
|
- **pythia8_card.dat** : Shower settings in Pythia8 (e.g. jet merging scales, MPI, ...)
|
|
|
|
- **madspin_card.dat** : Define the resonance decays for special cases (loop-induced or NLO processes)
|
|
|
|
|
|
|
|
Please checkout the launchpad [link](https://launchpad.net/mg5amcnlo) for more MadGraph specific details/instructions. To run Contur, **proc_card.dat** is necessary while others are only for (relatively) expert use cases. Here, we will first demonstrate how to run Contur steered MadGraph with **run_card.dat** and then explain the use case of **madspin_card.dat** as well.
|
|
|
|
Please checkout the launchpad [link](https://launchpad.net/mg5amcnlo) for more MadGraph specific details or instructions. To run Contur, **proc_card.dat** is necessary while others are only for (relatively) expert use cases. Here, we will first demonstrate how to run Contur steered MadGraph with **run_card.dat** and then explain the use case of **madspin_card.dat** as well.
|
|
|
|
|
|
|
|
Execute the commands below.
|
|
|
|
```
|
| ... | ... | @@ -73,39 +73,54 @@ Run Contur using command below. |
|
|
|
contur-batch -m madgraph -p zprime_mumu.dat -t zprime_mumu.sh -n 10000 -b 7TeV,8TeV,13TeV -r RunInfo
|
|
|
|
```
|
|
|
|
|
|
|
|
Result would probably look like below (a bit dependent on Contur release version one uses). We've excluded Z' mass up to 2000 GeV (for kl=1.0)!
|
|
|
|
[dominantPools0CLs.pdf](uploads/67f88c0444cbe12c89874b9ebd18336f/dominantPools0CLs.pdf)
|
|
|
|
Result would probably look like [dominantPools0CLs.pdf](uploads/67f88c0444cbe12c89874b9ebd18336f/dominantPools0CLs.pdf) (a bit dependent on Contur release version one uses).
|
|
|
|
We've excluded Z' mass up to 2000 GeV (for kl=1.0) using Contur!
|
|
|
|
|
|
|
|
|
|
|
|
# Relatively less easier example using `run_card.dat` and `madspin_card.dat`
|
|
|
|
|
|
|
|
MadSpin is used for resonance decays and we will not further explain what it is. For more details, please checkout MadGraph tutorials. Here, we will be using it for top quark decays which is decayed from Z'.
|
|
|
|
|
|
|
|
!! Provide other card examples (pythia8_card.dat, madspin_card.dat)
|
|
|
|
Execute the commands below.
|
|
|
|
```
|
|
|
|
source setupContur.sh
|
|
|
|
mkdir -p tutorial_with_madspin/RunInfo
|
|
|
|
cp -r data/Models/WZPrime/VPrime_NLO/ tutorial_with_madspin/RunInfo/
|
|
|
|
cd tutorial_with_madspin
|
|
|
|
cp -r RunInfo/VPrime_NLO/RunInfo_withMadSpin/*dat RunInfo/
|
|
|
|
cp RunInfo/VPrime_NLO/zprime_tt.* ./
|
|
|
|
```
|
|
|
|
|
|
|
|
Noticed `dat` and `sh` files have changed to `zprime_tt`? In `zprime_tt.sh`, only `generate` command has changed from the example above to generate top pair from Z' via s-channel in NLO.
|
|
|
|
```
|
|
|
|
import model <RunInfo/VPrime_NLO>
|
|
|
|
generate p p > zp > t t~ [QCD] # Z' to top pair via s-channel in NLO
|
|
|
|
output mgevents
|
|
|
|
launch
|
|
|
|
shower=PYTHIA8
|
|
|
|
set mzp {mzp}
|
|
|
|
set kl {kl}
|
|
|
|
set wzp Auto
|
|
|
|
```
|
|
|
|
|
|
|
|
[Run]
|
|
|
|
generator ='</Path/To/MadGraph/>/bin/mg5_aMC'
|
|
|
|
contur ='</Path/To/Contur/>/setupContur.sh'
|
|
|
|
environment = '</Path/To/Env/>/setup.sh'
|
|
|
|
Parameter scanning range has also changed as it's a different process. Take a look at `zprime_tt.dat` yourself. There is a new file named `madspin_card.dat` in `RunInfo` which defines the top quark decay modes. Let's take a look.
|
|
|
|
|
|
|
|
[Parameters]
|
|
|
|
[[mn1]]
|
|
|
|
mode = LOG
|
|
|
|
start = 10.0
|
|
|
|
stop = 1000.0
|
|
|
|
number = 13
|
|
|
|
[[vmun1]]
|
|
|
|
mode = LOG
|
|
|
|
start = 0.01
|
|
|
|
stop = 1.0
|
|
|
|
number = 13
|
|
|
|
[[sqvmun1]]
|
|
|
|
mode = REL
|
|
|
|
form = {vmun1}*{vmun1}
|
|
|
|
```
|
|
|
|
set Nevents_for_max_weight 250 # check MadGraph tutorial for details
|
|
|
|
set max_weight_ps_point 400 # check MadGraph tutorial for details
|
|
|
|
set max_running_process 1 # check MadGraph tutorial for details
|
|
|
|
|
|
|
|
define wdecay = u d c s u~ d~ c~ s~ e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~ # defining the particle container, wdecay now stands for all particles that can decay from W boson
|
|
|
|
|
|
|
|
# Running MadGraph with Contur steering commands
|
|
|
|
The set of MadGraph cards, run_card.dat, pythia8_card.dat, and madspin_card.dat (if one doesn't want to use default settings from MadGraph) should be put into RunInfo directory.
|
|
|
|
decay t > w+ b, w+ > wdecay wdecay # top quark decaying to all possible decay modes
|
|
|
|
decay t~ > w- b~, w- > wdecay wdecay # antitop quark decaying to all possible decay modes
|
|
|
|
|
|
|
|
launch # ! do not change !
|
|
|
|
```
|
|
|
|
contur-batch -r RunInfo -m madgraph -t proc_card.dat
|
|
|
|
|
|
|
|
Run Contur using command below.
|
|
|
|
```
|
|
|
|
contur-batch -m madgraph -p zprime_tt.dat -t zprime_tt.sh -n 10000 -b 7TeV,8TeV,13TeV -r RunInfo
|
|
|
|
```
|
|
|
|
|
|
|
|
Obviously, result would have changed as we are scanning totally different process. One peculiar point we can find is that low mass (where mass of Z' is smaller than 300 GeV) totally lost sensitivity. |
|
|
\ No newline at end of file |
