#include <Data.h>

Inheritance diagram for LArSamples::Data:

Collaboration diagram for LArSamples::Data:

Public Member Functions
	Data (const DataContainer &container, const EventData &eventData, const History *history, int index)
	Constructor.
	Data (std::unique_ptr< const DataContainer > container, const EventData &eventData, const History *history, int index)
	Data (const Data &other, const EventData eventData=0, const History history=0, int index=-1)
virtual	~Data ()
	Destructor.
CaloGain::CaloGain	gain () const
unsigned int	nSamples () const
const short &	sample (unsigned int i) const
float	pedestalSubtractedSample (unsigned int i) const
std::vector< float >	pedestalSubtractedSamples () const
const EventData &	eventData () const
int	run () const
int	event () const
int	lumiBlock () const
int	bunchId () const
double	energy () const
double	ofcTime () const
double	quality () const
double	pedestal () const
double	pedestalRMS () const
int	status () const
double	adcMax () const
double	meV2ADC () const
double	adc2MeV () const
bool	febError (LArFebErrorSummary::LArFebErrorType error) const
bool	isDisconnected () const
TString	problems (bool sayNone=false) const
TString	description (unsigned int verbosity=1) const
bool	isValid () const
double	footprint () const
short	maxSample () const
short	minSample () const
unsigned int	nPoints () const
double	time (unsigned int i) const
double	value (unsigned int i) const
double	covariance (unsigned int i, unsigned int j) const
unsigned int	nCorrs () const
double	noise () const
double	adcMax_new () const
double	covarianceMeV (unsigned int i, unsigned int j) const
const DataContainer &	container () const
const DataContainer *	dissolve ()
SimpleShape *	referenceShape () const
SimpleShape *	timeAlignedShape () const
const ScaledErrorData *	scaledErrorData () const
const History *	history () const
short	index () const
TVectorD	deltas () const
TVectorD	residuals () const
double	delta (short sample) const
double	residual (short sample) const
double	residualOffset (short sample, bool scale=true) const
double	resCorrN (ShapeErrorType shapeErrorType=BestShapeError) const
double	xi (short sample, ShapeErrorType shapeErrorType=BestShapeError, CaloGain::CaloGain g=CaloGain::UNKNOWNGAIN, bool xip=false) const
double	xiNorm (ShapeErrorType shapeErrorType=BestShapeError) const
double	xip (short sample, ShapeErrorType shapeErrorType=BestShapeError, CaloGain::CaloGain g=CaloGain::UNKNOWNGAIN) const
double	residualError (short sample1, short sample2=-1, bool scale=true) const
double	chi2 () const
bool	calcRefit (double &chi2, double &k, double &dT) const
bool	calcAdjust (double &k, double &dT) const
bool	isPassed (const TString &bitName) const
bool	goodForShapeCorr () const
double	upstreamEnergy () const
double	_sample (const DataFuncArgs &args) const
double	_pedestalSubstractedSample (const DataFuncArgs &args) const
double	_energy (const DataFuncArgs &) const
double	_ofcTime (const DataFuncArgs &) const
double	_quality (const DataFuncArgs &) const
double	_pedestal (const DataFuncArgs &) const
double	_pedestalRMS (const DataFuncArgs &) const
double	_adcMax (const DataFuncArgs &) const
double	_adcMax_new (const DataFuncArgs &) const
double	_noise (const DataFuncArgs &) const
double	_maxValue (const DataFuncArgs &) const
double	_minValue (const DataFuncArgs &) const
double	_null (const DataFuncArgs &) const
int	_run (const DataFuncArgs &) const
int	_event (const DataFuncArgs &) const
int	_lumiBlock (const DataFuncArgs &) const
int	_bunchId (const DataFuncArgs &) const
int	_index (const DataFuncArgs &) const
int	_gain (const DataFuncArgs &) const
double	_delta (const DataFuncArgs &args) const
double	_residual (const DataFuncArgs &args) const
double	_residualOffset (const DataFuncArgs &args) const
double	_resCorrN (const DataFuncArgs &) const
double	_xi (const DataFuncArgs &args) const
double	_xiCell (const DataFuncArgs &args) const
double	_xiRing (const DataFuncArgs &args) const
double	_xiNormCell (const DataFuncArgs &) const
double	_xiNormRing (const DataFuncArgs &) const
double	_xip (const DataFuncArgs &args) const
double	_residualError (const DataFuncArgs &args) const
double	_normResidualOffset (const DataFuncArgs &args) const
double	_normResidualError (const DataFuncArgs &args) const
double	_ofcSigma (const DataFuncArgs &args) const
double	_ofcGamma (const DataFuncArgs &args) const
int	_iEta (const DataFuncArgs &args) const
int	_iPhi (const DataFuncArgs &args) const
int	_region (const DataFuncArgs &args) const
int	_layer (const DataFuncArgs &args) const
int	_calo (const DataFuncArgs &args) const
int	_feedThrough (const DataFuncArgs &args) const
int	_slot (const DataFuncArgs &args) const
int	_channel (const DataFuncArgs &args) const
int	_hash (const DataFuncArgs &args) const
int	_ring (const DataFuncArgs &args) const
int	_nData (const DataFuncArgs &args) const
int	_maxPosition (const DataFuncArgs &) const
double	_count (const DataFuncArgs &) const
double	_chi2 (const DataFuncArgs &) const
double	_chi2_noCorr (const DataFuncArgs &args) const
double	_chi2_ringCorr (const DataFuncArgs &args) const
double	_chi2_cellCorr (const DataFuncArgs &args) const
double	_chi2_k (const DataFuncArgs &args) const
double	_refitDeltaT (const DataFuncArgs &args) const
double	_refitScale (const DataFuncArgs &args) const
double	_refitChi2 (const DataFuncArgs &args) const
double	_adjDeltaT (const DataFuncArgs &args) const
double	_adjScale (const DataFuncArgs &args) const
double	_x (const DataFuncArgs &args) const
double	_y (const DataFuncArgs &args) const
double	_z (const DataFuncArgs &args) const
double	_rt (const DataFuncArgs &args) const
double	_eta (const DataFuncArgs &args) const
double	_phi (const DataFuncArgs &args) const
double	_timeNoTOF (const DataFuncArgs &args) const
double	_timeForSplash1 (const DataFuncArgs &args) const
double	_timeForSplash2 (const DataFuncArgs &args) const
double	_peakSignif (const DataFuncArgs &args) const
double	_energySignif (const DataFuncArgs &args) const
int	_goodForShapeCorr (const DataFuncArgs &) const
double	_upstreamE (const DataFuncArgs &) const
double	_upstreamERatio (const DataFuncArgs &args) const
double	_chi2Anomaly (const DataFuncArgs &args) const
double	timeClCorr () const
double	timeECorr () const
double	_timeECorr (const DataFuncArgs &args) const
double	_pass (const DataFuncArgs &args) const
void	setCallBacks (const History *history, unsigned int index) const
virtual double	error (unsigned int i) const
TVectorD	values (int lwb, int upb) const
int	findTimeInterval (double time) const
int	interpolate (double time, double &value, double &error) const
bool	interpolate (const AbsShape &other, TVectorD &values, CovMatrix &errors, int lwb=-1, int upb=-1) const
int	interpolateDiff (double time, double &diff) const
bool	interpolateDiff (const AbsShape &other, TVectorD &diffs, int lwb=-1, int upb=-1) const
CovMatrix	covarianceMatrix (int lwb=-1, int upb=-1, const CovMatrix &refErr=CovMatrix(), bool withCorrs=true) const
CovMatrix	covarianceMatrix (unsigned int nPoints, bool withCorrs=true) const
CovMatrix	invCovarianceMatrix (int lwb=-1, int upb=-1, const CovMatrix &refErr=CovMatrix(), bool withCorrs=true) const
CovMatrix	invCovarianceMatrix (unsigned int nPoints, bool withCorrs=true) const
double	maxValue (bool withErrors=false) const
double	minValue (bool withErrors=false) const
int	maxPosition () const
int	minPosition () const
TGraphErrors *	graph (bool timeInUnitOfSamples=false) const
SimpleShape *	resample (unsigned int nPts) const

