diff --git a/cscripts/Makefile.am b/cscripts/Makefile.am index 3f01900..3a82f60 100644 --- a/cscripts/Makefile.am +++ b/cscripts/Makefile.am @@ -32,6 +32,7 @@ pythonscripts = $(top_builddir)/cscripts/__init__.py \ $(srcdir)/csobsdef.py \ $(srcdir)/csobsinfo.py \ $(srcdir)/csobsselect.py \ + $(srcdir)/csphasecrv.py \ $(srcdir)/cspull.py \ $(srcdir)/csresmap.py \ $(srcdir)/csroot2caldb.py \ @@ -63,6 +64,7 @@ parfiles = $(srcdir)/cscaldb.par \ $(srcdir)/csobsdef.par \ $(srcdir)/csobsinfo.par \ $(srcdir)/csobsselect.par \ + $(srcdir)/csphasecrv.par \ $(srcdir)/cspull.par \ $(srcdir)/csresmap.par \ $(srcdir)/csroot2caldb.par \ @@ -91,6 +93,7 @@ linkscripts = cscaldb \ csobsdef \ csobsinfo \ csobsselect \ + csphasecrv \ cspull \ csresmap \ csroot2caldb \ @@ -183,7 +186,8 @@ testdata = $(top_srcdir)/test/data/crab.xml \ $(top_srcdir)/test/data/iact_runlist_1.dat \ $(top_srcdir)/test/data/iact_runlist_2.dat \ $(top_srcdir)/test/data/model_temporal_phasecurve.xml \ - $(top_srcdir)/test/data/model_temporal_phasecurve.fits + $(top_srcdir)/test/data/model_temporal_phasecurve.fits \ + $(top_srcdir)/test/data/phased_events.fits.gz # Test data to be distributed in $(pyexecdir)/cscripts/tests/data/iactdata testiactdata = $(top_srcdir)/test/data/iactdata/hdu-index.fits.gz \ diff --git a/cscripts/__init__.py.in b/cscripts/__init__.py.in index 21554fa..e079b96 100644 --- a/cscripts/__init__.py.in +++ b/cscripts/__init__.py.in @@ -35,6 +35,7 @@ __all__ = [ "csobsdef", "csobsinfo", "csobsselect", + "csphasecrv", "cspull", "csresmap", "csroot2caldb", @@ -69,6 +70,7 @@ from cscripts.csobs2caldb import csobs2caldb from cscripts.csobsdef import csobsdef from cscripts.csobsinfo import csobsinfo from cscripts.csobsselect import csobsselect +from cscripts.csphasecrv import csphasecrv from cscripts.cspull import cspull from cscripts.csresmap import csresmap from cscripts.csroot2caldb import csroot2caldb diff --git a/cscripts/cscaldb.py b/cscripts/cscaldb.py index c6b94ab..d43f342 100755 --- a/cscripts/cscaldb.py +++ b/cscripts/cscaldb.py @@ -2,7 +2,7 @@ # ========================================================================== # Show the content of the ctools calibration database # -# Copyright (C) 2014-2016 Juergen Knoedlseder +# Copyright (C) 2014-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -55,7 +55,7 @@ class cscaldb(ctools.cscript): """ # Set name self._name = 'cscaldb' - self._version = '1.1.0' + self._version = ctools.__version__ # Initialise application by calling the appropriate class # constructor. diff --git a/cscripts/csebins.py b/cscripts/csebins.py index 9375529..b9507ca 100755 --- a/cscripts/csebins.py +++ b/cscripts/csebins.py @@ -40,7 +40,7 @@ class csebins(ctools.cscript): """ # Set name and version self._name = 'csebins' - self._version = '1.0.0' + self._version = ctools.__version__ # Set members self._ebounds = gammalib.GEbounds() diff --git a/cscripts/csfindobs.py b/cscripts/csfindobs.py index 7b6b426..23014cb 100755 --- a/cscripts/csfindobs.py +++ b/cscripts/csfindobs.py @@ -2,7 +2,7 @@ # ========================================================================== # Find observations from an IACT data store # -# Copyright (C) 2016 Michael Mayer +# Copyright (C) 2016-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -40,7 +40,7 @@ class csfindobs(ctools.cscript): """ # Set name self._name = 'csfindobs' - self._version = '1.2.0' + self._version = ctools.__version__ self._datapath = os.getenv('VHEFITS','') self._prodname = '' self._select_radec = True diff --git a/cscripts/csiactcopy.py b/cscripts/csiactcopy.py index 4ffa5e6..fe11a3f 100755 --- a/cscripts/csiactcopy.py +++ b/cscripts/csiactcopy.py @@ -44,7 +44,7 @@ class csiactcopy(ctools.cscript): """ # Set name self._name = 'csiactcopy' - self._version = '1.2.0' + self._version = ctools.__version__ self._datapath = os.getenv('VHEFITS','') self._remote_master = gammalib.GFilename() self._remote_base = '' diff --git a/cscripts/csiactdata.py b/cscripts/csiactdata.py index e021255..934d760 100755 --- a/cscripts/csiactdata.py +++ b/cscripts/csiactdata.py @@ -2,7 +2,7 @@ # ========================================================================== # Inspect an IACT data storage # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -45,7 +45,7 @@ class csiactdata(ctools.cscript): """ # Set name self._name = 'csiactdata' - self._version = '1.1.0' + self._version = ctools.__version__ self._datapath = os.getenv('VHEFITS','') self._prodnames = [] self._master_indx = '' diff --git a/cscripts/csiactobs.py b/cscripts/csiactobs.py index 1bf704a..3642d80 100755 --- a/cscripts/csiactobs.py +++ b/cscripts/csiactobs.py @@ -2,7 +2,7 @@ # ========================================================================== # Generates an IACT observation definition XML file. # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -47,7 +47,7 @@ class csiactobs(ctools.cscript): """ # Set name and version self._name = 'csiactobs' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise some members self._obs = gammalib.GObservations() diff --git a/cscripts/cslightcrv.par b/cscripts/cslightcrv.par index 597d5f3..825fcb9 100644 --- a/cscripts/cslightcrv.par +++ b/cscripts/cslightcrv.par @@ -36,7 +36,7 @@ emin, r, a, 0.1,,, "Lower energy limit of events (TeV)" emax, r, a, 100.0,,, "Upper energy limit of events (TeV)" enumbins, i, a, 0,0,100, "Number of energy bins for binned (0=unbinned)" coordsys, s, a, CEL,CEL|GAL,, "Coordinate System" -proj, s, a, CAR,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" nxpix, i, a, 200,,, "Size of the X axis in pixels" diff --git a/cscripts/cslightcrv.py b/cscripts/cslightcrv.py index 2a7c27f..401dfe2 100755 --- a/cscripts/cslightcrv.py +++ b/cscripts/cslightcrv.py @@ -69,7 +69,7 @@ class cslightcrv(ctools.cscript): """ # Set name self._name = 'cslightcrv' - self._version = '1.3.0' + self._version = ctools.__version__ # Initialise some members self._srcname = '' @@ -509,11 +509,7 @@ class cslightcrv(ctools.cscript): self._get_parameters() # Write observation into logger - if self._logTerse(): - self._log('\n') - self._log.header1(gammalib.number('Input observation',len(self._obs))) - self._log(str(self._obs)) - self._log('\n') + self._log_observations(gammalib.NORMAL, self._obs, 'Observation') # Get time boundaries tmin = self._tbins.tstart(0) diff --git a/cscripts/csmodelinfo.py b/cscripts/csmodelinfo.py index 7641b21..9d45f1f 100755 --- a/cscripts/csmodelinfo.py +++ b/cscripts/csmodelinfo.py @@ -2,7 +2,7 @@ # ========================================================================== # Shows the content of a model container # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -44,7 +44,7 @@ class csmodelinfo(ctools.cscript): """ # Set name self._name = 'csmodelinfo' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise class members self._models = gammalib.GModels() diff --git a/cscripts/csmodelmerge.py b/cscripts/csmodelmerge.py index ad285d9..151e4e0 100755 --- a/cscripts/csmodelmerge.py +++ b/cscripts/csmodelmerge.py @@ -2,7 +2,7 @@ # ========================================================================== # Merge model definition XML files # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -43,7 +43,7 @@ class csmodelmerge(ctools.cscript): """ # Set name self._name = 'csmodelmerge' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise class members self._files = None diff --git a/cscripts/csmodelselect.py b/cscripts/csmodelselect.py index 2f39643..e38e465 100755 --- a/cscripts/csmodelselect.py +++ b/cscripts/csmodelselect.py @@ -2,7 +2,7 @@ # ========================================================================== # Selects model from a model definition XML file # -# Copyright (C) 2016 Juergen Knoedlseder +# Copyright (C) 2016-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -38,7 +38,7 @@ class csmodelselect(ctools.cscript): """ # Set name self._name = 'csmodelselect' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise observation container from constructor arguments self._obs, argv = self._set_input_obs(argv) diff --git a/cscripts/csobs2caldb.py b/cscripts/csobs2caldb.py index 9c9d55a..3cde977 100755 --- a/cscripts/csobs2caldb.py +++ b/cscripts/csobs2caldb.py @@ -2,7 +2,7 @@ # ========================================================================== # Creates a calibration database entry for an IACT observation. # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -48,7 +48,7 @@ class csobs2caldb(ctools.cscript): """ # Set name and version self._name = 'csobs2caldb' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise members self._observation = gammalib.GCTAObservation() diff --git a/cscripts/csobsdef.py b/cscripts/csobsdef.py index d7d960b..05b3cd4 100755 --- a/cscripts/csobsdef.py +++ b/cscripts/csobsdef.py @@ -86,7 +86,7 @@ class csobsdef(ctools.cscript): """ # Set name and version self._name = 'csobsdef' - self._version = '1.3.0' + self._version = ctools.__version__ # Initialise class members self._obs = gammalib.GObservations() diff --git a/cscripts/csobsinfo.py b/cscripts/csobsinfo.py index eae9633..7249b4c 100755 --- a/cscripts/csobsinfo.py +++ b/cscripts/csobsinfo.py @@ -2,7 +2,7 @@ # ========================================================================== # Dump information about observation into log file # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -38,7 +38,7 @@ class csobsinfo(ctools.cscript): """ # Set name self._name = 'csobsinfo' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise class members self._obj_dir = None diff --git a/cscripts/csobsselect.py b/cscripts/csobsselect.py index a444184..bcc9969 100755 --- a/cscripts/csobsselect.py +++ b/cscripts/csobsselect.py @@ -2,7 +2,7 @@ # ========================================================================== # Selects observations from an observation definition XML file # -# Copyright (C) 2016 Juergen Knoedlseder +# Copyright (C) 2016-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -38,7 +38,7 @@ class csobsselect(ctools.cscript): """ # Set name self._name = 'csobsselect' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise observation container from constructor arguments self._obs, argv = self._set_input_obs(argv) diff --git a/cscripts/csphasecrv.par b/cscripts/csphasecrv.par index 7e08f68..fea2837 100644 --- a/cscripts/csphasecrv.par +++ b/cscripts/csphasecrv.par @@ -27,14 +27,14 @@ outfile, f, a, phasecrv.fits,,, "Output fits file name." # # Script parameters #================== -tbinalg, s, a, LIN,LIN|FILE,, "Algorithm for defining time bins" +phbinalg, s, a, LIN,LIN|FILE,, "Algorithm for defining phase bins" phbins, i, a, 5,1,10000, "Number of phase bins" phbinfile, f, a, NONE,,, "File defining the phase binning" emin, r, a, 0.1,,, "Lower energy limit of events (TeV)" emax, r, a, 100.0,,, "Upper energy limit of events (TeV)" enumbins, i, a, 0,0,100, "Number of energy bins for binned (0=unbinned)" coordsys, s, a, CEL,CEL|GAL,, "Coordinate System" -proj, s, a, CAR,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" nxpix, i, a, 200,,, "Size of the X axis in pixels" @@ -44,6 +44,7 @@ binsz, r, a, 0.02,,, "Pixel size (deg/pixel)" # # Standard parameters #==================== +publish, b, h, no,,, "Publish phase curve on VO Hub?" chatter, i, h, 2,0,4, "Chattiness of output" clobber, b, h, yes,,, "Overwrite existing output files with new output files?" debug, b, h, no,,, "Debugging mode activated" diff --git a/cscripts/csphasecrv.py b/cscripts/csphasecrv.py index fa3d769..695d57c 100755 --- a/cscripts/csphasecrv.py +++ b/cscripts/csphasecrv.py @@ -1,8 +1,8 @@ #! /usr/bin/env python # ========================================================================== -# Generates pulsar spectra as a function of phase +# Generates spectra as a function of phase # -# Copyright (C) 2014-2017 Rolf Buehler +# Copyright (C) 2017 Rolf Buehler # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -29,22 +29,22 @@ from cscripts import obsutils # ================ # class csphasecrv(ctools.cscript): """ - Generates a spectra in phase bins + Generates spectra in phase bins This script computes spectra by performing a maximum likelihood fit - using :doc:`ctlike` in a series of phase bins for pulsars. - The phase bins can be either specified in an ASCII file, as an interval - divided into equally sized phase bins. + using :doc:`ctlike` in a series of phase bins for pulsars. The phase bins + can be either specified in an ASCII file or as an interval divided into + equally sized phase bins. Examples: - >>> phrv = csphasecrv() - >>> phrv.run() + >>> phcrv = csphasecrv() + >>> phcrv.run() >>> ... (querying for parameters) ... - >>> phrv = phrv.phasecrv() + >>> phcrv = phrv.phasecrv() Generates phase fits and retrieves dictionary with best fit models. - >>> lcrv = csphasecrv(obs) - >>> lcrv.execute() + >>> phcrv = csphasecrv(obs) + >>> phcrv.execute() >>> ... (querying for parameters) ... Generates phase fits from the observations and saves results in XML files. @@ -57,13 +57,14 @@ class csphasecrv(ctools.cscript): """ # Set name self._name = 'csphasecrv' - self._version = '1.0.0' + self._version = ctools.__version__ - # Initialise some members + # Initialise some members. Phases are stored in a nested list + # [[ph1min,ph1max], [ph2min,ph2max],..] self._srcname = '' - # Phases are stored in a nested list [[ph1min,ph1max], [ph2min,ph2max],..] - self._phbins = [[0.0,1.0]] + self._phbins = [[0.0,1.0]] self._stacked = False + self._fits = gammalib.GFits() # Initialise observation container from constructor arguments. self._obs, argv = self._set_input_obs(argv) @@ -144,7 +145,7 @@ class csphasecrv(ctools.cscript): # Get algorithm to use for defining the time intervals. Possible # values are "FILE" or "LIN". This is enforced at # parameter file level, hence no checking is needed. - algorithm = self['tbinalg'].string() + algorithm = self['phbinalg'].string() # If the algorithm is "FILE" then handle a FITS or an ASCII file for # the time bin definition @@ -173,8 +174,8 @@ class csphasecrv(ctools.cscript): # Compute time step and setup the GTIs ph_step = 1.0 / float(nbins) for i in range(nbins): - phmin = i *ph_step - phmax = (i+1)*ph_step + phmin = i * ph_step + phmax = (i+1) * ph_step phbins.append([phmin,phmax]) # Return Good Time Intervals @@ -252,6 +253,168 @@ class csphasecrv(ctools.cscript): # Return new oberservation container return new_obs + def _get_free_par_names(self): + """ + Return list of free parameter names + + Returns + ------- + names : list of str + List of free parameter names. + """ + # Initialise empty list of free parameter names + names = [] + + # Collect list of free parameter names + for par in self._obs.models()[self._srcname]: + if par.is_free(): + names.append(par.name()) + + # Return names + return names + + def _create_fits_table(self, results): + """ + Creates FITS binary table containing light curve results + + Parameters + ---------- + results : list of dict + List of result dictionaries + + Returns + ------- + table : `~gammalib.GFitsBinTable` + FITS binary table containing light curve + """ + # Determine number of rows in FITS table + nrows = len(results) + + # Create FITS Table with extension "PHASECURVE" + table = gammalib.GFitsBinTable(nrows) + table.extname('PHASECURVE') + + # Append phase columns "PHASE_MIN" and "PHASE_MAX" + phmin = gammalib.GFitsTableDoubleCol('PHASE_MIN', nrows) + phmax = gammalib.GFitsTableDoubleCol('PHASE_MAX', nrows) + for i, result in enumerate(results): + phmin[i] = result['phmin'] + phmax[i] = result['phmax'] + table.append(phmin) + table.append(phmax) + + # Create parameter columns + for par in self._obs.models()[self._srcname]: + if par.is_free(): + name = par.name() + e_name = 'e_'+par.name() + col = gammalib.GFitsTableDoubleCol(name, nrows) + e_col = gammalib.GFitsTableDoubleCol(e_name, nrows) + col.unit(par.unit()) + e_col.unit(par.unit()) + for i, result in enumerate(results): + col[i] = result['values'][name] + e_col[i] = result['values'][e_name] + table.append(col) + table.append(e_col) + + # Return table + return table + + def _create_fits(self): + """ + Create FITS file object from fit results + """ + # Initialise list of result dictionaries + results = [] + + # Get source parameters + pars = self._get_free_par_names() + + # Loop over time bins + for i in range(len(self._phbins)): + + # Get time boundaries + phmin = self._phbins[i][0] + phmax = self._phbins[i][1] + + # Initialise result dictionary + result = {'phmin': phmin, + 'phmax': phmax, + 'pars': pars, + 'values': {}} + + # Store fit results + phname = str(phmin)+'-'+str(phmax) + source = self._fitmodels[phname][self._srcname] + for par in pars: + result['values'][par] = source[par].value() + result['values']['e_'+par] = source[par].error() + + # Append result to list of dictionaries + results.append(result) + + # Create FITS table from results + table = self._create_fits_table(results) + + # Create FITS file and append FITS table to FITS file + self._fits = gammalib.GFits() + self._fits.append(table) + + # Return + return + + def _save_fits(self): + """ + Saved phase fit results into a fits file + """ + # Write header + self._log_header1(gammalib.TERSE, 'Save phase FITS file') + + # Create FITS file + self._create_fits() + + # Get file name + outfile = self['outfile'].filename() + + # Log file name + self._log_value(gammalib.NORMAL, 'Phase curve file', outfile.url()) + + # Save to fits + self._fits.saveto(outfile, self._clobber()) + + # Return + return + + def _save_xml(self): + """ + Save phase fit results into one XML model each bin + """ + # Write header + self._log_header1(gammalib.TERSE, 'Save phase XML files') + + # Get file name + outname = self['outfile'].filename() + + # Loop over all phase bins + for phbin in self._phbins: + + # Set file name for phase bin + phname = str(phbin[0])+'-'+str(phbin[1]) + replace_name = '_phase_'+phname+'.xml' + outname_phase = (outname.file()).replace('.fits', replace_name) + outfile = outname.path() + outname_phase + + # Save XML file + self._fitmodels[phname].save(outfile) + + # Log file name + name = 'Phase [%s] file' % phname + self._log_value(gammalib.NORMAL, name, outfile) + + # Return + return + # Public methods def run(self): @@ -264,24 +427,26 @@ class csphasecrv(ctools.cscript): # Get parameters self._get_parameters() - + # Write observation into logger - if self._logTerse(): - self._log('\n') - self._log.header1(gammalib.number('Input observation',len(self._obs))) - self._log(str(self._obs)) - self._log('\n') - - #Store and clear observations to split them in phase + self._log_observations(gammalib.NORMAL, self._obs, 'Observation') + + # Store and clear observations to split them into phase bins orig_obs = self._obs.copy() - #self._obs.clear() - #Dictionary to save phase fitted models + # Dictionary to save phase fitted models self._fitmodels = {} - - #Loop over all phases + + # Write header + self._log_header1(gammalib.TERSE, 'Generate phase curve') + + # Loop over all phases for phbin in self._phbins: - + + # Write time bin into header + if self._logTerse(): + self._log.header2('PHASE '+str(phbin[0])+'-'+str(phbin[1])) + # Select events select = ctools.ctselect(orig_obs) select['emin'] = self['emin'].real() @@ -291,13 +456,13 @@ class csphasecrv(ctools.cscript): select['rad'] = 'UNDEFINED' select['ra'] = 'UNDEFINED' select['dec'] = 'UNDEFINED' - select['expr'] = "PHASE>"+str(phbin[0])+" && PHASE<"+str(phbin[1]) + select['expr'] = 'PHASE>'+str(phbin[0])+' && PHASE<'+str(phbin[1]) select.run() - #Add phase to observation id + # Add phase to observation id for i in range(0,select.obs().size()): oldid = select.obs()[i].id() - select.obs()[i].id(oldid+"_"+str(phbin[0])+"-"+str(phbin[1])) + select.obs()[i].id(oldid+'_'+str(phbin[0])+'-'+str(phbin[1])) obs = select.obs() # If a stacked analysis is requested then bin the events @@ -318,15 +483,18 @@ class csphasecrv(ctools.cscript): # Renormalize models to phase selection # TODO move the scaling from the temporal to the spectral component for model in like.obs().models(): - scaled_norm = model["Normalization"].value()/(phbin[1] - phbin[0]) - model["Normalization"].value(scaled_norm) + scaled_norm = model['Normalization'].value()/(phbin[1] - phbin[0]) + model['Normalization'].value(scaled_norm) - #Store fit model - self._fitmodels[str(phbin[0])+"-"+str(phbin[1])] = like.obs().models().copy() + # Store fit model + self._fitmodels[str(phbin[0])+'-'+str(phbin[1])] = like.obs().models().copy() + + # Create FITS file + self._create_fits() - ## Store observations - #for ob in select.obs(): - #self._obs.append(ob.copy()) + # Optionally publish light curve + if self['publish'].boolean(): + self.publish() # Return return @@ -350,171 +518,57 @@ class csphasecrv(ctools.cscript): # Return return - def saveXML(self): - """ - Save phase fits into one XML model each - """ - # Write header - self._log_header1(gammalib.TERSE, 'Save phase XML') - - # Get out filename - outfile = str() - - # Save XMLs - for phbin in self._phbins: - phname = str(phbin[0])+"-"+str(phbin[1]) - outname = self['outfile'].filename() - outname_ph = (outname.file()).replace(".fits","_ph_"+phname+".xml") - outfile = outname.path()+outname_ph - print outfile - self._fitmodels[phname].save(outfile) - - # Return - return - def save(self): """ Saves results to fits and XML """ - print "SAVING" - self.savefits() - self.saveXML() + # Save results to FITS + self._save_fits() + + # Save results to XML files (one per phase bin) + self._save_xml() # Return return - def phaserv(self): - """ - Return dictionary with best fit models - """ - # Return - return self._fitmodels - - def _get_free_par_names(self): + def publish(self, name=''): """ - Return list of free parameter names + Publish phase curve - Returns - ------- - names : list of str - List of free parameter names. - """ - # Initialise empty list of free parameter names - names = [] - - # Collect list of free parameter names - for par in self._obs.models()[self._srcname]: - if par.is_free(): - names.append(par.name()) - - # Return names - return names - - def savefits(self): - """ - Saved phase fit results to a fits file + Parameters + ---------- + name : str, optional + Name of phase curve """ - # Initialise list of result dictionaries - results = [] - - # Get source parameters - pars = self._get_free_par_names() - - # Loop over time bins - for i in range(len(self._phbins)): + # Write header + self._log_header1(gammalib.TERSE, 'Publish phase curve') - # Get time boundaries - phmin = self._phbins[i][0] - phmax = self._phbins[i][1] + # Continue only if light curve is valid + if self._fits.contains('PHASECURVE'): - # Write time bin into header - if self._logTerse(): - self._log.header2('PHASE '+str(phmin)+'-'+str(phmax)) + # Set default name is user name is empty + if not name: + user_name = self._name + else: + user_name = name - # Initialise result dictionary - result = {'phmin': phmin, - 'phmax': phmax, - 'pars': pars, - 'values': {}} - # Extract parameter values - - #Store fit results - phname = str(phmin)+"-"+str(phmax) - source = self._fitmodels[phname][self._srcname] - for par in pars: - result['values'][par] = source[par].value() - result['values']['e_'+par] = source[par].error() - - # Append result to list of dictionaries - results.append(result) - - # Create FITS table from results - table = self._create_fits_table(results) + # Log file name + self._log_value(gammalib.NORMAL, 'Phase curve name', user_name) - # Create FITS file and append FITS table to FITS file - self._fits = gammalib.GFits() - self._fits.append(table) + # Publish phase curve + self._fits.publish('PHASECURVE', user_name) - # Save to fits - outfile = self['outfile'].filename() - print outfile - self._fits.saveto(outfile, self._clobber()) - - #Return + # Return return - def _create_fits_table(self, results): + def phasecurve(self): """ - Creates FITS binary table containing light curve results - - Parameters - ---------- - results : list of dict - List of result dictionaries - - Returns - ------- - table : `~gammalib.GFitsBinTable` - FITS binary table containing light curve + Return dictionary with best fit models """ - # Determine number of rows in FITS table - nrows = len(results) - - # Create FITS Table with extension "LIGHTCURVE" - table = gammalib.GFitsBinTable(nrows) - table.extname('PHASESPEC') - - # Append time columns "MJD" and "e_MJD" - phmin = gammalib.GFitsTableDoubleCol('PHMIN', nrows) - phmax = gammalib.GFitsTableDoubleCol('PHMAX', nrows) - #mjd.unit('days') - #e_mjd.unit('days') - for i, result in enumerate(results): - phmin[i] = result['phmin'] - phmax[i] = result['phmax'] - table.append(phmin) - table.append(phmax) - - # Create parameter columns - for par in self._obs.models()[self._srcname]: - if par.is_free(): - name = par.name() - e_name = 'e_'+par.name() - col = gammalib.GFitsTableDoubleCol(name, nrows) - e_col = gammalib.GFitsTableDoubleCol(e_name, nrows) - col.unit(par.unit()) - e_col.unit(par.unit()) - for i, result in enumerate(results): - col[i] = result['values'][name] - e_col[i] = result['values'][e_name] - table.append(col) - table.append(e_col) - - # Return table - return table + # Return + return self._fitmodels - # ======================== # # Main routine entry point # # ======================== # diff --git a/cscripts/cspull.par b/cscripts/cspull.par index 61520cb..2f78313 100644 --- a/cscripts/cspull.par +++ b/cscripts/cspull.par @@ -31,7 +31,7 @@ ntrials, i, a, 10,,, "Number of trials" ra, r, a, 83.6331,0,360, "RA of pointing (deg)" dec, r, a, 22.0145,-90,90, "Dec of pointing (deg)" coordsys, s, a, CEL,CEL|GAL,, "Coordinate system (CEL - celestial, GAL - galactic)" -proj, s, a, TAN,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, TAN,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" emin, r, a, 0.1,,, "Lower energy limit (TeV)" emax, r, a, 100.0,,, "Upper energy limit (TeV)" enumbins, i, a, 0,,, "Number of energy bins (0=unbinned)" diff --git a/cscripts/cspull.py b/cscripts/cspull.py index e77857f..56eb1ac 100755 --- a/cscripts/cspull.py +++ b/cscripts/cspull.py @@ -2,7 +2,7 @@ # ========================================================================== # This script generates the pull distribution for all model parameters. # -# Copyright (C) 2011-2016 Juergen Knoedlseder +# Copyright (C) 2011-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -40,7 +40,7 @@ class cspull(ctools.cscript): """ # Set name self._name = 'cspull' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise some members self._edisp = False diff --git a/cscripts/csresmap.par b/cscripts/csresmap.par index 67564f2..ed0c42c 100644 --- a/cscripts/csresmap.par +++ b/cscripts/csresmap.par @@ -37,7 +37,7 @@ emax, r, a, 100.0,,, "Upper energy limit (TeV)" enumbins, i, h, 20,,, "Number of model cube energy bins" ebinfile, f, a, NONE,,, "Name of the file containing the energy bin definition" coordsys, s, a, CEL,CEL|GAL,, "Coordinate System" -proj, s, a, CAR,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" nxpix, i, a, 200,,, "Size of the X axis in pixels" diff --git a/cscripts/csresmap.py b/cscripts/csresmap.py index ff1e7cd..806c2bd 100755 --- a/cscripts/csresmap.py +++ b/cscripts/csresmap.py @@ -2,7 +2,7 @@ # ========================================================================== # Generates a residual map. # -# Copyright (C) 2014-2016 Michael Mayer +# Copyright (C) 2014-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -38,7 +38,7 @@ class csresmap(ctools.cscript): """ # Set name self._name = 'csresmap' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise class members self._resmap = None diff --git a/cscripts/csroot2caldb.py b/cscripts/csroot2caldb.py index 28c3df4..390beaf 100755 --- a/cscripts/csroot2caldb.py +++ b/cscripts/csroot2caldb.py @@ -50,7 +50,7 @@ class csroot2caldb(ctools.cscript): """ # Set name self._name = 'csroot2caldb' - self._version = '1.3.0' + self._version = ctools.__version__ # Initialise application by calling the appropriate class constructor self._init_cscript(argv) diff --git a/cscripts/cssens.py b/cscripts/cssens.py index f16590c..4165895 100755 --- a/cscripts/cssens.py +++ b/cscripts/cssens.py @@ -50,7 +50,7 @@ class cssens(ctools.cscript): """ # Set name self._name = 'cssens' - self._version = '1.3.0' + self._version = ctools.__version__ # Initialise class members self._ebounds = gammalib.GEbounds() diff --git a/cscripts/csspec.py b/cscripts/csspec.py index 64a5ebd..1e82724 100755 --- a/cscripts/csspec.py +++ b/cscripts/csspec.py @@ -2,7 +2,7 @@ # ========================================================================== # Generates a spectrum. # -# Copyright (C) 2014-2016 Michael Mayer +# Copyright (C) 2014-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -43,7 +43,7 @@ class csspec(ctools.cscript): """ # Set name self._name = 'csspec' - self._version = '1.1.0' + self._version = ctools.__version__ # Initialise some members self._ebounds = gammalib.GEbounds() diff --git a/cscripts/cssrcdetect.py b/cscripts/cssrcdetect.py index f053ddd..1cc0ad5 100755 --- a/cscripts/cssrcdetect.py +++ b/cscripts/cssrcdetect.py @@ -2,7 +2,7 @@ # ========================================================================== # Detects sources in a sky map # -# Copyright (C) 2016 Juergen Knoedlseder +# Copyright (C) 2016-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -40,7 +40,7 @@ class cssrcdetect(ctools.cscript): """ # Set name self._name = 'cssrcdetect' - self._version = '1.2.0' + self._version = ctools.__version__ # Set protected members self._models = gammalib.GModels() diff --git a/cscripts/cstsdist.par b/cscripts/cstsdist.par index 438b026..515742d 100644 --- a/cscripts/cstsdist.par +++ b/cscripts/cstsdist.par @@ -36,7 +36,7 @@ ntrials, i, a, 10,,, "Number of trials" ra, r, a, 83.63,0,360, "RA of pointing (deg)" dec, r, a, 22.01,-90,90, "Dec of pointing (deg)" coordsys, s, a, CEL,CEL|GAL,, "Coordinate system (CEL - celestial, GAL - galactic)" -proj, s, a, TAN,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, TAN,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" emin, r, a, 0.1,,, "Lower energy limit (TeV)" emax, r, a, 100.0,,, "Upper energy limit (TeV)" enumbins, i, a, 0,,, "Number of energy bins (0=unbinned)" diff --git a/cscripts/cstsdist.py b/cscripts/cstsdist.py index 2411f7e..f35029a 100755 --- a/cscripts/cstsdist.py +++ b/cscripts/cstsdist.py @@ -2,7 +2,7 @@ # ========================================================================== # Generates the TS distribution for a particular model # -# Copyright (C) 2011-2016 Juergen Knoedlseder +# Copyright (C) 2011-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -41,7 +41,7 @@ class cstsdist(ctools.cscript): """ # Set name self._name = 'cstsdist' - self._version = '1.2.0' + self._version = ctools.__version__ # Initialise some members self._srcname = '' diff --git a/cscripts/cstsmapmerge.py b/cscripts/cstsmapmerge.py index 9b12c8c..d385ee6 100755 --- a/cscripts/cstsmapmerge.py +++ b/cscripts/cstsmapmerge.py @@ -2,7 +2,7 @@ # ========================================================================== # Merge Test Statistic maps # -# Copyright (C) 2015-2016 Michael Mayer +# Copyright (C) 2015-2017 Michael Mayer # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -40,7 +40,7 @@ class cstsmapmerge(ctools.cscript): """ # Set name self._name = 'cstsmapmerge' - self._version = '1.1.0' + self._version = ctools.__version__ # Initialise class members self._files = None diff --git a/cscripts/cstsmapsplit.par b/cscripts/cstsmapsplit.par index b51bb67..bafeafc 100644 --- a/cscripts/cstsmapsplit.par +++ b/cscripts/cstsmapsplit.par @@ -41,9 +41,9 @@ nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, CAR, AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/cscripts/cstsmapsplit.py b/cscripts/cstsmapsplit.py index db736b6..48e8310 100755 --- a/cscripts/cstsmapsplit.py +++ b/cscripts/cstsmapsplit.py @@ -47,7 +47,7 @@ class cstsmapsplit(ctools.cscript): # Set name and version self._name = 'cstsmapsplit' - self._version = '1.2.0' + self._version = ctools.__version__ # Set data members self._outmap = gammalib.GFilename() diff --git a/cscripts/csviscube.py b/cscripts/csviscube.py index 92c9c13..9e3ae2a 100755 --- a/cscripts/csviscube.py +++ b/cscripts/csviscube.py @@ -2,7 +2,7 @@ # ========================================================================== # Compute a visibility cube # -# Copyright (C) 2016 Juergen Knoedlseder +# Copyright (C) 2016-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -50,7 +50,7 @@ class csviscube(ctools.cscript): """ # Set name self._name = 'csviscube' - self._version = '1.0.0' + self._version = ctools.__version__ # Initialise members self._cube = gammalib.GSkyMap() diff --git a/cscripts/csworkflow.py b/cscripts/csworkflow.py index 4ecdb91..98a9939 100755 --- a/cscripts/csworkflow.py +++ b/cscripts/csworkflow.py @@ -2,7 +2,7 @@ # ========================================================================== # Executes analysis workflow # -# Copyright (C) 2016 Juergen Knoedlseder +# Copyright (C) 2016-2017 Juergen Knoedlseder # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by @@ -53,7 +53,7 @@ class csworkflow(ctools.cscript): """ # Set name and version self._name = 'csworkflow' - self._version = '1.2.0' + self._version = ctools.__version__ # Set members self._workflow = gammalib.GXml() diff --git a/doc/source/users/reference_manual/cslightcrv.rst b/doc/source/users/reference_manual/cslightcrv.rst index d4c7509..d6a2f18 100644 --- a/doc/source/users/reference_manual/cslightcrv.rst +++ b/doc/source/users/reference_manual/cslightcrv.rst @@ -75,7 +75,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``xref [real]`` diff --git a/doc/source/users/reference_manual/csphasecrv.rst b/doc/source/users/reference_manual/csphasecrv.rst index ae1e624..8c2d944 100644 --- a/doc/source/users/reference_manual/csphasecrv.rst +++ b/doc/source/users/reference_manual/csphasecrv.rst @@ -9,13 +9,16 @@ Computes phase dependent spectra for a given source. Synopsis -------- -This script computes spectra by performing a maximum likelihood fit -using :doc:`ctlike` in a series of phase bins for pulsars. -The phase bins can be either specified in an ASCII file, as an interval -divided into equally sized time bins. -The format of the ASCII file is one row per time bin, each specifying the -start of stop value of the phase bin, separated by a whitespace. The phase -goes from 0.0 to 1.0 . +This script computes spectra by performing a maximum likelihood fit using +:doc:`ctlike` in a series of phase bins for pulsars. The phase bins can be +either specified in an ASCII file or as an interval divided into equally sized +phase bins. The format of the ASCII file is one row per phase bin, each +specifying the start of stop value of the phase bin, separated by a whitespace. +The phase goes from 0.0 to 1.0. + +On output the script writes the fitting results into a FITS file. The script +also produces one XML file per phase bin that contains the fitting results for +that bin. General parameters @@ -25,10 +28,10 @@ General parameters Input event list or observation definition XML file. ``inmodel [file]`` - Input model XML file. + Input model definition XML file. ``srcname [string]`` - Name of the periodic source in the source model XML file + Name of the periodic source in the source model XML file. ``caldb [string]`` Calibration database. @@ -40,10 +43,11 @@ General parameters Applies energy dispersion to response computation. ``outfile [file]`` - Name of the XML output file. The phase intervall will be automatically appended to the name. + Name of the XML output file. The phase interval will be automatically + appended to the name. -``(tbinalg = LIN) [string]`` - Algorithm for defining time bins. +``phbinalg [string]`` + Algorithm for defining phase bins. ``phbins [integer]`` Number of phase bins. @@ -63,7 +67,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``xref [real]`` @@ -81,9 +85,13 @@ General parameters ``binsz [real]`` Pixel size (in degrees/pixel). + Standard parameters ------------------- +``(publish = no) [boolean]`` + Specifies whether the phase curve should be published on VO Hub. + ``(chatter = 2) [integer]`` Verbosity of the executable: chatter = 0: no information will be logged diff --git a/doc/source/users/reference_manual/cspull.rst b/doc/source/users/reference_manual/cspull.rst index 7be3549..15190b9 100644 --- a/doc/source/users/reference_manual/cspull.rst +++ b/doc/source/users/reference_manual/cspull.rst @@ -71,7 +71,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``emin [real]`` diff --git a/doc/source/users/reference_manual/csresmap.rst b/doc/source/users/reference_manual/csresmap.rst index 35d6bd4..bdbed8e 100644 --- a/doc/source/users/reference_manual/csresmap.rst +++ b/doc/source/users/reference_manual/csresmap.rst @@ -76,7 +76,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``xref [real]`` diff --git a/doc/source/users/reference_manual/cstsdist.rst b/doc/source/users/reference_manual/cstsdist.rst index 9884871..90d37bd 100644 --- a/doc/source/users/reference_manual/cstsdist.rst +++ b/doc/source/users/reference_manual/cstsdist.rst @@ -78,7 +78,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``emin [real]`` diff --git a/doc/source/users/reference_manual/cstsmapsplit.rst b/doc/source/users/reference_manual/cstsmapsplit.rst index 27417d9..0660311 100644 --- a/doc/source/users/reference_manual/cstsmapsplit.rst +++ b/doc/source/users/reference_manual/cstsmapsplit.rst @@ -80,16 +80,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of map centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of map centre (J2000, in degrees). -``proj [string]`` - Projection method. - - Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctbin.rst b/doc/source/users/reference_manual/ctbin.rst index e0bb555..7314bdb 100644 --- a/doc/source/users/reference_manual/ctbin.rst +++ b/doc/source/users/reference_manual/ctbin.rst @@ -65,15 +65,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of cube centre (J2000, in degrees). -``proj [string]`` - Projection method. - Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctbkgcube.rst b/doc/source/users/reference_manual/ctbkgcube.rst index c6e2a27..dc9e34a 100644 --- a/doc/source/users/reference_manual/ctbkgcube.rst +++ b/doc/source/users/reference_manual/ctbkgcube.rst @@ -105,15 +105,15 @@ General parameters ``coordsysL [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of cube centre (J2000, in degrees). -``proj [string]`` - Projection method. - Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctedispcube.rst b/doc/source/users/reference_manual/ctedispcube.rst index bc75285..02f56c8 100644 --- a/doc/source/users/reference_manual/ctedispcube.rst +++ b/doc/source/users/reference_manual/ctedispcube.rst @@ -88,15 +88,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of cube centre (J2000, in degrees). -``proj [string]`` - Projection method. - ``(migramax = 2.0) [real]`` Upper bound of ratio between reconstructed and true photon energy. diff --git a/doc/source/users/reference_manual/ctexpcube.rst b/doc/source/users/reference_manual/ctexpcube.rst index 5a4dc6a..cbabc3d 100644 --- a/doc/source/users/reference_manual/ctexpcube.rst +++ b/doc/source/users/reference_manual/ctexpcube.rst @@ -95,15 +95,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` - Declination / Galactic latitude of cube centre (J2000, in degrees). - -``proj [string]`` - Projection method. - + Declination / Galactic latitude of cube centre (J2000, in degrees). + Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctmapcube.rst b/doc/source/users/reference_manual/ctmapcube.rst index f08ef8d..5ef17ec 100644 --- a/doc/source/users/reference_manual/ctmapcube.rst +++ b/doc/source/users/reference_manual/ctmapcube.rst @@ -61,15 +61,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of cube centre (J2000, in degrees). -``proj [string]`` - Projection method. - Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctmodel.rst b/doc/source/users/reference_manual/ctmodel.rst index 7a356c1..bdf3e76 100644 --- a/doc/source/users/reference_manual/ctmodel.rst +++ b/doc/source/users/reference_manual/ctmodel.rst @@ -125,15 +125,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of image centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of image centre (J2000, in degrees). -``proj [string]`` - Projection method. - Standard parameters ------------------- diff --git a/doc/source/users/reference_manual/ctpsfcube.rst b/doc/source/users/reference_manual/ctpsfcube.rst index abe25fa..03f159f 100644 --- a/doc/source/users/reference_manual/ctpsfcube.rst +++ b/doc/source/users/reference_manual/ctpsfcube.rst @@ -87,15 +87,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of cube centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of cube centre (J2000, in degrees). -``proj [string]`` - Projection method. - ``(amax = 0.3) [real]`` Upper bound of angular separation between true and measued photon direction (in degrees). diff --git a/doc/source/users/reference_manual/ctskymap.rst b/doc/source/users/reference_manual/ctskymap.rst index f4c3591..0c90583 100644 --- a/doc/source/users/reference_manual/ctskymap.rst +++ b/doc/source/users/reference_manual/ctskymap.rst @@ -65,7 +65,7 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). -``proj [string]`` +``proj [string]`` Projection method. ``xref [real]`` diff --git a/doc/source/users/reference_manual/cttsmap.rst b/doc/source/users/reference_manual/cttsmap.rst index 401125c..4a74e4b 100644 --- a/doc/source/users/reference_manual/cttsmap.rst +++ b/doc/source/users/reference_manual/cttsmap.rst @@ -78,15 +78,15 @@ General parameters ``coordsys [string]`` Coordinate system (CEL - celestial, GAL - galactic). +``proj [string]`` + Projection method. + ``xref [real]`` Right Ascension / Galactic longitude of map centre (J2000, in degrees). ``yref [real]`` Declination / Galactic latitude of map centre (J2000, in degrees). -``proj [string]`` - Projection method. - ``(binmin = -1) [integer]`` First bin to compute. diff --git a/doc/source/users/reference_manual/index.rst b/doc/source/users/reference_manual/index.rst index 50e8218..f81202a 100644 --- a/doc/source/users/reference_manual/index.rst +++ b/doc/source/users/reference_manual/index.rst @@ -43,13 +43,14 @@ cscripts cscaldb --- Lists available instrument response functions csebins --- Generates energy boundaries for stacked analysis - cslightcrv --- Computes lightcurve + cslightcrv --- Computes light curve csmodelinfo --- Shows model container content csmodelmerge --- Merges several model containers into one file csmodelselect --- Select models from model definition file csobsdef --- Generates observation definition file csobsinfo --- Shows observation container content csobsselect --- Select observations from observation definition file + csphasecrv --- Computes phase curve cspull --- Generates pull distribution csresmap --- Generates residual map cssens --- Computes CTA sensitivity diff --git a/src/ctbin/ctbin.par b/src/ctbin/ctbin.par index a55b087..394479e 100644 --- a/src/ctbin/ctbin.par +++ b/src/ctbin/ctbin.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctbin.par # -# Version: 1.1.0 # # # # Author: Juergen Knoedlseder # # Jurgen.Knodlseder@irap.omp.eu # @@ -37,9 +36,9 @@ nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" coordsys, s, a, CEL,CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, CAR,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/src/ctbkgcube/ctbkgcube.par b/src/ctbkgcube/ctbkgcube.par index 7b67403..0922af8 100644 --- a/src/ctbkgcube/ctbkgcube.par +++ b/src/ctbkgcube/ctbkgcube.par @@ -41,10 +41,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/src/ctedispcube/ctedispcube.par b/src/ctedispcube/ctedispcube.par index 602af07..3fcd398 100644 --- a/src/ctedispcube/ctedispcube.par +++ b/src/ctedispcube/ctedispcube.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctedispcube.par # -# Version: 1.2.0 # # # # Author: Maria Haupt # # maria.haupt@desy.de # @@ -40,10 +39,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 10,,, "Size of the X axis in pixels" nypix, i, a, 10,,, "Size of the Y axis in pixels" binsz, r, a, 1.0,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Angular separation binning parameters diff --git a/src/ctexpcube/ctexpcube.par b/src/ctexpcube/ctexpcube.par index d03845d..f938860 100644 --- a/src/ctexpcube/ctexpcube.par +++ b/src/ctexpcube/ctexpcube.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctexpcube.par # -# Version: 1.2.0 # # # # Author: Juergen Knoedlseder # # Jurgen.Knodlseder@irap.omp.eu # @@ -40,10 +39,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/src/ctmapcube/ctmapcube.par b/src/ctmapcube/ctmapcube.par index 08c82bd..c233337 100644 --- a/src/ctmapcube/ctmapcube.par +++ b/src/ctmapcube/ctmapcube.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctmapcube.par # -# Version: 1.0.0 # # # # Author: Juergen Knoedlseder # # Jurgen.Knodlseder@irap.omp.eu # @@ -37,9 +36,9 @@ nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" coordsys, s, a, CEL,CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR,AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, CAR,AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/src/ctmodel/ctmodel.par b/src/ctmodel/ctmodel.par index e4dcec1..03a3795 100644 --- a/src/ctmodel/ctmodel.par +++ b/src/ctmodel/ctmodel.par @@ -54,10 +54,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Standard parameters diff --git a/src/ctpsfcube/ctpsfcube.par b/src/ctpsfcube/ctpsfcube.par index d12bc49..6aab813 100644 --- a/src/ctpsfcube/ctpsfcube.par +++ b/src/ctpsfcube/ctpsfcube.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctpsfcube.par # -# Version: 1.2.0 # # # # Author: Chia-Chun Lu # # c.lu@mpi-hd.mpg.de # @@ -40,10 +39,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 10,,, "Size of the X axis in pixels" nypix, i, a, 10,,, "Size of the Y axis in pixels" binsz, r, a, 1.0,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Angular separation binning parameters diff --git a/src/ctselect/ctselect.par b/src/ctselect/ctselect.par index af15510..6240c80 100644 --- a/src/ctselect/ctselect.par +++ b/src/ctselect/ctselect.par @@ -5,7 +5,6 @@ # ------------------------------------------------------------------------- # # # # File: ctselect.par # -# Version: 1.2.0 # # # # Author: Jurgen Knodlseder # # Jurgen.Knodlseder@irap.omp.eu # @@ -32,8 +31,8 @@ tmin, r, a, 0.0,,, "Start time (CTA MET in seconds)" tmax, r, a, 0.0,,, "End time (CTA MET in seconds)" emin, r, a, 0.1,,, "Lower energy limit (TeV)" emax, r, a, 100.0,,, "Upper energy limit (TeV)" -expr, s, h, "",,, "Additional selection expression" -usethres,s, h, "NONE",NONE|DEFAULT|USER,, "Energy threshold type" +expr, s, h, ,,, "Additional selection expression" +usethres,s, h, NONE,NONE|DEFAULT|USER,, "Energy threshold type" # # Standard parameters diff --git a/src/ctskymap/ctskymap.par b/src/ctskymap/ctskymap.par index 6608657..e53f0ed 100644 --- a/src/ctskymap/ctskymap.par +++ b/src/ctskymap/ctskymap.par @@ -35,7 +35,7 @@ nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" -proj, s, a, CAR, AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" diff --git a/src/cttsmap/cttsmap.par b/src/cttsmap/cttsmap.par index 5a2c4e4..de9d3ad 100644 --- a/src/cttsmap/cttsmap.par +++ b/src/cttsmap/cttsmap.par @@ -40,10 +40,10 @@ usepnt, b, h, no,,, "Use pointing instead of xref/yref parameters?" nxpix, i, a, 200,,, "Size of the X axis in pixels" nypix, i, a, 200,,, "Size of the Y axis in pixels" binsz, r, a, 0.02,,, "Image scale (in degrees/pixel)" -coordsys, s, a, "CEL", CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +coordsys, s, a, CEL, CEL|GAL,,"Coordinate system (CEL - celestial, GAL - galactic)" +proj, s, a, CAR, AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN,, "Projection method" xref, r, a, 83.63,0,360, "First coordinate of image center in degrees (RA or galactic l)" yref, r, a, 22.01,-90,90, "Second coordinate of image center in degrees (DEC or galactic b)" -proj, s, a, "CAR", AIT|AZP|CAR|MER|MOL|STG|TAN,, "Projection method" # # Parameters for splitting and speed purpose diff --git a/test/data/phased_events.fits b/test/data/phased_events.fits deleted file mode 100644 index 1c2b9d4..0000000 --- a/test/data/phased_events.fits +++ /dev/null @@ -1,577 +0,0 @@ -SIMPLE = T / file does conform to FITS standard BITPIX = 8 / number of bits per data pixel NAXIS = 0 / number of data axes EXTEND = T / FITS dataset may contain extensions COMMENT FITS (Flexible Image Transport System) format is defined in 'AstronomyCOMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H COMMENT FITS (Flexible Image Transport System) format is defined in 'AstronomyCOMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H END XTENSION= 'BINTABLE' / binary table extension BITPIX = 8 / 8-bit bytes NAXIS = 2 / 2-dimensional binary table NAXIS1 = 36 / width of table in bytes NAXIS2 = 3813 / number of rows in table PCOUNT = 0 / size of special data area GCOUNT = 1 / one data group (required keyword) TFIELDS = 8 / number of fields in each row TTYPE1 = 'EVENT_ID' / label for field 1 TFORM1 = '1J ' / data format of field: 4-byte INTEGER TZERO1 = 2147483648 / offset for unsigned integers TSCAL1 = 1 / data are not scaled TTYPE2 = 'TIME ' / label for field 2 TFORM2 = '1D ' / data format of field: 8-byte DOUBLE TUNIT2 = 's ' / physical unit of field TTYPE3 = 'RA ' / label for field 3 TFORM3 = '1E ' / data format of field: 4-byte REAL TUNIT3 = 'deg ' / physical unit of field TTYPE4 = 'DEC ' / label for field 4 TFORM4 = '1E ' / data format of field: 4-byte REAL TUNIT4 = 'deg ' / physical unit of field TTYPE5 = 'ENERGY ' / label for field 5 TFORM5 = '1E ' / data format of field: 4-byte REAL TUNIT5 = 'TeV ' / physical unit of field TTYPE6 = 'DETX ' / label for field 6 TFORM6 = '1E ' / data format of field: 4-byte REAL TUNIT6 = 'deg ' / physical unit of field TTYPE7 = 'DETY ' / label for field 7 TFORM7 = '1E ' / data format of field: 4-byte REAL TUNIT7 = 'deg ' / physical unit of field TTYPE8 = 'PHASE ' / label for field 8 TFORM8 = '1E ' / data format of field: 4-byte REAL EXTNAME = 'EVENTS ' / name of this extension DSTYP1 = 'TIME ' / Data sub-space type DSUNI1 = 's ' / Data sub-space unit DSVAL1 = 'TABLE ' / Data sub-space value DSREF1 = ':GTI ' / Data sub-space reference DSTYP2 = 'ENERGY ' / Data sub-space type DSUNI2 = 'TeV ' / Data sub-space unit DSVAL2 = '0.1:100 ' / Data sub-space value DSTYP3 = 'POS(RA,DEC)' / Data sub-space type DSUNI3 = 'deg ' / Data sub-space unit DSVAL3 = 'CIRCLE(83.63,22.01,5)' / Data sub-space value NDSKEYS = 3 / Number of data sub-space keys CREATOR = 'GammaLib' / Program which created the file TELESCOP= 'CTA ' / Telescope OBS_ID = 0 / Observation identifier DATE_OBS= '2000-01-01' / Observation start date TIME_OBS= '11:58:56' / Observation start time DATE_END= '2000-01-01' / Observation end date TIME_END= '12:28:56' / Observation end time TSTART = 0.0000000000E+00 / [s] Mission time of start of observation TSTOP = 1.8000000000E+03 / [s] Mission time of end of observation MJDREFI = 51544 / [days] Integer part of time reference MJD MJDREFF = 5.0000000000E-01 / [days] Fractional part of time reference MJD TIMEUNIT= 's ' / Time unit TIMESYS = 'TT ' / Time system TIMEREF = 'LOCAL ' / Time reference TELAPSE = 1.8000000000E+03 / [s] Mission elapsed time ONTIME = 1.8000000000E+03 / [s] Total good time including deadtime LIVETIME= 1.7640000000E+03 / [s] Total livetime DEADC = 9.8000000000E-01 / Deadtime correction factor TIMEDEL = 1.0000000000E+00 / Time resolution OBJECT = ' ' / Observed object RA_OBJ = 0.0000000000E+00 / [deg] Target Right Ascension DEC_OBJ = 0.0000000000E+00 / [deg] Target Declination RA_PNT = 8.3630000000E+01 / [deg] Pointing Right Ascension DEC_PNT = 2.2010000000E+01 / [deg] Pointing Declination ALT_PNT = 9.0000000000E+01 / [deg] Average altitude of pointing AZ_PNT = 0.0000000000E+00 / [deg] Average azimuth of pointing RADECSYS= 'FK5 ' / Coordinate system EQUINOX = 2.0000000000E+03 / Epoch CONV_DEP= 0.0000000000E+00 / Convergence depth of telescopes CONV_RA = 0.0000000000E+00 / [deg] Convergence Right Ascension CONV_DEC= 0.0000000000E+00 / [deg] Convergence Declination OBSERVER= 'string ' / Observer N_TELS = 0 / Number of telescopes in event list TELLIST = 'string ' / Telescope IDs GEOLAT = 0.0000000000E+00 / [deg] Geographic latitude of array centre GEOLON = 0.0000000000E+00 / [deg] Geographic longitude of array centre ALTITUDE= 0.0000000000E+00 / [km] Altitude of array centre EUNIT = 'TeV ' / Energy unit EVTVER = 'draft1 ' / Event list version number END ?(BBA= -c:H=@BAs=Pm=R?]D(@BA> =yL?dFB@,?B- A>B?@`BAZ>"SS1?.>OgL@B!Ai=U?%@lhBDA4> )|;9?n@"Bw A^>4=EY=Œ=?NF @# ^B -A/i=>l=%?F= @$kbbB3A=B=Fݟ.>4Ѐ @'>B2A7 ->la| @*+BHAW?=T=<1|?Weo@-dB6iA">O@3xI>@1.~ BA>ώ섚>E&@2$B<AV{>u=@>\=@2iJ5BMAI>$P`<>*@5YJB\;Ak=`=.8=>T?lS@7hBXA7`>,ܽ=u>Q(@8ezBmA>(="=_?y'@<B)A>=6a>ݫv@f<B?y@@BaAA=[b$R?@A2BjA>/鼏fa=l>?,Ԁ@ABLA>>{U=|>X<þ>"8@CA!B.A (>.<k>@C9BoA4=m=Y=?CE@DaB6Aު>V=/==,&@EB/*AX>3=^lѷ>ن@FnBQA&>FU>ϑ<-?V@GaB3AL>(\e?ym @GWB'A=ٽ=ι?OÀ!@G!BUAc=$У=?` "@H@QQB A> -/=?1M#@IBAN>v>5>$@L?.B#bAt>P={g*>L%@MB,Aj=t=,Խ#M?_&@MRBwA(>BK>!/V>b'@NBSAs=9=>o =]?Q[)@N B54A=Pie(8>ù*@O5MB#AK>e>TZ(> +@OxDB Aw>2:ғ?;<,@P -@B'AJ>ṘG-?a"-@PMB3A=ܼ=W>=<.@PX[eB1@Q"B=A==?!h>Z2@Q7@BA*>pξU?=3@R BA>C9Dj?]4@RbxI@B-A!>;ú?=5@R:BA=X:R=ip6@S)BwA~>5q=:?)x7@SE~B1cA%==< 9?3L,8@S`8Bu'A>c=f޷=?QƲ9@SoBfbA4>Jp;@SBAq=}M=:1m><@T(6SBAM=铗<ܽ"?f=@T,aBE?AoJ>ݽ10;q?R4r>@TiЀB[A0>5|==;7?l?@UBGIAx>.J< @$?R@@V}BA>U>6O>b2?L%$A@W@BAO=y޽n7La?<9}B@W\@BAG=&(Jm">&C@XBAl> -q5[t?e|D@X=0B5A>7ϧ;g&&?F$ՀE@XxHBjwAS>u ;==^0F@X±BA>5Fٵ=?'G@Y~DBD0UI=;F?5=H@Y@BAr><==+?uI@YBfA =zX==pJ@Y0BCA=ܽ1;*>hK@Z tB}A=Hq;=e?;9L@Z)@BAl> ݽd>~U?l6M@[?B*AT}>C= W0Cb?iJˀN@[B@B?T,<"b%=[O@[(B8bA8=j"=ZYP@[zBmA\>-;=?%Q@\Bw*A>m;3 ==;|PR@\bB\QA>CU=?I?KS@\DBQ0A5=8-^<+>PT@\|B0:A =6nҼ?*$U@\PBiAF>u==Ȱ>\V@]w5B.A7>+NAԽS? OW@]1S@B7A=<=-ؼ>? QX@^\@BDAA> ef;I9 >YހY@^)BFAr> -w;J>ƴĀZ@^mBcAQZ=<y>"?#c[@^'B%AQ}>8G<Tm?Y[-\@^wB=A->&?Wˀ]@_GBTA;>X=?4s`@_oQB|'A>3O=՞=Ѐa@_gB6A>8z'?nb@`) BA>ƞQv> c@`!B`B tAF>+Kd@`)BCA =;$:=?~re@`@BisA>i=q'=(+>Of@`> B=A>#?sg@`h BvAX=ޢi=2?=z?wABh@`mBV=A -m>\q=>i@`BLA>;iCŽR?A"uj@aÇ?B?A#>?{;_>?rk@b @BkA͸>o r[=?8ɀl@b0eBT%AW=꽼h=?6vm@bf/ByZA> &:j=F~>#܀n@bmKBCA>!:T:ua=5o@bvw1BbAq> =ne>؈|p@bzdBAn><>9j?\熀q@b~eB^A> b>QRD>-i?u½r@b` BAV>7@N=ݑ =s@br`BAZ=`2Y=At@cB0A>'>Pu@cBHXA:T='*<+>v@c>\`BaA>0͜(>w0w@c` BMdA=>@<?Byx@c BmA->ü |=[?|y@dBkA39>7=>z@d~V B&A==iiL>^V{@ehŠBADy=ևKŽb/>@+Ѐ|@e2HBgA==齧`w?€}@e;T(@BKA>,q<%=Ѐ~@eF BJgA>85hWj60@es@B%A><=z W *?L̀@e| B]A:=VyO=H(?]@ecB~`A=Nb=m-=ݴJ=M@e^$BfA7>Ě=>p@e0BA>G潐?rH@f3+(BGA>?_0=i|<R>@f( Bt A>44ܪ=|?>€@f%>BdkA>

