THistSvc Class Reference

#include <src/THistSvc/THistSvc.h>

Inheritance diagram for THistSvc:

Collaboration diagram for THistSvc:

Classes
class	GlobalDirectoryRestore
struct	THistID

Public Member Functions
StatusCode	initialize () override
StatusCode	reinitialize () override
StatusCode	finalize () override
StatusCode	regHist (const std::string &name) override
StatusCode	regHist (const std::string &name, TH1 *) override
StatusCode	regHist (const std::string &name, TH2 *) override
StatusCode	regHist (const std::string &name, TH3 *) override
StatusCode	getHist (const std::string &name, TH1 *&) const override
StatusCode	getHist (const std::string &name, TH2 *&) const override
StatusCode	getHist (const std::string &name, TH3 *&) const override
StatusCode	regTree (const std::string &name) override
StatusCode	regTree (const std::string &name, TTree *) override
StatusCode	getTree (const std::string &name, TTree *&) const override
StatusCode	regGraph (const std::string &name) override
StatusCode	regGraph (const std::string &name, TGraph *) override
StatusCode	getGraph (const std::string &name, TGraph *&) const override
StatusCode	deReg (TObject *obj) override
StatusCode	deReg (const std::string &name) override
std::vector< std::string >	getHists () const override
std::vector< std::string >	getTrees () const override
std::vector< std::string >	getGraphs () const override
StatusCode	getTHists (TDirectory *td, TList &, bool recurse=false) const override
StatusCode	getTHists (const std::string &name, TList &, bool recurse=false) const override
StatusCode	getTHists (TDirectory *td, TList &tl, bool recurse=false, bool reg=false) override
StatusCode	getTHists (const std::string &name, TList &tl, bool recurse=false, bool reg=false) override
StatusCode	getTTrees (TDirectory *td, TList &, bool recurse=false) const override
StatusCode	getTTrees (const std::string &name, TList &, bool recurse=false) const override
StatusCode	getTTrees (TDirectory *td, TList &tl, bool recurse=false, bool reg=false) override
StatusCode	getTTrees (const std::string &name, TList &tl, bool recurse=false, bool reg=false) override
bool	exists (const std::string &name) const override
	THistSvc (const std::string &name, ISvcLocator *svc)