Static Public Member Functions
static TString	str (LArFebErrorSummary::LArFebErrorType error)
static TString	gainStr (CaloGain::CaloGain gain)
static double	timeShift ()
static void	setTimeShift (double shift=0)

Private Attributes
std::unique_ptr< const DataContainer >	m_container
const EventData &	m_eventData
const History *	m_history
unsigned int	m_index

Static Private Attributes
static std::atomic< double >	m_timeShift {0}

Friends
class	DataTweaker

Detailed Description

Definition at line 73 of file Data.h.

Constructor & Destructor Documentation

◆ Data() [1/3]

Data::Data	(	const DataContainer &	container,
		const EventData &	eventData,
		const History *	history,
		int	index )

Constructor.

Definition at line 32 of file Data.cxx.

   : AbsShape(), m_container(std::make_unique<DataContainer>(container)), m_eventData(eventData),
     m_history(history), m_index(index)
{
  ClassCounts::incrementInstanceCount("Data");
}

◆ Data() [2/3]

Data::Data	(	std::unique_ptr< const DataContainer >	container,
		const EventData &	eventData,
		const History *	history,
		int	index )

Definition at line 40 of file Data.cxx.

   : AbsShape(), m_container(std::move(container)), m_eventData(eventData),
     m_history(history), m_index(index) 
{
  ClassCounts::incrementInstanceCount("Data"); 
}

◆ Data() [3/3]

Data::Data	(	const Data &	other,
		const EventData *	eventData = 0,
		const History *	history = 0,
		int	index = -1 )

Definition at line 48 of file Data.cxx.

  : AbsShape(), m_container(new DataContainer(other.container())), 
    m_eventData(eventData ? *eventData : other.eventData()),
    m_history(history ? history : other.m_history), m_index(index >= 0 ? index : other.m_index)
{ 
  ClassCounts::incrementInstanceCount("Data"); 
}

◆ ~Data()

Data::~Data ( )

virtual

Destructor.

Definition at line 57 of file Data.cxx.

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

Member Function Documentation

◆ _adcMax()

double LArSamples::Data::_adcMax ( const DataFuncArgs & ) const

inline

Definition at line 198 of file Data.h.

198{ return adcMax(); }

LArSamples::Data::adcMax

double adcMax() const

Definition Data.h:126

◆ _adcMax_new()

double LArSamples::Data::_adcMax_new ( const DataFuncArgs & ) const

inline

Definition at line 199 of file Data.h.

199{ return adcMax_new(); }

LArSamples::Data::adcMax_new

double adcMax_new() const

Definition Data.cxx:76

◆ _adjDeltaT()

double Data::_adjDeltaT ( const DataFuncArgs & args ) const

Definition at line 288 of file Data.cxx.

{
  double k, dT;
  if (!calcAdjust(k, dT)) return -9999;
  return dT;
}

◆ _adjScale()

double Data::_adjScale ( const DataFuncArgs & args ) const

Definition at line 296 of file Data.cxx.

{
  double k, dT;
  if (!calcAdjust(k, dT)) return -1;
  return k;
}

◆ _bunchId()

int LArSamples::Data::_bunchId ( const DataFuncArgs & ) const

inline

Definition at line 208 of file Data.h.

208{ return bunchId(); }

LArSamples::Data::bunchId

int bunchId() const

Definition Data.cxx:30

◆ _calo()

int Data::_calo ( const DataFuncArgs & args ) const

Definition at line 325 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->calo();
}

◆ _channel()

int Data::_channel ( const DataFuncArgs & args ) const

Definition at line 353 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->channel();
}

◆ _chi2()

double LArSamples::Data::_chi2 ( const DataFuncArgs & ) const

inline

Definition at line 241 of file Data.h.

241{ return chi2(); }

LArSamples::Data::chi2

double chi2() const

Definition Data.cxx:204

◆ _chi2_cellCorr()

double Data::_chi2_cellCorr ( const DataFuncArgs & args ) const

Definition at line 225 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->chi2(m_index, -1, -1, DefaultChi2, CellShapeError);
}

◆ _chi2_k()

double Data::_chi2_k ( const DataFuncArgs & args ) const

Definition at line 232 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->chi2_k(m_index, args.x1);
}

◆ _chi2_noCorr()

double Data::_chi2_noCorr ( const DataFuncArgs & args ) const

Definition at line 211 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->chi2(m_index, -1, -1, DefaultChi2, NoShapeError);
}

◆ _chi2_ringCorr()

double Data::_chi2_ringCorr ( const DataFuncArgs & args ) const

Definition at line 218 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->chi2(m_index, -1, -1, DefaultChi2, RingShapeError);
}

◆ _chi2Anomaly()

double Data::_chi2Anomaly ( const DataFuncArgs & args ) const

Definition at line 596 of file Data.cxx.

{
  if (!m_history) return -1;
  return m_history->chi2Anomaly(args.x1);
}

◆ _count()

double LArSamples::Data::_count ( const DataFuncArgs & ) const

inline

Definition at line 239 of file Data.h.

239{ return 1; }

◆ _delta()

double LArSamples::Data::_delta ( const DataFuncArgs & args ) const

inline

Definition at line 211 of file Data.h.

211{ return delta(args.i1); }

LArSamples::Data::delta

double delta(short sample) const

Definition Data.cxx:430

◆ _energy()

double LArSamples::Data::_energy ( const DataFuncArgs & ) const

inline

Definition at line 193 of file Data.h.

193{ return energy(); }

LArSamples::Data::energy

double energy() const

Definition Data.h:110

◆ _energySignif()

double Data::_energySignif ( const DataFuncArgs & args ) const

Definition at line 408 of file Data.cxx.

{
  return energy()/noise();
}

◆ _eta()

double Data::_eta ( const DataFuncArgs & args ) const

Definition at line 502 of file Data.cxx.

502{ return m_history->cellInfo()->eta(); }

◆ _event()

int LArSamples::Data::_event ( const DataFuncArgs & ) const

inline

Definition at line 206 of file Data.h.

206{ return event(); }

LArSamples::Data::event

int event() const

Definition Data.cxx:28

◆ _feedThrough()

int Data::_feedThrough ( const DataFuncArgs & args ) const

Definition at line 339 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->feedThrough();
}

◆ _gain()

int LArSamples::Data::_gain ( const DataFuncArgs & ) const

inline

Definition at line 210 of file Data.h.

210{ return gain(); }

LArSamples::Data::gain

CaloGain::CaloGain gain() const

Definition Data.h:87

◆ _goodForShapeCorr()

int LArSamples::Data::_goodForShapeCorr ( const DataFuncArgs & ) const

inline

Definition at line 262 of file Data.h.

262{ return (goodForShapeCorr() ? 1 : 0); }

LArSamples::Data::goodForShapeCorr

bool goodForShapeCorr() const

Definition Data.cxx:536

◆ _hash()

int Data::_hash ( const DataFuncArgs & args ) const

Definition at line 360 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->hash();
}

◆ _iEta()

int Data::_iEta ( const DataFuncArgs & args ) const

Definition at line 304 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->iEta();
}

◆ _index()

int LArSamples::Data::_index ( const DataFuncArgs & ) const

inline

Definition at line 209 of file Data.h.

209{ return m_index; }

