LArSamples::History Class Reference

#include <History.h>

Collaboration diagram for LArSamples::History:

Public Member Functions
	History (const HistoryContainer &container, const std::vector< const EventData * > &eventData, unsigned int hash, const AbsShapeErrorGetter *shapeErrorGetter=0)
	Constructor.
	History (const std::vector< const Data * > &data, const CellInfo &info, const std::vector< const EventData * > &eventData, unsigned int hash, const AbsShapeErrorGetter *shapeErrorGetter=0)
virtual	~History ()
HistoryContainer *	dissolve ()
unsigned int	nData () const
const Data *	data (unsigned int i) const
const Data *	data_for_event (int event, int run=-1) const
const Data *	data_for_event (const EventData &eventData) const
const CellInfo *	cellInfo () const
const std::vector< const EventData * > &	eventData () const
bool	sum (SimpleShape *&sum, SimpleShape *&reference) const
bool	allShape (GraphShape *&allData, SimpleShape *&allRef) const
double	allChi2 (Chi2Params pars) const
double	chi2 (int i, int lwb=-1, int upb=-1, int chi2Params=DefaultChi2, ShapeErrorType shapeErrorType=BestShapeError, unsigned int *nDof=0) const
double	chi2_k (int i, double k, int lwb=-1, int upb=-1, int chi2Params=DefaultChi2) const
double	maxChi2 (int lwb=-1, int upb=-1, int chi2Params=DefaultChi2) const
OFC *	ofc (unsigned int i, int lwb=-1, int upb=-1, double time=Definitions::none, bool useCorrs=true) const
History *	refit (Chi2Params pars=DefaultChi2) const
History *	adjust () const
History *	filter (const TString &cuts) const
bool	refVal (unsigned int i, unsigned int sample, double &val, double &err) const
bool	delta (unsigned int i, unsigned int sample, double &del) const
TVectorD	deltas (unsigned int i, int lwb=-1, int upb=-1, bool correct=true) const
bool	isValid () const
TString	description (unsigned int verbosity=1) const
SimpleShape *	referenceShape (unsigned int k, double adcMax=-1, double time=Definitions::none, bool samplesOnly=false) const
SimpleShape *	deltaShape (unsigned int k, int lwb=-1, int upb=-1) const
unsigned int	hash () const
void	setShapeErrorGetter (const AbsShapeErrorGetter *err) const
const AbsShapeErrorGetter *	shapeErrorGetter () const
const ShapeErrorData *	shapeErrorData (CaloGain::CaloGain gain, ShapeErrorType shapeErrorType=BestShapeError, const Residual *res=0) const
const ScaledErrorData *	scaledErrorData (unsigned int i, double adcMax=-1, double time=Definitions::none, ShapeErrorType shapeErrorType=BestShapeError) const
bool	residualOffset (unsigned int i, short sample, double &offset, double adcMax=-1, double time=Definitions::none) const
bool	residualError (unsigned int i, short sample1, short sample2, double &offset, double adcMax=-1, double time=Definitions::none) const
bool	drawWithReference (int i, const TString &atlasTitle="") const
bool	drawSumWithReference () const
bool	drawAllWithReference (bool refit=false) const
bool	drawResiduals (int k=-1, bool errors=true, bool rescale=true) const
Averager *	calculatePedestal (int i) const
Residual *	residual (unsigned int k, bool correct=true, bool zeroTime=false) const
Residuals *	residuals (CaloGain::CaloGain gain=CaloGain::LARNGAIN, double absResTrunc=-1, bool correct=true, bool zeroTime=false) const
double	upstreamEnergy (unsigned int k) const
double	chi2Anomaly (double chi2Cut, unsigned int nDof=3) const
void	setInterface (const Interface *interface) const

Private Attributes
std::vector< const Data * >	m_data
CellInfo	m_cellInfo
std::vector< const EventData * >	m_eventData
unsigned int	m_hash
const AbsShapeErrorGetter *	m_shapeErrorGetter
const Interface *	m_interface = nullptr

Detailed Description

Definition at line 35 of file History.h.

Constructor & Destructor Documentation

History() [1/2]

