LArSamples::MultiTreeAccessor Class Reference

#include <MultiTreeAccessor.h>

Inheritance diagram for LArSamples::MultiTreeAccessor:

Collaboration diagram for LArSamples::MultiTreeAccessor:

Public Member Functions
	MultiTreeAccessor (std::vector< std::unique_ptr< const TreeAccessor > > &&accessors)
	Constructor.
virtual	~MultiTreeAccessor ()
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
bool	writeToFile (const TString &fileName) const
const TreeAccessor &	accessor (unsigned int i) const
unsigned int	nAccessors () const
std::vector< MultiTreeAccessor * >	filterComponents (const FilterList &filterList, const DataTweaker &tweaker) 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	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 std::unique_ptr< MultiTreeAccessor >	open (const std::vector< TString > &files)
static std::unique_ptr< MultiTreeAccessor >	openList (const TString &fileList)
static std::unique_ptr< MultiTreeAccessor >	openWild (const TString &wcName)

Protected Member Functions
const History *	getCellHistory (unsigned int i) const
const History *	getSCHistory (unsigned int i) const
const CellInfo *	getCellInfo (unsigned int i) const

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

Detailed Description

Definition at line 27 of file MultiTreeAccessor.h.

Constructor & Destructor Documentation

MultiTreeAccessor()

MultiTreeAccessor::MultiTreeAccessor

(

std::vector< std::unique_ptr< const TreeAccessor > > &&

accessors

)

Constructor.

Definition at line 26 of file MultiTreeAccessor.cxx.

  : m_accessors(std::move(accessors))
{
}

~MultiTreeAccessor()

MultiTreeAccessor::~MultiTreeAccessor ( )

virtual

Definition at line 101 of file MultiTreeAccessor.cxx.

{
}

Member Function Documentation

accessor()

const TreeAccessor & LArSamples::MultiTreeAccessor::accessor ( unsigned int i ) const

inline

Definition at line 51 of file MultiTreeAccessor.h.

{ return *m_accessors[i]; }

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

eventData()

const EventData * MultiTreeAccessor::eventData ( unsigned int i ) const

virtual

Implements LArSamples::Accessor.

Definition at line 106 of file MultiTreeAccessor.cxx.

{
  unsigned int nEventsSoFar = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    unsigned int n = accessor->nEvents();
    if (i < nEventsSoFar + n) return accessor->eventData(i - nEventsSoFar);
    nEventsSoFar += n;
  }
  return nullptr;
}

filterComponents()

std::vector< MultiTreeAccessor * > MultiTreeAccessor::filterComponents	(	const FilterList &	filterList,
		const DataTweaker &	tweaker ) const

Definition at line 268 of file MultiTreeAccessor.cxx.

{
  std::vector< std::vector<std::unique_ptr<const TreeAccessor> > > filteredAccessors(filterList.size());
  
  for (unsigned int i = 0; i < nAccessors(); i++) {
    const TreeAccessor* treeAcc = dynamic_cast<const TreeAccessor*>(&accessor(i));
    cout << "Processing data " << i << " of " << nAccessors();
    if (treeAcc) cout << " (fileName = " << treeAcc->fileName() << ")";
    cout << endl;
    FilterList thisFilterList;
    //
    for (unsigned int f = 0; f < filterList.size(); f++) {
      std::string pathname = (string)filterList.fileName(f);
      if( pathname.find("eos/atlas/") < pathname.length() ){
    int nslashes = 0, slpos = 0;
    for( int k1 = 0; k1 < (int)pathname.length(); k1++ ){
      if( nslashes > 2 )break;
      if( pathname[k1] != '/' )continue;
      nslashes++;
      slpos = k1;
    }
    pathname.resize( slpos );
      }
      TString thisFN = Form("%s_filter%d.root", pathname.c_str(), i );//filterList.fileName(f).Data(), i);
      thisFilterList.add(filterList.filterParams(f), thisFN);
    }
    std::vector<std::unique_ptr<TreeAccessor> > filteredTreeAccs = TreeAccessor::filter(accessor(i), thisFilterList, tweaker);
    if (filteredTreeAccs.size() != filterList.size()) {
      cout << "Filtering failed, exiting" << endl;
      return std::vector<MultiTreeAccessor*>();
    }
    for (unsigned int f = 0; f < filteredTreeAccs.size(); f++) filteredAccessors[f].push_back(std::move(filteredTreeAccs[f]));
  }
 
  std::vector<MultiTreeAccessor*> result;
  for (unsigned int f = 0; f < filteredAccessors.size(); f++) result.push_back(new MultiTreeAccessor(std::move(filteredAccessors[f])));
  return result;
}

getCellHistory()

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

protectedvirtual

Implements LArSamples::AbsLArCells.

Definition at line 171 of file MultiTreeAccessor.cxx.

