Learn R Programming

Hadron

An R implementation of fitting routines used in lattice QCD. It provides useful functions for extraction hadronic quantities and such like.

The license is GPL 3 or later, even though the DESCRIPTION only shows GPL-3.

master branch:

Copy Link

Version

Install

install.packages('hadron')

Monthly Downloads

39

Version

3.1.2

License

GPL-3

Issues

Pull Requests

Stars

Forks

Maintainer

Carsten Urbach

Last Published

September 10th, 2020

Functions in hadron (3.1.2)

analysis_online

analysis_online
addStat.raw_cf

Extend statistics of an existing raw_cf container
add.cf

Arithmetically adds two correlation functions
addConfIndex2cf

add a configuration index to an cf object
avg.cbt.cf

average close-by-times in a correlation function
alphas

compute alpha strong at given scale
add.raw_cf

add two raw_cf objects
analysis_gradient_flow

analysis_gradient_flow
CExp

Cosh Or Sinh Build Out Of Two Exps
addStat.cf

Combine statistics of two cf objects
bootstrap.gevp

perform a bootstrap analysis of a GEVP
bootstrap.hankel

GEVP method based on Hankel matrices.
c.cf

Concatenate correlation function objects
bootstrap.nlsfit

Bootstrap a non-linear least-squares fit
block.raw_cf

Block average correlation function data
bootstrap.analysis

Performs a Bootstrap with Blocking Analysis of a Timeseries
bootstrap.effectivemass

Computes effective masses with bootstrapping errors
bootstrap.cf

bootstrap a set of correlation functions
bootstrap.hankel_summed

GEVP method based on Hankel matrices.
bootstrap.meanerror

Compute the bootstrap error of the mean
cdhnew

finite size corrections a la Colangelo, Duerr, Haefeli, but re-expanded as series in the quark mass
cdh

finite size corrections a la Colangelo, Duerr, Haefeli
c.raw_cf

Concatenate raw_cf correlation function objects
cA2.09.48_3pi_I3_0_A1u_1_pc

A three pion correlator with significant thermal states.
cf_key_meson_3pt

Generate HDF5 key for CVC 'correlators' meson 3pt function with a local or derivative insertion The key for a meson three-point function has the form: /sud+-g-u-g/t10/dt12/gf5/pfx0pfy0pfz0/gc0/Ddim0_dir0/Ddim1_dir1/D.../gi5/pix0piy0piz0 where, from left to right:
  • 'u' is the flavour of the "backward" propagator

  • 'd' is the flavour of the "sequential" propagator

  • '+' indicates that 'sud' is daggered

  • 'g' indicates a gamma insertion

  • 'u' is the flavour of the foward propagator

  • 'g' indicates a Dirac structure at the source

  • 'tXX' is the source time slice

  • 'dtYY' is the source-sink separation

  • 'gfN' gamma structure at the sink in CVC indexing

  • 'pfxXpfyYpfzZ' is the sink momentum in CVC convention (sink and source phases are both e^ipx)

  • 'gcN' gamma structure at the current insertion point in CVC indexing

  • 'DdimJ_dirK' covariant displacement applied in dimension 'J', direction 'K' where it should be noted that this is. in operator notation, i.e., the right-most displacement is the one applied first.

  • ...

  • 'giN' gamma structure at the souce in CVC indexing

  • 'pixXpiyYpizZ' at the source in CVC convention

cf_key_meson_2pt

Generate key string to identify a meson 2pt function
cf_weighted

Weighted CF mixin constructor
cf_principal_correlator

Principal correlator CF mixin constructor
cf

Correlation function container
compute.plotlims

compute.plotlims
cf_boot

Bootstrapped CF mixin constructor
cf_shifted

Shifted CF mixin constructor
cf_smeared

Smeared CF mixin constructor
cf_subtracted

Subtracted CF mixin constructor
computeDisc

computes a disconnected correlation function from loops
computeacf

Computes The ACF and Integrated AC Time
correlators_key_meson_3pt