void	handle (const Incident &) override
StatusCode	io_reinit () override
	callback method to reinitialize the internal state of the component for I/O purposes (e.g. More...
Public Member Functions inherited from extends< BASE, Interfaces >
void *	i_cast (const InterfaceID &tid) const override
	Implementation of IInterface::i_cast. More...
StatusCode	queryInterface (const InterfaceID &ti, void **pp) override
	Implementation of IInterface::queryInterface. More...
std::vector< std::string >	getInterfaceNames () const override
	Implementation of IInterface::getInterfaceNames. More...
	~extends () override=default
	Virtual destructor. More...
void *	i_cast (const InterfaceID &tid) const override
	Implementation of IInterface::i_cast. More...
StatusCode	queryInterface (const InterfaceID &ti, void **pp) override
	Implementation of IInterface::queryInterface. More...
std::vector< std::string >	getInterfaceNames () const override
	Implementation of IInterface::getInterfaceNames. More...
	~extends () override=default
	Virtual destructor. More...
Public Member Functions inherited from extend_interfaces< Interfaces...>
	~extend_interfaces () override=default
	Virtual destructor. More...
	~extend_interfaces () override=default
	Virtual destructor. More...

Protected Member Functions
	~THistSvc () override=default

Private Types
enum	Mode {   READ, WRITE, UPDATE, APPEND,   SHARE, INVALID }
typedef std::map< std::string, THistID >	uidMap
typedef std::multimap< std::string, THistID >	idMap
typedef std::map< TObject *, THistID >	objMap
typedef std::multimap< std::string, std::string >	streamMap

Private Member Functions
template<typename T >
StatusCode	regHist_i (T *hist, const std::string &name)
template<typename T >
StatusCode	getHist_i (const std::string &name, T *&hist, bool quiet=false) const
template<typename T >
StatusCode	readHist_i (const std::string &name, T *&hist) const
StatusCode	readHist (const std::string &name, TH1 *&) const
StatusCode	readHist (const std::string &name, TH2 *&) const
StatusCode	readHist (const std::string &name, TH3 *&) const
StatusCode	readTree (const std::string &name, TTree *&) const
void	updateFiles ()
StatusCode	write ()
StatusCode	connect (const std::string &)
TDirectory *	changeDir (const THistSvc::THistID &hid) const
std::string	dirname (std::string &dir) const
void	removeDoubleSlash (std::string &) const
bool	browseTDir (TDirectory *dir) const
bool	findStream (const std::string &name, std::string &root, std::string &rem, TFile *&file) const
void	parseString (const std::string &id, std::string &root, std::string &rem) const
void	setupInputFile (Property &inputfile)
	call-back method to handle input stream property More...
void	setupOutputFile (Property &outputfile)
	call-back method to handle output stream property More...
void	setupCompressionLevel (Property &cmp)
void	copyFileLayout (TDirectory *, TDirectory *)
	helper function to recursively copy the layout of a TFile into a new TFile More...
void	MergeRootFile (TDirectory *target, TDirectory *source)
StatusCode	rootOpenAction (FILEMGR_CALLBACK_ARGS)
StatusCode	rootOpenErrAction (FILEMGR_CALLBACK_ARGS)

Private Attributes
MsgStream	m_log
StringArrayProperty	m_inputfile
StringArrayProperty	m_outputfile
std::vector< std::string >	m_Rstream
std::vector< std::string >	m_Wstream
IntegerProperty	m_autoSave
IntegerProperty	m_autoFlush
IntegerProperty	m_compressionLevel
IntegerProperty	m_maxFileSize
BooleanProperty	m_print
std::set< std::string >	m_alreadyConnectedInFiles
	list of already connected files. More...
std::set< std::string >	m_alreadyConnectedOutFiles
	list of already connected files. More...
uidMap	m_uids
idMap	m_ids
objMap	m_tobjs
std::map< std::string, std::pair< TFile *, Mode > >	m_files
streamMap	m_fileStreams
std::map< std::string, std::string >	m_sharedFiles
bool	signaledStop = false
bool	m_delayConnect = false
bool	m_okToConnect = false
std::string	m_curstream
IIncidentSvc *	p_incSvc = nullptr
IFileMgr *	p_fileMgr = nullptr

Additional Inherited Members
Public Types inherited from extends< BASE, Interfaces >
using	base_class = extends
	Typedef to this class. More...
using	extend_interfaces_base = extend_interfaces< Interfaces...>
	Typedef to the base of this class. More...
using	base_class = extends
	Typedef to this class. More...
using	extend_interfaces_base = extend_interfaces< Interfaces...>
	Typedef to the base of this class. More...
Public Types inherited from extend_interfaces< Interfaces...>
using	ext_iids = typename Gaudi::interface_list_cat< typename Interfaces::ext_iids...>::type
	take union of the ext_iids of all Interfaces... More...
using	ext_iids = typename Gaudi::interface_list_cat< typename Interfaces::ext_iids...>::type
	take union of the ext_iids of all Interfaces... More...

Detailed Description

Definition at line 26 of file THistSvc.h.

Member Typedef Documentation

typedef std::multimap<std::string, THistID> THistSvc::idMap

private

Definition at line 195 of file THistSvc.h.

typedef std::map<TObject*, THistID> THistSvc::objMap

private

Definition at line 196 of file THistSvc.h.

typedef std::multimap<std::string, std::string> THistSvc::streamMap

private

Definition at line 197 of file THistSvc.h.

typedef std::map<std::string, THistID> THistSvc::uidMap

private

Definition at line 194 of file THistSvc.h.

Member Enumeration Documentation

enum THistSvc::Mode

private

Enumerator
READ
WRITE
UPDATE
APPEND
SHARE
INVALID

Definition at line 105 of file THistSvc.h.

            {
    READ,
    WRITE,
    UPDATE,
    APPEND,
    SHARE,
    INVALID
  };

Constructor & Destructor Documentation

THistSvc::THistSvc	(	const std::string &	name,
		ISvcLocator *	svc
	)

Definition at line 62 of file THistSvc.cpp.

  : base_class(name, svc), m_log(msgSvc(), name )
{    
  declareProperty ("AutoSave", m_autoSave=0 );
  declareProperty ("AutoFlush", m_autoFlush=0 );
  declareProperty ("PrintAll", m_print=false);
  declareProperty ("MaxFileSize", m_maxFileSize=10240,
           "maximum file size in MB. if exceeded, will cause an abort. -1 to never check.");
  declareProperty ("CompressionLevel", m_compressionLevel=1 )->declareUpdateHandler( &THistSvc::setupCompressionLevel, this );
  declareProperty ("Output", m_outputfile )->declareUpdateHandler( &THistSvc::setupOutputFile, this );
  declareProperty ("Input", m_inputfile )->declareUpdateHandler ( &THistSvc::setupInputFile,  this );
}

THistSvc::~THistSvc ( )

overrideprotecteddefault

Member Function Documentation

bool THistSvc::browseTDir ( TDirectory *  dir ) const

private

Definition at line 405 of file THistSvc.cpp.

                                          {

  if (!dir) {
    std::cerr << "TDirectory == 0" << std::endl;
    return false;
  }

  GlobalDirectoryRestore restore;

  dir->cd();


  cout << "-> " << dir->GetPath() << "  "
       << dir->GetListOfKeys()->GetSize() << endl;

  //  TIter nextkey(dir->GetListOfKeys());
  TIter nextkey(dir->GetList());
  while (TKey *key = (TKey*)nextkey()) {

    TObject *obj = key->ReadObj();
    if (!obj) { cout << key->GetName() << " obj==0"<< endl; continue; }
    //    if (obj->IsA()->InheritsFrom("TDirectory")) {
      cout << "  Key: " << key->GetName() << "   "
           << " tit: " << obj->GetTitle() << "   "
           << " (" << key->GetClassName() << ")" << endl;
      //    }
  }

  nextkey = dir->GetListOfKeys();
  while (TKey *key = (TKey*)nextkey()) {

    TObject *obj = key->ReadObj();
    if (!obj) { cout << key->GetName() << " obj==0"<< endl; continue; }
    if (obj->IsA()->InheritsFrom("TDirectory")) {
      TDirectory *tt = dynamic_cast<TDirectory*>(obj);
      browseTDir(tt);
    }
  }

  return true;
}

TDirectory * THistSvc::changeDir ( const THistSvc::THistID &  hid ) const

private

Definition at line 1655 of file THistSvc.cpp.

                                                    {

  string uid = hid.id;
  TFile* file = hid.file;
  string stream, fdir, bdir, dir, id;

  if (file) {
    file->cd("/");
  } else {
    gROOT->cd();
  }

  fdir = uid;
  bdir = dirname(fdir);

  while ( (dir = dirname(fdir)) != "") {
    if (! gDirectory->GetKey(dir.c_str())) {
      gDirectory->mkdir(dir.c_str());
    }
    gDirectory->cd(dir.c_str());
  }

  return gDirectory;

}

StatusCode THistSvc::connect ( const std::string &  ident )

private

Definition at line 1431 of file THistSvc.cpp.

                                        {

  auto loc = ident.find(" ");
  string stream = ident.substr(0,loc);
  char typ(0);
  typedef std::pair<std::string,std::string>      Prop;
  std::vector<Prop> props;
  string filename, db_typ("ROOT");
  int cl(1);

  if (loc != string::npos) {
    using Parser = Gaudi::Utils::AttribStringParser;
    for (auto attrib: Parser(ident.substr(loc + 1))) {
      auto TAG = boost::algorithm::to_upper_copy(attrib.tag);
      auto VAL = boost::algorithm::to_upper_copy(attrib.value);

      if (TAG == "FILE" || TAG == "DATAFILE") {
        filename = attrib.value;
        removeDoubleSlash( filename );
      } else if ( TAG == "OPT" ) {
        if ( VAL == "APPEND" || VAL == "UPDATE" ) {
          typ = 'A';
        } else if ( VAL == "CREATE" || VAL == "NEW" || VAL == "WRITE" ) {
          typ = 'N';
        } else if ( VAL == "RECREATE" ) {
          typ = 'R';
        } else if (VAL == "SHARE") {
          typ = 'S';
        } else if ( VAL == "OLD" || VAL == "READ" ) {
          typ = 'O';
        } else {
          m_log << MSG::ERROR << "Unknown OPT: \"" << attrib.value << "\""
              << endmsg;
          typ = 0;
        }
      } else if (TAG == "TYP") {
        db_typ = std::move(attrib.value);
      } else if (TAG == "CL") {
        cl = std::stoi(attrib.value);
      } else {
        props.emplace_back( attrib.tag, attrib.value);
      }

    }
  }

  if (stream == "temp") {
    m_log << MSG::ERROR << "in JobOption \"" << ident
          << "\": stream name \"temp\" reserved."
          << endmsg;
    return StatusCode::FAILURE;
  }

  if (db_typ != "ROOT") {
    m_log << MSG::ERROR << "in JobOption \"" << ident
          << "\": technology type \"" << db_typ << "\" not supported."
          << endmsg;
    return StatusCode::FAILURE;
  }


  if (m_files.find(stream) != m_files.end()) {
    m_log << MSG::ERROR << "in JobOption \"" << ident
          << "\":\n stream \"" << stream << "\" already connected to file: \""
          << m_files[stream].first->GetName() << "\""
          << endmsg;
    return StatusCode::FAILURE;
  }

  Mode newMode;
  if (typ == 'O') {
    newMode = THistSvc::READ;
  } else if (typ == 'N') {
    newMode = THistSvc::WRITE;
  } else if (typ == 'A') {
    newMode = THistSvc::APPEND;
  } else if (typ == 'R') {
    newMode = THistSvc::UPDATE;
  } else if (typ == 'S') {
    newMode = THistSvc::SHARE;
  } else {
    // something else?
    m_log << MSG::ERROR << "No OPT= specified or unknown access mode in: "
          << ident << endmsg;
    return StatusCode::FAILURE;
  }

  // Is this file already connected to another stream?
  if (m_fileStreams.find(filename) != m_fileStreams.end()) {
    auto fitr = m_fileStreams.equal_range(filename);

    const std::string& oldstream = fitr.first->second;

    const auto& f_info = m_files[oldstream];

    if (newMode != f_info.second) {
      m_log << MSG::ERROR << "in JobOption \"" << ident
            << "\":\n file \"" << filename << "\" already opened by stream: \""
            << oldstream << "\" with different access mode."
            << endmsg;
      return StatusCode::FAILURE;
    } else {
      TFile *f2 = f_info.first;
      m_files[stream] = make_pair(f2,newMode);
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "Connecting stream: \"" << stream
              << "\" to previously opened TFile: \"" << filename << "\""
              << endmsg;
      return StatusCode::SUCCESS;
    }
  }


  IIncidentSvc *pi = nullptr;
  if (service("IncidentSvc",pi).isFailure()) {
    m_log << MSG::ERROR << "Unable to get the IncidentSvc" << endmsg;
    return StatusCode::FAILURE;
  }

  void* vf = nullptr;
  TFile *f = nullptr;

  if (newMode == THistSvc::READ) {
    // old file

    int r = p_fileMgr->open(Io::ROOT,name(), filename,Io::READ,vf,"HIST");

    if (r != 0) {
      m_log << "Unable to open ROOT file " << filename << " for reading"
        << endmsg;
      return StatusCode::FAILURE;
    }


    f = (TFile*) vf;

    // FIX ME!
    pi->fireIncident(FileIncident(name(), "BeginHistFile", filename));


  } else if (newMode == THistSvc::WRITE) {
    // new file. error if file exists

    int r = p_fileMgr->open(Io::ROOT,name(),filename, (Io::WRITE|Io::CREATE|Io::EXCL),
                vf,"HIST");

    if (r != 0) {
      m_log << "Unable to open ROOT file " << filename << " for writing"
        << endmsg;
      return StatusCode::FAILURE;
    }

    f = (TFile*)vf;

  } else if (newMode == THistSvc::APPEND) {
    // update file

    int r = p_fileMgr->open(Io::ROOT,name(),filename, (Io::WRITE | Io::APPEND),
                vf,"HIST");
    if (r != 0) {
      m_log << MSG::ERROR << "unable to open file \"" << filename
            << "\" for appending" << endmsg;
      return StatusCode::FAILURE;
    }

    f = (TFile*) vf;


  } else if (newMode == THistSvc::SHARE) {
    // SHARE file type
    //For SHARE files, all data will be stored in a temp file and will be merged into the target file
    //in write() when finalize(), this help to solve some confliction. e.g. with storegate

    static int ishared = 0;
    string realfilename=filename;
    filename = "tmp_THistSvc_"+ std::to_string(ishared++)+".root";

    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "Creating temp file \"" << filename
            << "\" and realfilename="<<realfilename << endmsg;
    m_sharedFiles[stream]=realfilename;


    int r = p_fileMgr->open(Io::ROOT,name(), filename, (Io::WRITE|Io::CREATE|Io::EXCL),
                vf,"HIST");

    if (r != 0) {
      m_log << "Unable to open ROOT file " << filename << " for writing"
        << endmsg;
      return StatusCode::FAILURE;
    }

    f = (TFile*)vf;

  } else if (newMode == THistSvc::UPDATE) {
    // update file

    int r = p_fileMgr->open(Io::ROOT,name(), filename, (Io::WRITE|Io::CREATE),
                vf, "HIST");

    if (r != 0) {
      m_log << "Unable to open ROOT file " << filename << " for appending"
        << endmsg;
      return StatusCode::FAILURE;
    }

    f = (TFile*)vf;

  }

  m_files[stream] = make_pair(f,newMode);
  m_fileStreams.insert(make_pair(filename,stream));

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "Opening TFile \"" << filename << "\"  stream: \""
          << stream << "\"  mode: \"" << typ << "\"" << " comp level: " << cl
          << endmsg;

  return StatusCode::SUCCESS;
}

void THistSvc::copyFileLayout ( TDirectory *  dst, TDirectory *  src  )

private

helper function to recursively copy the layout of a TFile into a new TFile

Definition at line 1872 of file THistSvc.cpp.

                                                          {

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG
      << "copyFileLayout() to dst path: " << dst->GetPath () << endmsg;

    // strip out URLs
    TString path ((char*)strstr (dst->GetPath(), ":"));
    path.Remove (0, 2);

    src->cd (path);
    TDirectory *cur_src_dir = gDirectory;

    // loop over all keys in this directory
    TList *key_list = cur_src_dir->GetListOfKeys ();
    int n = key_list->GetEntries ();
    for ( int j = 0; j < n; ++j ) {
      TKey *k = (TKey*)key_list->At (j);
      const std::string src_pathname = cur_src_dir->GetPath()
                                     + std::string("/")
                                     + k->GetName();
      TObject *o=src->Get (src_pathname.c_str());

    if ( o && o->IsA()->InheritsFrom ("TDirectory")) {
      if (m_log.level() <= MSG::VERBOSE)
    m_log << MSG::VERBOSE << " subdir [" << o->GetName() << "]..."
          << endmsg;
        dst->cd ();
        TDirectory * dst_dir = dst->mkdir (o->GetName(), o->GetTitle());
      copyFileLayout (dst_dir, src);
      }
    } // loop over keys
    return;
}

StatusCode THistSvc::deReg ( TObject *  obj )

override

Definition at line 913 of file THistSvc.cpp.

                            {

  auto itr = m_tobjs.find(obj);
  if (itr != m_tobjs.end()) {
    THistID hid = itr->second;

    auto itr2 = m_uids.find(hid.id);
    if (itr2 == m_uids.end()) {
      m_log << MSG::ERROR << "Problems deregistering TObject \""
            << obj->GetName()
            << "\" with id \"" << hid.id << "\"" << endmsg;
      return StatusCode::FAILURE;
    }

    std::string id,root,rem;
    parseString(hid.id, root, rem);

    auto mitr = m_ids.equal_range(rem);
    auto itr3 = std::find_if( mitr.first, mitr.second, [&](idMap::const_reference i) 
                                   { return i.second.obj == obj; } ) ;
    if (itr3 != mitr.second ) {
      m_log << MSG::ERROR << "Problems deregistering TObject \""
            << obj->GetName()
            << "\" with id \"" << hid.id << "\"" << endmsg;
      return StatusCode::FAILURE;
    }

    m_tobjs.erase(itr);
    m_uids.erase(itr2);
    m_ids.erase(itr3);

    return StatusCode::SUCCESS;

  } else {
    m_log << MSG::ERROR << "Cannot unregister TObject \"" << obj->GetName()
          << "\": not known to THistSvc" << endmsg;
    return StatusCode::FAILURE;
  }

}

StatusCode THistSvc::deReg ( const std::string &  name )

override

Definition at line 958 of file THistSvc.cpp.

                                   {

  auto itr = m_uids.find(id);
  if (itr == m_uids.end()) {
    m_log << MSG::ERROR << "Problems deregistering id \""
          << id << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  TObject* obj = itr->second.obj;
  return deReg(obj);
}

std::string THistSvc::dirname ( std::string &  dir ) const

private

Definition at line 1684 of file THistSvc.cpp.

                                      {


  string::size_type i = dir.find("/");

  if (i == string::npos) return {};

  if ( i == 0 ) {
    dir.erase(0,1);
    return dirname(dir);
  }

  string root = dir.substr(0,i);
  dir.erase(0,i);

  return root;

}

bool THistSvc::exists ( const std::string &  name ) const

override

Definition at line 1806 of file THistSvc.cpp.

                                              {

  TH1* h;
  return getHist_i(name,h,true).isSuccess();

}

StatusCode THistSvc::finalize ( )

override

Definition at line 241 of file THistSvc.cpp.

                   {

  GlobalDirectoryRestore restore;

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "THistSvc::finalize" << endmsg;

#ifndef NDEBUG
  if (m_log.level() <= MSG::DEBUG) {
  for (const auto& uid : m_uids ) {

    TObject* to = uid.second.obj;

    string dirname("none");
      if (to && to->IsA()->InheritsFrom("TTree")) {
      TTree* tr = dynamic_cast<TTree*>(to);
      if (tr->GetDirectory() != 0) {
        dirname = tr->GetDirectory()->GetPath();
      }
      } else if (to && to->IsA()->InheritsFrom("TGraph")) {
      if (!uid.second.temp) {
        dirname = uid.second.file->GetPath();
        string id2(uid.second.id);
        id2.erase(0,id2.find("/",1));
        id2.erase(id2.rfind("/"), id2.length());
        if (id2.find("/") == 0) {
          id2.erase(0,1);
        }
        dirname += id2;
      } else {
        dirname = "/tmp";
      }
      } else if (to && to->IsA()->InheritsFrom("TH1")) {
      TH1* th = dynamic_cast<TH1*>(to);
      if (th == 0) {
        m_log << MSG::ERROR << "Couldn't dcast: " << uid.first << endmsg;
      } else {
        if (th->GetDirectory() != 0) {
          dirname = th->GetDirectory()->GetPath();
        }
      }
      } else if (! to ) {
    m_log << MSG::WARNING << uid.first << " has NULL TObject ptr"
          << endmsg;
    }

      m_log << MSG::DEBUG << "uid: \"" << uid.first << "\"  temp: "
            << uid.second.temp << "  dir: " << dirname
            << endmsg;
  }
  }
#endif

  StatusCode sc = write();
  if (sc.isFailure()) {
    m_log << MSG::ERROR << "problems writing histograms" << endmsg;
  }

  if (m_print) {
    m_log << MSG::INFO << "Listing contents of ROOT files: " << endmsg;
  }
  vector<TFile*> deleted_files;
  for (auto& itr : m_files ) {

    if (find(deleted_files.begin(), deleted_files.end(), itr.second.first) ==
        deleted_files.end()) {
      deleted_files.push_back(itr.second.first);

#ifndef NDEBUG
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "finalizing stream/file " << itr.first << ":"
              << itr.second.first->GetName()
              << endmsg;
#endif
    } else {
#ifndef NDEBUG
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "already finalized stream " << itr.first << endmsg;
#endif
      continue;
    }


    if (m_print && m_log.level() <= MSG::INFO) {

      m_log << MSG::INFO;
      m_log << "==> File: " << itr.second.first->GetName()
            << "  stream: " << itr.first << endmsg;

      itr.second.first->Print("base");
    }

    string tmpfn=itr.second.first->GetName();

    p_fileMgr->close(itr.second.first, name());

    IIncidentSvc *pi = nullptr;
    if (service("IncidentSvc",pi).isFailure()) {
      m_log << MSG::ERROR << "Unable to get the IncidentSvc" << endmsg;
      return StatusCode::FAILURE;
    }

    if (itr.second.second==SHARE) {

      //Merge File
      void* vf = nullptr;
      int r = p_fileMgr->open(Io::ROOT,name(), m_sharedFiles[itr.first],
                  Io::WRITE|Io::APPEND,vf,"HIST");

      if (r) {
    m_log << MSG::ERROR << "unable to open Final Output File: \""
          << m_sharedFiles[itr.first] << "\" for merging"
          << endmsg;
        return StatusCode::FAILURE;
      }

      TFile *outputfile = (TFile*) vf;
      pi->fireIncident(FileIncident(name(), IncidentType::WroteToOutputFile,
                                     m_sharedFiles[itr.first]));

      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "THistSvc::write()::Merging Rootfile "<<endmsg;

      vf = nullptr;
      r = p_fileMgr->open(Io::ROOT,name(),tmpfn,Io::READ,vf,"HIST");

      if (r) {
    m_log << MSG::ERROR << "unable to open temporary file: \""
          << tmpfn << endmsg;
        return StatusCode::FAILURE;
      }

      TFile *inputfile = (TFile*) vf;

      outputfile->SetCompressionLevel( inputfile->GetCompressionLevel() );

      MergeRootFile(outputfile, inputfile);

      outputfile->Write();
      p_fileMgr->close(outputfile,name());
      p_fileMgr->close(inputfile,name());

      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "Trying to remove temporary file \"" << tmpfn
              << "\""<<endmsg;

      std::remove(tmpfn.c_str());
    }
    delete itr.second.first;
  }

  m_sharedFiles.clear();
  m_fileStreams.clear();
  m_files.clear();
  m_uids.clear();
  m_ids.clear();
  m_tobjs.clear();

  return Service::finalize();
}

bool THistSvc::findStream ( const std::string &  name, std::string &  root, std::string &  rem, TFile *&  file  ) const

private

Definition at line 1161 of file THistSvc.cpp.

                                       {

  auto pos = id.find("/");

  if (pos == string::npos) {
    stream = "temp";
    rem = id;
  } else if (pos != 0) {
    stream = "temp";
    rem = id;
  } else {

    auto pos2 = id.find("/",pos+1);

    if (pos2 == string::npos) {
      m_log << MSG::ERROR << "badly formed Hist/Tree id: \"" << id << "\""
            << endmsg;
      return false;
    }

    parseString(id,stream,rem);

  }

  if (stream == "temp") {
    file = nullptr;
    return true;
  }

  auto itr = m_files.find(stream);
  file = (itr != m_files.end() ? itr->second.first : nullptr );
  if (!file) {
    m_log << MSG::WARNING << "no stream \"" << stream
          << "\" associated with id: \"" << id << "\""
          << endmsg;
  }

  return true;
}

StatusCode THistSvc::getGraph ( const std::string &  name, TGraph *&  hist  ) const

override

Definition at line 1111 of file THistSvc.cpp.

                                                           {
  return getHist_i(id, hist);
}

std::vector< std::string > THistSvc::getGraphs ( ) const

override

Definition at line 1118 of file THistSvc.cpp.

                          {

  std::vector<std::string> names;
  names.reserve(m_uids.size());
  transform_if( std::begin(m_uids), std::end(m_uids),
                std::back_inserter(names), select1st,
                [](uidMap::const_reference i) { 
                    return i.second.obj->IsA()->InheritsFrom("TTree"); }
  );
  return names;

}

StatusCode THistSvc::getHist ( const std::string &  name, TH1 *&  hist  ) const

override

Definition at line 1052 of file THistSvc.cpp.

                                                       {
  return getHist_i(id, hist);
}

StatusCode THistSvc::getHist ( const std::string &  name, TH2 *&  hist  ) const

override

Definition at line 1059 of file THistSvc.cpp.

                                                       {
  return getHist_i(id, hist);
}

StatusCode THistSvc::getHist ( const std::string &  name, TH3 *&  hist  ) const

override

Definition at line 1066 of file THistSvc.cpp.

                                                       {
  return getHist_i(id, hist);
}

template<typename T >

StatusCode THistSvc::getHist_i ( const std::string &  name, T *&  hist, bool  quiet = false  ) const

private

Definition at line 142 of file THistSvc.icc.

                                                                   {
  // id starts with "/": unique

  GlobalDirectoryRestore restore;

  std::string idr(id);
  removeDoubleSlash( idr );

  if (idr.find("/") == 0) {
    auto itr = m_uids.find(id);
    if (itr == m_uids.end()) {
      if (!quiet) {
    m_log << MSG::ERROR << "Could not locate Hist with id \"" << idr << "\""
          << endmsg;
      }
      hist = nullptr;
      return StatusCode::FAILURE;
    }

    THistID hid = itr->second;
    if (!quiet) {
      hist = dynamic_cast<T*>(hid.obj);
      if ( !hist ) {
    m_log << MSG::ERROR << "dcast failed, Hist id: \"" << idr << "\"" 
          << endmsg;
    return StatusCode::FAILURE;
      }
      if (m_log.level() <= MSG::VERBOSE) {
    m_log << MSG::VERBOSE << "found unique Hist title: \"" 
          << hist->GetTitle()
          << "\"  id: \"" << idr << "\"" << endmsg;
      }
    } else {
      if (m_log.level() <= MSG::VERBOSE) {
    m_log << MSG::VERBOSE << "found unique Hist id: \"" << idr 
          << "\" type: \"" << hid.obj->IsA()->GetName() << "\""
          << endmsg;
      }
    }

    return StatusCode::SUCCESS;


    // not necessarily unique
  } else {

      auto mitr = m_ids.equal_range(idr);


    if (mitr.first == mitr.second) {
      m_log << MSG::ERROR << "Could not locate Hist with id \"" << idr << "\""
        << endmsg;
      hist = nullptr;
      return StatusCode::FAILURE;
    } else {

      if (distance(mitr.first,mitr.second) == 1) {
        THistID hid = mitr.first->second;
    if (!quiet) {
      hist = dynamic_cast<T*>(hid.obj);
      if (hist == 0) {
        m_log << MSG::ERROR << "dcast failed" << endmsg;
        return StatusCode::FAILURE;
      }
      if (m_log.level() <= MSG::VERBOSE) {
        m_log << MSG::VERBOSE << "found Hist title: \"" << hist->GetTitle()
          << "\"  id: \"" << idr << "\"" << endmsg;
      }
    } else {
      if (m_log.level() <= MSG::VERBOSE) {
        m_log << MSG::VERBOSE << "found Hist id: \"" << idr << "\" type: \""
          << hid.obj->IsA()->GetName() << "\""
          << endmsg;
      }
    }
        return StatusCode::SUCCESS;
      } else {
    if (!quiet) {
      // return failure if trying to GET a single hist
      m_log << MSG::ERROR << "Multiple matches with id \"" << idr << "\"."
        << " Further specifications required."
        << endmsg;
      hist = nullptr;
      return StatusCode::FAILURE;
    } else {
      // return a SUCCESS if just INQUIRING
      m_log << MSG::INFO << "Found multiple matches with id \"" << idr 
          << "\"" << endmsg;
      hist = nullptr;
      return StatusCode::SUCCESS;
    }
      }
    }
  }
}

std::vector< std::string > THistSvc::getHists ( ) const

override

Definition at line 1073 of file THistSvc.cpp.

                         {

  std::vector<std::string> names;
  names.reserve(m_uids.size());
  transform_if( std::begin(m_uids), std::end(m_uids),
                std::back_inserter(names), select1st,
                [](uidMap::const_reference i) { 
                    return i.second.obj->IsA()->InheritsFrom("TH11"); }
  );
  return names;

}

StatusCode THistSvc::getTHists ( TDirectory *  td, TList &  tl, bool  recurse = false  ) const

override

Definition at line 451 of file THistSvc.cpp.

                                                              {
  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  if (!td->cd()) {
    m_log << MSG::ERROR << "getTHists: No such TDirectory \"" << td->GetPath()
          << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "getTHists: \"" << td->GetPath() << "\": found "
          << td->GetListOfKeys()->GetSize() << " keys" << endmsg;

  TIter nextkey(td->GetListOfKeys());
  while (TKey *key = (TKey*)nextkey()) {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "  key: " << key->GetName();
    TObject *obj = key->ReadObj();
    if (obj != 0 && obj->IsA()->InheritsFrom("TDirectory")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
    } else if (obj != 0 && obj->IsA()->InheritsFrom("TH1")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
      tl.Add(obj);
    } else if (obj != 0) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " [" << obj->IsA()->GetName() << "]";
    }
    if (m_log.level() <= MSG::DEBUG)
      m_log << endmsg;
  }

  // operate recursively
  if (rcs) {
    nextkey = td->GetListOfKeys();
    while (TKey *key = (TKey*)nextkey()) {
      TObject *obj = key->ReadObj();
      if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
          TDirectory *tt = dynamic_cast<TDirectory*>(obj);
          getTHists(tt, tl, rcs);
      }
    }
  }

  return StatusCode::SUCCESS;


}

StatusCode THistSvc::getTHists ( const std::string &  name, TList &  tl, bool  recurse = false  ) const

override

Definition at line 506 of file THistSvc.cpp.

                                                                    {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

  std::string stream,rem,r2;
  parseString(dir,stream,rem);

  auto itr = m_files.find(stream);
  if (itr != m_files.end()) {
    r2 = itr->second.first->GetName();
    r2 += ":/";
    r2 += rem;

    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTHists: \"" << dir
            << "\" looks like a stream name."  << " associated TFile: \""
            << itr->second.first->GetName() << "\"" << endmsg;

    if (gDirectory->cd(r2.c_str())) {
      m_curstream = stream;
      sc = getTHists(gDirectory,tl,rcs);
      m_curstream = "";
      return sc;
    } else {
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "getTHists: no such TDirectory \""
              << r2 << "\"" << endmsg;
    }

  } else {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTHists: stream \"" << stream << "\" not found"
            << endmsg;
  }

  if (!gDirectory->cd(dir.c_str())) {
    m_log << MSG::ERROR << "getTHists: No such TDirectory/stream \"" << dir
          << "\"" << endmsg;
    sc = StatusCode::FAILURE;
  } else {
    sc = getTHists(gDirectory,tl,rcs);
  }

  return sc;

}

StatusCode THistSvc::getTHists ( TDirectory *  td, TList &  tl, bool  recurse = false, bool  reg = false  )

override

Definition at line 659 of file THistSvc.cpp.

                                                                  {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  if (!td->cd()) {
    m_log << MSG::ERROR << "getTHists: No such TDirectory \"" << td->GetPath()
          << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "getTHists: \"" << td->GetPath() << "\": found "
          << td->GetListOfKeys()->GetSize() << " keys" << endmsg;

  TIter nextkey(td->GetListOfKeys());
  while (TKey *key = (TKey*)nextkey()) {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "  key: " << key->GetName();
    TObject *obj = key->ReadObj();
    if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
    } else if (obj && obj->IsA()->InheritsFrom("TH1")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
      tl.Add(obj);
      if (reg && m_curstream != "") {
        string dir = td->GetPath();
        string fil = td->GetFile()->GetName();
        dir.erase(0,fil.length()+1);
        string id = "/" + m_curstream;
        if ( dir == "/" ) {
          id = id + "/" + key->GetName();
        } else {
          id = id + dir + "/" + key->GetName();
        }
        if (!exists(id)) {
          if (m_log.level() <= MSG::DEBUG)
            m_log << "  reg as \"" << id << "\"";
          regHist(id).ignore();
        } else {
          if (m_log.level() <= MSG::DEBUG)
            m_log << "  already registered";
        }
      }
    } else if (obj) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " [" << obj->IsA()->GetName() << "]";
    }
    if (m_log.level() <= MSG::DEBUG)
      m_log << endmsg;
  }

  // operate recursively
  if (rcs) {
    nextkey = td->GetListOfKeys();
    while (TKey *key = (TKey*)nextkey()) {
      TObject *obj = key->ReadObj();
      if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
          TDirectory *tt = dynamic_cast<TDirectory*>(obj);
          getTHists(tt, tl, rcs, reg);
      }
    }
  }

  return StatusCode::SUCCESS;

}

