Add a 2-particle RG-L EB test (#575)

* add rg-l generator

* add rg-l test

* 12 degrees for inbending

* fix column types

* fix RNG seed

* fix seed

* fix lund "index"

Author: baltzell

validation/advanced-tests/src/eb/scripts/gemc.sh

@@ -3,9 +3,10 @@
 usage() { echo "Usage: $0 [-g GEMC] [-n NEV] [-p PARTS] [-c GCARD] [-m]" 1>&2; exit $1; }
 
 run=11
-gemc=5.10
+gemc=dev
 nevents=100
 particles=()
+seed=1234567
 
 while getopts "g:n:c:p:mdh" o
 do
@@ -53,7 +54,7 @@ function run_gemc () {
     [ -e "$_output" ] && echo Output file already exists:  $_output && exit 3
     [ -z ${dryrun+x} ] && set -x
     $dryrun gemc \
-        $_gcard -RUNNO=$_run -USE_GUI=0 -N=$_nevents \
+        $_gcard -RANDOM=$seed -RUNNO=$_run -USE_GUI=0 -N=$_nevents \
         -SAVE_ALL_MOTHERS=1 -SKIPREJECTEDHITS=1 -INTEGRATEDRAW="*" -NGENP=50 \
         -INPUT_GEN_FILE="LUND, $_input" \
         -OUTPUT="hipo, $_output" &

validation/advanced-tests/src/eb/scripts/gen.py

@@ -0,0 +1,26 @@
+#!/usr/bin/env python
+import random,math,sys
+
+mp = 0.93827
+me = 0.00051
+fmt = ' '.join(6*['%d']+8*['%.4f'])
+
+def cartesian(p,t,h,m):
+    x = p * math.cos(h) * math.sin(t)
+    y = p * math.sin(h) * math.sin(t)
+    z = p * math.cos(t)
+    e = math.sqrt(p*p + m*m)
+    return [x,y,z,e]
+
+for i in range(int(sys.argv[1])):
+    p = 0.1+0.1*random.random()
+    e = math.sqrt(p*p+mp*mp)
+    t = math.radians(70+40*random.random())
+    h = 2*math.pi*random.random()
+    z = -5+10*random.random()
+    header = [2,1,1,0,0,0,0,0,0,0]
+    electron = [1,0,1,11,0,0] + cartesian(9.0,math.radians(12),0,me) + [me,0,0,z]
+    hadron = [2,0,1,2212,0,0] + cartesian(p,t,h,mp) + [mp,0,0,z]
+    print(' '.join([str(x) for x in header]))
+    print(fmt%tuple(electron))
+    print(fmt%tuple(hadron))

validation/advanced-tests/src/eb/scripts/list.txt

@@ -12,3 +12,4 @@ electrondeuteronC -pid 45 -cd
 electronFTpion -pid -211 -ft
 electronFTproton -pid -2212 -ft
 electronFTkaon -pid -321 -ft
+electronprotonA -pid -2212 -al