LArSamples::History::History	(	const HistoryContainer &	container,
		const std::vector< const EventData * > &	eventData,
		unsigned int	hash,
		const AbsShapeErrorGetter *	shapeErrorGetter = 0 )

Constructor.

History() [2/2]

History::History	(	const std::vector< const Data * > &	data,
		const CellInfo &	info,
		const std::vector< const EventData * > &	eventData,
		unsigned int	hash,
		const AbsShapeErrorGetter *	shapeErrorGetter = 0 )

Definition at line 56 of file History.cxx.

  : m_data(data), m_cellInfo(info), m_eventData(eventData), m_hash(hash), 
    m_shapeErrorGetter(shapeErrorGetter) 
{
  ClassCounts::incrementInstanceCount("History");
  unsigned int i = 0;
  for (const Data* pdata : m_data)
    pdata->setCallBacks(this, i);
}

~History()

History::~History ( )

virtual

Definition at line 69 of file History.cxx.

{
  ClassCounts::decrementInstanceCount("History");
  for (const Data* data : m_data)
    delete data;
  for (const EventData* eventData : m_eventData)
    delete eventData;  
}

Member Function Documentation

adjust()

History * History::adjust

(

)

const

Definition at line 259 of file History.cxx.

{
  std::vector<const Data*> datas;
  DataTweaker tw;
  tw.setAdjust(true);
 
  for (unsigned int j = 0; j < nData(); j++) {
    Data* newData = tw.tweak(*data(j));
    if (!newData) {
      for (unsigned int k = 0; k < datas.size(); k++) delete datas[k];
      return nullptr;
    }
    datas.push_back(newData);
  }
 
  std::vector<const EventData*> eventData;
  for (const EventData* event : m_eventData)
    eventData.push_back(new EventData(*event));
 
  return new History(datas, *cellInfo(), eventData, hash(), shapeErrorGetter());
}

allChi2()

double History::allChi2

(

Chi2Params

pars

)

const

Definition at line 439 of file History.cxx.

{
  Chi2Calc c2c(pars);
  
  GraphShape* allData = nullptr;
  SimpleShape* allRef = nullptr;
  if (!allShape(allData, allRef)) return -1;
  double chi2Value = c2c.chi2(*allData, *allRef);
  delete allData;
  delete allRef;
  return chi2Value;
}

allShape()

bool History::allShape	(	GraphShape *&	allData,
		SimpleShape *&	allRef ) const

Definition at line 413 of file History.cxx.

{
  allData = nullptr;
  allRef  = nullptr;
  
  for (unsigned int j = 0; j < nData(); j++) {
   if ( m_data[j]->isDisconnected()) continue; 
   if (m_data[j]->adcMax() < 1) continue;
    GraphShape* thisData = new GraphShape(*m_data[j], 1/m_data[j]->adcMax(), -m_data[j]->ofcTime());
    if (allData) {
      allData->add(*thisData);
      delete thisData;
    }
    else allData = thisData;
    SimpleShape* thisRef = referenceShape(j, 1); // normalized to 1, like the data
    if (allRef) {
      allRef->add(*thisRef);
      delete thisRef;
    }
    else allRef = thisRef;
  }
 
  return true;
}

calculatePedestal()

Averager * History::calculatePedestal

(

int

i

)

const

Definition at line 584 of file History.cxx.

{
  Averager* avg = new Averager(1);
  for (unsigned int k = 0; k < nData(); k++) {
    if (i >= 0 && i != (int)k) continue;
    SimpleShape* reference = referenceShape(k);
    if (!reference) continue;  
    double v,e;
    for (unsigned int l = 0; l < data(k)->nPoints(); l++)
      if (reference->interpolate(data(k)->time(l), v, e) == -1) { // we're before the ref shapeError
        TVectorD ped(1);
        ped(0) = data(k)->value(l);
        avg->fill(ped);
      }
    delete reference;
  }
  return avg;
}

cellInfo()

const CellInfo * LArSamples::History::cellInfo ( ) const

inline

Definition at line 56 of file History.h.

{ return &m_cellInfo; }

chi2()

