LArSamples::OccupancyMonitor Class Reference

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

#include <OccupancyMonitor.h>

Inheritance diagram for LArSamples::OccupancyMonitor:

Collaboration diagram for LArSamples::OccupancyMonitor:

Public Member Functions
	OccupancyMonitor (const Interface &interface)
	Constructor.
TH1I *	runOccupancyHistory () const
TH1I *	eventOccupancyHistory (unsigned int occForDump=0) const
TH1I *	eventOccupancy () const
TH1I *	cellOccupancy (int minForPrintout=-1) const
TH1I *	febOccupancy (CaloId calo) const
TH2D *	partitionOccupancy (PartitionId part) const
TH2D *	etaPhiOccupancy (CaloId calo, short layer, bool useEnergy=false) const
void	cellAndRingOccupancy (CaloId calo, unsigned int nMin=20) const
void	dump (short verbosity=2) 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
unsigned int	nChannels () const
const History *	getCellHistory (unsigned int i) const
const History *	getSCHistory (unsigned int i) const
const CellInfo *	getCellInfo (unsigned int i) const
virtual const History *	newCellHistory (unsigned int i) const
virtual const History *	cellHistory (unsigned int i) const
virtual 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
void	resetCellInfoCache ()

Static Public Member Functions
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
const History *	m_cellCache
std::vector< CellInfo * >	m_cellInfoCache

Detailed Description

storage of the time histories of all the cells

Definition at line 20 of file OccupancyMonitor.h.

Constructor & Destructor Documentation

OccupancyMonitor()

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

inline

Constructor.

Definition at line 26 of file OccupancyMonitor.h.

: MonitorBase(interface) { }

Member Function Documentation

cachePos()

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

inlineinherited

Definition at line 48 of file AbsLArCells.h.

{ return m_pos; }

cellAndRingOccupancy()

void OccupancyMonitor::cellAndRingOccupancy	(	CaloId	calo,
		unsigned int	nMin = 20 ) const

Definition at line 258 of file OccupancyMonitor.cxx.

{
  const unsigned int nPhiRings = 5516; // == Geo::nPhiRings()
  using block_t = unsigned int[4][3];
  block_t nCells, nRings;
  std::vector<block_t> ringOccupancy (nPhiRings);
  
  for (unsigned int i = 0; i < 4; i++)
    for (unsigned int j = 0; j < 3; j++) {
      for (unsigned int k = 0; k < nPhiRings; k++) ringOccupancy[k][i][j] = 0;
      nCells[i][j] = 0;
      nRings[i][j] = 0;
    }
    
  for (unsigned int i = 0; i < nChannels(); i++) {
    if (i%10000 == 0) cout << "Processing hash = " << i << endl;
    const History* history = cellHistory(i);
    if (!history) continue;
    if (!Id::matchCalo(history->cellInfo()->calo(), calo)) continue;
    unsigned int n[3] = {0,0,0};
    int layer = history->cellInfo()->layer();
    for (unsigned int k = 0; k < history->nData(); k++) {
      const auto globalPhiRing = history->cellInfo()->globalPhiRing();
      if (globalPhiRing<0) continue;
      ringOccupancy[globalPhiRing][layer][(int)history->data(k)->gain()]++;
      n[(int)history->data(k)->gain()]++;
    }
    for (unsigned int g = 0; g < 3; g++) if (n[g] >= nMin) nCells[layer][g]++;
    //cout << n[1] << " " << nCells[layer][1] << " " << layer << endl;
  }
  for (unsigned int i = 0; i < 4; i++)
    for (unsigned int j = 0; j < 3; j++) {
      for (unsigned int k = 0; k < nPhiRings; k++)
        if (ringOccupancy[k][i][j] >= nMin) nRings[i][j]++;
  }
  
  for (unsigned int i = 0; i < 4; i++) {
    cout << "layer " << i << endl;
    const auto nChann = Geo::nChannels(calo,i);
    const auto nEta = Geo::nEta(calo,i);
    for (unsigned int j = 0; j < 3; j++) {
      float nCellFract= infinity;
      if (nChann>0) nCellFract = nCells[i][j]*100.0/nChann;
      float nRingFract = infinity;
      if (nEta > 0) nRingFract = nRings[i][j]*100.0/nEta;
      cout << " gain " << j << " : " 
           << nCells[i][j] << " cells of " << nChann << " (" << nCellFract << " %), " 
           << nRings[i][j] << " rings of " << nEta << " (" << nRingFract << " %)." << endl;
    }
  }
  for (unsigned int i = 0; i < 4; i++) {
    const auto nChann = Geo::nChannels(calo,i);
    const auto nEta = Geo::nEta(calo,i);
    cout << i << " & " << nChann << " & " << nEta << " & ";
    for (unsigned int j = 0; j < 3; j++){
       float val = infinity;
       if (nChann>0) val = nCells[i][j]*100.0/nChann;
       cout << Form("%4.1f\\%% & ", val);
    }
    for (unsigned int j = 0; j < 3; j++){
       float val = infinity;
       if (nEta > 0) val = nRings[i][j]*100.0/nEta;
       cout << Form("%4.1f\\%% & ", val);
    }
    cout << endl;
  }
}