◆ _iPhi()

int Data::_iPhi ( const DataFuncArgs & args ) const

Definition at line 311 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->iPhi();
}

◆ _layer()

int Data::_layer ( const DataFuncArgs & args ) const

Definition at line 332 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->layer();
}

◆ _lumiBlock()

int LArSamples::Data::_lumiBlock ( const DataFuncArgs & ) const

inline

Definition at line 207 of file Data.h.

207{ return lumiBlock(); }

LArSamples::Data::lumiBlock

int lumiBlock() const

Definition Data.cxx:29

◆ _maxPosition()

int LArSamples::Data::_maxPosition ( const DataFuncArgs & ) const

inline

Definition at line 237 of file Data.h.

237{ return maxPosition(); }

LArSamples::AbsShape::maxPosition

int maxPosition() const

Definition AbsShape.cxx:52

◆ _maxValue()

double LArSamples::Data::_maxValue ( const DataFuncArgs & ) const

inline

Definition at line 201 of file Data.h.

201{ return maxValue(); }

maxValue

#define maxValue(current, test)

Definition CompoundLayerMaterialCreator.h:22

◆ _minValue()

double LArSamples::Data::_minValue ( const DataFuncArgs & ) const

inline

Definition at line 202 of file Data.h.

202{ return minValue(); }

minValue

#define minValue(current, test)

Definition CompoundLayerMaterialCreator.h:21

◆ _nData()

int Data::_nData ( const DataFuncArgs & args ) const

Definition at line 374 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->nData();
}

◆ _noise()

double LArSamples::Data::_noise ( const DataFuncArgs & ) const

inline

Definition at line 200 of file Data.h.

200{ return noise(); }

◆ _normResidualError()

double LArSamples::Data::_normResidualError ( const DataFuncArgs & args ) const

inline

Definition at line 223 of file Data.h.

223{ return residualError(args.i1, -1, false); }

LArSamples::Data::residualError

double residualError(short sample1, short sample2=-1, bool scale=true) const

Definition Data.cxx:490

◆ _normResidualOffset()

double LArSamples::Data::_normResidualOffset ( const DataFuncArgs & args ) const

inline

Definition at line 222 of file Data.h.

222{ return residualOffset(args.i1, false); }

LArSamples::Data::residualOffset

double residualOffset(short sample, bool scale=true) const

Definition Data.cxx:454

◆ _null()

double LArSamples::Data::_null ( const DataFuncArgs & ) const

inline

Definition at line 204 of file Data.h.

204{ return 0; }

◆ _ofcGamma()

double Data::_ofcGamma ( const DataFuncArgs & args ) const

Definition at line 577 of file Data.cxx.

{
  double result = -9999;
  if (!m_history) return result;
  std::unique_ptr<OFC> ofc (m_history->ofc(m_index));
  if (ofc->isInRange(args.i1) && ofc->isInRange(args.i2))
    result = ofc->Gamma()(args.i1, args.i2);
  return result;
}

◆ _ofcSigma()

double Data::_ofcSigma ( const DataFuncArgs & args ) const

Definition at line 566 of file Data.cxx.

{
  double result = -9999;
  if (!m_history) return result;
  std::unique_ptr<OFC> ofc (m_history->ofc(m_index));
  if (ofc->isInRange(args.i1))
    result = TMath::Sqrt(ofc->Gamma()(args.i1, args.i1));
  return result;
}

◆ _ofcTime()

double LArSamples::Data::_ofcTime ( const DataFuncArgs & ) const

inline

Definition at line 194 of file Data.h.

194{ return ofcTime(); }

LArSamples::Data::ofcTime

double ofcTime() const

Definition Data.h:113

◆ _pass()

double LArSamples::Data::_pass ( const DataFuncArgs & args ) const

inline

Definition at line 272 of file Data.h.

272{ return isPassed(args.str); }

LArSamples::Data::isPassed

bool isPassed(const TString &bitName) const

Definition Data.cxx:530

◆ _peakSignif()

double Data::_peakSignif ( const DataFuncArgs & args ) const

Definition at line 402 of file Data.cxx.

{
  return maxValue()/pedestalRMS();
}

◆ _pedestal()

double LArSamples::Data::_pedestal ( const DataFuncArgs & ) const

inline

Definition at line 196 of file Data.h.

196{ return pedestal(); }

LArSamples::Data::pedestal

double pedestal() const

Definition Data.h:119

◆ _pedestalRMS()

double LArSamples::Data::_pedestalRMS ( const DataFuncArgs & ) const

inline

Definition at line 197 of file Data.h.

197{ return pedestalRMS(); }

◆ _pedestalSubstractedSample()

double LArSamples::Data::_pedestalSubstractedSample ( const DataFuncArgs & args ) const

inline

Definition at line 192 of file Data.h.

192{ return pedestalSubtractedSample(args.i1); }

LArSamples::Data::pedestalSubtractedSample

float pedestalSubtractedSample(unsigned int i) const

Definition Data.h:94

◆ _phi()

double Data::_phi ( const DataFuncArgs & args ) const

Definition at line 503 of file Data.cxx.

503{ return m_history->cellInfo()->phi(); }

◆ _quality()

double LArSamples::Data::_quality ( const DataFuncArgs & ) const

inline

Definition at line 195 of file Data.h.

195{ return quality(); }

LArSamples::Data::quality

double quality() const

Definition Data.h:116

◆ _refitChi2()

double Data::_refitChi2 ( const DataFuncArgs & args ) const

Definition at line 268 of file Data.cxx.

{
  double chi2, k, dT;
  if (!calcRefit(chi2, k, dT)) return -1;
  return chi2;
}

◆ _refitDeltaT()

double Data::_refitDeltaT ( const DataFuncArgs & args ) const

Definition at line 252 of file Data.cxx.

{
  double chi2, k, dT;
  if (!calcRefit(chi2, k, dT)) return -9999;
  return dT;
}

◆ _refitScale()

double Data::_refitScale ( const DataFuncArgs & args ) const

Definition at line 260 of file Data.cxx.

{
  double chi2, k, dT;
  if (!calcRefit(chi2, k, dT)) return -1;
  return k;
}

◆ _region()

int Data::_region ( const DataFuncArgs & args ) const

Definition at line 318 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->region();
}

◆ _resCorrN()

double LArSamples::Data::_resCorrN ( const DataFuncArgs & ) const

inline

Definition at line 214 of file Data.h.

214{ return resCorrN(CellShapeError); }

LArSamples::Data::resCorrN

double resCorrN(ShapeErrorType shapeErrorType=BestShapeError) const

Definition Data.cxx:446

◆ _residual()

double LArSamples::Data::_residual ( const DataFuncArgs & args ) const

inline

Definition at line 212 of file Data.h.

212{ return residual(args.i1); }

LArSamples::Data::residual

double residual(short sample) const

Definition Data.cxx:438

◆ _residualError()

double LArSamples::Data::_residualError ( const DataFuncArgs & args ) const

inline

Definition at line 221 of file Data.h.

221{ return residualError(args.i1); }

◆ _residualOffset()

double LArSamples::Data::_residualOffset ( const DataFuncArgs & args ) const

inline

Definition at line 213 of file Data.h.

213{ return residualOffset(args.i1); }

◆ _ring()

int Data::_ring ( const DataFuncArgs & args ) const

Definition at line 367 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->globalPhiRing();
}

◆ _rt()

double Data::_rt ( const DataFuncArgs & args ) const

Definition at line 501 of file Data.cxx.

501{ return m_history->cellInfo()->rt(); }

◆ _run()

int LArSamples::Data::_run ( const DataFuncArgs & ) const

inline

Definition at line 205 of file Data.h.

205{ return run(); }