StatusCode THistSvc::getTHists ( const std::string &  name, TList &  tl, bool  recurse = false, bool  reg = false  )

override

Definition at line 733 of file THistSvc.cpp.

                                                                        {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

  std::string stream,rem,r2;
  parseString(dir,stream,rem);

  auto itr = m_files.find(stream);
  if (itr != m_files.end()) {
    r2 = itr->second.first->GetName();
    r2 += ":/";
    r2 += rem;

    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTHists: \"" << dir
            << "\" looks like a stream name."  << " associated TFile: \""
            << itr->second.first->GetName() << "\"" << endmsg;

    if (gDirectory->cd(r2.c_str())) {
      m_curstream = stream;
      sc = getTHists(gDirectory,tl,rcs,reg);
      m_curstream.clear();
      return sc;
    } 
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTHists: no such TDirectory \""
            << r2 << "\"" << endmsg;

  } else {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTHists: stream \"" << stream << "\" not found"
            << endmsg;
  }

  if (!gDirectory->cd(dir.c_str())) {
    m_log << MSG::ERROR << "getTHists: No such TDirectory/stream \"" << dir
          << "\"" << endmsg;
    sc = StatusCode::FAILURE;
  } else {
    if (reg) {
      m_log << MSG::WARNING << "Unable to register histograms automatically "
            << "without a valid stream name" << endmsg;
      reg = false;
    }
    sc = getTHists(gDirectory,tl,rcs,reg);
  }

  return sc;

}

