Commit 56537ca7 authored by Joany Manjarres's avatar Joany Manjarres
Browse files

See ChangeLog

parent 9091c5a5
2021-03-18 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-04-51
* add JO 451908-922 300471-2
2021-03-04 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-04-50
* add susy.399*slha
......
##############################################################
# Python snippet to generate EvtGen user decay file on the fly
# B+ -> J/psi(mu0.0mu0.0) K+: an analog of MC15.300404
# w/o cuts for muons and kaon
##############################################################
f = open("BU_JPSI_K_USER.DEC","w")
f.write("Particle B+ 5.27929 0.0\n") # PDG2014 mass
f.write("Particle B- 5.27929 0.0\n") # PDG2014 mass
f.write("Define dm_incohMix_B0 0.0\n") #disable neutral meson mixing
f.write("Define dm_incohMix_B_s0 0.0\n") #same, LEAVE AS B_s0
f.write("Alias myJ/psi J/psi\n")
f.write("Decay B+\n")
f.write("1.0000 myJ/psi K+ SVS;\n")
f.write("Enddecay\n")
f.write("Decay myJ/psi\n")
f.write("1.0000 mu+ mu- VLL;\n")
f.write("Enddecay\n")
f.write("End\n")
f.close()
##############################################################
evgenConfig.description = "Exclusive Bplus->Jpsi_mu3p5mu3p5_Kplus production"
evgenConfig.keywords = ["exclusive","Bplus","Jpsi","2muon"]
evgenConfig.minevents = 2000
include("MC15JobOptions/nonStandard/Pythia8B_A14_CTEQ6L1_EvtGen_Common.py")
### Put the content of MC15JobOptions/Pythia8B_exclusiveB_Common.py
### except actual closing B decays
# Set B+ mass in Pythia equal to that in EvtGen
genSeq.Pythia8B.Commands += ['521:m0 = 5.27929'] # PDG2014 mass
# Hard process
genSeq.Pythia8B.Commands += ['HardQCD:all = on'] # Equivalent of MSEL1
genSeq.Pythia8B.Commands += ['ParticleDecays:mixB = off']
genSeq.Pythia8B.Commands += ['HadronLevel:all = off']
# Event selection
genSeq.Pythia8B.SelectBQuarks = True
genSeq.Pythia8B.SelectCQuarks = False
genSeq.Pythia8B.VetoDoubleBEvents = True
# List of B-species
include("MC15JobOptions/Pythia8B_BPDGCodes.py")
###
###
#include("MC15JobOptions/BSignalFilter.py")
genSeq.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 7.']
genSeq.Pythia8B.QuarkPtCut = 0.0
genSeq.Pythia8B.AntiQuarkPtCut = 7.0
genSeq.Pythia8B.QuarkEtaCut = 102.5
genSeq.Pythia8B.AntiQuarkEtaCut = 2.6
genSeq.Pythia8B.RequireBothQuarksPassCuts = True
genSeq.Pythia8B.NHadronizationLoops = 1
# Make Pythia select the events with B+
genSeq.Pythia8B.SignalPDGCodes = [ 521 ]
genSeq.EvtInclusiveDecay.userDecayFile = "BU_JPSI_K_USER.DEC"
include("MC15JobOptions/ParentChildwStatusFilter.py")
filtSeq.ParentChildwStatusFilter.PDGParent = [521]
filtSeq.ParentChildwStatusFilter.StatusParent = [2]
filtSeq.ParentChildwStatusFilter.PtMinParent = 8900.
filtSeq.ParentChildwStatusFilter.EtaRangeParent = 2.5
filtSeq.ParentChildwStatusFilter.PDGChild = [443]
filtSeq.ParentChildwStatusFilter.PtMinChild = 0.
filtSeq.ParentChildwStatusFilter.EtaRangeChild = 999999.
# Final state selections
#filtSeq.BSignalFilter.LVL1MuonCutOn = True
#filtSeq.BSignalFilter.LVL2MuonCutOn = True
#filtSeq.BSignalFilter.LVL1MuonCutPT = 0
#filtSeq.BSignalFilter.LVL1MuonCutEta = 999999.
#filtSeq.BSignalFilter.LVL2MuonCutPT = 0
#filtSeq.BSignalFilter.LVL2MuonCutEta = 999999.
#filtSeq.BSignalFilter.B_PDGCode = 521
#filtSeq.BSignalFilter.Cuts_Final_hadrons_switch = True
#filtSeq.BSignalFilter.Cuts_Final_hadrons_pT = 0.0
#filtSeq.BSignalFilter.Cuts_Final_hadrons_eta = 999999.
##############################################################
# Python snippet to generate EvtGen user decay file on the fly
# B- -> J/psi(mu0.0mu0.0) K-: an analog of MC15.300405
# w/o cuts for muons and kaon
##############################################################
f = open("BU_JPSI_K_USER.DEC","w")
f.write("Particle B+ 5.27929 0.0\n") # PDG2014 mass
f.write("Particle B- 5.27929 0.0\n") # PDG2014 mass
f.write("Define dm_incohMix_B0 0.0\n") #disable neutral meson mixing
f.write("Define dm_incohMix_B_s0 0.0\n") #same, LEAVE AS B_s0
f.write("Alias myJ/psi J/psi\n")
f.write("Decay B-\n")
f.write("1.0000 myJ/psi K- SVS;\n")
f.write("Enddecay\n")
f.write("Decay myJ/psi\n")
f.write("1.0000 mu+ mu- VLL;\n")
f.write("Enddecay\n")
f.write("End\n")
f.close()
##############################################################
evgenConfig.description = "Exclusive Bminus->Jpsi_mu0p0mu0p0_Kminus production"
evgenConfig.keywords = ["exclusive","Bminus","Jpsi","2muon"]
evgenConfig.minevents = 2000
include("MC15JobOptions/nonStandard/Pythia8B_A14_CTEQ6L1_EvtGen_Common.py")
### Put the content of MC15JobOptions/Pythia8B_exclusiveB_Common.py
### except actual closing B decays
# Set B+ mass in Pythia equal to that in EvtGen
genSeq.Pythia8B.Commands += ['521:m0 = 5.27929'] # PDG2014 mass
# Hard process
genSeq.Pythia8B.Commands += ['HardQCD:all = on'] # Equivalent of MSEL1
genSeq.Pythia8B.Commands += ['ParticleDecays:mixB = off']
genSeq.Pythia8B.Commands += ['HadronLevel:all = off']
# Event selection
genSeq.Pythia8B.SelectBQuarks = True
genSeq.Pythia8B.SelectCQuarks = False
genSeq.Pythia8B.VetoDoubleBEvents = True
# List of B-species
include("MC15JobOptions/Pythia8B_BPDGCodes.py")
###
###
#include("MC15JobOptions/BSignalFilter.py")
genSeq.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 7.']
genSeq.Pythia8B.QuarkPtCut = 7.0
genSeq.Pythia8B.AntiQuarkPtCut = 0.0
genSeq.Pythia8B.QuarkEtaCut = 2.6
genSeq.Pythia8B.AntiQuarkEtaCut = 102.5
genSeq.Pythia8B.RequireBothQuarksPassCuts = True
genSeq.Pythia8B.NHadronizationLoops = 1
# Make Pythia select the events with B+
genSeq.Pythia8B.SignalPDGCodes = [ -521 ]
genSeq.EvtInclusiveDecay.userDecayFile = "BU_JPSI_K_USER.DEC"
include("MC15JobOptions/ParentChildwStatusFilter.py")
filtSeq.ParentChildwStatusFilter.PDGParent = [521]
filtSeq.ParentChildwStatusFilter.StatusParent = [2]
filtSeq.ParentChildwStatusFilter.PtMinParent = 8900.
filtSeq.ParentChildwStatusFilter.EtaRangeParent = 2.5
filtSeq.ParentChildwStatusFilter.PDGChild = [443]
filtSeq.ParentChildwStatusFilter.PtMinChild = 0.
filtSeq.ParentChildwStatusFilter.EtaRangeChild = 999999.
# Final state selections
#filtSeq.BSignalFilter.LVL1MuonCutOn = True
#filtSeq.BSignalFilter.LVL2MuonCutOn = True
#filtSeq.BSignalFilter.LVL1MuonCutPT = 0
#filtSeq.BSignalFilter.LVL1MuonCutEta = 999999.
#filtSeq.BSignalFilter.LVL2MuonCutPT = 0
#filtSeq.BSignalFilter.LVL2MuonCutEta = 999999.
#filtSeq.BSignalFilter.B_PDGCode = -521
#filtSeq.BSignalFilter.Cuts_Final_hadrons_switch = True
#filtSeq.BSignalFilter.Cuts_Final_hadrons_pT = 0.0
#filtSeq.BSignalFilter.Cuts_Final_hadrons_eta = 999999.
#Author: diallo boye <diallo.boye@cern.ch>
#This is a template joboption for generating events with the HAHM madgraph model
#To use it, copy the template and rename it by replacing the XXXXXX and YYYYYYYYY in the filename with:
# make XXXXXX match the runNumber you give to Generate_tf
# make the YYYYYYYY a short description
#Then modify:
# - proc_card (for the process),
# - param_card_extras (for the model parameters),
# - run_card_extras (for the generator-level cuts),
# - post_lhe_hook (a hook method for executing code after the lhe making (madgraph) but before showering (pythia)
#Run it with (for example):
# Generate_tf.py --ecmEnergy=13000 --runNumber=302076 --firstEvent=1 --asetup="" --maxEvents=5000 --randomSeed=1 --jobConfig="MC15.302076.MadGraphPythia8EvtGen_A14NNPDF23LO_HAHMggfZdZd4l_mZd5.py" --outputEVNTFile=my.evgen.root
#--------------------------
#<CHANGE THESE SETTINGS. keywords MUST BE IN LIST OF ALLOWED KEYWORDS>
evgenConfig.description="high mass higgs to ZZ"
evgenConfig.keywords+=['exotic','BSMHiggs']
evgenConfig.contact = ['diallo.boye@cern.ch']
evgenConfig.process = "HAHM_H_ZdZd_4l"
#<REPLACE WITH YOUR PROC_CARD CONTENT>
proc_card = """
import model HAHM_variableMW_v3_UFO
define l+ = e+ mu+
define l- = e- mu-
generate g g > h HIG=1 HIW=0 QED=0 QCD=0, (h > Zp Zp, Zp > l+ l-)"""
proc_name = evgenConfig.process #just used for the mg directory name
#modifications to the param_card.dat (generated from the proc_card i.e. the specific model)
#if you want to see the resulting param_card, run Generate_tf with this jobo, and look at the param_card.dat in the cwd
#If you want to see the auto-calculated values of the decay widths, look at the one in <proc_name>/Cards/param_card.dat (again, after running a Generate_tf)
param_card_extras = { "HIDDEN": { 'epsilon': '1e-4', #kinetic mixing parameter
'kap': '1e-4', #higgs mixing parameter
'mzdinput': '15', #Zd mass
'mhsinput':'1000.0' }, #dark higgs mass
"HIGGS": { 'mhinput':'240'}, #higgs mass
"DECAY": { 'wzp':'Auto', 'wh':'3.39', 'wt':'Auto' } #auto-calculate decay widths and BR of Zp, H, t
}
run_card_extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'ptj':'0',
'ptb':'0',
'pta':'0',
'ptl':'0',
'etaj':'-1',
'etab':'-1',
'etaa':'-1',
'etal':'-1',
'drjj':'0',
'drbb':'0',
'drll':'0',
'draa':'0',
'drbj':'0',
'draj':'0',
'drjl':'0',
'drab':'0',
'drbl':'0',
'dral':'0' }
def post_lhe_hook():
#You can put code here that will get called after the 'generate' method of MadGraphControl is called
#What follows is an example for modifying the lifetime of the zd:
#The 'add_lifetimes' method is defined in the file included below. It modifies the lhe file to include particle lifetimes on the particle with given pdgId
#add_lifetimes(pdgId=32, avgtau=5)
pass
include("MC15JobOptions/MadGraphControl_Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
#Author: diallo boye <diallo.boye@cern.ch>
#This is a template joboption for generating events with the HAHM madgraph model
#To use it, copy the template and rename it by replacing the XXXXXX and YYYYYYYYY in the filename with:
# make XXXXXX match the runNumber you give to Generate_tf
# make the YYYYYYYY a short description
#Then modify:
# - proc_card (for the process),
# - param_card_extras (for the model parameters),
# - run_card_extras (for the generator-level cuts),
# - post_lhe_hook (a hook method for executing code after the lhe making (madgraph) but before showering (pythia)
#Run it with (for example):
# Generate_tf.py --ecmEnergy=13000 --runNumber=302076 --firstEvent=1 --asetup="" --maxEvents=5000 --randomSeed=1 --jobConfig="MC15.302076.MadGraphPythia8EvtGen_A14NNPDF23LO_HAHMggfZdZd4l_mZd5.py" --outputEVNTFile=my.evgen.root
#--------------------------
#<CHANGE THESE SETTINGS. keywords MUST BE IN LIST OF ALLOWED KEYWORDS>
evgenConfig.description="high mass higgs to ZZ"
evgenConfig.keywords+=['exotic','BSMHiggs']
evgenConfig.contact = ['diallo.boye@cern.ch']
evgenConfig.process = "HAHM_H_ZdZd_4l"
#<REPLACE WITH YOUR PROC_CARD CONTENT>
proc_card = """
import model HAHM_variableMW_v3_UFO
define l+ = e+ mu+
define l- = e- mu-
generate g g > h HIG=1 HIW=0 QED=0 QCD=0, (h > Zp Zp, Zp > l+ l-)"""
proc_name = evgenConfig.process #just used for the mg directory name
#modifications to the param_card.dat (generated from the proc_card i.e. the specific model)
#if you want to see the resulting param_card, run Generate_tf with this jobo, and look at the param_card.dat in the cwd
#If you want to see the auto-calculated values of the decay widths, look at the one in <proc_name>/Cards/param_card.dat (again, after running a Generate_tf)
param_card_extras = { "HIDDEN": { 'epsilon': '1e-4', #kinetic mixing parameter
'kap': '1e-4', #higgs mixing parameter
'mzdinput': '15', #Zd mass
'mhsinput':'1000.0' }, #dark higgs mass
"HIGGS": { 'mhinput':'400'}, #higgs mass
"DECAY": { 'wzp':'Auto', 'wh':'29.16', 'wt':'Auto' } #auto-calculate decay widths and BR of Zp, H, t
}
run_card_extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'ptj':'0',
'ptb':'0',
'pta':'0',
'ptl':'0',
'etaj':'-1',
'etab':'-1',
'etaa':'-1',
'etal':'-1',
'drjj':'0',
'drbb':'0',
'drll':'0',
'draa':'0',
'drbj':'0',
'draj':'0',
'drjl':'0',
'drab':'0',
'drbl':'0',
'dral':'0' }
def post_lhe_hook():
#You can put code here that will get called after the 'generate' method of MadGraphControl is called
#What follows is an example for modifying the lifetime of the zd:
#The 'add_lifetimes' method is defined in the file included below. It modifies the lhe file to include particle lifetimes on the particle with given pdgId
#add_lifetimes(pdgId=32, avgtau=5)
pass
include("MC15JobOptions/MadGraphControl_Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
#Author: diallo boye <diallo.boye@cern.ch>
#This is a template joboption for generating events with the HAHM madgraph model
#To use it, copy the template and rename it by replacing the XXXXXX and YYYYYYYYY in the filename with:
# make XXXXXX match the runNumber you give to Generate_tf
# make the YYYYYYYY a short description
#Then modify:
# - proc_card (for the process),
# - param_card_extras (for the model parameters),
# - run_card_extras (for the generator-level cuts),
# - post_lhe_hook (a hook method for executing code after the lhe making (madgraph) but before showering (pythia)
#Run it with (for example):
# Generate_tf.py --ecmEnergy=13000 --runNumber=302076 --firstEvent=1 --asetup="" --maxEvents=5000 --randomSeed=1 --jobConfig="MC15.302076.MadGraphPythia8EvtGen_A14NNPDF23LO_HAHMggfZdZd4l_mZd5.py" --outputEVNTFile=my.evgen.root
#--------------------------
#<CHANGE THESE SETTINGS. keywords MUST BE IN LIST OF ALLOWED KEYWORDS>
evgenConfig.description="high mass higgs to ZZ"
evgenConfig.keywords+=['exotic','BSMHiggs']
evgenConfig.contact = ['diallo.boye@cern.ch']
evgenConfig.process = "HAHM_H_ZdZd_4l"
#<REPLACE WITH YOUR PROC_CARD CONTENT>
proc_card = """
import model HAHM_variableMW_v3_UFO
define l+ = e+ mu+
define l- = e- mu-
generate g g > h HIG=1 HIW=0 QED=0 QCD=0, (h > Zp Zp, Zp > l+ l-)"""
proc_name = evgenConfig.process #just used for the mg directory name
#modifications to the param_card.dat (generated from the proc_card i.e. the specific model)
#if you want to see the resulting param_card, run Generate_tf with this jobo, and look at the param_card.dat in the cwd
#If you want to see the auto-calculated values of the decay widths, look at the one in <proc_name>/Cards/param_card.dat (again, after running a Generate_tf)
param_card_extras = { "HIDDEN": { 'epsilon': '1e-4', #kinetic mixing parameter
'kap': '1e-4', #higgs mixing parameter
'mzdinput': '80', #Zd mass
'mhsinput':'1000.0' }, #dark higgs mass
"HIGGS": { 'mhinput':'400'}, #higgs mass
"DECAY": { 'wzp':'Auto', 'wh':'29.16', 'wt':'Auto' } #auto-calculate decay widths and BR of Zp, H, t
}
run_card_extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'ptj':'0',
'ptb':'0',
'pta':'0',
'ptl':'0',
'etaj':'-1',
'etab':'-1',
'etaa':'-1',
'etal':'-1',
'drjj':'0',
'drbb':'0',
'drll':'0',
'draa':'0',
'drbj':'0',
'draj':'0',
'drjl':'0',
'drab':'0',
'drbl':'0',
'dral':'0' }
def post_lhe_hook():
#You can put code here that will get called after the 'generate' method of MadGraphControl is called
#What follows is an example for modifying the lifetime of the zd:
#The 'add_lifetimes' method is defined in the file included below. It modifies the lhe file to include particle lifetimes on the particle with given pdgId
#add_lifetimes(pdgId=32, avgtau=5)
pass
include("MC15JobOptions/MadGraphControl_Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
#Author: diallo boye <diallo.boye@cern.ch>
#This is a template joboption for generating events with the HAHM madgraph model
#To use it, copy the template and rename it by replacing the XXXXXX and YYYYYYYYY in the filename with:
# make XXXXXX match the runNumber you give to Generate_tf
# make the YYYYYYYY a short description
#Then modify:
# - proc_card (for the process),
# - param_card_extras (for the model parameters),
# - run_card_extras (for the generator-level cuts),
# - post_lhe_hook (a hook method for executing code after the lhe making (madgraph) but before showering (pythia)
#Run it with (for example):
# Generate_tf.py --ecmEnergy=13000 --runNumber=302076 --firstEvent=1 --asetup="" --maxEvents=5000 --randomSeed=1 --jobConfig="MC15.302076.MadGraphPythia8EvtGen_A14NNPDF23LO_HAHMggfZdZd4l_mZd5.py" --outputEVNTFile=my.evgen.root
#--------------------------
#<CHANGE THESE SETTINGS. keywords MUST BE IN LIST OF ALLOWED KEYWORDS>
evgenConfig.description="high mass higgs to ZZ"
evgenConfig.keywords+=['exotic','BSMHiggs']
evgenConfig.contact = ['diallo.boye@cern.ch']
evgenConfig.process = "HAHM_H_ZdZd_4l"
#<REPLACE WITH YOUR PROC_CARD CONTENT>
proc_card = """
import model HAHM_variableMW_v3_UFO
define l+ = e+ mu+
define l- = e- mu-
generate g g > h HIG=1 HIW=0 QED=0 QCD=0, (h > Zp Zp, Zp > l+ l-)"""
proc_name = evgenConfig.process #just used for the mg directory name
#modifications to the param_card.dat (generated from the proc_card i.e. the specific model)
#if you want to see the resulting param_card, run Generate_tf with this jobo, and look at the param_card.dat in the cwd
#If you want to see the auto-calculated values of the decay widths, look at the one in <proc_name>/Cards/param_card.dat (again, after running a Generate_tf)
param_card_extras = { "HIDDEN": { 'epsilon': '1e-4', #kinetic mixing parameter
'kap': '1e-4', #higgs mixing parameter
'mzdinput': '125', #Zd mass
'mhsinput':'1000.0' }, #dark higgs mass
"HIGGS": { 'mhinput':'450'}, #higgs mass
"DECAY": { 'wzp':'Auto', 'wh':'46.82', 'wt':'Auto' } #auto-calculate decay widths and BR of Zp, H, t
}
run_card_extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'ptj':'0',
'ptb':'0',
'pta':'0',
'ptl':'0',
'etaj':'-1',
'etab':'-1',
'etaa':'-1',
'etal':'-1',
'drjj':'0',
'drbb':'0',
'drll':'0',
'draa':'0',
'drbj':'0',
'draj':'0',
'drjl':'0',
'drab':'0',
'drbl':'0',
'dral':'0' }
def post_lhe_hook():
#You can put code here that will get called after the 'generate' method of MadGraphControl is called
#What follows is an example for modifying the lifetime of the zd:
#The 'add_lifetimes' method is defined in the file included below. It modifies the lhe file to include particle lifetimes on the particle with given pdgId
#add_lifetimes(pdgId=32, avgtau=5)
pass
include("MC15JobOptions/MadGraphControl_Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
#Author: diallo boye <diallo.boye@cern.ch>
#This is a template joboption for generating events with the HAHM madgraph model
#To use it, copy the template and rename it by replacing the XXXXXX and YYYYYYYYY in the filename with:
# make XXXXXX match the runNumber you give to Generate_tf