LArSamples::Data::run

int run() const

Definition Data.cxx:27

◆ _sample()

double LArSamples::Data::_sample ( const DataFuncArgs & args ) const

inline

Definition at line 191 of file Data.h.

191{ return sample(args.i1); }

LArSamples::Data::sample

const short & sample(unsigned int i) const

Definition Data.h:93

◆ _slot()

int Data::_slot ( const DataFuncArgs & args ) const

Definition at line 346 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->cellInfo()->slot();
}

◆ _timeECorr()

double Data::_timeECorr ( const DataFuncArgs & args ) const

Definition at line 641 of file Data.cxx.

{
  // Modify this to add Isabelle's correction
  return ofcTime();
}

◆ _timeForSplash1()

double Data::_timeForSplash1 ( const DataFuncArgs & args ) const

Definition at line 388 of file Data.cxx.

{
  if (!m_history) return -9999;
  return _timeNoTOF(args) + m_history->cellInfo()->position().z()/(TMath::C()*1E-6); // 10^-6 to convert to mm/ns
}

◆ _timeForSplash2()

double Data::_timeForSplash2 ( const DataFuncArgs & args ) const

Definition at line 395 of file Data.cxx.

{
  if (!m_history) return -9999;
  return _timeNoTOF(args) - m_history->cellInfo()->position().z()/(TMath::C()*1E-6); // 10^-6 to convert to mm/ns
}

◆ _timeNoTOF()

double Data::_timeNoTOF ( const DataFuncArgs & args ) const

Definition at line 381 of file Data.cxx.

{
  if (!m_history) return -9999;
  return ofcTime() + m_history->cellInfo()->position().Mag()/(TMath::C()*1E-6); // 10^-6 to convert to mm/ns
}

◆ _upstreamE()

double LArSamples::Data::_upstreamE ( const DataFuncArgs & ) const

inline

Definition at line 263 of file Data.h.

263{ return upstreamEnergy(); }

LArSamples::Data::upstreamEnergy

double upstreamEnergy() const

Definition Data.cxx:560

◆ _upstreamERatio()

double Data::_upstreamERatio ( const DataFuncArgs & args ) const

Definition at line 588 of file Data.cxx.

{ 
  double upe = upstreamEnergy();
  if (upe < 0 || energy() == 0) return -1;
  return upe/energy(); 
}

◆ _x()

double Data::_x ( const DataFuncArgs & args ) const

Definition at line 498 of file Data.cxx.

498{ return m_history->cellInfo()->position().X(); }

◆ _xi()

double LArSamples::Data::_xi ( const DataFuncArgs & args ) const

inline

Definition at line 215 of file Data.h.

215{ return xi(args.i1); }

LArSamples::Data::xi

double xi(short sample, ShapeErrorType shapeErrorType=BestShapeError, CaloGain::CaloGain g=CaloGain::UNKNOWNGAIN, bool xip=false) const

Definition Data.cxx:462

◆ _xiCell()

double LArSamples::Data::_xiCell ( const DataFuncArgs & args ) const

inline

Definition at line 216 of file Data.h.

216{ return xi(args.i1, CellShapeError, args.i2 >= 0 ? (CaloGain::CaloGain)args.i2 : CaloGain::UNKNOWNGAIN); }

CaloGain::CaloGain

CaloGain

Definition CaloGain.h:11

CaloGain::UNKNOWNGAIN

@ UNKNOWNGAIN

Definition CaloGain.h:21

◆ _xiNormCell()

double LArSamples::Data::_xiNormCell ( const DataFuncArgs & ) const

inline

Definition at line 218 of file Data.h.

218{ return xiNorm(CellShapeError); }

LArSamples::Data::xiNorm

double xiNorm(ShapeErrorType shapeErrorType=BestShapeError) const

Definition Data.cxx:474

◆ _xiNormRing()

double LArSamples::Data::_xiNormRing ( const DataFuncArgs & ) const

inline

Definition at line 219 of file Data.h.

219{ return xiNorm(RingShapeError); }

◆ _xip()

double LArSamples::Data::_xip ( const DataFuncArgs & args ) const

inline

Definition at line 220 of file Data.h.

220{ return xip(args.i1); }

LArSamples::Data::xip

double xip(short sample, ShapeErrorType shapeErrorType=BestShapeError, CaloGain::CaloGain g=CaloGain::UNKNOWNGAIN) const

Definition Data.cxx:484

◆ _xiRing()

double LArSamples::Data::_xiRing ( const DataFuncArgs & args ) const

inline

Definition at line 217 of file Data.h.

217{ return xi(args.i1, RingShapeError, args.i2 >= 0 ? (CaloGain::CaloGain)args.i2 : CaloGain::UNKNOWNGAIN); }

◆ _y()

double Data::_y ( const DataFuncArgs & args ) const

Definition at line 499 of file Data.cxx.

499{ return m_history->cellInfo()->position().Y(); }

◆ _z()

double Data::_z ( const DataFuncArgs & args ) const

Definition at line 500 of file Data.cxx.

500{ return m_history->cellInfo()->position().Z(); }

◆ adc2MeV()

double LArSamples::Data::adc2MeV ( ) const

inline

Definition at line 128 of file Data.h.

128{ return (adcMax() > 0 && energy() > 0 ? energy()/adcMax() : -1); }

◆ adcMax()

double LArSamples::Data::adcMax ( ) const

inline

Returns: OFC peak parameters

Definition at line 126 of file Data.h.

126{ return container().adcMax(); }

LArSamples::DataContainer::adcMax

float adcMax() const

Definition DataContainer.h:86

◆ adcMax_new()

double Data::adcMax_new ( ) const

Definition at line 76 of file Data.cxx.

{
  return ((energy()/noise())*pedestalRMS());
 
}

◆ bunchId()

int Data::bunchId ( ) const

Definition at line 30 of file Data.cxx.

30{ return eventData().bunchId(); }

LArSamples::EventData::bunchId

unsigned int bunchId() const

Definition LArCalorimeter/LArCafJobs/LArCafJobs/EventData.h:50

◆ calcAdjust()

bool Data::calcAdjust	(	double &	k,
		double &	dT ) const

Definition at line 276 of file Data.cxx.

{
  if (!m_history) return false;
  std::unique_ptr<OFC> ofc (m_history->ofc(m_index));
  if (!ofc) return false;
  
  k = ofc->A(*this)/adcMax();
  dT = ofc->time(*this);
  return true;
}

◆ calcRefit()

bool Data::calcRefit	(	double &	chi2,
		double &	k,
		double &	dT ) const

Definition at line 239 of file Data.cxx.

{
  ShapeFitter fitter;
  std::unique_ptr<SimpleShape> reference (referenceShape());
  if (!reference){
   return false;
  }
  std::unique_ptr<const ScaledErrorData> sed (scaledErrorData());
  bool result = fitter.fit(*this, *reference, k, dT, chi2, sed.get());
  return result;
}

◆ chi2()

double Data::chi2 ( ) const

Definition at line 204 of file Data.cxx.

{
  if (!m_history) return 0;
  return m_history->chi2(m_index);
}

◆ container()

const DataContainer & LArSamples::Data::container ( ) const

inline

Definition at line 158 of file Data.h.

158{ return *m_container; }

◆ covariance()

double Data::covariance	(	unsigned int	i,
		unsigned int	j ) const

virtual

Implements LArSamples::AbsShape.

Definition at line 82 of file Data.cxx.

{
  if (i >= nSamples() || j >= nSamples()) return 0;
  const std::vector<float>& corrs = m_container->corrs();
  unsigned int index = TMath::Abs(int(i) - int(j));  // we use time-translation-invariant correlations...
  if (index >= corrs.size()) return 0;
  double corr = corrs[index];
  return corr*pedestalRMS()*pedestalRMS()/(noise()*noise()); // hack to correct the noise per sample/ADC_max noise error
}