StatusCode THistSvc::getTree ( const std::string &  name, TTree *&  hist  ) const

override

Definition at line 1088 of file THistSvc.cpp.

                                                         {
  return getHist_i(id, hist);
}

std::vector< std::string > THistSvc::getTrees ( ) const

override

Definition at line 1095 of file THistSvc.cpp.

                         {

  std::vector<std::string> names;
  names.reserve(m_uids.size());
  transform_if( std::begin(m_uids), std::end(m_uids),
                std::back_inserter(names),
                select1st,
                [](uidMap::const_reference i) { 
                    return i.second.obj->IsA()->InheritsFrom("TTree"); }
  );
  return names;

}

StatusCode THistSvc::getTTrees ( TDirectory *  td, TList &  tl, bool  recurse = false  ) const

override

Definition at line 559 of file THistSvc.cpp.

                                                              {
  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  if (!td->cd()) {
    m_log << MSG::ERROR << "getTTrees: No such TDirectory \""
          << td->GetPath() << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "getTHists: \"" << td->GetPath() << "\": found "
          << td->GetListOfKeys()->GetSize() << " keys" << endmsg;

  TIter nextkey(td->GetListOfKeys());
  while (TKey *key = (TKey*)nextkey()) {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "  key: " << key->GetName();
    TObject *obj = key->ReadObj();
    if (obj != 0 && obj->IsA()->InheritsFrom("TDirectory")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
    } else if (obj != 0 && obj->IsA()->InheritsFrom("TTree")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
      tl.Add(obj);
    } else if (obj != 0) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " [" << obj->IsA()->GetName() << "]";
    }
    m_log << endmsg;
  }

  // operate recursively
  if (rcs) {
    nextkey = td->GetListOfKeys();
    while (TKey *key = (TKey*)nextkey()) {
      TObject *obj = key->ReadObj();
      if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
          TDirectory *tt = dynamic_cast<TDirectory*>(obj);
          getTTrees(tt, tl, rcs);
      }
    }
  }

  return StatusCode::SUCCESS;

}

