LArSamples::OFC Class Reference

#include <OFC.h>

Inheritance diagram for LArSamples::OFC:

Collaboration diagram for LArSamples::OFC:

Public Member Functions
	OFC (const AbsShape &shape, const AbsShape &data, int lwb=-1, int upb=-1, const ShapeErrorData *sed=0, bool withAutoCorr=true)
	Constructor.
	OFC (const AbsShape &shape, const AbsShape &data, const CovMatrix &refErr, int lwb=-1, int upb=-1, const ShapeErrorData *sed=0, bool withAutoCorr=true)
bool	initGValues (const AbsShape &shape, const AbsShape &data, const ShapeErrorData *sed)
bool	initOFCs (const AbsShape &data, const CovMatrix &refErr, bool useCorrs)
bool	initRVectors ()
	OFC (const OFC &other)
virtual	~OFC ()
unsigned int	nSamples () const
int	lwb () const
int	upb () const
double	g (unsigned int i) const
double	gp (unsigned int i) const
const TVectorD &	g () const
const TVectorD &	gp () const
double	G () const
double	Gp () const
double	Gpp () const
double	a (unsigned int i) const
double	b (unsigned int i) const
const TVectorD &	a () const
const TVectorD &	b () const
const TVectorD &	r (unsigned int i) const
unsigned int	rIdx (unsigned int i) const
const TMatrixD &	resProj () const
double	A (const AbsShape &data) const
double	B (const AbsShape &data) const
double	time (const AbsShape &data) const
const CovMatrix &	invGamma () const
const CovMatrix &	Gamma () const
double	dot (const TVectorD &form, const TVectorD &v) const
TVectorD	residual (const TVectorD &v) const
bool	isInRange (int i) const
bool	hasSameRange (int lw, int up) const
bool	hasSameRange (const IndexRange &other) const
bool	hasSameRange (const TVectorD &v) const
bool	providesRange (int lw, int up) const
bool	providesRange (const IndexRange &other) const
bool	providesRange (const TVectorD &v) const
TString	rangeStr () const
int	commonLwb (const IndexRange &other) const
int	commonUpb (const IndexRange &other) const
bool	checkRange (int &l, int &h) const

Private Attributes
TVectorD	m_g
TVectorD	m_gp
TVectorD	m_a
TVectorD	m_b
double	m_G = 0.0
double	m_Gp = 0.0
double	m_Gpp = 0.0
CovMatrix	m_invGamma
CovMatrix	m_Gamma
TMatrixD	m_resProj
std::vector< TVectorD >	m_r
std::vector< unsigned int >	m_rIdx

Detailed Description

Definition at line 27 of file OFC.h.

Constructor & Destructor Documentation

OFC() [1/3]

OFC::OFC	(	const AbsShape &	shape,
		const AbsShape &	data,
		int	lwb = -1,
		int	upb = -1,
		const ShapeErrorData *	sed = 0,
		bool	withAutoCorr = true )

Constructor.

Definition at line 21 of file OFC.cxx.

  : m_g(lw, up), m_gp(lw, up)
{
  ClassCounts::incrementInstanceCount("OFC"); 
  if (!initGValues(shape, data, sed)) return;
  if (!initOFCs(data, CovMatrix(), withAutoCorr)) return;
}

OFC() [2/3]

OFC::OFC	(	const AbsShape &	shape,
		const AbsShape &	data,
		const CovMatrix &	refErr,
		int	lwb = -1,
		int	upb = -1,
		const ShapeErrorData *	sed = 0,
		bool	withAutoCorr = true )

Definition at line 30 of file OFC.cxx.

  : m_g(lw, up), m_gp(lw, up)
{
  ClassCounts::incrementInstanceCount("OFC"); 
  if (!initGValues(shape, data, sed)) return;
  if (!initOFCs(data, refErr, withAutoCorr)) return;
}

OFC() [3/3]

LArSamples::OFC::OFC ( const OFC & other )

inline

Definition at line 41 of file OFC.h.

                          :
    IndexRange(),
        m_g(other.m_g), m_gp(other.m_gp), m_a(other.m_a), m_b(other.m_b),
        m_G(other.m_G), m_Gp(other.m_Gp), m_Gpp(other.m_Gpp),
        m_invGamma(other.m_invGamma) { }

~OFC()

OFC::~OFC ( )

virtual

Definition at line 40 of file OFC.cxx.