double History::chi2	(	int	i,
		int	lwb = -1,
		int	upb = -1,
		int	chi2Params = DefaultChi2,
		ShapeErrorType	shapeErrorType = BestShapeError,
		unsigned int *	nDof = 0 ) const

Definition at line 164 of file History.cxx.

{
  SimpleShape* reference = referenceShape(i);
  if (!reference) return -1;
  Chi2Calc c2c(chi2Params);
  if ( m_data[i]->isDisconnected()) return -1;
  const ScaledErrorData* sea = scaledErrorData(i, -1, Definitions::none, shapeErrorType);
  if (!sea && shapeErrorType != NoShapeError && shapeErrorType != BestShapeError) return -1;
  double chi2Value = c2c.chi2(*m_data[i], *reference, sea, lwb, upb);
  if (nDof) *nDof = c2c.nDof();
  if (sea) delete sea;
  delete reference;
  return chi2Value;
}

chi2_k()

double History::chi2_k	(	int	i,
		double	k,
		int	lwb = -1,
		int	upb = -1,
		int	chi2Params = DefaultChi2 ) const

Definition at line 180 of file History.cxx.

{
  const AbsShapeErrorGetter* oldGetter = shapeErrorGetter();
  UniformShapeErrorGetter* uniGetter = new UniformShapeErrorGetter(k);
  CombinedShapeErrorGetter* combGetter = new CombinedShapeErrorGetter();
  if (oldGetter) combGetter->add(*oldGetter);
  combGetter->add(*uniGetter);
  setShapeErrorGetter(combGetter);
  double chi2Value = chi2(i, lwb, upb, chi2Params);
  setShapeErrorGetter(oldGetter);
  delete combGetter;
  delete uniGetter;
  return chi2Value;
}

chi2Anomaly()

double History::chi2Anomaly	(	double	chi2Cut,
		unsigned int	nDof = 3 ) const

Definition at line 658 of file History.cxx.

{
  if (nData() == 0 || chi2Cut <= 0) return -1;
  double nBadChi2 = 0;
  for (unsigned int k = 0; k < nData(); k++) if (chi2(k) > chi2Cut) nBadChi2++;
  double effBad = nBadChi2/nData();
  double effRef = TMath::Prob(chi2Cut, nDof);
  double dEff = sqrt(effRef*(1 - effRef))/sqrt(nData());
  if (dEff == 0) dEff = 1/sqrt(nData());
  return (effBad - effRef)/dEff;
}

data()

const Data * History::data

(

unsigned int

i

)

const

Definition at line 91 of file History.cxx.

{ 
  if (i >= nData()) return nullptr;
  return m_data[i]; 
}

data_for_event() [1/2]

const Data * History::data_for_event

(

const EventData &

eventData

)

const

Definition at line 109 of file History.cxx.

{
  // turn off caching
 
  //if (m_dataForEvent.size() == 0) 
  //  for (unsigned int k = 0; k < nData(); k++)
  //    m_dataForEvent[data(k)->event()] = std::make_pair(data(k)->run(), data(k));
 
  //std::unordered_map<int, std::pair<int, const Data*> >::const_iterator
  //  data = m_dataForEvent.find(eventData.event());
  //if (data == m_dataForEvent.end()) return 0;
  //if (data->second.first != eventData.run()) return data_for_event(eventData.event(), eventData.run()); // wrong run: revert to pedestrian method...
  //return data->second.second;
 
  return data_for_event(eventData.event(), eventData.run());
}

data_for_event() [2/2]

const Data * History::data_for_event	(	int	event,
		int	run = -1 ) const

Definition at line 98 of file History.cxx.

{
  for (unsigned int k = 0; k < nData(); k++) {
    if (data(k)->event() != event) continue;
    if (run > 0 && data(k)->run() != run) continue;
    return data(k);
  }
  return nullptr;
}

delta()

bool History::delta	(	unsigned int	i,
		unsigned int	sample,
		double &	del ) const

Definition at line 366 of file History.cxx.

{
  if (k >= nData()) return 0;
  TVectorD dv = deltas(k, sample, sample);
  if (dv.GetNrows() == 0) return false;
  del = dv[sample];
  return true;
}

