LArSamples::Interface Class Reference

#include <Interface.h>

Inheritance diagram for LArSamples::Interface:

Collaboration diagram for LArSamples::Interface:

Public Member Functions
	Interface (const Accessor &accessor)
	Constructor.
virtual	~Interface ()
unsigned int	nEvents () const
unsigned int	nRuns () const
const EventData *	eventData (unsigned int i) const
const RunData *	runData (unsigned int i) const
unsigned int	historySize (unsigned int i) const
unsigned int	historySizeSC (unsigned int i) const
HistoryIterator	begin (unsigned int pos=0, double eMin=-1, double adcMaxMin=-1) const
unsigned int	end () const
unsigned int	size () const
unsigned int	nFilledChannels () const
void	printFilledRanges (unsigned int skip=0) const
bool	highEData (double eCut, TArrayI &hashes, TArrayI &indices) const
HistoryIterator	findEtaPhi (CaloId calo, short layer, short iEta, short iPhi, short region=0) const
HistoryIterator	findFebChannel (CaloId calo, short feb, short channel) const
HistoryIterator	findFTSlotChannel (CaloId calo, short ft, short slot, short channel) const
bool	neighbors (const CellInfo &cell, double dRCut, std::vector< unsigned int > &hashes) const
bool	firstNeighbors (unsigned int hash, std::vector< unsigned int > &hashes, short layer=-2) const
bool	data (const std::vector< unsigned int > &hashes, const EventData &event, std::vector< const Data * > &data) const
Interface *	merge (const Interface &other, const TString &fileName) const
Interface *	merge (const Interface &other, const TString &fileName, const TString &LBFile) const
Interface *	filter (const TString &sel, const TString &fileName, const TString &tweaks="") const
Interface *	filter (const FilterParams &filterParams, const DataTweaker &tweaker, const TString &fileName) const
Interface *	makeTemplate (const TString &fileName) const
Interface *	refit (const TString &newFileName, Chi2Params pars=DefaultChi2) const
const Accessor &	accessor () const
bool	isValid () const
TH1D *	Draw (const TString &var, int nBins, double xMin, double xMax, const TString &sel="", const TString &opt="") const
TH2D *	Draw (const TString &varList, int nBinsX, double xMin, double xMax, int nBinsY, double yMin, double yMax, const TString &sel="", const TString &opt="") const
TH2D *	DrawPartition (PartitionId partition, const TString &var, const TString &sel="", const TString &opt="", CombinationType comb=TotalValue) const
TH2D *	DrawEtaPhi (CaloId calo, short layer, const TString &var, const TString &sel="", const TString &opt="", CombinationType comb=TotalValue) const
bool	Scan (const TString &vars, const TString &sel="", unsigned int verbosity=1) const
bool	Scan (const TString &vars, CombinationType comb, const TString &sel="", const TString &ranges="", unsigned int verbosity=1) const
bool	Show (unsigned int hash, unsigned int verbosity=1) const
bool	Show (const TString &sel="", unsigned int verbosity=1) const
bool	ShowEvents (const TString &sel="", unsigned int verbosity=1) const
bool	ShowRuns (unsigned int verbosity=1) const
bool	ShowStats (const TString &varList, const TString &sel="", bool withErrors=false) const
OccupancyMonitor	occupancyMonitor () const
DigitMonitor	digitMonitor () const
FitMonitor	fitMonitor () const
bool	dumpEventTuple (const TString &variables, const TString &fileName) const
void	setShapeErrorGetter (const AbsShapeErrorGetter *err)
void	setShapeError (double k)
void	setShapeError (const TString &fileName)
const AbsShapeErrorGetter *	shapeErrorGetter () const
bool	writeToFile (const TString &fileName) const
const History *	cellHistory (unsigned int i) 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 CellInfo *	cellInfo (unsigned int i) const
virtual unsigned int	nChannels () 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 Interface *	open (const TString &fileName)
static Interface *	open (const std::vector< TString > &fileNames)
static Interface *	openList (const TString &fileList)
static Interface *	openWild (const TString &wcName)
static Interface *	merge (const std::vector< const Interface * > &interfaces, const TString &fileName)
static Interface *	merge (const std::vector< const Interface * > &interfaces, const TString &fileName, const TString &LBFile)
static Interface *	merge (const TString &listFileName, const TString &fileName)
static Interface *	merge (const TString &listFileName, const TString &fileName, const TString &LBFile)
static bool	filterAndMerge (const TString &listFileName, const TString &outFile, const TString &filters, const TString &tweaks="")
static TString	addSuffix (const TString &fileName, const TString &suffix)

Private Attributes
const Accessor *	m_accessor
const AbsShapeErrorGetter *	m_shapeErrorGetter
bool	m_ownShapeErrorGetter
std::vector< std::vector< unsigned int > * >	m_neighborCache
std::vector< unsigned int >	m_neighborHistoryPos
std::vector< const History * >	m_neighborHistories
unsigned int	m_pos
const History *	m_cellCache
std::vector< CellInfo * >	m_cellInfoCache

Detailed Description

Definition at line 36 of file Interface.h.

Constructor & Destructor Documentation

Interface()