cellCache()

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

inlineinherited

Definition at line 47 of file AbsLArCells.h.

{ return m_cellCache; }

cellHistory()

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

virtualinherited

Reimplemented in LArSamples::Interface.

Definition at line 59 of file AbsLArCells.cxx.

{ 
  if (m_pos == i) return m_cellCache;
  resetCache();
  const History* history = newCellHistory(i);
  if (!history) return nullptr;
  m_cellCache = history;
  m_pos = i;
  return m_cellCache;
}

cellInfo()

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

virtualinherited

Definition at line 71 of file AbsLArCells.cxx.

{
  const CellInfo* info = cellInfoCache(i);
  if (info) return (info->isValid() ? new CellInfo(*info) : nullptr);
  info = getCellInfo(i);
  if (info)  m_cellInfoCache[i] =  new CellInfo(*info, false); 
  //m_cellInfoCache[i] = (info ? new CellInfo(*info, false) : new CellInfo());
  return info;
}

cellInfoCache()

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

inherited

Definition at line 82 of file AbsLArCells.cxx.

{
  return m_cellInfoCache[i];
}

cellOccupancy()

TH1I * OccupancyMonitor::cellOccupancy

(

int

minForPrintout = -1

)

const

Definition at line 120 of file OccupancyMonitor.cxx.

{
  std::vector<unsigned int> occ(nChannels());
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);  
    if (history) {
      occ[i] = occ[i] + history->nData();
    }
  }
  
  unsigned int oMin = 999999999, oMax = 0;
  
  for (std::vector<unsigned int>::const_iterator cell = occ.begin();
       cell != occ.end(); ++cell) 
  {
    if (*cell < oMin) oMin = *cell;
    if (*cell > oMax) oMax = *cell;
  }
  if (oMin > oMax) {
    oMin = 0;
    oMax = 1;
  }
  
  //std::map<std::pair<unsigned int, std::string>, unsigned int> names;
  std::map<std::pair<unsigned int, unsigned int>, std::string> names;
  TH1I* h = new TH1I("occ", "Cell occupancies", oMax - oMin + 1, oMin - 0.5, oMax + 0.5);
  h->GetXaxis()->SetTitle("Occupancy");
  h->GetYaxis()->SetTitle("Number of cells");
 
  unsigned int i = 0;
  for (std::vector<unsigned int>::const_iterator cellOcc = occ.begin();
       cellOcc != occ.end(); ++cellOcc, ++i) {
    h->Fill(*cellOcc);
    const History* history = cellHistory(i);
    if (!history) continue;
    if (minForPrintout >= 0) {
      TString heading = history->cellInfo()->location(2);
      if (history->nData() > 0) heading += " " + history->data(0)->problems();
      names[ std::pair<unsigned int, unsigned int>(*cellOcc, i) ] = heading.Data();
    }
  }
  
  if (minForPrintout >= 0) {
/*    for (std::map<std::pair<unsigned int, std::string>, unsigned int>::reverse_iterator cell = names.rbegin();
    cell != names.rend(); ++cell)
    {
    if (cell->first.first < (unsigned int)minForPrintout) break;
    cout << cell->first.first << " : " << cell->second << ", " << cell->first.second << endl;
  }*/
    for (std::map<std::pair<unsigned int, unsigned int>, std::string>::reverse_iterator cell = names.rbegin();
         cell != names.rend(); ++cell)
    {
      if (cell->first.first < (unsigned int)minForPrintout) break;
      cout << cell->first.first << " : " << cell->first.second << ", " << cell->second << endl;
    }
  }
 
  return h;
}

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/3]

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/3]

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;
}

