add 377954-63, 412116-7,21, 346589-94

ChangeLog

-2019-05-16 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
+2019-06-05 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
+	* tagging MC15JobOptions-01-00-38
+	* add 377954-63, 412116-7,21, 346589-94
+	* add MadGraphControl_SimplifiedModel_GG_direct_LongLived_RHadron.py
+
+2019-06-04 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
 	* tagging MC15JobOptions-01-00-37
 	* add 450694-6, 450697-8, 312243, 449616-772, 346587, update 346551
 

common/MadGraph/MadGraphControl_SimplifiedModel_GG_direct_LongLived_RHadron.py

+include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
+
+infoStrings = runArgs.jobConfig[0].split("/")[-1].split(".")[2].split("_")
+
+masses['1000021'] = float(infoStrings[5])
+masses['1000022'] = float(infoStrings[6].split('.py')[0])
+if "stab" in infoStrings[7] :
+  lifetime = -1
+else :
+  lifetime = float(infoStrings[7].replace("ns","").split('.py')[0].replace('p','.') ) #in ns
+spectrum = int(infoStrings[8].split('.py')[0].replace('sp',''))
+if '_gl' in runArgs.jobConfig[0]:
+    gluinoballProb = float(infoStrings[-1].split('.')[0].replace("gl",""))*0.01
+else:
+    gluinoballProb = -1
+
+
+if masses['1000022']<0.5: masses['1000022']=0.5
+
+gentype   = str(infoStrings[2])
+decaytype = str(infoStrings[3])
+process   = '''
+generate p p > go go $ susysq susysq~ @1
+add process p p > go go j $ susysq susysq~ @2
+add process p p > go go j j $ susysq susysq~ @3
+'''
+njets = 2
+evt_multiplier = 6
+evgenLog.info('Registered generation of gluino grid '+str(runArgs.runNumber))
+
+evgenConfig.contact  = [ "lawrence.lee.jr@cern.ch" ]
+evgenConfig.keywords += ['simplifiedModel','gluino','longLived']
+evgenConfig.description = 'Gluino-gluino production, glu->qq+LSP in simplified model, m_glu = %d GeV, m_N1 = %d GeV. Gluino lifetime of %f, including RHadron setup. gluino-ball probability: %f, '%(masses['1000021'],masses['1000022'],lifetime,gluinoballProb)
+
+include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
+
+if njets>0:
+    genSeq.Pythia8.Commands += ["Merging:Process = pp>{go,1000021}{go,1000021}"]
+
+genSeq.Pythia8.Commands += ["Init:showChangedSettings = on"]
+genSeq.Pythia8.Commands += ["Rhadrons:allow = on"]
+genSeq.Pythia8.Commands += ["RHadrons:allowDecay = off"]
+# Check if we are explicitly setting the gluinoball fraction; if not use the default value
+if gluinoballProb!=-1:
+    genSeq.Pythia8.Commands += ["RHadrons:probGluinoball = %f"%gluinoballProb]
+else:
+    genSeq.Pythia8.Commands += ["RHadrons:probGluinoball = 0.1"]
+genSeq.Pythia8.Commands += ["Next:showScaleAndVertex = on"]
+genSeq.Pythia8.Commands += ["Check:nErrList = 2"]
+
+# Make sure that the RHadrons are allowed by TestHepMC
+if 'testSeq' in dir():
+    extra_pdgids_f = open('extra_pdgids.txt','w')
+    from RHadrons.RHadronMasses import offset_options
+    # Just allow all of them to have an anti-particle.  It's more future-proof.
+    for a in offset_options:
+        extra_pdgids_f.write(str(a)+'\n')
+        extra_pdgids_f.write('-'+str(a)+'\n')
+    extra_pdgids_f.close()
+    testSeq.TestHepMC.G4ExtraWhiteFile='extra_pdgids.txt'
+
+evgenConfig.specialConfig = "LIFETIME={lifetime};preInclude=SimulationJobOptions/preInclude.RHadronsPythia8.py".format(lifetime=lifetime)
+
+# Set up R-hadron masses in Pythia8
+import os
+if os.access('param_card.dat',os.R_OK):
+    from RHadrons.RHadronMasses import get_Pythia8_commands
+    genSeq.Pythia8.Commands += get_Pythia8_commands('param_card.dat',spectrum)
+# Otherwise this is going to be done a different way by the simulation pre-include
+