LArSamples::Interface::Interface ( const Accessor & accessor )

inline

Constructor.

Definition at line 41 of file Interface.h.

        : m_accessor(&accessor), m_shapeErrorGetter(0), m_ownShapeErrorGetter(0), m_neighborCache(nChannels()) { }

~Interface()

Interface::~Interface ( )

virtual

Definition at line 63 of file Interface.cxx.

{ 
  if (m_ownShapeErrorGetter) delete m_shapeErrorGetter;
  delete m_accessor; 
 
  for (std::vector<unsigned int>* neighbors : m_neighborCache)
    delete neighbors;
}

Member Function Documentation

accessor()

const Accessor & LArSamples::Interface::accessor ( ) const

inline

Definition at line 94 of file Interface.h.

{ return *m_accessor; }

addSuffix()

TString Interface::addSuffix ( const TString & fileName, const TString & suffix )

static

Definition at line 382 of file Interface.cxx.

{
  TString rootName = fileName;
  if (fileName(fileName.Length() - 5, 5) == ".root") rootName = fileName(0, fileName.Length() - 5);
  
  return  rootName + "_" + suffix + ".root";
}

begin()

HistoryIterator Interface::begin	(	unsigned int	pos = 0,
		double	eMin = -1,
		double	adcMaxMin = -1 ) const

Definition at line 182 of file Interface.cxx.

{
  return HistoryIterator(*this, pos, eMin, adcMaxMin);
}

cachePos()

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

inlineinherited

Definition at line 48 of file AbsLArCells.h.

{ return m_pos; }

cellCache()

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

inlineinherited

Definition at line 47 of file AbsLArCells.h.

{ return m_cellCache; }

cellHistory()

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

virtual

Reimplemented from LArSamples::AbsLArCells.

Definition at line 125 of file Interface.cxx.

{ 
  const History* history = accessor().cellHistory(i);
  if (history) {
    history->setShapeErrorGetter(m_shapeErrorGetter);
    history->setInterface(this);
  }
  return history;
}

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

data()

bool Interface::data	(	const std::vector< unsigned int > &	hashes,
		const EventData &	event,
		std::vector< const Data * > &	data ) const

Definition at line 770 of file Interface.cxx.

{
  if (hashes != m_neighborHistoryPos) {
    for (const History* history : m_neighborHistories)
      delete history;
    m_neighborHistories.clear();
    m_neighborHistoryPos.clear();
    for (std::vector<unsigned int>::const_iterator hash = hashes.begin(); hash != hashes.end(); ++hash) {
      const History* history = AbsLArCells::newCellHistory(*hash);// bypasses history caching in order not to invalidate cell
      m_neighborHistories.push_back(history);
      m_neighborHistoryPos.push_back(*hash);
    }
  }
 
  for (const History* history : m_neighborHistories) {
    if (!history) continue;
    const Data* dataForEvent = history->data_for_event(event);
    if (dataForEvent) data.push_back(new Data(*dataForEvent));
  }
  return true;
}

digitMonitor()

DigitMonitor LArSamples::Interface::digitMonitor ( ) const

inline

Definition at line 115 of file Interface.h.

{ return DigitMonitor(*this); }

Draw() [1/2]

TH1D * Interface::Draw	(	const TString &	var,
		int	nBins,
		double	xMin,
		double	xMax,
		const TString &	sel = "",
		const TString &	opt = "" ) const

Definition at line 463 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return nullptr;
  
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(var, vars, funcs, args) || funcs.size() != 1) {
    cout << "Invalid variable specification " << var << endl;
    return nullptr;
  }
 
  TString title = vars[0];
  if (TString(sel) != "") title = title + ", " + sel;
  
  TH1D* h = m.dist(funcs[0], args[0], vars[0], nBins, xMin, xMax, 
           title, vars[0], "digits", f);
  if (!h) return nullptr;
  h->Draw(opt);
  return h;
}

Draw() [2/2]

TH2D * Interface::Draw	(	const TString &	varList,
		int	nBinsX,
		double	xMin,
		double	xMax,
		int	nBinsY,
		double	yMin,
		double	yMax,
		const TString &	sel = "",
		const TString &	opt = "" ) const

Definition at line 488 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return nullptr;
  
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(varList, vars, funcs, args) || funcs.size() != 2) {
    cout << "Invalid variable specification " << varList << endl;
    return nullptr;
  }
 
  TString title = vars[1] + " vs. " + vars[0];
 
  if (TString(sel) != "") title = title + ", " + sel;
  TH2D* h = m.dist(funcs[0], args[0], funcs[1], args[1], vars[0] + "_" + vars[1], 
           nBinsX, xMin, xMax,  nBinsY, yMin, yMax, 
           title, vars[0], vars[1], f);
  if (!h) return nullptr;
  h->Draw(opt);
  return h;
}

DrawEtaPhi()

TH2D * Interface::DrawEtaPhi	(	CaloId	calo,
		short	layer,
		const TString &	var,
		const TString &	sel = "",
		const TString &	opt = "",
		CombinationType	comb = TotalValue ) const

