1 # Using the magic encoding
   2 # -*- coding: utf-8 -*-
   3 # FLAG plot for mc_pdg
   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	= "$\mathrm{\overline{m}_c(\overline{m}_c)}$"	        		# plot title
  17 plotnamestring	= "mc"			# filename for plot
  18 #plotxticks	= ([[0.95, 1.00, 1.05, 1.10]])# locations for x-ticks
  19 plotxticks	= ([[1.25, 1.30, 1.35, 1.40]])# locations for x-ticks
  20 yaxisstrings	= [				# x,y-location for y-strings
  21 		   [1.07,13.,"$\\rm N_f=2+1+1$"],
  22 		   [1.07,5.0,"$\\rm N_f=2+1$"],
  23                    #[1.07,2.0,"$\\rm N_f=2$"]
  24 		   ]
  25 		   
  26 LABEL=1
  27 xaxisstrings	= [[1.45,-2.420,"GeV"]]		# x,y-location for x-strings
  28 
  29 		   
  30 xlimits		= [1.20,1.55]			# plot's xlimits
  31 logo		= 'upper left'			# either of 'upper right'
  32 					  	#           'upper left'
  33 					  	#           'lower left'
  34 					  	#           'lower right'
  35 tpos = 1.42					# 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 # 1		position on the y-axis
  45 # 2		neg. error
  46 # 3		pos. error
  47 # 4		Collaboration string
  48 # 5		this column contains layout parameters for the plot-symbol and
  49 #  		and the collaboration text:
  50 #		column	item
  51 #		0 	marker style (see bottom of this file for a full list)
  52 #		1	marker face color r=red, b=blue, k=black, w=white
  53 #					  g=green
  54 #		2	marker color (errorbar and frame), color coding as 
  55 #			for face color
  56 #		3	color intensity 0=full, 1=soso, 2=bleak
  57 dat2p1p1=[
  58 [1.348, 0.046, 0.046, 0.0, 0.0, "ETM 14" 	,['s','g','g',0,tpos]],
  59 [1.3478,0.0027, 0.0027, 0.0195, 0.0195, "ETM 14A" 	,['s','g','g',0,tpos]],
  60 [1.2715, 0.0095, 0.0095, 0.0, 0.0, "HPQCD 14A " 	,['s','g','g',0,tpos]],
  61 [1.273, 0.011, 0.011,0,0, "FNAL/MILC/TUMQCD 18" 	,['s','g','g',0,tpos]],
  62 [1.2757, 0.0084, 0.0084,0,0, "HPQCD 18"         ,['s','g','g',0,tpos]],
  63 ] 
  64 dat2p1=[
  65 [1.268, 0.009, 0.009,  0.0,     0.0,     "HPQCD 08B" ,['s','l','g',0,tpos]],
  66 [1.273, 0.006,0.006,  0.0,     0.0,     "HPQCD 10"  ,['s','g','g',0,tpos]],
  67 [1.304, 0.005, 0.005,  0.02,   0.02,   "$\chi$QCD 14"   ,['s','g','g',0,tpos]],
  68 [1.2871, 0.0123, 0.0123, 0.00, 0.00, "JLQCD 16"  ,['s','g','g',0,tpos]],
  69 [1.267, 0.012, 0.012, 0.00, 0.00, "Maezawa 16"  ,['s','w','r',0,tpos]],
  70 ]
  71 datapheno=[
  72 [1.275, 0.035, 0.025, 0, 0, "PDG",['^','k','k',0,tpos]],
  73 ]
  74 # The color coding for the FLAG-average (below) is as for the data itself; 
  75 # the additional list at the end is the RGB-code for the errorband.
  76 # Please note that the order inside the list FLAGaverage (if there is more 
  77 # than one) # is significant, the last item will be plotted last and will 
  78 # therefore be plotted on top of all others.
  79 
  80 
  81 FLAGaverage=[ # there should be as many entries for average as there are data sets
  82 	      # in the case where no average is provided for a given data set
  83 	      # just replace the central value by "NaN"
  84 	      # 
  85 	      # with respect to the above there is another column at the end specifying
  86 	      # the type of FLAG-band
  87 	      # 0 -> dashed lines marking the edge of the interval
  88 	      # 1 -> grey area bounded by solid black lines
  89 [np.nan, 0, 0, 0, 0, "",['^','k','k',0,tpos],1],
  90 #[np.nan, 0, 0, 0, 0, "",['^','k','k',0,tpos],1],
  91 #[np.nan, 2.5, 2.5, 0, 0, "",['^','k','k',0,tpos],1]]
  92 #[np.nan, 0, 0, 0, 0, "",['^','k','k',0,tpos],1],
  93 #[1.03, 0.04, 0.04, 0.0, 0.0, "FLAG average for $\\rm N_f=2$",['s','k','k',0,tpos],1],
  94 [1.275, 0.005, 0.005, 0.0, 0.0, "FLAG average for $\\rm N_f=2+1$",['s','k','k',0,tpos],1],
  95 [1.280, 0.013, 0.013, 0.0, 0.0, "FLAG average for $\\rm N_f=2+1+1$",['s','k','k',0,tpos],1]]
  96 # The follwing list should contain the list names of the previous DATA-blocks
  97 datasets=[datapheno,dat2p1,dat2p1p1]
  98 
  99 ################################################################################
 100 # DO NOT EDIT THE FOLLOWING ####################################################
 101 ################################################################################
 102 
 103 
 104 # Now that all the relevant has been specified, generate the plot.
 105 # Please make sure that you have a directory called "plots" in you current 
 106 # working directory. That's where all plots will be stored 
 107 #execfile('FLAGplot.py')
 108 #import FLAGplot.py
 109 #exec(open(FLAGplot.py).read())
 110 exec(compile(open('FLAGplot.py').read(), 'flagplot.py', 'exec'))
 111 
 112 # additional documentation:
 113 # '.' 	point marker
 114 # ',' 	pixel marker
 115 # 'o' 	circle marker
 116 # 'v' 	triangle_down marker
 117 # '^' 	triangle_up marker
 118 # '<' 	triangle_left marker
 119 # '>' 	triangle_right marker
 120 # '1' 	tri_down marker
 121 # '2' 	tri_up marker
 122 # '3' 	tri_left marker
 123 # '4' 	tri_right marker
 124 # 's' 	square marker
 125 # 'p' 	pentagon marker
 126 # '*' 	star marker
 127 # 'h' 	hexagon1 marker
 128 # 'H' 	hexagon2 marker
 129 # '+' 	plus marker
 130 # 'x' 	x marker
 131 # 'D' 	diamond marker
 132 # 'd' 	thin_diamond marker
 133 # '|' 	vline marker
 134 # '_' 	hline marker