1 # Using the magic encoding
   2 # -*- coding: utf-8 -*-
   3 # FLAG plot for ms
   4 # check whether numpy library is present
   5 try: 
   6  import numpy
   7 except ImportError:
   8  print ("numpy library not found. If you want the FLAG-logo to be added to the")
   9  print ("plot, please install this library")
  10 import numpy as np
  11 ################################################################################
  12 # Please edit the following blocks #############################################
  13 ################################################################################
  14 
  15 # layout specifications
  16 titlestring	= "$\\rm m_s$"			# plot title
  17 plotnamestring	= "ms"			# filename for plot
  18 plotxticks	= ([[80., 90., 100.]])# locations for x-ticks
  19 yaxisstrings	= [				# x,y-location for y-strings
  20 		   [+70.,25.5,"$\\rm N_f=2+1+1$"],
  21 		   [+70.,13.5,"$\\rm N_f=2+1$"],
  22                    [+70.,3,"pheno."]
  23 		   ]
  24 		   
  25 LABEL=1
  26 xaxisstrings	= [				# x,y-location for x-strings
  27 		   [+121.,-3.3,"MeV"] 	 		    
  28                    ]		   
  29 		   
  30 xlimits		= [80.,140.]			# plot's xlimits
  31 logo		= 'upper left'			# either of 'upper right'
  32 					  	#           'upper left'
  33 					  	#           'lower left'
  34 					  	#           'lower right'
  35 tpos = 120.					# x-position for the data labels
  36 	
  37 
  38 # the following blocks contain the list of 
  39 # lattice results that will be plotted, one block for 2, 2+1 and 2+1+1, 
  40 # respectively
  41 # 
  42 # column	item
  43 # 0		central value
  44 # 2		neg. stat. error
  45 # 3		pos. stat. error
  46 # 4		neg. tot. error
  47 # 5		pos. tot. error
  48 # 6		Collaboration string
  49 # 7		this column contains layout parameters for the plot-symbol and
  50 #  		and the collaboration text:
  51 #		column	item
  52 #		0 	marker style (see bottom of this file for a full list)
  53 #		1	marker face color r=red, b=blue, k=black, w=white
  54 #					  g=green
  55 #		2	marker color (errorbar and frame), color coding as 
  56 #			for face color
  57 #		3	color intensity 0=full, 1=soso, 2=bleak
  58 #               4       position of label
  59 dat2p1p1=[
  60 [99.6, 4.3,4.3,3.6, 3.6, "ETM 14"		,['s','g','g',0,tpos]],
  61 [93.7, 0.8, 0.8, 0., 0., "HPQCD 14A "		,['s','g','g',0,tpos]],
  62 [92.52, 0.69,0.69,0.,0., "FNAL/MILC/TUMQCD 18"		,['s','g','g',0,tpos]],
  63 [94.49, 0.96,0.96,0.,0., "HPQCD 18  "		,['s','g','g',0,tpos]],
  64 ] 
  65 dat2p1=[
  66 # [76, 7.6, 7.6, 0.0, 0.0, "MILC 04, HPQCD/MILC/UKQCD 04" 	,['D','w','r',0,tpos]],
  67 # [87, 5.7, 5.7, 0.0, 0.0, "HPQCD 05"   	    	,['D','w','g',0,tpos]],
  68 # [90.1, 6.1, 17.2, 4.3, 4.3, "CP-PACS/JLQCD 07"  ,['D','w','r',0,tpos]],
  69 # [107.3, 11.7, 11.7, 4.4, 4.4, "RBC/UKQCD 08"	,['s','w','r',0,tpos]],
  70 # [72.72, 0.78, 0.78, 0.0, 0.0, "PACS-CS 08"	,['D','w','r',0,tpos]],
  71 # [88, 5.0, 5.0, 0.0, 0.0, "MILC 09"   	    	,['D','w','g',0,tpos]],
  72 [89, 4.78, 4.78, 0.2, 0.2, "MILC 09A"  	    	,['D','g','g',0,tpos]],
  73 [92.4, 1.5, 1.5, 0.0, 0.0, "HPQCD 09A" 		,['s','l','g',0,tpos]],
  74 [92.75, 1.11, 1.11, 0.58, 0.58, "PACS-CS 09"	,['s','w','r',0,tpos]],
  75 [97.6, 6.22, 6.22, 2.9, 2.9, "Blum 10"		,['s','w','r',0,tpos]],
  76 [96.2, 2.6, 2.6, 1.6, 1.6, "RBC/UKQCD 10A"	,['s','l','g',0,tpos]],
  77 [92.2, 1.3, 1.3, 0.0, 0.0, "HPQCD 10"		,['s','g','g',0,tpos]],
  78 [86.7, 2.3, 2.3, 0.0, 0.0, "PACS-CS 10"		,['s','w','r',0,tpos]],
  79 [95.5, 1.86, 1.86, 1.1, 1.1, "BMW 10A, 10B"	,['s','g','g',0,tpos]],
  80 #[94.2, 5.86, 5.86, 1.4, 1.4, "Laiho 11"		,['D','w','g',0,tpos]],
  81 [83.60, 2.30, 2.30, 0.58, 0.58, "PACS-CS 12"	,['s','w','r',0,tpos]],
  82 [92.31, 1.88, 1.88, 1.26, 1.26,  "RBC/UKQCD 12"	,['s','l','g',0,tpos]],
  83 [90.3, 1.35, 1.35, 0.9, 0.9,  "RBC/UKQCD 14B"	,['s','g','g',0,tpos]],
  84 [92.0, 1.7, 1.7, 0.0, 0.0,  "Maezawa 16"	,['s','g','g',0,tpos]],
  85 ]
  86 datapheno=[
  87 [102, 0,0,0,0, "Vainshtein 78",['o','b','b',0,tpos]],
  88 [104, 15,15,0,0, "Narison 06",['o','b','b',0,tpos]],
  89 [92, 9,9,0,0, "Jamin 06",['o','b','b',0,tpos]],
  90 [105, 9.2,9.2,0,0, "Chetyrkin 06",['o','b','b',0,tpos]],
  91 [103, 8,8,0,0, "Dominguez 09",['o','b','b',0,tpos]],
  92 #[93.5, 2.5, 2.5, 0, 0, "PDG",['^','k','k',0,tpos]],
  93 [95, 5, 5, 0, 0, "PDG",['^','k','k',0,tpos]],
  94 ]
  95 # datapheno=[
  96 # [95, 5, 5, 5, 5, "PDG",['s','w','k',0,tpos],1]
  97 # ]
  98 # The color coding for the FLAG-average (below) is as for the data itself; 
  99 # the additional list at the end is the RGB-code for the errorband.
 100 # Please note that the order inside the list FLAGaverage (if there is more 
 101 # than one) # is significant, the last item will be plotted last and will 
 102 # therefore be plotted on top of all others.
 103 
 104 FLAGaverage=[ # there should be as many entries for average as there are data sets
 105 	      # in the case where no average is provided for a given data set
 106 	      # just replace the central value by "NaN"
 107 	      # 
 108 	      # with respect to the above there is another column at the end specifying
 109 	      # the type of FLAG-band
 110 	      # 0 -> dashed lines marking the edge of the interval
 111 	      # 1 -> grey area bounded by solid black lines
 112 #[95, 6.32, 6.32, 2, 2, "our estimate for $N_f=2^x$",['s','w','k',0,tpos],1],
 113 #[95, 5, 5, 0, 0, "",['^','k','k',0,tpos],1],
 114 [np.nan, 2.5, 2.5, 0, 0, "",['^','k','k',0,tpos],1],
 115 #[0,0,0,0.,0., "dummy",['s','k','k',0,0],0],
 116  [92.03,0.88,0.88,0.,0., "FLAG average for $\\rm N_f=2+1$",['s','k','k',0,tpos],1],
 117 [93.12,0.69,0.69,0.,0., "FLAG average for $\\rm N_f=2+1+1$",['s','k','k',0,tpos],1]]
 118 
 119 # The follwing list should contain the list names of the previous DATA-blocks
 120 datasets=[datapheno,dat2p1,dat2p1p1]
 121 
 122 ################################################################################
 123 # DO NOT EDIT THE FOLLOWING ####################################################
 124 ################################################################################
 125 
 126 
 127 # Now that all the relevant has been specified, generate the plot.
 128 # Please make sure that you have a directory called "plots" in you current 
 129 # working directory. That's where all plots will be stored 
 130 #execfile('FLAGplot.py')
 131 exec(compile(open('FLAGplot.py').read(), 'FLAGplot.py', 'exec'))
 132 
 133 # additional documentation:
 134 # '.' 	point marker
 135 # ',' 	pixel marker
 136 # 'o' 	circle marker
 137 # 'v' 	triangle_down marker
 138 # '^' 	triangle_up marker
 139 # '<' 	triangle_left marker
 140 # '>' 	triangle_right marker
 141 # '1' 	tri_down marker
 142 # '2' 	tri_up marker
 143 # '3' 	tri_left marker
 144 # '4' 	tri_right marker
 145 # 's' 	square marker
 146 # 'p' 	pentagon marker
 147 # '*' 	star marker
 148 # 'h' 	hexagon1 marker
 149 # 'H' 	hexagon2 marker
 150 # '+' 	plus marker
 151 # 'x' 	x marker
 152 # 'D' 	diamond marker
 153 # 'd' 	thin_diamond marker
 154 # '|' 	vline marker
 155 # '_' 	hline marker