Definition at line 543 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return nullptr;
 
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(var, vars, funcs, args) || funcs.size() != 1) {
    cout << "Invalid variable specification " << var << endl;
    return nullptr;
  }
  TString title = var;
  title += Form(", %s, layer %d", Id::str(calo).Data(), layer);
  if (TString(sel) != "") title = title + ", " + sel;
 
  TH2D* h = m.etaPhiMap(funcs[0], args[0], var, calo, layer, title, comb, f);
  if (!h) return nullptr;
  h->Draw(opt);
  return h;
}

DrawPartition()

TH2D * Interface::DrawPartition	(	PartitionId	partition,
		const TString &	var,
		const TString &	sel = "",
		const TString &	opt = "",
		CombinationType	comb = TotalValue ) const

Definition at line 516 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return nullptr;
 
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(var, vars, funcs, args) || funcs.size() != 1) {
    cout << "Invalid variable specification " << var << endl;
    return nullptr;
  }
 
  TString title = var;
  title = title + ", " + Id::str(partition);
  if (TString(sel) != "") title = title + ", " + sel;
 
  TH2D* h = m.partitionMap(funcs[0], args[0], var, partition, title, comb, f);
  if (!h) return nullptr;
  h->Draw(opt);
  return h;
}

dumpEventTuple()

bool Interface::dumpEventTuple	(	const TString &	variables,
		const TString &	fileName ) const

Definition at line 793 of file Interface.cxx.

{
  std::vector<float*> floatVars;
  std::vector<int*> intVars;
  std::vector<std::vector<float>*> floatVects;
  std::vector<std::vector<int>*> intVects;
  std::map<TString, unsigned int> varIndex;
 
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(variables, vars, funcs, args)) return false;
 
  cout << "Making trees..." << endl;
 
  TFile* flatFile = TFile::Open(fileName + "_tmpFlatFile.root", "RECREATE");
  TTree* flatTree = new TTree("flatTree", "Flat tree");
  
  TFile* eventFile = TFile::Open(fileName, "RECREATE");
  TTree* eventTree = new TTree("eventTree", "Event tree");
  
  for (unsigned int j = 0; j < vars.size(); j++) {
    unsigned int index = 0;
    if (funcs[j].isNull()) return false;
    if (funcs[j].isInt()) {
      int* varCont = new int(0);
      std::vector<int>* vectCont = new std::vector<int>();      
      index = intVars.size();
      intVars.push_back(varCont);
      intVects.push_back(vectCont);
      flatTree->Branch(vars[j], varCont);
      eventTree->Branch(vars[j], vectCont);
    }
    else {
      float* varCont = new float(0);
      std::vector<float>* vectCont = new std::vector<float>();      
      index = floatVars.size();
      floatVars.push_back(varCont);
      floatVects.push_back(vectCont);
      flatTree->Branch(vars[j], varCont);
      eventTree->Branch(vars[j], vectCont);
    }
    varIndex[vars[j]] = index;
    cout << vars[j] << " -> " << index << endl;
  }
  std::map< unsigned int, std::map< unsigned int, std::vector<long long> > > runEventIndices;
  cout << "Making flat ntuple" << endl;
  unsigned int count = 0;
  for (HistoryIterator iter = begin(); iter.isValid(); iter.next()) {
    const History* hist = iter.history();
    count++;
    if (count % 100 == 0) cout << "Processing entry " << count << " (hash = " << iter.pos() << "), size = " << hist->nData() << endl;
    for (unsigned int k = 0; k < hist->nData(); k++) {
      for (unsigned int j = 0; j < vars.size(); j++) {
        if (funcs[j].isInt())
          *intVars[varIndex[vars[j]]] = int(funcs[j].intVal(*hist->data(k), args[j]));
        else
          *floatVars[varIndex[vars[j]]] = funcs[j].doubleVal(*hist->data(k), args[j]);
      }
      runEventIndices[hist->data(k)->run()][hist->data(k)->event()].push_back(flatTree->GetEntries());
      flatTree->Fill();
    }
  }
  
  cout << "Making event tuple" << endl;
  unsigned int runCount = 0;
  for (const auto& run : runEventIndices) {
    runCount++;
    cout << "Processing run " << run.first << " (" << runCount << " of " << runEventIndices.size() << ")" << endl;
    unsigned int eventCount = 0;
    for (const auto& event : run.second) {
      eventCount++;
      if (eventCount % 1000 == 0) 
        cout << "  processing event " << event.first << " (" << eventCount << " of " << run.second.size() << "), size = " << event.second.size() << endl;
      for (unsigned int j = 0; j < intVects.size(); j++) intVects[j]->clear();
      for (unsigned int j = 0; j < floatVects.size(); j++) floatVects[j]->clear();
      for (long long index : event.second) {
        flatTree->GetEntry(index);
       for (unsigned int j = 0; j < vars.size(); j++) {
         if (funcs[j].isInt())
           intVects[varIndex[vars[j]]]->push_back(*intVars[varIndex[vars[j]]]);
         else
           floatVects[varIndex[vars[j]]]->push_back(*floatVars[varIndex[vars[j]]]);
        }
      }
      eventTree->Fill();
    }
  }
  
  cout << "Writing data..." << endl;  
  flatFile->cd();
  flatTree->Write();
  eventFile->cd();
  eventTree->Write();  
  
  cout << "Cleaning up..." << endl;
  delete eventTree;
  delete eventFile;
  
  delete flatTree;
  delete flatFile;
  for (unsigned int j = 0; j < intVects.size(); j++) delete intVects[j];
  for (unsigned int j = 0; j < floatVects.size(); j++) delete floatVects[j];      
  for (unsigned int j = 0; j < intVars.size(); j++) delete intVars[j];
  for (unsigned int j = 0; j < floatVars.size(); j++) delete floatVars[j];      
  
  cout << "Done!" << endl;
  return true;
}

