mue C++ API Reference

struct cpdf : public intfunction

Public Functions

inline void init(double wv)
inline double cpandel(double d, double t)
inline void destroy()
inline virtual ~cpdf()
inline virtual double ifunction(double d, double t)
inline cpdf(double l, double la, double tau, double jit, double lf, double ls, double laf, double las, double wv, bool iflag, const char *rname)
inline void lx(double x)
inline void ll(double x1, double x2)

Public Members

double l
double la
double tau
double jit
double wv
double np
double ng
int size
double dh
double hmin
double hmax
double *abs
double *sca
double *abs_sum
double *sca_sum
double lf
double ls
double laf
double las
bool kurt
bool cpandel_flag
interpol *cpandel_intc

Public Static Attributes

static const int verbose = 0
static double udf = 1.e-5
static double thmax = 100
static bool updown = true
static bool invert = false
static const double le = 0.3
static const double lr = 0.6
static const double rlen = 0
static const double elen = 0
static const double th0 = 90
static const double dd0 = 0
static const double c = 0.299792458
static const double lm = 0.5e-6
static const double a1 = 0.84
static const double a00 = 3.1
static const double a01 = -3.9
static const double a02 = 4.6
static const double a2 = 0.3
static const double at = 1.
static const double ts1 = 5
static const double ts2 = 25
static const double pmin = 1.e-50
static const double pre = 1.e-2
static const double dmin = -150
static const double dmax = 1050
static const double admin = -100
static const double admax = 1000
static const double tmin = -150
static const double tmax = 3050
static const double atmin = -100
static const double atmax = 3000
static const bool fixdir = true
static const bool lcflag = true
static const bool mnball = true
static const bool nnhit = false
static const int om_pad = 1
static const int dnum = 1000
static const int tnum = 1000
struct cpdf

Public Functions

inline cpdf(double ita, double sig, double ef, double es, double af, double as)
inline double lnscat(double d, double t, double Q)
inline void relax()
inline void restore()

Private Functions

inline double gaus(double x, double sig)
inline double student(double x, double sig)
inline double addn(double x, double t)
inline void set(double d)
inline double pandel(double t)
inline double cpandel(double t)
inline double cf(double t)
inline double dcf(double t, double s)
inline double dcf(double t, double s, double Q)
inline double npandel(double t, double Q)

Private Members

double ita
double sig
double ef
double es
double af
double as
double xi
double rho
struct dir

Subclassed by MUEX::vert

Public Functions

inline dir()
inline dir(const double n[])
inline bool operator<(const dir &rhs) const

Public Members

double n[3]
struct epar

Public Members

int num
double val
double sto
void *par
gsl_vector *x
struct event
#include <reader.h>

Public Functions

inline bool operator<(const event &rhs) const
inline event(wform *wf)
inline void copy(const event &e)
inline void destroy()
inline event(const event &e)
inline const event &operator=(const event &e)
inline ~event()
inline event()

Public Members

unsigned long long id
unsigned long long gt
double Q
double W
char fi
int ijs
wform *wf
struct evless

Public Functions

inline bool operator()(const event *lhs, const event *rhs) const
struct face

Public Functions

inline face(const dir &x, const dir &y, const dir &z)
inline void half(dir h[])

Public Members

dir v[3]
struct hit

Public Functions

inline hit(double t, double q)
inline bool operator<(const hit &rhs) const

Public Members

double t
double q
class I3mue : public I3ConditionalModule
#include <I3mue.h>

This module performs reconstruction of IceCube events.

class: I3mue

Version $Id$

Date: Tue November 1 2005

(c) 2005 IceCube Collaboration

Author

dima

Public Functions

SET_LOGGER ("I3mue")
I3mue(const I3Context &ctx)
virtual ~I3mue()
void Configure()
void Physics(I3FramePtr frame)
void Geometry(I3FramePtr frame)
void Calibration(I3FramePtr frame)
void DetectorStatus(I3FramePtr frame)

Private Members

bool fInitialized
unsigned int fEvent
unsigned long long fToffset
std::vector<std::string> fRecoPulseSeriesNames
std::string fDataReadoutName
std::string fRecoResult
int fRecoRnum
int fRecoIntr
std::string fOutputParticle
int fVerbose
double fDetCenterDepth
std::map<OMKey, unsigned long long> fIds
struct Ice

Public Functions

