Commit c3107690 authored by Joany Manjarres's avatar Joany Manjarres
Browse files

See ChangeLog

parent a7dfe649
2019-10-16 Joany Manjarres <joany@cern.ch>a
* tag MC15JobOptions-01-01-70
* add param_card.SM.N2N3.ZZ.dat param_card.SM.N2N3.Zh.dat param_card.SM.N2N3.hh.dat MadGraphControl_SimplifiedModel_N2N3.py
* update cvs 346707-714
* add JO 346788-89
2019-10-15 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-01-69
* add JOs 312580-598 395372-365 397207-216 398551-564
......
#--------------------------------------------------------------
# Standard pre-include
#
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
#--------------------------------------------------------------
# Some options for local testing. Shouldn't hurt anything in production.
#
# Nominal configuration for production: MadSpin decays, with 0,1,2+ partom emissions in the matrix element
#
#--------------------------------------------------------------
# Interpret the name of the JO file to figure out the mass spectrum
#
def MassToFloat(s):
if "p" in s:
return float(s.replace("p", "."))
return float(s)
#--------------------------------------------------------------
# split up the JO file input name to interpret it
# e.g. jobConfig: MadGraphControl_SimplifiedModel_N2N3_ZZ_800_200_MadSpin.py
splitConfig = runArgs.jobConfig[0].rstrip('.py').split('_')
#C1/N2/N3 degenerate
masses['1000025'] = MassToFloat(splitConfig[4])
masses['1000023'] = MassToFloat(splitConfig[4])
masses['1000024'] = MassToFloat(splitConfig[4])
#N1
masses['1000022'] = MassToFloat(splitConfig[5])
if masses['1000022']<0.5: masses['1000022']=0.5
# interpret the generation type, so we know which processes to run.
gentype = splitConfig[2] #will be N2N3
# decaytype is fixed. This really only tells MGC which param card to use.
# In principle this does nothing since it gets overwritten, but we leave
# it in for clarity! The JO needs to have "_HinoZh50_" in the proper position.
decaytype = splitConfig[3] #e.g. ZZ,Zh,hh,HinoZh50
madspindecays=False
if (decaytype == 'ZZ' or decaytype == 'Zh' or decaytype == 'Zh50') and ('MadSpin' in runArgs.jobConfig[0]) :
madspindecays = True;
print "gentype", gentype
print "decaytype", decaytype
print "madspindecays", madspindecays
# max number of jets will be two, unless otherwise specified.
njets = 2
#--------------------------------------------------------------
# MadGraph options
#
process = '''
define w = w+ w-
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define lv = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
define f = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~ u u~ d d~ c c~ s s~ b b~
generate p p > n2 n3 $ susystrong @1
add process p p > n2 n3 j $ susystrong @2
add process p p > n2 n3 j j $ susystrong @3
'''
mergeproc="{n2,1000023}{n3,1000025}"
msdecaystring=""
if madspindecays == True:
if decaytype == 'ZZ':
msdecaystring="decay n2 > z n1, z > f f \ndecay n3 > z n1, z > f f\n"
elif decaytype == 'Zh':
msdecaystring="decay n2 > z n1, z > f f \ndecay n3 > h01 n1\n"
elif decaytype == 'Zh50':
msdecaystring="decay n2 > z n1, z > f f \ndecay n2 > h01 n1\n decay n3 > z n1, z > f f \ndecay n3 > h01 n1\n"
# print the process, just to confirm we got everything right
print "Final process card:"
print process
#--------------------------------------------------------------
# Madspin configuration
#
if madspindecays==True:
if msdecaystring=="":
raise RuntimeError("Asking for MadSpin decays, but no decay string provided!")
madspin_card='madspin_card_N2N3_'+decaytype+'.dat'
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 Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
set BW_cut 15 # default for onshell
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
set spinmode none
# specify the decay for the final state particles
%s
# running the actual code
launch"""%(runArgs.randomSeed,msdecaystring))
mscard.close()
mergeproc+="LEPTONS,NEUTRINOS"
#--------------------------------------------------------------
# Pythia options
#
pythia = genSeq.Pythia8
pythia.Commands += ["23:mMin = 0.2"]
pythia.Commands += ["24:mMin = 0.2"]
# information about this generation
evgenLog.info('Registered generation of ~chi20 ~chi30 production, decaying into ~chi10 via ZZ/Zh/hh; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000025']) + ' ' + str(masses['1000022']))
evgenConfig.contact = [ "shion.chen@cern.ch" ]
evgenConfig.keywords += ['gaugino', 'chargino', 'neutralino']
evgenConfig.description = '~chi20 ~chi30 production, decaying into ~chi10 via ZZ/Zh/hh in simplified model. Br(N2->ZN1)=Br(N3->ZN1)=Br(N2->hN1)=Br(N3->hN1)=0.5, m_N2N3 = %s GeV, m_N1 = %s GeV'%(masses['1000025'],masses['1000022'])
#--------------------------------------------------------------
# No filter at the moment
evt_multiplier=2
evgenLog.info('inclusive processes will be generated')
#--------------------------------------------------------------
# Standard post-include
#
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
#--------------------------------------------------------------
# Merging options
#
if njets>0:
genSeq.Pythia8.Commands += [ "Merging:Process = pp>{n2,1000023}{n3,1000025}",
"1000025:spinType = 1",
"1000023:spinType = 1" ]
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
# JO for Pythia 8 jet jet JZ2 slice with fat jet filter of 160 GeV
evgenConfig.description = "Dijet truth jet slice JZ2 with fat jet filter of 160 GeV, with the A14 NNPDF23 LO tune"
evgenConfig.keywords = ["QCD", "jets", "SM"]
evgenConfig.contact = ["amoroso@cern.ch","christopher.young@cern.ch"]
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_ShowerWeights.py")
genSeq.Pythia8.Commands += ["HardQCD:all = on",
"PhaseSpace:Bias2Selection=on",
"PhaseSpace:pTHatMin = 15."]
include("MC15JobOptions/JetFilter_JZ2.py")
include("MC15JobOptions/AntiKt10TruthJets.py")
from GeneratorFilters.GeneratorFiltersConf import QCDTruthJetFilter
if "QCDTruthJetFilterLargeR" not in filtSeq:
filtSeq += QCDTruthJetFilter("QCDTruthJetFilterLargeR")
filtSeq.QCDTruthJetFilterLargeR.TruthJetContainer = "AntiKt10TruthJets"
filtSeq.QCDTruthJetFilterLargeR.MinPt = 160*GeV
evgenConfig.minevents = 1000
# JO for Pythia 8 jet jet JZ3 slice with fat jet filter of 400 GeV
evgenConfig.description = "Dijet truth jet slice JZ3 with fat jet filter of 400 GeV, with the A14 NNPDF23 LO tune"
evgenConfig.keywords = ["QCD", "jets", "SM"]
evgenConfig.contact = ["amoroso@cern.ch","christopher.young@cern.ch"]
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_ShowerWeights.py")
genSeq.Pythia8.Commands += ["HardQCD:all = on",
"PhaseSpace:Bias2Selection=on",
"PhaseSpace:pTHatMin = 50."]
include("MC15JobOptions/JetFilter_JZ3.py")
include("MC15JobOptions/AntiKt10TruthJets.py")
from GeneratorFilters.GeneratorFiltersConf import QCDTruthJetFilter
if "QCDTruthJetFilterLargeR" not in filtSeq:
filtSeq += QCDTruthJetFilter("QCDTruthJetFilterLargeR")
filtSeq.QCDTruthJetFilterLargeR.TruthJetContainer = "AntiKt10TruthJets"
filtSeq.QCDTruthJetFilterLargeR.MinPt = 400*GeV
evgenConfig.minevents = 1000
......@@ -2979,14 +2979,14 @@ Dsid,energy,inputeventfile,(inputconffile)
346704, 13000, mc15_13TeV.345038.PowhegPythia8EvtGen_NNPDF30_AZNLO_ZH125J_Zincl_MINLO.evgen.TXT.e5590
346705, 13000, mc15_13TeV.345039.PowhegPythia8EvtGen_NNPDF30_AZNLO_WpH125J_Wincl_MINLO.evgen.TXT.e5590
346706, 13000, mc15_13TeV.345040.PowhegPythia8EvtGen_NNPDF30_AZNLO_WmH125J_Wincl_MINLO.evgen.TXT.e5590
346707, 13000, ,group.phys-gener.MG5_aMCatNLO260.346707.tHjb110_4fl.TXT.mc15_v1
346708, 13000, ,group.phys-gener.MG5_aMCatNLO260.346708.tHjb122_4fl.TXT.mc15_v1
346709, 13000, ,group.phys-gener.MG5_aMCatNLO260.346709.tHjb123_4fl.TXT.mc15_v1
346710, 13000, ,group.phys-gener.MG5_aMCatNLO260.346710.tHjb124_4fl.TXT.mc15_v1
346711, 13000, ,group.phys-gener.MG5_aMCatNLO260.346711.tHjb126_4fl.TXT.mc15_v1
346712, 13000, ,group.phys-gener.MG5_aMCatNLO260.346712.tHjb127_4fl.TXT.mc15_v1
346713, 13000, ,group.phys-gener.MG5_aMCatNLO260.346713.tHjb130_4fl.TXT.mc15_v1
346714, 13000, ,group.phys-gener.MG5_aMCatNLO260.346714.tHjb140_4fl.TXT.mc15_v1
346707, 13000, ,group.phys-gener.MG5_aMCatNLO260.346707.tHjb110_4fl.TXT.mc15_v2
346708, 13000, ,group.phys-gener.MG5_aMCatNLO260.346708.tHjb122_4fl.TXT.mc15_v2
346709, 13000, ,group.phys-gener.MG5_aMCatNLO260.346709.tHjb123_4fl.TXT.mc15_v2
346710, 13000, ,group.phys-gener.MG5_aMCatNLO260.346710.tHjb124_4fl.TXT.mc15_v2
346711, 13000, ,group.phys-gener.MG5_aMCatNLO260.346711.tHjb126_4fl.TXT.mc15_v2
346712, 13000, ,group.phys-gener.MG5_aMCatNLO260.346712.tHjb127_4fl.TXT.mc15_v2
346713, 13000, ,group.phys-gener.MG5_aMCatNLO260.346713.tHjb130_4fl.TXT.mc15_v2
346714, 13000, ,group.phys-gener.MG5_aMCatNLO260.346714.tHjb140_4fl.TXT.mc15_v2
346729, 13000, mc15_13TeV.345054.PowhegPythia8EvtGen_NNPDF3_AZNLO_WpH125J_MINLO_lvbb_VpT.evgen.TXT.e5706
346730, 13000, mc15_13TeV.345053.PowhegPythia8EvtGen_NNPDF3_AZNLO_WmH125J_MINLO_lvbb_VpT.evgen.TXT.e5706
346764, 13000,, group.phys-gener.MG5_aMCatNLO262.346764.aMcAtNloPythia8_NNPDF30ME_tHjb_4fl_yt_minus1_CPalpha_0.TXT.mc15_v1
......
MC15JobOptions-01-01-69/
\ No newline at end of file
MC15JobOptions-01-01-70/
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment