LArSamples::DigitMonitor Class Reference

storage of the time histories of all the cells More...

#include <DigitMonitor.h>

Inheritance diagram for LArSamples::DigitMonitor:

Collaboration diagram for LArSamples::DigitMonitor:

Public Member Functions
	DigitMonitor (const Interface &interface)
	Constructor.
TH1D *	chi2Dist (const TString &name, int nBins, double max, double kFactor=0, double fitMax=-1, int lwb=-1, int upb=-1, unsigned int chi2Pars=DefaultChi2, ShapeErrorType shapeErrorType=BestShapeError, unsigned int nDof=0) const
TH1D *	bestChi2Dist (const TString &name, int nBins, double max, double refErrMin, double refErrMax, unsigned int refErrNBins, int lwb=-1, int upb=-1, unsigned int chi2Pars=DefaultChi2, unsigned int nDof=0) const
TH1D *	gainDist (const TString &name) const
TH1D *	layerDist (const TString &name) const
TH1D *	energyDist (const TString &name, int nBins, double eMax=10000) const
TH1D *	timeDist (const TString &name, int nBins, double tMin=-25, double tMax=25) const
TH1D *	adcMaxDist (const TString &name, int nBins, double aMin=0, double aMax=4096) const
TH1D *	noiseDist (const TString &name, int nBins, double max) const
TH2D *	maxValueMap (TString name, PartitionId partition) const
TH2D *	minValueMap (TString name, PartitionId partition) const
TH1D *	residualDist (unsigned int k, const TString &name, int nBins, double rMin, double rMax, bool norm=false) const
bool	residualPlots (CaloId calo, unsigned int layer, CaloGain::CaloGain gain=CaloGain::LARHIGHGAIN, bool xip=false, bool ring=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
bool	residualPlotsRingComp (CaloId calo, unsigned int layer, CaloGain::CaloGain gain=CaloGain::LARHIGHGAIN, bool xip=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
bool	residualPlotsGainComp (CaloId calo, unsigned int layer, bool ring=false, bool xip=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
TH1D *	shapeErrorDist (unsigned int k, const TString &name, int nBins, double rMin, double rMax, double mean=0) const
double	residualCorr (unsigned int k1, unsigned int k2) const
bool	makeResidualCorrections (const TString &outputFile, short resTrunc=-1, short timeTrunc=-1, double absResTrunc=-1, unsigned int minSize=0, bool weigh=false, bool adjust=false, bool zeroTime=false) const
std::unique_ptr< Residuals >	getResiduals (unsigned int hash, CaloGain::CaloGain gain, double absResTrunc=-1, bool adjust=false, bool zeroTime=false) const
bool	residualParams (int lwb, int upb, CovMatrix &k, TVectorD &means) const
CovMatrix	kMatrix (int lwb, int upb)
TVectorD	means (int lwb, int upb)
int	combine (SimpleShape *&shape, SimpleShape *&ref, const TString &selection="", bool timeAligned=true) const
TH1D *	dist (const DataFuncSet &func, const DataFuncArgs &args, const TString &name, int nBins, double xMin, double xMax, const TString &title="", const TString &xTitle="", const TString &yTitle="", const FilterParams &f=FilterParams()) const
TH2D *	dist (const DataFuncSet &funcX, const DataFuncArgs &argsX, const DataFuncSet &funcY, const DataFuncArgs &argsY, const TString &name, int nBinsX, double xMin, double xMax, int nBinsY, double yMin, double yMax, const TString &title="", const TString &xTitle="", const TString &yTitle="", const FilterParams &f=FilterParams()) const
TH2D *	partitionMap (const DataFuncSet &func, const DataFuncArgs &args, TString name, PartitionId partition, const TString &title="", CombinationType comb=AverageValue, const FilterParams &f=FilterParams()) const
TH2D *	etaPhiMap (const DataFuncSet &func, const DataFuncArgs &args, const TString &name, CaloId calo, short layer, TString title="", CombinationType comb=AverageValue, const FilterParams &f=FilterParams()) const
bool	prepareDumpParams (const TString &vars, int verbosity, std::vector< TString > &variables, std::vector< DataFuncSet > &funcs, std::vector< DataFuncArgs > &args, std::vector< TString > &formats, TString &locFormat, TString &locHeader, TString &varHeader) const
bool	dump (const TString &vars, const FilterParams &f=FilterParams(), unsigned int verbosity=1) const
bool	dump (const TString &vars, CombinationType comb=AverageValue, const FilterParams &f=FilterParams(), const TString &ranges="", unsigned int verbosity=1) const
bool	statParams (const std::vector< DataFuncSet > &funcs, const std::vector< DataFuncArgs > &args, TVectorD &mean, TVectorD &meanErr, TMatrixD &covMatrix, TMatrixD &covMatrixErr, const FilterParams &f=FilterParams()) const
const Interface &	interface () const
virtual unsigned int	nChannels () const override
virtual std::unique_ptr< const History >	getCellHistory (unsigned int i) const override
virtual std::unique_ptr< const History >	getSCHistory (unsigned int i) const override
virtual std::unique_ptr< const CellInfo >	getCellInfo (unsigned int i) const override
virtual std::unique_ptr< const History >	newCellHistory (unsigned int i) const
virtual const History *	cellHistory (unsigned int i) const
virtual std::unique_ptr< const CellInfo >	cellInfo (unsigned int i) const
virtual unsigned int	nChannelsSC () const
virtual void	resetCache () const
const CellInfo *	cellInfoCache (unsigned int i) const
const History *	pass (unsigned int i, const FilterParams &f) const
const History *	cellCache () const
unsigned int	cachePos () const

Static Public Member Functions
static TF1 *	chi2Func (const char *name, double xMin, double xMax)
static TF1 *	fitChi2 (TH1D &h, const char *name, double xMin, double xMax, double nDof=-1, double fitMin=Definitions::none, double fitMax=Definitions::none)
static double	history_value (const History &history, const DataFuncSet &func, const DataFuncArgs &args, CombinationType comb, const FilterParams &f, unsigned int &nValues)
static DataFuncSet	func (const TString &var)
static bool	parseVariables (TString varStr, std::vector< TString > &vars, std::vector< DataFuncSet > &funcs, std::vector< DataFuncArgs > &args)
static TString	str (CombinationType comb)

Static Public Attributes
static const int	printPeriodicity = 10000

Private Attributes
const Interface *	m_interface
unsigned int	m_pos
std::unique_ptr< const History >	m_cellCache
std::vector< std::unique_ptr< CellInfo > >	m_cellInfoCache

Detailed Description

storage of the time histories of all the cells

Definition at line 28 of file DigitMonitor.h.

Constructor & Destructor Documentation

DigitMonitor()

LArSamples::DigitMonitor::DigitMonitor ( const Interface & interface )

inline

Constructor.

Definition at line 34 of file DigitMonitor.h.

: MonitorBase(interface) { } 

Member Function Documentation

adcMaxDist()

TH1D * DigitMonitor::adcMaxDist	(	const TString &	name,
		int	nBins,
		double	aMin = 0,
		double	aMax = 4096 ) const

Definition at line 47 of file DigitMonitor.cxx.

{ return dist(&Data::_adcMax, DataFuncArgs(), name, nBins, aMin, aMax, "ADCMax distribution", "ADCMax (counts)"); }

bestChi2Dist()

TH1D * DigitMonitor::bestChi2Dist	(	const TString &	name,
		int	nBins,
		double	max,
		double	refErrMin,
		double	refErrMax,
		unsigned int	refErrNBins,
		int	lwb = -1,
		int	upb = -1,
		unsigned int	chi2Pars = DefaultChi2,
		unsigned int	nDof = 0 ) const

Definition at line 125 of file DigitMonitor.cxx.

{
  TF1* fChi2 = chi2Func("fChi2", 0, max);
  if (nDof > 0) fChi2->FixParameter(1, nDof);
 
  double delta = (refErrMax - refErrMin)/refErrNBins;
  TH1D* fitResults = new TH1D("fitResults", "", refErrNBins, refErrMin - delta/2, refErrMax - delta/2);
  
  for (unsigned int k = 0; k < refErrNBins; k++) {
    double refErr = refErrMin + delta*k;
    TH1D* h = new TH1D(name, "#chi^{2} distribution", nBins, 0, max);
    UniformShapeErrorGetter kFactorGetter(refErr);
    CombinedShapeErrorGetter shapeErrorGetter;
    if (interface().shapeErrorGetter()) shapeErrorGetter.add(*interface().shapeErrorGetter());
    shapeErrorGetter.add(kFactorGetter);
    for (unsigned int i = 0; i < nChannels(); i++) {
      const History* history = cellHistory(i);
      if (!history) continue;
      history->setShapeErrorGetter(&shapeErrorGetter);
      for (unsigned int j = 0; j < history->nData(); j++) {
        h->Fill(history->chi2(j, lwb, upb, chi2Pars));
      }
    }
    h->Fit(fChi2); 
    fitResults->SetBinContent(k+1, fChi2->GetChisquare());
    delete h;
  }
      
  TF1* fPol2 = new TF1("fPol2", "[0]*(x - [1])*(x - [1]) + [2]", refErrMin, refErrMax);
  
  double center = fitResults->GetBinCenter(fitResults->GetMinimumBin());
  fPol2->SetParameter(0, (fitResults->GetMaximum() - fitResults->GetMinimum())/center/center);
  fPol2->SetParameter(1, center);
  fPol2->SetParameter(2, fitResults->GetMinimum());
  
  fitResults->Draw();
  fitResults->Fit(fPol2);
  
  TPaveText* p = new TPaveText(0.6, 0.85, 0.85, 0.65, "NDC");
  p->AddText(Form("Fit k_{best} = %.1f%%", fPol2->GetParameter(1)*100));
  p->AddText(Form("Fit #delta k_{best} = %.2f%%", 1/TMath::Sqrt(fPol2->GetParameter(0))*100));
  p->AddText(Form("Fit #chi^{2}_{min} = %.2f", fPol2->GetParameter(2)));
  p->Draw();
  
  return fitResults;
}

cachePos()

unsigned int LArSamples::AbsLArCells::cachePos ( ) const

inlineinherited

Definition at line 49 of file AbsLArCells.h.

{ return m_pos; }

cellCache()

const History * LArSamples::AbsLArCells::cellCache ( ) const

inlineinherited

Definition at line 48 of file AbsLArCells.h.

{ return m_cellCache.get(); }

cellHistory()

const History * AbsLArCells::cellHistory ( unsigned int i ) const

virtualinherited

Reimplemented in LArSamples::Interface.

Definition at line 54 of file AbsLArCells.cxx.

{ 
  if (m_pos == i) return m_cellCache.get();
  resetCache();
  std::unique_ptr<const History> history = newCellHistory(i);
  if (!history) return nullptr;
  m_cellCache = std::move(history);
  m_pos = i;
  return m_cellCache.get();
}

cellInfo()

std::unique_ptr< const CellInfo > AbsLArCells::cellInfo ( unsigned int i ) const

virtualinherited

Definition at line 66 of file AbsLArCells.cxx.

{
  const CellInfo* info = cellInfoCache(i);
  if (info) {
    if (info->isValid()) {
      return std::make_unique<CellInfo> (*info);
    }
    return nullptr;
  }
  std::unique_ptr<const CellInfo> infop = getCellInfo(i);
  if (infop)  m_cellInfoCache[i] =  std::make_unique<CellInfo>(*infop, false);
  return infop;
}

cellInfoCache()

const CellInfo * AbsLArCells::cellInfoCache ( unsigned int i ) const

inherited

Definition at line 81 of file AbsLArCells.cxx.

{
  return m_cellInfoCache[i].get();
}

chi2Dist()

TH1D * DigitMonitor::chi2Dist	(	const TString &	name,
		int	nBins,
		double	max,
		double	kFactor = 0,
		double	fitMax = -1,
		int	lwb = -1,
		int	upb = -1,
		unsigned int	chi2Pars = DefaultChi2,
		ShapeErrorType	shapeErrorType = BestShapeError,
		unsigned int	nDof = 0 ) const

Definition at line 60 of file DigitMonitor.cxx.

{
  TH1D* h = new TH1D(name, "#chi^{2} distribution", nBins, 0, max);
  unsigned int nDofEff = 0, nDofLast = 0;
  std::unique_ptr<UniformShapeErrorGetter> kFactorGetter;
  std::unique_ptr<CombinedShapeErrorGetter> shapeErrorGetter;
  if (kFactor > 0) {
    kFactorGetter = std::make_unique<UniformShapeErrorGetter>(kFactor);
    shapeErrorGetter = std::make_unique<CombinedShapeErrorGetter>();
    if (interface().shapeErrorGetter()) shapeErrorGetter->add(*interface().shapeErrorGetter());
    shapeErrorGetter->add(*kFactorGetter);
  }
  for (unsigned int i = 0; i < nChannels(); i++) {
    if ((i+1) % 10000 == 0) cout << "Processing entry # " << i+1 << endl;
    const History* history = cellHistory(i);
    if (!history) continue;
    if (shapeErrorGetter) history->setShapeErrorGetter(shapeErrorGetter.get()); // Do it at history level so it doesn't "stick" afterwards,,,
    for (unsigned int j = 0; j < history->nData(); j++) {
      h->Fill(history->chi2(j, lwb, upb, chi2Pars, shapeErrorType, &nDofEff));
      if (nDofLast > 0 && nDofLast != nDofEff) cout << "WARNING @ hash = " << i << ", index = " << j << " : nDof varied from " << nDofLast << " to " << nDofEff << endl;
      nDofLast = nDofEff;
    }
  }
  
  if (nDof > 0 && nDof != nDofEff) cout << "WARNING : fixing nDof = " << nDof << ", but computer nDofEff = " << nDofEff << endl;
  h->Draw();
  if (fitMax < 0) return h;
  TF1* fChi2 = fitChi2(*h, "chi2", 0, max, nDof, 0, fitMax);  
  TPaveText* p = new TPaveText(0.6, 0.85, 0.85, 0.65, "NDC");
  p->AddText(Form("k = %.1f%%", kFactor*100));
  p->AddText(Form("Fit #chi^{2} = %.1f", fChi2->GetChisquare()));
  p->AddText(Form("Fit nDOF = %.1f", fChi2->GetParameter(1)));
  p->Draw();
  return h;
}

chi2Func()

TF1 * DigitMonitor::chi2Func ( const char * name, double xMin, double xMax )

static

Definition at line 99 of file DigitMonitor.cxx.

{
  TF1* fChi2 = new TF1(name, "[0]*TMath::Power(x/2, [1]/2)/x*TMath::Exp(-x/2)/TMath::Gamma([1]/2)", xMin, xMax);
  fChi2->SetParameter(0, 100);
  fChi2->SetParameter(1, 3);
  fChi2->SetParLimits(1, 1, 7);
  return fChi2;
}

combine()

int DigitMonitor::combine	(	SimpleShape *&	shape,
		SimpleShape *&	ref,
		const TString &	selection = "",
		bool	timeAligned = true ) const

Definition at line 520 of file DigitMonitor.cxx.

{
  FilterParams f;
  if (!f.set(selection)) return 0;
  int n = 0;
  pshape = pref = nullptr;
  //don't change existing interface, use these internally
  std::unique_ptr<SimpleShape> shape;
  std::unique_ptr<SimpleShape> ref;
  double maxSum = 0;
  
  for (unsigned int i = 0; i < nChannels(); i++) {  
    if ((i+1) % 50000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
    const History* history = pass(i, f);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (!f.passEvent(*history->data(j))) continue;
      maxSum += history->data(j)->maxValue();
      cout << "Adding pulse (" << n+1 << ") at hash " << i << ", index " << j << ", max = " << maxSum/(n+1) << endl;
      //cout << i << " " << j << endl;
      auto thisData = std::make_unique<SimpleShape>(*history->data(j));
      auto thisRef =  history->referenceShape(j);
      if (timeAligned) {
        if (!SimpleShape::scaleAndShift(thisData, 1, -history->data(j)->ofcTime())) return -1;
        if (!SimpleShape::scaleAndShift(thisRef,  1, -history->data(j)->ofcTime())) return -1;
      }
      if (!SimpleShape::add(shape, *thisData)) return -1;
      if (!SimpleShape::add(ref, *thisRef)) return -1;      
      n++;
    }
  }
 
  if (n==0) return -1;
  if (!SimpleShape::scaleAndShift(shape, 1.0/n)) return -1;
  if (!SimpleShape::scaleAndShift(ref,   1.0/n)) return -1;
   //hand these back in existing interface
  pshape = shape.release();
  pref = ref.release();
  return n;
}

dist() [1/2]

TH1D * MonitorBase::dist ( const DataFuncSet & func, const DataFuncArgs & args, const TString & name, int nBins, double xMin, double xMax, const TString & title = "", const TString & xTitle = "", const TString & yTitle = "", const FilterParams & f = FilterParams() ) const

inherited

Definition at line 69 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  TH1D* h = new TH1D(name, title, nBins, xMin, xMax);
  if (xTitle != "") h->GetXaxis()->SetTitle(xTitle);
  if (yTitle != "") h->GetYaxis()->SetTitle(yTitle);
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    if ((i+1) % printPeriodicity == 0) cout << "Cell # " << i+1 << "/" << m_interface->nChannels() << endl;
    const History* history = pass(i, f);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (!f.passEvent(*history->data(j))) continue;
      h->Fill(func.val(*history->data(j), args));
    }
  }
  return h;
}

