lhef::H5Handler Class Reference

#include <LH5Reader.h>

Public Member Functions
void	counter (int, int)
std::vector< lheh5::Particle >	getEvent ()
std::pair< lheh5::EventHeader, std::vector< lheh5::Particle > >	getEventProperties ()
lheh5::EventHeader	getHeader ()
	H5Handler (const std::string &fileNameIn)
void	readBlock ()
virtual	~H5Handler ()

Public Attributes
HighFive::Group	_event
lheh5::Events	_events1
lheh5::Events2	_events2
unsigned int long	_eventsRead
HighFive::Group	_index
HighFive::Group	_init
HighFive::Group	_particle
HighFive::Group	_procInfo
std::unique_ptr< HighFive::File >	h5file
bool	indexStatus
int	npLO
int	npNLO

Private Attributes
int	_blocksRead
int	_eventsInBlock
unsigned int long	_eventsTotal
int	_formatType

Detailed Description

Definition at line 16 of file LH5Reader.h.

Constructor & Destructor Documentation

H5Handler()

lhef::H5Handler::H5Handler

(

const std::string &

fileNameIn

)

Definition at line 70 of file LH5Reader.cc.

References _blocksRead, _eventsInBlock, _eventsRead, _eventsTotal, _formatType, _index, _procInfo, Exception, h5file, indexStatus, npLO, npNLO, and mps_update::status.

      : h5file(new HighFive::File(fileNameIn)),
        indexStatus(h5file->exist("/index")),
        _index(h5file->getGroup(indexStatus ? "index" : "event")),
        _particle(h5file->getGroup("particle")),
        _event(h5file->getGroup("event")),
        _init(h5file->getGroup("init")),
        _procInfo(h5file->getGroup("procInfo")) {
    hid_t dspace;
    _formatType = 1;

    if (indexStatus) {
      dspace = H5Dget_space(h5file->getDataSet("index/start").getId());
      _formatType = 2;
    } else {
      _index = h5file->getGroup("event");
      dspace = H5Dget_space(h5file->getDataSet("event/start").getId());
      _formatType = 1;
    }

    _eventsTotal = H5Sget_simple_extent_npoints(dspace);
    _eventsRead = 0;
    _eventsInBlock = 6400 / 4;  // configurable parameter??
    if (_eventsTotal % _eventsInBlock != 0)
      throw cms::Exception("ReadError") << "block size does not match HDF5 file" << std::endl;
    _blocksRead = 0;

    // Check if the file contains the npLO, npNLO information
    bool status = h5file->exist("/procInfo/npLO");
    if (status) {
      HighFive::DataSet _npLO = _procInfo.getDataSet("npLO");
      _npLO.read(npLO);
      HighFive::DataSet _npNLO = _procInfo.getDataSet("npNLO");
      _npNLO.read(npNLO);
    } else {
      npLO = 1;
      npNLO = 0;
    }
  };

~H5Handler()

virtual lhef::H5Handler::~H5Handler ( )

inlinevirtual

Definition at line 19 of file LH5Reader.h.

{}

Member Function Documentation

counter()

void lhef::H5Handler::counter	(	int	firstEventIn,
		int	maxEventsIn
	)

Definition at line 110 of file LH5Reader.cc.

References _blocksRead, _eventsInBlock, _eventsRead, _eventsTotal, Exception, and createfilelist::int.

Referenced by counter.Counters::__getitem__().

                                                                       {
    if (maxEventsIn > 0 && firstEventIn > maxEventsIn)
      throw cms::Exception("ConfigurationError") << "\" firstEvent > maxEvents \"" << std::endl;
    // Reset maximum number of events to read
    if (_blocksRead == 0 && maxEventsIn >= 0 && (unsigned long int)maxEventsIn < _eventsTotal)
      _eventsTotal = (unsigned long int)maxEventsIn;
    // If there are multiple files, assume you only want to jump through the first file
    if (firstEventIn > 0 && _blocksRead > 0)
      _eventsRead = firstEventIn - 1;
    // Set blocks read to be in the correct place
    _blocksRead = _eventsRead / _eventsInBlock;
  }

getEvent()

