temporarily work around mpi problems and run sherpa n times ^^

prog/runcards/pp/Sherpa.yaml

@@ -24,6 +24,10 @@ ME_GENERATORS:
 
 SCALES: 'VAR{1/2*Abs2(p[0]+p[1])}'
 
+PARTICLE_DATA:
+  94:
+    Massive: false
+
 # the main scattering process
 PROCESSES:
 - 94 -94 -> 22 22:

prog/runcards/pp_phaeno/Sherpa_Base.yaml

@@ -9,7 +9,7 @@ PDF_SET_VERSIONS: [0, 0]
 ME_GENERATORS:
 - Comix
 
-SCALES: 'VAR{1/2*Abs2(p[0]+p[1])}'
+#SCALES: 'VAR{1/2*Abs2(p[0]+p[1])}'
 
 # cuts
 SELECTORS:

prog/runcards/pp_phaeno/makefile

@@ -18,7 +18,7 @@ OUT_DIR = out
 EVENT_COUNT = 100000
 
 # List of Runcards, see RUNCARD_FOLDER
-RUNCARDS = basic with_pT
+RUNCARDS = basic with_pT with_pT_and_fragmentation
 
 # Base Runcard
 BASE_RUNCARD = Sherpa_Base.yaml
@@ -29,9 +29,11 @@ SHERPA_EXT_RUNCARD = Sherpa_Ext.yaml
 #                                    Macros                                   #
 ###############################################################################
 
+NUM_CPUS = $(shell nproc)
+EVENT_COUNT_CPU = $(shell echo "$(EVENT_COUNT) / $(NUM_CPUS)" | bc)
 THIS_DIR := $(dir $(abspath $(firstword $(MAKEFILE_LIST))))
 SET_RIVET := RIVET_ANALYSIS_PATH="$(THIS_DIR)$(ANALYSIS_DIR)"
-SHERPA_CMD := $(SHERPA) -e $(EVENT_COUNT) "RESULT_DIRECTORY: Results" "ANALYSIS_OUTPUT: analysis"
+SHERPA_ARGS := -e $(EVENT_COUNT_CPU) "RESULT_DIRECTORY: Results" "ANALYSIS_OUTPUT: analysis"
 OUT_YODA_FILES := $(foreach card,$(RUNCARDS),$(OUT_DIR)/$(card).yoda)
 
 # check if outdir is list
@@ -71,10 +73,9 @@ $(OUT_YODA_FILES): $(OUT_DIR)/%.yoda: $(RUNCARD_DIR)/%/$(SHERPA_EXT_RUNCARD) $(A
 	cat $(BASE_RUNCARD) $(R_DIR)$(SHERPA_EXT_RUNCARD) > $(R_DIR)Sherpa.yaml
 
 	@printf "\033[1;32mRunning Sherpa\033[0m\n"
-	cd $(R_DIR) && $(SHERPA_CMD) -e 0
+	cd $(R_DIR) && $(SHERPA) $(SHERPA_ARGS) -e 0
 	cd $(R_DIR) && \
-		$(SET_RIVET) mpirun --use-hwthread-cpus \
-			--use-hwthread-cpus $(SHERPA_CMD) "EVT_OUTPUT: 0"
+		$(SET_RIVET) $(THIS_DIR)/run_n_sherpas.fish $(NUM_CPUS) $(SHERPA_ARGS) "EVT_OUTPUT: 0"
 	cd $(R_DIR) && yodamerge *.yoda -o $(YODA_FILE) \
 		&& mv $(YODA_FILE) $(THIS_DIR)$(OUT_DIR)/
 

prog/runcards/pp_phaeno/run_n_sherpas.fish

@@ -0,0 +1,6 @@
+#!/bin/fish
+
+for i in (seq $argv[1])
+    Sherpa -R $i $argv[2..-1] "ANALYSIS_OUTPUT: analysis_$i" &
+end
+wait

prog/runcards/pp_phaeno/runcards/basic/Sherpa_Ext.yaml

@@ -14,3 +14,5 @@ MI_HANDLER: None
 
 # no transverse impulses
 BEAM_REMNANTS: false
+
+SCALES: 'VAR{1/2*Abs2(p[0]+p[1])}'

prog/runcards/pp_phaeno/runcards/with_pT/Sherpa_Ext.yaml

@@ -15,3 +15,5 @@ MI_HANDLER: None
 # showers are enabled as well
 INTRINSIC_KPERP:
   ENABLED: true
+
+SHOWER_GENERATOR: CSS

prog/runcards/pp_phaeno/runcards/with_pT_and_fragmentation/.gitignore

@@ -0,0 +1,3 @@
+*
+!Sherpa_Ext.yaml
+!.gitignore

prog/runcards/pp_phaeno/runcards/with_pT_and_fragmentation/Sherpa_Ext.yaml

@@ -0,0 +1,20 @@
+# the main scattering process
+PROCESSES:
+- 94 -94 -> 22 22:
+    Order: {QCD: 0, EW: 2}
+    Integration_Error: 0.001
+
+# use pure QED
+EW_SCHEME: 0
+
+# Enable fragmentation
+FRAGMENTATION: Ahadic
+
+MI_HANDLER: None
+
+# we explicitly enable transverse momenta
+# showers are enabled as well
+INTRINSIC_KPERP:
+  ENABLED: true
+
+SHOWER_GENERATOR: CSS