StatusCode THistSvc::getTTrees ( const std::string &  name, TList &  tl, bool  recurse = false  ) const

override

Definition at line 612 of file THistSvc.cpp.

                                                                    {
  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

  std::string stream,rem,r2;
  parseString(dir,stream,rem);

  auto itr = m_files.find(stream);
  if (itr != m_files.end()) {
    r2 = itr->second.first->GetName();
    r2 += ":/";
    r2 += rem;

    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTTrees: \"" << dir
            << "\" looks like a stream name."  << " associated TFile: \""
            << itr->second.first->GetName() << "\"" << endmsg;

    if (gDirectory->cd(r2.c_str())) {
      return getTTrees(gDirectory,tl,rcs);
    }
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTTrees: no such TDirectory \""
            << r2 << "\"" << endmsg;
  } else {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTTrees: stream \"" << stream << "\" not found"
            << endmsg;
  }

  if (!gDirectory->cd(dir.c_str())) {
    m_log << MSG::ERROR << "getTTrees: No such TDirectory/stream \"" << dir
          << "\"" << endmsg;
    sc = StatusCode::FAILURE;
  } else {
    sc = getTTrees(gDirectory,tl,rcs);
  }
  return sc;
}

StatusCode THistSvc::getTTrees ( TDirectory *  td, TList &  tl, bool  recurse = false, bool  reg = false  )

override

Definition at line 790 of file THistSvc.cpp.

                                                                  {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  if (!td->cd()) {
    m_log << MSG::ERROR << "getTTrees: No such TDirectory \""
          << td->GetPath() << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "getTHists: \"" << td->GetPath() << "\": found "
          << td->GetListOfKeys()->GetSize() << " keys" << endmsg;

  TIter nextkey(td->GetListOfKeys());
  while (TKey *key = (TKey*)nextkey()) {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "  key: " << key->GetName();
    TObject *obj = key->ReadObj();
    if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
    } else if (obj && obj->IsA()->InheritsFrom("TTree")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " (" << obj->IsA()->GetName() << ")";
      tl.Add(obj);
      if (reg && m_curstream != "") {
        string dir = td->GetPath();
        string fil = td->GetFile()->GetName();
        dir.erase(0,fil.length()+1);
        string id = "/" + m_curstream;
        if ( dir == "/" ) {
          id = id + "/" + key->GetName();
        } else {
          id = id + dir + "/" + key->GetName();
        }
        if (!exists(id)) {
          if (m_log.level() <= MSG::DEBUG)
            m_log << "  reg as \"" << id << "\"";
          regHist(id).ignore();
        } else {
          if (m_log.level() <= MSG::DEBUG)
            m_log << "  already registered";
        }
      }
    } else if (obj != 0) {
      if (m_log.level() <= MSG::DEBUG)
        m_log << " [" << obj->IsA()->GetName() << "]";
    }
    if (m_log.level() <= MSG::DEBUG)
      m_log << endmsg;
  }

  // operate recursively
  if (rcs) {
    nextkey = td->GetListOfKeys();
    while (TKey *key = (TKey*)nextkey()) {
      TObject *obj = key->ReadObj();
      if (obj && obj->IsA()->InheritsFrom("TDirectory")) {
          TDirectory *tt = dynamic_cast<TDirectory*>(obj);
          getTTrees(tt, tl, rcs, reg);
      }
    }
  }

  return StatusCode::SUCCESS;

}

StatusCode THistSvc::getTTrees ( const std::string &  name, TList &  tl, bool  recurse = false, bool  reg = false  )

override

Definition at line 864 of file THistSvc.cpp.

                                                                        {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

  std::string stream,rem,r2;
  parseString(dir,stream,rem);

  auto itr = m_files.find(stream);
  if (itr != m_files.end()) {
    r2 = itr->second.first->GetName();
    r2 += ":/";
    r2 += rem;

    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTTrees: \"" << dir
            << "\" looks like a stream name."  << " associated TFile: \""
            << itr->second.first->GetName() << "\"" << endmsg;

    if (gDirectory->cd(r2.c_str())) {
      return getTTrees(gDirectory,tl,rcs,reg);
    } else {
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "getTTrees: no such TDirectory \""
              << r2 << "\"" << endmsg;
    }

  } else {
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "getTTrees: stream \"" << stream << "\" not found"
            << endmsg;
  }

  if (!gDirectory->cd(dir.c_str())) {
    m_log << MSG::ERROR << "getTTrees: No such TDirectory/stream \"" << dir
          << "\"" << endmsg;
    return StatusCode::FAILURE;
  } 

  return getTTrees(gDirectory,tl,rcs,reg);

}

void THistSvc::handle ( const Incident &  )

override

Definition at line 1816 of file THistSvc.cpp.

                                    {

  if (signaledStop) return ;

  if (m_maxFileSize.value() == -1) return;

  // convert to bytes.
  Long64_t mfs = (Long64_t)m_maxFileSize.value() * (Long64_t)1048576;
  Long64_t mfs_warn = mfs * 95 / 100;

  updateFiles();

  map<string, pair<TFile*,Mode> >::const_iterator itr;
  for (const auto& f : m_files) {
    TFile* tf = f.second.first;

#ifndef NDEBUG
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "stream: " << f.first << "  name: "
            << tf->GetName() << "  size: " << tf->GetSize()
            << endmsg;
#endif

    // Signal job to terminate if output file is too large
    if (tf->GetSize() > mfs) {

      signaledStop = true;

      m_log << MSG::FATAL << "file \"" << tf->GetName()
            << "\" associated with stream \"" << f.first
            << "\" has exceeded the max file size of "
            << m_maxFileSize.value() << "MB. Terminating Job."
            << endmsg;

      IEventProcessor* evt = nullptr;
      if (service("ApplicationMgr", evt, true).isSuccess()) {
        evt->stopRun();
        evt->release();
      } else {
        abort();
      }
    } else if (tf->GetSize() > mfs_warn) {
      m_log << MSG::WARNING << "file \"" << tf->GetName()
            << "\" associated with stream \"" << f.first
            << "\" is at 95% of its maximum allowable file size of "
            << m_maxFileSize.value() << "MB"
            << endmsg;
    }
  }
}

StatusCode THistSvc::initialize ( )

override

Definition at line 78 of file THistSvc.cpp.

                     {
  GlobalDirectoryRestore restore;

  // Super ugly hack to make sure we have the OutputLevel set first, so we
  // can see DEBUG printouts in update handlers.
  auto jos = serviceLocator()->service<IJobOptionsSvc>( "JobOptionsSvc", true );
  if( jos ) { 
    const auto props = jos->getProperties( name() );

    if ( props ) {
      auto i = std::find_if(std::begin(*props),std::end(*props),
                            [](const Property* p) { return p->name() == "OutputLevel"; } );
      if ( i != std::end(*props)) {
          setProperty( **i ).ignore();
          m_log.setLevel( m_outputLevel.value() );
      }
    }
  }


  StatusCode status = Service::initialize();
  m_log.setLevel( m_outputLevel.value() );

  if (status.isFailure()) {
    m_log << MSG::ERROR << "initializing service" << endmsg;
    return status;
  }

  StatusCode st(StatusCode::SUCCESS);

  try {
    setupOutputFile( m_outputfile );
  } catch ( GaudiException& err ) {
    m_log << MSG::ERROR
          << "Caught: " << err << endmsg;
    st = StatusCode::FAILURE;
  }

  try {
    setupInputFile( m_inputfile );
  } catch ( GaudiException& err ) {
    m_log << MSG::ERROR
          << "Caught: " << err << endmsg;
    st = StatusCode::FAILURE;
  }

  // Protect against multiple instances of TROOT
  if ( !gROOT )   {
    static TROOT root("root","ROOT I/O");
    //    gDebug = 99;
  } else {
    if (m_log.level() <= MSG::VERBOSE)
      m_log << MSG::VERBOSE << "ROOT already initialized, debug = "
            << gDebug<< endmsg;
  }

  if (service("IncidentSvc", p_incSvc, true).isFailure()) {
    m_log << MSG::ERROR << "unable to get the IncidentSvc" << endmsg;
    st = StatusCode::FAILURE;
  } else {
    p_incSvc->addListener( this, "EndEvent", 100, true);
  }

  if (service("FileMgr",p_fileMgr,true).isFailure()) {
    m_log << MSG::ERROR << "unable to get the FileMgr" << endmsg;
    st = StatusCode::FAILURE;
  } else {
    m_log << MSG::DEBUG << "got the FileMgr" << endmsg;
  }


  // Register open/close callback actions

  using namespace std::placeholders;
  auto boa = [this](const Io::FileAttr* fa, const std::string& caller) { return this->rootOpenAction(fa,caller); };
  if (p_fileMgr->regAction(boa, Io::OPEN, Io::ROOT).isFailure()) {
    m_log << MSG::ERROR
      << "unable to register ROOT file open action with FileMgr"
      << endmsg;
  }
  auto bea = [this](const Io::FileAttr* fa, const std::string& caller) { return this->rootOpenErrAction(fa,caller); };
  if (p_fileMgr->regAction(bea, Io::OPEN_ERR, Io::ROOT).isFailure()) {
    m_log << MSG::ERROR
      << "unable to register ROOT file open Error action with FileMgr"
      << endmsg;
  }


  m_okToConnect = true;

  if (m_delayConnect) {
    if (!m_inputfile.value().empty()) { setupInputFile(m_inputfile); }
    if (!m_outputfile.value().empty()) { setupOutputFile(m_outputfile); }

    m_delayConnect = false;

  }
  m_alreadyConnectedOutFiles.clear();
  m_alreadyConnectedInFiles.clear();


  IIoComponentMgr* iomgr = nullptr;

  if (service("IoComponentMgr", iomgr, true).isFailure()) {
    m_log << MSG::ERROR << "unable to get the IoComponentMgr" << endmsg;
    st = StatusCode::FAILURE;
  } else {

    if ( !iomgr->io_register (this).isSuccess() ) {
      m_log << MSG::ERROR
            << "could not register with the I/O component manager !"
            << endmsg;
      st = StatusCode::FAILURE;
    } else {
      bool all_good = true;
      // register input/output files...
      for ( const auto& reg : m_files ) { 
        const std::string& fname = reg.second.first->GetName();
        const IIoComponentMgr::IoMode::Type iomode =
          ( reg.second.second==THistSvc::READ
        ? IIoComponentMgr::IoMode::READ
        : IIoComponentMgr::IoMode::WRITE );
        if ( !iomgr->io_register (this, iomode, fname).isSuccess () ) {
          m_log << MSG::WARNING << "could not register file ["
                << fname << "] with the I/O component manager..." << endmsg;
          all_good = false;
        } else {
          m_log << MSG::INFO << "registered file [" << fname << "]... [ok]"
              << endmsg;
        }
      }
      if (!all_good) {
        m_log << MSG::ERROR
              << "problem while registering input/output files with "
              << "the I/O component manager !" << endmsg;
        st = StatusCode::FAILURE;
      }
    }

  }

  if (st.isFailure()) {
    m_log << MSG::FATAL << "Unable to initialize THistSvc" << endmsg;
  }

  return st;

}

