THistSvc Class Reference

#include <src/THistSvc/THistSvc.h>

Inheritance diagram for THistSvc:

Collaboration diagram for THistSvc:

Classes
class	GlobalDirectoryRestore
struct	THistID

Public Member Functions
virtual StatusCode	initialize ()
	Initialization (from CONFIGURED to INITIALIZED). More...
virtual StatusCode	reinitialize ()
	Initialization (from INITIALIZED or RUNNING to INITIALIZED, via CONFIGURED). More...
virtual StatusCode	finalize ()
	Finalize (from INITIALIZED to CONFIGURED). More...
virtual StatusCode	regHist (const std::string &name)
virtual StatusCode	regHist (const std::string &name, TH1 *)
virtual StatusCode	regHist (const std::string &name, TH2 *)
virtual StatusCode	regHist (const std::string &name, TH3 *)
virtual StatusCode	getHist (const std::string &name, TH1 *&) const
virtual StatusCode	getHist (const std::string &name, TH2 *&) const
virtual StatusCode	getHist (const std::string &name, TH3 *&) const
virtual StatusCode	regTree (const std::string &name)
virtual StatusCode	regTree (const std::string &name, TTree *)
virtual StatusCode	getTree (const std::string &name, TTree *&) const
virtual StatusCode	regGraph (const std::string &name)
virtual StatusCode	regGraph (const std::string &name, TGraph *)
virtual StatusCode	getGraph (const std::string &name, TGraph *&) const
virtual StatusCode	deReg (TObject *obj)
virtual StatusCode	deReg (const std::string &name)
virtual std::vector< std::string >	getHists () const
virtual std::vector< std::string >	getTrees () const
virtual std::vector< std::string >	getGraphs () const
virtual StatusCode	getTHists (TDirectory *td, TList &, bool recurse=false) const
virtual StatusCode	getTHists (const std::string &name, TList &, bool recurse=false) const
virtual StatusCode	getTHists (TDirectory *td, TList &tl, bool recurse=false, bool reg=false)
virtual StatusCode	getTHists (const std::string &name, TList &tl, bool recurse=false, bool reg=false)
virtual StatusCode	getTTrees (TDirectory *td, TList &, bool recurse=false) const
virtual StatusCode	getTTrees (const std::string &name, TList &, bool recurse=false) const
virtual StatusCode	getTTrees (TDirectory *td, TList &tl, bool recurse=false, bool reg=false)
virtual StatusCode	getTTrees (const std::string &name, TList &tl, bool recurse=false, bool reg=false)
virtual bool	exists (const std::string &name) const
	THistSvc (const std::string &name, ISvcLocator *svc)
void	handle (const Incident &)
	Inform that a new incident has occurred. More...