Generate HDF5 key for CVC 'correlators' meson 3pt function with a local or derivative insertion The key for a meson three-point function has the form: /sud+-g-u-g/t10/dt12/gf5/pfx0pfy0pfz0/gc0/Ddim0_dir0/Ddim1_dir1/D.../gi5/pix0piy0piz0 where, from left to right:
  • 'u' is the flavour of the "backward" propagator

  • 'd' is the flavour of the "sequential" propagator

  • '+' indicates that 'sud' is daggered

  • 'g' indicates a gamma insertion

  • 'u' is the flavour of the foward propagator

  • 'g' indicates a Dirac structure at the source

  • 'tXX' is the source time slice

  • 'dtYY' is the source-sink separation

  • 'gfN' gamma structure at the sink in CVC indexing

  • 'pfxXpfyYpfzZ' is the sink momentum in CVC convention (sink and source phases are both e^ipx)

  • 'gcN' gamma structure at the current insertion point in CVC indexing

  • 'DdimJ_dirK' covariant displacement applied in dimension 'J', direction 'K' where it should be noted that this is. in operator notation, i.e., the right-most displacement is the one applied first.

  • ...

  • 'giN' gamma structure at the souce in CVC indexing

  • 'pixXpiyYpizZ' at the source in CVC convention

computefpsOS

Computes the pseudoscalar decay constant for the Osterwalder Seiler case from the pseudoscalar amplitude and mass
cf_orig

Original data CF mixin constructor
conj_raw_cf

Take the complex conjugate of a raw_cf object
cf_meta

CF metadata mixin constructor
create_displ_chains

create list of chains of displacements Multilpe covariant displacements, when applied in order, form a list of displacments. Each consists of a direction and a dimension.
computefps

Computes the pseudoscalar decay constant for the twisted mass case from the pseudoscalar amplitude and mass
cvc_local_loop_key

Generate HDF5 key for a momentum and spin-projected CVC loop
cvc_read_loops

read HDF5 loop files in the CVC loop format
effectivemass

effectivemass
effectivemass.cf

Computes effective mass values for a correlation function
correlators_key_meson_2pt

Generate HDF5 key for CVC 'correlators' meson 2pt function
correlatormatrix

Sample correlator matrix
effmass

effmass
dispersion_relation

Continuum dispersion relation for CM to lattice frame
disc_3pt

disconnected contribution to current insertion three-point function
cvc_to_raw_cf

Convert correlation function read from CVC HDF5 or AFF format to 'raw_cf'
cyprus_make_key_scalar

HDF5 key for Cyprus CalcLoops scalar-type loops
extract.loop

Extract a single loop from an object of class cmiloop
escapeLatexSpecials

Escape special LaTeX characters for use in LaTeX labels
concat.cf

Concatenate two correlation function objects
extractSingleCor.cf

extract one single correlator object as cf object from a large cf object.
fs.mpia0

Finite Size Corrections to \(q\cot\delta\) for I=2 \(\pi\pi\) near threshold
foldr1

Folds the non-empty list with the binary function
concat.raw_cf

Concatenate two raw_cf correlation function objects
extract.obs

Extract One or More Gamma Combinations from am CMI Correlator
-.cf

Arithmetically subtract correlators
effmass2

effmass2
extract_mass.effectivemassfit

specialisation of extract_mass to objects of type effectivemassfit
extract_mass

generic function to extract a fitted mass
-.raw_cf

add two raw_cf objects
getorderedconfignumbers

Creates an ordered vector of gauge config file numbers
g1

g1
fit.effectivemass

Fits a constant to effective mass data
getorderedfilelist

Creates an ordered filelist from a basename and a path
construct_onlinemeas_rundir

cyprus_make_key_vector

HDF5 key for Cyprus CalcLoops derivative-type loops
h5_names_exist

check if group names exist in HDF5 file
fit.plateau2cf

fits a plateau to an object of class cf
get_plotdata_raw_cf