deltas()

TVectorD History::deltas	(	unsigned int	i,
		int	lwb = -1,
		int	upb = -1,
		bool	correct = true ) const

Definition at line 376 of file History.cxx.

{
  if (k >= nData()) return TVectorD();
  SimpleShape* reference = referenceShape(k);
  if (!reference) return TVectorD();  
  Chi2Calc c2c;
  CovMatrix errors;
  const ScaledErrorData* sea = (correct ? scaledErrorData(k) : nullptr);
  TVectorD dv = c2c.deltas(*data(k), *reference, errors, sea, lwb, upb);
  if (sea) delete sea;
  delete reference;
  return dv;
}

deltaShape()

SimpleShape * History::deltaShape	(	unsigned int	k,
		int	lwb = -1,
		int	upb = -1 ) const

Definition at line 556 of file History.cxx.

{
  if (k >= nData()) return nullptr;
  SimpleShape* reference = referenceShape(k);
  if (!reference) return nullptr;  
  Chi2Calc c2c;
  CovMatrix errors;
  const ScaledErrorData* sea = scaledErrorData(k);
  TVectorD dv = c2c.deltas(*data(k), *reference, errors, sea, lwb, upb);
  if (sea) delete sea;
  delete reference;
  if (dv.GetNrows() == 0) return nullptr;
  SimpleShape* shape = new SimpleShape(dv.GetNrows(), Definitions::samplingInterval, data(k)->time(c2c.lwb()) + data(k)->ofcTime());
  for (int l = c2c.lwb(); l <= c2c.upb(); l++) shape->set(l - c2c.lwb(), dv(l), TMath::Sqrt(errors(l, l)));
  return shape;
}

description()

TString History::description

(

unsigned int

verbosity = 1

)

const

Definition at line 574 of file History.cxx.

{
  TString desc = "";
  for (unsigned int i = 0; i < nData(); i++)
    desc += Form("  #%-2d : ", i) + data(i)->description(verbosity) + "\n";
  if (desc == "") return desc;
  return cellInfo()->location(2) + "\n" + desc;
}

dissolve()

HistoryContainer * History::dissolve

(

)

Definition at line 79 of file History.cxx.

{
  HistoryContainer* histCont = new HistoryContainer(new CellInfo(*cellInfo()));
  for (unsigned int k = 0; k < nData(); k++) {
    Data* newData = new Data(*data(k));
    histCont->add(newData->dissolve());
  }
  delete this;
  return histCont;
}

drawAllWithReference()

bool History::drawAllWithReference

(

bool

refit = false

)

const

Definition at line 484 of file History.cxx.

{  
  if (!nData()) {
    cout << "No data" << endl;
    return false;
  }
  
  if (doRefit) {
    History* refitted_history = refit(DefaultChi2);
    if (!refitted_history) return false;
    bool result = refitted_history->drawAllWithReference(false);
    delete refitted_history;
    return result;
  }
  
  // Use the shape of the gain of the first pulse...
  SimpleShape* refShape = referenceShape(0, 1000, 0);
  SimpleShape* smpShape = referenceShape(0, 1000, 0, true);
  std::vector<const AbsShape*> shapes;
  
  for (unsigned int i = 0; i < nData(); i++) {
    if (m_data[i]->adcMax() < 1) continue;
    SimpleShape* shape = new SimpleShape(*m_data[i], 1000/m_data[i]->adcMax(), -m_data[i]->ofcTime());
    const ScaledErrorData* sed = scaledErrorData(i, 1000);
    if (sed) {
      for (unsigned int k = 0; k < shape->nPoints(); k++)
        if (sed->isInRange(k))
          shape->set(k, shape->value(k) - sed->offsets()(k));
      delete sed;
    }
    shapes.push_back(shape);
  }
  
  ShapeDrawer drawer(Legend);
  return drawer.drawAndDelete(Form("%s (Normalized Shape, max at 1000)", cellInfo()->location(2).Data()), shapes, refShape, smpShape);
}

drawResiduals()

bool History::drawResiduals	(	int	k = -1,
		bool	errors = true,
		bool	rescale = true ) const