dist() [2/2]

TH2D * MonitorBase::dist ( const DataFuncSet & funcX, const DataFuncArgs & argsX, const DataFuncSet & funcY, const DataFuncArgs & argsY, const TString & name, int nBinsX, double xMin, double xMax, int nBinsY, double yMin, double yMax, const TString & title = "", const TString & xTitle = "", const TString & yTitle = "", const FilterParams & f = FilterParams() ) const

inherited

Definition at line 89 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  TH2D* h = new TH2D(name, title, nBinsX, xMin, xMax, nBinsY, yMin, yMax);
  if (xTitle != "") h->GetXaxis()->SetTitle(xTitle);
  if (yTitle != "") h->GetYaxis()->SetTitle(yTitle);
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    if ((i+1) % printPeriodicity == 0) cout << "Cell # " << i+1 << "/" << m_interface->nChannels() << endl;
    const History* history = pass(i, f);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (!f.passEvent(*history->data(j))) continue;
      h->Fill(funcX.val(*history->data(j), argsX), 
          funcY.val(*history->data(j), argsY));
    }
  }
  return h;
}

dump() [1/2]

bool MonitorBase::dump ( const TString & vars, CombinationType comb = AverageValue, const FilterParams & f = FilterParams(), const TString & ranges = "", unsigned int verbosity = 1 ) const