StatusCode THistSvc::io_reinit ( )

override

callback method to reinitialize the internal state of the component for I/O purposes (e.g.

upon fork(2))

Definition at line 1912 of file THistSvc.cpp.

{
  bool all_good = true;
  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "reinitializing I/O..." << endmsg;

  // retrieve the I/O component manager...

  IIoComponentMgr* iomgr = nullptr;

  if (service("IoComponentMgr", iomgr, true).isFailure()) {
    m_log << MSG::ERROR << "could not retrieve I/O component manager !"
          << endmsg;
    return StatusCode::FAILURE;
  }

  GlobalDirectoryRestore restore;
  // to hide the expected errors upon closing the files whose
  // file descriptors have been swept under the rug...
  gErrorIgnoreLevel = kFatal;

  typedef std::map<std::string, std::pair<TFile*,Mode> > FileReg_t;

  for (auto & ifile : m_files ) {
    TFile *f = ifile.second.first;
    std::string fname = f->GetName();
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "file [" << fname << "] mode: ["
          << f->GetOption() << "] r:"
          << f->GetFileBytesRead()
          << " w:" << f->GetFileBytesWritten()
          << " cnt:" << f->GetFileCounter()
          << endmsg;

    if ( ifile.second.second == READ ) {
      if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG
            << "  TFile opened in READ mode: not reassigning names" << endmsg;
      continue;
    }

    if ( !iomgr->io_retrieve (this, fname).isSuccess () ) {
      m_log << MSG::ERROR << "could not retrieve new name for [" << fname
            << "] !!" << endmsg;
      all_good = false;
      continue;
    } else {
      if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "got a new name [" << fname << "]..." << endmsg;
    }
    // create a new TFile
    // TFile *newfile = TFile::Open (fname.c_str(), f->GetOption());

    void* vf;
    Option_t *opts = f->GetOption();
    int r = p_fileMgr->open(Io::ROOT,name(),fname,Io::WRITE,vf,"HIST");
    if (r != 0) {
      m_log << MSG::ERROR << "unable to open file \"" << fname
        << "\" for writing" << endmsg;
      return StatusCode::FAILURE;
    }
    TFile *newfile = (TFile*) vf;
    newfile->SetOption(opts);


    if (ifile.second.second != THistSvc::READ) {
      copyFileLayout (newfile, f);
      ifile.second.first = newfile;
    }

    // loop over all uids and migrate them to the new file
    // XXX FIXME: this double loop sucks...
    for ( auto& uid : m_uids ) {
      THistID& hid = uid.second;
      if ( hid.file != f ) continue;
      TDirectory *olddir = this->changeDir (hid);
      hid.file = newfile;
      // side-effect: create needed directories...
      TDirectory *newdir = this->changeDir (hid);
      TClass *cl = hid.obj->IsA();

      // migrate the objects to the new file.
      // thanks to the object model of ROOT, it is super easy.
      if (cl->InheritsFrom ("TTree")) {
        dynamic_cast<TTree*> (hid.obj)->SetDirectory (newdir);
        dynamic_cast<TTree*> (hid.obj)->Reset();
      }
      else if (cl->InheritsFrom ("TH1")) {
        dynamic_cast<TH1*> (hid.obj)->SetDirectory (newdir);
        dynamic_cast<TH1*> (hid.obj)->Reset();
      }
      else if (cl->InheritsFrom ("TGraph")) {
        olddir->Remove (hid.obj);
        newdir->Append (hid.obj);
      } else {
        m_log << MSG::ERROR
              << "id: \"" << hid.id << "\" is not a inheriting from a class "
              << "we know how to handle (received [" << cl->GetName()
              << "], " << "expected [TTree, TH1 or TGraph]) !"
              << endmsg
              << "attaching to current dir [" << newdir->GetPath() << "] "
              << "nonetheless..." << endmsg;
        olddir->Remove (hid.obj);
        newdir->Append (hid.obj);
      }
    }
    f->ReOpen ("READ");
    p_fileMgr->close(f,name());
    f = newfile;
  }

  return all_good ? StatusCode::SUCCESS : StatusCode::FAILURE;
}

void THistSvc::MergeRootFile ( TDirectory *  target, TDirectory *  source  )

private

Definition at line 1730 of file THistSvc.cpp.

                                                                   {

  if (m_log.level() <= MSG::DEBUG)
    m_log <<MSG::DEBUG << "Target path: " << target->GetPath() << endmsg;
  TString path( (char*)strstr(target->GetPath(), ":") );
  path.Remove( 0, 2);

  source->cd(path);
  TDirectory *current_sourcedir = gDirectory;

  // loop over all keys in this directory
  TList *lkeys=current_sourcedir->GetListOfKeys();
  int nkeys=lkeys->GetEntries();
  TKey *key;
  for (int jj=0; jj<nkeys; jj++) {
    key=(TKey*) lkeys->At(jj);
    string pathnameinsource=current_sourcedir->GetPath()+string("/")+key->GetName();
    if (m_log.level() <= MSG::DEBUG)
      m_log <<MSG::DEBUG << "Reading Key:" << pathnameinsource << endmsg;
    //key->Dump();
    //TObject *obj=key->ReadObj();
    TObject *obj=source->Get(pathnameinsource.c_str());

    if (obj) {
    if (obj->IsA()->InheritsFrom("TDirectory") ) {
      // it's a subdirectory

      if (m_log.level() <= MSG::DEBUG)
        m_log <<MSG::DEBUG << "Found subdirectory " << obj->GetName()
              << endmsg;

      // create a new subdir of same name and title in the target file
      target->cd();
      TDirectory *newtargetdir =
        target->mkdir(obj->GetName(), obj->GetTitle() );

      MergeRootFile(newtargetdir, source);

    } else if (obj->IsA()->InheritsFrom("TTree")) {
      if (m_log.level() <= MSG::DEBUG)
        m_log <<MSG::DEBUG << "Found TTree " << obj->GetName() << endmsg;
      TTree *mytree=dynamic_cast<TTree*>(obj);
      int nentries=(int) mytree->GetEntries();
      mytree->SetBranchStatus("*",1);

      if (m_log.level() <= MSG::DEBUG)
        m_log <<MSG::DEBUG << "Dumping TTree " << nentries <<" entries"
              << endmsg;
      //mytree->Print();
      //for (int ij=0; ij<nentries; ij++) {
      //m_log <<MSG::DEBUG << "Dumping TTree Show( " << ij <<" )"
      //<< endmsg;
      //mytree->Show(ij);
      //}
      target->cd();
      mytree->CloneTree();

      //m_log <<MSG::DEBUG << "Writing TTree to target file: ( "
      //<< mycopiedtree->Write(key->GetName()) <<" ) bytes written"
      //<< endmsg;

    } else {
      target->cd();
      obj->Write(key->GetName() );
    }
    }

  } // while ( ( TKey *key = (TKey*)nextkey() ) )

  // save modifications to target file

}

void THistSvc::parseString ( const std::string &  id, std::string &  root, std::string &  rem  ) const

private

Definition at line 1205 of file THistSvc.cpp.

                                                                       {
  auto pos = id.find("/");

  if (pos == string::npos) {
    root.clear();
    rem = id;
  } else if (pos == 0) {
    parseString(id.substr(1),root,rem);
  } else {
    root = id.substr(0,pos);
    rem = id.substr(pos+1);
  }

}

StatusCode THistSvc::readHist ( const std::string &  name, TH1 *&  hist  ) const

private

Definition at line 1133 of file THistSvc.cpp.

                                                        {
  return readHist_i(id, hist);
}

StatusCode THistSvc::readHist ( const std::string &  name, TH2 *&  hist  ) const

private

Definition at line 1140 of file THistSvc.cpp.

                                                        {
  return readHist_i(id, hist);
}

StatusCode THistSvc::readHist ( const std::string &  name, TH3 *&  hist  ) const

private

Definition at line 1147 of file THistSvc.cpp.

                                                        {
  return readHist_i(id, hist);
}

template<typename T >

StatusCode THistSvc::readHist_i ( const std::string &  name, T *&  hist  ) const

private