{
  ClassCounts::decrementInstanceCount("OFC"); 
}

Member Function Documentation

A()

double OFC::A

(

const AbsShape &

data

)

const

Definition at line 145 of file OFC.cxx.

{
  TVectorD vals(lwb(), upb());
 
  for (int i = lwb(); i <= upb(); i++)
    vals[i] = data.value(i);
 
  return dot(a(), vals);
}

a() [1/2]

const TVectorD & LArSamples::OFC::a ( ) const

inline

Definition at line 67 of file OFC.h.

{ return m_a; }

a() [2/2]

double LArSamples::OFC::a ( unsigned int i ) const

inline

Definition at line 64 of file OFC.h.

{ return m_a[i]; }

B()

double OFC::B

(

const AbsShape &

data

)

const

Definition at line 156 of file OFC.cxx.

{
  TVectorD vals(lwb(), upb());
 
  for (int i = lwb(); i <= upb(); i++)
    vals[i] = data.value(i);
 
  return dot(b(), vals);
}

b() [1/2]

const TVectorD & LArSamples::OFC::b ( ) const

inline

Definition at line 68 of file OFC.h.

{ return m_b; }

b() [2/2]

double LArSamples::OFC::b ( unsigned int i ) const

inline

Definition at line 65 of file OFC.h.

{ return m_b[i]; }

checkRange()

bool IndexRange::checkRange ( int & l, int & h ) const

inherited

Definition at line 14 of file IndexRange.cxx.

{  
  if (l < 0) 
    l = lwb();
  else if (!isInRange(l)) {
    cout << "IndexRange::checkRange : lower bound " << l << " is out of bounds" << endl;
    return false;
  }
 
  if (h < 0) 
    h = upb();
  else if (!isInRange(h)) {
    cout << "IndexRange::checkRange : upper bound " << h << " is out of bounds" << endl;
    return false;
  }
 
  return true;
}

commonLwb()

int LArSamples::IndexRange::commonLwb ( const IndexRange & other ) const

inlineinherited

Definition at line 39 of file IndexRange.h.

{ return (other.lwb() > lwb() ? other.lwb() : lwb()); }

commonUpb()

int LArSamples::IndexRange::commonUpb ( const IndexRange & other ) const

inlineinherited

Definition at line 40 of file IndexRange.h.

{ return (other.upb() < upb() ? other.upb() : upb()); }

dot()

double OFC::dot	(	const TVectorD &	form,
		const TVectorD &	v ) const

Definition at line 173 of file OFC.cxx.

{
  double s = 0;
   for (int i = lwb(); i <= upb(); i++) s += form[i]*v[i];
  return s;
}

G()

double LArSamples::OFC::G ( ) const

inline

Definition at line 60 of file OFC.h.

{ return m_G; }

g() [1/2]

const TVectorD & LArSamples::OFC::g ( ) const

inline

Definition at line 57 of file OFC.h.

{ return m_g; }

g() [2/2]

double LArSamples::OFC::g ( unsigned int i ) const

inline

Definition at line 54 of file OFC.h.

{ return m_g[i]; }

Gamma()

const CovMatrix & LArSamples::OFC::Gamma ( ) const

inline

Definition at line 83 of file OFC.h.

{ return m_Gamma; }

Gp()

double LArSamples::OFC::Gp ( ) const

inline

Definition at line 61 of file OFC.h.

{ return m_Gp; }

gp() [1/2]

const TVectorD & LArSamples::OFC::gp ( ) const

inline

Definition at line 58 of file OFC.h.

{ return m_gp; }

gp() [2/2]

double LArSamples::OFC::gp ( unsigned int i ) const

inline

Definition at line 55 of file OFC.h.

{ return m_gp[i]; }

Gpp()

double LArSamples::OFC::Gpp ( ) const

inline

Definition at line 62 of file OFC.h.

{ return m_Gpp; }

hasSameRange() [1/3]

bool LArSamples::IndexRange::hasSameRange ( const IndexRange & other ) const

inlineinherited

Definition at line 30 of file IndexRange.h.

{ return hasSameRange(other.lwb(), other.upb()); }

hasSameRange() [2/3]

bool LArSamples::IndexRange::hasSameRange ( const TVectorD & v ) const

inlineinherited

Definition at line 31 of file IndexRange.h.

{ return hasSameRange(v.GetLwb(), v.GetUpb()); }