◆ covarianceMatrix() [1/2]

CovMatrix AbsShape::covarianceMatrix	(	int	lwb = -1,
		int	upb = -1,
		const CovMatrix &	refErr = CovMatrix(),
		bool	withCorrs = true ) const

inherited

Definition at line 198 of file AbsShape.cxx.

{
  if (lwb < 0 || upb < 0) { lwb = 0; upb = nPoints() - 1; }
  CovMatrix c(lwb, upb);
  for (int i = lwb; i <= upb; i++) {
    for (int j = lwb; j <= upb; j++) {
      double cov = (withCorrs || i == j ? covariance(i, j) : 0);
      if (refErr.GetNrows() > 0) cov += refErr(i, j);
      c(i, j) = cov;
    }
  }
  return c;
}

◆ covarianceMatrix() [2/2]

CovMatrix AbsShape::covarianceMatrix	(	unsigned int	nPoints,
		bool	withCorrs = true ) const

inherited

Definition at line 223 of file AbsShape.cxx.

{
  return covarianceMatrix(0, n - 1, CovMatrix(), withCorrs);
}

◆ covarianceMeV()

double Data::covarianceMeV	(	unsigned int	i,
		unsigned int	j ) const

Definition at line 93 of file Data.cxx.

{
  double adc2mev = adc2MeV();
  if (adc2mev < 0) return 0;
  return covariance(i,j)*adc2mev*adc2mev; // we use time-translation-invariant correlations...
  //return (i == j ? m_corrs[0] : 0); // temporary
}

◆ delta()

double Data::delta ( short sample ) const

Definition at line 430 of file Data.cxx.

{
  double val = 0;
  if (!m_history || !m_history->delta(m_index, sample, val)) return -9999;
  return val;
}

◆ deltas()

TVectorD Data::deltas ( ) const

Definition at line 414 of file Data.cxx.

{
  if (!m_history) return TVectorD();
  return m_history->deltas(m_index);
}

◆ description()

TString Data::description ( unsigned int verbosity = 1 ) const

Definition at line 515 of file Data.cxx.

{
  TString desc = "";
  desc += Form("%d samples, ", nSamples());
  desc += "gain = " + gainStr(gain()) + ", ";
  desc += Form("energy = %6.1f MeV, ", energy());
  if (verbosity >= 2) desc += Form("ADCMax = %6.1f, ", adcMax());
  desc += Form("time = %6.1f ns, ", ofcTime());
  if (verbosity >= 2) desc += Form("chi^2/n = %6.1f, ", chi2());
  desc += Form("run = %d, event = %d, ", run(), event());
  desc += Form("problems = %s", problems(true).Data());
  return desc;
}

◆ dissolve()

const DataContainer * Data::dissolve ( )

Definition at line 63 of file Data.cxx.

{
  const DataContainer* container = m_container.release();
  delete this;
  return container;
}

◆ energy()

double LArSamples::Data::energy ( ) const

inline

Returns: energy

Definition at line 110 of file Data.h.

110{ return container().energy(); }

LArSamples::DataContainer::energy

float energy() const

Definition DataContainer.h:67

◆ error()

double AbsShape::error ( unsigned int i ) const

virtualinherited

Reimplemented in LArSamples::ScaledShiftedShape.

Definition at line 24 of file AbsShape.cxx.

{ 
  return TMath::Sqrt(covariance(i, i)); 
}

◆ event()

int Data::event ( ) const

Returns: event

Definition at line 28 of file Data.cxx.

28{ return eventData().event(); }

LArSamples::EventData::event

int event() const

Definition LArCalorimeter/LArCafJobs/LArCafJobs/EventData.h:44

◆ eventData()

const EventData & LArSamples::Data::eventData ( ) const

inline

Definition at line 97 of file Data.h.

97{ return m_eventData; }

◆ febError()

bool Data::febError ( LArFebErrorSummary::LArFebErrorType error ) const

Definition at line 118 of file Data.cxx.

{
  return ((status() >> (error + 17)) & 0x1) != 0;
}

◆ findTimeInterval()

int AbsShape::findTimeInterval ( double time ) const

inherited

Definition at line 71 of file AbsShape.cxx.

{
  if (t < time(0)) return -1;
  if (t > time(nPoints() - 1)) return int(nPoints()) - 1;
  
  for (unsigned int i = 0; i < nPoints() - 1; i++)
    if (time(i) <= t && t < time(i+1)) return i;
 
  return nPoints() - 1;
}

◆ footprint()

double LArSamples::Data::footprint ( ) const

inline

Definition at line 140 of file Data.h.

140{ return sizeof(*this) + nSamples()*sizeof(short); }

xAOD::short

short

Definition Vertex_v1.cxx:165

◆ gain()

CaloGain::CaloGain LArSamples::Data::gain ( ) const

inline

Returns: gain value

Definition at line 87 of file Data.h.

87{ return container().gain(); }

LArSamples::DataContainer::gain

CaloGain::CaloGain gain() const

Definition DataContainer.h:54

◆ gainStr()

TString Data::gainStr ( CaloGain::CaloGain gain )

static

Definition at line 506 of file Data.cxx.

{
  if (gain == CaloGain::LARHIGHGAIN) return "HIGH";
  if (gain == CaloGain::LARMEDIUMGAIN) return "MEDIUM";
  if (gain == CaloGain::LARLOWGAIN) return "LOW";
  return "";
}

◆ goodForShapeCorr()

bool Data::goodForShapeCorr ( ) const

Definition at line 536 of file Data.cxx.

{
  if (adcMax() < 0) return false;
 
  CaloId caloId = (m_history ? m_history->cellInfo()->calo()  : UNKNOWN_CALO);
  //short   layer = (m_history ? m_history->cellInfo()->layer() : -1);
 
  //if (Id::matchCalo(caloId, FCAL)       && energy()/noise() < 250) return false;
  if (Id::matchCalo(caloId, FCAL)) return false;
  if (Id::matchCalo(caloId, HEC)        && energy()/noise() < 20) return false;
  if (Id::matchCalo(caloId, EMEC_INNER) && energy()/noise() < 250) return false;
  if (Id::matchCalo(caloId, EMEC_OUTER) && energy()/noise() < 50) return false;
  if (Id::matchCalo(caloId, EMB)        && energy()/noise() < 30) return false;
 
  if (Id::matchCalo(caloId, HEC)) {
    if (TMath::Abs(ofcTime()) > 10) return false;
  }
  else {
    if (TMath::Abs(ofcTime()) > 5) return false;
  }
 
  return true;
}

◆ graph()

TGraphErrors * AbsShape::graph ( bool timeInUnitOfSamples = false ) const

inherited

Definition at line 186 of file AbsShape.cxx.

{
  TGraphErrors* graph = new TGraphErrors(nPoints());
  double unit = (timeInUnitOfSamples ? Definitions::samplingInterval : 1);
  for (unsigned int i = 0; i < nPoints(); i++) {
    graph->SetPoint(i, time(i)/unit, value(i));    
    graph->SetPointError(i, 0, error(i));
  }
  return graph;
}

◆ history()

const History * LArSamples::Data::history ( ) const

inline

Definition at line 167 of file Data.h.

167{ return m_history; }

◆ index()

short LArSamples::Data::index ( ) const

inline

Definition at line 168 of file Data.h.

168{ return m_index; }

◆ interpolate() [1/2]

bool AbsShape::interpolate	(	const AbsShape &	other,
		TVectorD &	values,
		CovMatrix &	errors,
		int	lwb = -1,
		int	upb = -1 ) const