inherited

Definition at line 317 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  std::vector<TString> variables, formats;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  TString locFormat, locHeader, varHeader;
  if (!prepareDumpParams(vars, verbosity, variables, funcs, args,
                         formats, locFormat, locHeader, varHeader)) return false;
 
  std::vector<double> minVals, maxVals;
  TObjArray* rangeList = ranges.Tokenize(":");
  for (unsigned int k = 0; k < variables.size(); k++) {
    TString range = (k < (unsigned int)rangeList->GetEntries() ? ((TObjString*)rangeList->At(k))->String() : "");
    if (range == "") {
      minVals.push_back(-DBL_MAX);
      maxVals.push_back(+DBL_MAX);
      continue;
    }
    TObjArray* minMaxList = range.Tokenize(",");
    if (minMaxList->GetEntries() == 2) {
      minVals.push_back(((TObjString*)minMaxList->At(0))->String().Atof());
      maxVals.push_back(((TObjString*)minMaxList->At(1))->String().Atof());
    }
    else {
      cout << "Invalid range specification " << range << endl;
      return false;
    }
    delete minMaxList;
  }
  delete rangeList;
 
 TString header = locHeader + varHeader;
  TString separator = "";
  for (int i = 0; i < header.Length(); i++) separator += "-";
  cout << separator << endl << header << endl << separator << endl;
  
  unsigned int nEntries = 0, nValues = 0;
  
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    const History* history = pass(i, f);
    if (!history) continue;
    TString line = Form(locFormat, i, history->cellInfo()->location(verbosity).Data());
    bool ok = true;
    for (unsigned int k = 0; k < variables.size(); k++) {
      double val = history_value(*history, funcs[k], args[k], comb, f, nValues);
      if (nValues == 0 || val < minVals[k] || val >= maxVals[k]) { ok = false; break; }
      line += Form(formats[k].Data(), val);
    }
    if (!ok) continue;
    cout << line << endl;
    nEntries++;
  }
  cout << nEntries << " entries selected" << endl;
  return true;
}

dump() [2/2]

bool MonitorBase::dump ( const TString & vars, const FilterParams & f = FilterParams(), unsigned int verbosity = 1 ) const

inherited

Definition at line 284 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  std::vector<TString> variables, formats;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  TString locFormat, locHeader, varHeader;
  if (!prepareDumpParams(vars, verbosity, variables, funcs, args,
                         formats, locFormat, locHeader, varHeader)) return false;
  TString header = locHeader + " idx |" + varHeader;
  TString separator = "";
  for (int i = 0; i < header.Length(); i++) separator += "-";
  cout << separator << endl << header << endl << separator << endl;
  
  unsigned int nEntries = 0;
  
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    const History* history = pass(i, f);
    if (!history) continue;
    TString part1 = Form(locFormat, i, history->cellInfo()->location(verbosity).Data());
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (!f.passEvent(*history->data(j))) continue;
      TString line = part1 + Form(" %3d |", j);
      for (unsigned int k = 0; k < variables.size(); k++)
        line += Form(formats[k].Data(), funcs[k].val(*history->data(j), args[k]));
      cout << line << endl;
      nEntries++;
    }
  }
  cout << nEntries << " entries selected" << endl;
  return true;
}

energyDist()

TH1D * DigitMonitor::energyDist	(	const TString &	name,
		int	nBins,
		double	eMax = 10000 ) const

Definition at line 41 of file DigitMonitor.cxx.

{return dist(&Data::_energy, DataFuncArgs(), name, nBins, 0, eMax, "Energy distribution", "Energy (MeV)"); }

etaPhiMap()

TH2D * MonitorBase::etaPhiMap ( const DataFuncSet & func, const DataFuncArgs & args, const TString & name, CaloId calo, short layer, TString title = "", CombinationType comb = AverageValue, const FilterParams & f = FilterParams() ) const

inherited

Definition at line 131 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  if (title == "") title = Form("%s %s (%s, layer %d)", name.Data(), str(comb).Data(), Id::str(calo).Data(), layer);
  TH2D* h = Geo::etaPhiHist(calo, layer, name, title);
  FilterParams ff(f);
  ff.addCalo(calo);
  ff.addLayer(layer);
  unsigned int nValues = 0;
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    if ((i+1) % printPeriodicity == 0) cout << "Cell # " << i+1 << "/" << m_interface->nChannels() << endl;
    const History* history = pass(i, ff);
    if (!history) continue;
    h->Fill(history->cellInfo()->eta(), history->cellInfo()->phi(), 
        history_value(*history, func, args, comb, f, nValues));
  }
  return h;
}

fitChi2()

TF1 * DigitMonitor::fitChi2 ( TH1D & h, const char * name, double xMin, double xMax, double nDof = -1, double fitMin = Definitions::none, double fitMax = Definitions::none )

static

Definition at line 109 of file DigitMonitor.cxx.

{
  if (fitMin == Definitions::none) fitMin = xMin;
  if (fitMax == Definitions::none) fitMin = xMax;
 
  TF1* fChi2 = chi2Func(name, xMin, xMax);
  fChi2->SetParameter(0, 4*h.GetMaximum());
  if (nDof > 0) fChi2->FixParameter(1, nDof);
  
  h.Fit(fChi2, "", "", fitMin, fitMax);
  fChi2->SetRange(xMin, xMax); 
  return fChi2;
}

func()

DataFuncSet MonitorBase::func ( const TString & var )

staticinherited