{ 
  std::unique_ptr<CellInfo> cellInfo;
  std::vector<const Data*> allData;
  std::vector<const EventData*> allEventData;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    std::unique_ptr<const History> thisHistory (accessor->getCellHistory(i));
    if (!thisHistory) continue;
    if (!cellInfo) {
      cellInfo = std::make_unique<CellInfo>(*thisHistory->cellInfo());
    }
    const std::vector<const EventData*>& thisEventData = thisHistory->eventData();
    std::map<const EventData*, const EventData*> eventMap;
    for (const EventData* event : thisEventData) {
      if (eventMap.find(event) != eventMap.end()) continue;
      auto newED = new EventData(*event);
      eventMap[event] = newED;
      allEventData.push_back(newED);
    }
   
    for (unsigned int j = 0; j < thisHistory->nData(); j++) {
      allData.push_back(new Data(*thisHistory->data(j), eventMap[&thisHistory->data(j)->eventData()], nullptr, -1));
      if (!cellInfo->shape(thisHistory->data(j)->gain())) {
        const ShapeInfo* thisShape = thisHistory->cellInfo()->shape(thisHistory->data(j)->gain());
        cellInfo->setShape(thisHistory->data(j)->gain(), thisShape ? new ShapeInfo(*thisShape) : nullptr);
      }
    }
  }
  //data are copied from cellInfo into History member variable
  auto * h  = cellInfo ? new History(allData, *cellInfo, allEventData, i): nullptr;
  return h; 
}

getCellInfo()

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

protectedvirtual

Reimplemented from LArSamples::AbsLArCells.

Definition at line 246 of file MultiTreeAccessor.cxx.

{
  resetCache();
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    const HistoryContainer* cont = accessor->historyContainer(i);
    if (cont && cont->cellInfo()) return new CellInfo(*cont->cellInfo());
  }
  return nullptr;
}

getSCHistory()

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

protectedvirtual

Implements LArSamples::AbsLArCells.

Definition at line 204 of file MultiTreeAccessor.cxx.

{ 
  std::unique_ptr<CellInfo> cellInfo;
  std::vector<const Data*> allData;
  std::vector<const EventData*> allEventData;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    //cout << "---> Getting history for a treeAccessor..." << endl; 
    std::unique_ptr<const History> thisHistory (accessor->getSCHistory(i));
    //cout << "---> done Getting history for a treeAccessor..." << endl; 
    if (!thisHistory) continue;
    if (!cellInfo) {
      cellInfo = std::make_unique<CellInfo>(*thisHistory->cellInfo());
    //cout << "---> done new cell info" << endl; 
    }
    //cout << "---> Creating new event data N = " << thisHistory->eventData().size() << endl; 
    const std::vector<const EventData*>& thisEventData = thisHistory->eventData();
    std::map<const EventData*, const EventData*> eventMap;
    for (const EventData* event : thisEventData) {
      if (eventMap.find(event) != eventMap.end()) continue;
      auto newED = std::make_unique<EventData>(*event);
      eventMap[event] = newED.get();
      allEventData.push_back(newED.release());
    }
    //cout << "---> Creating new data N = " << thisHistory->nData() << endl; 
   
    for (unsigned int ii = 0; ii < thisHistory->nData(); ii++) {
     //cout << "------> Creating new data " << i << endl; 
      allData.push_back(new Data(*thisHistory->data(ii), eventMap[&thisHistory->data(ii)->eventData()], nullptr, -1));
     //cout << "------> done Creating new data " << i << endl; 
      if (!cellInfo->shape(thisHistory->data(ii)->gain())) {
        const ShapeInfo* thisShape = thisHistory->cellInfo()->shape(thisHistory->data(ii)->gain());
        cellInfo->setShape(thisHistory->data(ii)->gain(), thisShape ? new ShapeInfo(*thisShape) : nullptr);
      }
     //cout << "------> done shape " << i << endl; 
    }
     //cout << "---> done Creating new data, deleting treeAcc history" << endl; 
  }
  //cout << "--->returning new history..." << endl; 
  return (cellInfo ? new History(allData, *cellInfo, allEventData, i) : nullptr); 
}

historySize()

unsigned int MultiTreeAccessor::historySize ( unsigned int i ) const

virtual

Implements LArSamples::Accessor.

Definition at line 148 of file MultiTreeAccessor.cxx.

{
  resetCache();
  unsigned int size = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    const HistoryContainer* cont = accessor->historyContainer(i);
    if (cont) size += cont->nDataContainers();
  }
  return size;
}

historySizeSC()

unsigned int MultiTreeAccessor::historySizeSC ( unsigned int i ) const

virtual

Implements LArSamples::Accessor.

Definition at line 159 of file MultiTreeAccessor.cxx.

{
  resetCache();
  unsigned int size = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    const HistoryContainer* cont = accessor->historyContainerSC(i);
    if (cont) size += cont->nDataContainers();
  }
  return size;
}

nAccessors()

unsigned int LArSamples::MultiTreeAccessor::nAccessors ( ) const

inline

Definition at line 52 of file MultiTreeAccessor.h.

