gp_rdiff: fix plot for missing medium

tschaume · Apr 25, 2014 · c4b126e · c4b126e
1 parent da62909
commit c4b126e
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/pyana/examples/gp_rdiff.py b/pyana/examples/gp_rdiff.py
@@ -92,6 +92,9 @@ def gp_rdiff(version, nomed, noxerr):
         # build list of data points
         if key in dataOrdered: dataOrdered[key].append(dp)
         else: dataOrdered[key] = [ dp ]
+    if energy not in medium: # add pseudo-data for missing medium (27GeV)
+      key = ' '.join([energy, 'GeV (Med.)'])
+      dataOrdered[key] = [ [0.1,-1,0,0,0], [1,-1,0,0,0] ]
 
   # make plot
   nSets = len(dataOrdered)
@@ -170,4 +173,4 @@ def gp_rdiff(version, nomed, noxerr):
   logging.basicConfig(
     format='%(message)s', level=getattr(logging, loglevel)
   )
-  print gp_rdiff(version, nomed, noxerr)
+  print gp_rdiff(args.version, args.nomed, args.noxerr)