Definition at line 431 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  if (var == "sample")         return DataFuncSet(&Data::_sample);
  if (var == "pedSubSample")   return DataFuncSet(&Data::_pedestalSubstractedSample);
  if (var == "energy")         return DataFuncSet(&Data::_energy);
  if (var == "adcMax")         return DataFuncSet(&Data::_adcMax);
  if (var == "adcMax_new")     return DataFuncSet(&Data::_adcMax_new);
  if (var == "gain")           return DataFuncSet(&Data::_gain);
  if (var == "peakSignif")     return DataFuncSet(&Data::_peakSignif);
  if (var == "energySignif")   return DataFuncSet(&Data::_energySignif);
  if (var == "ofcTime")        return DataFuncSet(&Data::_ofcTime);
  if (var == "sampleMax")      return DataFuncSet(&Data::_maxValue);
  if (var == "maxPosition")    return DataFuncSet(&Data::_maxPosition);
  if (var == "quality")        return DataFuncSet(&Data::_quality);
  if (var == "noise")          return DataFuncSet(&Data::_noise);
  if (var == "pedestal")       return DataFuncSet(&Data::_pedestal);
  if (var == "pedestalRMS")    return DataFuncSet(&Data::_pedestalRMS);
  if (var == "maxValue")       return DataFuncSet(&Data::_maxValue);
  if (var == "minValue")       return DataFuncSet(&Data::_minValue);
  if (var == "count")          return DataFuncSet(&Data::_count);
  if (var == "chi2")           return DataFuncSet(&Data::_chi2);
  if (var == "chi2_noCorr")    return DataFuncSet(&Data::_chi2_noCorr);
  if (var == "chi2_ringCorr")  return DataFuncSet(&Data::_chi2_ringCorr);
  if (var == "chi2_cellCorr")  return DataFuncSet(&Data::_chi2_cellCorr);
  if (var == "chi2_k")         return DataFuncSet(&Data::_chi2_k);
  if (var == "refitScale")     return DataFuncSet(&Data::_refitScale);
  if (var == "refitDeltaT")    return DataFuncSet(&Data::_refitDeltaT);
  if (var == "refitChi2")      return DataFuncSet(&Data::_refitChi2);
  if (var == "adjScale")       return DataFuncSet(&Data::_adjScale);
  if (var == "adjDeltaT")      return DataFuncSet(&Data::_adjDeltaT);
  if (var == "run")            return DataFuncSet(&Data::_run);
  if (var == "event")          return DataFuncSet(&Data::_event);
  if (var == "lumiBlock")      return DataFuncSet(&Data::_lumiBlock);
  if (var == "bunchId")        return DataFuncSet(&Data::_bunchId);
  if (var == "x")              return DataFuncSet(&Data::_x);
  if (var == "y")              return DataFuncSet(&Data::_y);
  if (var == "z")              return DataFuncSet(&Data::_z);
  if (var == "rt")             return DataFuncSet(&Data::_rt);
  if (var == "eta")            return DataFuncSet(&Data::_eta);
  if (var == "phi")            return DataFuncSet(&Data::_phi);
  if (var == "iEta")           return DataFuncSet(&Data::_iEta);
  if (var == "iPhi")           return DataFuncSet(&Data::_iPhi);
  if (var == "region")         return DataFuncSet(&Data::_region);
  if (var == "calo")           return DataFuncSet(&Data::_calo);
  if (var == "layer")          return DataFuncSet(&Data::_layer);
  if (var == "feedthrough")    return DataFuncSet(&Data::_feedThrough);
  if (var == "FT")             return DataFuncSet(&Data::_feedThrough);
  if (var == "slot")           return DataFuncSet(&Data::_slot);
  if (var == "channel")        return DataFuncSet(&Data::_channel);
  if (var == "ring")           return DataFuncSet(&Data::_ring);
  if (var == "index")          return DataFuncSet(&Data::_index);
  if (var == "hash")           return DataFuncSet(&Data::_hash);
  if (var == "nData")          return DataFuncSet(&Data::_nData);
  if (var == "timeNoTOF")      return DataFuncSet(&Data::_timeNoTOF);
  if (var == "timeForSplash1") return DataFuncSet(&Data::_timeForSplash1);
  if (var == "timeForSplash2") return DataFuncSet(&Data::_timeForSplash2);
  if (var == "delta")          return DataFuncSet(&Data::_delta);
  if (var == "residual")       return DataFuncSet(&Data::_residual);
  if (var == "residualOffset") return DataFuncSet(&Data::_residualOffset);
  if (var == "resCorrN")       return DataFuncSet(&Data::_resCorrN);
  if (var == "xi")             return DataFuncSet(&Data::_xi);
  if (var == "xiCell")         return DataFuncSet(&Data::_xiCell);
  if (var == "xiRing")         return DataFuncSet(&Data::_xiRing);
  if (var == "xiNormCell")     return DataFuncSet(&Data::_xiNormCell);
  if (var == "xiNormRing")     return DataFuncSet(&Data::_xiNormRing);
  if (var == "xip")            return DataFuncSet(&Data::_xip);
  if (var == "residualError")  return DataFuncSet(&Data::_residualError);
  if (var == "ofcSigma")       return DataFuncSet(&Data::_ofcSigma);
  if (var == "ofcGamma")       return DataFuncSet(&Data::_ofcGamma);
  if (var == "normResidualOffset") return DataFuncSet(&Data::_normResidualOffset);
  if (var == "normResidualError")  return DataFuncSet(&Data::_normResidualError);
  if (var == "pass")               return DataFuncSet(&Data::_pass);
  if (var == "goodForShapeCorr")   return DataFuncSet(&Data::_goodForShapeCorr);
  if (var == "upstreamE")          return DataFuncSet(&Data::_upstreamE);
  if (var == "upstreamERatio")     return DataFuncSet(&Data::_upstreamERatio);
  if (var == "chi2Anomaly")        return DataFuncSet(&Data::_chi2Anomaly);
  if (var == "timeECorr")      return DataFuncSet(&Data::_timeECorr);
 
  return DataFuncSet();
}

gainDist()

TH1D * DigitMonitor::gainDist

(

const TString &

name

)

const

Definition at line 176 of file DigitMonitor.cxx.

{
  TH1D* h = new TH1D(name, "gain distribution", 3, -0.5, 2.5);
  h->GetXaxis()->SetBinLabel(1, "Low");
  h->GetXaxis()->SetBinLabel(2, "Medium");
  h->GetXaxis()->SetBinLabel(3, "High");
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) h->Fill(history->data(j)->gain());
  }
  return h;
}

getCellHistory()

std::unique_ptr< const History > MonitorBase::getCellHistory ( unsigned int i ) const

overridevirtualinherited

Implements LArSamples::AbsLArCells.

Definition at line 34 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{ 
  return interface().getCellHistory(i); 
}

getCellInfo()

std::unique_ptr< const CellInfo > MonitorBase::getCellInfo ( unsigned int i ) const

overridevirtualinherited

Reimplemented from LArSamples::AbsLArCells.

Definition at line 44 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{ 
  return interface().getCellInfo(i);
}

getResiduals()

std::unique_ptr< Residuals > DigitMonitor::getResiduals	(	unsigned int	hash,
		CaloGain::CaloGain	gain,
		double	absResTrunc = -1,
		bool	adjust = false,
		bool	zeroTime = false ) const

Definition at line 561 of file DigitMonitor.cxx.

{
  const History* history = cellHistory(hash);
  if (!history) return nullptr;
  std::unique_ptr<const History> historyptr;
  if (adjust) {
    historyptr = history->adjust();
    history = historyptr.get();
  }
  return history->residuals(gain, absResTrunc, false, zeroTime);
}

getSCHistory()

std::unique_ptr< const History > MonitorBase::getSCHistory ( unsigned int i ) const

overridevirtualinherited

Implements LArSamples::AbsLArCells.

Definition at line 39 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{ 
  return interface().getSCHistory(i); 
}

history_value()

double MonitorBase::history_value ( const History & history, const DataFuncSet & func, const DataFuncArgs & args, CombinationType comb, const FilterParams & f, unsigned int & nValues )

staticinherited