common/MadGraph/MadGraphControl_tWH_CP_NLO.py

@@ -4,41 +4,28 @@ from MadGraphControl.MadGraphUtils import *
 ## number of events to generate + safety margin
 nevents=1.1*runArgs.maxEvents
 runName='run_01'
-###generate events
-#gridpack_dir=None
-#gridpack_mode=False
-#cluster_type='condor'
-#cluster_queue='tomorrow'
-#mode=0
-#njobs=1
-###gridpack production
 gridpack_dir='madevent/'
 gridpack_mode=False
 mode=2 #mode=0 for single-core run, mode=1 if access to acluster, mode=2 for multicore production
 njobs=132
 cluster_type='lsf'
 cluster_queue='8nh'
-#mode=1#NO
-#njobs=20#NO
-#cluster_type='lsf'#NO
-#cluster_queue='1nd'#NO
-
 
 ######
 
 ## map DSID to process settings
 # select diagram removal scheme: DR1 or DR2?
 tWH_DR1 = []
-tWH_DR2 = [ 346486, 346511, 346536, 346537, 346545, 346546, 346547, 346548, 346549, 346550]
+tWH_DR2 = [ 346486, 346511, 346545, 346546, 346547, 346548, 346549, 346550, 346589, 346590, 346591, 346592, 346593, 346594]
 
 # select any BSM top Yukawa couplings (default: SM):
-#tWH_CPalpha_0 = [ 346486, 346511 ]
-tWH_CPalpha_15 = [ 346545 ]
-tWH_CPalpha_30 = [ 346546 ]
-tWH_CPalpha_45 = [ 346547, 346536 ]
-tWH_CPalpha_60 = [ 346548 ]
-tWH_CPalpha_75 = [ 346549 ]
-tWH_CPalpha_90 = [ 346550, 346537 ]
+#tWH_CPalpha_0 = [ 346486]
+tWH_CPalpha_15 = [ 346545, 346589]
+tWH_CPalpha_30 = [ 346546, 346590]
+tWH_CPalpha_45 = [ 346547, 346591]
+tWH_CPalpha_60 = [ 346548, 346592]
+tWH_CPalpha_75 = [ 346549, 346593]
+tWH_CPalpha_90 = [ 346550, 346594]
 
 tWH_yt_minus1 = []
 tWH_yt_plus2 = []

common/MadGraph/MadGraphControl_ttbar.py

@@ -17,11 +17,13 @@ MET200_filter     = [407048]
 pT300_cut     = [410368]
 
 # LHE production
-tt_inclusive_LHE = [410440]
+tt_inclusive_LHE = [410440, 412121]
 
+# for matching with H7
+herwig_dsids   = [412121]
 
-enhanced_bwcut = [410440]
-thisDSID     = runArgs.runNumber
+enhanced_bwcut = [410440, 412121]
+thisDSID       = runArgs.runNumber
 
 # --------------------------------------------------------------
 #  Some global production settings
@@ -46,6 +48,9 @@ dyn_scale = '10'    # user-defined scale -> Dominic's definition of mt+1/2*(pt^2
 lhe_version=3
 bwcut = 15
 
+if thisDSID in herwig_dsids:
+    parton_shower='HERWIGPP'
+
 # --------------------------------------------------------------
 #  Setting up the process
 # --------------------------------------------------------------