end()

unsigned int LArSamples::Interface::end ( ) const

inline

Definition at line 61 of file Interface.h.

{ return nChannels(); }

eventData()

const EventData * LArSamples::Interface::eventData ( unsigned int i ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 54 of file Interface.h.

{ return accessor().eventData(i); }

filter() [1/2]

Interface * Interface::filter	(	const FilterParams &	filterParams,
		const DataTweaker &	tweaker,
		const TString &	fileName ) const

Definition at line 391 of file Interface.cxx.

{
  TreeAccessor* newAcc = TreeAccessor::filter(accessor(), filterParams, fileName, tweaker);
  return new Interface(*newAcc);
}

filter() [2/2]

Interface * Interface::filter	(	const TString &	sel,
		const TString &	fileName,
		const TString &	tweaks = "" ) const

Definition at line 363 of file Interface.cxx.

{
  FilterParams f;
  if (!f.set(sel)) return nullptr;
 
  DataTweaker tweak;
  if (!tweak.set(tweaks)) return nullptr;
 
  TString thisFN = fileName;
  if (thisFN.Index(".root") < 0 && dynamic_cast<const TreeAccessor*>(&accessor())) {
    auto pAccess = dynamic_cast<const TreeAccessor*>(&accessor());
    if (not pAccess) return nullptr;
    TString newFN = addSuffix(pAccess->fileName(), fileName);
    if (newFN != "") thisFN = std::move(newFN);
  }
  return filter(f, tweak, thisFN);
}

filterAndMerge()

bool Interface::filterAndMerge ( const TString & listFileName, const TString & outFile, const TString & filters, const TString & tweaks = "" )

static

Definition at line 293 of file Interface.cxx.

{
  FilterList filterList;
  
  TObjArray* list = filters.Tokenize(",;");
  if (list->GetEntries() == 0) {
    cout << "No filtering specified, exiting.";
    delete list;
    return 0;
  }
  
  for (int k = 0; k < list->GetEntries(); k++) {
    TObjString* tobs = (TObjString*)(list->At(k));
    TObjArray* items = tobs->String().Tokenize(":");
    if (items->GetEntries() != 2) {
      cout << "Invalid filter entry " << tobs->String() << ", exiting." << endl;
      delete list;
      delete items;
      return 0;
    }
    TString params = ((TObjString*)(items->At(0)))->String();
    TString suffix = ((TObjString*)(items->At(1)))->String();    
    FilterParams f;
    if (!f.set(params)) return 0;
    cout << "---" << endl;
    filterList.add(f, addSuffix(outFile, suffix));
    delete items;
  }
  delete list;
 
  
  DataTweaker tweak;
  if (!tweak.set(tweaks)) return 0;
 
  Interface* multi = openList(listFileName);
  if (!multi) return 0;
  const MultiTreeAccessor* mt = dynamic_cast<const MultiTreeAccessor*>(&multi->accessor());
  if (!mt){
    delete multi;
    return 0;
  } 
  std::vector<MultiTreeAccessor*> filtered_mts = mt->filterComponents(filterList, tweak);
  if (filtered_mts.size() != filterList.size()){
    delete multi;
    return 0;
  } 
  delete multi;
  cout << "Component filtering done!" << endl;
  // The following line should work, but doesn't... so the block of code below replaces it.
  //Interface* filtered_multi = new Interface(*filtered_mt);
  //
  for (unsigned int f = 0; f < filtered_mts.size(); f++) {
    std::vector<TString> files; 
    for (unsigned int i = 0; i < filtered_mts[f]->nAccessors(); i++) {
      files.push_back(((const TreeAccessor*)&filtered_mts[f]->accessor(i))->fileName());
      cout << "Added " << files.back() << endl;
    }
    delete filtered_mts[f];
    Interface* filtered_multi = open(files);
    //
    std::vector<const Interface*> justOne;
    justOne.push_back(filtered_multi);
    Interface* interface = merge(justOne, filterList.fileName(f));
    delete interface;
  }
  
  return true;
}

findEtaPhi()

HistoryIterator Interface::findEtaPhi	(	CaloId	calo,
		short	layer,
		short	iEta,
		short	iPhi,
		short	region = 0 ) const

Definition at line 415 of file Interface.cxx.

{
  for (unsigned int i = 0; i < nChannels(); i++) {
    const CellInfo* info = cellInfo(i);
    if (!info) continue;
    if (!Id::matchCalo(info->calo(), calo)) continue;
    if (info->layer()  != layer) continue;
    if (info->iEta()   != iEta) continue;
    if (info->iPhi()   != iPhi) continue;
    if (info->region() != region) continue;
    return HistoryIterator(*this, i);
  }
  
  return HistoryIterator(*this, end());
}

findFebChannel()

HistoryIterator Interface::findFebChannel	(	CaloId	calo,
		short	feb,
		short	channel ) const

Definition at line 432 of file Interface.cxx.

{
  for (unsigned int i = 0; i < nChannels(); i++) {
    const CellInfo* info = cellInfo(i);
    if (!info) continue;
    if (!Id::matchCalo(info->calo(), calo)) continue;
    if (info->feb()     != feb)     continue;
    if (info->channel() != channel) continue;
    return HistoryIterator(*this, i);
  }
  
  return HistoryIterator(*this, end());
}

findFTSlotChannel()

HistoryIterator Interface::findFTSlotChannel	(	CaloId	calo,
		short	ft,
		short	slot,
		short	channel ) const

Definition at line 447 of file Interface.cxx.

{
  for (unsigned int i = 0; i < nChannels(); i++) {
    const CellInfo* info = cellInfo(i);
    if (!info) continue;
    if (!Id::matchCalo(info->calo(), calo))    continue;
    if (ft >= 0      && info->feedThrough() != ft)      continue;
    if (slot >= 0    && info->slot()        != slot)    continue;
    if (channel >= 0 && info->channel()     != channel) continue;
    return HistoryIterator(*this, i);
  }
  
  return HistoryIterator(*this, end());
}

firstNeighbors()

bool Interface::firstNeighbors	(	unsigned int	hash,
		std::vector< unsigned int > &	hashes,
		short	layer = -2 ) const

Definition at line 746 of file Interface.cxx.

{
  const CellInfo* cell = cellInfo(hash);
  if (!cell) return true;
  if (!Id::matchCalo(cell->calo(), HEC)) return false; // for now!  
  if (layer < 0) return true;
  std::vector<unsigned int> allHashes;
  const std::vector<unsigned int>* cache = m_neighborCache[hash];
  if (cache) 
    allHashes = *cache;
  else {
    if (!neighbors(*cell, 0.15, allHashes)) return false;
    m_neighborCache[hash] = new std::vector<unsigned int>(allHashes);
  }
  for (unsigned int h : allHashes) {
    const CellInfo* info = cellInfo(h);
    if (!info) continue;
    if (info->layer() == layer) hashes.push_back(h);
    delete info;
  }
  return true;
}

fitMonitor()

FitMonitor LArSamples::Interface::fitMonitor ( ) const

inline

Definition at line 116 of file Interface.h.

{ return FitMonitor(*this); }

getCellHistory()

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

virtual

Implements LArSamples::AbsLArCells.

Definition at line 105 of file Interface.cxx.

{ 
  const History* history = accessor().getCellHistory(i);
  if (history) {
    history->setShapeErrorGetter(m_shapeErrorGetter);
    history->setInterface(this);
  }
  return history;
}

getCellInfo()

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

virtual

Reimplemented from LArSamples::AbsLArCells.

Definition at line 136 of file Interface.cxx.

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

getSCHistory()

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

virtual

Implements LArSamples::AbsLArCells.

Definition at line 115 of file Interface.cxx.

{ 
  const History* history = accessor().getSCHistory(i);
  if (history) {
    history->setInterface(this);
  }
  return history;
}

highEData()

bool Interface::highEData	(	double	eCut,
		TArrayI &	hashes,
		TArrayI &	indices ) const

Definition at line 188 of file Interface.cxx.

{
  std::vector<unsigned int> hashV, indexV;
  unsigned int nTot = 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++) {
      const Data& data = *history->data(j);
      nTot++;
      if (data.isDisconnected()) continue;
      if (nTot % 10000 == 0) cout << nTot << endl;
      if (data.energy() > eCut) {
        cout << "E = " << data.energy() << " " << i << " " << j << endl;
        hashV.push_back(i);
        indexV.push_back(j);
      }
    }
  }
  
  hashes.Set(hashV.size());
  indices.Set(indexV.size());
  
  for (unsigned int i = 0; i < hashV.size(); i++) {
    hashes[i] = hashV[i];
    indices[i] = indexV[i];
  }
  
  cout << hashV.size() << "/" << nTot << endl;  
  return true;
}