Definition at line 404 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  double val = 0, sum1 = 0, sum2 = 0;
  nValues = 0;
  if (comb == MaxValue) val = DBL_MIN;
  if (comb == MinValue) val = DBL_MAX;
  for (unsigned int j = 0; j < history.nData(); j++) {
    if (!f.passEvent(*history.data(j))) continue;
    double thisVal = func.val(*history.data(j), args);
    nValues++;
    switch (comb) {
      case AverageValue : val += thisVal; break;
      case MaxValue     : if (thisVal > val) val = thisVal; break;
      case MinValue     : if (thisVal < val) val = thisVal; break;
      case TotalValue   : val += thisVal; break;
      case RMSValue     : sum1 += thisVal; sum2 += thisVal*thisVal; break; 
    }
  }
  if (comb == AverageValue && nValues > 0) val /= nValues;
  if (comb == RMSValue) 
    val = (nValues > 0 ? TMath::Sqrt((sum2 - sum1*sum1/nValues)/nValues) : -1);
  if (TMath::Abs(val) > 0.9*DBL_MAX) val = 0;
  return val;
}

interface()

const Interface & LArSamples::MonitorBase::interface ( ) const

inlineinherited

Definition at line 74 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h.

{ return *m_interface; }

kMatrix()

CovMatrix LArSamples::DigitMonitor::kMatrix ( int lwb, int upb )

inline

Definition at line 78 of file DigitMonitor.h.

{ CovMatrix k; TVectorD means; residualParams(lwb, upb, k, means); return k; }

layerDist()

TH1D * DigitMonitor::layerDist

(

const TString &

name

)

const

Definition at line 191 of file DigitMonitor.cxx.

{
  TH1D* h = new TH1D(name, "layer distribution", 4, -0.5, 3.5);
  h->GetXaxis()->SetBinLabel(1, "PS");
  h->GetXaxis()->SetBinLabel(2, "Strips");
  h->GetXaxis()->SetBinLabel(3, "Middle");
  h->GetXaxis()->SetBinLabel(4, "Back");
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) h->Fill(history->cellInfo()->layer());
  }
  return h;
}

makeResidualCorrections()

bool DigitMonitor::makeResidualCorrections	(	const TString &	outputFile,
		short	resTrunc = -1,
		short	timeTrunc = -1,
		double	absResTrunc = -1,
		unsigned int	minSize = 0,
		bool	weigh = false,
		bool	adjust = false,
		bool	zeroTime = false ) const

Definition at line 574 of file DigitMonitor.cxx.

{
  TreeShapeErrorGetter* shapeError = new TreeShapeErrorGetter(outputFile, true);
  std::vector< std::vector<ResidualCalculator> > ringCalcs;
  std::vector<ResidualCalculator> gainRingCalcs(Geo::nPhiRings());
  for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) ringCalcs.push_back(gainRingCalcs);
  
  cout << "Processing cells" << endl;
  
  for (unsigned int i = 0; i < nChannels(); i++) {
    if ((i+1) % 10000 == 0) cout << "Processing hash = " << (i+1) << endl;
    for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) {
      std::unique_ptr<Residuals> residuals = getResiduals(i, (CaloGain::CaloGain)g, absResTrunc, adjust, zeroTime);
      if (residuals && (resTrunc > 0 || timeTrunc > 0)) {
        residuals = residuals->truncate(resTrunc, timeTrunc);
      }
 
      if (residuals && residuals->size() < minSize) {
        residuals = nullptr;
      }
 
      ResidualCalculator* resCalc = nullptr;
      if (residuals) {
        resCalc = residuals->calculator(weigh);
      }     
 
      if (!resCalc) {
        shapeError->addCell(ResidualCalculator(), (CaloGain::CaloGain)g);
        continue;
      }
 
      cout << i << " : Phi ring " << cellHistory(i)->cellInfo()->globalPhiRing() << " : adding " << resCalc->size() << endl;
      shapeError->addCell(*resCalc, (CaloGain::CaloGain)g);
      short iring = cellHistory(i)->cellInfo()->globalPhiRing();
      if (iring >= 0) {
        ringCalcs[g][iring].append(*resCalc);
        if (iring == 1142) cout << ringCalcs[g][iring].regresser()->mean(0) << endl;
      }
      delete resCalc;
    }
  }
  
  cout << "Processing phi-symmetric data" << endl;
  
  for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) {
    for (short i = 0; i < Geo::nPhiRings(); i++) {
      if (ringCalcs[g][i].size() > 0) 
        cout << "Saving data for gain " << g << ", phi ring " << i << ", n = " << ringCalcs[g][i].size() << ", range = " << ringCalcs[g][i].rangeStr() << endl;
      shapeError->addRing(ringCalcs[g][i], (CaloGain::CaloGain)g);
    }
  }
  
  cout << "Done!" << endl;
  
  delete shapeError;
  return true;
}

maxValueMap()

TH2D * DigitMonitor::maxValueMap	(	TString	name,
		PartitionId	partition ) const

Definition at line 53 of file DigitMonitor.cxx.

{ return partitionMap(&Data::_maxValue, DataFuncArgs(), std::move(name), partition, "Max sample value (ADC counts)", MaxValue); }

means()

TVectorD LArSamples::DigitMonitor::means ( int lwb, int upb )

inline

Definition at line 79 of file DigitMonitor.h.

{ CovMatrix k; TVectorD means; residualParams(lwb, upb, k, means); return means; }

minValueMap()

TH2D * DigitMonitor::minValueMap	(	TString	name,
		PartitionId	partition ) const

Definition at line 56 of file DigitMonitor.cxx.

{ return partitionMap(&Data::_minValue, DataFuncArgs(), std::move(name), partition, "Min sample value (ADC counts)", MinValue); }

nChannels()

unsigned int MonitorBase::nChannels ( ) const

overridevirtualinherited

Reimplemented from LArSamples::AbsLArCells.

Definition at line 50 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{ 
  return interface().nChannels(); 
}

nChannelsSC()

virtual unsigned int LArSamples::AbsLArCells::nChannelsSC ( ) const

inlinevirtualinherited

Definition at line 36 of file AbsLArCells.h.

{ return Definitions::nChannelsSC; }

newCellHistory()

std::unique_ptr< const History > AbsLArCells::newCellHistory ( unsigned int i ) const

virtualinherited

Definition at line 39 of file AbsLArCells.cxx.

{ 
  std::unique_ptr<const History> history = getCellHistory(i);
  if (!history) return nullptr;
  if (!m_cellInfoCache[i]) {
    const CellInfo* ci=history->cellInfo();
    if (ci) {
      m_cellInfoCache[i]=std::make_unique<CellInfo>(*ci,false);
    }
  }
  //  m_cellInfoCache[i] = (history->cellInfo() ? new CellInfo(*history->cellInfo(), false) : new CellInfo());
  return history;
}

noiseDist()

TH1D * DigitMonitor::noiseDist	(	const TString &	name,
		int	nBins,
		double	max ) const

Definition at line 50 of file DigitMonitor.cxx.

{ return dist(&Data::_noise, DataFuncArgs(), name, nBins, 0, max, "Noise distribution", "Noise (ADC counts)"); }

parseVariables()

bool MonitorBase::parseVariables ( TString varStr, std::vector< TString > & vars, std::vector< DataFuncSet > & funcs, std::vector< DataFuncArgs > & args )

staticinherited