inherited

Definition at line 147 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = other.nPoints() - 1;
  values.ResizeTo(lwb, upb);
  errors.ResizeTo(lwb, upb, lwb, upb);  
  int actualLwb = -1, actualUpb = upb;
  for (int i = lwb; i <= upb; i++) {
    double val, err;
    int stat = interpolate(other.time(i), val, err);
    if (stat ==  0 && actualLwb < 0) actualLwb = i;
    if (stat == +1) { actualUpb = i - 1; break; }
    values(i) = val;
    errors(i, i) = err*err;
  }
  values.ResizeTo(actualLwb, actualUpb);
  errors.ResizeTo(actualLwb, actualUpb, actualLwb, actualUpb);  
  return (actualLwb >= 0);
}

◆ interpolate() [2/2]

int AbsShape::interpolate	(	double	time,
		double &	value,
		double &	error ) const

inherited

Definition at line 83 of file AbsShape.cxx.

{
  if (nPoints() < 2) return -1;
  
  int status = 0;
  int i1 = findTimeInterval(t);
  
  if (i1 < 0) {
    status = -1;
    i1 = 0;
  }
 
  if (i1 > int(nPoints()) - 2) {
    status = 1;
    i1 = nPoints() - 2;
  }
 
  double delta = (t - time(i1))/(time(i1+1) - time(i1));
  val = value(i1) + delta*(value(i1 + 1) - value(i1));
  err = error(i1) + delta*(error(i1 + 1) - error(i1));
  //cout << t << " " << i1 << " " << time(i1) << " " << delta << " " << value(i1) << " " << value(i1+1) << " " << val << " " << status << endl;
  return status;
}

◆ interpolateDiff() [1/2]

bool AbsShape::interpolateDiff	(	const AbsShape &	other,
		TVectorD &	diffs,
		int	lwb = -1,
		int	upb = -1 ) const

inherited

Definition at line 168 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = other.nPoints() - 1;
  diffs.ResizeTo(lwb, upb);
  int actualLwb = -1, actualUpb = upb;  
  for (int i = lwb; i <= upb; i++) {
    double diff;
    int stat = interpolateDiff(other.time(i), diff);
    if (stat ==  0 && actualLwb < 0) actualLwb = i;
    if (stat == +1) { actualUpb = i - 1; break; }
    diffs(i) = diff;
  }
  diffs.ResizeTo(actualLwb, actualUpb);
  return (actualLwb >= 0);
}

◆ interpolateDiff() [2/2]

int AbsShape::interpolateDiff	(	double	time,
		double &	diff ) const

inherited

Definition at line 108 of file AbsShape.cxx.

{
  int status = 0;
  int i1 = findTimeInterval(t);  
  if (i1 < 0) {
    i1 = 0;
    status = -1;
  }  
  if (i1 > int(nPoints()) - 1) {
    i1 = nPoints() - 1;
    status = 1;
  }  
  if (i1 > 1 && i1 < int(nPoints()) - 2)
    diff = (value(i1 - 2) - 8*value(i1 - 1) + 8*value(i1 + 1) - value(i1 + 2))/(time(i1 + 2) - time(i1 - 2))/3;
  else
    if (i1 > 0) 
      if (i1 <  int(nPoints()) - 1)
        diff = (value(i1 + 1) - value(i1 - 1))/(time(i1 + 1) - time(i1 - 1));
      else
        diff = (value(i1) - value(i1 - 1))/(time(i1) - time(i1 - 1));
    else 
      diff = (value(i1 + 1) - value(i1))/(time(i1 + 1) - time(i1));
 
  return status;
}

◆ invCovarianceMatrix() [1/2]

CovMatrix AbsShape::invCovarianceMatrix	(	int	lwb = -1,
		int	upb = -1,
		const CovMatrix &	refErr = CovMatrix(),
		bool	withCorrs = true ) const

inherited

Definition at line 215 of file AbsShape.cxx.

{ 
  return covarianceMatrix(lwb, upb, refErr, withCorrs).Invert(); 
}

◆ invCovarianceMatrix() [2/2]

CovMatrix AbsShape::invCovarianceMatrix	(	unsigned int	nPoints,
		bool	withCorrs = true ) const

inherited

Definition at line 230 of file AbsShape.cxx.

{
  return covarianceMatrix(n, withCorrs).Invert();
}

◆ isDisconnected()

bool Data::isDisconnected ( ) const

Definition at line 113 of file Data.cxx.

{
  return ((status() >> 16) & 0x1) != 0;
}

◆ isPassed()

bool Data::isPassed ( const TString & bitName ) const

Definition at line 530 of file Data.cxx.

{ 
  return eventData().isPassed(bitName);
}

◆ isValid()

bool Data::isValid ( ) const

Definition at line 161 of file Data.cxx.

{
  if (!m_container) return false;
  if (!container().isValid()) return false;
  return true;
}

◆ lumiBlock()

int Data::lumiBlock ( ) const

Returns: lumi block

Definition at line 29 of file Data.cxx.

29{ return eventData().lumiBlock(); }

LArSamples::EventData::lumiBlock

unsigned int lumiBlock() const

Definition LArCalorimeter/LArCafJobs/LArCafJobs/EventData.h:47

◆ maxPosition()

int AbsShape::maxPosition ( ) const

inherited

Definition at line 52 of file AbsShape.cxx.

{
  int maxPosition = -1;
  double maxValue = -1;
  for (unsigned int i = 0; i < nPoints(); i++)
    if (value(i) > maxValue) { maxPosition = i; maxValue = value(i); }
  return maxPosition;
}

◆ maxSample()

short Data::maxSample ( ) const

Definition at line 169 of file Data.cxx.

{
  int maxPos = maxPosition();
  if (maxPos >= 0) return sample(maxPos);
  return -1;
}

◆ maxValue()

double AbsShape::maxValue ( bool withErrors = false ) const

inherited

Returns: sample max parameters

Definition at line 30 of file AbsShape.cxx.

{
  double maxValue = -DBL_MAX;
  for (unsigned int i = 0; i < nPoints(); i++) {
    double val = value(i) + (withErrors ? error(i) : 0);
    if (val > maxValue) maxValue = val;
  }
  return maxValue;
}

◆ meV2ADC()

double LArSamples::Data::meV2ADC ( ) const

inline

Definition at line 127 of file Data.h.

127{ return (adcMax() > 0 && energy() > 0 ? adcMax()/energy() : -1); }

◆ minPosition()

int AbsShape::minPosition ( ) const

inherited

Definition at line 61 of file AbsShape.cxx.

{
  int minPosition = -1;
  double minValue = Definitions::none;
  for (unsigned int i = 0; i < nPoints(); i++)
    if (value(i) < minValue) { minPosition = i; minValue = value(i); }
  return minPosition;
}

◆ minSample()

short Data::minSample ( ) const

Definition at line 176 of file Data.cxx.

{
  int minPos = minPosition();
  if (minPos >= 0) return sample(minPos);
  return -1;
}

◆ minValue()

double AbsShape::minValue ( bool withErrors = false ) const

inherited

Definition at line 41 of file AbsShape.cxx.

{
  double minValue = DBL_MAX;
  for (unsigned int i = 0; i < nPoints(); i++) {
    double val = value(i) - (withErrors ? error(i) : 0);
    if (val < minValue) minValue = val;
  }
  return minValue;
}

◆ nCorrs()

unsigned int LArSamples::Data::nCorrs ( ) const

inline

Definition at line 151 of file Data.h.

151{ return container().corrs().size(); }

LArSamples::DataContainer::corrs

const std::vector< float > & corrs() const

Definition DataContainer.h:83

◆ noise()

double Data::noise ( ) const

Definition at line 71 of file Data.cxx.