historySize()

unsigned int LArSamples::Interface::historySize ( unsigned int i ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 57 of file Interface.h.

{ return accessor().historySize(i); }

historySizeSC()

unsigned int LArSamples::Interface::historySizeSC ( unsigned int i ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 58 of file Interface.h.

{ return accessor().historySizeSC(i); }

isValid()

bool Interface::isValid

(

)

const

Definition at line 222 of file Interface.cxx.

{  
  for (unsigned int i = 0; i < nChannels(); i++) {
    const History* history = cellHistory(i);
    if (!history) continue;
    if (!history->isValid()) { 
      cout << "Invalid LArSamplesHistory at hash = " << i << endl;
      return false;
    }
  }
  
  return true;
}

makeTemplate()

Interface * Interface::makeTemplate

(

const TString &

fileName

)

const

Definition at line 398 of file Interface.cxx.

{
  TreeAccessor* newAcc = TreeAccessor::makeTemplate(accessor(), fileName);
  return new Interface(*newAcc);
}

merge() [1/6]

Interface * Interface::merge	(	const Interface &	other,
		const TString &	fileName ) const

Definition at line 237 of file Interface.cxx.

{
  std::vector<const Interface*> interfaces;
  interfaces.push_back(this);
  interfaces.push_back(&other);
  return merge(interfaces, fileName);
}

merge() [2/6]

Interface * Interface::merge	(	const Interface &	other,
		const TString &	fileName,
		const TString &	LBFile ) const

Definition at line 255 of file Interface.cxx.

{
  std::vector<const Interface*> interfaces;
  interfaces.push_back(this);
  interfaces.push_back(&other);
  return merge(interfaces, fileName, LBFile);
}

merge() [3/6]

Interface * Interface::merge ( const std::vector< const Interface * > & interfaces, const TString & fileName )

static

Definition at line 246 of file Interface.cxx.

{  
  std::vector<const Accessor*> accessors;
  for (unsigned int i = 0; i < interfaces.size(); i++)
    accessors.push_back(&interfaces[i]->accessor());
  TreeAccessor* newAccessor = TreeAccessor::merge(accessors, fileName);
  return new Interface(*newAccessor);
}

merge() [4/6]

Interface * Interface::merge ( const std::vector< const Interface * > & interfaces, const TString & fileName, const TString & LBFile )

static

Definition at line 264 of file Interface.cxx.

{  
  std::vector<const Accessor*> accessors;
  for (unsigned int i = 0; i < interfaces.size(); i++)
    accessors.push_back(&interfaces[i]->accessor());
  TreeAccessor* newAccessor = TreeAccessor::merge(accessors, fileName, LBFile);
  return new Interface(*newAccessor);
}

merge() [5/6]

Interface * Interface::merge ( const TString & listFileName, const TString & fileName )

static

Definition at line 274 of file Interface.cxx.

{
  Interface* multi = openList(listFileName);
  if (!multi) return nullptr;
  std::vector<const Interface*> justOne;
  justOne.push_back(multi);
  return merge(justOne, fileName);
}

merge() [6/6]

Interface * Interface::merge ( const TString & listFileName, const TString & fileName, const TString & LBFile )

static

Definition at line 283 of file Interface.cxx.

{
  Interface* multi = openList(listFileName);
  if (!multi) return nullptr;
  std::vector<const Interface*> justOne;
  justOne.push_back(multi);
  return merge(justOne, fileName, LBFile);
}

nChannels()

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

inlinevirtualinherited

Reimplemented in LArSamples::MonitorBase.

Definition at line 34 of file AbsLArCells.h.

{ return Definitions::nChannels; }

nChannelsSC()

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

inlinevirtualinherited

Definition at line 35 of file AbsLArCells.h.

{ return Definitions::nChannelsSC; }

neighbors()

bool Interface::neighbors	(	const CellInfo &	cell,
		double	dRCut,
		std::vector< unsigned int > &	hashes ) const

Definition at line 732 of file Interface.cxx.

{
  for (unsigned int i = 0; i < nChannels(); i++) {
    const CellInfo* otherCell = cellInfo(i);
    if (!otherCell) continue;
    if (cell.position().DeltaR(otherCell->position()) > dRCut) { delete otherCell; continue; }
    //cout << "Adding hash = " << i << " " << otherCell->location(3) << endl;
    hashes.push_back(i);
    delete otherCell;
  }  
  return true;
}

nEvents()

unsigned int LArSamples::Interface::nEvents ( ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 51 of file Interface.h.

{ return accessor().nEvents(); }

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

nFilledChannels()

unsigned int Interface::nFilledChannels

(

)

const

Definition at line 142 of file Interface.cxx.

{
  unsigned int n = 0;
  for (unsigned int i = 0; i < nChannels(); i++)
    n += (accessor().historySize(i) ? 1 : 0);
  return n;
}

nRuns()

unsigned int LArSamples::Interface::nRuns ( ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 52 of file Interface.h.

{ return accessor().nRuns(); }

occupancyMonitor()

OccupancyMonitor LArSamples::Interface::occupancyMonitor ( ) const

inline

Definition at line 114 of file Interface.h.

{ return OccupancyMonitor(*this); }

open() [1/2]

Interface * Interface::open ( const std::vector< TString > & fileNames )

static

Definition at line 42 of file Interface.cxx.

{
  MultiTreeAccessor* accessor = MultiTreeAccessor::open(fileNames);
  return (accessor ? new Interface(*accessor) : nullptr);
}

open() [2/2]

Interface * Interface::open ( const TString & fileName )

static

Definition at line 35 of file Interface.cxx.

{
  TreeAccessor* accessor = TreeAccessor::open(fileName);
  return (accessor ? new Interface(*accessor) : nullptr);
}

openList()

Interface * Interface::openList ( const TString & fileList )

static

Definition at line 49 of file Interface.cxx.

{
  MultiTreeAccessor* accessor = MultiTreeAccessor::openList(fileList);
  return (accessor ? new Interface(*accessor) : nullptr);
}

openWild()

Interface * Interface::openWild ( const TString & wcName )

static

Definition at line 56 of file Interface.cxx.

{
  MultiTreeAccessor* accessor = MultiTreeAccessor::openWild(wcName);
  return (accessor ? new Interface(*accessor) : nullptr);
}

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

printFilledRanges()

void Interface::printFilledRanges

(

unsigned int

skip = 0

)

const

Definition at line 151 of file Interface.cxx.

{
  int i1 = 0, i2 = 0;
  unsigned int nNull = 0;
  unsigned int size = accessor().historySize(0);
  bool ok = (size > 0);
  bool inInterval = ok;
 
  for (unsigned int i = 0; i < nChannels(); i++) {
    size = accessor().historySize(i);
    if (inInterval && nNull == 0) i2 = i; 
    bool ok = (size > 0);
    if (!ok && inInterval) {
      nNull++;
      if (nNull > skip) {
        nNull = 0;
        inInterval = false;
        cout << i1 << "-" << i2 << endl;
      }
    }
    if (ok) nNull = 0;
    if (ok && !inInterval) {
      inInterval = true;
      i1 = i;
    }
  }
 
  if (inInterval) cout << i1 << "-" << nChannels() << endl;
}

refit()

Interface * Interface::refit	(	const TString &	newFileName,
		Chi2Params	pars = DefaultChi2 ) const

Definition at line 405 of file Interface.cxx.

{
  FilterParams f;
  DataTweaker tw;
  tw.setRefit(true);
  tw.setFitParams(pars);
  return filter(f, tw, newFileName);
}

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

runData()

const RunData * LArSamples::Interface::runData ( unsigned int i ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 55 of file Interface.h.

{ return accessor().runData(i); }

Scan() [1/2]

bool Interface::Scan	(	const TString &	vars,
		CombinationType	comb,
		const TString &	sel = "",
		const TString &	ranges = "",
		unsigned int	verbosity = 1 ) const

Definition at line 578 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return 0;
  return m.dump(vars, comb, f, ranges, verbosity);
}

Scan() [2/2]

bool Interface::Scan	(	const TString &	vars,
		const TString &	sel = "",
		unsigned int	verbosity = 1 ) const

Definition at line 569 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return 0;
  return m.dump(vars, f, verbosity);
}

setShapeError() [1/2]

void Interface::setShapeError

(

const TString &

fileName

)

Definition at line 90 of file Interface.cxx.

{
  setShapeErrorGetter(new TreeShapeErrorGetter(fileName));
}

setShapeError() [2/2]

void Interface::setShapeError

(

double

k

)

Definition at line 81 of file Interface.cxx.

{
  if (m_ownShapeErrorGetter) delete m_shapeErrorGetter;
  if (k == 0) { m_shapeErrorGetter = nullptr; m_ownShapeErrorGetter = false; return; }
  m_shapeErrorGetter = new UniformShapeErrorGetter(k);
  m_ownShapeErrorGetter = true;
}

setShapeErrorGetter()

void Interface::setShapeErrorGetter

(

const AbsShapeErrorGetter *

err

)

Definition at line 73 of file Interface.cxx.

{
  if (m_ownShapeErrorGetter) delete m_shapeErrorGetter;
  m_shapeErrorGetter = err;
  m_ownShapeErrorGetter = false;
}

shapeErrorGetter()

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

inline

Definition at line 123 of file Interface.h.

{ return m_shapeErrorGetter; }

Show() [1/2]

bool Interface::Show	(	const TString &	sel = "",
		unsigned int	verbosity = 1 ) const

Definition at line 595 of file Interface.cxx.

{
  FilterParams f;
  if (!f.set(sel)) return false;
  for (unsigned int i = 0; i < nChannels(); i++) {  
    const History* history = pass(i, f);
    if (!history) continue;
    History* filtered = history->filter(sel);
    TString hDesc = filtered->description(verbosity);
    delete filtered;
    if (hDesc == "") continue;
    cout << Form("Hash = %-5d : ", i) << hDesc
         << "-----------------------------------------------------------------------------" 
         << endl;
  }
  return true;
}

Show() [2/2]

bool Interface::Show	(	unsigned int	hash,
		unsigned int	verbosity = 1 ) const

Definition at line 587 of file Interface.cxx.

{
  const History* history = cellHistory(hash);
  if (!history) return false;
  cout << history->description(verbosity) << endl;
  return true;
}

ShowEvents()

bool Interface::ShowEvents	(	const TString &	sel = "",
		unsigned int	verbosity = 1 ) const

Definition at line 614 of file Interface.cxx.

{
  FilterParams f;
  if (!f.set(sel)) return false;
  std::map< std::pair<unsigned int, unsigned int>, unsigned int > eventCells;
  std::map< std::pair<unsigned int, unsigned int>, double > eventEnergy;
 
  if (verbosity & 8) {
    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++) {
        std::pair<unsigned int, unsigned int> ev = std::make_pair(history->data(j)->run(), history->data(j)->event());
        eventCells[ev]++;
        eventEnergy[ev] += history->data(j)->energy();
      }
    }
  }
  for (unsigned int i = 0; i < nEvents(); i++) {
    const EventData* evtData = eventData(i);
    std::pair<unsigned int, unsigned int> id(evtData->run(), evtData->event());
    TString printout = Form("%d : ", i) + evtData->description(verbosity);
    if (verbosity & 8)
      printout += Form(" : %6d LAr hits, %7.1f MeV", eventCells[id], eventEnergy[id]);
    cout << printout << endl;
  }
  
  return true;
}

ShowRuns()

bool Interface::ShowRuns

(

unsigned int

verbosity = 1

)

const

Definition at line 646 of file Interface.cxx.

{
  std::map<unsigned int, unsigned int> events;
 
  if (verbosity & 8) {
    for (unsigned int i = 0; i < nEvents(); i++)
      events[eventData(i)->run()]++;
  }
  
  for (unsigned int i = 0; i < nRuns(); i++) {
    const RunData* rData = runData(i);
    TString printout = rData->description(verbosity);
    if (verbosity & 8)
      printout += Form(" : %6d events", events[rData->run()]);
    cout << printout << endl;
  }
  
  return true;
}

ShowStats()

bool Interface::ShowStats	(	const TString &	varList,
		const TString &	sel = "",
		bool	withErrors = false ) const

Definition at line 667 of file Interface.cxx.

{
  MonitorBase m(*this);
  FilterParams f;
  if (!f.set(sel)) return 0;
 
  std::vector<TString> vars;
  std::vector<DataFuncSet> funcs;
  std::vector<DataFuncArgs> args;
  if (!MonitorBase::parseVariables(varList, vars, funcs, args)) {
    cout << "Invalid variable specification " << varList << endl;
    return 0;
  }
  
  TVectorD mean(vars.size()), meanErr(vars.size());
  TMatrixD covMatrix(vars.size(), vars.size()), covMatrixErr(vars.size(), vars.size());
  if (!m.statParams(funcs, args, mean, meanErr, covMatrix, covMatrixErr, f)) return false;
  
  if (!withErrors) {
    cout << "---------------------------" << endl;
    for (unsigned int i = 0; i < vars.size(); i++)
      cout << Form("| %10s | %-9.4g |", vars[i].Data(), mean(i)) << endl;
    cout << "---------------------------" << endl << endl;
 
    cout << "|            |";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "            |";
    cout << endl << "--------------";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "-------------";
    cout << endl;
    for (unsigned int i1 = 0; i1 < vars.size(); i1++) {
      cout << Form("| %10s |", vars[i1].Data());
      for (unsigned int i2 = 0; i2 < vars.size(); i2++)
        cout << Form(" %-9.4g |", covMatrix(i1, i2));
      cout << endl;
    }
    cout << "--------------";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "-------------";
    cout << endl;
  }
  else {
    cout << "---------------------------" << endl;
    for (unsigned int i = 0; i < vars.size(); i++)
      cout << Form("| %10s | %-9.4g +/- %-9.4g |", vars[i].Data(), mean(i), meanErr(i)) << endl;
    cout << "---------------------------" << endl << endl;
 
    cout << "|            |";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "                          |";
    cout << endl << "--------------";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "---------------------------";
    cout << endl;
    for (unsigned int i1 = 0; i1 < vars.size(); i1++) {
      cout << Form("| %10s |", vars[i1].Data());
      for (unsigned int i2 = 0; i2 < vars.size(); i2++)
        cout << Form(" %-9.4g +/- %-9.4g  |", covMatrix(i1, i2), covMatrixErr(i1, i2));
      cout << endl;
    }
    cout << "--------------";
    for (unsigned int i = 0; i < vars.size(); i++) cout << "---------------------------";
    cout << endl;
  }
  
  return true;
}

size()

unsigned int Interface::size

(

)

const

Definition at line 96 of file Interface.cxx.

{
  unsigned int n = 0;
  for (unsigned int i = 0; i < nChannels(); i++)
    n += accessor().historySize(i);
  return n;
}

writeToFile()

bool LArSamples::Interface::writeToFile ( const TString & fileName ) const

inlinevirtual

Implements LArSamples::Accessor.

Definition at line 125 of file Interface.h.

{ return accessor().writeToFile(fileName); }

Member Data Documentation

m_accessor

const Accessor* LArSamples::Interface::m_accessor

private

Definition at line 136 of file Interface.h.

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_neighborCache

std::vector<std::vector<unsigned int>*> LArSamples::Interface::m_neighborCache

mutableprivate

Definition at line 140 of file Interface.h.

m_neighborHistories

std::vector<const History*> LArSamples::Interface::m_neighborHistories

mutableprivate

Definition at line 142 of file Interface.h.

m_neighborHistoryPos

std::vector<unsigned int> LArSamples::Interface::m_neighborHistoryPos

mutableprivate

Definition at line 141 of file Interface.h.

m_ownShapeErrorGetter

bool LArSamples::Interface::m_ownShapeErrorGetter

private

Definition at line 138 of file Interface.h.

m_pos

unsigned int LArSamples::AbsLArCells::m_pos

mutableprivateinherited

Definition at line 54 of file AbsLArCells.h.

m_shapeErrorGetter

const AbsShapeErrorGetter* LArSamples::Interface::m_shapeErrorGetter

private

Definition at line 137 of file Interface.h.

---

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

• Interface.h
• Interface.cxx