Extended statistics for single-top + photon NLO samples

Make sure that all the items in the checklists are duly followed:
follow the instructions of each checklist item and report potential issues to the repository mainteners.

JIRA ticket

ATLMCPROD-10765

Checklist for approval

This checklist should be normally filled by MC contacts (but it can also be done by the person opening the merge request).
If something is unclear get in touch with the repository maintainers.

• I have followed the guidelines in McSampleRequestProcedure
• I have provided the ATLMCPROD JIRA ticket where the request has been discussed
• The request has been approved by the PA/CP group
• The pipeline has run and the status is green (still check the CI job output - if you see any errors report to the repository maintainers)
• If the pipeline is orange: check the CI job output to sort out the issues. If the issue is unclear contact the repository maintainers
• CI jobs have not been skipped (if jobs have been skipped, the requesters should provide a reason for skipping and confirmation from the PMG conveners is necessary before merging)
• Check that no log.generate.short files are included in the commit to be merged to master. If such files are present, it indicates that something went wrong in the pipeline. Check the job logs and contact the package maintainers if needed.
• Check that the title of the MR is descriptive enough (if not ask the requester to modify it)
• Check that the DSID range assigned is correct (500000-999999) and corresponds to right generator block

523xxx/523588/mc.aMCPy8_tqgammaSM_tchan_4fl_NLO.py

+import MadGraphControl.MadGraph_NNPDF30NLOnf4_Base_Fragment
+from MadGraphControl.MadGraphUtils import *
+import os,subprocess,fileinput
+
+# General settings 
+
+evgenConfig.nEventsPerJob = 500
+nevents = runArgs.maxEvents*1.1 if runArgs.maxEvents>0 else evgenConfig.nEventsPerJob*1.1
+mode = 0
+#gridpack_mode=True
+#gridpack_dir='madevent/'
+
+name = 'tqgammaSM_tchan_4fl_NLO'
+runName='mc.MadGraph_'+str(name)
+
+process = """
+import model loop_sm
+define p = g u c d s u~ c~ d~ s~
+define j = g u c d s u~ c~ d~ s~ 
+define l+ = e+ mu+ ta+ 
+define l- = e- mu- ta-
+define vl = ve vm vt 
+define vl~ = ve~ vm~ vt~
+generate p p > t b~ j a $$ w+ w- [QCD]
+add process p p > t~ b j a $$ w+ w- [QCD] 
+output -f"""
+
+beamEnergy=-999
+if hasattr(runArgs,'ecmEnergy'):
+    beamEnergy = runArgs.ecmEnergy / 2.
+else:
+    raise RuntimeError("No center of mass energy found.")
+
+
+extras = { 'lhe_version'  : '3.0',
+           'parton_shower' :'PYTHIA8',
+           'ptgmin'        :10,
+           'r0gamma'       :0.2, 
+           'maxjetflavor'  :4,
+           'dynamical_scale_choice': '3', #sum of the transverse mass divided by 2
+           'ptl'           :0.0,
+           'ptj'           :0.0,
+           'etal'          :5.0,
+           'etagamma'      :5.0,
+           'etaj'          :-1,
+           'drll'          :0.0,
+           'nevents'       : nevents,
+           }
+
+
+process_dir = new_process(process)
+
+modify_run_card(process_dir=process_dir,runArgs=runArgs,settings=extras)      
+               
+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                                                                                                                          
+decay t > w+ b, w+ > l+ vl                                                                                                                                               
+decay t~ > w- b~, w- > l- vl~                                                                                                                                                                                                                                                                                                      
+# running the actual code                                                                                                                                                  
+launch"""%runArgs.randomSeed)                                                                                                                                              
+mscard.close()
+print_cards()
+
+
+generate(runArgs=runArgs,process_dir=process_dir)
+outputDS = arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=True)
+
+
+#### Shower
+evgenConfig.description = 'aMC@NLO_'+str(name)
+evgenConfig.generators += ["aMcAtNlo", "Pythia8"]
+evgenConfig.keywords+= ['SM', 'top',  'photon','singleTop','lepton'] 
+evgenConfig.contact = ['nils.julius.abicht@cern.ch']
+runArgs.inputGeneratorFile=outputDS
+
+
+include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
+
+
+
+if 'ATHENA_PROC_NUMBER' in os.environ:
+    evgenLog.info('Noticed that you have run with an athena MP-like whole-node setup.  Will re-configure now to make sure that the remainder of the job runs serially.')
+    njobs = os.environ.pop('ATHENA_PROC_NUMBER')
+    # Try to modify the opts underfoot
+    if not hasattr(opts,'nprocs'): 
+        mglog.warning('Did not see option!')
+    else: opts.nprocs = 0
+    print(opts)
+
+include("Pythia8_i/Pythia8_aMcAtNlo.py")
+#include("Pythia8_i/Pythia8_ShowerWeights.py")