Definition at line 242 of file THistSvc.icc.

                                                        {

  GlobalDirectoryRestore restore;

  std::string idr(id);
  removeDoubleSlash( idr );

  std::string stream, rem, dir, fdir, bdir, fdir2;
  TFile *file;

  if (!findStream(idr, stream, rem, file) ) {
    return StatusCode::FAILURE;
  }

  if ( !file ) {
    m_log << MSG::ERROR << "no associated file found" << endmsg;
    return StatusCode::FAILURE;
  }

  file->cd("/");

  fdir = idr;
  bdir = dirname(fdir);
  fdir2 = fdir;
  while ( (dir=dirname(fdir)) != "" ) {
    if (! gDirectory->GetKey(dir.c_str())) {
      m_log << MSG::ERROR << "Directory \"" << fdir2 << "\" doesnt exist in "
        << file->GetName() << endmsg;
      return StatusCode::FAILURE;
    }
    gDirectory->cd(dir.c_str());
  }

  TObject *to=nullptr;
  gDirectory->GetObject(fdir.c_str(), to);

  if ( !to ) {
    m_log << MSG::ERROR << "Could not get obj \"" << fdir << "\" in "
      << gDirectory->GetPath() << endmsg;
    return StatusCode::FAILURE;
  }



  hist = dynamic_cast<T*>(to);
  if ( !hist ) {
    m_log << MSG::ERROR << "Could not convert \"" << idr << "\" to a "
      << System::typeinfoName(typeid(*hist)) << " as is a "
      << to->IsA()->GetName()
      << endmsg;
    return StatusCode::FAILURE;
  }


  if (m_log.level() <= MSG::DEBUG) {
    m_log << MSG::DEBUG << "Read in " << hist->IsA()->GetName() << "  \""
      << hist->GetName() << "\" from file "
      << file->GetName() << endmsg;
    hist->Print();
  }

  return StatusCode::SUCCESS;

}

StatusCode THistSvc::readTree ( const std::string &  name, TTree *&  hist  ) const

private

Definition at line 1154 of file THistSvc.cpp.

                                                          {
  return readHist_i(id, hist);
}

StatusCode THistSvc::regGraph ( const std::string &  name )

override

Definition at line 1024 of file THistSvc.cpp.

                                      {
  TGraph *hist = nullptr;
  return regHist_i(hist, id);
}

StatusCode THistSvc::regGraph ( const std::string &  name, TGraph *  hist  )

override

Definition at line 1032 of file THistSvc.cpp.

                                                    {
  if ( strcmp(hist->GetName(),"Graph") == 0 ) {

    std::string id2(id);
    string::size_type i = id2.rfind("/");
    if (i != string::npos) {
      id2.erase(0,i+1);
    }

    m_log << MSG::INFO << "setting name of TGraph id: \"" << id << "\" to \""
          << id2 << "\" since it is unset" << endmsg;
    hist->SetName(id2.c_str());
  }

  return regHist_i(hist, id);
}

StatusCode THistSvc::regHist ( const std::string &  name )

override

Definition at line 974 of file THistSvc.cpp.

                                     {

  TH1 *hist = nullptr;
  return regHist_i(hist, id);
}

StatusCode THistSvc::regHist ( const std::string &  name, TH1 *  hist  )

override

Definition at line 983 of file THistSvc.cpp.

                                                {
  return regHist_i(hist, id);
}

StatusCode THistSvc::regHist ( const std::string &  name, TH2 *  hist  )

override

Definition at line 990 of file THistSvc.cpp.

                                                {
  return regHist_i(hist, id);
}

StatusCode THistSvc::regHist ( const std::string &  name, TH3 *  hist  )

override

Definition at line 997 of file THistSvc.cpp.

                                                {
  return regHist_i(hist, id);
}

template<typename T >

StatusCode THistSvc::regHist_i ( T *  hist, const std::string &  name  )

private

Definition at line 17 of file THistSvc.icc.

                                                           {

  GlobalDirectoryRestore restore;

  std::string idr(id);
  removeDoubleSlash( idr );

  if (idr.find("/") == idr.length()) {
    m_log << MSG::ERROR << "Badly formed identifier \"" << idr << "\": "
        << "Must not end with a /" << endmsg;
    return StatusCode::FAILURE;
  }


  TFile *f = nullptr;
  std::string stream,rem;
  if (!findStream(idr, stream, rem, f)) {
    m_log << MSG::ERROR << "Could not register id: \"" << idr << "\""
        << endmsg;
    return StatusCode::FAILURE;
  }

  std::string uid = "/" + stream + "/" + rem;
  auto itr = m_uids.find(uid);
  if (itr != m_uids.end()) {
    m_log << MSG::ERROR << "already registered an object with identifier \""
      << idr << "\"" << endmsg;
    return StatusCode::FAILURE;
  }


  bool temp = false;
  if ( !f ) {
    temp = true;
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "Historgram with id \"" << idr << "\" is temporary"
        << endmsg;
  }


  TObject *to = nullptr;
  THistID hid;

  // check to see if this hist is to be read in;
  if (!temp && m_files.find(stream)->second.second == READ) {

    if (hist != 0) {
      m_log << MSG::WARNING <<  "Registering id: \"" << idr
        << "\" with non zero pointer!" << endmsg;
    }

    if (readHist_i(idr,hist).isFailure()) {
      m_log << MSG::ERROR <<  "Unable to read in hist" << endmsg;
      return StatusCode::FAILURE;
    }
    to = dynamic_cast<TObject*>(hist);
    hid = THistID(uid,temp,to,f,m_files.find(stream)->second.second);

  } else if (hist == 0) {
    m_log << MSG::ERROR << "Unable to read in hist with id: \""
      << idr << "\"" << endmsg;
    return StatusCode::FAILURE;

  } else {

    to = dynamic_cast<TObject*>(hist);
    if (to == 0) {
      m_log << MSG::ERROR << "Could not dcast to TObject. id: \"" << idr
        << "\"" << endmsg;
      return StatusCode::FAILURE;
    }

    auto oitr = m_tobjs.find(to);
    if (oitr != m_tobjs.end()) {
      m_log << MSG::ERROR << "already registered id: \"" << idr
        << "\" with identifier \"" << oitr->second.id << "\"" << endmsg;
      return StatusCode::FAILURE;
    }

    hid = THistID(uid,temp,to,f,m_files.find(stream)->second.second);
    TDirectory* dir = changeDir(hid);

    if ( dynamic_cast<TTree*>(hist) ) {
      dynamic_cast<TTree*>(hist)->SetDirectory(dir);
    } else if ( dynamic_cast<TH1*>(hist) ) {
      dynamic_cast<TH1*>(hist)->SetDirectory(dir);
    } else if ( dynamic_cast<TGraph*>(hist) ) {
      dir->Append(hist);
    } else {
      m_log << MSG::ERROR << "id: \"" << idr
        << "\" is not a TH, TTree, or TGraph. Attaching it to current dir."
        << endmsg;
      dir->Append(hist);
    }

  }

  std::string fname;
  if ( !f ) {
    fname = "none";
  } else {
    fname = f->GetName();
  }

  if (m_log.level() <= MSG::DEBUG) 
    m_log << MSG::DEBUG << "Registering " << System::typeinfoName(typeid(*hist))
      << " title: \"" << hist->GetTitle()
      << "\"  id: \"" << uid << "\"  dir: "
      //      << hist->GetDirectory()->GetPath() << "  "
      << changeDir(hid)->GetPath()
      << "  file: " << fname
      << endmsg;

  m_ids.emplace(rem, hid);
  m_uids[uid] = hid;
  m_tobjs[to] = hid;

  return StatusCode::SUCCESS;

}

StatusCode THistSvc::regTree ( const std::string &  name )

override

Definition at line 1004 of file THistSvc.cpp.

                                     {
  TTree *hist = nullptr;
  return regHist_i(hist, id);
}

StatusCode THistSvc::regTree ( const std::string &  name, TTree *  hist  )

override

Definition at line 1012 of file THistSvc.cpp.

                                                  {
  StatusCode sc = regHist_i(hist, id);
  if (hist && sc.isSuccess()) {
    if (m_autoSave != 0) hist->SetAutoSave(m_autoSave);
    hist->SetAutoFlush(m_autoFlush);
  }
  return sc;
}

StatusCode THistSvc::reinitialize ( )

override

Definition at line 230 of file THistSvc.cpp.

                       {

  GlobalDirectoryRestore restore;
  m_log << MSG::WARNING << "reinitialize not implemented" << endmsg;
  return StatusCode::SUCCESS;

}

void THistSvc::removeDoubleSlash ( std::string &  id ) const

private

Definition at line 1720 of file THistSvc.cpp.

                                               {

  while (id.find("//") != std::string::npos) {
    id.replace(id.find("//"),2,"/");
  }

}

StatusCode THistSvc::rootOpenAction ( FILEMGR_CALLBACK_ARGS  )

private

Definition at line 2030 of file THistSvc.cpp.

                                                                       {

  if (fa->tech() != Io::ROOT) {
    // This should never happen
    return StatusCode::SUCCESS;
  }

  if (fa->desc() != "HIST") {
    return StatusCode::SUCCESS;
  }

  p_incSvc->fireIncident(FileIncident(caller, "OpenHistFile", fa->name()));

  if ( fa->flags().isRead() ) {
    p_incSvc->fireIncident(FileIncident(caller, "BeginHistFile", fa->name()));
  } else if ( fa->flags().isWrite() ) {
    p_incSvc->fireIncident(FileIncident(caller, IncidentType::BeginOutputFile,
                    fa->name()));
  } else {
    // for Io::RW
    p_incSvc->fireIncident(FileIncident(caller, IncidentType::BeginOutputFile,
                    fa->name()));
  }

  return StatusCode::SUCCESS;


}

StatusCode THistSvc::rootOpenErrAction ( FILEMGR_CALLBACK_ARGS  )

private

Definition at line 2062 of file THistSvc.cpp.

                                                                          {

  if (fa->tech() != Io::ROOT) {
    // This should never happen
    return StatusCode::SUCCESS;
  }

  if (fa->desc() != "HIST") {
    return StatusCode::SUCCESS;
  }

  if ( fa->flags().isRead() ) {
    p_incSvc->fireIncident(FileIncident(caller, IncidentType::FailInputFile,
                    fa->name()));
  } else if ( fa->flags().isWrite() ) {
    p_incSvc->fireIncident(FileIncident(caller, IncidentType::FailOutputFile,
                    fa->name()));
  } else {
    // for Io::RW
    p_incSvc->fireIncident(FileIncident(caller, "FailRWFile", fa->name()));
  }

  return StatusCode::SUCCESS;
}