Definition at line 522 of file History.cxx.

{  
  std::vector<const AbsShape*> shapes;
  for (unsigned int i = 0; i < nData(); i++) {
    if (m_data[i]->adcMax() < 1) continue;
    if (k >= 0 && k != (int)i) continue;
    SimpleShape* shape = deltaShape(i);
    if (!errors) 
      for (unsigned int idx = 0; idx < shape->nPoints(); idx++) shape->setError(idx, 0);
    if (!shape) continue;
    if (rescale) {
      SimpleShape* scaled = new SimpleShape(*shape, 1/m_data[i]->adcMax(), 9.0*i/nData() - m_data[i]->ofcTime());
      delete shape;
      shape = scaled;
    }
    shapes.push_back(shape);
  }
  ShapeDrawer drawer(DataFirst);
  TString title = (rescale ? "Normalized " : "") + TString("residuals for %s");
  return drawer.drawAndDelete(Form(title.Data(), cellInfo()->location(2).Data()), shapes);  
}

drawSumWithReference()

bool History::drawSumWithReference

(

)

const

Definition at line 474 of file History.cxx.

{
  SimpleShape* dataShape, *refShape;
  if (!sum(dataShape, refShape)) return false;
    
  ShapeDrawer drawer(DataFirst | Legend);
  return drawer.drawAndDelete("", dataShape, refShape);
}

drawWithReference()

bool History::drawWithReference	(	int	i,
		const TString &	atlasTitle = "" ) const

Definition at line 453 of file History.cxx.

{
  if ((unsigned int)k >= nData()) return false;
  auto refShape = std::unique_ptr<SimpleShape>(referenceShape(k));
  auto smpShape = std::unique_ptr<SimpleShape>(referenceShape(k, -1, Definitions::none, true));
  
  if (!refShape || !smpShape) return false;
  int pars = DataFirst | Legend;
  TString title = "";
  if (atlasTitle != "") { // Make sure AtlasStyle is set in ROOT, otherwise will not be as pretty...
    pars |= AtlasStyle;
    title = atlasTitle;
  }
  else 
    title = Form("%s, run %d, event %d", cellInfo()->location(1).Data(), m_data[k]->run(), m_data[k]->event());
  ShapeDrawer drawer(pars);
  bool result = drawer.draw(title, m_data[k], refShape.get(), smpShape.get());
  return result;
}

eventData()

const std::vector< const EventData * > & LArSamples::History::eventData ( ) const

inline

Definition at line 58 of file History.h.

{ return m_eventData; }

filter()

History * History::filter

(

const TString &

cuts

)

const

Definition at line 282 of file History.cxx.

{
  FilterParams f;
  if (!f.set(cuts)) return nullptr;
  
  std::vector<const Data*> datas;
 
  for (unsigned int j = 0; j < nData(); j++) {
    if (!f.pass(hash(), *this, j)) continue;
    datas.push_back(new Data(*data(j)));
  }
 
  std::vector<const EventData*> eventData;
  for (const EventData* event : m_eventData)
    eventData.push_back(new EventData(*event));
 
  return new History(datas, *cellInfo(), eventData, hash(), shapeErrorGetter());
}

hash()

unsigned int LArSamples::History::hash ( ) const

inline

Definition at line 86 of file History.h.

{ return m_hash; }

isValid()

bool History::isValid

(

)

const

Definition at line 152 of file History.cxx.

{
  if (!m_cellInfo.isValid()) return false;
  if (nData() == 0) return false;
  
  for (const Data* data : m_data)
    if (!data->isValid()) return false;
  
  return true;
}

maxChi2()

double History::maxChi2	(	int	lwb = -1,
		int	upb = -1,
		int	chi2Params = DefaultChi2 ) const

Definition at line 196 of file History.cxx.

{
  double maxChi2 = -2;  
  for (unsigned int i = 0; i < nData(); i++) {
    double chi2Value = chi2(i, lwb, upb, chi2Params);
    if (chi2Value > maxChi2) maxChi2 = chi2Value;
  }
  return maxChi2;
}

nData()

unsigned int LArSamples::History::nData ( ) const