{ return m_accessors.size(); }

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

nEvents()

unsigned int MultiTreeAccessor::nEvents ( ) const

virtual

Implements LArSamples::Accessor.

Definition at line 130 of file MultiTreeAccessor.cxx.

{
  unsigned int n = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors)
     n += accessor->nEvents();
  return n;
}

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

nRuns()

unsigned int MultiTreeAccessor::nRuns ( ) const

virtual

Implements LArSamples::Accessor.

Definition at line 139 of file MultiTreeAccessor.cxx.

{
  unsigned int n = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors)
     n += accessor->nRuns();
  return n;
}

open()

std::unique_ptr< MultiTreeAccessor > MultiTreeAccessor::open ( const std::vector< TString > & files )

static

Definition at line 32 of file MultiTreeAccessor.cxx.

{
  std::vector<std::unique_ptr<const TreeAccessor> > accessors;
  for (const TString& fileName : files) {
    std::unique_ptr<const TreeAccessor> accessor = TreeAccessor::open(fileName);
    if (!accessor) {
      cout << "Skipping invalid file " << fileName << endl;
      continue;
    }
    accessors.push_back(std::move(accessor));
  }
  
  if (accessors.empty()) return nullptr;
  return std::make_unique<MultiTreeAccessor>(std::move(accessors));
}

openList()

std::unique_ptr< MultiTreeAccessor > MultiTreeAccessor::openList ( const TString & fileList )

static

Definition at line 49 of file MultiTreeAccessor.cxx.

{
  std::ifstream f(fileList);
  if (!f) {
    cout << "file " << fileList << " not accessible" << endl;
    return nullptr;
  }
 
  std::string fileName;  
  unsigned int i = 0;
 
  std::vector<std::unique_ptr<const TreeAccessor> > accessors;
  
  while (f >> fileName) {
    std::unique_ptr<const TreeAccessor> accessor = TreeAccessor::open(fileName.c_str());
    if (!accessor) {
      cout << "Skipping invalid file " << fileName << endl;
      continue;
    }
    cout << std::setw(2) << ++i << " - " << fileName << endl;
    accessors.push_back(std::move(accessor));
  }
  
  if (accessors.empty()) return nullptr;
  return std::make_unique<MultiTreeAccessor>(std::move(accessors));
}

openWild()

std::unique_ptr< MultiTreeAccessor > MultiTreeAccessor::openWild ( const TString & wcName )

static

Definition at line 77 of file MultiTreeAccessor.cxx.

{
  // Piggyback on TChain wildcarding feature...
  TChain chain("");
  chain.Add(wcName);
  
  std::vector<std::unique_ptr<const TreeAccessor> > accessors;
  
  for (int i = 0; i < chain.GetListOfFiles()->GetEntries(); i++) {
    std::string fileName = chain.GetListOfFiles()->At(i)->GetTitle();
    std::unique_ptr<const TreeAccessor> accessor = TreeAccessor::open(fileName.c_str());
    if (!accessor) {
      cout << "Skipping invalid file " << fileName << endl;
      continue;
    }
    cout << std::setw(2) << i+1 << " - " << fileName <<  " , nEvents = " << accessor->nEvents() << ", nRuns = " << accessor->nRuns() << endl;
    accessors.push_back(std::move(accessor));
  }
  
  if (accessors.empty()) return nullptr;
  return std::make_unique<MultiTreeAccessor>(std::move(accessors));
}

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

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.

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

runData()

const RunData * MultiTreeAccessor::runData ( unsigned int i ) const

virtual

Implements LArSamples::Accessor.

Definition at line 118 of file MultiTreeAccessor.cxx.

{ 
  unsigned int nRunsSoFar = 0;
  for (const std::unique_ptr<const TreeAccessor>& accessor : m_accessors) {
    unsigned int n = accessor->nRuns();
    if (i < nRunsSoFar + n) return accessor->runData(i - nRunsSoFar);
    nRunsSoFar += n;
  }
  return nullptr;
}

writeToFile()

bool MultiTreeAccessor::writeToFile ( const TString & fileName ) const

virtual

Implements LArSamples::Accessor.

Definition at line 257 of file MultiTreeAccessor.cxx.

{
  std::vector<const Accessor*> accessors;
  for (unsigned int i = 0; i < m_accessors.size(); i++)
    accessors.push_back(m_accessors[i].get());
  cout << "Merging data..." << endl;
  std::unique_ptr<TreeAccessor> singleContainer = TreeAccessor::merge(accessors, fileName);
  return !!singleContainer;
}

Member Data Documentation

m_accessors

std::vector<std::unique_ptr<const TreeAccessor> > LArSamples::MultiTreeAccessor::m_accessors

private

Definition at line 64 of file MultiTreeAccessor.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_pos

unsigned int LArSamples::AbsLArCells::m_pos

mutableprivateinherited

Definition at line 54 of file AbsLArCells.h.

---

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

• MultiTreeAccessor.h
• MultiTreeAccessor.cxx