inline Ice()
inline void ini(string dir)
inline ~Ice()
inline void lx(double x, double y, double z)
inline void ll(double x1, double y1, double z1, double x2, double y2, double z2)
inline void set(double le, double la)
inline double flux(double d, double coseta, bool type)
inline double aprof(double cospsi, double dd)

Public Members

double wv
double np
double ng
double c
double cm
double na
double le
double la
double coschr
double sinchr
double tanchr

Private Members

bool kurt
int size
double dh
double hmin
double hmax
double *abs
double *sca
double *abs_sum
double *sca_sum
struct Ico

Public Functions

inline Ico(int num)

Public Members

int vn
int fn
vector<dir> v
vector<vector<int>> f
class interpol : private intfunction
#include <interpol.h>

Public Functions

interpol(int, double, double, int, double, double, double (*)(double, double))
interpol(int, double, double, int, double, double, intfunction*)
virtual ~interpol()
double interpolate(double, double)

Private Functions

inline virtual double ifunction(double x, double y)
void init(int, double, double, int, double, double, intfunction*)
double interpolate(double, int)
double interpolate(double)

Private Members

double *x
double *y
double **f
double *iX
double *iY
double *iF
double *c
double *d
int starti
int romberg
bool rational
int nx
int ny
double xmin
double xmax
double xstep
double ymin
double ymax
double ystep
double (*func)(double, double)
struct intfunction
#include <interpol.h>

Subclassed by I3FRreco::cpdf, interpol

Public Functions

virtual double ifunction(double, double) = 0
inline virtual ~intfunction()

Public Members

const char *rname
struct line

Subclassed by MUEX::vtemp

Public Functions

inline line()
inline line(double n[3], double x, double y, double z, double t)
inline line(double th, double ph, double x, double y, double z, double t)

Public Members

double nx
double ny
double nz
double x0
double y0
double z0
double t0
template<int d>
class Miniball
#include <miniball.h>

Public Types

typedef std::list<Point<d>>::iterator It
typedef std::list<Point<d>>::const_iterator Cit
typedef std::list<Point<d>>::iterator It
typedef std::list<Point<d>>::const_iterator Cit

Public Functions

inline Miniball()
void check_in(const Point<d> &p)
void build()
Point<d> center() const
double squared_radius() const
int nr_points() const
Cit points_begin() const
Cit points_end() const
int nr_support_points() const
Cit support_points_begin() const
Cit support_points_end() const
double accuracy(double &slack) const
bool is_valid(double tolerance = 1e-15) const
inline Miniball()
void check_in(const Point<d> &p)
void build()
Point<d> center() const
double squared_radius() const
int nr_points() const
Cit points_begin() const
Cit points_end() const
int nr_support_points() const
Cit support_points_begin() const
Cit support_points_end() const
double accuracy(double &slack) const
bool is_valid(double tolerance = 1e-15) const

Private Functions

void mtf_mb(It k)
void pivot_mb(It k)
void move_to_front(It j)
double max_excess(It t, It i, It &pivot) const
void mtf_mb(It k)
void pivot_mb(It k)
void move_to_front(It j)
double max_excess(It t, It i, It &pivot) const

Private Members

std::list<Point<d>> L
Miniball_b<d> B
It support_end
template<int d>
class Miniball

Public Types

typedef std::list<Point<d>>::iterator It
typedef std::list<Point<d>>::const_iterator Cit

Public Functions

inline Miniball()
void check_in(const Point<d> &p)
void build()
Point<d> center() const
double squared_radius() const
int nr_points() const
Cit points_begin() const
Cit points_end() const
int nr_support_points() const
Cit support_points_begin() const
Cit support_points_end() const
double accuracy(double &slack) const
bool is_valid(double tolerance = 1e-15) const

Private Functions

void mtf_mb(It k)
void pivot_mb(It k)
void move_to_front(It j)
double max_excess(It t, It i, It &pivot) const

Private Members

std::list<Point<d>> L
Miniball_b<d> B
It support_end
template<int d>
class Miniball_b
#include <miniball.h>

Public Functions

inline Miniball_b()
const double *center() const
double squared_radius() const
int size() const
int support_size() const
double excess(const Point<d> &p) const
void reset()
bool push(const Point<d> &p)
void pop()
double slack() const
inline Miniball_b()
const double *center() const
double squared_radius() const
int size() const
int support_size() const
double excess(const Point<d> &p) const
void reset()
bool push(const Point<d> &p)
void pop()
double slack() const

Private Members