dump() [3/3]

void OccupancyMonitor::dump

(

short

verbosity = 2

)

const

Definition at line 241 of file OccupancyMonitor.cxx.

{
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    cout << history->cellInfo()->location(verbosity);
    if (verbosity < 5) 
      cout << " : " << history->nData() << " pulse" << endl;
    else {
      cout << " : ";
      for (unsigned int i = 0; i < history->nData(); i++)
    cout << "E = " << history->data(i)->energy() << " MeV, " << endl;
    }
  }
}

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;
}

etaPhiOccupancy()

TH2D * OccupancyMonitor::etaPhiOccupancy	(	CaloId	calo,
		short	layer,
		bool	useEnergy = false ) const

Definition at line 231 of file OccupancyMonitor.cxx.

{
  DataFuncSet func = (useEnergy ? MonitorBase::func("energy") : MonitorBase::func("count"));
  return etaPhiMap(func, DataFuncArgs(), Form("occupancy_%s_%d", Id::str(calo).Data(), layer), 
                   calo, layer, 
                   Form("Occupancies for partition %s, layer %d", Id::str(calo).Data(), layer),
                   TotalValue);                   
}

eventOccupancy()

TH1I * OccupancyMonitor::eventOccupancy

(

)

const

Definition at line 86 of file OccupancyMonitor.cxx.

{
  std::map<unsigned int, unsigned int> occ;
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int k = 0; k < history->nData(); k++)
      occ[history->data(k)->event()]++;
  }
  
  unsigned int oMin = 999999999, oMax = 0;
  
  for (std::map<unsigned int, unsigned int>::const_iterator event = occ.begin();
       event != occ.end(); ++event) 
  {
    if (event->second < oMin) oMin = event->second;
    if (event->second > oMax) oMax = event->second;
  }
  if (oMin > oMax) {
    oMin = 0;
    oMax = 1;
  }
  
  TH1I* h = new TH1I("occ", "Event occupancies", oMax - oMin + 1, oMin - 0.5, oMax + 0.5);
  h->GetXaxis()->SetTitle("Occupancy");
  h->GetYaxis()->SetTitle("Number of events");
  for (std::map<unsigned int, unsigned int>::const_iterator event = occ.begin();
       event != occ.end(); ++event) 
    h->Fill(event->second);
  
  return h;
}

eventOccupancyHistory()

TH1I * OccupancyMonitor::eventOccupancyHistory

(

unsigned int

occForDump = 0

)

const

Definition at line 60 of file OccupancyMonitor.cxx.

{
  std::map<std::pair<unsigned int, unsigned int>, unsigned int> occ;
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int k = 0; k < history->nData(); k++)
      occ[std::make_pair(history->data(k)->run(),history->data(k)->event())]++;
  }
  
  TH1I* h = new TH1I("occ", "Event occupancies", occ.size(), 0, occ.size());
  h->GetXaxis()->SetTitle("Recorded Event");
  h->GetYaxis()->SetTitle("Occupancy");
  
  unsigned int bin = 1;
  for (std::map<std::pair<unsigned int, unsigned int>, unsigned int>::const_iterator event = occ.begin();
       event != occ.end(); ++event, ++bin) 
  {
    if (occForDump > 0 && occForDump <= event->second) 
      cout << "run = " << event->first.first << " event = " << event->first.second << " occ = " << event->second << endl;
    h->SetBinContent(bin, event->second);
  }
  return h;
}

febOccupancy()

TH1I * OccupancyMonitor::febOccupancy

(

CaloId

calo

)

const

Definition at line 181 of file OccupancyMonitor.cxx.

