file /home/anarendran/Documents/temp/rivet/pyext/build/lib.linux-x86_64-3.10/rivet/plotinfo.py

/home/anarendran/Documents/temp/rivet/pyext/build/lib.linux-x86_64-3.10/rivet/plotinfo.py

Namespaces

Name
rivet
rivet::plotinfo

Classes

Name
classrivet::plotinfo::PlotParser

Source code

from __future__ import print_function
import os, re
from .util import texpand
import sys

class PlotParser(object):
    """
    Reads Rivet's .plot files and determines which attributes to apply to each histo path.
    """

    pat_begin_block = re.compile(r'^(#*\s*)?BEGIN (\w+) ?(\S+)?')
    pat_begin_name_block = re.compile(r'^(#*\s*)?BEGIN (\w+) ?(\S+)? ?(\w+)?')
    pat_end_block =   re.compile(r'^(#*\s*)?END (\w+)')
    pat_comment = re.compile(r'^\s*#|^\s*$')
    pat_property = re.compile(r'^(\w+?)\s*=\s*(.*)$')
    pat_property_opt = re.compile('^ReplaceOption\[(\w+=\w+)\]=(.*)$')
    pat_path_property  = re.compile(r'^(\S+?)::(\w+?)=(.*)$')
    pat_paths = {}

    def __init__(self, plotpaths=None, addfiles=[]):
        """
        Parameters
        ----------
        plotpaths : list of str, optional
            The directories to search for .plot files.
            The default is to call the rivet.getAnalysisPlotPaths() function to get
            the directory where the .plot files can be found. (Usually equivalent to calling :command:`rivet-config --datadir`)

        Raises
        ------
        ValueError
            If `plotpaths` is not specified and calling
            :command:`rivet-config` fails.
        """
        self.addfiles = addfiles

        self.plotpaths = plotpaths
        if not self.plotpaths:
            try:
                import rivet
                self.plotpaths = rivet.getAnalysisPlotPaths()
            except Exception as e:
                sys.stderr.write("Failed to load Rivet analysis plot paths: %s\n" % e)
                raise ValueError("No plot paths given and the rivet module could not be loaded!")


    def getSection(self, section, hpath):
        """Get a section for a histogram from a .plot file.

        Parameters
        ----------
        section : ('PLOT'|'SPECIAL'|'HISTOGRAM')
            The section that should be extracted.
        hpath : str
            The histogram path, i.e. /AnalysisID/HistogramID .

        TODO:
         * Caching! The result of the lookup is not cached so every call requires a file to be searched for and opened.
        """
        if section not in ['PLOT', 'SPECIAL', 'HISTOGRAM']:
            raise ValueError("Can't parse section \'%s\'" % section)

        
        from rivet.aopaths import AOPath
        try:
            aop = AOPath(hpath)
        except ValueError:
            print("Found analysis object with non-standard path structure:", hpath, "... skipping")
            return None

        
        plotfile = aop.basepathparts()[0] + ".plot"
        ret = {'PLOT': {}, 'SPECIAL': None, 'HISTOGRAM': {}}
        for pidir in self.plotpaths:
            plotpath = os.path.join(pidir, plotfile)
            self._readHeadersFromFile_readHeadersFromFile_readHeadersFromFile(plotpath, ret, section, aop.basepath())
            
            if ret[section]: #< neatly excludes both empty dicts and None, used as null defaults above
                break

        
        for extrafile in self.addfiles:
            self._readHeadersFromFile_readHeadersFromFile_readHeadersFromFile(extrafile, ret, section, hpath)
        return ret[section]


    def getSections(self, sections, hpath):
        """Get all sections for a histogram from a .plot file.

        Parameters
        ----------
        section : ('SPECIAL')
            The section that should be extracted. Only Specials allowed to occur multiple times.
        hpath : str
            The histogram path, i.e. /AnalysisID/HistogramID .

        TODO:
         * Caching! The result of the lookup is not cached so every call requires a file to be searched for and opened.
        """
        if sections not in ['SPECIAL']:
            raise ValueError("Can't parse section \'%s\'" % section)

        
        from rivet.aopaths import AOPath
        try:
            aop = AOPath(hpath)
        except ValueError:
            print("Found analysis object with non-standard path structure:", hpath, "... skipping")
            return None

        
        plotfile = aop.basepathparts()[0] + ".plot"
        ret = {'SPECIAL': {}}
        for pidir in self.plotpaths:
            plotpath = os.path.join(pidir, plotfile)
            self._readNamedHeadersFromFile_readNamedHeadersFromFile_readNamedHeadersFromFile(plotpath, ret, sections, aop.basepath())
            
            if ret[sections]: #< neatly excludes both empty dicts and None, used as null defaults above
                break

        
        for extrafile in self.addfiles:
            self._readNamedHeadersFromFile_readNamedHeadersFromFile_readNamedHeadersFromFile(extrafile, ret, sections, hpath)
        return ret[sections]


    def _readHeadersFromFile(self, plotfile, ret, section, hpath):
        """Get a section for a histogram from a .plot file."""
        if not os.access(plotfile, os.R_OK):
            return
        startreading = False
        f = open(plotfile)
        msec = None
        for line in f:
            m = self.pat_begin_block.match(line)
            if m:
                tag, pathpat = m.group(2,3)
                # pathpat could be a regex
                if pathpat not in self.pat_paths:
                    try:
                        self.pat_paths[pathpat] = re.compile(pathpat)
                    except TypeError:
                        print("Error reading plot file for {}. Skipping.".format(plotfile))
                        return
                if tag == section:
                    m2 = self.pat_paths[pathpat].match(hpath)
                    if m2:
                        msec = m2
                        startreading = True
                        if section in ['SPECIAL']:
                            ret[section] = ''
                        continue
            if not startreading:
                continue
            if self.isEndMarkerisEndMarkerisEndMarker(line, section):
                startreading = False
                continue
            elif self.isCommentisCommentisComment(line):
                continue
            if section in ['PLOT', 'HISTOGRAM']:
                vm = self.pat_property.match(line)
                if vm:
                    prop, value = vm.group(1,2)
                    if msec:
                        oldval = value
                        try:
                            
                            value = value.encode("string-escape")
                            #print(value)
                            value = re.sub("(\\\\)(\\d)", "\\2", value) #< r-strings actually made this harder, since the \) is still treated as an escape!
                            #print(value)
                            value = msec.expand(value)
                            #print(value)
                        except Exception as e:
                            #print(e)
                            value = oldval #< roll back escapes if it goes wrong
                    ret[section][prop] = texpand(value) #< expand TeX shorthands
                vm = self.pat_property_opt.match(line)
                if vm:
                    prop, value = vm.group(1,2)
                    ret[section]['ReplaceOption[' + prop + ']'] = texpand(value)
            elif section in ['SPECIAL']:
                ret[section] += line
        f.close()


    def _readNamedHeadersFromFile(self, plotfile, ret, section, hpath):
        """Get a section for a histogram from a .plot file."""
        if not os.access(plotfile, os.R_OK):
            return
        startreading = False
        name = ''
        f = open(plotfile)
        for line in f:
            m = self.pat_begin_name_block.match(line)
            if m:
                tag, pathpat, name = m.group(2,3,4)
                if name == None:
                    continue
                # pathpat could be a regex
                if pathpat not in self.pat_paths:
                    self.pat_paths[pathpat] = re.compile(pathpat)
                if tag == section:
                    if self.pat_paths[pathpat].match(hpath):
                        startreading = True
                        if section in ['SPECIAL']:
                            ret[section][name] = ''
                        continue
            if not startreading:
                continue
            if self.isEndMarkerisEndMarkerisEndMarker(line, section):
                startreading = False
                continue
            elif self.isCommentisCommentisComment(line):
                continue
            if section in ['SPECIAL']:
                ret[section][name] += line
        f.close()


    def getHeaders(self, hpath):
        """Get the plot headers for histogram hpath.

        This returns the PLOT section.

        Parameters
        ----------
        hpath : str
            The histogram path, i.e. /AnalysisID/HistogramID .

        Returns
        -------
        plot_section : dict
            The dictionary usually contains the 'Title', 'XLabel' and
            'YLabel' properties of the respective plot.

        See also
        --------
        :meth:`getSection`
        """
        return self.getSectiongetSectiongetSection('PLOT', hpath)
    
    getPlot = getHeaders


    def getSpecial(self, hpath):
        """Get a SPECIAL section for histogram hpath.

        The SPECIAL section is only available in a few analyses.

        Parameters
        ----------
        hpath : str
            Histogram path. Must have the form /AnalysisID/HistogramID .

        See also
        --------
        :meth:`getSection`
        """
        return self.getSectiongetSectiongetSection('SPECIAL', hpath)

    def getSpecials(self, hpath):
        """Get all SPECIAL sections for histogram hpath.

        The SPECIAL section is only available in a few analyses.

        Parameters
        ----------
        hpath : str
            Histogram path. Must have the form /AnalysisID/HistogramID .

        See also
        --------
        :meth:`getSections`
        """
        return self.getSectionsgetSectionsgetSections('SPECIAL', hpath)


    def getHistogramOptions(self, hpath):
        """Get a HISTOGRAM section for histogram hpath.

        The HISTOGRAM section is only available in a few analyses.

        Parameters
        ----------
        hpath : str
            Histogram path. Must have the form /AnalysisID/HistogramID .

        See also
        --------
        :meth:`getSection`
        """
        return self.getSectiongetSectiongetSection('HISTOGRAM', hpath)


    def isEndMarker(self, line, blockname):
        m = self.pat_end_block.match(line)
        return m and m.group(2) == blockname


    def isComment(self, line):
        return self.pat_comment.match(line) is not None


    def updateHistoHeaders(self, hist):
        headers = self.getHeadersgetHeadersgetHeaders(hist.histopath)
        if "Title" in headers:
            hist.title = headers["Title"]
        if "XLabel" in headers:
            hist.xlabel = headers["XLabel"]
        if "YLabel" in headers:
            hist.ylabel = headers["YLabel"]


def mkStdPlotParser(dirs=None, addfiles=[]):
    """
    Make a PlotParser with the standard Rivet .plot locations automatically added to
    the manually set plot info dirs and additional files.
    """
    if dirs is None:
        dirs = []
    from .core import getAnalysisPlotPaths
    dirs += getAnalysisPlotPaths()
    seen = set()
    dirs = [d for d in dirs if d not in seen and not seen.add(d)]
    return PlotParser(dirs, addfiles)

Updated on 2022-08-07 at 20:46:08 +0100