int m
int s
double q0[d]
double z[d + 1]
double f[d + 1]
double v[d + 1][d]
double a[d + 1][d]
double c[d + 1][d]
double sqr_r[d + 1]
double *current_c
double current_sqr_r
template<int d>
class Miniball_b

Public Functions

inline Miniball_b()
const double *center() const
double squared_radius() const
int size() const
int support_size() const
double excess(const Point<d> &p) const
void reset()
bool push(const Point<d> &p)
void pop()
double slack() const

Private Members

int m
int s
double q0[d]
double z[d + 1]
double f[d + 1]
double v[d + 1][d]
double a[d + 1][d]
double c[d + 1][d]
double sqr_r[d + 1]
double *current_c
double current_sqr_r
class muex : public I3ConditionalModule

Public Functions

inline muex(const I3Context &ctx)
inline void Configure()
inline void Physics(I3FramePtr frame)
inline void Geometry(I3FramePtr frame)
inline void Calibration(I3FramePtr frame)
inline void DetectorStatus(I3FramePtr frame)

Private Functions

SET_LOGGER ("muex")

Private Members

bool ini
std::string pulses
std::string badoms
std::string rectrk
std::string result
std::string elistn
std::string icedir
unsigned int lcsp
unsigned int repeat
bool rectyp
bool usempe
bool detail
bool energy
bool compat
struct multe
#include <reader.h>

Public Members

int evnum
int max
int mlt
int tot
unsigned long long t0
unsigned long long ts
int inic
int itop
int icsn
int itsn
std::vector<preco> tracks
std::vector<event> events
struct mydom
#include <reader.h>

Public Functions

inline mydom()

Public Members

int num
unsigned long long id
double hv
double st
double vem
double x
double y
double z
int om
int str
struct mydom

Public Members

double x
double y
double z
double hv
double eff
double st
struct mykey

Public Functions

inline mykey(int str, int om)
inline bool operator<(const mykey &rhs) const

Public Members

int str
int om
struct mykey
#include <reader.h>

Public Functions

inline bool operator<(const mykey &rhs) const

Public Members

int om
int str
struct otemp

Public Functions

inline otemp()
inline void norm(int i)
inline void norm()

Public Members

double tc
double t0
double th
double ph
double x0
double y0
double z0
double a0
double k0
short dth
short dir
template<int d>
class Point

Public Functions

inline Point()
inline Point(const Point &p)
inline Point(const double *p)
inline Point &operator=(const Point &p)
inline double &operator[](int i)
inline const double &operator[](int i) const
inline const double *begin() const
inline const double *end() const

Private Members

double coord[d]
template<int d>
class Point
#include <miniball.h>

Public Functions

inline Point()
inline Point(const Point &p)
inline Point(const double *p)
inline Point &operator=(const Point &p)
inline double &operator[](int i)
inline const double &operator[](int i) const
inline const double *begin() const
inline const double *end() const
inline Point()
inline Point(const Point &p)
inline Point(const double *p)
inline Point &operator=(const Point &p)
inline double &operator[](int i)
inline const double &operator[](int i) const
inline const double *begin() const
inline const double *end() const

Private Members

double coord[d]
struct preco
#include <reader.h>

Public Members

int type
int hits
int wforms
int hitdoms
int strings
double rllh
double t0
double z0
double th
double dd
double ph
double az
double x0
double y0
double n0
double ne
double nl
double t0_min
double z0_min
double th_min
double dd_min
double n0_min
double ph_min
double az_min
double x0_min
double y0_min
double t0_max
double z0_max
double th_max
double dd_max
double n0_max
double ph_max
double az_max
double x0_max
double y0_max
short dir
short ddf
short en
short sp
struct rrt

Public Functions

inline void ini(gsl_rng *G = NULL)
inline double *rot(double *m)
inline double *tor(double *m)
inline void ini(double *n)
inline double *ste(double x, double y)

Public Members

gsl_rng *G
double A[3][3]
double n[3]
struct sens

Public Functions

inline double s(double coseta)
inline double s(double coseta, double y)
inline sens()

Public Members

double ave
double a[11]
struct sens

Public Functions

inline double s(double coseta)
inline sens()

Public Members

double ave
double ped_pmt[11]
struct set

Subclassed by MUEX::vtemp::setq, MUEX::vtemp::sett

Public Functions

inline set(const mydom &dom)
inline void loc() const
inline double deff() const

Public Members

double x
double y
double z
double q
double eff
double rn
double dd
double tdel
double coseta
double le
double la
struct setq : public MUEX::vtemp::set

Public Functions