extract data from 'raw_cf' in format convenient to plot
invalidate.samples.cf

Invalidate samples
invcosh

numerically invert the cosh function for the mass
fs.a0

Finite Size Corrections to \(q\cot\delta\) for I=2 \(\pi\pi\) near threshold
hadron

The Hadron Package
fs.qcotdelta

Finite Size Corrections to \(q\cot\delta\) for I=2 \(\pi\pi\) near threshold
gevp2cf

Extracts a principle correlator from a GEVEP
invertCovMatrix

Inverts the covariance matrix for noisy data
is.cf

Checks whether an object is a cf
h5_get_dataset

get dataset from HDF5 file
gm_mu

is_empty.raw_cf

check if an obect is of class raw_cf and empty otherwise
mul.cf

Arithmetically scale a correlator by a scalar a
jackknife.cf

jackknife a set of correlation functions
cyprus_read_loops

read HDF5 loop files in the Cyprus CalcLoops format
fit.cosh

Fits a sum of several cosh-functions
extract_mass.matrixfit

specialisation of extract_mass to objects of type matrixfit
gevp.hankel

GEVP method based on Hankel matrices.
mul.raw_cf

scale raw_cf data
hankeldensity2effectivemass

hankeldensity2effectivemass
is.raw_cf

check if an object is of class raw_cf
gevp

solve GEVP for correlator matrix
hankel2cf

hankel2cf
gm

List of arrays of gamma structures
gevp.hankel_summed

GEVP method based on Hankel matrices.
hankel2effectivemass

hankel2effectivemass
gevp2amplitude

Extracts physical amplitudes from a GEVP
has_icf

Checks whether the cf object contains an imaginary part
plot.effmass

plot.effmass
loop_stochav

average over stochastic samples of loop
parametric.nlsfit

NLS fit with parametric bootstrap
matrixfit

Routine For A Factorising Matrix Fit
is_empty.cf

Checks whether the cf object contains no data
mom_combinations

Generate table of momentum component combinations
parametric.nlsfit.cov

parametric.nlsfit.cov
matrixModel

Correlator matrix model.
make_parlist

Create a parameter list for matrixfit
loop_vev_subtract

subtract vev from loop data
plot.gevp.amplitude

plot.gevp.amplitude
parametric.bootstrap

Parametric bootstrap
plot_hankel_spectrum

plot_hankel_spectrum
+.cf

Arithmetically add correlators
pscor.sample

Sample pseudoscalar correlator
plot_eigenvalue_timeseries

plot_eigenvalue_timeseries
plotwitherror

Plot Command For XY Plots With Error Bars
jackknife_error

Estimates error from jackknife samples
int_idx_matrix.raw_cf

Construct tensor index set for the internal degrees of freedom
jackknife_cov

jackknife_cov
idx_matrix.raw_cf

Construct the tensor index set for the entire raw correlator
raw_cf

Container for raw correlation functions
loopdata

Sample loop data
pcModel

Principal correlator two state model.
make_parind

Create a parameter index matrix for matrixfit
readhlcor

readhlcor
readgradflow

Read Gradient Flow Output Files in tmLQCD format
plot.averx

Plots averx data
print.raw_cf

Print summary of data contained in raw_cf container
print.ofit

print.ofit
parametric.bootstrap.cov

Parametric bootstrap with covariance
plot.bootstrapfit

Plot a bootstrap NLS fit
readcmidisc

reads disconnected loops in cmi format
pcac

Computes the pcac mass
plot.cf

Plot a correlation function
readcmifiles

Read Single Data Files in Chris Michael Format
plot.outputdata

Plot Command For Class Ouputdata
resampling_is_concatenable

Checks whether the resampling of two cf objects is concatenable
loop_spin_project

spin projection of quark loop data
new_matrixfit

perform a factorising fit of a matrix of correlation functions
loop_2pt

compute two-point correlation function between quark loops
old_removeTemporal.cf

Remove temporal states
pcacfit

pcacfit
/.raw_cf