Definition at line 152 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  varStr.ReplaceAll(" ", "");
  TObjArray* varList = varStr.Tokenize(":");
  for (int i = 0; i < varList->GetEntries(); i++) {
    TString varAndSample = ((TObjString*)varList->At(i))->String();
    TString var = varAndSample;
    short sample1 = -1, sample2 = -1;
    double par = Definitions::none;
    TString str = "";
    if (varAndSample.CountChar('[') + varAndSample.CountChar(']') > 0) {
      if (varAndSample.CountChar('[') != 1 || varAndSample.CountChar(']') != 1 || varAndSample[varAndSample.Length() - 1] != ']') { 
    cout << "Invalid use of [ ] delimiters" << endl;
    delete varList; 
    return false; 
      }
      TObjArray* varSampList = varAndSample.Tokenize("[,]");
      if (varSampList->GetEntries() == 2) { 
        var = ((TObjString*)varSampList->At(0))->String();
        TString sampStr = ((TObjString*)varSampList->At(1))->String();
        delete varSampList;
        if (!sampStr.IsDigit()) { 
          cout << "Sample argument [sample] should be an integer" << endl;        
          delete varList; 
          return false; 
        }
        sample1 = sampStr.Atoi();
      }
      else if (varSampList->GetEntries() == 3) { 
        var = ((TObjString*)varSampList->At(0))->String();
        TString sampStr1 = ((TObjString*)varSampList->At(1))->String();
        TString sampStr2 = ((TObjString*)varSampList->At(2))->String();
        delete varSampList;
        if (!sampStr1.IsDigit() || !sampStr2.IsDigit()) { 
          cout << "Sample arguments [sample1, sample2] should be an integer" << endl;
          delete varList; 
          return false; 
        }
        sample1 = sampStr1.Atoi();
        sample2 = sampStr2.Atoi();
      }
      else {
    cout << "Invalid specification of sample argument (expected var[sample])" << endl;  
    delete varSampList; 
    delete varList; 
    return false; 
      }
    }
    if (varAndSample.CountChar('(') + varAndSample.CountChar(')') > 0) {
      if (varAndSample.CountChar('(') != 1 || varAndSample.CountChar(')') != 1 || varAndSample[varAndSample.Length() - 1] != ')') { 
        cout << "Invalid use of ( ) delimiters" << endl;
        delete varList; 
        return false; 
      }
      TObjArray* varSampList = varAndSample.Tokenize("(,)");
      if (varSampList->GetEntries() == 2) { 
        var = ((TObjString*)varSampList->At(0))->String();
        TString sampStr = ((TObjString*)varSampList->At(1))->String();
        delete varSampList;
        if (!sampStr.IsFloat()) { 
          cout << "floating-point argument expected between ()" << endl;        
          delete varList; 
          return false; 
        }
        par = sampStr.Atof();
      }
      else {
        cout << "Invalid specification of floating-point argument (expected func(arg))" << endl;      
        delete varSampList; 
        delete varList; 
        return false; 
      }
    }
    if (varAndSample(0, 5) == "pass_") {
      str = varAndSample(5, varAndSample.Length() - 5);
      var = "pass";
    }
    DataFuncSet fcn = func(var);
    if (fcn.isNull()) {
      cout << "Unknown variable " << var << endl;
      return false;
    }
    vars.push_back(varAndSample);
    funcs.push_back(fcn);
    args.emplace_back(sample1, sample2, par, Definitions::none, str);
  }
  delete varList;
  return true;
}

partitionMap()

TH2D * MonitorBase::partitionMap ( const DataFuncSet & func, const DataFuncArgs & args, TString name, PartitionId partition, const TString & title = "", CombinationType comb = AverageValue, const FilterParams & f = FilterParams() ) const

inherited

Definition at line 113 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  TH2D* h = Geo::partitionHist(partition, std::move(name), title + " " + Id::str(partition));
  unsigned int nValues = 0;
  FilterParams ff(f);
  ff.addPartition(partition);
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    if ((i+1) % printPeriodicity == 0) cout << "Cell # " << i+1 << "/" << m_interface->nChannels() << endl;
    const History* history = pass(i, ff);
    if (!history) continue;
    h->SetBinContent(history->cellInfo()->feb() + 1, history->cellInfo()->channel() + 1, 
             history_value(*history, func, args, comb, f, nValues));
  }
  return h;
}

pass()

const History * AbsLArCells::pass ( unsigned int i, const FilterParams & f ) const

inherited

Definition at line 96 of file AbsLArCells.cxx.

{ 
  //std::cout << "Called AbsLArCells with hash " << i  << std::endl;
  if (!f.passHash(i)) return nullptr;
  std::unique_ptr<const CellInfo> info = cellInfo(i);
  if (!info) {
    return nullptr;
  }
  //std::cout << "Called AbsLArCells::pass on a cell belonging to " << Id::str(info->calo()) << std::endl;
  bool result = f.passCell(*info);
  return result ? cellHistory(i) : nullptr;
}

prepareDumpParams()

bool MonitorBase::prepareDumpParams ( const TString & vars, int verbosity, std::vector< TString > & variables, std::vector< DataFuncSet > & funcs, std::vector< DataFuncArgs > & args, std::vector< TString > & formats, TString & locFormat, TString & locHeader, TString & varHeader ) const

inherited

Definition at line 244 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  if (!parseVariables(vars, variables, funcs, args)) {
    cout << "Invalid variable specification " << vars << endl;
    return false;
  }
 
  locFormat = "| %6d | %-";
  if (verbosity <= 1)
    locFormat += "26";
  else if (verbosity == 2)
    locFormat += "80";
  else if (verbosity == 3)
    locFormat += "108";
  else
    locFormat += "120";
  locFormat += "s |";
  
  TString locFormatHeader = locFormat; locFormatHeader.ReplaceAll("d", "s");
  locHeader = Form(locFormatHeader, "hash", "location");
    
  for (unsigned int i = 0; i < variables.size(); i++) {
    formats.emplace_back(" %-9.8g |");
    TString fH = formats[i]; fH.ReplaceAll("g", "s");
    TString varName = variables[i];
    if (args[i].i2 != -1) varName += Form("[%d, %d]", args[i].i1, args[i].i2);
    else if (args[i].i1 != -1) varName += Form("[%d]", args[i].i1);
    if (!Definitions::isNone(args[i].x1)) varName += Form("[%.5f]", args[i].x1);
    if (args[i].str != "") varName += Form("_%s", args[i].str.Data());
    varHeader += Form(fH, varName.Data());
  }
  return true;
}

resetCache()

void AbsLArCells::resetCache ( ) const

virtualinherited

Definition at line 32 of file AbsLArCells.cxx.

{
  m_cellCache.reset();
  m_pos = nChannels() + 1;
}

residualCorr()

double DigitMonitor::residualCorr	(	unsigned int	k1,
		unsigned int	k2 ) const

Definition at line 446 of file DigitMonitor.cxx.

{
  double sum = 0;
  int n = 0;
  
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      double residual1 = history->data(j)->delta(k1);
      double residual2 = history->data(j)->delta(k2);
      sum += residual1*residual2;
      n++;
    }  
  }
  
  return n==0 ? 0 : sum/n;
}

residualDist()

TH1D * DigitMonitor::residualDist	(	unsigned int	k,
		const TString &	name,
		int	nBins,
		double	rMin,
		double	rMax,
		bool	norm = false ) const

Definition at line 207 of file DigitMonitor.cxx.

{
  TH1D* h = new TH1D(name, "Residual distribution", nBins, rMin, rMax);
  h->GetXaxis()->SetTitle(norm ? "Residual/ADCMax (%)" : "Residual (ADC counts)");
    
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) 
      h->Fill(history->data(j)->delta(k)*(norm ? 100/history->data(j)->adcMax() : 1));
  }
  return h;
}

residualParams()

bool DigitMonitor::residualParams	(	int	lwb,
		int	upb,
		CovMatrix &	k,
		TVectorD &	means ) const

Definition at line 466 of file DigitMonitor.cxx.

{
  CovMatrix sum2N(lwb, upb), sumCN(lwb, upb);
  TVectorD sum(lwb, upb), sumN(lwb, upb);
  double sumA = 0, n = 0;
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (history->data(j)->adcMax() < 0) continue;
      std::unique_ptr<OFC> ofc = history->ofc(j);
      if (!ofc || ofc->hasSameRange(lwb, upb)) { cout << "Invalid index bounds!" << endl; return false; }
      for (int i1 = lwb; i1 <= upb; i1++) {
        sum(i1)   = sum(i1)  + history->data(j)->delta(i1);
        sumN(i1)  = sumN(i1) + history->data(j)->delta(i1)/TMath::Power(history->data(j)->adcMax(), 2);
        for (int i2 = lwb; i2 <= upb; i2++) {
          sumCN(i1, i2) = sumCN(i1, i2) + ofc->Gamma()(i1, i2)/TMath::Power(history->data(j)->adcMax(), 2);
          sum2N(i1, i2) = sum2N(i1, i2) + history->data(j)->delta(i1)*history->data(j)->delta(i2)/TMath::Power(history->data(j)->adcMax(),2);
        }
      }
      sumA +=  1/TMath::Power(history->data(j)->adcMax(), 2);
      n++;
    }
  }
 
  if (n > 0) {
    sumA /= n;
  }
  k.ResizeTo(lwb, upb);
  means.ResizeTo(lwb, upb);
  
  for (int i1 = lwb; i1 <= upb; i1++) {
    if (n > 0) {
      sum(i1)   = sum(i1)/n;
      sumN(i1)  = sumN(i1)/n;
    }
    means(i1) = sum(i1);
  }
 
  for (int i1 = lwb; i1 <= upb; i1++) {
    for (int i2 = lwb; i2 <= upb; i2++) {
      if (n != 0) {
        sumCN(i1, i2) = sumCN(i1, i2)/n;
        sum2N(i1, i2) = sum2N(i1, i2)/n;
      }
      k(i1, i2) = sum2N(i1, i2) - sumN(i1)*sum(i2) - sum(i1)*sumN(i2)
                + sum(i1)*sum(i2)*sumA - sumCN(i1,i2);
      k(i1, i2) = (k(i1, i2) >= 0 ? sqrt(k(i1, i2)) : -sqrt(-k(i1, i2)));
    }
  }
  return true;
}