{ 
  return (m_container->corrs().empty() ? 0 : TMath::Sqrt(m_container->corrs()[0]));
}

◆ nPoints()

unsigned int LArSamples::Data::nPoints ( ) const

inlinevirtual

Returns: data points

Implements LArSamples::AbsShape.

Definition at line 146 of file Data.h.

146{ return nSamples(); }

◆ nSamples()

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

inline

Returns: number of samples

Definition at line 90 of file Data.h.

90{ return container().nSamples(); }

LArSamples::DataContainer::nSamples

unsigned int nSamples() const

Definition DataContainer.h:57

◆ ofcTime()

double LArSamples::Data::ofcTime ( ) const

inline

Returns: time

Definition at line 113 of file Data.h.

113{ return container().ofcTime() + m_timeShift; }

LArSamples::DataContainer::ofcTime

float ofcTime() const

Definition DataContainer.h:70

LArSamples::Data::m_timeShift

static std::atomic< double > m_timeShift

Definition Data.h:287

◆ pedestal()

double LArSamples::Data::pedestal ( ) const

inline

Returns: pedestal

Definition at line 119 of file Data.h.

119{ return container().pedestal(); }

LArSamples::DataContainer::pedestal

float pedestal() const

Definition DataContainer.h:76

◆ pedestalRMS()

double LArSamples::Data::pedestalRMS ( ) const

inline

Definition at line 120 of file Data.h.

120{ return container().pedestalRMS(); }

LArSamples::DataContainer::pedestalRMS

float pedestalRMS() const

Definition DataContainer.h:77

◆ pedestalSubtractedSample()

float LArSamples::Data::pedestalSubtractedSample ( unsigned int i ) const

inline

Definition at line 94 of file Data.h.

94{ return sample(i) - pedestal(); }

◆ pedestalSubtractedSamples()

std::vector< float > Data::pedestalSubtractedSamples ( ) const

Definition at line 102 of file Data.cxx.

{
  std::vector<float> samples;
  
  for (unsigned int i = 0; i < nSamples(); i++)
    samples.push_back(pedestalSubtractedSample(i));
  
  return samples;
}

◆ problems()

TString Data::problems ( bool sayNone = false ) const

Definition at line 142 of file Data.cxx.

{
  if (isDisconnected()) return "DISCONNECTED";
 
/*   
  TString problems = "";
  for (CaloBadChannel::CALOPROBLEM pb = CaloBadChannel::DEAD; pb < CaloBadChannel::NPB; pb = CaloBadChannel::CALOPROBLEM(pb + 1))
    if (problem(pb)) {
      if (problems != "") problems += " ";
      problems += str(pb);
    }
*/  
 
  TString problems = "";
  if (problems =="" && sayNone) problems = "None";
  return problems;
}

◆ quality()

double LArSamples::Data::quality ( ) const

inline

Returns: pedestal

Definition at line 116 of file Data.h.

116{ return container().quality(); }

LArSamples::DataContainer::quality

float quality() const

Definition DataContainer.h:73

◆ referenceShape()

SimpleShape * Data::referenceShape ( ) const

Definition at line 190 of file Data.cxx.

{
  if (!m_history) return nullptr;
  return m_history->referenceShape(m_index);
}

◆ resample()

SimpleShape * AbsShape::resample ( unsigned int nPts ) const

inherited

Definition at line 237 of file AbsShape.cxx.

{
  std::vector<double> values, errors;
  double t0 = time(0);
  double t1 = time(nPoints() - 1);
  double dT = (t1 - t0)/nPts;
  double t = t0;
  for (unsigned int i = 0; i < nPts; i++, t += dT) {
    double val, err;
    int inRange = interpolate(t, val, err);
    if (inRange != 0) return nullptr;
    values.push_back(val);
    errors.push_back(err);
  }
  
  return new SimpleShape(values, errors, dT, t0);
}

◆ resCorrN()

double Data::resCorrN ( ShapeErrorType shapeErrorType = BestShapeError ) const

Definition at line 446 of file Data.cxx.

{
  if (!m_history) return -999;
  const ShapeErrorData* sed = m_history->shapeErrorData(gain(), type);
  return (sed ? sed->n() : 0);
}

◆ residual()

double Data::residual ( short sample ) const

Definition at line 438 of file Data.cxx.

{
  double val = 0;
  if (!m_history || adcMax() < 0 || !m_history->delta(m_index, sample, val)) return -9999;
  return val/adcMax();
}

◆ residualError()

double Data::residualError	(	short	sample1,
		short	sample2 = -1,
		bool	scale = true ) const

Definition at line 490 of file Data.cxx.

{
  double val = 0;
  if (!m_history || !m_history->residualError(m_index, sample1, sample2, val, scale)) return -9999;
  return val;
}

◆ residualOffset()

double Data::residualOffset	(	short	sample,
		bool	scale = true ) const

Definition at line 454 of file Data.cxx.

{
  double val = 0;
  if (!m_history || !m_history->residualOffset(m_index, sample, val, (scale ? -1 : 1))) return -9999;
  return val;
}

◆ residuals()

TVectorD Data::residuals ( ) const

Definition at line 421 of file Data.cxx.

{
  if (!m_history || adcMax() < 0) return TVectorD();
  TVectorD resv = m_history->deltas(m_index);
  for (int i = 0; i < resv.GetNrows(); i++) resv(i) /= adcMax();
  return resv;
}

◆ run()

int Data::run ( ) const

Returns: run

Definition at line 27 of file Data.cxx.

27{ return eventData().run(); }

LArSamples::EventData::run

int run() const

Definition EventData.cxx:59

◆ sample()

const short & LArSamples::Data::sample ( unsigned int i ) const

inline

Returns: a sample value

Definition at line 93 of file Data.h.

93{ return container().sample(i); }

LArSamples::DataContainer::sample

const short & sample(unsigned int i) const

Definition DataContainer.h:61

◆ scaledErrorData()

const ScaledErrorData * Data::scaledErrorData ( ) const

Definition at line 197 of file Data.cxx.

{
  if (!m_history) return nullptr;
  return m_history->scaledErrorData(m_index);
}

◆ setCallBacks()

void LArSamples::Data::setCallBacks	(	const History *	history,
		unsigned int	index ) const

inline

Definition at line 274 of file Data.h.

274{ m_history = history; m_index = index; }

◆ setTimeShift()

void LArSamples::Data::setTimeShift ( double shift = 0 )

inlinestatic

Definition at line 278 of file Data.h.

278{ m_timeShift = shift; }

◆ status()

int LArSamples::Data::status ( ) const

inline

Returns: status

Definition at line 123 of file Data.h.

123{ return m_container->status(); }

◆ str()

TString Data::str ( LArFebErrorSummary::LArFebErrorType error )

static

Definition at line 123 of file Data.cxx.

{
  if (error == LArFebErrorSummary::Parity)         return "Parity";
  if (error == LArFebErrorSummary::BCID)           return "BCID";
  if (error == LArFebErrorSummary::SampleHeader)   return "SampleHeader";
  if (error == LArFebErrorSummary::EVTID)          return "EVTID";
  if (error == LArFebErrorSummary::ScacStatus)     return "ScacStatus";
  if (error == LArFebErrorSummary::ScaOutOfRange)  return "ScaOutOfRange";
  if (error == LArFebErrorSummary::GainMismatch)   return "GainMismatch";
  if (error == LArFebErrorSummary::TypeMismatch)   return "TypeMismatch";
  if (error == LArFebErrorSummary::NumOfSamples)   return "NumOfSamples";
  if (error == LArFebErrorSummary::EmptyDataBlock) return "EmptyDataBlock";
  if (error == LArFebErrorSummary::DspBlockSize)   return "DspBlockSize";
  if (error == LArFebErrorSummary::CheckSum)       return "CheckSum";
  if (error == LArFebErrorSummary::MissingHeader)  return "MissingHeader";
  return "";
}