{
  std::map<unsigned int, unsigned int> occ;
  FilterParams f; 
  f.addCalo(calo);
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = pass(i, f);
    if (!history) continue;
    for (unsigned int k = 0; k < history->nData(); k++)
      occ[history->cellInfo()->feb()]++;
  }
  
  unsigned int oMin = 999999999, oMax = 0;
  
  for (std::map<unsigned int, unsigned int>::const_iterator feb = occ.begin();
       feb != occ.end(); ++feb) 
  {
    if (feb->second < oMin) oMin = feb->second;
    if (feb->second > oMax) oMax = feb->second;
  }
  if (oMin > oMax) {
    oMin = 0;
    oMax = 1;
  }
  
  TH1I* h = new TH1I("occ", "FEB occupancies", oMax - oMin + 1, oMin - 0.5, oMax + 0.5);
  h->GetXaxis()->SetTitle("Occupancy");
  h->GetYaxis()->SetTitle("Number of FEBs");
  for (std::map<unsigned int, unsigned int>::const_iterator feb = occ.begin();
       feb != occ.end(); ++feb) 
    h->Fill(feb->second);
  
  return h;
}

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();
}

getCellHistory()

const History * MonitorBase::getCellHistory ( unsigned int i ) const

virtualinherited

Implements LArSamples::AbsLArCells.

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

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

getCellInfo()

const CellInfo * MonitorBase::getCellInfo ( unsigned int i ) const

virtualinherited

Reimplemented from LArSamples::AbsLArCells.

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

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

getSCHistory()

const History * MonitorBase::getSCHistory ( unsigned int i ) const

virtualinherited

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; }

nChannels()

unsigned int MonitorBase::nChannels ( ) const

virtualinherited

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 35 of file AbsLArCells.h.

{ return Definitions::nChannelsSC; }

newCellHistory()

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

virtualinherited

Definition at line 44 of file AbsLArCells.cxx.

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

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;
}

partitionOccupancy()

TH2D * OccupancyMonitor::partitionOccupancy

(

PartitionId

part

)

const

Definition at line 217 of file OccupancyMonitor.cxx.

{
  TH2D* occ = Geo::partitionHist(part, TString("occupancy_") + Id::str(part), "Occupancies for partition " +  Id::str(part));
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    if (Id::partition(history->cellInfo()->calo()) != part) continue;
    occ->Fill(history->cellInfo()->feb(), history->cellInfo()->channel(), history->nData());
  }
 
  return occ;
}

pass()

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

inherited

Definition at line 99 of file AbsLArCells.cxx.

{ 
  //std::cout << "Called AbsLArCells with hash " << i  << std::endl;
  if (!f.passHash(i)) return nullptr;
  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);
  delete 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 34 of file AbsLArCells.cxx.

{
  if (m_cellCache) {
    delete m_cellCache;
    m_cellCache = nullptr;
  }
  m_pos = nChannels() + 1;
}

resetCellInfoCache()

void AbsLArCells::resetCellInfoCache ( )

inherited

Definition at line 114 of file AbsLArCells.cxx.

{
  unsigned int i = 0;
  for (std::vector<CellInfo*>::iterator cellInfo = m_cellInfoCache.begin();
       cellInfo != m_cellInfoCache.end(); ++cellInfo, i++)
    if (*cellInfo) {
      delete *cellInfo;
      *cellInfo = 0;
    }
}

runOccupancyHistory()

TH1I * OccupancyMonitor::runOccupancyHistory

(

)

const

Definition at line 28 of file OccupancyMonitor.cxx.

{
  std::map<unsigned int, unsigned int> occ;
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    for (unsigned int k = 0; k < history->nData(); k++)
      occ[history->data(k)->run()]++;
  }
  
  unsigned int runMin = 0;
  unsigned int runMax = 1;
  if (!occ.empty()) {
    runMin = occ.begin()->first;
    runMax = occ.rbegin()->first;
  }
 
  TH1I* h = new TH1I("occ", "Run occupancies", runMax - runMin + 1, runMin - 0.5, runMax + 0.5);
  h->GetXaxis()->SetTitle("Run number");
  h->GetYaxis()->SetTitle("Occupancy");
  for (std::map<unsigned int, unsigned int>::const_iterator run = occ.begin();
       run != occ.end(); ++run) 
  {
    //cout << event->first << " " << event->second << endl;
    int bin = run->first - runMin + 1;
    h->SetBinContent(bin, run->second);
    h->GetXaxis()->SetBinLabel(bin, Form("%d", run->first));
  }
  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 "";
}

Member Data Documentation

m_cellCache

const History* LArSamples::AbsLArCells::m_cellCache

mutableprivateinherited

Definition at line 55 of file AbsLArCells.h.

m_cellInfoCache

std::vector<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:

• OccupancyMonitor.h
• OccupancyMonitor.cxx