residualPlots()

bool DigitMonitor::residualPlots	(	CaloId	calo,
		unsigned int	layer,
		CaloGain::CaloGain	gain = CaloGain::LARHIGHGAIN,
		bool	xip = false,
		bool	ring = false,
		double	lo = -0.05,
		double	hi = 0.05,
		const TString &	fileName = "residuals" )

Definition at line 222 of file DigitMonitor.cxx.

{
  TH1D* h1[5];
  TH2D* hEtaPhi[5];
  TString globalTitle =  TString(xip ? "Derivative" : "Offset") + " correction";
  for (unsigned short k = 0; k < 5; k++) {
    TString title =  globalTitle + " for sampling " + Form("%d", k);
    h1[k] = new TH1D(Form("h%d", k), title, 100, lo, hi);
    h1[k]->GetXaxis()->SetTitle("correction");
    hEtaPhi[k] = Geo::etaPhiHist(calo, layer, Form("hEtaPhi%d", k), title);
    hEtaPhi[k]->SetMinimum(lo);
    hEtaPhi[k]->SetMaximum(hi);
  }
  FilterParams f;
  f.addCalo(calo);
  f.addLayer(layer);
  TH2D* hAll = new TH2D("hAll", globalTitle + " vs. sampling", 5, 0, 5, 100, lo, hi);
  hAll->GetXaxis()->SetTitle("sampling");
  hAll->GetYaxis()->SetTitle("correction");    
  for (unsigned int i = 0; i < nChannels(); i++) {
    if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
    const History* history = pass(i, f);
    if (!history || history->nData() == 0) continue;
    for (unsigned short k = 0; k < 5; k++) {
      ShapeErrorType type = (ring ? RingShapeError : CellShapeError);
      double xi = history->data(0)->xi(k, type, gain, xip);
      h1[k]->Fill(xi);
      hEtaPhi[k]->Fill(history->cellInfo()->eta(), history->cellInfo()->phi(), xi);
      if (Id::matchCalo(calo, HEC) && history->cellInfo()->region() == 1) // For HEC, take care of larger phi size in region 1
        hEtaPhi[k]->Fill(history->cellInfo()->eta(), history->cellInfo()->phi() - Geo::phiSize(HEC_A, 1, 1), xi);      
    }
  }
  for (unsigned int i = 1; i <= 100; i++) 
    for (unsigned short k = 0; k < 5; k++) 
      hAll->SetBinContent(k+1, i, h1[k]->GetBinContent(i));
  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
  for (unsigned short k = 0; k < 5; k++) {
    h1[k]->Write();
    hEtaPhi[k]->Write();
  }
  hAll->Write();
  ff->Close();
  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 2000);
  c1->Divide(2,3);
  for (unsigned short k = 0; k < 5; k++) {
    c1->cd(k+1); 
    hEtaPhi[k]->Draw("COLZ");
  }
  c1->cd(6);
  hAll->Draw();
  c1->Print(fileName + ".png");
  c1->Print(fileName + ".eps");
  return true;
}

residualPlotsGainComp()

bool DigitMonitor::residualPlotsGainComp	(	CaloId	calo,
		unsigned int	layer,
		bool	ring = false,
		bool	xip = false,
		double	lo = -0.05,
		double	hi = 0.05,
		const TString &	fileName = "residuals" )

Definition at line 350 of file DigitMonitor.cxx.