inline setq(const set &dom, bool sat)
inline setq(const mydom &dom)

Public Members

double p
bool sat
struct sett : public MUEX::vtemp::set

Public Functions

inline sett(const mydom &dom)

Public Members

vector<hit> h
double t
double ave
bool lc
struct spq

Public Members

double p
double q
double r
struct Tilt

Public Functions

inline Tilt()
inline void ini(string dir)
inline void set_r0()
inline double zshift(double rx, double ry, double rz)

Public Members

bool tilt
int lnum
int lpts
int l0
double lmin
double lrdz
double r0
double lnx
double lny
double lr[LMAX]
double lp[LMAX][LYRS]
struct vert : public MUEX::dir

Public Functions

inline vert()
inline vert(const dir n, double f)

Public Members

double f
struct vtemp

Public Functions

inline vtemp()
inline void norm(int i)
inline void norm()

Public Members

double t0
double z0
double th
double dd
double ph
double az
double nn
double x0
double y0
short dth
short ddf
struct vtemp : public MUEX::line

Public Functions

inline vtemp(I3RecoPulseSeriesMapConstPtr pmap, bool type, unsigned int lcsp, unsigned int unc, bool mpe, bool en, gsl_rng *G = NULL)
inline void ini(set &dom, bool type = true)
inline double illt(bool reset = true, bool all = true)
inline void set_t0(double n[])
inline double llh(double n[])
inline void fastini()
inline void recoini()
inline void advance(double r)
inline double llh()
inline void set_gsl(const gsl_vector *x)
inline double solve(int r)
inline double Solve(int ini)
inline double solve(double *xi, const double num)
inline double rllt()
inline double lle(double n)
inline double sign(double x, double y)
inline double square(double x)
inline void g1(double a, double b, double &cterm, double &sterm, double &sig)
inline void h12(int mode, int lpivot, int l1, int m, double u[], int iue, double &up, double c[], int ice, int icv, int ncv)
inline int nnls(double a[], int mda, int m, int n, double b[], double x[], double &rnorm, double w[], double zz[], int index[], int &nsetp)
inline double dlle(bool flag)
inline double ille(bool type)
inline void setParticle(double shift, double energy)
inline double rlle(bool type)
inline void set_track(const line &l, bool track = true)
inline bool set_track(const I3ParticleConstPtr &track)
inline bool reconstruct(I3RecoPulseSeriesMapConstPtr pmap, I3FramePtr frame, std::string result)

Public Members

bool type
bool en
bool mpe
unsigned int lcsp
unsigned int unc
double nn
double na
double xc
double yc
double zc
double tc
double dc
double rn
double qtot
double norm
map<mykey, sett> hits
map<mykey, setq> hnhs
double step
double ares
double dres
bool useall
unsigned int nhits
unsigned int hitdoms
unsigned int strings
rrt rt
int rmode
struct MUEX::vtemp::XL xl
I3VectorI3Particle elist
I3Particle muon
struct wform
#include <reader.h>

Public Functions

inline wform()
inline void copy(const wform &e)
inline void destroy()
inline wform(const wform &e)
inline const wform &operator=(const wform &e)
inline ~wform()

Public Members

unsigned long long id
int events
struct XL

Public Functions

inline void ini(bool type)
inline double gfx(double x)
inline double dif(double x)
inline double Ld(double x)
inline double f(double mu, double Q)

Public Members

int nu
double sig
double s2
double b
double w
double W0
double zx
double D
double q
double mx
double x
double y
double z
namespace [anonymous]
namespace I3FRreco

Functions