◆ time()

double LArSamples::Data::time ( unsigned int i ) const

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 147 of file Data.h.

147{ return Definitions::samplingTime(i); }

LArSamples::Definitions::samplingTime

static double samplingTime(unsigned int i)

Definition LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:17

◆ timeAlignedShape()

SimpleShape * Data::timeAlignedShape ( ) const

Definition at line 184 of file Data.cxx.

{
  return new SimpleShape(*this, 1, -ofcTime());
}

◆ timeClCorr()

double Data::timeClCorr ( ) const

Definition at line 602 of file Data.cxx.

{  
  if (!m_history) return -9999.;
  Double_t timeClCorr=-99.;
  
  if(run()==203636) timeClCorr = ofcTime() + 0.04;
  else if(run()==203680) timeClCorr = ofcTime() + 0.109155;
  else if(run()==203719) timeClCorr = ofcTime() + 0.141762;
  else if(run()==203745) timeClCorr = ofcTime() - 0.107597;
  else if(run()==203779) timeClCorr = ofcTime() - 0.0032581;
  else if(run()==203876) timeClCorr = ofcTime() + 0.109328;
  else if(run()==203934) timeClCorr = ofcTime() - 0.0404719;
  else if(run()==204073) timeClCorr = ofcTime() - 0.0244807;
  else if(run()==204153) timeClCorr = ofcTime() + 0.185949;
  else if(run()==204158) timeClCorr = ofcTime() + 0.0383347;
  else if(run()==204240) timeClCorr = ofcTime() + 0.0635586;
  else if(run()==204265) timeClCorr = ofcTime() + 0.13898;
  else if(run()==204416) timeClCorr = ofcTime() + 0.1074;
  else if(run()==204564) timeClCorr = ofcTime() + 0.0369925;
  else if(run()==210302) timeClCorr = ofcTime() - 0.2;
  else if(run()==210308) timeClCorr = ofcTime() - 0.2;
  else timeClCorr = ofcTime();
 
  return timeClCorr;  
}  

◆ timeECorr()

double Data::timeECorr ( ) const

Definition at line 628 of file Data.cxx.

{  
  if (!m_history) return -9999.;
  Double_t timeCorr=-99.;  
      
  //   Int_t calo =  m_history->cellInfo()->calo();
  //   Double_t p0, p1, p2, p3;
  
  timeCorr=ofcTime();
  
  return timeCorr;  
}  

◆ timeShift()

double LArSamples::Data::timeShift ( )

inlinestatic

Definition at line 277 of file Data.h.

277{ return m_timeShift; }

◆ upstreamEnergy()

double Data::upstreamEnergy ( ) const

Definition at line 560 of file Data.cxx.

{
  if (!m_history) return -1;
  return m_history->upstreamEnergy(m_index);
}

◆ value()

double LArSamples::Data::value ( unsigned int i ) const

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 148 of file Data.h.

148{ return pedestalSubtractedSample(i); }

◆ values()

TVectorD AbsShape::values	(	int	lwb,
		int	upb ) const

inherited

Definition at line 135 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = (int)nPoints() - 1;
  if (upb >= (int)nPoints()) upb = (int)nPoints() - 1;
  TVectorD vals(lwb, upb);
  for (int i = lwb; i <= upb; i++) 
    vals(i) = value(i);
  return vals;
}

◆ xi()

double Data::xi	(	short	sample,
		ShapeErrorType	shapeErrorType = BestShapeError,
		CaloGain::CaloGain	g = CaloGain::UNKNOWNGAIN,
		bool	xip = false ) const

Definition at line 462 of file Data.cxx.

{
  if (!m_history) return -999;
  //njpb don't see the point of next line, commenting out
  //if (g != CaloGain::UNKNOWNGAIN && type != RingShapeError) type = CellShapeError;
  std::unique_ptr<const ShapeErrorData> sed (m_history->shapeErrorData(g == CaloGain::UNKNOWNGAIN ? gain() : g, type));
  if (!sed) return -999;
  double val = (sed->isInRange(sample) ? (xip ? sed->xip()(sample) : sed->xi()(sample)) : -999);
  return val;
}

◆ xiNorm()

double Data::xiNorm ( ShapeErrorType shapeErrorType = BestShapeError ) const

Definition at line 474 of file Data.cxx.

{
  if (!m_history) return -1;
  std::unique_ptr<const ShapeErrorData> sed (m_history->shapeErrorData(gain(), type));
  if (!sed) return -1;
  double val = TMath::Sqrt(sed->xi().Norm2Sqr());
  return val;
}

◆ xip()

double Data::xip	(	short	sample,
		ShapeErrorType	shapeErrorType = BestShapeError,
		CaloGain::CaloGain	g = CaloGain::UNKNOWNGAIN ) const

Definition at line 484 of file Data.cxx.

{
  return xi(sample, type, g, true);
}

◆ DataTweaker

friend class DataTweaker

friend

Definition at line 275 of file Data.h.

Member Data Documentation

◆ m_container

std::unique_ptr<const DataContainer> LArSamples::Data::m_container

private

Definition at line 283 of file Data.h.

◆ m_eventData

const EventData& LArSamples::Data::m_eventData

private

Definition at line 284 of file Data.h.

◆ m_history

const History* LArSamples::Data::m_history

mutableprivate

Definition at line 285 of file Data.h.

◆ m_index

unsigned int LArSamples::Data::m_index

mutableprivate

Definition at line 286 of file Data.h.

◆ m_timeShift

std::atomic<double> LArSamples::Data::m_timeShift {0}

inlinestaticprivate

Definition at line 287 of file Data.h.

287{0}; // specify a global time shift between first sample time and reported ofc time

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

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ Data() [1/3]

◆ Data() [2/3]

◆ Data() [3/3]

◆ ~Data()

Member Function Documentation

◆ _adcMax()

◆ _adcMax_new()

◆ _adjDeltaT()

◆ _adjScale()

◆ _bunchId()

◆ _calo()

◆ _channel()

◆ _chi2()

◆ _chi2_cellCorr()

◆ _chi2_k()

◆ _chi2_noCorr()

◆ _chi2_ringCorr()

◆ _chi2Anomaly()

◆ _count()

◆ _delta()

◆ _energy()

◆ _energySignif()

◆ _eta()

◆ _event()

◆ _feedThrough()

◆ _gain()

◆ _goodForShapeCorr()

◆ _hash()

◆ _iEta()

◆ _index()

◆ _iPhi()

◆ _layer()

◆ _lumiBlock()

◆ _maxPosition()

◆ _maxValue()

◆ _minValue()

◆ _nData()

◆ _noise()

◆ _normResidualError()

◆ _normResidualOffset()

◆ _null()

◆ _ofcGamma()

◆ _ofcSigma()

◆ _ofcTime()

◆ _pass()

◆ _peakSignif()

◆ _pedestal()

◆ _pedestalRMS()

◆ _pedestalSubstractedSample()

◆ _phi()

◆ _quality()

◆ _refitChi2()

◆ _refitDeltaT()

◆ _refitScale()

◆ _region()

◆ _resCorrN()

◆ _residual()

◆ _residualError()

◆ _residualOffset()

◆ _ring()

◆ _rt()

◆ _run()

◆ _sample()

◆ _slot()

◆ _timeECorr()

◆ _timeForSplash1()

◆ _timeForSplash2()

◆ _timeNoTOF()

◆ _upstreamE()

◆ _upstreamERatio()

◆ _x()

◆ _xi()

◆ _xiCell()