@@ -277,14 +282,32 @@ outputDS=arrange_output(run_name=name,
 # --------------------------------------------------------------
 # Run Pythia 8 Showering
 # --------------------------------------------------------------
-evgenConfig.generators    += ["aMcAtNlo", "Pythia8"]
-evgenConfig.description    = 'MG5_aMC@NLO+Pythia8+EvtGen '+name+' OTF, A14 NNPDF 2.3 LO, ME NNPDF 3.0 NLO, using scale sqrt(sum_i mT(i)**2/2)), for i = top quarks'
-evgenConfig.keywords      += [ 'SM', 'top']
-evgenConfig.contact        = [ 'ian.connelly@cern.ch', 'steffen.henkelmann@cern.ch' ]
-runArgs.inputGeneratorFile=outputDS
-
-include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
-include("MC15JobOptions/Pythia8_aMcAtNlo.py")
+
+if not runArgs.runNumber in herwig_dsids:
+    evgenConfig.generators    += ["aMcAtNlo", "Pythia8"]
+    evgenConfig.description    = 'MG5_aMC@NLO+Pythia8+EvtGen '+name+' OTF, A14 NNPDF 2.3 LO, ME NNPDF 3.0 NLO, using scale sqrt(sum_i mT(i)**2/2)), for i = top quarks'
+    evgenConfig.keywords      += [ 'SM', 'top']
+    evgenConfig.contact        = [ 'ian.connelly@cern.ch', 'steffen.henkelmann@cern.ch' ]
+    runArgs.inputGeneratorFile=outputDS    
+    include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+    include("MC15JobOptions/Pythia8_aMcAtNlo.py")
+elif runArgs.runNumber in herwig_dsids:
+    runArgs.inputGeneratorFile=outputDS
+    evgenConfig.generators += ["aMcAtNlo", "Herwig7"]
+    evgenConfig.description = 'MG5_aMC@NLO+Herwig7+EvtGen '+name+' OTF, H7p1 default tune, ME NNPDF 3.0 NLO, using scale sqrt(sum_i mT(i)**2/2)), for i = top quarks'
+    evgenConfig.tune = "H7.1-Default"
+    include("MC15JobOptions/Herwig7_LHEF.py")
+    Herwig7Config.me_pdf_commands(order="NLO", name="NNPDF30_nlo_as_0118")
+    Herwig7Config.tune_commands()
+    Herwig7Config.lhef_mg5amc_commands(lhe_filename=runArgs.inputGeneratorFile, me_pdf_order="NLO")
+    include("MC15JobOptions/Herwig71_EvtGen.py")
+    Herwig7Config.run()
+else:
+    theApp.finalize()
+    theApp.exit()
+
+
+
 
 # --------------------------------------------------------------
 # Apply TTbarWToLeptonFilter

share/DSID346xxx/MC15.346589.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_15.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID346xxx/MC15.346590.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_30.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID346xxx/MC15.346591.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_45.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID346xxx/MC15.346592.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_60.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID346xxx/MC15.346593.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_75.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID346xxx/MC15.346594.aMcAtNloPythia8EvtGen_tWH125_4fl_CPalpha_90.py

+include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')

share/DSID377xxx/MC15.377954.MGPy8EG_A14N23LO_GG_direct_1200_1100_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=1000

share/DSID377xxx/MC15.377955.MGPy8EG_A14N23LO_GG_direct_1200_1150_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=500

share/DSID377xxx/MC15.377956.MGPy8EG_A14N23LO_GG_direct_1300_1200_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=1000

share/DSID377xxx/MC15.377957.MGPy8EG_A14N23LO_GG_direct_1300_1250_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=500

share/DSID377xxx/MC15.377958.MGPy8EG_A14N23LO_GG_direct_1400_1300_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=1000

share/DSID377xxx/MC15.377959.MGPy8EG_A14N23LO_GG_direct_1400_1350_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=500

share/DSID377xxx/MC15.377960.MGPy8EG_A14N23LO_GG_direct_1500_1400_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=1000

share/DSID377xxx/MC15.377961.MGPy8EG_A14N23LO_GG_direct_1500_1450_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=500

share/DSID377xxx/MC15.377962.MGPy8EG_A14N23LO_GG_direct_1600_1500_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=1000

share/DSID377xxx/MC15.377963.MGPy8EG_A14N23LO_GG_direct_1600_1550_MET100.py

+include('MC15JobOptions/MadGraphControl_SimplifiedModel_GG_direct.py')
+
+evgenConfig.minevents=500