virtual StatusCode	io_reinit ()
	callback method to reinitialize the internal state of the component for I/O purposes (e.g. More...
Public Member Functions inherited from extends3< Service, ITHistSvc, IIncidentListener, IIoComponent >
	extends3 (A1 a1, A2 a2, A3 a3)
	Templated constructor with 3 arguments. More...
	extends3 (A1 a1, A2 a2)
	Templated constructor with 2 arguments. More...
	extends3 (A1 a1)
	Templated constructor with 1 argument. More...
	extends3 ()
	Default constructor. More...
virtual void *	i_cast (const InterfaceID &tid) const
	Implementation of IInterface::i_cast. More...
virtual StatusCode	queryInterface (const InterfaceID &ti, void **pp)
	Implementation of IInterface::queryInterface. More...
virtual std::vector< std::string >	getInterfaceNames () const
	Implementation of IInterface::getInterfaceNames. More...
virtual	~extends3 ()
	Virtual destructor. More...
Public Member Functions inherited from Service
virtual unsigned long	release ()
	Release Interface instance. More...
virtual const std::string &	name () const
	Retrieve name of the service. More...
virtual StatusCode	configure ()
	Configuration (from OFFLINE to CONFIGURED). More...
virtual StatusCode	start ()
	Start (from INITIALIZED to RUNNING). More...
virtual StatusCode	stop ()
	Stop (from RUNNING to INITIALIZED). More...
virtual StatusCode	terminate ()
	Initialization (from CONFIGURED to OFFLINE). More...
virtual Gaudi::StateMachine::State	FSMState () const
	Get the current state. More...
virtual Gaudi::StateMachine::State	targetFSMState () const
	When we are in the middle of a transition, get the state where the transition is leading us. More...
virtual StatusCode	restart ()
	Initialization (from RUNNING to RUNNING, via INITIALIZED). More...
virtual StatusCode	sysInitialize ()
	Initialize Service. More...
virtual StatusCode	sysStart ()
	Initialize Service. More...
virtual StatusCode	sysStop ()
	Initialize Service. More...
virtual StatusCode	sysFinalize ()
	Finalize Service. More...
virtual StatusCode	sysReinitialize ()
	Re-initialize the Service. More...
virtual StatusCode	sysRestart ()
	Re-initialize the Service. More...
virtual StatusCode	setProperty (const Property &p)
	Set the property by property. More...
virtual StatusCode	setProperty (const std::string &s)
	Set the property by string. More...
virtual StatusCode	setProperty (const std::string &n, const std::string &v)
	Set the property by std::string. More...
virtual StatusCode	getProperty (Property *p) const
	Get the property by property. More...
virtual const Property &	getProperty (const std::string &name) const
	Get the property by name. More...
virtual StatusCode	getProperty (const std::string &n, std::string &v) const
	Get the property by std::string. More...
virtual const std::vector < Property * > &	getProperties () const
	Get list of properties. More...
template<class TYPE >
StatusCode	setProperty (const std::string &name, const TYPE &value)
	set the property form the value More...
	Service (const std::string &name, ISvcLocator *svcloc)
	Standard Constructor. More...
SmartIF< ISvcLocator > &	serviceLocator () const
	Retrieve pointer to service locator. More...
StatusCode	setProperties ()
	Method for setting declared properties to the values specified for the job. More...
template<class T >
StatusCode	service (const std::string &name, const T *&psvc, bool createIf=true) const
	Access a service by name, creating it if it doesn't already exist. More...
template<class T >
StatusCode	service (const std::string &name, T *&psvc, bool createIf=true) const
template<class T >
StatusCode	service (const std::string &svcType, const std::string &svcName, T *&psvc) const
	Access a service by name and type, creating it if it doesn't already exist. More...
template<class T >
Property *	declareProperty (const std::string &name, T &property, const std::string &doc="none") const
	Declare the named property. More...
Property *	declareRemoteProperty (const std::string &name, IProperty *rsvc, const std::string &rname="") const
	Declare remote named properties. More...
SmartIF< IAuditorSvc > &	auditorSvc () const
	The standard auditor service.May not be invoked before sysInitialize() has been invoked. More...
Public Member Functions inherited from CommonMessaging< implements3< IService, IProperty, IStateful > >
	CommonMessaging (const A1 &a1, const A2 &a2, const A3 &a3)
	Templated constructor with 3 arguments. More...
	CommonMessaging (const A1 &a1, const A2 &a2)
	Templated constructor with 2 arguments. More...
	CommonMessaging (const A1 &a1)
	Templated constructor with 1 argument. More...
	CommonMessaging ()
	Default constructor. More...
virtual	~CommonMessaging ()
	Virtual destructor. More...
SmartIF< IMessageSvc > &	msgSvc () const
	The standard message service. More...
MsgStream &	msgStream () const
	Return an uninitialized MsgStream. More...
MsgStream &	msgStream (const MSG::Level level) const
	Predefined configurable message stream for the efficient printouts. More...
MsgStream &	always () const
	shortcut for the method msgStream(MSG::ALWAYS) More...
MsgStream &	fatal () const
	shortcut for the method msgStream(MSG::FATAL) More...
MsgStream &	err () const
	shortcut for the method msgStream(MSG::ERROR) More...
MsgStream &	error () const
	shortcut for the method msgStream(MSG::ERROR) More...
MsgStream &	warning () const
	shortcut for the method msgStream(MSG::WARNING) More...
MsgStream &	info () const
	shortcut for the method msgStream(MSG::INFO) More...
MsgStream &	debug () const
	shortcut for the method msgStream(MSG::DEBUG) More...
MsgStream &	verbose () const
	shortcut for the method msgStream(MSG::VERBOSE) More...
MsgStream &	msg () const
	shortcut for the method msgStream(MSG::INFO) More...
MSG::Level	msgLevel () const
	get the output level from the embedded MsgStream More...
bool	msgLevel (MSG::Level lvl) const
	get the output level from the embedded MsgStream More...
Public Member Functions inherited from extend_interfaces3< IService, IProperty, IStateful >
virtual	~extend_interfaces3 ()
	Virtual destructor. More...
Public Member Functions inherited from IService
	DeclareInterfaceID (IService, 3, 0)
	InterfaceID. More...
virtual	~IService ()
	virtual destructor More...
Public Member Functions inherited from INamedInterface
	DeclareInterfaceID (INamedInterface, 1, 0)
	InterfaceID. More...
virtual	~INamedInterface ()
	Virtual destructor (always needed for abstract classes). More...
Public Member Functions inherited from IInterface
virtual unsigned long	addRef ()=0
	Increment the reference count of Interface instance. More...
virtual unsigned long	refCount () const =0
	Current reference count. More...
virtual	~IInterface ()
	Virtual destructor. More...
Public Member Functions inherited from IProperty
	DeclareInterfaceID (IProperty, 2, 0)
	InterfaceID. More...
Public Member Functions inherited from IStateful
	DeclareInterfaceID (IStateful, 1, 0)
	InterfaceID. More...
virtual	~IStateful ()
Public Member Functions inherited from extend_interfaces3< ITHistSvc, IIncidentListener, IIoComponent >
virtual	~extend_interfaces3 ()
	Virtual destructor. More...
Public Member Functions inherited from ITHistSvc
	DeclareInterfaceID (ITHistSvc, 2, 0)
	InterfaceID. More...
virtual	~ITHistSvc ()
	virtual destructor More...
Public Member Functions inherited from IIncidentListener
	DeclareInterfaceID (IIncidentListener, 2, 0)
	InterfaceID. More...
Public Member Functions inherited from IIoComponent
	DeclareInterfaceID (IIoComponent, 1, 0)
virtual	~IIoComponent ()
	Destructor: More...
virtual StatusCode	io_finalize ()

Protected Member Functions
virtual	~THistSvc ()
Protected Member Functions inherited from Service
virtual	~Service ()
	Standard Destructor. More...
int	outputLevel () const
	get the Service's output level More...
Protected Member Functions inherited from CommonMessaging< implements3< IService, IProperty, IStateful > >
void	updateMsgStreamOutputLevel (int level)
	Update the output level of the cached MsgStream. More...

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
bool	m_delayConnect
bool	m_okToConnect
std::string	m_curstream
IIncidentSvc *	p_incSvc
IFileMgr *	p_fileMgr

Additional Inherited Members
Public Types inherited from extends3< Service, ITHistSvc, IIncidentListener, IIoComponent >
typedef extends3	base_class
	Typedef to this class. More...
typedef extend_interfaces3 < ITHistSvc, IIncidentListener, IIoComponent >	extend_interfaces_base
	Typedef to the base of this class. More...
typedef extend_interfaces_base::ext_iids	interfaces
	MPL set of all the implemented interfaces. More...
Public Types inherited from Service
typedef Gaudi::PluginService::Factory2 < IService *, const std::string &, ISvcLocator * >	Factory
Public Types inherited from CommonMessaging< implements3< IService, IProperty, IStateful > >
typedef CommonMessaging	base_class
Public Types inherited from implements3< IService, IProperty, IStateful >
typedef implements3	base_class
	Typedef to this class. More...
typedef extend_interfaces3 < IService, IProperty, IStateful >	extend_interfaces_base
	Typedef to the base of this class. More...
typedef extend_interfaces_base::ext_iids	interfaces
	MPL set of all the implemented interfaces. More...
Public Types inherited from extend_interfaces3< IService, IProperty, IStateful >
typedef mpl::fold< typename IService::iid::iids::type, typename mpl::fold< typename IProperty::iid::iids::type, typename IStateful::iid::iids::type, mpl::insert< mpl::_1, mpl::_2 > >::type, mpl::insert< mpl::_1, mpl::_2 > >::type	ext_iids
	MPL set of interfaces extended by this one. More...
Public Types inherited from IInterface
enum	Status { SUCCESS = 1, NO_INTERFACE, VERSMISMATCH, LAST_ERROR }
	Return status. More...
typedef Gaudi::InterfaceId < IInterface, 0, 0 >	iid
	Interface ID. More...
typedef mpl::set1< iid >	ext_iids
	Extra interfaces. More...
Public Types inherited from extend_interfaces3< ITHistSvc, IIncidentListener, IIoComponent >
typedef mpl::fold< typename ITHistSvc::iid::iids::type, typename mpl::fold< typename IIncidentListener::iid::iids::type, typename IIoComponent::iid::iids::type, mpl::insert< mpl::_1, mpl::_2 > >::type, mpl::insert< mpl::_1, mpl::_2 > >::type	ext_iids
	MPL set of interfaces extended by this one. More...
Static Public Member Functions inherited from IInterface
static const InterfaceID &	interfaceID ()
	Return an instance of InterfaceID identifying the interface. More...
Protected Attributes inherited from Service
IntegerProperty	m_outputLevel
	Service output level. More...
Gaudi::StateMachine::State	m_state
	Service state. More...
Gaudi::StateMachine::State	m_targetState
	Service state. More...
Protected Attributes inherited from CommonMessaging< implements3< IService, IProperty, IStateful > >
SmartIF< IMessageSvc >	m_msgsvc
	Pointer to the message service;. More...
std::auto_ptr< MsgStream >	m_msgStream
	The predefined message stream. More...
bool	m_streamWithService
	Flag to create a new MsgStream if it was created without the message service. 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 48 of file THistSvc.cpp.

  : base_class(name, svc), m_log(msgSvc(), name ), signaledStop(false),
    m_delayConnect(false),m_okToConnect(false),
    p_incSvc(0), p_fileMgr(0) {

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

protectedvirtual

Definition at line 67 of file THistSvc.cpp.

                    {

}

Member Function Documentation

bool THistSvc::browseTDir ( TDirectory *  dir ) const

private

Definition at line 413 of file THistSvc.cpp.

                                          {

  if (dir == 0) {
    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 == 0) { 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 == 0) { 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 1746 of file THistSvc.cpp.

                                                    {

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

  if (file != 0) {
    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 1512 of file THistSvc.cpp.

                                        {

  Tokenizer tok(true);

  string::size_type 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 val,VAL,TAG,filename,db_typ("ROOT");
  int cl(1);

  tok.analyse(ident.substr(loc+1,ident.length()), " ", "", "", "=", "'", "'");

  for ( Tokenizer::Items::iterator i = tok.items().begin(); i != tok.items().end(); i++)    {
    const std::string& tag = (*i).tag();
    TAG = tag;
    toupper(TAG);

    val = (*i).value();
    VAL = val;
    toupper(VAL);

    if (TAG == "FILE" || TAG == "DATAFILE") {
      filename = val;
      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: \"" << (*i).value() << "\""
            << endmsg;
        typ = 0;
      }
    } else if (TAG == "TYP") {
      db_typ = (*i).value();
    } else if (TAG == "CL") {
      cl = atoi(val.c_str());
    } else {
      props.push_back( Prop((*i).tag(), (*i).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()) {
    std::pair<streamMap::iterator, streamMap::iterator> fitr =
      m_fileStreams.equal_range(filename);

    std::string oldstream = (fitr.first)->second;

    std::pair<TFile*,Mode> 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(0);
  if (service("IncidentSvc",pi).isFailure()) {
    m_log << MSG::ERROR << "Unable to get the IncidentSvc" << endmsg;
    return StatusCode::FAILURE;
  }

  void* vf(0);
  TFile *f(0);

  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;
    stringstream out;
    string realfilename=filename;
    out << ishared++;
    filename = string("tmp_THistSvc_")+out.str()+string(".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 1969 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 != NULL && 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 )

virtual

Implements ITHistSvc.

Definition at line 929 of file THistSvc.cpp.

                            {

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

    uidMap::iterator 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);

    idMap::iterator itr3;
    bool found(false);

    std::pair<idMap::iterator, idMap::iterator> mitr = m_ids.equal_range(rem);
    if (mitr.first == mitr.second) {
      m_log << MSG::ERROR << "Problems deregistering TObject \""
            << obj->GetName()
            << "\" with id \"" << hid.id << "\"" << endmsg;
      return StatusCode::FAILURE;
    } else {
      for (itr3 = mitr.first; itr3 != mitr.second; ++itr3) {
        if (itr3->second.obj == obj) {
          found = true;
          break;
        }
      }
      if (!found) {
        m_log << MSG::ERROR << "Problems deregistering TObject \""
              << obj->GetName()
              << "\" with id \"" << hid.id << "\"" << endmsg;
      }
    }

    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 )

virtual

Implements ITHistSvc.

Definition at line 987 of file THistSvc.cpp.

                                   {

  uidMap::iterator 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 1775 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

virtual

Implements ITHistSvc.

Definition at line 1899 of file THistSvc.cpp.

                                              {

  TH1* h;

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


}

StatusCode THistSvc::finalize ( )

virtual

Finalize (from INITIALIZED to CONFIGURED).

Reimplemented from Service.

Definition at line 247 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) {
  uidMap::const_iterator uitr;
  for (uitr=m_uids.begin(); uitr != m_uids.end(); ++uitr) {

    TObject* to = uitr->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 (!uitr->second.temp) {
        dirname = uitr->second.file->GetPath();
        string id2(uitr->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: " << uitr->first << endmsg;
      } else {
        if (th->GetDirectory() != 0) {
          dirname = th->GetDirectory()->GetPath();
        }
      }
      } else if (! to ) {
    m_log << MSG::WARNING << uitr->first << " has NULL TObject ptr"
          << endmsg;
    }

      m_log << MSG::DEBUG << "uid: \"" << uitr->first << "\"  temp: "
            << uitr->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;
  map<string, pair<TFile*,Mode> >::const_iterator itr;
  for (itr = m_files.begin(); itr != m_files.end(); ++itr) {

    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(0);
    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(0);
      int r = p_fileMgr->open(Io::ROOT,name(), m_sharedFiles[itr->first],
                  Io::WRITE|Io::APPEND,vf,"HIST");

      if (r != 0 ) {
    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 = 0;
      r = p_fileMgr->open(Io::ROOT,name(),tmpfn,Io::READ,vf,"HIST");

      if (r != 0) {
    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 1210 of file THistSvc.cpp.

                                       {

  string::size_type pos = id.find("/");

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

    string::size_type 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 = 0;
    return true;
  }

  map< string,pair<TFile*,Mode> >::const_iterator itr = m_files.find(stream);
  if (itr != m_files.end()) {
    file = itr->second.first;
  } else {
    file = 0;
    m_log << MSG::WARNING << "no stream \"" << stream
          << "\" associated with id: \"" << id << "\""
          << endmsg;
  }

  return true;

}

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

virtual

Implements ITHistSvc.

Definition at line 1154 of file THistSvc.cpp.

                                                           {
  return getHist_i(id, hist);
}

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

virtual

Implements ITHistSvc.

Definition at line 1161 of file THistSvc.cpp.

                          {

  std::vector<std::string> names;
  names.reserve(m_uids.size());

  uidMap::const_iterator itr;
  for (itr = m_uids.begin(); itr != m_uids.end(); ++itr) {
    THistID tid = itr->second;

    if (tid.obj->IsA()->InheritsFrom("TGraph")) {
      names.push_back(itr->first);
    }

  }

  return names;

}

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

virtual

Implements ITHistSvc.

Definition at line 1084 of file THistSvc.cpp.

                                                       {
  return getHist_i(id, hist);
}

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

virtual

Implements ITHistSvc.

Definition at line 1091 of file THistSvc.cpp.

                                                       {
  return getHist_i(id, hist);
}

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

virtual

Implements ITHistSvc.

Definition at line 1098 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) {
    std::map<std::string, THistID>::const_iterator 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 = 0;
      return StatusCode::FAILURE;
    }

    THistID hid = itr->second;
    if (!quiet) {
      hist = dynamic_cast<T*>(hid.obj);
      if (hist == 0) {
    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 {

    std::pair< std::multimap<std::string,THistID>::const_iterator,
      std::multimap<std::string,THistID>::const_iterator > mitr = m_ids.equal_range(idr);


    if (mitr.first == mitr.second) {
      m_log << MSG::ERROR << "Could not locate Hist with id \"" << idr << "\""
        << endmsg;
      hist = 0;
      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 = 0;
      return StatusCode::FAILURE;
    } else {
      // return a SUCCESS if just INQUIRING
      m_log << MSG::INFO << "Found multiple matches with id \"" << idr 
          << "\"" << endmsg;
      hist = 0;
      return StatusCode::SUCCESS;
    }
      }
    }
  }
}

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

virtual

Implements ITHistSvc.

Definition at line 1105 of file THistSvc.cpp.

                         {

  std::vector<std::string> names;
  names.reserve(m_uids.size());

  uidMap::const_iterator itr;
  for (itr = m_uids.begin(); itr != m_uids.end(); ++itr) {
    THistID tid = itr->second;

    if (tid.obj->IsA()->InheritsFrom("TH1")) {
      names.push_back(itr->first);
    }

  }

  return names;

}

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

virtual

Implements ITHistSvc.

Definition at line 459 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 != 0 && 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

virtual

Implements ITHistSvc.

Definition at line 514 of file THistSvc.cpp.

                                                                    {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

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

  map< string,pair<TFile*,Mode> >::const_iterator 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  )

virtual

Implements ITHistSvc.

Definition at line 672 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);
      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 != 0 && 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  )

virtual

Implements ITHistSvc.

Definition at line 746 of file THistSvc.cpp.

                                                                        {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

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

  map< string,pair<TFile*,Mode> >::const_iterator 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 = "";
      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 {
    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

virtual

Implements ITHistSvc.

Definition at line 1126 of file THistSvc.cpp.

                                                         {
  return getHist_i(id, hist);
}

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

virtual

Implements ITHistSvc.

Definition at line 1133 of file THistSvc.cpp.

                         {

  std::vector<std::string> names;
  names.reserve(m_uids.size());

  uidMap::const_iterator itr;
  for (itr = m_uids.begin(); itr != m_uids.end(); ++itr) {
    THistID tid = itr->second;

    if (tid.obj->IsA()->InheritsFrom("TTree")) {
      names.push_back(itr->first);
    }

  }

  return names;

}

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

virtual

Implements ITHistSvc.

Definition at line 567 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 != 0 && 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

virtual

Implements ITHistSvc.

Definition at line 620 of file THistSvc.cpp.

                                                                    {
  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

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

  map< string,pair<TFile*,Mode> >::const_iterator 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);
    } 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;
    sc = StatusCode::FAILURE;
  } else {
    sc = getTTrees(gDirectory,tl,rcs);
  }

  return sc;


}

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

virtual

Implements ITHistSvc.

Definition at line 804 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);
      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 != 0 && 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  )

virtual

Implements ITHistSvc.

Definition at line 878 of file THistSvc.cpp.

                                                                        {

  GlobalDirectoryRestore restore;

  gErrorIgnoreLevel = kBreak;

  StatusCode sc;

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

  map< string,pair<TFile*,Mode> >::const_iterator 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;
    sc = StatusCode::FAILURE;
  } else {
    sc = getTTrees(gDirectory,tl,rcs,reg);
  }

  return sc;

}

void THistSvc::handle ( const Incident &  )

virtual

Inform that a new incident has occurred.

Implements IIncidentListener.

Definition at line 1911 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 (itr=m_files.begin(); itr!= m_files.end(); ++itr) {
    TFile* tf = itr->second.first;

#ifndef NDEBUG
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "stream: " << itr->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 \"" << itr->first
            << "\" has exceeded the max file size of "
            << m_maxFileSize.value() << "MB. Terminating Job."
            << endmsg;

      IEventProcessor* evt(0);
      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 \"" << itr->first
            << "\" is at 95% of its maximum allowable file size of "
            << m_maxFileSize.value() << "MB"
            << endmsg;
    }
  }
}

StatusCode THistSvc::initialize ( )

virtual

Initialization (from CONFIGURED to INITIALIZED).

Reimplemented from Service.

Definition at line 74 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.
  IJobOptionsSvc* jos(0);
  if( serviceLocator()->service( "JobOptionsSvc", jos, true ).isSuccess() ) {
    const std::vector<const Property*> *props = jos->getProperties( name() );

    if (props != NULL) {
      for ( std::vector<const Property*>::const_iterator cur = props->begin();
            cur != props->end(); cur++) {
        if ( (*cur)->name() == "OutputLevel" ) {
          setProperty( **cur ).ignore();
          m_log.setLevel( m_outputLevel.value() );
          break;
        }
      }
    }
  }


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

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

  vector<string>::const_iterator itr;
  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 ( 0 == 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

  Io::bfcn_action_t boa = boost::bind(&THistSvc::rootOpenAction, this, _1,_2);
  if (p_fileMgr->regAction(boa, Io::OPEN, Io::ROOT).isFailure()) {
    m_log << MSG::ERROR
      << "unable to register ROOT file open action with FileMgr"
      << endmsg;
  }
  Io::bfcn_action_t bea = boost::bind(&THistSvc::rootOpenErrAction, this, _1,_2);
  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 == true) {
    if (m_inputfile.value().size() > 0) { setupInputFile(m_inputfile); }
    if (m_outputfile.value().size() > 0) { setupOutputFile(m_outputfile); }

    m_delayConnect = false;

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


  IIoComponentMgr* iomgr(0);

  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;
      typedef std::map<std::string, std::pair<TFile*,Mode> > Registry_t;
      // register input/output files...
      for ( Registry_t::const_iterator
              ireg = m_files.begin(),
              iend = m_files.end();
            ireg != iend;
            ++ireg ) {
        const std::string fname = ireg->second.first->GetName();
        const IIoComponentMgr::IoMode::Type iomode =
          ( ireg->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 ( )

virtual

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

upon fork(2))

Implements IIoComponent.

Definition at line 2009 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(0);

  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 (FileReg_t::iterator ifile = m_files.begin(), iend=m_files.end();
       ifile != iend; ++ifile) {
    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 ( uidMap::iterator uid = m_uids.begin(), uid_end = m_uids.end();
          uid != uid_end;
          ++uid ) {
      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 1823 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 1257 of file THistSvc.cpp.

                                                                       {
  string::size_type pos = id.find("/");

  if (pos == string::npos) {
    root = "";
    rem = id;
    return;
  }

  if (pos == 0) {
    parseString(id.substr(1,id.length()),root,rem);
  } else {
    root = id.substr(0,pos);
    rem = id.substr(pos+1,id.length());
  }

}

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

private

Definition at line 1182 of file THistSvc.cpp.

                                                        {
  return readHist_i(id, hist);
}

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

private

Definition at line 1189 of file THistSvc.cpp.

                                                        {
  return readHist_i(id, hist);
}

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

private

Definition at line 1196 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 243 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 == 0) {
    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;
  gDirectory->GetObject(fdir.c_str(), to);

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



  hist = dynamic_cast<T*>(to);
  if (hist == 0) {
    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 1203 of file THistSvc.cpp.

                                                          {
  return readHist_i(id, hist);
}

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

virtual

Implements ITHistSvc.

Definition at line 1056 of file THistSvc.cpp.

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

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

virtual

Implements ITHistSvc.

Definition at line 1064 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 )

virtual

Implements ITHistSvc.

Definition at line 1004 of file THistSvc.cpp.

                                     {

  TH1 *hist(0);

  return regHist_i(hist, id);
}

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

virtual

Implements ITHistSvc.

Definition at line 1014 of file THistSvc.cpp.

                                                {
  return regHist_i(hist, id);
}

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

virtual

Implements ITHistSvc.

Definition at line 1021 of file THistSvc.cpp.

                                                {
  return regHist_i(hist, id);
}

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

virtual

Implements ITHistSvc.

Definition at line 1028 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(0);
  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;
  uidMap::const_iterator 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 == 0) {
    temp = true;
    if (m_log.level() <= MSG::DEBUG)
      m_log << MSG::DEBUG << "Historgram with id \"" << idr << "\" is temporary"
        << endmsg;
  }


  TObject *to;
  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;
    }

    objMap::const_iterator 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) != 0 ) {
      dynamic_cast<TTree*>(hist)->SetDirectory(dir);
    } else if ( dynamic_cast<TH1*>(hist) != 0 ) {
      dynamic_cast<TH1*>(hist)->SetDirectory(dir);
    } else if ( dynamic_cast<TGraph*>(hist) != 0 ) {
      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 == 0) {
    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.insert(std::pair<std::string,THistID>(rem, hid));
  m_uids[uid] = hid;
  m_tobjs[to] = hid;

  return StatusCode::SUCCESS;

}

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

virtual

Implements ITHistSvc.

Definition at line 1035 of file THistSvc.cpp.

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

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

virtual

Implements ITHistSvc.

Definition at line 1043 of file THistSvc.cpp.

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

StatusCode THistSvc::reinitialize ( )

virtual

Initialization (from INITIALIZED or RUNNING to INITIALIZED, via CONFIGURED).

Reimplemented from Service.

Definition at line 233 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 1813 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 2132 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 2164 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 1278 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 1291 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;

  typedef std::vector<std::string> Strings_t;
  for ( Strings_t::const_iterator
          itr  = m_inputfile.value().begin(),
          iEnd = m_inputfile.value().end();
        itr != iEnd;
        ++itr ) {
    if ( m_alreadyConnectedInFiles.end() ==
         m_alreadyConnectedInFiles.find( *itr ) ) {
      if ( connect(*itr).isFailure() ) {
        sc = StatusCode::FAILURE;
      } else {
        m_alreadyConnectedInFiles.insert( *itr );
      }
    }
  }

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

  }

  return;
}

void THistSvc::setupOutputFile ( Property &  outputfile )

private

call-back method to handle output stream property

Definition at line 1337 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;

  typedef std::vector<std::string> Strings_t;
  for ( Strings_t::const_iterator
          itr  = m_outputfile.value().begin(),
          iEnd = m_outputfile.value().end();
        itr != iEnd;
        ++itr ) {
    if ( m_alreadyConnectedOutFiles.end() ==
         m_alreadyConnectedOutFiles.find( *itr ) ) {
      if ( connect(*itr).isFailure() ) {
        sc = StatusCode::FAILURE;
      } else {
        m_alreadyConnectedOutFiles.insert( *itr );
      }
    }
  }

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

void THistSvc::updateFiles ( )

private

Definition at line 1375 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;


  uidMap::iterator uitr, uitr2;
  for (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;
        findStream(uitr->second.id, streamName, rem, dummy);

        map<string, pair<TFile*,Mode> >::iterator itr;
        for (itr=m_files.begin(); itr!= m_files.end(); ++itr) {
          if (itr->second.first == oldFile) {
            itr->second.first = newFile;

          }
        }

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

        streamMap::iterator sitr;
        for (sitr = m_fileStreams.begin(); sitr!=m_fileStreams.end(); ++sitr) {
          if (sitr->second == streamName) {
            oldFileName = sitr->first;
            break;
          }
        }


#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 != "") {
          while ( (sitr=m_fileStreams.find(oldFileName)) != m_fileStreams.end() ) {

#ifndef NDEBUG
            if (m_log.level() <= MSG::DEBUG)
              m_log << MSG::DEBUG << "changing filename \"" << oldFileName
                    << "\" to \"" << newFileName << "\" for stream \""
                    << sitr->second << "\"" << endmsg;
#endif
            m_fileStreams.erase(sitr);
        m_fileStreams.insert( make_pair(newFileName,streamName) );
          }


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

      }

    }
  }

}

StatusCode THistSvc::write ( )

private

Definition at line 1478 of file THistSvc.cpp.

                {

  updateFiles();

  map<string, pair<TFile*,Mode> >::const_iterator itr;
  for (itr=m_files.begin(); itr!= m_files.end(); ++itr) {
    if (itr->second.second == WRITE || itr->second.second == UPDATE
        ||itr->second.second==SHARE) {
      itr->second.first->Write("",TObject::kOverwrite);
    } else if (itr->second.second == APPEND) {
      itr->second.first->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++) {
    if (m_log.level() <= MSG::DEBUG)
      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

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

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

private

Definition at line 214 of file THistSvc.h.

IIncidentSvc* THistSvc::p_incSvc

private

Definition at line 213 of file THistSvc.h.

bool THistSvc::signaledStop

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