inline

Definition at line 51 of file History.h.

{ return m_data.size(); }

ofc()

OFC * History::ofc	(	unsigned int	i,
		int	lwb = -1,
		int	upb = -1,
		double	time = Definitions::none,
		bool	useCorrs = true ) const

Definition at line 207 of file History.cxx.

{
  if (k >= nData()) return nullptr;
  if (Definitions::isNone(time)) {
    time = data(k)->ofcTime();
    //cout << "Using reference time = " << time << endl;
  }
  SimpleShape* reference = referenceShape(k, 1, time); // ADC=1 : we need a normalized shape
  if (!reference) return nullptr;
  OFC* result = new OFC(*reference, *m_data[k], lwb, upb, shapeErrorData(CaloGain::LARHIGHGAIN), withAutoCorr); // FixMe
  delete reference;
  if (!result) return nullptr;
  if (result->g().GetNrows() == 0) { delete result; result = nullptr; }
  return result;
}

referenceShape()

SimpleShape * History::referenceShape	(	unsigned int	k,
		double	adcMax = -1,
		double	time = Definitions::none,
		bool	samplesOnly = false ) const

Definition at line 545 of file History.cxx.

{
  if (!cellInfo()->shape(m_data[k]->gain())) return nullptr;
  if (adcMax < 0) adcMax = m_data[k]->adcMax();
  if (Definitions::isNone(time)) time = m_data[k]->ofcTime();
  SimpleShape* shape = new SimpleShape(*cellInfo()->shape(m_data[k]->gain()), adcMax, time, samplesOnly);
  return shape;
}

refit()

History * History::refit

(

Chi2Params

pars = DefaultChi2

)

const

Definition at line 235 of file History.cxx.

{
  std::vector<const Data*> datas;
  DataTweaker tw;
  tw.setRefit(true);
  tw.setFitParams(pars);
 
  for (unsigned int j = 0; j < nData(); j++) {
    Data* refitData = tw.tweak(*data(j));
    if (!refitData) {
      for (unsigned int k = 0; k < datas.size(); k++) delete datas[k];
      return nullptr;
    }
    datas.push_back(refitData);
  }
 
  std::vector<const EventData*> eventData;
  for (const EventData* event : m_eventData)
    eventData.push_back(new EventData(*event));
 
  return new History(datas, *cellInfo(), eventData, hash(), shapeErrorGetter());
}

refVal()

bool History::refVal	(	unsigned int	i,
		unsigned int	sample,
		double &	val,
		double &	err ) const

Definition at line 224 of file History.cxx.

{
  if (k >= nData()) return false;
  SimpleShape* reference = referenceShape(k);
  if (!reference) return false;  
  if (reference->interpolate(m_data[k]->time(sample), val, err) != 0) { delete reference; return false; }
  delete reference;
  return true;
}

residual()

Residual * History::residual	(	unsigned int	k,
		bool	correct = true,
		bool	zeroTime = false ) const

Definition at line 604 of file History.cxx.

{
  if (k >= nData()) return nullptr;
  TVectorD del = deltas(k, -1, -1, correct);
  return new Residual(del, data(k)->run(), data(k)->event(), data(k)->adcMax(), 
                      (zeroTime ? 0 : data(k)->ofcTime() /*- Definitions::samplingTime(del.GetLwb()) */));
}

residualError()

bool History::residualError	(	unsigned int	i,
		short	sample1,
		short	sample2,
		double &	offset,
		double	adcMax = -1,
		double	time = Definitions::none ) const

Definition at line 402 of file History.cxx.

{
  const ScaledErrorData* sea = scaledErrorData(k, adcMax, time);
  if (!sea) return false;
  if (!sea->isInRange(sample1) || !sea->isInRange(sample2)) { delete sea; return false; }
  error = sea->errors()(sample1, sample2);
  delete sea;
  return true;
}

residualOffset()

bool History::residualOffset	(	unsigned int	i,
		short	sample,
		double &	offset,
		double	adcMax = -1,
		double	time = Definitions::none ) const

Definition at line 391 of file History.cxx.