@f`B0LA4=l)?#d@fy B&.Ab>";}߽R{=?΀@fq1 BR/Ar>/I= zO K3>Uf@gM{ B-Ao=2E=f?t@gUB$A=b#>@g@y B A#O>%;Kh>i@gg BEAUT>,w';?|H̀@guBA=#'=̽?=@gk3BM;A@>D~ս:.t@hBkAAe>%tV=y=C> 3<@h&UB1FA>>ս|ѥF >^@hA6BKAD>fk5!ս>^x @h@BA~>\>^?Bր@hB:A=4Lm |Gx=[ @i$)BJeA>)'ؾ D9=B@iS BA>ŌG2T?{@iYtP@BPA>?8=<0{= @idFB.A=cmh>p@it}B4mAJ=T9*@i"BjAX>">5_=G=>B@i`B.$Ay>Ffjt/?N@i`BKA>sf<.?L-@j@BxnAr>==ǯ>@j[BAZ>+={&e?e/@jvnB#7A+> > ch@=Y`@jDӀBcAr=䧅C=wo=>@jB\A0>= =C?~΀@j;B}A>7׸;s=я?9@k B1A>H<=t-)?@k#KB%AGo>sblZ>Ͳ@k BA=I6(ޕ?2Gu@kĹ7B*xAQ>w<2>@kBAk=V=/ϾJr>gր@k@B)(A>>eḀ@lB$BAO>L6]׽=zP@lD BvA>"<=^>b@laݕB?AH`>$<ϔL}?'р@lBCA>=`@l|B(A=&n=Gý@=@ls`B[A=@a/=:N\?N@miBJA#> A_;?.Ǹ=+E>϶@n -"1B>At>%=> B8?jZ}@n.pހB>AM> Z[M>$@n>B&A ->"q6HآL>oҘ@na@BPA= w+O>wjp@nC BGA> C3e<?uÀ@n@B|A#X>)ꤽ F?Iր@o+B`A>Ab=_z> @oB,A=l=ЉT?>}@pAkpBR^Af>I4=$ -l@p$t_PBDTAD>G <ϗ;R>ڀ@p5.!BrAz>]9=R?pހ@p7 7B A==w1y>>j@pc_p0B[Aj>t=,=9&>@p~~BA> -=jU>.@p/РB+ A=ǽƱ.y?6@p(@BwA>) >O=ĥ=@pIBIUA2E>+ ׃M=]tl>pP@pIBdA1=2=az={|?i[@q\BBNA=CQt?2!@q}X@BaUAH=34>DiU=c>|O$@q8PBeA==}=<'@qB6A>,μ}n?/@qfPB$AB=dy0So="><:@r k<`B[Ahy=/&??x@r-/EBBA=Y=vKN9F?E?@r1aB6A.=*Y o`@r73BaAy>-O1=e[?C!@rVJBRNA!=ֺ׾K<Ȑ>O2@rjB\5A+>6z%<;[=>1?@rBڐB$A>( \">AȀ@rVBSAG=եxI@rJQBMyA|>-<,>|@rN@BAי>?Լ>?.@r\B -Ab={&?n@r@B67A=e=UC>@sȃB A>=}J>€@suPB+A>;-=b*?@s ŀBLA>St= ?<V?:T@s;B{KAt>R=Ҁ?}@s!7BIA=!~@s9@BlA6A>9+/?uc@sYKBInA'-=x> Y` Cq K=B?RQI@s͍BV;A>0kZ=T=̑8@sԮB5ZA"F=T48?t⒀@seBAy>!һ;">z>ӌ@sQB?A>|ڻA?}'@t B,1A>Hz9%?lo@ti -BIAΠ>0 ʽDua=\F=? @t!ֶB%A=ݺY½U=>@t^`BA=6=ZM>d@thc;0B+A.>Jbn.?78@tkk`B A&>0L0<nDz=-Ѐ@tv5pBAz> - =K9=?O9@tO`B9Aq==a>ƀ@t̊BA> "T=W^>Zx@tXBhDA>==>ዀ@uNtB)\A>'F:NM;%t?J,@u֚`BrAc3===8>&@u6BAK>*)<W?O@ux BRA>+ @y<x?]9π@u"0BCA>=zz;z_>e@u8\BhAj=>q4=>ql@upB(IAt(> aNTB >@u;@B6VAF(>*j>μ!+>r@ud BA|>*=d̽=h@v@BA<=3-m=C?pw@vHZ!B ?A|=}w_v'>Ȁ@vR"BA6= =U%8?U@vXFB,ZA>;=Iƽ$? -@v]z.BJAb=ؘ` " ,Gz@v>`BYAP> z<=,?J@vB>A(>:$߻\>c@vUpBfUA>Ͻ=?cր@vBAkA^>Ξ"8ֻs?cXA @vxgB)A->GG> 6N>B -@vBktAs>PŽ=N?@N @vw/BZApd>< ӽ/= @w B'PAo=`=I>P @wd)BhcA>HV=G>>Z@wz BC.A5==:?*@w!TBn4꼲px=;.9<>i@xZ BAen=Rӽ!?ڀ@xMKB0AU\=D9?>&@xM= B)AA>E۪8?|w@xyBA*=J<1@>2>L@xBKA=ܽ}i5@xBMAͦ>)>\2x@y*0BJnA>G~%A?@y- B%A&> =>F]@y6ٰBtzAh>NV=)=>𼸀@y66B A==<="?7@yj(B[A.>Y$=H@yF/@BhAE+>5<=r>*Ȁ@y B$A> 6;=w^Ou?#@yBzAH=u/>F=ІF>& @z2P@B"WA4>-bdo>!@zD<BiA'=}Ļ I=m|>"@zqB'AYR=_= eG_d=#@zz4Bp.A>#v~=C >Q%$@zKB]cA>'V:Nҫ=F<~@%@zق0BA=r=mi=&@z߷*BqA=6==(?B'@z BuA=LCW0=?7H(@z+B"Ax>HO;vjy?T)@{1pB7A>7.?HtZ*@{SBCSA>"#:?Td+@{BAr>3Dwyf>Fʀ,@{+mpB1AoP>&e)Q>Ԁ-@{E0BWA:=ʽR{p=@?6.@{%BA>7%= r>j/@{RB5A>4%h 5?ha0@{y0BK4AA=Z