void THistSvc::setupCompressionLevel ( Property &  cmp )

private

Definition at line 1223 of file THistSvc.cpp.

{

  m_log << MSG::WARNING << "\"CompressionLevel\" Property has been deprecated. "
        << "Set it via the \"CL=\" parameter in the \"Output\" Property"
        << endmsg;

}

void THistSvc::setupInputFile ( Property &  inputfile )

private

call-back method to handle input stream property

Definition at line 1236 of file THistSvc.cpp.

{

  if (FSMState() < Gaudi::StateMachine::CONFIGURED || !m_okToConnect ) {

    m_log <<MSG::DEBUG << "Delaying connection of Input Files until Initialize"
      << ". now in " << FSMState()
      << endmsg;

    m_delayConnect = true;
  } else {

    m_log <<MSG::DEBUG << "Now connecting of Input Files"
      << endmsg;

    StatusCode sc = StatusCode::SUCCESS;

    for ( const auto& itr : m_inputfile.value() ) {
      if ( m_alreadyConnectedInFiles.end() != 
           m_alreadyConnectedInFiles.find( itr ) ) continue;
      if ( connect(itr).isFailure() ) {
        sc = StatusCode::FAILURE;
      } else {
        m_alreadyConnectedInFiles.insert( itr );
      }
      
    }

    if ( !sc.isSuccess() ) {
      throw GaudiException( "Problem connecting inputfile !!", name(),
                            StatusCode::FAILURE );
    }

  }
}

void THistSvc::setupOutputFile ( Property &  outputfile )

private

call-back method to handle output stream property

Definition at line 1275 of file THistSvc.cpp.

{
  if (FSMState() < Gaudi::StateMachine::CONFIGURED || !m_okToConnect) {
    m_log <<MSG::DEBUG << "Delaying connection of Input Files until Initialize"
      << ". now in " << FSMState()
      << endmsg;
    m_delayConnect = true;
  } else {

    StatusCode sc = StatusCode::SUCCESS;
    for ( const auto & itr : m_outputfile.value() ) {
      if ( m_alreadyConnectedOutFiles.end() !=
           m_alreadyConnectedOutFiles.find( itr ) ) continue;
      if ( connect(itr).isFailure() ) {
        sc = StatusCode::FAILURE;
      } else {
        m_alreadyConnectedOutFiles.insert( itr );
      }
    }

    if ( !sc.isSuccess() ) {
      throw GaudiException( "Problem connecting outputfile !!", name(),
                            StatusCode::FAILURE );
    }
  }
}

void THistSvc::updateFiles ( )

private

Definition at line 1305 of file THistSvc.cpp.

                      {

  // If TTrees grow beyond TTree::fgMaxTreeSize, a new file is
  // automatically created by root, and the old one closed. We
  // need to migrate all the UIDs over to show the correct file
  // pointer. This is ugly.

  if (m_log.level() <= MSG::DEBUG)
    m_log << MSG::DEBUG << "updateFiles()" << endmsg;


  for (auto uitr=m_uids.begin(); uitr != m_uids.end(); ++uitr) {
#ifndef NDEBUG
    if (m_log.level() <= MSG::VERBOSE)
      m_log << MSG::VERBOSE << " update: " << uitr->first << " "
        << uitr->second.id << " " << uitr->second.mode << endmsg;
#endif
    TObject* to = uitr->second.obj;
    TFile* oldFile = uitr->second.file;
    if (!to) {
      m_log << MSG::WARNING << uitr->first << ": TObject == 0" << endmsg;
    } else if ( uitr->second.temp || uitr->second.mode == READ ) {
      // do nothing - no need to check how big the file is since we
      // are just reading it.
#ifndef NDEBUG
    if (m_log.level() <= MSG::VERBOSE)
      m_log << MSG::VERBOSE << "     skipping" << endmsg;
#endif

    } else if (to->IsA()->InheritsFrom("TTree")) {
      TTree* tr = dynamic_cast<TTree*>(to);
      TFile* newFile = tr->GetCurrentFile();

      if (oldFile != newFile) {
        std::string newFileName = newFile->GetName();
        std::string oldFileName, streamName, rem;
        TFile* dummy = nullptr;
        findStream(uitr->second.id, streamName, rem, dummy);

        for (auto& itr : m_files ) {
          if (itr.second.first == oldFile) itr.second.first = newFile;
        }

        for (auto uitr2 = uitr; uitr2 != m_uids.end(); ++uitr2) {
          if (uitr2->second.file == oldFile) {
            uitr2->second.file = newFile;
          }
        }

        auto sitr = std::find_if( std::begin(m_fileStreams), std::end(m_fileStreams),
                                  [&](streamMap::const_reference s) {
                                      return s.second == streamName;
        });
        if (sitr!=std::end(m_fileStreams)) oldFileName = sitr->first;

#ifndef NDEBUG
      if (m_log.level() <= MSG::DEBUG)
        m_log << MSG::DEBUG << "migrating uid: " << uitr->second.id
              << "   stream: " << streamName
              << "   oldFile: " << oldFileName
              << "   newFile: " << newFileName
              << endmsg;
#endif


        if (!oldFileName.empty()) {
          auto i = m_fileStreams.lower_bound(oldFileName);
          while (i != std::end(m_fileStreams) && i->first == oldFileName) {

#ifndef NDEBUG
            if (m_log.level() <= MSG::DEBUG)
              m_log << MSG::DEBUG << "changing filename \"" << i->first
                    << "\" to \"" << newFileName << "\" for stream \""
                    << i->second << "\"" << endmsg;
#endif
            std::string nm = std::move(i->second);
            i = m_fileStreams.erase(i);
            m_fileStreams.emplace( newFileName, std::move(nm) ); 
          }


        } else {
          m_log << MSG::ERROR
                << "Problems updating fileStreams with new file name" << endmsg;
        }

      }
    }
  }
}

StatusCode THistSvc::write ( )

private

Definition at line 1399 of file THistSvc.cpp.

                {

  updateFiles();

  std::for_each( m_files.begin(), m_files.end(), [](std::pair<const std::string,std::pair<TFile*,Mode>>& i) {
      auto mode = i.second.second;
      auto file = i.second.first;
      if ( mode == WRITE || mode == UPDATE || mode == SHARE ) {
          file->Write("",TObject::kOverwrite);
      } else if ( mode == APPEND ) {
          file->Write("");
      }
  } );

  if (m_log.level() <= MSG::DEBUG) {
    m_log << MSG::DEBUG << "THistSvc::write()::List of Files connected in ROOT "
          << endmsg;
    TSeqCollection *filelist=gROOT->GetListOfFiles();
    for (int ii=0; ii<filelist->GetEntries(); ii++) {
        m_log << MSG::DEBUG
              << "THistSvc::write()::List of Files connected in ROOT: \""
              << filelist->At(ii)->GetName()<<"\""<<endmsg;
    }
  }

  return StatusCode::SUCCESS;

}

Member Data Documentation

std::set<std::string> THistSvc::m_alreadyConnectedInFiles

private

list of already connected files.

This is to keep track of files registered by the setupInputFile callback method

Definition at line 187 of file THistSvc.h.

std::set<std::string> THistSvc::m_alreadyConnectedOutFiles

private

list of already connected files.

This is to keep track of files registered by the setupOutputFile callback method

Definition at line 191 of file THistSvc.h.

IntegerProperty THistSvc::m_autoFlush

private

Definition at line 182 of file THistSvc.h.

IntegerProperty THistSvc::m_autoSave

private

Definition at line 182 of file THistSvc.h.

IntegerProperty THistSvc::m_compressionLevel

private

Definition at line 182 of file THistSvc.h.

std::string THistSvc::m_curstream

mutableprivate

Definition at line 211 of file THistSvc.h.

bool THistSvc::m_delayConnect = false

private

Definition at line 209 of file THistSvc.h.

std::map<std::string, std::pair<TFile*,Mode> > THistSvc::m_files

private

Definition at line 203 of file THistSvc.h.

streamMap THistSvc::m_fileStreams

private

Definition at line 204 of file THistSvc.h.

idMap THistSvc::m_ids

private

Definition at line 200 of file THistSvc.h.

StringArrayProperty THistSvc::m_inputfile

private

Definition at line 180 of file THistSvc.h.

MsgStream THistSvc::m_log

mutableprivate

Definition at line 178 of file THistSvc.h.

IntegerProperty THistSvc::m_maxFileSize

private

Definition at line 182 of file THistSvc.h.

bool THistSvc::m_okToConnect = false

private

Definition at line 209 of file THistSvc.h.

StringArrayProperty THistSvc::m_outputfile

private

Definition at line 180 of file THistSvc.h.

BooleanProperty THistSvc::m_print

private

Definition at line 183 of file THistSvc.h.

std::vector<std::string> THistSvc::m_Rstream

private

Definition at line 181 of file THistSvc.h.

std::map<std::string, std::string > THistSvc::m_sharedFiles

private

Definition at line 206 of file THistSvc.h.

objMap THistSvc::m_tobjs

private

Definition at line 201 of file THistSvc.h.

uidMap THistSvc::m_uids

private

Definition at line 199 of file THistSvc.h.

std::vector<std::string> THistSvc::m_Wstream

private

Definition at line 181 of file THistSvc.h.

IFileMgr* THistSvc::p_fileMgr = nullptr

private

Definition at line 214 of file THistSvc.h.

IIncidentSvc* THistSvc::p_incSvc = nullptr

private

Definition at line 213 of file THistSvc.h.

bool THistSvc::signaledStop = false

private

Definition at line 208 of file THistSvc.h.

---

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

• GaudiSvc/src/THistSvc/THistSvc.h
• GaudiSvc/src/THistSvc/THistSvc.cpp
• GaudiSvc/src/THistSvc/THistSvc.icc