void pp_start(int h)
void pp_stop()
void pp_setpars(double udf, double thmax, bool invert)
double tfit(double t, double z, double th, double dd, double ddp)
double cfit(double t, double z, double dd, double ddp)
double vector_get(const void *par, const gsl_vector *x, int i)
vtemp &gettmp_inic(const vtemp &par, const gsl_vector *x)
int settmp_inic(gsl_vector *&x, const vtemp &vtmp)
otemp &gettmp_itop(const otemp &par, const gsl_vector *x)
int settmp_itop(gsl_vector *&x, const otemp &otmp)
void reco_e(vtemp&, preco*, double)
void reco_e(otemp&, preco*, double)
void reco_e(multe &e, preco &r, double n)
double illh(const gsl_vector *x, gsl_vector *df, void *par, bool flag)
static const double pepp(const mydom &dom)
double pdens(const mydom &dom, double q)
double tllh(const gsl_vector *x, gsl_vector *df, void *par, bool flag)
double llh_f(const gsl_vector *x, void *par)
void llh_df(const gsl_vector *x, void *par, gsl_vector *df)
void llh_fdf(const gsl_vector *x, void *par, double *f, gsl_vector *df)
double llh_fx(double x, void *par)
void errors(void *par, double val, gsl_vector *x, gsl_vector *xmin, gsl_vector *xmax)
void order(double &x, double &y)
void reconstruct(multe &e, bool do_reco)
void inimctrack(multe &evts, preco &track)
ostream &operator<<(ostream &out, const preco &it)
int indoms()
const map<unsigned long long, mydom> &idoms()
const map<mykey, unsigned long long> &mykeys()
const map<int, int> &mystrings()
static void new_dom(unsigned long long cid)
ostream &operator<<(ostream &out, const event &it)
deque<event*> &ievents()
deque<event*>::iterator get_eitf()
static void everase(deque<event*>::iterator eiti, deque<event*>::iterator eitf)
inline bool evsort(const event *eva, const event *evb)
static void output()
int mdata()
int data()
deque<event*>::iterator next()
deque<multe> &itevents()
ostream &operator<<(ostream &out, const multe &it)
int getinterpo()
void amain(unsigned int rnuma, int intra)
void setgeo(unsigned long long domid, int om, int str, double x, double y, double z, double hv)
int get_itop()
int get_inic()
int get_itsn()
int get_icsn()
unsigned int get_rmsk()
int get_strn()

Variables

static cpdf *ppf[2]
static bool ppstf = false
static int freco
static int vnum
static int vnux
static bool rspe
static bool tfix = false
static bool rfx = true
static double rfx_th = 0
static double rfx_ph = 0
static double rfx_ct = 0
static multe *evts
static double lle
static double dlle
static double nllr
double (*llh)(const gsl_vector*, gsl_vector*, void*, bool)
static bool extf = true
map<unsigned long long, double> qdoms
bool parinic = false
bool reco_reset = true
static int itopreco = 2
static bool itopdok0 = false
static double tAmean
static double tArms
bool paritop = false
struct I3FRreco::sens asens
static int ndoms = 0
static map<unsigned long long, mydom> doms
static map<mykey, unsigned long long> ikeys
static map<int, int> istrings
static unsigned long long ev_t0 = 0
static const long long ONESEC = 10000000000LL
static long long TSTEP = ONESEC
static deque<event*> events
static deque<event*>::iterator eit = events.begin()
static deque<event*>::iterator eiti = eit
static deque<event*>::iterator eitf = eit
static bool nnd = false
unsigned long maxevents = 0
static unsigned long evoutcount
static unsigned long fullevents
static vector<unsigned long long> dreq
static unsigned long xnum = (unsigned long)1.e6
static unsigned long long gt0 = 0
static unsigned long long Ggt = 0
static bool efl = false
static deque<multe> tevents
static unsigned int rnum
static int interpo = 0
namespace MUEX

Functions

ostream &operator<<(ostream &out, const dir &n)
struct MUEX::Ico ico (2)
void check(double &x)
struct MUEX::cpdf pdf (0, 15, 1.0, 35., 1.0, 100)
inline double mb_sqr(double r)
template<int d>
std::ostream &operator<<(std::ostream &os, const Point<d> &p)
void normalize(double n[3])
double flle(double, void*)
double fllh(const gsl_vector*, void*)

Variables

bool compat = false
set<mykey> boms
map<mykey, mydom> doms
static const double dr = 0.16510
static const double pmin = 1e-9
static const double dura = 2.e3
static const double nois = 0.5e-6
static const double umin = nois * dura
static const double efcr = 0.75 * 0.9
static const double zoff = 1948.07
static const double rho = 0.9216
static const double arho = 0.26 * rho
static const double brho = 0.36e-3 * rho
static const double Cf = 2450.08 * (M_PI * dr * dr)
static const double Ff = 5.21 * (0.924 / rho)
struct MUEX::sens asens
struct MUEX::Tilt tilt
struct MUEX::Ice ice
static const double SQRT2PI = M_SQRT2 * M_SQRTPI
class MUEX::Miniball pdf
namespace std

STL namespace.

file I3mue.cxx
#include “mue/I3mue.h
#include “reader.h
#include “llhreco.h
#include “dataclasses/geometry/I3Geometry.h”
#include “dataclasses/status/I3DetectorStatus.h”
#include “dataclasses/physics/I3RecoPulse.h”
#include “dataclasses/physics/I3Particle.h”
#include <sstream>

