Commit 30824b86 authored by Nishu Nishu's avatar Nishu Nishu
Browse files

2HDM+a signal for a 4top signature (Dsids 510250-510314)

parent 80a3d8ab
Pipeline #3469934 passed with stages
in 53 seconds
import MadGraphControl.MadGraph_NNPDF30NLO_Base_Fragment
from MadGraphControl.MadGraphUtils import *
import re
# to check lepton filter
import os
### get MC job-options filename
FIRST_DIR = (os.environ['JOBOPTSEARCHPATH']).split(":")[0]
jofiles = [f for f in os.listdir(FIRST_DIR) if (f.startswith('mc') and f.endswith('.py'))]
filter_string=jofiles[0].split('_')[-1].replace(".py","")
#----------------------------------------------------------------------------
# Process definition
#----------------------------------------------------------------------------
initialGluons=False
#
# Specical syntax: see the 2HDM+a white paper (https://arxiv.org/pdf/1810.09420.pdf)
# QCD<=2: Four-top production from NP processes, including interference terms amongH,A,a
# QED<=2: Four-top production including both SM and NP contributions and their interference.
# QED<=0: Four-top production within the SM.
#
process="""
import model Pseudoscalar_2HDM -modelname
define p = g d u s c b d~ u~ s~ c~ b~
define j = g d u s c b d~ u~ s~ c~ b~
generate p p > t t~ t t~ / a z h1 QCD<=2
output -f
"""
#----------------------------------------------------------------------------
# Beam energy
#----------------------------------------------------------------------------
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
#---------------------------------------------------------------------------
# Number of events to generate
#---------------------------------------------------------------------------
nevents=evgenConfig.nEventsPerJob
if LHE_EventMultiplier>0:
nevents=runArgs.maxEvents*LHE_EventMultiplier
#---------------------------------------------------------------------------
# MG5 Run Card
#---------------------------------------------------------------------------
extras = {}
if initialGluons:
extras = {
'maxjetflavor' : 4,
'asrwgtflavor' : 4,
'lhe_version' : '3.0',
'cut_decays' : 'F',
'nevents' : nevents,
}
else:
extras = {
'maxjetflavor' : 5,
'asrwgtflavor' : 5,
'lhe_version' : '3.0',
'cut_decays' : 'F',
'nevents' : nevents,
}
# Build run_card
process_dir = new_process(process)
modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=extras)
#---------------------------------------------------------------------------
# MG5 parameter Card
#---------------------------------------------------------------------------
# Build param_card.dat
params={}
## blocks might be modified
dict_blocks={
"mass": ["MB", "MXd", "mh2", "mh3", "mhc", "mh4"],
"DMINPUTS" : ["gPXd", ],
"FRBlock": ["tanbeta", "sinbma", ],
"Higgs": ["lam3", "laP1", "laP2", "sinp"],
}
for bl in dict_blocks.keys():
for pa in dict_blocks[bl]:
if pa in THDMparams.keys():
if bl not in params: params[bl]={}
if pa=="MB":
params[bl]["5"]=THDMparams[pa]
else:
params[bl][pa]=THDMparams[pa]
## auto calculation of decay width
THDMparams_decay={
"25": "Auto",
"35": "Auto",
"36": "Auto",
"37": "Auto",
"55": "Auto",
}
params["decay"]=THDMparams_decay
print("Updating parameters:")
print(params)
modify_param_card(process_dir=process_dir,params=params)
# Build reweight_card.dat
if reweight:
# Create reweighting card
reweight_card_loc=process_dir+'/Cards/reweight_card.dat'
rwcard = open(reweight_card_loc,'w')
for rw_name in reweights:
params_rwt = params.copy()
for param in rw_name.split('-'):
param_name, value = param.split('_')
if param_name == "SINP":
params_rwt['HIGGS']['SINP'] = value
elif param_name == "TANB":
params_rwt['FRBLOCK']['TANBETA'] = value
param_card_reweight = process_dir+'/Cards/param_card_reweight.dat'
#shutil.copy(process_dir+'/Cards/param_card_default.dat', param_card_reweight)
shutil.copy(process_dir+'/Cards/param_card.dat', param_card_reweight)
param_card_rwt_new=process_dir+'/Cards/param_card_rwt_%s.dat' % rw_name
modify_param_card(param_card_input=param_card_reweight, process_dir=process_dir,params=params_rwt, output_location=param_card_rwt_new)
rwcard.write("launch --rwgt_info=%s\n" % rw_name)
rwcard.write("%s\n" % param_card_rwt_new)
rwcard.close()
print_cards()
# Build MadSpin card
madspin_card=process_dir+'/Cards/madspin_card.dat'
if os.access(madspin_card,os.R_OK):
os.unlink(madspin_card)
mscard = open(madspin_card,'w')
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
# specify the decay for the final state particles
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
decay t > w+ b
decay t~ > w- b~
decay w+ > all all
decay w- > all all
# running the actual code
launch"""%runArgs.randomSeed)
mscard.close()
#---------------------------------------------------------------------------
# Generate the events
#---------------------------------------------------------------------------
generate(process_dir=process_dir,runArgs=runArgs)
arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=True)
#---------------------------------------------------------------------------
# Metadata
#---------------------------------------------------------------------------
if initialGluons:
evgenConfig.process = "p p > t t~ t t~ / a z h1 QCD<=2"
initialStateString = "gluon fusion"
evgenConfig.description = "Pseudoscalar 2HDMa model for 4-top events"
evgenConfig.keywords = ["exotic","BSM","WIMP"]
evgenConfig.contact = ["Johanna Gramling <jgramlin@cern.ch>", "Lailin Xu <lailin.xu@cern.ch>"]
#---------------------------------------------------------------------------
# Shower
#---------------------------------------------------------------------------
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_MadGraph.py")
# Teach pythia about the dark matter particle
genSeq.Pythia8.Commands += ["SLHA:allowUserOverride = on",
"1000022:all = xd xd~ 2 0 0 %d 0.0 0.0 0.0 0.0" % (int(THDMparams['MXd'])),
"1000022:isVisible = false",
"1000022:mayDecay = off"
]
#particle data = name antiname spin=2s+1 3xcharge colour mass width (left out, so set to 0: mMin mMax tau0)
genSeq.Pythia8.Commands += ["1000022:all = xd xd~ 2 0 0 %d 0.0 0.0 0.0 0.0" % (int(THDMparams['MXd'])),
"1000022:isVisible = false"]
#---------------------------------------------------------------------------
# Filters
#---------------------------------------------------------------------------
# Generator filter
if (filter_string == "1L"):
evgenLog.info('1lepton filter applied')
include('GeneratorFilters/LeptonFilter.py' )
filtSeq.LeptonFilter.Ptcut = 20000.
filtSeq.LeptonFilter.Etacut = 2.8
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 200 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 200 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 200 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 200 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 400 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 400 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 400 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 200 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 500 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 500 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 500 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 200 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.50 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 600.0 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 600.0 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 600.0 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 200 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.10-TANB_1.0',
'SINP_0.20-TANB_1.0',
'SINP_0.30-TANB_1.0',
'SINP_0.40-TANB_1.0',
'SINP_0.60-TANB_1.0',
'SINP_0.70-TANB_1.0',
'SINP_0.80-TANB_1.0',
'SINP_0.90-TANB_1.0',
'SINP_0.10-TANB_0.5',
'SINP_0.20-TANB_0.5',
'SINP_0.30-TANB_0.5',
'SINP_0.40-TANB_0.5',
'SINP_0.50-TANB_0.5',
'SINP_0.60-TANB_0.5',
'SINP_0.70-TANB_0.5',
'SINP_0.80-TANB_0.5',
'SINP_0.90-TANB_0.5',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 200 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 200 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 200 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 250 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 400 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 400 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 400 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 250 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 500 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 500 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 500 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 250 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")
../../510xxx/510250/MadGraphControl_Py8G_2HDMa_4t_MadSpin.py
\ No newline at end of file
THDMparams = {}
THDMparams['gPXd'] = 1.0 # The coupling of the additional pseudoscalar mediator to dark matter (DM). This coupling is called $y_\chi$ in (2.5) of arXiv:1701.07427.
THDMparams['tanbeta'] = 1.0 # The ratio of the vacuum expectation values $ an eta = v_2/v_1$ of the Higgs doublets $H_2$ and $H_1$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['sinbma'] = 1.0 # The sine of the difference of the mixing angles $\sin (eta - lpha)$ in the scalar potential containing only the Higgs doublets. This quantity is defined in Section 3.1 of arXiv:1701.07427.
THDMparams['lam3'] = 3.0 # The quartic coupling of the scalar doublets $H_1$ and $H_2$. This parameter corresponds to the coefficient $\lambda_3$ in (2.1) of arXiv:1701.07427.
THDMparams['laP1'] = 3.0 # The quartic coupling between the scalar doublets $H_1$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P1}$ in (2.2) of arXiv:1701.07427.
THDMparams['laP2'] = 3.0 # The quartic coupling between the scalar doublets $H_2$ and the pseudoscalar $P$. This parameter corresponds to the coefficient $\lambda_{P2}$ in (2.2) of arXiv:1701.07427.
THDMparams['sinp'] = 0.35 # The sine of the mixing angle $ heta$, as defined in Section 2.1 of arXiv:1701.07427.
THDMparams['MXd'] = 10.0 # The mass of the fermionic DM candidate denoted by $m_\chi$ in arXiv:1701.07427.
THDMparams['mh1'] = 125. # The mass of the lightest scalar mass eigenstate $h$, which is identified in arXiv:1701.07427 with the Higgs-like resonance found at the LHC.
THDMparams['mh2'] = 600 # The mass of the heavy scalar mass eigenstate $H$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh3'] = 600 # The mass of the heavy pseudoscalar mass eigenstate $A$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mhc'] = 600 # The mass of the charged scalar eigenstate $H^\pm$. See Section 2.1 of arXiv:1701.07427 for further details.
THDMparams['mh4'] = 250 # The mass of the pseudoscalar mass eigenstate $a$ that decouples for $\sin heta = 0$. See Section 2.1 of arXiv:1701.07427 for further details.
initialGluons = True # Determines initial state to generate from. True --> top loop induced production (gluon Fusion). False --> b b~ -initiated production (b-anti-b annihilation). See Sections 5.4 and 6.5 of arXiv:1701.07427 for further details
reweight = True # Determines whether to store alternative weights for tanb and sinp
reweights=[
'SINP_0.7-TANB_0.3',
'SINP_0.7-TANB_0.5',
'SINP_0.7-TANB_1.0',
'SINP_0.7-TANB_3.0',
'SINP_0.7-TANB_5.0',
'SINP_0.35-TANB_0.3',
'SINP_0.35-TANB_0.5',
'SINP_0.35-TANB_3.0',
'SINP_0.35-TANB_5.0',
]
# Number of events to generate
evgenConfig.nEventsPerJob=1000
# Event multiplier
LHE_EventMultiplier = 5 #For Filter, Multiplier on maxEvents
include("MadGraphControl_Py8G_2HDMa_4t_MadSpin.py")