std::vector< lheh5::Particle > lhef::H5Handler::getEvent

(

)

Definition at line 133 of file LH5Reader.cc.

References _blocksRead, _events1, _eventsInBlock, _eventsRead, _eventsTotal, lheh5::Events::mkEvent(), and readBlock().

                                               {
    std::vector<lheh5::Particle> _vE;
    if (_eventsRead > _eventsTotal - 1)
      return std::vector<lheh5::Particle>();
    int checkEvents = (_blocksRead)*_eventsInBlock - 1;
    _vE = _events1.mkEvent(_eventsRead);
    _eventsRead++;
    if (checkEvents >= 0 && _eventsRead > (unsigned long int)checkEvents)
      readBlock();
    return _vE;
  }

getEventProperties()

std::pair< lheh5::EventHeader, std::vector< lheh5::Particle > > lhef::H5Handler::getEventProperties

(

)

Definition at line 150 of file LH5Reader.cc.

References _blocksRead, _events1, _eventsInBlock, _eventsRead, _eventsTotal, lheh5::Events::mkEvent(), lheh5::Events::mkEventHeader(), lheh5::EventHeader::nparticles, and readBlock().

                                                                                     {
    std::vector<lheh5::Particle> _vE;
    lheh5::EventHeader _h;
    _h.nparticles = -1;
    if (_eventsRead > _eventsTotal - 1)
      return std::make_pair(_h, _vE);
    int checkEvents = (_blocksRead)*_eventsInBlock - 1;
    _vE = _events1.mkEvent(_eventsRead % _eventsInBlock);
    _h = _events1.mkEventHeader(_eventsRead % _eventsInBlock);
    _eventsRead++;
    if (checkEvents >= 0 && _eventsRead > (unsigned long int)checkEvents)
      readBlock();
    return std::make_pair(_h, _vE);
  }

getHeader()

lheh5::EventHeader lhef::H5Handler::getHeader

(

)

Definition at line 145 of file LH5Reader.cc.

References _events1, _eventsRead, and lheh5::Events::mkEventHeader().

                                        {
    // fragile, must be called before getEvent
    return _events1.mkEventHeader(_eventsRead);
  }

readBlock()

void lhef::H5Handler::readBlock

(

)

Definition at line 123 of file LH5Reader.cc.

References _blocksRead, _event, _events1, _eventsInBlock, _eventsTotal, _index, _particle, createIOVlist::nEvents, and lheh5::readEvents().

Referenced by getEvent(), and getEventProperties().

                            {
    // Start counting at 0
    size_t iStart = _blocksRead * _eventsInBlock;
    size_t nEvents = _eventsInBlock;
    if ((unsigned long int)(iStart + nEvents) > _eventsTotal)
      return;
    _events1 = readEvents(_index, _particle, _event, iStart, nEvents);
    _blocksRead++;
  }

Member Data Documentation

_blocksRead

int lhef::H5Handler::_blocksRead

private

Definition at line 37 of file LH5Reader.h.

Referenced by counter(), getEvent(), getEventProperties(), H5Handler(), and readBlock().

_event

HighFive::Group lhef::H5Handler::_event

Definition at line 24 of file LH5Reader.h.

Referenced by readBlock().

_events1

lheh5::Events lhef::H5Handler::_events1

Definition at line 28 of file LH5Reader.h.

Referenced by getEvent(), getEventProperties(), getHeader(), and readBlock().

_events2

lheh5::Events2 lhef::H5Handler::_events2

Definition at line 27 of file LH5Reader.h.

_eventsInBlock

int lhef::H5Handler::_eventsInBlock

private

Definition at line 35 of file LH5Reader.h.

Referenced by counter(), getEvent(), getEventProperties(), H5Handler(), and readBlock().

_eventsRead

unsigned int long lhef::H5Handler::_eventsRead

Definition at line 26 of file LH5Reader.h.

Referenced by counter(), getEvent(), getEventProperties(), getHeader(), and H5Handler().

_eventsTotal

unsigned int long lhef::H5Handler::_eventsTotal

private

Definition at line 34 of file LH5Reader.h.

Referenced by counter(), getEvent(), getEventProperties(), H5Handler(), and readBlock().

