Skip to content
Snippets Groups Projects
README.md 4.96 KiB
Newer Older
Yu-Heng Chen's avatar
Yu-Heng Chen committed
# 2HDM A/H > tt~ (+Interf) event generation/reweighting
The code is used to produce the signal-plus-interference pattern of BSM Higgs produced in the top-antitop quark pair final states in either a generic 2HDM or 2HDM+a model.
Yu-Heng Chen's avatar
Yu-Heng Chen committed

Yu-Heng Chen's avatar
Yu-Heng Chen committed
Prerequisites
-------------
#### Necessary
*   MadGraph >= 2.6.0
*   AanalysisTop,21.2.87 or similar
    - SLC6 environment is recommended to avoid conflicts due to the python2.7->3.0 transition.
*   HTCondor or similar batch systems
*   UFO model files:
    -    Generic 2HDM: ![Higgs_Effective_Couplings_FormFact]
    -    2HDM+a: ![Pseudoscalar_2HDM_FormFactor]

[Higgs_Effective_Couplings_FormFact]: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Models/ggHFullLoop
[Pseudoscalar_2HDM_FormFactor]: https://gitlab.cern.ch/yuchen/Pseudoscalar_2HDM_FormFactor
Yu-Heng Chen's avatar
Yu-Heng Chen committed

Yu-Heng Chen's avatar
Yu-Heng Chen committed
#### Optional
*   IPython (for debugging)
*   coloredlogs (for debugging)

Setup
-----
1.  setupATLAS
    ```bash
    setupATLAS
    ```
Yu-Heng Chen's avatar
Yu-Heng Chen committed
2.  `git clone` this package
    ```bash
Yizhou Cai's avatar
Yizhou Cai committed
    git clone --recurse-submodules https://gitlab.cern.ch/caiyi/2HDM_SI.git
    git checkout MS_new
Yu-Heng Chen's avatar
Yu-Heng Chen committed
    ```
3.  Modify `config.py` to adapt to your needs
4.  source `setup.sh` every time whenever you log in.
Yu-Heng Chen's avatar
Yu-Heng Chen committed
    ```bash
    source ./setup.sh
    ```

Usage
-----
`./main.py --help` to show this help message:
Yu-Heng Chen's avatar
Yu-Heng Chen committed
```bash
Yu-Heng Chen's avatar
Yu-Heng Chen committed
usage: main.py [-h] [-i ISEED] [-e ECM] [--param_card PARAM_CARD]
               [-p [sim_opt [sim_opt ...]]] [-pv [sim_opt [sim_opt ...]]]
               [-t [sim_opt [sim_opt ...]]] [-tv [sim_opt [sim_opt ...]]]
               [-c CLUSTER] [--ffmt FFMT] [--model {2HDM,2HDM+a}]
Yizhou Cai's avatar
Yizhou Cai committed
               [-s SYSTEMATICS]
               {parameters} ...
Yu-Heng Chen's avatar
Yu-Heng Chen committed

Yu-Heng Chen's avatar
Yu-Heng Chen committed
MadgraphPythia Type-II 2HDM gg > H/A > tt~ (+ Interf) production.

positional arguments:
Yizhou Cai's avatar
Yizhou Cai committed
  {parameters}
Yu-Heng Chen's avatar
Yu-Heng Chen committed
    parameters          Specify a process using physical parameters.

optional arguments:
  -h, --help            show this help message and exit
  -i ISEED, --iseed ISEED
Yu-Heng Chen's avatar
Yu-Heng Chen committed
                        Random seed number (default: 21)
  -e ECM, --ecm ECM     Central Mass Energy sqrt(s). (default: None)
Yu-Heng Chen's avatar
Yu-Heng Chen committed
  --param_card PARAM_CARD
                        Instead of directly using the parameter set from run
                        number/parameters, import a MG5 parameter card. Only
Yizhou Cai's avatar
Yizhou Cai committed
                        used when generating event (default: ) (Suggested to be auto)
Yu-Heng Chen's avatar
Yu-Heng Chen committed
  -p [sim_opt [sim_opt ...]], --parton_level [sim_opt [sim_opt ...]]
Yizhou Cai's avatar
Yizhou Cai committed
                        Parton-level simulation. Powered by MG5_aMC@NLOv2.6.7
Yu-Heng Chen's avatar
Yu-Heng Chen committed
                        (default: None)
  -pv [sim_opt [sim_opt ...]], --parton_level_validation [sim_opt [sim_opt ...]]
                        Produce a ntuple that has all the neccessary parton-
                        level information inside. The Matrix Element will be
                        calculated as well. If <FFMT> is given, then the file
                        matching will be done and stored into the ntuples.
                        (default: None)
  -t [sim_opt [sim_opt ...]], --truth_level [sim_opt [sim_opt ...]]
                        Truth-level simulation. Powered by Pythia6 (default:
                        None)
  -tv [sim_opt [sim_opt ...]], --truth_level_validation [sim_opt [sim_opt ...]]
                        Produce a ntuple that has all the neccessary truth-
                        level information. RWGT will be done at the procedure.
                        (default: None)
  -c CLUSTER, --cluster CLUSTER
                        Batch system to use. (default: condor)
  --ffmt FFMT           File format that will be produced or matched (default:
                        DAOD_TRUTH1)
  --model {2HDM,2HDM+a}
                        Model to run. (default: None)
Yizhou Cai's avatar
Yizhou Cai committed
  -s SYSTEMATICS, --systematics SYSTEMATICS
                        Caculate systematics uncertainty (default: False)
Yu-Heng Chen's avatar
Yu-Heng Chen committed
```
Yu-Heng Chen's avatar
Yu-Heng Chen committed

Yu-Heng Chen's avatar
Yu-Heng Chen committed
MC Generation
-------------
Yu-Heng Chen's avatar
Yu-Heng Chen committed
LHE/DAOD_TRUTH1 Production. __You may also need this for reweighting.__
Yu-Heng Chen's avatar
Yu-Heng Chen committed
1.  For example,
    ```bash
Yizhou Cai's avatar
Yizhou Cai committed
    ./main.py parameters A 500 0.4 SI "[:1000000:10000]" 100 0 --model 2HDM+a -e 13000 -p (or -t or -pv or -tv) --param_card auto
Yu-Heng Chen's avatar
Yu-Heng Chen committed
    ```
Yizhou Cai's avatar
Yizhou Cai committed
    *   The `parameters` is used to control the model parameter, following with mediator, mass of heavy Higgs, tan\beta, signal type, slicer, mass of a, sin\theta in order.
Yu-Heng Chen's avatar
Yu-Heng Chen committed
    *   The slicer is used to control the sample size and number of events per sub-LHE file. `[:1000000:10000]` means one million events are produced and ten thousand events are stored in each sub-LHE file.
Yizhou Cai's avatar
Yizhou Cai committed
    *   With `-p`, LHE files of ME-level events are produced. With `-s`, systematics uncertainty will be considered. Add `--param_card [File-Path]` if you want to specify a parameter card.  
Yu-Heng Chen's avatar
Yu-Heng Chen committed
        with `-pv`, flat ntuples of ME-level events with reweighting-related information included are produced.  
        with `-t`, DAOD_TRUTH1 files of events after showering are produced.  
        with `-tv`, flat ntuples of events after showering with reweighting information (and matching information if any) included are produced.
Yizhou Cai's avatar
Yizhou Cai committed
    
Yu-Heng Chen's avatar
Yu-Heng Chen committed
2.   Please remember the slicing information and use the same throughout the whole event generation/reweighting procedure.