Functions

I3_MODULE(I3mue)
double res(double x, double mod)
file I3mue.h
#include <icetray/OMKey.h>
#include <icetray/I3ConditionalModule.h>
#include <icetray/I3Frame.h>
#include <icetray/I3Logging.h>
file interpol.cxx
#include <cmath>
#include <cstdlib>
#include “reader.h
#include “interpol.h
#include <fcntl.h>
#include <unistd.h>
#include <iostream>
#include <cstdio>
#include <sys/stat.h>

Defines

FILESAVE

This code is a c++ version of the simplified java class Interpolate of mmc.

file interpol.h
file llhreco.cxx
#include <cmath>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_sf_gamma.h>
#include <gsl/gsl_sf_hyperg.h>
#include <gsl/gsl_multimin.h>
#include <gsl/gsl_roots.h>
#include <fstream>
#include “reader.h
#include “llhreco.h
#include “interpol.h
#include “miniball.h
file /llhreco.cxx
#include “miniball.h
#include “nnls.cxx

Defines

LMAX
LYRS

Functions

Ico(int num)
cpdf(double ita, double sig, double ef, double es, double af, double as)
double lnscat(double d, double t, double Q)
void relax()
void restore()
double gaus(double x, double sig)
double student(double x, double sig)
double addn(double x, double t)
void set(double d)
double pandel(double t)
double cpandel(double t)
double cf(double t)
double dcf(double t, double s)
double dcf(double t, double s, double Q)
double npandel(double t, double Q)

Variables

int vn
int fn
vector<dir> v
vector<vector<int>> f
double ita
double sig
double ef
double es
double af
double as
double xi
double rho
file llhreco.h
file miniball.h
#include <cstdlib>
#include <cassert>
#include <cmath>
#include <iostream>
#include <list>

Functions

inline double mb_sqr(double r)
template<int d>
std::ostream &operator<<(std::ostream &os, const Point<d> &p)
file /miniball.h

Defines

MINIBALL_H_INCLUDED

Functions

inline double mb_sqr(double r)
template<int d>
std::ostream &operator<<(std::ostream &os, const Point<d> &p)
file muex.cxx
#include <icetray/I3ConditionalModule.h>
#include <dataclasses/I3Double.h>
#include <dataclasses/physics/I3MCTree.h>
#include <dataclasses/physics/I3MCTreeUtils.h>
#include <dataclasses/physics/I3RecoPulse.h>
#include <dataclasses/physics/I3Particle.h>
#include <dataclasses/geometry/I3Geometry.h>
#include <dataclasses/status/I3DetectorStatus.h>
#include <dataclasses/calibration/I3Calibration.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_sf_gamma.h>
#include <gsl/gsl_sf_hyperg.h>
#include <gsl/gsl_sf_lambert.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_multimin.h>
#include <fstream>
#include <limits.h>
#include “llhreco.cxx

Typedefs

typedef std::vector<OMKey> OMKeySeries

Functions

I3_POINTER_TYPEDEFS(OMKeySeries)
I3_MODULE(muex)
file nnls.cxx

Defines

c(x)
u(x, y)
b(x)
x(y)
w(x)
zz(x)
index(x)
a(x, y)
c(x)
u(x, y)
b(x)
x(y)
w(x)
zz(x)
index(x)
a(x, y)

Functions

double sign(double x, double y)
double square(double x)
void g1(double a, double b, double &cterm, double &sterm, double &sig)
void h12(int mode, int lpivot, int l1, int m, double u[], int iue, double &up, double c[], int ice, int icv, int ncv)
int nnls(double a[], int mda, int m, int n, double b[], double x[], double &rnorm, double w[], double zz[], int index[], int &nsetp)
file reader.cxx
#include <cstdio>
#include <fcntl.h>
#include <unistd.h>
#include <getopt.h>
#include <cerrno>
#include <cstring>
#include <cstdlib>
#include <sys/stat.h>
#include <sys/mman.h>
#include <sys/types.h>
#include <cmath>
#include <fstream>
#include <algorithm>
#include “reader.h
#include “llhreco.h
#include <gsl/gsl_sf_erf.h>
#include <gsl/gsl_multimin.h>
file reader.h
#include <map>
#include <deque>
#include <vector>
#include <iostream>
#include <cstring>
dir icetray
dir mue
dir mue
dir mue
dir muex
dir private
dir public