{
  const ScaledErrorData* sea = scaledErrorData(k, adcMax, time);
  if (!sea) return false;
  if (!sea->isInRange(sample)) { delete sea; return false; }
  offset = sea->offsets()(sample);
  delete sea;
  return true;
}

residuals()

Residuals * History::residuals	(	CaloGain::CaloGain	gain = CaloGain::LARNGAIN,
		double	absResTrunc = -1,
		bool	correct = true,
		bool	zeroTime = false ) const

Definition at line 613 of file History.cxx.

{
  Chi2Calc c2c;
  CovMatrix dummyErrors;
  auto residuals = std::make_unique<Residuals>();
  for (unsigned int k = 0; k < nData(); k++) {
    if (gain != CaloGain::LARNGAIN && data(k)->gain() != gain) continue;
    auto res = std::unique_ptr<Residual>(residual(k, correct, zeroTime));
    if (!res) { cout << "Error calculating residual for hash = " << m_hash << ", index = " << k << endl; return nullptr; }
    if (residuals->size() > 0 && !residuals->hasSameRange(*res)) {
      cout << "Warning for hash = " << m_hash << ", index = " << k << " : index interval changed from [" 
           << residuals->lwb() << ", " << residuals->upb() << "] to " << res->rangeStr() << endl;      
      return nullptr;
    }
    bool pass = true;
    for (int i = res->lwb(); i < res->upb(); i++) {
      if (absResCut > 0 && TMath::Abs(res->scaledDelta(i)) > absResCut) {
        pass = false;
        break;
      }
    }
    if (pass) residuals->add(*res);
  }
  return residuals.release();
}

scaledErrorData()

const ScaledErrorData * History::scaledErrorData	(	unsigned int	i,
		double	adcMax = -1,
		double	time = Definitions::none,
		ShapeErrorType	shapeErrorType = BestShapeError ) const

Definition at line 349 of file History.cxx.

{
  if (shapeErrorType == NoShapeError || !shapeErrorGetter()) return nullptr;
  if (k >= nData() || data(k)->adcMax() <= 0) return nullptr;
 
  const ShapeErrorData* sed = shapeErrorData(data(k)->gain(), shapeErrorType, residual(k, false));
  if (!sed) return nullptr;
 
  double sf = (adcMax < 0 ? data(k)->adcMax() : adcMax);
  double ts = (Definitions::isNone(time) ? data(k)->ofcTime() : time);  
  ScaledErrorData* sced = new ScaledErrorData(*sed, sf, ts);
  delete sed;
  return sced;
}

setInterface()

void LArSamples::History::setInterface ( const Interface * interface ) const

inline

Definition at line 107 of file History.h.

{ m_interface = interface; }

setShapeErrorGetter()

void LArSamples::History::setShapeErrorGetter ( const AbsShapeErrorGetter * err ) const

inline

Definition at line 87 of file History.h.

{ m_shapeErrorGetter = err; }

shapeErrorData()

const ShapeErrorData * History::shapeErrorData	(	CaloGain::CaloGain	gain,
		ShapeErrorType	shapeErrorType = BestShapeError,
		const Residual *	res = 0 ) const

Definition at line 302 of file History.cxx.