hasSameRange() [3/3]

bool LArSamples::IndexRange::hasSameRange ( int lw, int up ) const

inlineinherited

Definition at line 29 of file IndexRange.h.

{ return (lwb() == lw && upb() == up); }

initGValues()

bool OFC::initGValues	(	const AbsShape &	shape,
		const AbsShape &	data,
		const ShapeErrorData *	sed )

Definition at line 46 of file OFC.cxx.

{
  CovMatrix errors;
  if (!shape.interpolate(data, m_g, errors, lwb(), upb())) {
    cout <<"ERROR in OFC::initGValues : indices were invalid!" << endl;
    m_g.ResizeTo(0);
    return false;
  }  
  if (!shape.interpolateDiff(data, m_gp, lwb(), upb())) { // lwb and upb come from g, so set above correctly...
    cout <<"ERROR in OFC::initGValues : (diff) indices were invalid!" << endl;
    m_g.ResizeTo(0);
    m_gp.ResizeTo(0);
    return false;
  }
  m_g.ResizeTo(m_gp.GetLwb(), m_gp.GetUpb()); // in case gp was more restricted
  if (sed) {
    TVectorD xi = sed->xi(m_gp.GetLwb(), m_gp.GetUpb());
    TVectorD xip = sed->xip(m_gp.GetLwb(), m_gp.GetUpb());
    m_g += xi;
    m_gp += xip;
  }
  return true;
}

initOFCs()

bool OFC::initOFCs	(	const AbsShape &	data,
		const CovMatrix &	refErr,
		bool	useCorrs )

Definition at line 71 of file OFC.cxx.

{
  m_a.ResizeTo(lwb(), upb());
  m_b.ResizeTo(lwb(), upb());
 
  //cout << nSamples() << endl;
  //data.covarianceMatrix().Print();
  
  CovMatrix noiseCov = data.covarianceMatrix(lwb(), upb(), refErr, withAutoCorr);
  CovMatrix noiseInvCov = data.invCovarianceMatrix(lwb(), upb(), refErr, withAutoCorr);
  //noiseCov.Print();
  
  TVectorD gam = m_g;
  TVectorD gamp = m_gp;
  gam *= noiseInvCov;
  gamp *= noiseInvCov;
 
  m_G   = m_g*gam;
  m_Gp  = m_g*gamp;
  m_Gpp = m_gp*gamp;
 
  double denom = G()*Gpp() - Gp()*Gp();
  //cout << m_G   << " " << m_g*gam << endl;
  //cout << m_Gp  << " " << m_g*gamp << " " << m_gp*gam << endl;
  //cout << m_Gpp << " " << m_gp*gamp << endl;
 
  for (int i = lwb(); i <= upb(); i++) {
    m_a(i) = (Gpp()*gam(i) - Gp()*gamp(i))/denom;
    m_b(i) = (Gp()*gam(i) - G()*gamp(i))/denom;
  }
 
  m_resProj.ResizeTo(lwb(), upb(), lwb(), upb());
  for (int i = lwb(); i <= upb(); i++) {
    for (int j = lwb(); j <= upb(); j++) {
      m_resProj(i,j) = (i == j ? 1 : 0) - g(i)*a(j) + gp(i)*b(j);
    }
  }
 
  //m_a.Print();
  //cout << m_a*m_g << " " << m_a*m_gp << " " << m_b*m_g << " "<< m_b*m_gp << endl;
  
  CovMatrix resNorm(lwb(), upb());
  for (int i = lwb(); i <= upb(); i++) {
    for (int j = lwb(); j <= upb(); j++) {
      double rn = noiseCov(i, j) - (Gpp()*g(i)*g(j) - Gp()*(g(i)*gp(j) + g(j)*gp(i)) + G()*gp(i)*gp(j))/denom;
      resNorm(i, j) = rn;
    }
  }
  
  m_Gamma.ResizeTo(lwb(), upb(), lwb(), upb());
  m_Gamma = noiseCov;
  m_Gamma.Similarity(m_resProj);
 
  m_invGamma.ResizeTo(lwb(), upb(), lwb(), upb());
  m_invGamma = noiseInvCov;  
  m_invGamma.Similarity(m_resProj);
  return true;
}

initRVectors()

bool OFC::initRVectors

(

)

Definition at line 131 of file OFC.cxx.