divide two raw_cf objects
summary.hadronacf

summary.hadronacf
store_correl

Store a 'raw_cf' correlator in an associative array together with a description The object cf will be stored as an element of cmap under key out_key in the member obj of cmap. The data frame passed via desc will be appended as a row to cmap[[out_key]]$map. If out_key does not exist as a key in cmap, a new element will be created. If it already exists, addStat.raw_cf is called to add statistics to the existing raw_cf. Requires the 'hash' package.
plot.coshfit

Plot a cosh-fit
plaq.sample

Sample plaquette time series
plot.cfit

plot.c1fit
summary.coshfit

Summarize a cosh-fit
simple.nlsfit

NLS fit with without bootstrap
*.cf

Divide two cf objects by each other measurement by measurement
residual_plot

residual_plot
zetazp

Computes the running of Z_P from scale mu0 to scale mu2
weight_shift_reweight.cf

Weight-shift-reweight a correlation function
summary.hankel_summed

summary.hankel_summed
plot.effectivemass

plot.effectivemass
plot.uwerr

Plot Command For Class UWerr
plot.raw_cf

plot all correlators in raw_cf object
plot_timeseries

plot_timeseries
raw_cf_data

Original data mixin constructor for raw_cf
plot.pionff

plot.pionff
raw_cf_meta

raw_cf metadata mixin constructor
summary.effectivemass

summary.effectivemass
print.effectivemassfit

print.effectivemassfit
print.cf

print.cf
raw_cf_to_cf

Extract a particular internal component of a 'raw_cf' into a 'cf'
onlinemeas

determines pion mass and pcac mass from online measured correlator of the HMC code
overview_plot_raw_cf

create convenient overview plots for a raw_cf object
plothlinewitherror

plothlinewitherror
restore_seed

Restore random number generator state
readnissatextcf

reader for Nissa text format correlation functions
readoutputdata

Read Data In output.data Format of tmLQCD
readbinarycf

read correlation function from binary files
summary.matrixfit

summary.matrixfit
samplecf

Sample cf data
takeTimeDiff.cf

Take time difference
plot.massfit

plot.massfit
plot.hadronacf

plot.hadronacf
pointswithslantederror

pointswithslantederror
tex.catwitherror

paste a number with error in tex-ready format
+.raw_cf

add two raw_cf objects
plot.matrixfit

Plot a matrixfit
summary.ofit

summary.ofit
tikz.finalize

tikz.finalize
readbinarysamples

Read binary correlation function by sample
readtextcf

Read correlator data from single file
readbinarydisc

read disconnected loops from binary files
resample_hankel

Resample bootstrap samples in Hankel effmass
shift.cf

shift a correlation function by 'places' time-slices
removeTemporal.cf

Remove Thermal States by Weighting and Shifting
shift.raw_cf

shift a raw_cf correlation function by 'places' time-slices
plot.ofit

plot.ofit
print.bootstrapfit

Print a bootstrap NLS fit
predict.bootstrapfit

Predict values for bootstrapfit
tikz.init

tikz.init
summary.bootstrapfit

Summarize a bootstrap NLS fit
summary.cf

summary.cf
summary.effectivemassfit

summary.effectivemassfit
subtract.excitedstates

Substract excited states.
string2error

string2error
unsymmetrise.cf

Unfold a correlation function which has been symmetrised
*.raw_cf

multiply two raw_cf objects
uwerr

Time Series Analysis With Gamma Method
uwerr.cf

uwerr.cf
summary.raw_cf

Print summary of data contained in raw_cf container
summary.uwerr

summary.uwerr
resampling_is_compatible

Checks whether the resampling of two cf objects is compatible
summary.gevp.amplitude

summary.gevp.amplitude
symmetrise.cf

Average backward and forward-dominated parts of the correlation function
swap_seed

Set seed and store a seed which can be used to reset the random number generator
uwerr.raw_cf

Gamma method analysis on all time-slices in a 'raw_cf' object
weight.cf

Weight a correlation function