{
  if (shapeErrorType == NoShapeError || !shapeErrorGetter()) return nullptr;
  if (shapeErrorType == BestShapeError) {
    for (unsigned int i = 0; i < NShapeErrorTypes; i++) {
      const ShapeErrorData* sed = shapeErrorData(gain, (ShapeErrorType)i, res);
      if (sed) return sed;
    }
    return nullptr;
  }
  
  if (shapeErrorType == CellShapeError) {
    ShapeErrorData* sed = shapeErrorGetter()->shapeErrorData(hash(), gain, res);
    if (!sed) return nullptr;
    sed->setShapeErrorType(CellShapeError);
    return sed;
  }
  
  if (shapeErrorType == LowGainCellShapeError || shapeErrorType == MedGainCellShapeError || shapeErrorType == HighGainCellShapeError) {
    CaloGain::CaloGain fbGain = (shapeErrorType == LowGainCellShapeError ? CaloGain::LARLOWGAIN :
                                (shapeErrorType == MedGainCellShapeError ? CaloGain::LARMEDIUMGAIN :CaloGain::LARHIGHGAIN));
    ShapeErrorData* sed = shapeErrorGetter()->shapeErrorData(hash(), fbGain, res);
    if (!sed) return nullptr;
    sed->setShapeErrorType(shapeErrorType);
    return sed;
  }
  
  if (shapeErrorType == RingShapeError) {
    ShapeErrorData* sed = shapeErrorGetter()->phiSymShapeErrorData(cellInfo()->globalPhiRing(), gain, res);
    if (!sed) return nullptr;
    sed->setShapeErrorType(RingShapeError);
    return sed;
  }
 
  if (shapeErrorType == LowGainRingShapeError || shapeErrorType == MedGainRingShapeError || shapeErrorType == HighGainRingShapeError) {
    CaloGain::CaloGain fbGain = (shapeErrorType == LowGainRingShapeError ? CaloGain::LARLOWGAIN :
                                (shapeErrorType == MedGainRingShapeError ? CaloGain::LARMEDIUMGAIN :CaloGain::LARHIGHGAIN));
    ShapeErrorData* sed = shapeErrorGetter()->phiSymShapeErrorData(cellInfo()->globalPhiRing(), fbGain, res);
    if (!sed) return nullptr;
    sed->setShapeErrorType(shapeErrorType);
    return sed;
  }
 
  return nullptr;
}

shapeErrorGetter()

const AbsShapeErrorGetter * LArSamples::History::shapeErrorGetter ( ) const

inline

Definition at line 88 of file History.h.

{ return m_shapeErrorGetter; }

sum()

bool History::sum	(	SimpleShape *&	sum,
		SimpleShape *&	reference ) const

Definition at line 127 of file History.cxx.

{
  reference = sum = nullptr;
  if (nData() == 0) return false;
  
  sum = new SimpleShape(m_data[0]->nSamples());
  reference = nullptr;
 
  for (unsigned int j = 0; j < nData(); j++) {
    if (sum->nPoints() !=  m_data[j]->nSamples()) return false;
    for (unsigned int k = 0; k < m_data[j]->nSamples(); k++)
      sum->add(k, m_data[j]->pedestalSubtractedSample(k), m_data[j]->error(k));
    SimpleShape* thisRef = referenceShape(j);
    if (!thisRef) { delete sum; return false; }
    if (reference) {
      reference->add(*thisRef);
      delete thisRef;
    }
    else reference = thisRef;
  }
    
  return true;
}

upstreamEnergy()

double History::upstreamEnergy

(

unsigned int

k

)

const

Definition at line 640 of file History.cxx.

{
  if (!m_interface || !cellInfo()) return -1;
  if (k >= nData() || !data(k)->eventData()) return -1;
  std::vector<unsigned int> upstreamNeighbors;
  if (!m_interface->firstNeighbors(hash(), upstreamNeighbors, cellInfo()->layer() - 1)) return -1;
  if (upstreamNeighbors.empty()) return -1;
  std::vector<const Data*> unData;
  if (!m_interface->data(upstreamNeighbors, *data(k)->eventData(), unData)) return -1;
  double upstreamE = 0;
  for (const Data* data : unData) {
    upstreamE += data->energy();
    delete data;
  }
  return upstreamE;
}

Member Data Documentation

m_cellInfo

CellInfo LArSamples::History::m_cellInfo

private

Definition at line 112 of file History.h.

m_data

std::vector<const Data*> LArSamples::History::m_data

private

Definition at line 111 of file History.h.

m_eventData

std::vector<const EventData*> LArSamples::History::m_eventData

private

Definition at line 113 of file History.h.

m_hash

unsigned int LArSamples::History::m_hash

private

Definition at line 114 of file History.h.

m_interface

const Interface* LArSamples::History::m_interface = nullptr

mutableprivate

Definition at line 116 of file History.h.

m_shapeErrorGetter

const AbsShapeErrorGetter* LArSamples::History::m_shapeErrorGetter

mutableprivate

Definition at line 115 of file History.h.

---

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

• History.h
• History.cxx