{
  m_r.clear();
  m_rIdx.clear();
  for (unsigned int k = 0; k < nSamples() - 2; k++) {
    TVectorD res(nSamples());
    res[k] = 1;
    m_r.push_back(residual(res));
    m_rIdx.push_back(k);
  }
  return true;
}

invGamma()

const CovMatrix & LArSamples::OFC::invGamma ( ) const

inline

Definition at line 82 of file OFC.h.

{ return m_invGamma; }

isInRange()

bool LArSamples::IndexRange::isInRange ( int i ) const

inlineinherited

Definition at line 27 of file IndexRange.h.

{ return (i >= lwb() && i <= upb()); }

lwb()

int LArSamples::OFC::lwb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 51 of file OFC.h.

{ return g().GetLwb(); }

nSamples()

unsigned int LArSamples::OFC::nSamples ( ) const

inline

Definition at line 49 of file OFC.h.

{ return m_g.GetNoElements(); }

providesRange() [1/3]

bool LArSamples::IndexRange::providesRange ( const IndexRange & other ) const

inlineinherited

Definition at line 34 of file IndexRange.h.

{ return providesRange(other.lwb(), other.upb()); }

providesRange() [2/3]

bool LArSamples::IndexRange::providesRange ( const TVectorD & v ) const

inlineinherited

Definition at line 35 of file IndexRange.h.

{ return providesRange(v.GetLwb(), v.GetUpb()); }

providesRange() [3/3]

bool LArSamples::IndexRange::providesRange ( int lw, int up ) const

inlineinherited

Definition at line 33 of file IndexRange.h.

{ return (lwb() <= lw && upb() >= up); }

r()

const TVectorD & LArSamples::OFC::r ( unsigned int i ) const

inline

Definition at line 70 of file OFC.h.

{ return m_r[i]; }

rangeStr()

TString LArSamples::IndexRange::rangeStr ( ) const

inlineinherited

Definition at line 37 of file IndexRange.h.

{ return Form("[%d, %d]", lwb(), upb()); }

residual()

TVectorD OFC::residual

(

const TVectorD &

v

)

const

Definition at line 167 of file OFC.cxx.

{
  return v + (-1)*dot(a(), v)*g() + (-1)*dot(b(), v)*gp();
}

resProj()

const TMatrixD & LArSamples::OFC::resProj ( ) const

inline

Definition at line 73 of file OFC.h.

{ return m_resProj; }

rIdx()

unsigned int LArSamples::OFC::rIdx ( unsigned int i ) const

inline

Definition at line 71 of file OFC.h.

{ return m_rIdx[i]; }

time()

double LArSamples::OFC::time ( const AbsShape & data ) const

inline

Definition at line 77 of file OFC.h.

                                            {
      double den = A(data);
      return den == 0 ? 0 : B(data)/den;
    }

upb()

int LArSamples::OFC::upb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 52 of file OFC.h.

{ return g().GetUpb(); }

Member Data Documentation

m_a

TVectorD LArSamples::OFC::m_a

private

Definition at line 90 of file OFC.h.

m_b

TVectorD LArSamples::OFC::m_b

private

Definition at line 90 of file OFC.h.

m_G

double LArSamples::OFC::m_G = 0.0

private

Definition at line 91 of file OFC.h.

m_g

TVectorD LArSamples::OFC::m_g

private

Definition at line 90 of file OFC.h.

m_Gamma

CovMatrix LArSamples::OFC::m_Gamma

private

Definition at line 92 of file OFC.h.

m_Gp

double LArSamples::OFC::m_Gp = 0.0

private

Definition at line 91 of file OFC.h.

m_gp

TVectorD LArSamples::OFC::m_gp

private

Definition at line 90 of file OFC.h.

m_Gpp

double LArSamples::OFC::m_Gpp = 0.0

private

Definition at line 91 of file OFC.h.

m_invGamma

CovMatrix LArSamples::OFC::m_invGamma

private

Definition at line 92 of file OFC.h.

m_r

std::vector<TVectorD> LArSamples::OFC::m_r

private

Definition at line 94 of file OFC.h.

m_resProj

TMatrixD LArSamples::OFC::m_resProj

private

Definition at line 93 of file OFC.h.

m_rIdx

std::vector<unsigned int> LArSamples::OFC::m_rIdx

private

Definition at line 95 of file OFC.h.

---

The documentation for this class was generated from the following files:

• OFC.h
• OFC.cxx