8(1@{ZB.A:>}G?T2@|`B7 A>0)W?Q 3@|BpA"==g=I?(!؀4@|wBA>JG= ?T55@|w>0BKA >/lh<%9?3O6@|"BNAGO>k5-n8@}-BnAZ>P= 8d>:9@} -< BX)A\R>,T== .=:@}9@BGA>62>>p>m;@}K3PB{A'=?=8=`=X<@}OõBZA > fز=59?>;=@}[ BDAO>k<;9E?R:;>@}jPB4zAv>{j<`?xǀ?@})PB\A>%B \=D4>C@@~]PB%AF>צ;~d)>B@~NW@BA>=y=ׯ?b C@~R jBWGA<>4=?D@~Y@BLCAv>L/=D<=(E@~_BVA*>-<2 -=?g9F@~gcpB)A>9_=iF8R>G@~iHPB$A>BK=Rٽl>H@~N `BA_>~r [W?݀I@~jBA=戏-?J@~wBKA/=5=4<~$>*ȀK@*RpB AFp>' }<NY>%L@0KBoYAXk>V== >s)M@3.B A'@=> fx>ϰN@< RBFAN>FYp<=O@<PB.A&\> l<?4P@F2pB=AY>w @?hQ@PB -A=旣n"~>R@DpB7A>=_2,s >.S@BnkA>GT?=0=?i -T@RB=A>>0% ?zU@2BY A>6N/="=&>(V@ BEA>3?o=;dz>+t W@vB2A*>IK8]>0X@'B0hA>$n/Ƽ>۠Y@BWAP>P=H=4i>1ƀZ@B=At>q=@ >dƀ[@ʸB0A >H,[i>|}(\@.0BA >#e\(=Ꝼ>?}T]@# BhA==>=.?eD^@,B^A>LHԼ͡=L,#?SE__@.B%}A">1`H;nW?hb(`@30BHAK>8n<8<:N>va@=G7(B_*A>&o==T?=Tb@=}BT=A=ж .G=%>Rc@_hB5IA>%m>\6ZĚ;πd@{BRAW>,d<2>e@)#(BA=~.??Rf@r8BE~A>4;;?ހg@.0B-AQK==Zh@&`BAר>6g>#<2@i@J8B&A'> #ĽLa J= j@9BsAl5>oTY=>oeĀk@BWVAW=^q=0>ײl@vTHB=DA>11?:m@jFpBLA>9;c<?K$n@ 0BA>mӼL>J>.o@5B%:AV> 1=RYv?ep@G/0BAi=ޙ=*+eT>q@B]A =ڻj=D1>6r@1UxBA>D9p@&?Ss@7ۻB*{A D=1=۽2q<t@>A(BfKA3><<|`=~?׀u@>8B3AA=ͿW=p/0>]v@@BA=H >W >w@AC BbAl>&/==mB?QDx@AhB!Az>bq;Hy@E)B#A>CWme>k0z@T~Bm+A>$<=&=?/{@o9pBDA|>~e=;x>h|@rfB Ah>I17X?Z}@tfxBAA~N>,ӽ-QU~>E>~@u(B'AC>GԽdgE?Q@@z%BdA>Kx=|l?>M@{QXB)Ao1=i=5sr7A?L*Հ@`BKA=֓=W<`>^U|@܎xBKJA>CőRAH{5U<;->^P@-BCA >> ݹ;&&?*3@C@BxzAa==g==)}@*SPBheA>'=g =A)?ˀ@BAQ>@y,.=g?OD@ң}pBPA+^>z <џ8>!Ԁ@}0B +ApR=XBp?uU@8BHA>1J1<>?i@,Bb!AY)>} K=jJ>\@" hBQxA=K <->@)hBz+Av?>:=C=%p>[@M@Bf!A8 >15r@O&sBVAKs>>}=?WJ@YB_AO>}x<~=W>i@\BR)A;A>/<@nhB;:Ai>;=>Ya?{@oBsA=>$$=>Qؤ@"0BxA=푶=߿?@dB>A s=r=t>Ԁ@%BAZ>@= 9)><^@ԞBaOA^=ѾJ=;=c>fd@FB(JA,]>KI> Bl>V>@ާhB+MA=J@d ,u<տ @BBA><߼zy?V3@9BAֆ> i>L@{BA(=X=Ł?\Q@8BmA7><0>b@ 1*B:;A>M=)v?^ȁ@'BAA \>)b?Db@BFA=a= ,;>N,@ BvAQ=8}PǮ=fM?$+?@"B,`A>^8T$K?i@*XBA=BBھvPB?@/=pB+QAt=5=@",X;?Cq@0EBNAA(=Ǫ@=wB+A > {*{>∀@JB*zALX>v=2?L@QBZXA~>=t=0N)?й@ZBKDA=tH<.?~p@m!B;A>Ҿ%6TQS"> k@w9AB_A>2v==V??ŀ@HB?A1>g_=c>Tp@BMA>cԽ}T<=c?r@BA>7Bˋ>6?iZy@ȰPBF@Af=䗨;/N>ȡ@?B<A+==@2>㊀@HB#A>!#i>q@PB/AP>0J<@B?SP@oXB#A>C= de=ċȀ@nBAm> 4 =2|>2C=P@]B0A>@YBy?' @,BHA$=ߥ<3k >ي*@+HB?A>,2ʼU4ˊ?{@`BHA>F<9=?V΀@BWA#> @(B9Au=@;W?w0p@)B[A7==>0>^(@27 XB=bAt>6!/Ѽ?@3\>(BQ A=$\=l7G<äY=Va>K@9Q7HBiAxL=lc=H129>€@9x;B9e9=@BBUQA͈= ,= ( ?l&Ā@DRcB|A|>4A=PR=؇>Ā@P~xB,A#>ސ; j>< @TڐBkA > W8=g?i@`B A6>`@G>|B?)@i_\Bd?AZ=^x=z}?j@lPBtA">"W==7?oڀ@lBPUAv>(J=C/<ۻ@pBXlAˆ=h#K="0?e@XBNAb= =qσ0BA>Ap>&=n*?\+΀@" B:A>"EhE?U@AB9A`=Ƽ?c@íB-A>\ƽ=8@#?BCA(>@Zn:]O_?y2@BsA>rK&f==0@b Bb%AL.>s<=j7 ?Aʀ@,B>A

A#?M"@B'OA=>oȽI9?dV7@-BEA>5{P*;>\ʀ@#BuA =;"1=- =A @ 8B^A> !=l=Q%>@DPBA/>0.>T>Pop>z@B KAo>' -~Hp?"q@{BA3o=J<>4m?h7Ȁ@BA=,=*@DB9A>=J~X?1D@S2B,A>(u #-?5 z@Y1B̤A>F@ǾZ%>m@taBaA>6@Ƃ=fh?V@v6BO#A>F wY<>9Ā@wB[Ah==(/=:]>P@6BA&1=\}<-?`;@gB`Aɶ=A=@=Z\>;;@BQAm>) -@$BdfA =C}={3:>*@yB$(==C?A@@Е0BVA>@¼>=U?.S@ӛBoAg>JGX=L>&@BcA=^/>CL=uDl>zW@TB]lA=@=2=G!?g@BAN>$\? > r>K@+ּ= :?@)@HdBA =+VW*_>̀@KgިB<AH> <Έ@=6@P&B7Au=Gh>8K =b@Z0BDA)>BDuG;xSC=Uh@aPBEA>Gּ`;3[?xɱ@u5BqTA}=k<=d=_=zQ@* hB9A>*==_?IL@>B#&A}l= - i??@͉@B)hA>=g:gt? "@ԟOB֛Ac`=彰@HLL?,@ B4A/>RG=ְֵ6>`@B-AM>.H> ]?GDu@ -lxBdA}a>བྷ+*=X@ՈBZAh> -=(r=.t>\*@$ҀB(A>{Dx>@(=BLAx=L=Y=>?sP@Di`B:FA,=?> /un=l@TܐBgA>J=:f?x35@ZpǠB^ AS>;u>Y=Ky=[Ȁ@^YBe]A>)]Ō=>n8@j BAA>vݾHL>5 @7BcAGO>%y`<ˁ|>N9D @8B)CAH=&˽ڽ;c??G~ -@BGA=|= <ؖ> @F}BLAf>(=q*p @HBA>)=>I>4@HBdA>lX>=~?:@H<BMA@>I@<(<3{>ư@@U0ByA> =O=!>$@;B3Ay>4?r@̑qBztA>2^V=N>6@BhA$=NV9==ӤP@CBA>\>D-.?@>B@FA)> j=?]d@ 9BA9=IV=M?~Ih@|BA=`2>%=FЀ@0tBFÁ=K ;lC=NЀ@IB\A>6PE%p=A?yw_@JB"A~{>I~*=T lK?,j@S\0B>iAb>hr2 -:2>@`]DBeAb> q7=Ic?]B@|޶(BRA=K=<R?A5 @xBcAx>=HCq? @ K`BJ4Ae@=H=!fd>ۺ!@:BKA=7=A<?N5"@q8BFA>'M[;F=#@dpBvA >r=UR?Ax$@ BA[>r=h=n?V)%@%j^ BA>(8=jC=>&@2B:A>)P`#`}g? -蒀'@5>/@?b\Ҁ(@B׶BA> <=ts=ݽ!>h+@n\0BAԔ>9&Mp>(%,@vBsA=]Z0=>-@w?@BWA=K4B?.@|gBuAK= <+=?/@}JGXBKAO>AVBA;=Z铻m>$1@XBaAJ>9%=h? ـ2@BLA*><3<>K3@PBR9&<ޙ<肊>6 4@t"pBaAI>>=hY?,5@B9A>˼M?pB6@BFA=ڼ=r?lR?G=7@,PB6AO='1ZǼo?38@%B.A>.cz=q>]>9@r(B;LA-i>w:@FBVA>=n==Ѐ;@P8BcA=ᢽv<<@.BhBA>5̽Pǽ==@:.pBm[AA>P]=? C>@FFwBfA=p =>=kU?5?@KZBjAXP> -X=,={>z@@MB.A>̲?i A@S0BR[AU>8=DDB@S B -A>&wս5B=pC@^\B-A>#~载(a?HkD@dBzBA>G -I=g=OfE@dUBA-A>inh$r?\F@f6 BAY>:w=0 a?0#G@iBPAS>GnǸH@n8B#A>u=e>èI@ojB"A=a=#j'?SJ@p2:BkK@t\BeQAI=`<ՙP=M>L@.7hBvA4=ש<^ɽC?cM@^YBQAi=њ\=k"<ި?DaN@VEBX9Ad=̓% = P> O@zBwA:>: 1=">P@؀EPB: A=]ƙ|3?H~DQ@ۘrBEAo=Vu ->]R@bhBgTA>Z#=a?% -7S@pBeA>]=|>YT@DB@Au>/'>0og>ʀU@5B'Aӝ=̘ES>V@BAL>F~u=d?|W@^ߘBTJA]>0==w= X@3HBA>½}փo?菀Y@B23A>S}q?oZ@2>BeA==mb=v>*[@4F8B4=Aa=c`1IԈd?*hF\@>Y(BA>>>!v H?]@?РBvA`>7ʌb"=>WH^@^+BYA6=T2 =*%>_@r%XBq>DԼ6v>Ā`@ŠBQArt=x<טg?1a@ BA6=0=G>Ҁb@<BDxA0=5RU;b?< <̀B>#e@ZB[SA]==|O=7=0f@,`BsA>w=?x€g@IIBwRAŔ>- @K=÷u?eUbh@ZB&A^=*=ŽN n?ci@ܳBvAX>%%=C>jj@˒8B,vA>;Rx=a#;fk@CBA~= ?gl@@B$A>ځ=ΠybN)>*m@ޒXB=A6> wƀn@dB=PA)>.<%#>'΀o@MBmAB>=D?l!sp@/B5{AU> -= "^>5Ԁq@BA=6k|(y?+rr@#\BA== o=s@S0BKA==|:>-? (t@V,BDA>Hʼn;{O>Yu@_L0BA5C=̀v@e -@BIAd5=>'F<@k=Ew@jExB2A> -x+?x@sHB_EAٓ=X==T>6y@BAs=ْ=?0?&?zz@B5AP> ByüS<1V{@pBj -A=ȵ=A=Pq?@|@Z8BA> =¬l#>}@6BhAn=U<>- =?.~@2BSA==W<*?N@"`B3A>eu=S6=렀@B.ZAT=9 ߙ?c~s@8B_iA=Ȍ=IU=U>c⬀@ӗBAc1=3y>5> p@} B@A>53|#5pE?[׀@B_Aە>+`==W5>Y@4" B@Am@>E?r=1[?'o@QBA>TT> <@@x`B\AX>3pٽI'==XP@B2ZAB ><>b>@@BV&A,w>@0u@+z BAt>=?N4?{@.{ B^=r=Mk? =@1 ܈BfAZ>,2O= =,H>K@<BNzA"=Ҝyž<>@MYhXBgAJ=>@QhBC@A>"C:>z^@S=PBA4>9~^-@T`Bc/Ai=\,>*=qr?d4@VB6A>1:f`>6@[0B& A=T;~S?m@~dhBKA>3(=ZE<=V@;BQrA5>2UF"<"e>1@BQlAQ> >FM<4?]=π@ XBA2=ѷ˘>@7f0BiAb=W==H?"@2B A>6-C=zq>;@SBA >S(>! @ BGFN1=l(@EBGSA_i=~=R< Y>$@+B LA=/=e?IÙ@쩐BCTAfZ>=##:?d@@\BfAW>HF=on=?0\@G>HXBtAC=׽*l>#?UW@OB%A >+s:uU>O4@V BA>{'8>IQ<@VlBAa=gѽ ?.}0@X<@BA3>,y =P@d@B?A9=SI<<@E pB?A=?=a8V>1T@Bu-A=6Z=Q?n)"@;eB4|A>E벽cۼz?k@wB{A-=@BJsAhD>2$>)n:@B24Am=>1:=1]|? -A@ՊBlA>> <=?,@jBJA>2ˁ}@hBdA%>w< =<>,@ߒ`BA > =̽>@B@AU=SKDz>OL@ψBA=ϲ=WUU>3@(Bd AM>=w=x)? x@B>xB|A=|=6S@MB]Ac+>CC=h=H)?;@#=,BeAR>3UNB=c?P:@/fBGA]=Ql=?\J@>hBA>:9B* ?88@DjJB:A(=>Wk?(A@t=UB+qAl>%nXZ+\>@~ BJA=ۈ=Ҁ@xBFA>"5W;K>䬀@BQA>+9V9=M=Gp@fBJQA>$:UDR"\'dԾ>\@BJAA>+-<$^<_3>@K B6XA¦>/rY#aJ>ӽP@+8B~`A> :>=?\.@r8B!^A(>Fx뜞va>/Ԁ@] BG4A(>"J=T< "?c$@lBjAN>̽G=V;@B;AO=6̯\L&>(@j`BA-=$=y>¨ր@;<BIAA>c P?Hŀ@-`BKA> =h@w:T1=w;k@/àBgA= =e>h>d>LZ@4B A0=$BGrAA>e=<>܀@EBKBl)Ag=~>54=?Te݀@FtBkA>'}*=v=> @IphB`AH> Y=_?\G@STBRA$><T$<շ>OL@Uk)(B:Aô>>!sc>䈀@Y=\BAG=>Y>5?|L@ZYB5A=4>c켻> @]BIIY?`>@]EBA> =wwx?@^B{#A~=V'q=;>_@b;ԸBAQX=񚚽N?}-@fzUlBJA>,м3<_7>@h NlBA#>;K@[?xVq@j(B>pAR> +/>@l?BqA===S>gҀ@t\DB@Ai>=)C?e:? @yR$BSA>6=Z= :>C@z B>A1>!:>S?g?z<@{KDBj6A=n=0=?$0@#B5Ay>Jf=JG>n@`pBA6=e>KN.>RV|@BAa>Fٺ?G@YB%AU=+H>=`@(BJVA@>1nPNesq?f@TB!A>J2Er P?.@XBA =j{? J@vHBM(A>&3.<4>Y@:lBnAI>+뼔w)>f@ CBtLA~>3$=}6?4_@iSBUAx== 1?5l@H,BXA6>:?=W|=$>Ҁ@mPBA>2R=bi=@ÅBFA?>9tHR@ bBVzA/&>?Ј9@ɓ)DB0+Ad&=צ?R~X?hm@啔B\A~f>8RK}=CV>nvЀ@|B:AO>8/v<߼e]?tIԀ@ܠaBMkA2= "A<c>&@(BHtAg>q<.p?">@BjAw> =r>&@_BAc=W">I?@q@HsBEA*>.q =ݮ;@> 2@mBcA=k[=v=@%LBCVA>!v=p:0?U@hBhA >U_$=J>]^d@ q}B5A9>0}=>_?U@ x:BPuA>?V;"<vh@. BAW=+=-S=:t@*DBKA="}8!==4e =ݱ>> -@/`XB8UA>.м 5= @2B.-A>'ȽPN @9hFB8A=_?* @T3BZA~=e=TQ=.m<@V$BA >>]˻pp?'x@XBjgAk8>+b9J=׶?oЀ@[,BvA>7=ȉo=I>y@biXBAg=!k(>6?OQ>@gzBaA> e9=fF>G@l -BA+>(8+=@mV`B_qA=꟔P=V; ?`@ppBAAr==<'w?nؙ@sNBPAU>76=<>@uB#7A8=8Ty}hX9N@@u|DBQA{=$Y<8 ?@yvB)A7>>>A=J?)Uc@& BoA1`=}I@HdB:A)=UQ[{E=h@BrAb>@6==)=>' @fB@ A>F4:7:?b!@B UA>.:ɽJ?l"@oB%!AHi=q<ڽZ>?u #@B3AA>=u%?i$@E#BKAAH>H?R˛<4>U,%@`B\WA q>ӹ==? m?M&@iHB@Aq>- `RF\?-p'@,,BA>.Wq>\(@cpBEAӣ>zݾ k;c< )@BAy> $=I"> -><j?6+@B/rAgq>=% W>`D,@T B/Aq>9w;^ ?--@^BQA>IS$4L<:'dU? v/@DFBNAI=z=&=zv;E?sـ1@L`B33A4>';l<üc?gπ2@BhA -=< -M=$ -?}u3@ |B<^AK==V7?8+4@ j,B7~At> A)?E5@"#BA>}=D?Aef6@3"8B*A7>>6A > 4m=o"q?F8@6B Ar;>83G=;j?{9@A B<A>vd?G?]^:@F;̬BMA>5<6>ʀ;@FmB|Aw>:Y=}?ۀ<@Gé8BoAF>N=}? =@YpBnEA{>f_q=.>>@[BUAw>;== ">?Ro2?@]^hBFAT5=<"e?C@@^QdB*Ac=LȽ]Ͻ>A@^rm,BAA4x>->ک?VXB@dyB*A=\Qn0??C@i3BAm==1ML>D@sQBlAo>b~X=M?j\ME@wg BZAM8>I0<=2_\?WF@xUBVAj>q=>RހG@~oBt@A=2TrM=p=+`H@?.dBA=@=ˏ*=??mI@B9A[>4=9 ~~>N\J@B5uAEi>6<þlZL>ܜK@BoAn>A7d?M$L@4BA>*:?œM@BA>>D=Ϛnz?2N@vLBRrAe -=ܯo<?XFnO@ B{A=UH=ќ?6P@\B[A>!=9=9?T{Q@xBXA#>/ F;=%g?rrR@B{AE=D_=a>ɺS@lB[AE=90>1T@- B&)AK>?ܼIXMR?(U@eBrA.==.=z?a V@BtAܤ>2I O= ?]SW@8BA=3$=zl>6X@CBA>(]QP1>+B?FRȀY@B A->p5='`>@x[@ BgAA=1z=?N\@B'A-=7F=D?|L]@b BqA=P=|=TA?`5^@HBAV"=ս=>!b؀_@BDA=吁=` ; -?ր`@*dBA=–=҉[>>a@+BOA=㾁e<=LPb@8BjAx>溽=i<`c@:BPAf> =.<ѭ?jzd@FPXBkA=s=Cr=\? e@K:WHBlAs>%==%~?vjf@KB_TA>(=S=UA&>슀g@R0BA>듾H?vqh@RPB Am>3ƽ5y? i@]6$BfjA>W=)=ڍ?oj@aEBA>=";k@kB_Au> 绽;=V==z l@pRBOA[=rX<Ã??ƿm@sBLAc=ޖܼwG>?0n@uV8B34AB>WV<ɼ?[o@uBPA> 4===E>zq@_pB[A>#ؼL=:9?jlr@eDBA=[YN4?;es@0B A=$-<-t@{BKAp>,<|w?ngu@ļBqOA:N>>f=/?6 v@^IB;?Ao"==5NԼX5=$Ww@_@BFA@={;5?..x@>4BGLAC>5ٽ|og< =y@{B#cAR=bFkgPr?X]z@-BIAhA>A='D2 H/>c|@ ZBKA7=y=<5?Bo}@sB*A==vɽ/.>uzp~@hBsAQ>ev==LL?C@RBQRA= -;q B @ݵ8BQA>=jn6@#B8Ai>9!=*/=v@+ BDAAb=>s;+C?tN@aBAA4>ۨ>aL? - -@ }`tBA > -*g>L@vB4yA{>+˼b?~@BKA_>B}!=me]>߷@$BFA2=e'"<h>@!ґxBT7A > EY=I?@# ?BmA=>9|@.̆BWA=ՃEf?Rao@57֝Ȁ@:|BoA(>A@IBTlA\=VU=>]̀@JBL]A6>$:!B:X?]Q@Y]hB+AAe>'n>.@Z0B\ A==f=<4? l@Z'0B$A>b:r?k߀@\3 B̾&G?D@_.DBqAIw>Fi<(8==@a -]0BbjA>)}R=lKv?4f@mk8B( -A> > K1KDj=7`@v0BAk> o=/=&@vBhCAM=p˽ƭ=> -Ȁ@8BHA > A<>u?@ B-A> g?B݀@UB1A?=O#?y@gdB(aA*>4<3B,> @tDB A>>_R^>^@B: A/=f=|Hu?@sBeA> н=Xh?_|@{BblA%=2< M=lP??D! @싴B`A/ >FzBOPˁ;1?09 -@wLB]At>u>~=It>o<@@BGA=fJesP<<@РBRA.z={@1ՌB5A=Uм]VY>|@XBA=A½k=m>&',@^BA(=G -?'@4BAm=KV=2齗r>@PθB~A= >1&?N@ -B A4>I vl>.@/OB2Ar=H9@=Ѱ؀@$Bb:AP>5Gy'F=j?#N@ B+}AZ>%a= +>@&oB4aAݷ>4Wg -Ck? -D|@-BA>#@)A>5@@34BAh>IO=(pUc? 7@?2bLB`!A2> dr =[\>cml@@?BUA>(%2ws=?H@JXB8A|=^I=PO>@OOPB1Ag>=:=%6rF?F:@QBDA=cREw='?Z<@SqlB9A>=>@XE1BH"A =W=<%2s?Sŀ@XuB1zA>eu>.L@X߼BCA> i= :]vQ?aȀ@\hB?A>#*w?U;@iBbA=ވ9=su=of=Ϊ@n Be AO=c=i={?E8@owLBAh=Q,?e@yh0B4A"Q>;g?%>CB@zהBSpA|>,=<)>i@zϠB;A`d>m_?Tˀ@{j@B[A=4 -9*=9.?8@b$$BqAC=k;=9?oi@ BV-A**=M<-=y?@qkTB3A̽>B&xd-?ra@SpBpA̒=NF=D3?TQ@BA%> Q=׺wT?U@78B4Au>;<>$P@GB)fA>y=PN:@?qCd@B]AHS>Zo>Ā@^0B6AU==a({>@BA~=c=1?8,@۔BX:A>?^== >n @BVEA>8b=f =?_@0B0An=̀L>@ÑB~3A8>4a<=%=Ѐ@|@B2VA?8>L>_^d><@˝lBA>%z=^y]>v@xBXA= ?}{B=%5?X(@Վ BQA>=n <:>@|BAr">H=;v֜>@rB&AX> =kA}(?y@BU6A>|[>F= -"?)$@B Ag>! -=O?-{@$BbA}=H(=o>,@"8B{A>"$=fuk=?$|p@{|BA Y=>o܀@ BQ Ap=ž,z=d<B>@& HB"A7=$[=—kq?΀@/VR B{kA>&[|@:8PBTA>Ϛ@>\B%Ab=KԽU>'~@B_XBA>Y:{c =@C@BN-Ab&=E=W<;?T4@@C&SpB*A=Ӯ>9{2;d?@ˀ@G;B/hA>)bgG0}=`?o@OL BZA/={=5$s?$$À@OذBA E>aI>?2+i@PQ4BfpAH=gO=->@Z10BEA=ji>*=ɷp@\}>B`Ag>3'-='=_S>Eǹ@_BA=>$qHq>t@k!BJA?=Zֽv@qLBAe=s>)>F@yb"(BAT=ٱk/8?!v@y̽BAT~>Ś+>Jx@LB=Ab>(Z1}l>@n@B A>3> -_[>CX@1B8Aq>= o>8!=q<@ȍB1$AS>ڥ<(U>4@BXAER=&S=#L?|jw@dB"A><5=VӴp=᐀@ɩB9AY> Hl4 <:@pB9A>)=>~>`@VXBcAt>/j@ -=wyL>@ϭ@Bj A>þ;Jt=bW>8@խ LBuA&L=^>&=?@վBA5=VjG> 8@BA>=O =@1LBF6A9=>c;im>,@㻍XBTAP=u6>=TK>>d@\B)A6>ݩy9&>΀@ -4B2A=9Ӽ>=#80@ BeAl=>1>+=ф<` @&YBA>@s>s?+>h -@ jHBp+A>*-)=6>nj @ IBBKAj|>%=,5<,>+ @v84B -NA>'Ž'_Э7?0 @a BA>.&2ӽ)>@bB0A -=}>C6>~B@a<B$$A=-~=yZan,=EX@RjBiAo>==8@PB-A<>d)>@/BAk>e=Yý`>@3BRA9>uڡ6= R?B @3MB7A>=><#@@9 B_A9J=ͭs<|=Pq>@={B$aADs=>_o>@AxDBF?Aw>GE==;ژ>L/x@OBPA1=ඥ>E'fH=B+y? -@XbBQzA>Y)g D!>z;-?5@_BA>,=Đb-l>bP@f^xB&A>HȽo}O?x+@hBB_AD>$<`>h|@pRxBjA)> Pkxk?B0M @x*v4BoAB=")t=?;d7!@}B/A=ycyqA =-?"7#@u0B3A>m H ?냀$@XhBpAW>>O>=`+?%@PB/$A?>'<= &@9BTAA> -h5=Dc?0'@!ElBAb=̀>)?p(@BX>Al>$=~= 1?[5Ѐ)@~xBXAo>='=$v]?8ڱ*@‡4BH}A=ͼ9 +@~B%A1>9Ƭ,@3BPA=7%~'<Ѵk>Ԁ-@Ñ0BA*=l=s_>;\?I0.@В8B'A>+'>>=쐀/@ԚFBnA=P=>ew0@0,BdA+=ϙ =F=s>O1@8BA{I==M >O?l 2@<B ^AƐ>>=)b? -o3@BHA>;ɼ39`=84@dWBX5A`K=E-.= ?m{5@ޣ<BRAg> +<|?U6@ B9A(>)@&r==Fl>T9@ \$B$Ao>C=6l^gJ>z:@LpBNA=`0=ܓ?;@%BVA>)M/K=q<@ dB AK><ݻ8?%=@`Bp=Aq==9@=Y? >@.`BA'>$<߽r>j?@/;BA==K?q݀@@;0<BC,AmY= !:r&=5xA@>BpA@>I.G%?89C@C0܈BJKA>hp;N.06D@DKBA>+ؼ7>"E@aP tB#A>=jIi^>Ö4F@btBeA*> hZ<3U=> -,G@rBFcAb=@n80;>ЀH@B5AW>K^Q>HI@"BbAS==iΓ?1J@.|BWA==f=?FVdK@BAtl>RA?0F.L@ѵBIA=Dqgn_tM@ B9A͹>=TٽQ? ɗN@TBEWAUK>B =;N>LO@JBA=d=]>\>NP@}|Be3A> 5K=~u>q%Q@XB*pAC>=2ϵ>%R@ *BB#A_>hdoH=~\S@(BJA=ܣ@=vi%Y)n=d?QU@BZA]>@ֽR5=4>k@V@ľO Bi,A=+'A$==}?F{W@Ŭ$BHA>*O;,4?{5X@X`BqA>C<k?N`wY@ BA=AJ>¾Z@טBqA[> (=|=?P;[@BAg=%x> k?8柀\@BAd>a>, 6?X]@?B,A[=X>#G!C=Ѐ^@_BA/o=K,_@ \B:A\== ɼjR>w`@XBO$Ar> 51=<=@a@B~Abp>43=-=e?v$b@CXBxAu>o<=BLtsB>d=P<+ d@BWA\P>=J,< e@%BAA9=b6/ ?$@f@* BA=Fl>?2 g@,BA>H3$=W?y&h@.@BAAO>`=úG?v9i@9װBIA7=(:2\Hj@CXB|A>KtU!= U?*k@HB AM>k=co"?4l@Q%HBcZAD=j ag=sE0?| -6m@]|B3PA>pNHC>6n@bikdBLAoy>6D޽<{?!;o@hxBA>9~=f?Wp@mBHAW>,dG<7h? -ꒀq@v|BO|AՖ=u>`~|=l?g rs@dBe'Aw*=K=ElB==%>5Lt@BEA.==wf$;>Oju@ՔB sANl>&< ;ȡ}?oVv@BLAA>3=Á=fh=|>ޥ^x@BA=U=>#?,d<y@BZJAR>#~<=/?W@z@$NB:A>!U=h}< {@#$B]A>'-9=J/?"|@^yBdAɨ>a_=|fB>}@BHA(><":? =~@ (BNiA =&=<'tJܽ>@%ܐB8VA%>*"< .t;?D@!B&A_>l8zS=0ֽ>n5>`@ò<BovA=W_=>jZ@ǛB,BFA1==c?P_@Ǧ BHA)=~=*> -x?)YZ@`BA-=g=Ѐ@BXA-8> +>l@GjBBAq#>W:+S?!@ܞpBsA=>#'pܼ=>@IS|BAe>5D=u7>H@֤BfA>D(ۼS=>Զ@BA0>#p@xBejACJ=7#e@\BVA= j=qG?!@ HB,AT=O) ͽ"`>@n<BAF=)]fC?r߀@ $B|A4>#"==>?L@)SBYA0>)>*=,=߸@0NBA&=t}{>f@GByARs>8@<+*=4< o@HBEA=5=m;5>Z@J}B9Ae|>>(?+@J?BA:=ׇ}/p?#z@LGB8A>>@Sv?l@Onh|BtlAǢ>-=6C=ʝ?N߀@_BObA==<I?B@ekBmA1> J>W=>uX@vǸB?A=>%僚>e`@yBDA~h>GI=S|;+B>%@BIAq>E0=:|BqGA,>8$"c=L>ߧ@ph`BHA>>u =V J@v@BAI>=Zn0=}>@v@HyBspAy=QM= =)Z>@B%Ad=3=TI;YL@B7AF@>I4<5my/?Dv@+-|BRA>v="罣?j߀@õpB!AC>!E<"q{>c@SB8SAӨ>R<Ӽt?ITƀ@[BBqA==L]c>@TBMAH>+1l@=5!<%?FBǽY<̷?\@e0B8AO>Aw\PVPZ= 0@`BUmȡ*4|?w@;BVA3>)(n=(?Y @jB>A="=we?|~@5(BpA=N:=^?g|@ByA>:==˵?u눀@nXByAYH>(t= Z=˵?\Ā@ uw`BLGA>IQ<.?Q>|@[BKA/>6c=z> H@ pB=A`>Ӽ?4@%BހA==`Ǿ9o?/@0,BNAX>f<>s<@:;4BA=J~2齐0 ?A@B{ B,tAs=FS, $>:@MB9A>;X\?;*@SB1=}<@VC9B|A ->;4n=ج>4@WNdBA> )>>G@WBTAJ"=<֥%=>Ȁ@XzBF -A > =5;3V?#U@^8B:|A8>fd[(@e1~B+"A_>(~d=-\?:U@keBnA?=?X=n?0@qBg7A=4r= ?H^@yBthAd>4+%=-'=l=4@ze BLfAyX==IS<>I`@wlB5#A>!D˼T=(@k(Bb5A*>޽K'=j?\R@B8A=44X缛>x@߲4B%A->+8=f=C=@B[A=*ĸF>{>@Z^\B5tA"==ZK>5n@pB=WA=>"{|{ż?@BB0A.>I2{Q!=ׄX>|*@!B"A+S>#!Jk?6@ǥBNAx> ј=<= /@@@B(|Ae=i}2A>@B8A>/=R２?A@!4B$A>$[?3À@ӌTBA>\k>z@BVgAfd>9ӣ=#=>Ѐ@hBA#> a=Y>>&>ىҀ@LBA">Aj=>%?P @qBfnA-l> -I$@B'AI>A =YFE?yՀ@;B4A)o>Hk<'`=@-BWVA1*=x>X=>[ @BuA>Ͻ=,>]@eBMoA">){U@?Ř>Ԁ@K|BHAH=S:f<5=M@hB?A=׷4>ڜ@ -B+A2>Rހ@8BA -S>2e+K?%Ѐ@8B@A=:;I:>>?I@BbA)>μ/=p? -?F BA}>d?g #? BFAn>vJ=s;>D@B(A>7R>\>)@6PB?KA>@» >@B? AX>׀=ݘ,<<@ۀBEAY>= -Hb;$>~@&-BQA>)3<.>܀@'B-Ah>X,=(=d< @(RBf1ANb>}<缌=>C@,<BEA@>O(u;?N@@.BcA>x;0-=u>@0BnA>Z=D[>^@1BGeA>WҾ:<r?+'@4X*5B$dA>qAѽ_?K̀@9kBCfAf<[:=X@:(BI5Ai#>ۤ=)Ol0=sg]?N@;PBodA>3V=p=0;?A7@>H%+B.{A>ʺ=@??Tʀ@@HBWA`>R_1H=>|@@ԀBbAr>>*½6`=n2?r@ABAeAov>z=5 -[?Հ@B+BDA!>m̽;~?<+.@B]B?AE>w2S=a-?i@B.ĀBRCAe>*n<?o@CB@A>ټZ= P@CB:Ak>[=- l0Q><@E0BA3>>b+T @GB@`A>=K> -@G.BA>]8z> ->:Ӡ @I<BAO>1a<?zr @IxBzAT>%=RQ?5 @K -ߡBl.AA)>8<=Q)??@L vLB\AM>Wy<:=B?g@M=OBA>ʼsZ>>n@MMB@AXL>tT?.h?I@NNJBBAY>k= -9=?]w@NaBrA\S>ʠ=N=?kׁ@Qb@B4A-3>, @RBKA4>{=ùhZI@R`IBPA>iPW=vE>Hj@R׿@BXbA9>UPX@S2#B_'Au>Z=B}o=Sy?OO@Uki:8<>l@WB|BZA(2>_<ן=4?=@WxB>A`>a> :>D@WWBQNA >eBW=<ڛ>&@XXBtA><=0=0?z@X@B4AZ>ߦ2?@YB.SA0>֫?wQ׀@Y+@BQA >)мs<0!>Ȁ @YO(@BSAf> >]<b?Zl !@Y]GBLA#X>H;ޏw =z?">ڀ#@Z~B;Aj>Q2=fC=$$@[BB/A>~Ehat3Uf?۠%@\BTA9>\=P>R"&@]JBAf>_=$> >>ƀ'@]s~B83A -V> <:4>:j(@^C BA>YLg?MbЀ)@^׹YB;A4>o[=+ސҽz(?O,@`^BJA6F>\ES^.@a'lB`sAl>ͻw=]?@/@aBB8AO>wx<좞>}0@ay@BqA>g_<|=U>z1@bB'Az>;H>02@b՜BA> ->ty>E23@cW獠B4A >z'Mֳ%>耀4@cs诀BkA.>$%5@cȌaBSAn>0=4V== ?(7@dCBlaA>l>OB-=3Q?<Ӆ8@dBjA>ZO =x_?bct9@d̠BQA>ht]r8BZ>;;@e} BG(A>6K?c<f{><@e%B'A>_]a:}D>>?=@f[JBmAj>2"=,{˽?[5Q>@fBLA >Ҋ><3=?@f0B<"AT>>? ހ@@g(B>rA\>=ML>ĀA@g} `B]A>Z:w=H>B@gB&UA>=l{QG(?@ZC@hX@BA>Q$(9?RD@hB$Acy>c=+=mb`?C݀E@hgB0A>S><F@hrIB^A'>)x<p=P/>G@h1Bm A>Zj=>RH@hoB%=A>UYg?isxI@iaBWA\>Q=V?q*J@ioB.HAW'>p=\`?kK@iBK,A>j5<?hL@i2BCA>J=eA:?eM@j( B0Ab>= - &?%ӍN@jQ`BcFAmw>W=2=r=ȀO@k5rB:AC_>aT=?1GQ@kTٳBYA>P10>=,?R@kX.`BAy>`=K =Z??sS@kBB-gA*r>UT@k -BKSA)>> -4l=<?xV@k`BgoA>yz. n0=W>RW@lvBTA>xۍ=f=?up>X@lzژB:SA>;=tW?dY@l6BIA/>X}Z@lBbEA+>+~1=kL>Ā[@l B6A>b,漰>?7a\@lBLA>]|Z=0Gp]@l1BOAO>OT2i=w<@_@meB70Ac>^t弨k?.z`@mS&#@B-'Aj>yE=,q?:a@mWf`B@AA>L=1ד?eb@mB:5ADe>33c@m[EBA>BV+>nd@mB^A > XdQz=e@m̻BOqAN>i<藹< %>Z`f@n[B{FA>oJM=cm?Ig@n`B)AD>V$ʊ81?|$h@oG@B<-AO|>3U==瀀i@oJ BqAA>nFҟ=?Uj@olBynAqz>N];=:"=i> k@oBRAb>=4< ->l@o$B.A:>yrG>m@p -B%A6>hSZKn@p0BM*A>NWtN9 r=:D> p@pKjBJA>OMTP%⽇>Zdr@pqtB$A>D½a=`Hs@pd`BDUAW>}Cݽ0;R4>t@p2X@B8A)>mC=ϣ=؀u@pfB!A>=js߀?H@v@pАBSA>z=bM<J=ݲw@pGPBAA>':dv?m8x@q$_B3A|>m[= >ːfy@qGB3DAp|>xZ> hz@qP(BA!>kU4;i>$?0`{@qe+ޠBrAr>p֒=='=|@qB$A![>\;X\l=Dx}@qBz|A}>lb==G?I~@qbBLAk>Z`u<=Ȁ@r B9A6>yq@<<> @r"B-Aa>m=zν?F>@r7"pB\{A>|==@%>M.@r?BUA{s>b=M=1?VD@rMsB%A3>,Ｑ T?k @rYzrB\AO>WWPq==?#@r@B=A>MҴ]> -B>Vh@rB%AQ> <}UM> -@rE@B"A'>V;vB7q`A>qZ@rɨ@BhAm>#cZ -=?zY@rqBWAW> =( -=Z7>At@r߽BBA>!>/<]9ǔ?u\@rPBgA>E=>5@s:BHA>9p<7)?Ca@sv. B#kA >Jg>Y@sB-A_>&єT>|@s@BSAA>= 0u=݁=`V=Ѐ@s>*`BF"A%>lό<Z;?L@tv.0B:A7>u<2o=>yt@tB/A2>p$=5 #=j `@t0} -B&A>>j<.ĽO?_+!@t4k BPA>\k=<Ѣ)>@tQBmUAݚ>Q&X =F>@tc$B A>=+(>q4@tlPBRgA(>\<"e=e= D?"@t$a0B94A>19i<>g@tèB8AQ>jll=xH? +À@tBB/LA>਽>I%@u7BRAv>w=M?Հ@u  BKrA̬>y5i@u. _B4@A>Tl=lܘ,k>ƀ@ud;BNAR> <@?qL̀@up`BbAk>QJ= ,??OU8@u B: AV>8=%z$?"@uB5^Am>N2õz?.@uMB*A>Sʼҽ1<@uapBOA$d>+;@uRBA>cn?pd@uByA>=1>@v 'BVA7><$>=?cɀ@v2fBf@A%><=Z>-_h@vL: B>A>|9&e{A>[Ā@vOBCiA>w:o?>g@vrBKRA>}=PPVM>N@v B5AY>b@ܼD>~fT@vyFPBTAR>>==>@vXpBI"A>^]C=A@vB(1#)>D@v>pB]A>o>?܀@v|BAN>ihsa?z@w*sBAA> _>Z@waʁBvFA>=d=?Qr̀@wp6BjAX>:=N>@w-B1A>g1>Ԩ@wqBBA>_2=dpg{>n@w縙 B4yA>ȧ=d9Ǽߜ>fĀ@w<7忕?*v@x!B(HA>z="Bɉ?@x!яB'aAI>g<]TI>6@x(Q BoA>O^=u?"ʀ@x+!B=A>>|=[ټ ?Xb@xANBbA>м=n>@x#B$A->:,b %=K@x -B.DA>VȽ˜?Jh@xا)BqA|[>L?f@y`B?Aί>Ny>]8@"?[Ѐ@y x@B+A>н(}u?L@y&,`BMA>\Mk=U~31?6\@y`BOA>=O<ƻ?.π@y@B Aܤ>޼s$y>@yRPBbAR>uj1;V={@yyrBWxAJ>ڱи==@y!bB^1A:~>u< =L?r@zBcA>Ũ=X=v?QX܀@zBgA> X=1=t@z#* BFA >].;>=P@z-B,A>T Iq>"@z;D~>p@zPi BAAuH>h3:<@z^@BAC>Xj%<=f?u@zhWBgAR>bk=#ꓽH->d€@{Gz@BGAU>RpD<5?p>@{u3k@BUA>Q==ª?}`@{ pB`YAe>c(k1=]+I<@{BEAd>= W;?I->@{BBH7Ai>ʽ<'M>V@{޼0BI_A:>T<|QT=g֙>u@|ABcdA>M=U:=s_?)J׀@|lj@B$A >gJh;7=^R?@|ڞfBA>Y e`=ŊH@}-B5A>jhZTC? р@}z)BE A>>=%;9?) - @}BZ2Aa>S:fZ=/fM>M@}WBwAae>,=ѽ{=?T@}yPBJA><]> @}BPeAMk>^I>x<)?5)d@}RBDA>mؼn;1N=(@}BS{A>m6=)CD<ћ@? 3@}Ʉv`BASA>;=xbr>Ҁ@}Ӱ BcA@M>m5YMPBvA >\/S=={ @@}ZB^(A>Ƚһm=L?32@}B0Aɚ>XӪ= Q=@~3 BG|A>=G<?Up@~!ZBGA>_9<r??@~+BOoAÒ>\=<>@~2HB*VA>{~13=[P@~lIBsAl>=0(s=7]>@~JBfzA>84=;>^€@~q] BXfȂ>="><`@~1B/A>`C: ><4@~l@B3Ae>,|)Լ%<@@~bB"A>=xWma=-Ȁ@~ٱ?pBQkA(}>?> -<4p?WY@ z̠BXAyo>;ཌྷB>jf4@AB8A>Y=@R=Ѐ@ -(0BGA>𽁢Ƌ=k<؏;<@$BA@>5<']"=Ѐ@kBXAey>yN="<=%R>.@=@Bz A>R[n.=̓;@h?BDSA>T*;QQ?p@)BRA>`)=B>2;„Y>:>z@MM+pBIA>W]@U@B.A:>b<`C>p@YhBuAp>P=82=U?_k@fQHB8A>L=~=n @kVB!Ao> =UgrC>D -@k8BA/>|=^?n^  @o:jHB'3A15>?( t_6ݹ>=P@!BK A>%L=`Kx<{^?ak @SBQA/>a|6Gda~M?W<@F xB:QA>k%at=-@ {hBVA>6=J>->@ xB5uA>y5is?oN@*ɛB3A[>z噭? -@5sBbMA8o>_<=ka?VRx@Ix BlLA>k=f?.\@_ -PB:XA[>gX>s?@p*`BmnAd>=!}==d=L?e@-7B7A>=7&3>r@ѸB!A?$>k8$? -Υ@ΛqB0A>>"Y>IV@᧺B^FAY>c = -=M>ր@BfArm>a=;=? @ʴBTA,0>T<==*>7o(!@ =B;A>ʼ㒼I?BL"@OzB0AAn>bm>#@+BRA>< `=K ->}%@8cBR8A>xS=\ -=cҼ*|(@aB+bAH>Xo<+y?).)@B@Ar>ŽẻV?x̀*@x@B_A\T> ->#=S>+@ؗBAiAi~>=*cs?+n,@eBsAU>W=X' ?||-@(Bi A:B>>)=d?6$.@BFA>/"E;?-sn/@xBAA/>j׺ ?c)?0@ -:BQhA >Tdk=6K?Dz2@lBx'0?K{3@—@B5fA>=\ '>^bd4@ bBmAG>YA>=v?5@E @B?Aq>s$=:ܻ?=6@LBXAn>oUXN=CF?/+7@MNBSA>԰G=p|?jI8@m4XB8A>pKY;V9@o(B+A >>e=(J?#:@oXB_Ak>~?==Y >*;@rTBA>M/{Q>|>2a܀<@s@B;AU>gM0=!G/?z2=@BXkAA>wvN<؆="3?0|>@QXBKBA89><<m?I?@(øB&0AO>M=SZ+=O=+}? wA@. BBA ->1S"=B@XB9$AX>Y=> -C@us`BVfAy>[|=J =.?\D@mBCA>}Oxo;O>zE@̥ xB{A>[=Bf==?`F@B+Aߤ>P)i>^G@HB`A]>=?=`}>ҀH@BBAW>W=)'?,I@ BVAQ>1=U?&hJ@4{BWA>l =?-ـK@5U\8BSAT> #==>Z7ȀL@@?BPA4>UϴԀM@H(B;A:(>o}<]?IpN@JHB $Az>Z0=L4=yO@jjhpB LA>V|Y?Bz P@LiBSVAű>w=7<=oQ@QB9A>Rv4h7s?6R@݌xBAb>i=?!S@\xB%IAњ>[{=$X?$1T@ǷB8A>Ph;~&L<=$U>ʀV@4}BPA>^=a+h?W@d\xBXoA>0="=">dހX@*B!Ar>%½̻s|?HʀY@]0B"A >S[ ;*jsq?J+{Z@BIA>=T|iнM:XJ[>n\@KxB_0A>fH{=&#=T,(?ʀ]@M,_B[A>`?=O=5Qg=`^@{dWBA ->VKH>릀_@ҭhB'YA7>f=IE>W4`@hBWKA>L!7=ү?a2a@0B`XAu>x=\=\>ZDb@4BZA3&>V mSj?aʀd@ BA:>^C?7_Fe@]yBHA><7:?X8ʀf@kBTA>)=1=U?sfcg@BXA>:w=$+=,bph@(BA=>捺==?hi@\BA G> ;R>7j@ B)0AT>>i{<1?k@pBE2A>I_O;b[>g8l@:`BTDA>R$6=Sf>*m@RjHB:xAG><,&pm;_n@X0B5QA>} }>Ⱦo@bT(B|A4>&x?Tp@~w8BuA|> =y=>.Kq@}BXA?%>F;Pqcꂽ85:Hs@ݺBB A+>\ @!y??9!Jt@[BEfA>h=j8;B??u@VBAu>=B@=aJ>ev@FB!A>9u㣽t?CRw@RBBVmA<>WMA=Z? tx@O(B9AgC>=%?yʀy@iXB3A>=r0廊>̀z@ -tBHDAQ>˛==ssrJ$?)"|@"B8AG>v<=%3>a}@;{BA>a%6>e~@=BvAk>Qz?K>N$@O~B:Ac>;b*/,k%?{"@YBLA6>kQXx<?-7@^(BjA>[c]S=?h@BGA^V>!=<?ESGs=ML?w82@$ Bq&A)>B<+X=P?w㬀@BS}A:>U<<>O<@TB`A@$>0#<=^7? @iBgxAO>‡!<=<@@wπBwAƹ>ga=G>о@% -&BV1A>漽=\=@&pB)4A>ŝ<.=0@*@BDA4>W=;?@;/3B +A>vٽty?kŀ@=㹨BA>ɜ~s3P>(@irBMAY >L潻oR<#>؀@BBj~As>6= t?D@{8B/Aa>>=- -?vv[@1R(BoA>싼У=ß> \@`BJA>:Sr_= *:lUqXF= =>&@DB7A4>i4@?f@5B/A>NLg -R?O@nB jA>=l}?@0BT9A@>T7=FkL?* -@@HBFfA>uϞ>{;Э>@JcB12AO>yqĭ7ތ?,im@f9BGRAe>G=""-< #E?Vj@KB2QAVB>_0=?5À@"B& A>=Rˎ=@͇ B`A>zrL>u=^?w==Λ?@<0B}A$>t<=ڝP>T@$B>Av.>=Cf]S=@ ƨB5IA ->X/@=Ѐ@'B.AjE>=+ڽ?@CĈBV[A>:==@E.Bm0Au ->s4=A3=?6@I8M8BZAb>g= -=1M?"@JB2AŨ>G]*=@S9@B5oAߑ>NмӤ>v@lBFA>T V+;Y>70@o/BcA@>rm7Ӷ=u?b@roHB:Ad>瀽,gc9?=^1@zZHBWAo>=`=>@|(BXA^>n=z=$-?0@9BDAfO>a8O=#;5S?Pa7@@BD3A>=;BLs=@TpB A>=Ez>@z?BtAw>O=Fu=>e\@BvA>KU==S@o"hBjAK>˻$v>} @µXBpA>O+U4#=g>u@BHjA> \ <-<=@@0BJA&X>U<LEv:s(V>sô@ 8BHA=>{Iה>B\@`pBFA~>]=T;_?TL8@ -B6E nL@pBAR>JT5b>@$hB3A_>V45H9? HZ@;OBp-A=>i劼q-=/? @<(B A"k> -;;i?o@RsB=A>=XlI"=0Ȁ@h6B4A>Lϼ\>C@r`B A>,>6̀@ry(BX;A>u/= Ɓ?t @wtBtA>Yr=O?~d@IhBYAi>I@>=-?H@{0B)Az>Ƥ9t?T -@yPBAu>a=BmQ?A @`B AS>t[XƑ;>\:x@HPBYAeF>LO=!=+>׆@E*B:A>6E=(Ci>~<@RRB$[A>w;`=&@@߷HBm>A:>]U=%=*;>>@-rPBAx><1=G*_;9@OBQcA:>^m<<>y@赨BSA>B<;@1t0B\wA>q=C=?~>RN@HfB+ A>4榽.?@KB;A8>M<7\O=Ȁ@L(Bs AG>e <ʜ$=z?`@M.B<+A>fI=?5R@V9pBZ8AM>V<@e=/t>E@hg B*9A~>iW4?N#@hBtAF><\=w>=?@yH8BA>`"۽t?J@I`B2Az>=ʼI?P}@fxBA>hJp=cɾ.D?RIP@@BPA+v>L\,<7?<6>l"@VXB5AS> <#g?s܀@OVBRqAk>W9t<랤>St@5pBD A>f:{v;9?I"@֚vBA>ur=҄s%?s@B\AL>;ԫ=={>̀@f_BAR\>!<߁%? 0@hBzA>NS=?CՀ@ƅB@AM>O,HuE?0(@ DBHA>\8<<~r?)׀@BHA>x}9&<3u?I8@6 xB1A/>z=#?[ۀ@$@XB3:A1>d/L@%X PBGA6b>\<0<4>b@D.[HB A>W=}jp>&@PBHA>N:7Ϛ<1 e?!Ԁ@5BwWAƣ>1=92=Ý+?=r{@ʘBwA>y-f=u?u@;HB8uA\>G =&Tͪ?xɀ@VB.Ay>Q=J]&P>@ZB=A>gg"?"d@&AhB23Ap>2=7J? ǀ@#EBjAq>n =V4?,@ypB=Anb>W=3͌?lH@ BIIA>A9*Q,@[bBHAk> <$?@BLA<>fu\@15BA܍>T=l< `@79GBGA>=y <>@NB8%A>yOT\q?{%@\tB-A>}KlGƗ?^@_BNA6~>/<<>l@tBF^A?>\%@{AhB[A{>f=9z=f0@lB-A8>yWR}<&1@@B4FAk6>L4Dr=K@ BsAi>csHW=_6=Rf @b8Be/A|>(/==S=>I @e8BAb>sڝ> =x@wuPBR+A}>]Rüw=UP>\@ nByAg>O8=&=J?i؀@&BY`AY>"ud=)M?9V@B6A>_U>Iƀ @80BFA_>=;?}R -@"B2A#4>m;뮠`?hԯ @3YB4VA6>>B߼>{& @DpB>AX>`L=Lj?; @NuB+A>D(=\C'>UV@vB>tA>&>'ݻ^*>3{@w^BHA>cka4_}<#Z?^ݱ@{0PBh?A >,ɼK=K>.@|bB"Aқ>ckkQ>i@¸BA>k٭(s%>d@7xB7cA>̾zj0>Ԁ@NƐB|A{>kw=N7=i >Ȁ@އhB:A^>uz:a?Kd@~HB\(Aw>c=F+==Ŷ?7+@rIBzmA1;>yؾ@ƟxBAA>=7’>w @pxB4vA:>X)==m<ߡ@`BRACB>߀Q@[BP?AX/>;=i<%?z窀@OB@BSEAQ>Ӆ<5<\?d"@YBQA*V>f =V?T ̀@(B(AN>TK=װDE>$ @zBAܛ>K3ޤ/^S>o!@ JBO2Ao>+=^j߰=(8N>q#@=$B/A>O=NȽ=a0$@KBA>Hmu+F?}%@Qt@Bt A(> 1ѽL3>'@BA>`GMx?ly(@ &XB]>A>w=t==Eq?1m)@!B1A/>$*@#CB+AQ>&ѽ(`?+@'B}ZA>Ui$H>=?7,@/BAA?}>~񢽤o>-@J BV;A\>r]=e=܈? ;.@LLB A>m82_>!/@L B/A1>@T;b?f1@euB*hA?O>yrn<^53 ?_2@r|.BAt>HfŽ1:?~d3@y -(BBAK>{<9>μ4@~BEuA[>|&@@j; -?t5@"3TBPyA&>ƣ=\<jpKo~=(?Umi8@. ,B"A|>)9l>R9@|B4hAF>!Π >–:@\B[A>z=^=5d>񠆀;@BGA>m<O?J<@OB2A]y>=%괃?S>)=@ByAa>X=c=T?U>@dB0A4,><~XD=Z?,?@7sBFA>gbh<Mu>@@9\,B -A+>T=vη>A@ۺBXA>en9i-=<>JB@ -xBtAW>[==؉?O_GC@,B4Ac>F½P.3L?TD@<B]A.>CҀE@pBLCA>a <>CyF@h&pBsAo>==y>NaxG@ -sxB0.A(!>c]<@ɽN>H@{B%|A8>MBY!.ؼ”>J@ S֨B8A>9~V?{܀K@%Bg=3?g,L@+ސBqXAݜ>^ļF6=1?mV"M@/1B2A>~ -;eJy> 0N@74BtAo1>u=5=?#O@7֘BbAb>l%=|?x7P@@#BLAB>~<Q@DیB9A>Cv b?;9R@ESBZA3>`L$彆O>R\S@iy(BAU>)=]?oЀT@mHBS%A#>;;<==!U@pf3BN A>>x=]>}6ge~>㚀W@u ƀBEOA1>k+=; =8X@DB^)Ar>[b=<5=Lw?MoY@LB-KA>TTz|>P9Z@pBYA>F>b >e=-_?"@[@VBQA4>\<_<|?4\@UBzA>=߭=Vd>0j]@BFAU>0=<? ^@ BMA|>=tH<^>^_@3BuAA\>uS< =$>x`@9LB,qAF>4a@ B;gA>8׼mTX?:+_b@@B:Ai>=*b˼`Ý?Q'pc@9uB64A>Yӣ;aO?jd@@BGA>P.R=L2=%a?vOf@*B$A>=~<[]=hg@ݨB0AqO>b=9tT?ɀh@NBkA?w>d<]=>TxЀi@XB%HA >o,S tY-">j@ -tB>AJ>^0>!k@HBgAA>"ŽVbn=w?Al@5lB.WA>Oh=U? gm@7BRA>n;d<?6ߧn@_tBCA>{S췽ju?mڀo@ B"AH>ZYfEmE?.2p@ BDA>| ci;z=q@cB-#A>>5?4r@ BlA#>aq;8=f?րs@SBAJ>müc=C>럎t@)tBVA>i;=jfB=UE>0Lu@1EB?5A>ox?nv@3|BTpAJ>/=1>w@8UxB]AfF>`=#=D?nx@9*B=0A>ǦQy>9ny@<BwA`\>v=5=>z@=OBA`>=p3< {@>D4B:}AE>Mοo>|@E;LBB?A>G=?E^}@KÎBXAO>e{N<.=l>[<~@OBAxA>q û>~@RsBAA>k[80~?S9@ZPB=A>QlBE5A>\=;?A߀@uB?VAO>m?WTa?5(y@|SBfcAJ>&QU=b?@4BE=A>Z?\;?<̀@ -B|A>v~=04@?v+!@ԬBnA>l=+?w1@2BcA>>==v8?e @$B4cA:_>>m?"꜀@PB> -A>Sn@e#=?zЀ@tB1MA>cwf>/4=᠀@BkBWTA>W=w=:D@R4BIA>>g4<``h<ڣݾ8[>ϵ@`B=A>bڼQq=`@2ҜB6nA+>pQN ]?;Ѐ@70B$Aq>nͽV^%?}L@`4BOfAVC>=NH}p@hB6,A>m(QＷ5=@ -IRB4`A5>dj{?;8@ _BZjA{>Yg=0>@BA=>R<>%?B[w@B!A>рx]=0@+ȔBrAp>n=˰>b@-1`BBA.> -iCM>n|@:lB]A>zv=ҡ=K/ ?tq@?B(SAj>e:mBX?<@?7BWA>=_@)=.>ª@DTB0A>a溼}̈́ޓ="b?mԀ@N~$B9A6><\>T@P BAA>.[C=@\lBnAs>f=>YT?5@bpB45AK>wb / -U>@nhB$?A#y>a?=@uQBZ_A>-m;=0?n*@}xB=yAf>P=#(>UqЀ@\BVA>]ٻ3===j`?&q@ZBOAK\>m<ۈ*<7D>QĀ@=9B05Ax>=H}?nK?@џBA>X0>l@`BAG> O$>\@{iB:Ag<>l#=%qgX?Z@K3xB\&A ->;2f==>x@J8BjAle> 2=>F@}ƬBAŠ>WB>٦@ĞBcAy>;=J=t?1@ BfA>}n0C>>-?瀀@Ыk|BfSA>Tc={>&X@ՑV<BHARo> ~<6>@2 8BkA>n =d>lf@ۖBD`;tW=吀@B;Ay\>R=IhG>}@ B-A\b>S=Ҡ>z@%_B4A%>< -М>z@' B,AU>/=E&b> $@.^BF'Aȷ>F˽Y;=Q>@1BFA> =pb<?Z8@8`B=rAC>tlB(:A(z>Y< ٽCI?k?@J(xBMA>mj<=L@K0Bx)A>];f=?g#@NB>AJq>< ۍ>5@SgBYJAi>Vo=(>fm@U!rBA>=gɽ/?>@fHdBVA)>*<*s=?#p@i?B\Aݣ>XdN=AA=ucЀ@lʧB=AQ?>PX<r>Z̀@x ~Bi^}>S%j>ҷ@x"~TB9A]o>Lk=荼>@z,8B AU><%@>m@ըBGAg>Yp&<б>Zl@fHB]As>K6=D>JV@?4B]iAQ>k 藵=G0t=iЀ@}BpAW>kC=0="=,@BRAA8>s<>xr@BjAw>k=F'=P?#@@BJYA>$9 ̀@<BS7Ai>`i=]ި<$?Y%@TB3A@>W[<>X@BjA">2ٗ=3>@@@BDA`A>~=;r>L@$BjA> ]=F?v@N@BA!>s>{>S~ @,B5SA2>[' f\k?sT@PBA>Ė=m̽H>@BFAf>`=$ ;6>j @BUAA> _^ @t$B;A>v -`Sj^m<˽&ܟU>Ѐ@OOBA|>ϼ?b@׸4BkA!>u=̰`=.?Dy׀@o tBc,AxY>=G=qʇ>P@9*\B1A>Fd>F@XB`A|>=W=`C>L@BNAs>4g=3hz?-@B]A>ÿ=E?R<@8$BAY> -Rp`';q@"XB0A>㩽X'>fĀ@$B+%Av>u=Ev-<y@зB.A-Z>S@BN9Ȁ>ծ< -?)(4@kBB9A>=޽a9=IP@"--HBkYA>c7 =L>^x@&C9`BvA>=> @(%`B$bA>.輾{_i?#@-1DBTA0>@,B;AY>Xf= -E?hKu@IHB0A->̽Lg>P@`E\DB+A+>w5>@auBJ AR>cl<]ֹ?D߀@h BZA>_IR;&=2?:z@zBD*A3}>G4?@ TBaA>_=5&X/>B@<=B^iA->zt]=< ?F盀@nlB^A/>䚅=OS>@B*A{>[W=Nf/Z}=b@ȺN`B4A9>bƍ<|D>vB@tBA">&;Ὀ>2@xB6qA>>S =ý>f@ӤBHA>=xk<$=@AlBDA>`H=;ux?1 @|B~A>W:%k>!1 ?e@pBPA>y\uI<#>1 @BlIA<>.|h?=?D>@cbBA͑>L="\b?EΌ@QB6A >O=a%g>@ -%BNA>> <9<9?kp@ 1BGAף>ļ`a<>@uB*A;>h.ߦ?V@[B@A>9F?kA @ XB(A>dһ ?>6 -@Z8BMA>TH%%O @)3BShA^>z=8#<ҡ?/gQ @2>9dB8A>fGX|?B @4B1rA4>Z=yuk?Lj@55n BXA>=< =$? -@6YBAy>^b=JNW?ŀ@7BNA`>o?[=<<9>@;pBaA;@>:@5 =h*w?I@>#BJA>=X)==C?>@?4B^Ai>@=)7=R)>@I/ B*A->yhP@L B~A?>6$BTA>s~=:@WDBEAi>^`R;˸=Ѐ@gdB[A >*^=;?5@2@g8B!A>6O`>c>@hXBKA\>E=<>p@l49 BB`A>==6@mJBMuA{}>=NW<6?k@r^:pBUA >>ui= `?N@vYBA>><=l >P @wa<BLA4+>&<}C<>wԀ!@yBSA">Z;H=,?t"@zB\AW>w~*=CI^?+A&#@kpBFnAr>h;l;坝>|$@cB%A֬>zɼ%U3?U%@UByA0>C4:q=?_I&@LB+Ac> =;-"?)'@hBAM>bE1?E(@BI,AE$>w<¦[<)@WB>=;+x'[8Y?s2,@$BAw>=F>>ʀ-@B.vA>GP"c?^.@aBwA_>]76m=(>$€/@&BbAb>g[=i>10@ơB9VAgF>n=%>G1@U;xBQcAȧ>d=-I<·=/02@RBctA>p^=s?X3@ɯ4BYOAr>&޽=(m?o:v4@͏<B\A>T=R =?5@64BoA$p>s< -=?=̀6@҂BOqAx>[=H< -?[7@R$BaAJ>;H=er?g8@ BQJASB>\\7X<ڤ?3?9@ BiA>pؼ =.Y?\r:@'H0BeA0>xT<܁ ><@NDBMAY(>}= [<ۨ=@B40Aԕ>7Y=?׀>@iBX%AJ>a==?2,?@lBJA>=O<^f?W<@@`@BXA|>k=%>"bA@ǥBO!A/I>cW(6O?l}ҀC@2^0BLA=>fA<>fD@3wdBVA%0>ja0<k=v>7pE@>"B5A>_hWP=IHF@?ٓ@BRA[>i=NF=2vR:BDLJ@TDB_Aa>O=d=XV>r=$𽬊=&>,M@f.BE8A>To;?|DN@j8B=;A>k=;+>AȀO@o^B>A̸>!L>-P@tlBb Aj>=,{=i^=lQ@}e1BAS\>*>>z?6%R@Q B=A>>Ji"?!S@VhBl{Aa>˧T==jW=pT@[lBK0A1>f.4U@B&(A>uKCR?NV@$PBAe>n=#L?xW@B`7A0b>ҷ<_a=[+?BэX@ƜBMvAf>1Q<>Y@{B7A,l>? LZ@* B'AN>;{E4>Ew[@VBPA>j?7&dm=9t<\?-y]@TBmAp>=D< q^@BCA46>YD;t?IȀ_@0B"A&X>RL<j?${`@I B5A%>>ټ^=a@Bi@A>SJB m=>0^b@0@B>`AQ>F<0D ?0Yd@SBRA?{>N<,<\,>ee@$B,AQd>j<[!)?P.zf@jBG{A8># <0P<y?G՚g@#B-A:>l<)=h@ƇB6SA>rZ۪J?c-~i@-n B )Al>}X0>\|j@.dB5A>{rm5jv?P)Hk@/ylBjA}%>b=Qg=Z=wXl@6MBAF>ỉ<>?m@<8BVA6m>ep<=?^mQn@UNB=XA2>L{f51لp@f$BUA>Z==oq@lB_ZA٣>Z=-=Ud=SBr@apBKA>}d*|[?t@XU$B[AXQ>@==;j>զu@olB}A[>p=1=͸=먀v@YxBF'A >$i=c;#>.娀w@$BUAb>z߽2-%= ?c=΀x@>dBRAN>x~y@ۻy,B<A9>bڹB]?z@RBQA>mu<߰o?j@{@NBFAv>dE2 ;=D>ހ|@B&mA>QP?)}@YB=IAw>i3?'ـ~@BOAZ>Yn(<,*=@0ŜB%A4t>}$L<<=? 'ɀ@:B`A4>e - \=`? \@:LB2AT>o4= -I?SҀ@@DB;A1>a I]\.>:@@ʰB_|A >U}=6=VX?M@NLHB7As>U=wݼ?y@V/HBR-A2>a=<)=|@YTBDZA>kϮMm;T?@^BEAP>ć>=w;B>7@` BQAq>g׾J5@hqLBQA'>L<<?7'@|tҴBbAZ>!j=n_?Bfe@BBbqA#7>W;;*=luV=Ɇ؀@g`BGcAL>b=}<?e{@4B)WA}>[W== ?r@B)iA>Wc:?#*k@STB&A>PE= N>@ؗB A9I>+?5@mPB]ATp><=J>f@xB\$A>ty===?aX@xB@Acg>g#=QPP?^\Q@SBsA>RL==Ec>H@ RBAe%>i]P?I=W,@VB^A>=O=Q M>;:@G0B-A>80!Y?QC@٦B=A/(>L^Rɼ ?#i@A:BRA -> l=F<>&T@+ Bd{Au>„'={X=l`@Kv(BDA'>aJ;o? 髀@BI@A->|h7I<u>;@~B-Ax>MK^զR=(@BBA>}I8ֶ>^@9,BA(>!<qӨ>>z@y~BWjA'>fhp=ˠ>׼\@ ;BQAd)>Sڙ=^yWqg>"k'="@BA7>~=罞e?CB@ -.eBuA'>w9b^&>ך@WBJOA6>'(<u<[@$yxB[A>k.=:m?QD@*BIGA>~;0@ @?E0@,BTeA>$a=T?O@1DBTPA>=zC=?kAɀ@5v|BKEA2>ު|@7ɤBU Ae>핽8= >a@>&`BA>|½? @Da B=&A_>`2=qr {?a܀@UlB">A>ea$=Jѽo?HL -@^ZpB2WA>wM=>Łd@hLB*A>Z딽"͹1V=W@@}1ݸBmA?>TQ4==!@EBHPAa>]<*.>+_@BNA>[ Q<Q?~@JB]AAȠ>|=)=E߬?Tx @@BR_U(?r/@@BIAЬ>\4֒Ԁ@Փ> BN?A(>Pz<"ʩXU>;X|?r@$B$A>]=P^?K+A@ʙPB-AM>{DQ<}hR>@'BbOA>uwn=k?@ý@20BNA>8:=XW?ހ@aOB fAz>Sl=qe?FY@+uPB!ZAw>T3޽Wbvg>Z@ 9BiA+>K6<:}w=9P>>X@䓨B4A >ǻxЂ>̀@8BJVA/W>þ@?TB=lAk>sq}>C@-lB'Aq>xr'J >ٻ@|B]A>z;=I>>ePĀ@'-B4AG>]<γΤ?T@)|B#A1>J$%ea?; @*aB;rA >| S?/>@+SxBsA>s,=a=0>Ld@/lB7 AWK>=%>@11BdWA':><1=z;>)Ȁ@:'BA@>=ߡ=G?@S!=a_?!@?)BHAu>WɌ=<7]?ET@@>BrYA>έo -3=<>(2,@GB>VA@)>W8<>Z`@Sq0BRA͕>= PAg>:vCR>a/<@rֶB -ALS><߼H=?6y@wbBvA>ijJ===^@}[wB A>ev=yj>FT@~vBKA%>f< -Ey ӽ~>~@B +A>N3t>Ɇ >d\@8PBAF>'<ӽ{K> @BRAH>$mk9`@B_2AHc>E<ϳ/=TY>S؀@)B*OAoN>=5)3>_@YBN1A">y;<>Z@@LBSAD>f:Tj&ڻ=Q2%>@0BO6A>:=<>zD@=(BONAOb>Lr+ݽx?|@B:AAzT>#&v>0r@܋BHA>nK<8?q@VwLBRA>WX= -V^;̇>4@$B@SAi ->|pqɞ?f+@nBe{A>B-=pށ=q1>t@1`B5A*o?>.4@[1B-À>߈7Xv>!@ MPBW0A#5?;8;ͩ=?]]@%LBKA >hd((<%>a4@(:B7A?b² ->zc @(B5 8@2weBAwAF>z;ym4)?/2c@2V9B_A>UYu=X> @7BB%Ah_?!D<;@:#BVDA:?\*|=0? @<8)B&A?L?w@@/]B&Ag?#5|(wνN癢תz>38 @CB)AH6>̼a0 @H(BX|A?=M="`>ݢ @I]|BRA0W>]<^|^ @L.B'A> e]E?T @M$BJ A.>Fs<6. c!>ހ @QȀ @QΐU@BYA>ҼG=\=*>$ @R#_A@B8A9>T4 @R0|1B -IAsw?3\=>@О>čt @TBlA'?JxQ ?T @TiBE_A?GE/;?pS @Tn5BNAa? -~>sh @WkBAFN>S<=S0 @WB4NA3>6\ @XB)Az$?D=KYy7R-?/%& @X'>BKAmS?c=1<4"?E݇ @Y=5E<; @ZB;A>=V -~> T @^@BCArY?`=;:?R) @^9BSA> ֻ&@@ @^N@BpA?F-{=>9ƀ @_9΀B9A@?9$Ju˼`?T; !@__)B]A>Cz=IG??P "@`Zp BFAW?!t=lGP;Ԯ?O/ #@`/B5PAP?%?=d?6Ā $@`"`BWAא?I=> %@a(B/0A߸>H3>ހ &@b `BTA ?;;=> '@b/ BFA?#u=fY\;M>: (@b B3E )@cH*y@B2Af?YQ=Ѐ *@csa BP|A3?X +@cB"A(?;qa>]Ѐ ,@csBDAMP>7? ᧀ -@dѠB.A?+T  )?KL .@d~BA}A7L>՛O9A /@e=ĠBAa>ü-'"= 0@ed#8BsA7?a=?c 1@e|B0A? =ρh/? ^ 2@em@B?A;? 5@fBiAb>[=mp2<`?E 6@f!B4A>?wE?a 7@f BH>Am?g2<(>I 8@fC0B*AH>-zt^?1\?_1 9@g7UBA>ڥ=xڽ\`>~T :@gV8aB4A.?"=Z~ =@h`ߠB6A?9=xAս= >@hQBGWAB?@ 3Ā ?@h' B'"AA?3<Kd?jg @@i`vBBA>X!Y89Ψ=Sb A@iE B=UA?g=Z߼?XB* B@i BzzA?=ge>A C@j>aBA:?*<ؽt>l D@j|B3AF%? <ƭR&?~ E@jr BIA?>f=oO;?`0 H@lH̹B@A?;G_̓>[ I@lVMBAow?=6 Q$@?;j J@lmGB.AAx?y=H{>?QP K@m0BEAC^? qR_㼹v!>3؀ M@nBAA> ׼A?bs N@n -B;BAl >}=/??\ O@n `B:SAQ4>0k< -tOj?*"& P@n$/BL4Ai?69C<< Q@n}B]A">?;Ḟ=H(>4 R@nBcA>ҏO==?aGE S@n4@BDA6>Ns<;*?UY T@o B+A?3q=n@*C>?:݀ U@o$ϠB^BAaE>l==M^|>^: V@oiBWA>= =_?UF\ W@ohBilAMO?3A\?= -FB>̀ [@phB2'AO>#ę> \@pBF(A>Ϯ0/;}>3Ā ]@pd`BFeA? ::H;t!?hs ^@q7$ʀBZAz>#u=.V?8 _@q^B8qA?=u[)?0 } `@qoPBG1A/K>lQȼU_nP>9( b@q~sBOA>*?8B<<)I?n}A c@q%#BEuA -?Ռj;?k\ d@qӒB2@A ?ۉ4?w e@q?B5QAj? ;>o=>]K f@q{pB(GA<>9ۻó -BK?0& g@r5BGA;>rP~<>̀ h@r -=BZA?l}=2?`: i@r"B@PA?1FMD?/ j@r/BTMA>'vp=?Bve k@rl+BAUA? --ܽu|[? l@r1BVA?>2q~=qL?9a m@s\:10BQA?#}[;V<ޗ?8 n@s|v&B1vA?i;+M/t?t o@s RPBgA>6=J7=?{y p@sxBlAd?=L=>U^ q@tw@BH!A?.F<%=>F r@tBcA8?;<=v=e@ s@t'B2A>Npύ "->p!l t@tގ`B<A~2? ՓH@Sn>x u@u`B>CA9>,"<>F< v@uIBBAn$>齦WN9暩?9| w@u%"`B(1AC:>< #CA=h x@uN`B#A ? l,{eDI? y@u{g#BRpAԯ?:U+<덝?Oh| z@u3`B qA:?<ܽb;? {@vAiB;AF\?pg<ǃLk`>r$Ȁ |@vYPB+JA ?&ooy,y=( }@vbB6HA> ;'<Z ~@wS+B@YA > WA= @wTB0mAo?7U;cs{= @wVB,fA8>y?L+$]? @xBLA?#Q=S(wj=+ 3> @x+GPB@|A!?%EfŻv?9 - @xgwpBE!A>)˽x;tY> @x3B/zA? w=S z0?In @xBXFA&>νQ=!Z>ݔ @xR0BMtAJ>T?,=I>I. @y<BNZA>fIϽJY=>v( @yaEBRA> ;@V<>ک @y`BAA)I>i<&tKk=t|0 @ygqpBhA#? 1=~=> @zXF00BC>A >;T:n>R @zl`B|#Ad?#ި= 3=uw?%: @zu'JBOA>AԽm1h+;f>i @ze B1AjE?x<=+^=X @z(OBhA?z.7B=`>ٖ @{BeuA>g4=td>uԄ @{ BA;U>7 > @{Lsi BWAh?w@==j؀ @{PB)A%>[!; 8v>.P0 @{IB,A#>>J;v#> @{BEA@>;ȗ? @{pB6A?hɮ?R @|NB Ab>@s> @|1$'Bo AD9>ᤔ<'=>=P @|VB!A?[&Ss[>v @|X{BHnA?y<.,>L @|mBd A#>WsJ!=xh?8 @|tB+AV? _=_'<N@ @|B[A?-Z=9C/> @|BBvA>4:~:d=Фp @|`B_A">;\=V ?_ @|&BZA>ؼG=4S=\ @}@B𕫼,Ab?8 @}\\BDfA?(vF;ZV?fjĀ @}[ BEFA+??"<6;r?*O @~*楠BSLA>6;Ү @~:`BPAR=?H<,M@?+n @~ڜBT#AO>m<`=Y?E @~ǨB@Ac? = S> @~B+A>$=\S'> @B;A&>ڼMf= @CPBVA>%t=>ď @&QBuAr>R =0?Tٟ @EI0B-A߀?+:X5>\ @^@B1JA??<|s%N>q @m`B?AC>Eb<ػ]?z3 @z-BoA!>M=x?{< @xBGA|?v;-<(>DŽ @6~BjAf>=$=B>٭p @c˰B#SA[>+E;[bg> @hBCA>[>;{>|< @i%^BLAkb?`A4<ɨ>% @pdXBQ/ASe>n<<ݹ?-d @qXBA?%͙>"C ?+V @|8BQ/A>S|<> @ hB]kA'}?$;g?^ @0B%A{>g=NY2T>b @$|(BA> =?>#QI>* @57PBFA;.>b<<!n? @BHzAF?(ڽ4]ҽ&,? -? @8B5AD?S꼌׵=2.p @ .xB1kANH?&S<<ʼ=V?2 @!59(B~LA?==k?~Ny @7KHB6}AQ~> <M?|tӀ @ZQ*B$AtA>{4E^F? @vB)YAK>׻ڽ;?`P̀ @E@B9A?E.='3> @{@B^A^B?= =Rj>K4 @hB8A9?/t @BJA><Ÿ=b=GT= @9`PBM5A*?Gp<,Ā @FpޘB0A]? -J=2B?J# @O250B"Av8? P=C,m =dX=P>8 @4HBIA?H}6NMeN>گ - @H0BaA{?F==el>T @x+BI8AC?r<ϱ<> @BA2>ګ;!˼Z=帀 @WY pB,$Ahd?.='ب&>?O @dB:Av?瘽ct>FP @rBAP>w e\ @5˜B.A ?;Ƶy?e @BVA>"L=G=N?_ @IBNnA?9d:nr<>% @MB9AX>֊=輁= @uB(wA&?HK漑AAj?:W> @!18BnAV> -=<>d @ nB;OAܜ?rxW#?NZl @B4A>.$*> @4DB.A4>A">>ֽ8<]>6hl @B7BSA?Ҽr<7>S @ׇB'A}?@{ZԽLS>ր @e8BTGA/?BS<["|=c?4 @ B"jAK>= -n3> @3ѽBKfA?#Gּ<P?)ʠ @BPA+? Rz`r|k&>gt @fXBDA7Y?# @;B!iAW ->H=+u=a5 @=pBMAK?#<]+< =Ð @oB;vAS\? @se,hBdAQ?c<=~!= @wBH.AA?Hx='r=>Ά @7BA>րt+=洞>? @8B<*A.?L=:?n> -@_xBEA?:M=`~;™> - -@9B5^A#o?*%;Oï>` -@ (BlA? u=r*=a?#} -@}BnA?aLqb >0ހ -@ {B3(A?7{=D?3i -@T[@BDA>=;>T -@CʸB^NAI?\Z=M*?l-T -@2#BRA?&w+}<4?{Ā -@6$BH -A5><D<"?$ɀ - @HsXB<A)d>֤#Ѽ@>D - -@ZWBHHA}?"<)0=G( - @B>(A">;ź}= - @3hB'A!?";ҽK? - @ZBRA K>e`=U; Rl? -@"B0A?L=dr'> -@-BTA7?&<^==} -@9fxB>A>?k ?c€ -@?B.,A?2ν>>B -@lDшB@Ax>_i -i"8> T -@B*bAd?/BrJ&3e?I0 -@L(BWA?Q6=@?_. -@#B@fA?(^gaC>Fƀ -@BdUA(X?ͩ<=zt`?(=; -@ B#*A,?===h?~ -@LBNPAo? =6B` -@vB6|A8? ~,=5 ?J' -@$BbA$h?$X;2=mj^?@2 -@tVmB&At?=0M>X -@KBQA[J>&= <-?,v -@OHBA&_?H<> -@vB>AnQ?ؚ - @kxBIA?=J) -"@RB$ -A'?$9<ޛ -#@ A%B,`A> $u?XHԀ -$@*xB A;? '<ؽz<詀 -%@-ɏBWAg?=&n=ٿ?K| -&@F0BIAc?|o=K\=/lmA=,v -+@ mzBGwA|>==[<?p! -,@+D)XBMAV$?=R -/@W\BGA>q<>I& -0@bTBjA1L?qd̀ -1@k@BEJAg? -?=&1Q;$>8 -2@oWBA>> l<$7>K?v~M -3@|hB gA?2I6<&?Z5 -4@B1A?D|Ee >f -5@یBLvA%?@~:M<ޠ=v' -6@B5gAyl>=Ix,>ʀ@ -7@cxB<A0? =;k9>œ -9@3ZB,A$>ۻ"<Ƚ!b=00 -:@=bBS)A Y>΃<> lqv?e%Ԁ -<@IH8B@AD>kν(n_=ȬȀ -=@VFB2PA? -2nL?YU ->@p@pBQA0>ڗf;PD -?@ (B;A5&?;ׄ -@@WB(AB?$RQ<˽== -A@qB@AI>K0 -C@}BPAG?=)<ѐG>t -D@ hB=A?.qN.?h" -E@5 PPBX5A!?= = m>D -F@:B=Az?60XbQ>R -G@qhBFA?gbsmX<f?zn -H@r,8B6AE_?<Òؼ?R+ -I@v*(BPA?>䜻Ͷn<ҕi?em -J@{fBorA >`;c=t?m -K@eB3A~>-%ۼá>G -L@h0BENA$?J<;?~â -M@lV(BKA&{?<*hP -N@nGPBG.A";?,O;<k -P@Z@BgHA/?+|.8=9>6\>ah -R@؞XBQAF?%d'*fW<[r>* -S@DHBhAs>̦i=%T>c| -T@BxA8? -І<7=ɐ> -U@ABMA>]O<=( -V@YB(A ->=a&>0> -W@"~V\BrnA?'=|^=m>? -X@"\BEMA&>{=٫;{> -Y@*HB'OA6>8NPJ2>~ -Z@6BUA8?<&=[?~ -[@=uLBVA?h=?T -\@YT BLAX?2A.;;X =@ -^@`mB,TA?Ӥ½$?'W7 -_@phB5HA, ?%'<<^s?  -`@q$fBKA,k?㴽S?i -a@TBW A8?c=?Vŀ -b@BVBZA9>Ju<=1i>G< -c@%BK:Adw?BA=9<#?Nw -d@x=LBJA:? - -nK>n -e@19BJAwh>֔=Ej<=?' -7 -h@]B2A >===M -i@BgAp>˂x=Q> -j@UBM)A?-;F<8k?jw -k@0]B9A>ׁaļ? -l@ DBIAWz>z=<\A>7 -m@]BJA>ýGB<_\>mH -n@`B6As>䴾Ｍ2#>tz -o@&B"~A?F@Omi?6V -p@T1B=A>:W:?'? -q@"}BUPA?v= $?|3@ -r@AbBOA>\@mv<5}=Z -s@BB BAP8>їZ>+R?@\0 -t@\+hB$Af4? -=#|Ž\6?& -u@\BMAv?)=g<% ->`2Ѐ -v@n,B<)A`M?2==f>7" -w@t6PBVAl&>'6=/x2A=qP -x@~rB$Aw?t=FĽ\?.uR -y@BTfAE?Gg2=W?$m -z@?PB[A->ؕ@o鍰?gp -}@pBA3c?t= S؀ -@M4B@5.>=L+? -@ -]DB/AS?69=J -@ mBQA>Jؼ<": -@$"<BiA-B>)㯷mm=簠 -@EB(TA?0o -/B=ڕ0 -@MBEA?kK;]> -@RlB1A>>o=2pA> -@cTBWTAr> -=!=<=8 -@dtB\{A >ʔ; =@?>N -@gdB%AÝ>=!W*=@ -@|BAPA$?:X=sQջ*\> -@+YBAP>cȀ -@XtB7A?"C=k=h -@skBP0A>D?<)<= -@4BWA͐?0V=d=]= -@qPBVlA?2;"qK=R:?s_ -@B,~Ah[?IO='{#?߀ -@ݧ-=p=?<4 -@0#BR+A[O>?*΀ -@g?DB7jA/?Hd<\ TW?j -@B*MA -?F]jٳڽ3?kZ -@40B]A?Dq^M=H?JW -@+3`B4 Al>Ӆ=SXe> -@-TBVAL><*=u=] -@5BQ7AP?5*A -@EHB9A?N=3漃M?_ -@QK_hB0AM?/`=g=r5 -@ BRAT>+&3<>ƀ -@B{XB?PA1>N=z?YKC -@<B^9A?(;=M&`?Dx4 -@ ;B?A3>=}p໧V|?#> -@R(B'"A?1q;fKm?} -@q BYoA8>h=)h> -@$HBFxAt?=v - %;[?WU -@1BXjA>j.="h?bL -@\GB(A6? W=yAD/?#ր -@iB7KA>X=Z M? -@j/B%AQ?Q7nT^>e\ -@BBA#? fX=}Pg9Q>@S| -@f+4B`AZ ?2,= *=^>?p -@BsSA(?1UR =ڮ>W -@BA>V&O2?a. -@WB*'A`?1 Ͻ5&?, -@BmKA7?^=Ǡ=_> -@lWxBH,A?=x]<&nC>> -@ۦBBAG>pNg%:<>xn -@(jdBGA(>W=wZ<?Q -@3GxB8A ?JYļy>Ķb -@:NBA><=nt>, -@="B=A1-? [k -@DqB?lA5>̨ -<>v>r, -@DvB|A>3E=Oo>8\ -@E=xBOA#^?$U;C;<;> -@S"BGA?!Fan<MU=cȀ -@n GB'Al? =0oGޙ?_@ -@nBe.An? -=iN?cOƀ -@w.lB]AL?<<7=J>Ҁ -@xB' A/?9ǂL=B -@{BcA.#> -@},B1AD{?<<<?g۹ -@CDBAZ>ο= =$>ڀ -@tBBeAf?J5?#b -@(BJA? :<^Z?U= -@SHBWA ?;G=\>Ȁ -@ԄB^ A7>r; `=K?h -@XB5A?"85ڼ4?:՝ -@jBAAr? >LJ -@BIA?(뿽;@ -@NB;AO?lWݼIR>k -@B7A?ج=h뼪?; -@B0Ax>ٖ=Hq>H -@pQB;AJ>]<1 Gn>Ȉ -@ 9@B#%A>;Dfi> -@$B,Aݣ?R#y}= -@)%hB-A{>$=N$">JN -@C#|B"AJ>9<`l`>v -@H&=B#eA?0=t̽g >Z -@Yi-B@߆=v?Ȁ -@\J dBZAL> ݼa\=3$?z| -@e^B$cAK?T<Ƕ_?1^ -@i BITA5>ۂ>ѻ=>% -@sФB?ABW>S\=Oik=>ϓ̀ -@vbB/A9?;, A=yn -@wB8An>຤? -@BSA4?i)<~;=A ?FW -@6IBA?JaGCL>N -@?BM]A ?;]AC -@(hBJ}A|K>]q=Oz؀ -@ݛBWAM?0<=V > -@8BAE*?E<ͽ=? -@XB:A>ν=?cd@=4 -@+=(Bj"A?)xi=#> -@/3PBGJA#:?Cd;< I?gK -@/:qB1~A?R2o>o< -@L!HBA> -𭽣Q? -@QeBl!AR>=:>p -@X`BLA?- -ԼE<|>* -@\湸BSAn??XZ%=?E孀 -@mӔB1Al{??B -@y<BehA?. -=@?~ -@BVA7><7y= -? -@lBxA&?H;W>=>A -@,%BPA>?/=E;H?+p -@dB]AT? i-=D=Hi?BE -@BNA?CpF<>Tv>Ԁ -@(BVAw?&=D>V -@\tBujAV? -j==s+;c -@CxB9Alp>i ->Zڀ -@'B8A7.?45gU8W=:?A; @aLB(AU?1L'=>?L @bBCA/> -=; n=Xv @{lTlB-A?PZJ?% @|#|B;A50>XQ<ͼHL=H @}'8BA]?=Ygږ? 3 @~8ԼB&AA>w3%5=H;E @WB{,Aտ?QN=I?E! @BF/Au?Q.hQ;7?M2 @PBFA??qc;ZV?S -@BDBAg>;.Ay?rH @g'8B@4A >aG:: @$B#A'>I<|e? @xB[AV>fD=m1=;L?Bt @۟BjA"?@vW=VB @+4hB\IAj?"ܼ=>?S/ @>UB+eA?=E6+ސ?gŀ @eGHB;ZAD??G `m7U&?6s. @fmB&A?-ȼjO?A Q @pu<BHA?h<6hm>pZ @ŠB+A -?Ņ=/(?U> @ґBAHA"?'S; 旊=е<,2?ڶ @G\B;AB>N<üOΥ=^ @C B@A^> rSet?!D @اlB?0AY?(=ދ>) @8XB4A1>⿼#ʼ/?q7 !@UUBFAi??.yb;,?a. "@h@B Ab?9l=FC?V;ǀ #@BPHA?FR='h?pˀ %@1BFA>9;=Z?? &@?~E<B+%A;g>̜b<-B?Q_ '@OBTA>䨻H=_=P (@YB@jA>=cd~=IԠ )@f0B5A?*n> Vt *@jiB/pA?ν ?no +@mBEA>ٔ?;d? , ,@{$B:A?|Rj?]u -@ BwAo>va= >Ƙ .@UCDBaAPT>K<{i=gb> /@BbAW1?D.=w=i`>] 0@SHB2A ?ͩ+je?3H9 1@B5 Apr>2zȿ?$R 2@sB/2Ak>>=}m?` 3@Y4 BXtA> ¼uFG="d>1N\ 4@<B6jA?R=x w?OU} 5@* BA>J=wrY= 6@,B/Ao?ɽ -,>Ҁ 7@DBFAi>>򉽫 9@ wOpBOƀ ;@eBYeA!?$8=R=)J{>ހ <@3\BC=Ag?w==&Ce:?a΀ =@= B)A?b>H<7>$L >@=BJlAo?!=6m4 ?@\<BIeA ?ChK{BQ< B@u6B$A]-?@=n^a?! C@tB0A{?+M=No6?CeM D@`BAAQQ?$[ F@.sBMA ?C;B=2N?Hyp I@PB$A ?!\>y J@BSj` K@WB!Aa>BrȽrBG?oz L@ pBLA>cрT N@kBMAV?;M= (<>AԀ O@8BAq>e=:L罏?G#c P@!>DB2%A؞>־ۼmؼ?]7 Q@"{XB! A>t9:x>F R@%G,B-A?=!=Q?X? S@]3BRAn?=5 <9> T@nB,9AC)>y* U@}xB%A>ۀ½U> V@=LB_A>j=WS?$! W@2BIAK>c򼼾n@H X@%TB?A ?sym>V Y@B? Avb>}M@?s Z@|BQpAm>;><ܨ|?h_ [@.B -A?=&[ -> \@BTUA?t =ا?(A ]@WBCA R?/c!:|> L ^@0B#A>P;QHb?Z' _@kBgMAt?'[=S=l `@BYA~?/;e=-?E; a@4BDA?=K;e>$ b@,B9mA>=*Oo%?:% c?40B&}A?=o>P>>YP d@mB*nA?-]3?t e@{B9A7w?<#=p@ f@#"ݰBLAa?<=>/ g@$B)hA^? x3 :>v`Ԁ h@)azBDwA8?R4Ox i@5}B@|A?R; -vw?{ j@8]B4A[k?Rsʼ9>r k@=@gB[A -??X=5d%=S l@=KBUA;?Bَ m@@HBBA}?t-9f7>l n@CQBH^A?W;<,U?G o@D%yBQ A0?w=;vW<ּO=ͥ0 p@HсBBDA?i ^+;`=p q@IfBAPn?X<ཏ0?J r@KxB=An3?K[=3p "?o s@O=B-A"?;u;>, t@Qyy@BMADg?jOzަ^ w@W,`;BBAO?vX<:,_?$ x@WE@B8A#?^;?F>?B y@WV(B1AN?<n>f z@W@BA -?kyi!\? {@Y@BFA(?Db. |@[ aBWOAV?=5=ی>> }@[G1BLAT& ~@\Z -@BA'?<t꽦&)?\% @]ʮBdA\?R=sw @_KB+zA?k;/k+6>F; @_ß@BJAl=?k=/| @_4@BALA?vO+b=<؀ @`ׇBTHAR?VDVbT @c1+8@BA?tR.U;ɀ @c6X B?CAt??N½;?F @c_R BZA ?SLy=2A=ָ @d쀈@BBAS?IJ=w89=^ @f󮘀BOA?=r$P<> @g BYA ?g|;R =*H? @gBSAz?ef<#>*J @heJ0@B(pAa?u0=Aj>rCЀ @iܿB4A?åSɗ*? - @kݠB[RAr?O-K=7E={ @k^qB+FA=`?^<(,?cĀ @l[BedA?P-t3C=, > - @l]}B,AK?B<ː#H;?h_ @m B*aA?;=3D>2TD @n_`B!Aا?yY3AtH>F @ntBrkAcA?\=L=n?TЀ @oiܮBCA ?at=`; >Ȉ @o&BRA?F=q3<'?" @p% -0B,A7t? @p:pPB;AIT?~&3GS?qae @pONBh8AV?|6==?RȀ @pmBHA?`#<7q1>X$ @qCBFAl?'ݼ:;?i8 @r -qtBUxA^?h`5)= D9?3\| @rsBIcA?}sz~  @rmBGA7}?<<>z @r B3AU?MEHI D?I @sSB6dA?s0 ?^€ @sB0AY?Wz:V -?gN @tJBTA?s:=*?\o @tHpB_dAZ?hE"o=U?la @t]EBIA?_z=r{` @u iSB1AE?Ǥh]m>* @u0BOAf? ֻD<4=E @uߐB!A?%BmOs?_f: @vB5Ac?%=>D @vO%@B^AL?cA @v1BB=GAL?=N^A<4@ @w%*BJANS?DԀ @w.1BDMA?je^:;N?,E @wB7As?fi> @wr~B6hA?/'aM?H` @w:`BFYA#?}+H;> - @wD`B=A$?A<A ->Ԁ @xF?BAW?VX9?=! @xSB=TA7?y'".>iP @xyBIEAn?_1G @xBFAq?¼(K#;N>a @x_ BQbA}?U=RL @ya -BdCAY?v= -=y>l @y!RB(Aif?~=)R@=3 @yBdxA -?Sɦ={?nԙ @y$@Bi6 @zB]vAa?j<=.=Gx,>I^ @{ kBWNA#?w6=j!=?N$ @{ӠBA$?E,;-w>z @|)RB5APT?~bG @|F@BXOA5?E5< =!G?o/ @}/:cB+A$?=ͬ*(>X( @}9BOoA?"=b<݊>)4 @}|=@BlA?[u=Q?^J @}cW0BGA2?}9V<#?' @}Ű0B]DA?%/*=F?װ @}]pBHAl~?W=0<<=^ @~70B< A+?ȼA ?^ @~B#A?joex:? @PB_TA?vd =Ug^>Zy @xmBT!A@?xV$? @IPBKA ?nȭ=& @"DOB%A?^\<%U>r @O|B2A?Go뤷>$L @YOHB]Aժ?{z0=KM>N @fB+A)?z<*(t>\ @@B9A?Rn☼h?4X @(DBSAn?-׀ @XapBrAB?93<=> @ PB-A?Y-==>h8 @%!B9 Ap?c;SƼG8?M€ @ xPXBWQA?;=> @*A`BO@As?6=><3A?O½Y⍨> @LBT5A?J ==>/x @BBHIAޠ?Fj <)? @SBG>A%?0&< -l>y @nHB0A7?G<]$?9>ڀ @EhxBV$#>L @ȌBxAW?M=KG=Ɋ> @@`HBdA?p˼wP={I?6| @^!BCA ?LSJD:>) @"0B:;Am?:mw @BXA|?=4=P=#?-љ @ BRA*{?}<0> @1``BAvn?b=C+3=p @b@B*AL?Kǎ66?Hì @jB0A2?v>-{ @c0BWA?e=\YŽI=x @B6A?XJBp\?|O @(F0BdA?LIv=~?( @MBT}A?=>m @\'BLdAY?^= -duSZ @wBgOA7?{==K ? f @B MAx?)~(>$ @pBOUAi?j| -"r @oߔpBA'Aa?R= q'>h @{XB]A}?7=RP<=H0=. @|EBqA$?8˼(=?VW @XB>GA?ռ\nV$= -@BOA?d(;J<ŵ?l @сqBZAO\?K]<=3=#=4@ @((BLAl?w=/@Ԁ @0BTAe?ml="7G=:. @:c`B:A?~gr?h @Q8B_UAh#?d[='ZK=US9> Ȁ @aٸB^A]?_*=O>A& @vhB/A?':b -&? @TȐB.A?d=rI.> -J @)`B^AS?I<>=az @PB(A?w*F@ڳ>ƅP @B^NA?5`t=M?oZ} @3dhBr @nBS(Af?c=$-H< =^ @}hBHA?@bA<5?V!4 @ӃBA#?j .-> !@T(0BIA.?W8i $@B;@A ?km= ʼX>8 %@`BKA?:=k4$ )@\8B+(A?ǘ=`-r>,%x *@de B-A\?A=?^ +@BJARY?Ҥ 4@rdBh{A?,ļ]J=?%N 5@i3BTA3?|<{P:=|>`> 6@ BA? Nd&͏?Uw 7@\BD`A1%?, ;@iB&AS1?_E Nf?y9Հ <@B6A?mgӽ?Ѹ>Q% =@[B-Ah?@?GI >@:!`B^AU?c2=Q=R)">d4 ?@`BA?Y A@L7hBHiA$H?LX;9<-G>?H B@JB]*A?䉽S!N=Eg> C@\(B4A`j?Z_M=߮ɇ?b* D@4@B(kA+?<::A??Rfd E@AM?=c=\H> I@_MBI#À J@|ݨBWAF?/W L@~ -BfAb?=i\?+e M@%B,AN?fTq#"|?z K N@%BLOA1?":~ R@y2~pBOA ?U$h”H X@GB(A@D?<3>`?: Y@OBKxAMH?E/<5< Z@xKB6A? e5):> [@^pBEA?$z; -i;ƕ>q \@Ȟ -B;A"?wV;W=?dc ]@BDAL?;,?_hǀ ^@BdA?sPe - ={4?? _@! PBfA ?Zͨ=Q>gF `@&IBFhAsr?{3==;.q?1 a@=uBOBA 0?G;d<`?G+[ b@HB/.A٦?渼H=680 c@f>xB4AC?lռW? g d@llBSA?ٵWJ=`? e@AB6A ?(|?W f@yBIAa?i`C=U0?lN h@;NhB=A)?G6=P i@TB)Ae?f̭JR9TE?: j@ -B2gA1?do! k@PB*YA{?)F~A3>܀ l@ιBA?^:膽&k>-߈ m@SHB71AUp??=? n@#BQaA/*?4C t@(_BBCApE?mML=7F =S@ u@-o,B:A%?|vJ7|=TX v@1B$B:Aw?pa6k?% w@BGBJA!?WI]zDL {@ߝB2AƔ?Oۼ?h |@±B@A?y()i?B }@WqBTAQ?aBO=Dl=it@ ~@1LB*AAS?ߑ)S, @KňBHAQ?`P<$ < @1FB;A-t?Y}=I_ @=R @_O$B[hA?seֽQ1=8W?=L @kB?AF? -. @yBb9AZ?ɽ6-=jA?9" @ B(bAm\?_fB@; @BYAs?=,>:Z8 @wBJA:?dؽ@n/̀ @۰pBZA+?Y2=3t> @4SB6A?Ƽ-l'=/ @BjA$[?;A= ?a P @LB@A?kic=]=@ @4ITBA? W @`tBXA?sǽ ZQ="O?QH @rQ)B4Aq:?=9τy> @BlkARf?f<=3?,x @-B=A1?^CX @ KBkA?ʽՅ=Y?MQa @.WXBGA?Zk'< >֕Ҁ @PB6AG$?:<ʦ?yR( @YMB+A`?=m (5b\ @ B?{AA?.Vt&>f` @oB.A3"?dKAk̀ @\BF@A}?_=R;Ε=5 @$BY`A?K\=)?e>E @B:!AT?l=PzR>( @q$BS$A?g<? @BYA?͚T=+>Q - @q$BLA+^?ky<7m @*e$B4A?hTӒ>. @,lBBA"?l.;r?O @8lB>A?ֻػ?\π @D -BDA?1[f;l= @F B\Axt?tV=GN===?8̥ @PXB]AN?S-<7=H>A: @WϮBAAq?v_&(ƒ>| @`B5A\?h=Լ>R @ud"B7A?a \> @&B?A=?< =Ÿ?AD @}B/AL2?c.<ɽ s?k؀ @٫dB4A?iXν -FϨ&?dp @$B%A?4Uy= @`B.A5]?^<(?7 @|BMA$.?`;-`a= @BlA?ݤ'=? @>BF>A,?W rx;ql= @GB2oAiO?4=)R#?nY] @B=Ao"?vY=5Mʼ$_H T @ 3BjAKN?0 @9eB[A?ҽA_8=:&~ @v`B@4A?W">?au6 @pB_NAU?_f==U%?8 @Bg~A?ɼ?=?xv @nN|B')A0?<\VK6?:Ê @QsB;nAJ?J@=Sc>ƭԀ @Bf,A_?vL'=&?e- @qBErA?0W;U?t @ GtB7KA?KQ>(r @.rxBEA`0?b͎=j;?'@ @LBCA6O?o<#?%s @sɔBXAo?O==%e9?"] @~eWBXA?Q=2=#>. @\@B9A?M꼍O>dI @İB+A?SU;ٽ.>8d @gBBA?blM=9?^ @ܩB8bA?!TG? -_ @BC;A?3$N: >ƀ @]JB2VA ?OPüM?-g @7BXA?f==#>Ⱦ @% -@B:A?.9{;>J @5KHB$A_?=ڵ\?v2Ȁ @Nx#B0A?O輋 >]"̀ @Y}BD -A~?L=;/V>u1@ @^ҽ B]bA64?U<=F?]ƀ @qBWA%?=r?{ @srBJAIz?K% @'zBVA*t?aB. @u\B!A?sɟٽuGH?d @x,BXA?c:="?J# @BFAq?r=P<sc<& @/OB;Aa?YfG43?1N @BHA?'=m<<> @BC>A`F?VLd=f:4>윀 @9BA)A?zp&M?^HI @ABPZARy?ќ_0 @ -$xBb?A@?5э=k?y׀ @ B>A,,?{n<=P=Ѐ @ MmBbA?N=ip=| @$ B=A?ű|k<,@ @K4B2)AMa?<㝷:?-R @NBA?gm.>| @UBNA{?Г=N>)<>ߎ @tPBA AkX?8=-)3f>΀ @PBr_ @BkpA'?z;=>Hd @\B AA?_F<*yq>찀 @B7A -?8=^;?CFԀ @ǡB9A?f'p?#I @Ǽ$BS -A @ΡB4At?PL?oP @+BdA6?4<&=}X?y߀ @z$B`RAB?O/f=\>#F @IB#RA?0QIg>؀ @ iB6A?C(弭gD?PEs @B@Aȇ?sYC?O @vBRfA?<?H4 @-}:LBEAt?4{N ;$?+- @2RBJA\?Ǥ<uhB?Aސ?<#ת?H @C4P B0}A?M ?Xc> @I]p$BLA2?Q!ỳ @x9BLA}?U=QP<`>Ԁ @{B7AY-?1Mr>H΀ @tB?+AP?I<^F=[ @BY A? w:-=&? ֌ @,B5A"?p_HW>/ @hBHA ?gn<%Q?RĀ @|?B7AH?n꼟?= @x# B5A?;P8 -= 8 @3BDA5?{;|?m @2hBUhA?D9N=  @A BU5Ai@BV |= -X= @JB6A?젽'ܼ|?&@ @NwB@AM@HC @P@BL3A[?&=F+ @U B&Aj@-;?!M> !@WBA/?׌x` "@Y= -BIAB$@:O<)F (@bµB\A(@^0<"=& -@drBPgAa@^J=wp /@fl BYA}@ =R<=,?S`8 0@jKHBJ=A9?ZY: n뒀 4@m&Bb0A?\yg=j>.d 5@nB6A2|@ ? 6@n{QB4fA,@:(P 7@nB%JA@ Z8Y*=( 8@o#kBD#A}@;;;h?ƀ 9@pBZAiC?̊6=)=3(?؀ :@pD= BBA@0?><ҳl03?w] ;@p^eB;A)@=]\I> ̀ <@pBA A$?ڒ4<}ʻ3B.?@8߀ =@pЇ@B7uAޞ?qE/=a >@p3BFBAf@.=$;R?C% ?@q|| PBbA&?ǽ:=n}7?}, @@rXF BFA@[;%;3?T% A@r8TBMQA@ -9oV F@suB> AЙ?{Ҽ;p? G@sVBTAb@Ge;}"=!?yϺ H@t\BMA)@ռz8 L@wH@BH{A@ b=W~ M@wʀBOA7@'@< O@{NPB0A@=V298A P@{pc BAѺ@|be'> - Q@{ϿBSAp@Ee=8g= .?}Z R@|B*lAݥ@&RGѽ3>t S@~BGAC@"w<+3?# T@~B@ A8?ғ=T ?P蚀 U@~ԁBGA@k6`?<x?6c V@azBGA@w\<R?(U W@pBMA?ݺ;<?% X@N?BhAN?ӟ=O?  Y@rB;A#@;Fe?NI Z@v#B(A?O+= ?Q [@yYBLAA.( ^@ړXB8A@u? ?/e _@-pBIA@~ -B2AB+?B<Ř3>4 a@V3-BS A\@=<X?1 b@BNA?K7+ c@pnBPA*u@744 d@C%BZVA*?̯;,s=0Y?E e@ŌB\Z= f@]B0WA[@Lx0*F?/z g@IqnBUAE?=#Z? h@_B&AX?&"<սO?zz i@+AB+ABX?ك<.>Zz j@B+A9?lx'̿*|>Zƀ k@BQA @V<?ߧ=)u?N o@BHAQ?V^S<8>jD p@눢B>AŢ?b> q@B8A2@Ae:\"?W` r@4j`B+A#@<*cg'9>4 s@75BA7A+@WIJ|=x t@Ώ0BZA@aۣ=.Y+=n - u@GhB9A@('[&?U_׀ v@ܮ;pB`aA?.=) L=]8?Xz w@ t0BPAn@;.=3F<ӟ4?3 x@(BLzA?ݷԺ?<a?܀ y@7B.A5@ 8G<>v z@"BEA @IYu;W= {@]BGA@4(; -#< >' |@YB3DA!?ޛ6zL?yR }@ - BQzAw@E T?n @4oB@kA?a=f~&S?u @PBGAE?!f-8 @B5A@A:69w>t @p@HB1A@tcI> @:pBL;A -?TS]T<8 @XBVA'?勞<j=?nx @BQA@K}I` @ABK}A?=(<]B?w @ BcA?ҩGvB=u,? @^mB:A@+oO:yb?  @RB'A@˪G9G>| @2HBmAN?<=Q&? =b @P9BJAc?=x= @JB:A=?2z<}k?^k) @ 8BAP?ߦ @NBBA @"_W4>I @uBGAu?|X%<?T>Ke @2B: Ah?=~='|c?ޠ @"hBLfAN@>#<jƀ @oBBAMW@7t, @BbA?==mro?9@ @BRA?;f}տ, @-6xBJrAu@D׼6 @>RBBHAg@.;b> @DBJA_@{=A<} @m\B@A@.eJi!?Az @BHA7 -@A.<>x<0 ?&F @ˑ(BNAG@9x @ fBgA&L@3i <=·?a @:DBTA@*E*=j> - @s XBVA}@==& @4BbA@g=<=m~i<[ @B:AM@1otjJ> @BnAQ?$;=+?>O @HHB.wA"@ -;\ =EmЀ @}BQYA)@9"<, Ÿ8 @sB6WA1X@eAT| @oB2A6? 2> @7@BCSA@:8k?0. @(<BE}Az@3ڈ=L;7?_S @BLAl@1=0Ƹ @DBOAn'?3v=3Z'ƀ @_,B2A?魼'G>N^ @`d7,BK#A\s?^G=o<~dF?pɀ @aFl@BUA@Xc=~= j?etǀ @x7(B.oAU?Ԅ=g,>?W_2 @:B@EA@iNǻ> @2ZB[A}? @=5eH?r' @-}BRAvR@1=C<.? ɭ @ɫ,B< -ArS?A'NA>B" @YlBEAWD??=,;c?f/ @!BCYAk@3ZI2:?o# @.0B<#A8@>x<^>{?Q @hB5A?h{+>, @|\BF-A@/QY;s= @~xBNAU@g=$<?9WC @ͼB4A@;H̨>]G @B;tAG@v@><&7 ? \n @n[ BeAm@y == @?8B\A@+oj=@?C -L @ͺlB[A4U@*<~4=5?x @PBdA?@i:e=P>` @{˸B4=AR@>6<s?K @aDBDAK?IA @ʩTBE,A?;; ->ʀ @wJBO'Ak?=.<[>T< @*AdBUsA@-f%= ?N̂ @fLB5A?j9?:1Ԁ @|?B\uAh@,P='?=@U>_ @slBJQA@,=f@d @mB;4A8@5ZW>ׂ @QBBA?ұ=9?lO @> B.As@ JX(s?ц @!pB>A$e@";X ?h @,pBVA3@++%;7Kb=w?4X @6B7A=Z?<?Gi @]BB[A -4@1~i;m?}l @ BrAR%@FGI<=?s+~ @ɴB=VA5K@ @ BchA?Y==s{p= @ձFtB\~A4/@J <}=@N> @(B>A -?ҽd╻?) @BZA@:%E=1=Ӹ @zBGAW @6*=!`<? @hܤBD9A_@Tf=;Ea;/ @>!XBDAo@Ѽ_;t? @Re\BY Až?ЃU#N='l>D΀ @jd'BdjAܢ@BмIJ={'?i @K0B\vA&@|D=@? @m\B6Aw@=n q0> @'B1AU@ =H?π @+DB0EAj@9սl?t @dppB8PAT@ J=}> -p @BMA?@-sp @7.rBLA]@,'=<>l(@@FBaAV5?==eU?j禀@v! tBCA@GŶ:t>6b@B6tA{&?-=MX׼?ހ@HBHAm@)=2<#>.@B6hA:@$i<ռ>>斀@'B?A@T b >g@1BBAL@ ˘o:,A!?5Q@1a/B*bAaK@P-=3F>o@Ob%B=Ā @b̷BmJ@r"ǀBHwA@cTHǁǥƀ@tYBBcA9@w^'<ea?*i@t蚴BHAm@Ӽ6<<>V@v&<B-AF@`/<R>^.@v( BBA @J:+<=`@v-B#Am@};~{em%?A@wgB9?A?3@ܬ<*"?=@xY BCAp@YE=8V;>@y&B5jA8@e<@ӼC> @z$PB;6A@yVZ?Ph @{ `BjAލ@P0ۼ׏=?E9!@{FBXA@2= =?h"@}8B=wA@aeU/K>c,D#@~6B%A@]=m8vZ_@=Ē@$@XFhBPA(@cesHcp&@KUB;qA@n] -VS??o'@rD BB3Aa@?Nh*>2j(@~BGAG@ZּI<4t?6v)@U+VBIIAX@˼*@kBUAc@oI= ?W+@ BL A@=w0 -d-@GB7AJi@<׺W>.@HB*AqB@f =9_08? n/@;BXA @=%f?PG0@ BLA}<@j6=Qƞ5@18MB@6A@@l<^J>F6@gZ0Bj>Ae`@~=!܀= >+7@E5`BJA+4@ZP<5ʤzA9@<cRZ>9w,:@.MhB7A@l(9N!G={;@2ֆBDoAa@~=V3;^Z>=$<@M;(BGA%@{< -3+< `=?Fo߀=@^o8B82A@a_;0>>@3`B9TAu@b -G{?EĀ?@lӜB&IA@t,Q>@@ B?A@O!Cw?A@;֚ BHAԸ@c -w<2"t?#WB@0B8+A҈@F㯼4?{C@bBGAi@OF<>< D@0B A*n@j6E@BevA@aQ=p#>:FF@`BYA2@=%WG@.BO-A&@<<<z?_RH@W.B>A@SFs ?-AI@ -pB@AE@d<F>J@*T(BWA%>@?<,=?2ɀK@3BLdA,C@M<>C<ʺ>) L@BYA@7=+_q>M@JBOA@Z-Z<˗?΀N@/BTEA[@`:=Wn?WFO@sBSA -@PS|P<T>spP@FojTBK^Ao@Sj8[[@ hB.A!/@wc;l=+8\@dB0A@6Xm=]@ۃH@B>AC@Nfc@3qB;ACJ@Q:<=vG$?Q(Ed@aBWAv@JU:mJ=ag@|BNA@;a3+i@)dBbxA9@]<)=l>`Ȁj@;B9Ah@Ҽ@T>k@EZDB8[A'e@<Tͼb<5Ѐl@BIgA+@E'<;|ypp@RSPB3AB@@U;<N>w”q@B5AR@:UQN?"r@`bB/AV:@q= >`s@$$B.yA&@lZD ?87t@kө|BPAK@|H<~zu@ylB-An3@6=3v-?Fv@GB8)AV@e˽4I;L?kw@ BJ|A&$@< Lɪ{@? B_AE@S `#=Yh>؀|@A'^BOCACZ@N<ZC}@B/A8Q@C~@ XBLA*@d<,N{@,B/Am@=2} -:?R|@[BckA @;=sX=V@~BCA_@;W~@B`"A!l@;b=[Q8?S8@rB8Aw@eQ=Fռ[?ۀ@1B9A}@]\s'>Et@3BJA)@<+(@K4BQA@e:ʀ@sB.|AC@xż?Ξ@٦@B?UA@CC@l4XB[A @lUgi=5Ч?4@XBWA=@X<=b ?I*@SBEA;-@?z<;c?˰@`&˿BZzAA3:1=1b\>c؀@g -`BSAR@ˋ|<_?K΀@m#BQABl@8<$<޺?lԀ@r:+B[SA"%A4';3=7?4@vuBGAn@߼it<>ހ@x50BSZAәA9<'=@{`BIAfA=$>A@@ Td@tKB6)AJ@pE>t@WBEAA*/_;?@@B@fA1A R'A?5A@XBQTAA(t\R<?Q@BGAd@=  @{exBTdA@; =>G?@%[@YkfB@%AEA ̻BN?:@/BLAY@ۏȽA<>@X/BWAA9:+=?Lـ@U8BVdAv@ ?=p>{j@BCJAi@G":7w= -mP@0B?"Ac_@=7>%@BPzA)xA*P><'@GgB=AB[@3<5 - >@0cBE)AiA";;M?\6@VB-FAA6R#>5 @B>AmI@J=1&;?f@XB6A>A@<缹>hL܀@zyTBVAAf=9=V?>9@BψB?HAzA՘=LDԻv>5p@fB0A$A<߽>@!B:A*RAĬ<.Ƽcga>@8ݘBMA!A%S1;F <:> -J@XB5AAJ7uSGU?Dd@$Z8BQA%A8S2<ゕ?5U׀@*;dB* -A@腒5W?J{@Q4,BH AA!+E*<"h=Ĩ@ʸBVAT6ApMĀ@"B;4AJ#@,<֣FZ3>߮B@EL+B75A@ ;~Od!<@@iNhBGA+A w<;u,<?nl@0bBG@+AhB%AAH8p4GSf?+&@hBS AI A*<ҙ<}?9@#Bm AAǼm==>H*@zBOAT@>1=W1<>T/@m0B:&A@ث%?y?䜀@5BKAEA)E @dULBB"A@J>U@B+A,'@Ц<=z'>!H@t?B4AA|ݻ/͇=aȀ@~B;APA~qT>I>>8@B1rAS7A|<խ>,f@8BPAAL-<<[@CpBNA AŠt`JS8@:HxBCAJAl3ABAI<3>ܩ@EDB8AAV?; 3?\@y8B9QAA\!ɼm>=@d<B;AAT%\^?`ƀ@DXB?A-As.*@P߈BCA0B8@dڀBKA-A_r5}UA7UBэp@pBSA_B*=8ƀ@?0BBAU^BjS=:?atk@rBHABV*<2f?fɌ@לjBBABt)69!r>̀@BDKAOB-Qya?;MƓ>pXTENSION= 'BINTABLE' / binary table extension BITPIX = 8 / 8-bit bytes NAXIS = 2 / 2-dimensional binary table NAXIS1 = 16 / width of table in bytes NAXIS2 = 1 / number of rows in table PCOUNT = 0 / size of special data area GCOUNT = 1 / one data group (required keyword) TFIELDS = 2 / number of fields in each row TTYPE1 = 'START ' / label for field 1 TFORM1 = '1D ' / data format of field: 8-byte DOUBLE TTYPE2 = 'STOP ' / label for field 2 TFORM2 = '1D ' / data format of field: 8-byte DOUBLE EXTNAME = 'GTI ' / name of this extension MJDREFI = 51544 / [days] Integer part of time reference MJD MJDREFF = 5.0000000000E-01 / [days] Fractional part of time reference MJD TIMEUNIT= 's ' / Time unit TIMESYS = 'TT ' / Time system TIMEREF = 'LOCAL ' / Time reference END @ \ No newline at end of file diff --git a/test/data/phased_events.fits.gz b/test/data/phased_events.fits.gz new file mode 100644 index 0000000..2c7ffee --- /dev/null +++ b/test/data/phased_events.fits.gz diff --git a/test/test_csphasecrv.py b/test/test_csphasecrv.py new file mode 100755 index 0000000..50f9b9d --- /dev/null +++ b/test/test_csphasecrv.py @@ -0,0 +1,271 @@ +#! /usr/bin/env python +# ========================================================================== +# This scripts performs unit tests for the csphasecrv script. +# +# Copyright (C) 2017 Juergen Knoedlseder +# +# This program is free software: you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation, either version 3 of the License, or +# (at your option) any later version. +# +# This program is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program. If not, see . +# +# ========================================================================== +import gammalib +import cscripts +from testing import test + + +# ================================ # +# Test class for csphasecrv script # +# ================================ # +class Test(test): + """ + Test class for csphasecrv script + + This test class makes unit tests for the csphasecrv script by using it + from the command line and from Python. + """ + + # Constructor + def __init__(self): + """ + Constructor + """ + # Call base class constructor + test.__init__(self) + + # Set test data + self._phased_events = self._datadir + '/phased_events.fits' + + # Return + return + + # Set test functions + def set(self): + """ + Set all test functions. + """ + # Set test name + self.name('csphasecrv') + + # Append tests + self.append(self._test_cmd, 'Test csphasecrv on command line') + self.append(self._test_python, 'Test csphasecrv from Python') + + # Return + return + + # Test csphasecrv on command line + def _test_cmd(self): + """ + Test csphasecrv on the command line. + """ + # Set script name + csphasecrv = self._script('csphasecrv') + + # Setup csphasecrv command + cmd = csphasecrv+' inobs="'+self._phased_events+'"'+ \ + ' inmodel="'+self._model+'" srcname="Crab"'+ \ + ' caldb="'+self._caldb+'" irf="'+self._irf+'"'+ \ + ' phbinalg="LIN" phbins=5'+ \ + ' enumbins=0 emin=0.1 emax=100.0'+ \ + ' outfile="csphasecrv_cmd1.fits"'+ \ + ' logfile="csphasecrv_cmd1.log" chatter=1' + + # Check if execution of wrong command fails + self.test_assert(self._execute('command_that_does_not_exist') != 0, + 'Self test of test script') + + # Check if execution was successful + self.test_assert(self._execute(cmd) == 0, + 'Check successful execution from command line') + + # Check phase curve + self._check_phase_curve('csphasecrv_cmd1.fits', 5) + + # Setup csphasecrv command + cmd = csphasecrv+' inobs="events_that_do_not_exist.fits"'+ \ + ' inmodel="'+self._model+'" srcname="Crab"'+ \ + ' caldb="'+self._caldb+'" irf="'+self._irf+'"'+ \ + ' phbinalg="LIN" phbins=5'+ \ + ' enumbins=0 emin=0.1 emax=100.0'+ \ + ' outfile="csphasecrv_cmd2.fits"'+ \ + ' logfile="csphasecrv_cmd2.log" chatter=1' + + # Check if execution failed + self.test_assert(self._execute(cmd) != 0, + 'Check invalid input file when executed from command line') + + # Return + return + + # Test csphasecrv from Python + def _test_python(self): + """ + Test csphasecrv from Python + """ + # Set-up unbinned csphasecrv + pcrv = cscripts.csphasecrv() + pcrv['inobs'] = self._phased_events + pcrv['inmodel'] = self._model + pcrv['srcname'] = 'Crab' + pcrv['caldb'] = self._caldb + pcrv['irf'] = self._irf + pcrv['phbinalg'] = 'LIN' + pcrv['phbins'] = 5 + pcrv['enumbins'] = 0 + pcrv['emin'] = 0.1 + pcrv['emax'] = 100.0 + pcrv['outfile'] = 'csphasecrv_py1.fits' + pcrv['logfile'] = 'csphasecrv_py1.log' + pcrv['chatter'] = 2 + pcrv['publish'] = True + + # Run csphasecrv script and save phase curve + pcrv.logFileOpen() # Make sure we get a log file + pcrv.run() + pcrv.save() + + # Check phase curve + self._check_phase_curve('csphasecrv_py1.fits', 5) + + # Now use FILE as time bin algorithm. For this we need first to + # create an ASCII file. We use now 2 phase bins. The ASCII file + # is saved into the file "csphasecrv_py2.dat". + csv = gammalib.GCsv(2,2) + phmin = 0.0 + phdelta = 1.0 / float(csv.nrows()) + for i in range(csv.nrows()): + csv[i,0] = '%.5f' % (phmin + i * phdelta) + csv[i,1] = '%.5f' % (phmin + (i+1) * phdelta) + csv.save('csphasecrv_py2.dat', ' ', True) + + # Set-up unbinned csphasecrv + pcrv = cscripts.csphasecrv() + pcrv['inobs'] = self._phased_events + pcrv['inmodel'] = self._model + pcrv['srcname'] = 'Crab' + pcrv['caldb'] = self._caldb + pcrv['irf'] = self._irf + pcrv['phbinalg'] = 'FILE' + pcrv['phbinfile'] = 'csphasecrv_py2.dat' + pcrv['enumbins'] = 0 + pcrv['emin'] = 0.1 + pcrv['emax'] = 100.0 + pcrv['outfile'] = 'csphasecrv_py2.fits' + pcrv['logfile'] = 'csphasecrv_py2.log' + pcrv['chatter'] = 3 + + # Execute csphasecrv script + pcrv.logFileOpen() # Make sure we get a log file + pcrv.execute() + + # Check phase curve + self._check_phase_curve('csphasecrv_py2.fits', csv.nrows()) + + # Now we setup an observation container on input. We attached the + # model to the observation container so that csphasecrv should + # no longer query for the parameter. + cta = gammalib.GCTAObservation(self._phased_events) + obs = gammalib.GObservations() + obs.append(cta) + obs.models(self._model) + + # Set-up unbinned csphasecrv from observation container. Now use + # the GTI algorithm so that we test all timing algorithms. + pcrv = cscripts.csphasecrv(obs) + pcrv['srcname'] = 'Crab' + pcrv['caldb'] = self._caldb + pcrv['irf'] = self._irf + pcrv['phbinalg'] = 'LIN' + pcrv['phbins'] = 5 + pcrv['enumbins'] = 0 + pcrv['emin'] = 0.1 + pcrv['emax'] = 100.0 + pcrv['outfile'] = 'csphasecrv_py3.fits' + pcrv['logfile'] = 'csphasecrv_py3.log' + pcrv['chatter'] = 4 + + # Execute csphasecrv script + pcrv.logFileOpen() # Make sure we get a log file + pcrv.execute() + + # Check phase curve + self._check_phase_curve('csphasecrv_py3.fits', 5) + + # Finally we set-up a binned csphasecrv + pcrv = cscripts.csphasecrv() + pcrv['inobs'] = self._phased_events + pcrv['inmodel'] = self._model + pcrv['srcname'] = 'Crab' + pcrv['caldb'] = self._caldb + pcrv['irf'] = self._irf + pcrv['phbinalg'] = 'LIN' + pcrv['phbins'] = 2 + pcrv['emin'] = 0.1 + pcrv['emax'] = 100.0 + pcrv['enumbins'] = 10 + pcrv['coordsys'] = 'CEL' + pcrv['proj'] = 'TAN' + pcrv['xref'] = 83.63 + pcrv['yref'] = 22.01 + pcrv['nxpix'] = 20 + pcrv['nypix'] = 20 + pcrv['binsz'] = 0.02 + pcrv['outfile'] = 'csphasecrv_py4.fits' + pcrv['logfile'] = 'csphasecrv_py4.log' + pcrv['chatter'] = 4 + + # Execute csphasecrv script + pcrv.logFileOpen() # Make sure we get a log file + pcrv.execute() + + # Check phase curve + self._check_phase_curve('csphasecrv_py4.fits', 2) + + # Return + return + + # Check phase curve result file + def _check_phase_curve(self, filename, bins): + """ + Check phase curve file + """ + # Expected column names + cols = ['PHASE_MIN', 'PHASE_MAX', 'Prefactor', 'e_Prefactor', + 'Index', 'e_Index'] + + # Open FITS file + fits = gammalib.GFits(filename) + + # Check FITS file structure + self.test_value(fits.size(), 2, + 'Check for 2 extensions in phase curve FITS file') + self.test_assert(fits.contains('PHASECURVE'), + 'FITS file contains "PHASECURVE" extension') + + # Get PHASECURVE table + table = fits['PHASECURVE'] + + # Check FITS table structure + self.test_value(table.ncols(), len(cols), + 'Check for %d columns in phase curve FITS table' % len(cols)) + self.test_value(table.nrows(), bins, + 'Check for %d rows in phase curve FITS table' % bins) + for col in cols: + self.test_assert(table.contains(col), + 'FITS file contains "'+col+'" column') + + # Close FITS file + fits.close() + + # Return + return diff --git a/test/test_python_cscripts.py b/test/test_python_cscripts.py index ac1c055..1daaf4d 100755 --- a/test/test_python_cscripts.py +++ b/test/test_python_cscripts.py @@ -36,6 +36,7 @@ import test_csobs2caldb import test_csobsdef import test_csobsinfo import test_csobsselect +import test_csphasecrv import test_cspull import test_csresmap import test_csroot2caldb @@ -124,6 +125,7 @@ def test(installed=False, debug=False): test_csobsdef.Test(), test_csobsinfo.Test(), test_csobsselect.Test(), + test_csphasecrv.Test(), test_cspull.Test(), test_csresmap.Test(), test_csroot2caldb.Test(),