_formatType

int lhef::H5Handler::_formatType

private

Definition at line 36 of file LH5Reader.h.

Referenced by H5Handler().

_index

HighFive::Group lhef::H5Handler::_index

Definition at line 24 of file LH5Reader.h.

Referenced by ntupleDataFormat._Object::__getattr__(), python.cmstools.EventTree::__getitem__(), python.cmstools.EventTree::__iter__(), python.cmstools.EventBranch::__readData(), python.cmstools.EventTree::__setBranchIndicies(), ntupleDataFormat.TrackingParticle::_nMatchedSeeds(), ntupleDataFormat._SimHitMatchAdaptor::_nMatchedSimHits(), ntupleDataFormat._TrackingParticleMatchAdaptor::_nMatchedTrackingParticles(), ntupleDataFormat.TrackingParticle::_nMatchedTracks(), html.HtmlReport::addNote(), ntupleDataFormat.TrackingVertex::daughterTrackingParticles(), ntupleDataFormat.TrackingParticle::decayVertices(), ntupleDataFormat._DetIdStrAdaptor::detIdStr(), H5Handler(), ntupleDataFormat._Object::index(), python.cmstools.EventTree::index(), ntupleDataFormat.Seed::indexWithinAlgo(), ntupleDataFormat._Object::isValid(), ntupleDataFormat._DetIdStrAdaptor::layerStr(), ntupleDataFormat.InvalidHit::layerStr(), ntupleDataFormat.TrackingParticle::matchedSeedInfos(), ntupleDataFormat._SimHitMatchAdaptor::matchedSimHitInfos(), ntupleDataFormat.TrackingParticle::matchedTrackInfos(), ntupleDataFormat._TrackingParticleMatchAdaptor::matchedTrackingParticleInfos(), ntupleDataFormat.GluedHit::monoHit(), ntupleDataFormat.TrackingVertex::nDaughterTrackingParticles(), ntupleDataFormat.SimHit::nRecHits(), ntupleDataFormat._HitObject::nseeds(), ntupleDataFormat.GluedHit::nseeds(), ntupleDataFormat.TrackingVertex::nSourceTrackingParticles(), ntupleDataFormat._HitObject::ntracks(), ntupleDataFormat.Vertex::nTracks(), ntupleDataFormat.TrackingParticle::parentVertex(), readBlock(), ntupleDataFormat.SeedMatchInfo::seed(), ntupleDataFormat.Track::seed(), ntupleDataFormat._HitObject::seeds(), ntupleDataFormat.GluedHit::seeds(), python.cmstools.EventBranch::setIndex(), ntupleDataFormat.SimHitMatchInfo::simHit(), ntupleDataFormat.TrackingVertex::sourceTrackingParticles(), ntupleDataFormat.GluedHit::stereoHit(), ntupleDataFormat.Seed::track(), ntupleDataFormat.SimHit::trackingParticle(), ntupleDataFormat._HitObject::tracks(), ntupleDataFormat.Vertex::tracks(), ntupleDataFormat.Track::vertex(), and html.HtmlReport::write().

_init

HighFive::Group lhef::H5Handler::_init

Definition at line 24 of file LH5Reader.h.

Referenced by ElectronMVAID.ElectronMVAID::__call__(), and FWLite.ElectronMVAID::__call__().

_particle

HighFive::Group lhef::H5Handler::_particle

Definition at line 24 of file LH5Reader.h.

Referenced by readBlock().

_procInfo

HighFive::Group lhef::H5Handler::_procInfo

Definition at line 24 of file LH5Reader.h.

Referenced by H5Handler().

h5file

std::unique_ptr<HighFive::File> lhef::H5Handler::h5file

Definition at line 22 of file LH5Reader.h.

Referenced by H5Handler().

indexStatus

bool lhef::H5Handler::indexStatus

Definition at line 23 of file LH5Reader.h.

Referenced by H5Handler().

npLO

int lhef::H5Handler::npLO

Definition at line 25 of file LH5Reader.h.

Referenced by H5Handler().

npNLO

int lhef::H5Handler::npNLO

Definition at line 25 of file LH5Reader.h.

Referenced by H5Handler().