{
  unsigned int nBins = 100;
  TH1D* h[5][5];
  TString globalTitle =  TString(xip ? "Derivative" : "Offset") + " correction";
  for (unsigned short l = 0; l < 5; l++)
    for (unsigned short k = 0; k < 5; k++)
      h[l][k] = new TH1D(Form("h%d%d", l, k), "", nBins, l < 2 ? -15 : lo, l < 2 ? 15 : hi);
  FilterParams f;
  f.addCalo(calo);
  f.addLayer(layer);
  TH2D* hAll[5];
  hAll[0] = new TH2D("h0All", "#xi_{high}/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
  hAll[1] = new TH2D("hDAll", "(#xi_{high} - #xi_{low})/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
  hAll[2] = new TH2D("hEAll", globalTitle + " higher-gain value vs. sampling", 5, 0, 5, nBins, 0, hi);
  hAll[3] = new TH2D("hRAll", globalTitle + " lower-gain value vs. sampling", 5, 0, 5, nBins, lo, hi);
  hAll[4] = new TH2D("hCAll", globalTitle + " higher-gain error vs. sampling", 5, 0, 5, nBins, lo, hi);
  TH2D* hCor = new TH2D("hCor", globalTitle + " lower-gain vs. higher-gain values", nBins, lo, hi, nBins, lo, hi);
  hCor->GetXaxis()->SetTitle("#xi_{high}");
  hCor->GetYaxis()->SetTitle("#xi_{low}");    
  for (unsigned short l = 0; l < 4; l++) {
    hAll[l]->GetXaxis()->SetTitle("sampling");
    hAll[l]->GetYaxis()->SetTitle("significance");    
  }
  
  CaloGain::CaloGain hiGain = (Id::matchCalo(calo, HEC) || Id::matchCalo(calo, FCAL) ? CaloGain::LARMEDIUMGAIN : CaloGain::LARHIGHGAIN);
  CaloGain::CaloGain loGain = (Id::matchCalo(calo, HEC) || Id::matchCalo(calo, FCAL) ? CaloGain::LARLOWGAIN : CaloGain::LARMEDIUMGAIN);
 
  for (unsigned int i = 0; i < nChannels(); i++) {
    if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
    const History* history = pass(i, f);
    if (!history || history->nData() == 0) continue;
    for (unsigned short k = 0; k < 5; k++) {
      std::unique_ptr<const ShapeErrorData> hgSED = history->shapeErrorData(hiGain, ring ? RingShapeError : CellShapeError);
      std::unique_ptr<const ShapeErrorData> lgSED = history->shapeErrorData(loGain, ring ? RingShapeError : CellShapeError);
      if (!hgSED || !lgSED) {
        continue;
      }
      double hgVal = (xip ? hgSED->xi()(k) : hgSED->xip()(k));
      double lgVal = (xip ? lgSED->xi()(k) : lgSED->xip()(k));
      double hgErr = TMath::Sqrt(xip ? hgSED->xiErr()(k,k) : hgSED->xipErr()(k,k));
      double error = TMath::Sqrt(xip ? hgSED->xiErr()(k,k) + lgSED->xiErr()(k,k) : hgSED->xipErr()(k,k) + lgSED->xipErr()(k,k));
      h[0][k]->Fill(hgVal/hgErr);
      h[1][k]->Fill((hgVal - lgVal)/error);
      h[2][k]->Fill(hgVal);
      h[3][k]->Fill(lgVal);
      h[4][k]->Fill(hgErr);
      hCor->Fill(hgVal, lgVal);
    }
  }
  for (unsigned int i = 1; i <= nBins; i++) 
    for (unsigned short l = 0; l < 5; l++)
      for (unsigned short k = 0; k < 5; k++)
        hAll[l]->SetBinContent(k+1, i, h[l][k]->GetBinContent(i));
  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
  for (unsigned short l = 0; l < 5; l++) {
    for (unsigned short k = 0; k < 5; k++)
      h[l][k]->Write();
    hAll[l]->Write();
  }
  ff->Close();
  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 500);
  c1->Divide(3,1);
  c1->cd(1); hCor->Draw();
  TLatex l; l.SetNDC();
  TString corrStr = Form("Correlation : %.1f%%", 100*hCor->GetCorrelationFactor());
  l.DrawLatex(0.2, 0.8, corrStr);
  c1->cd(2); hAll[0]->Draw("COL");
  c1->cd(3); hAll[1]->Draw("COL");
  c1->Print(fileName + ".png");
  c1->Print(fileName + ".eps");
  return true;
}

residualPlotsRingComp()

bool DigitMonitor::residualPlotsRingComp	(	CaloId	calo,
		unsigned int	layer,
		CaloGain::CaloGain	gain = CaloGain::LARHIGHGAIN,
		bool	xip = false,
		double	lo = -0.05,
		double	hi = 0.05,
		const TString &	fileName = "residuals" )

Definition at line 279 of file DigitMonitor.cxx.

{
  unsigned int nBins = 100;
  TH1D* h[5][5];
  TString globalTitle =  TString(xip ? "Derivative" : "Offset") + " correction";
  for (unsigned short l = 0; l < 5; l++)
    for (unsigned short k = 0; k < 5; k++)
      h[l][k] = new TH1D(Form("h%d%d", l, k), "", nBins, l < 2 ? -15 : lo, l < 2 ? 15 : hi);
  FilterParams f;
  f.addCalo(calo);
  f.addLayer(layer);
  TH2D* hAll[5];
  hAll[0] = new TH2D("h0All", "#xi_{cell}/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
  hAll[1] = new TH2D("hDAll", "(#xi_{cell} - #xi_{ring})/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
  hAll[2] = new TH2D("hEAll", globalTitle + " per-cell value vs. sampling", 5, 0, 5, nBins, 0, hi);
  hAll[3] = new TH2D("hRAll", globalTitle + " per-ring value vs. sampling", 5, 0, 5, nBins, lo, hi);
  hAll[4] = new TH2D("hCAll", globalTitle + " per-cell error vs. sampling", 5, 0, 5, nBins, lo, hi);
  TH2D* hCor = new TH2D("hCor", globalTitle + " per-ring vs. per-cell values", nBins, lo, hi, nBins, lo, hi);
  hCor->GetXaxis()->SetTitle("#xi_{cell}");
  hCor->GetYaxis()->SetTitle("#xi_{ring}");    
  for (unsigned short l = 0; l < 4; l++) {
    hAll[l]->GetXaxis()->SetTitle("sampling");
    hAll[l]->GetYaxis()->SetTitle("significance");    
  }
  for (unsigned int i = 0; i < nChannels(); i++) {
    if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
    const History* history = pass(i, f);
    if (!history || history->nData() == 0) continue;
    for (unsigned short k = 0; k < 5; k++) {
      std::unique_ptr<const ShapeErrorData> cellSED = history->shapeErrorData(gain, CellShapeError);
      std::unique_ptr<const ShapeErrorData> ringSED = history->shapeErrorData(gain, RingShapeError);
      if (!cellSED || !ringSED) {
        continue;
      }
      double cellVal = (xip ? cellSED->xi()(k) : cellSED->xip()(k));
      double ringVal = (xip ? ringSED->xi()(k) : ringSED->xip()(k));
      double cellErr = TMath::Sqrt(xip ? cellSED->xiErr()(k,k) : cellSED->xipErr()(k,k));
      h[0][k]->Fill(cellVal/cellErr);
      h[1][k]->Fill((cellVal - ringVal)/cellErr);
      h[2][k]->Fill(cellVal);
      h[3][k]->Fill(ringVal);
      h[4][k]->Fill(cellErr);
      hCor->Fill(cellVal, ringVal);
    }
  }
  for (unsigned int i = 1; i <= nBins; i++) 
    for (unsigned short l = 0; l < 5; l++)
      for (unsigned short k = 0; k < 5; k++)
        hAll[l]->SetBinContent(k+1, i, h[l][k]->GetBinContent(i));
  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
  for (unsigned short l = 0; l < 5; l++) {
    for (unsigned short k = 0; k < 5; k++)
      h[l][k]->Write();
    hAll[l]->Write();
  }
  ff->Close();
  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 500);
  c1->Divide(3,1);
  c1->cd(1); hCor->Draw();
  TLatex l; l.SetNDC();
  TString corrStr = Form("Correlation : %.1f%%", 100*hCor->GetCorrelationFactor());
  l.DrawLatex(0.2, 0.8, corrStr);
  c1->cd(2); hAll[0]->Draw("COL");
  c1->cd(3); hAll[1]->Draw("COL");
  c1->Print(fileName + ".png");
  c1->Print(fileName + ".eps");
  return true;
}

shapeErrorDist()

TH1D * DigitMonitor::shapeErrorDist	(	unsigned int	k,
		const TString &	name,
		int	nBins,
		double	rMin,
		double	rMax,
		double	mean = 0 ) const

Definition at line 426 of file DigitMonitor.cxx.

{
  TH1D* h = new TH1D(name, "Shape error", nBins, rMin, rMax);
  h->GetXaxis()->SetTitle("#epsilon^{2}");
  
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) { 
      std::unique_ptr<OFC> ofc = history->ofc(j);
      double delta = history->data(j)->delta(k) - mean;
      double adcMax = history->data(j)->adcMax();
      h->Fill((delta*delta - ofc->Gamma()(k, k))/(adcMax*adcMax));
    }
  }
  cout << sqrt(h->GetMean())*100 << " +/- " << sqrt(h->GetMeanError())*100 << endl;
  return h;
}

statParams()

bool MonitorBase::statParams ( const std::vector< DataFuncSet > & funcs, const std::vector< DataFuncArgs > & args, TVectorD & mean, TVectorD & meanErr, TMatrixD & covMatrix, TMatrixD & covMatrixErr, const FilterParams & f = FilterParams() ) const

inherited

Definition at line 375 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  Averager avg(funcs.size());
    
  for (unsigned int i = 0; i < m_interface->nChannels(); i++) {
    if ((i+1) % printPeriodicity == 0) cout << "Cell # " << i+1 << "/" << m_interface->nChannels() << endl;
    const History* history = pass(i, f);
    if (!history) continue;
    for (unsigned int j = 0; j < history->nData(); j++) {
      if (!f.passEvent(*history->data(j))) continue;
      TVectorD vals(funcs.size());
      for (unsigned int k = 0; k < funcs.size(); k++)
        vals(k) = funcs[k].val(*history->data(j), args[k]);
      avg.fill(vals);
    }
  }
  
  mean = avg.means();
  meanErr = avg.meanErrors();
  covMatrix = avg.covarianceMatrix();
  covMatrixErr = avg.covarianceMatrixErrors();
 
  return true;
}

str()

TString MonitorBase::str ( CombinationType comb )

staticinherited

Definition at line 56 of file LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx.

{
  switch (comb) {
    case AverageValue : return "average";
    case MaxValue     : return "max value";
    case MinValue     : return "min value";
    case TotalValue   : return "total value";
    case RMSValue     : return "RMS";
  }
  return "";
}

timeDist()

TH1D * DigitMonitor::timeDist	(	const TString &	name,
		int	nBins,
		double	tMin = -25,
		double	tMax = 25 ) const

Definition at line 44 of file DigitMonitor.cxx.

{ return dist(&Data::_ofcTime, DataFuncArgs(), name, nBins, tMin, tMax, "Time distribution", "Time (ns)"); }

Member Data Documentation

m_cellCache

std::unique_ptr<const History> LArSamples::AbsLArCells::m_cellCache

mutableprivateinherited

Definition at line 55 of file AbsLArCells.h.

m_cellInfoCache

std::vector<std::unique_ptr<CellInfo> > LArSamples::AbsLArCells::m_cellInfoCache

mutableprivateinherited

Definition at line 56 of file AbsLArCells.h.

m_interface

const Interface* LArSamples::MonitorBase::m_interface

privateinherited

Definition at line 87 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h.

m_pos

unsigned int LArSamples::AbsLArCells::m_pos

mutableprivateinherited

Definition at line 54 of file AbsLArCells.h.

printPeriodicity

const int LArSamples::MonitorBase::printPeriodicity = 10000

staticinherited

Definition at line 79 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h.

---

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

• DigitMonitor.h
• DigitMonitor.cxx