Anp::HistMan Class Reference

#include <HistMan.h>

List of all members.

Public Types
typedef std::map< std::string, HistInfo >	InfoMap
typedef InfoMap::const_iterator	InfoIter
typedef std::map< std::string, InfoMap >	DirMap
typedef DirMap::const_iterator	DirIter
Public Member Functions
bool	ReadFile (const std::string file)
bool	KeyExists (int key, const std::string &dir=".") const
bool	KeyExists (const std::string &key, const std::string &dir=".") const
unsigned int	NMiss () const
const Hist1d< double >	CreateHist1d (int key, const std::string &dir=".")
const Hist1d< double >	CreateHist1d (const std::string &key, const std::string &dir=".")
TH1 *	CreateTH1 (int key, const std::string &dir=".")
TH2 *	CreateTH2 (int key, const std::string &dir=".")
TH1 *	CreateTH1 (const std::string &key, const std::string &dir=".")
TH2 *	CreateTH2 (const std::string &key, const std::string &dir=".")
TH2 *	GetTH2 (int xkey, int ykey, const std::string &dir=".")
TH2 *	GetTH2 (const std::string &xkey, const std::string &ykey, const std::string &dir=".")
void	Print (std::ostream &o=std::cout) const
Static Public Member Functions
HistMan &	Instance ()
Private Member Functions
	HistMan ()
	~HistMan ()
	HistMan (const HistMan &)
const HistMan &	operator= (const HistMan &)
const Hist1d< double >	CreateHist1d (const HistInfo &info)
TH1 *	CreateTH1 (const HistInfo &info, const std::string &opt="double")
TH2 *	CreateTH2 (const HistInfo &info, const std::string &opt="float")
bool	ReadHistogramBlock (TiXmlElement &hist)
bool	ReadHistogram (const TiXmlElement &hist, HistInfo &info)
bool	ReadAxis (const TiXmlElement &hist, AxisInfo &info)
const std::vector< int >	GetIntVec (const std::string &list) const
const std::vector< std::string >	GetStringVec (const std::string &list) const
const std::string	Convert (int key, int width=0) const
unsigned short	ComputeWidth (int key) const
Private Attributes
unsigned int	fNMiss
unsigned short	fKeyWidth
DirMap	fDirMap
Static Private Attributes
HistMan *	fgInstance = 0
	do not stream
Mutex	fMutex
	do not stream

---

Member Typedef Documentation

typedef DirMap::const_iterator Anp::HistMan::DirIter

Definition at line 95 of file PhysicsNtuple/Hist/HistMan.h. Referenced by CreateHist1d(), CreateTH1(), CreateTH2(), GetTH2(), KeyExists(), and Print().

typedef std::map<std::string, InfoMap> Anp::HistMan::DirMap

Definition at line 94 of file PhysicsNtuple/Hist/HistMan.h.

typedef InfoMap::const_iterator Anp::HistMan::InfoIter

Definition at line 93 of file PhysicsNtuple/Hist/HistMan.h. Referenced by CreateHist1d(), CreateTH1(), CreateTH2(), GetTH2(), and Print().

typedef std::map<std::string, HistInfo> Anp::HistMan::InfoMap

Definition at line 92 of file PhysicsNtuple/Hist/HistMan.h. Referenced by CreateHist1d(), CreateTH1(), CreateTH2(), GetTH2(), KeyExists(), Print(), and ReadHistogramBlock().

---

Constructor & Destructor Documentation

Anp::HistMan::HistMan (   )  [private]

Definition at line 56 of file PhysicsNtuple/Hist/HistMan.cxx. :fNMiss(0), fKeyWidth(4) { }

Anp::HistMan::~HistMan (   )  [private]

Definition at line 63 of file PhysicsNtuple/Hist/HistMan.cxx. { }

Anp::HistMan::HistMan (  const HistMan &   )  [private]

---

Member Function Documentation

unsigned short Anp::HistMan::ComputeWidth (  int  key  )  const [private]

Definition at line 865 of file PhysicsNtuple/Hist/HistMan.cxx. Referenced by Convert(), and ReadHistogramBlock(). { // // Compute key width // return static_cast<unsigned int>(std::floor(std::log10(double(std::abs(key) + 1)))) + 1; }

const string Anp::HistMan::Convert (  int  key, int  width = 0 )  const [private]

Definition at line 849 of file PhysicsNtuple/Hist/HistMan.cxx. References ComputeWidth(). Referenced by CreateTH1(), CreateTH2(), and GetTH2(). { stringstream strkey; if(width == 0) { strkey << setw(ComputeWidth(key)) << key; } else { strkey << setfill('0') << setw(width) << key; } return strkey.str(); }

const Anp::Hist1d< double > Anp::HistMan::CreateHist1d (  const HistInfo &  info  )  [private]

Definition at line 643 of file PhysicsNtuple/Hist/HistMan.cxx. References Anp::AxisInfo::GetBins(), Anp::AxisInfo::GetMax(), Anp::AxisInfo::GetMin(), Anp::AxisInfo::GetNbins(), and Anp::HistInfo::GetXaxis(). { const AxisInfo &xaxis = info.GetXaxis(); if(!xaxis.GetBins().empty()) { return Hist1d<double>(xaxis.GetBins()); } return Hist1d<double>(xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax()); }

const Anp::Hist1d< double > Anp::HistMan::CreateHist1d (  const std::string &  key, const std::string &  dir = "." )

Definition at line 466 of file PhysicsNtuple/Hist/HistMan.cxx. References DirIter, fDirMap, fMutex, InfoIter, and InfoMap. { Anp::Lock<Mutex> lock(fMutex); DirIter dit = fDirMap.find(dir); if(dit == fDirMap.end()) { ++fNMiss; return Hist1d<double>(0, 0.0, 0.0); } const InfoMap &imap = dit -> second; InfoIter hit = imap.find(key); if(hit == imap.end()) { ++fNMiss; return Hist1d<double>(0, 0.0, 0.0); } return HistMan::CreateHist1d(hit -> second); }

const Anp::Hist1d< double > Anp::HistMan::CreateHist1d (  int  key, const std::string &  dir = "." )

Definition at line 460 of file PhysicsNtuple/Hist/HistMan.cxx. { return HistMan::CreateHist1d(HistMan::Convert(key), dir); }

TH1 * Anp::HistMan::CreateTH1 (  const HistInfo &  info, const std::string &  opt = "double" )  [private]

Definition at line 656 of file PhysicsNtuple/Hist/HistMan.cxx. References Convert(), fKeyWidth, Anp::AxisInfo::GetBins(), Anp::HistInfo::GetKey(), Anp::AxisInfo::GetMax(), Anp::AxisInfo::GetMin(), Anp::HistInfo::GetName(), Anp::AxisInfo::GetNbins(), Anp::AxisInfo::GetTitle(), Anp::HistInfo::GetTitle(), Anp::HistInfo::GetXaxis(), Anp::HistInfo::GetYaxis(), and Anp::SetDir(). { string strkey = info.GetKey(); string name = info.GetName(); string title = info.GetTitle(); const int intkey = std::atoi(strkey.c_str()); if(intkey != 0) { strkey = "key_" + Convert(intkey, fKeyWidth); } if(name.empty()) { name = strkey; } if(title.empty()) { title = info.GetName(); } const AxisInfo &xaxis = info.GetXaxis(); const AxisInfo &yaxis = info.GetYaxis(); TH1 *h = 0; if(xaxis.GetBins().empty()) { Anp::Lock<Anp::Mutex> lock(Anp::GetROOTMutex()); if(opt.find("profile") != string::npos) { h = new TProfile(name.c_str(), title.c_str(), xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax()); } else if(opt.find("float") != string::npos) { h = new TH1F(name.c_str(), title.c_str(), xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax()); } else { h = new TH1D(name.c_str(), title.c_str(), xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax()); } } else { Anp::Lock<Anp::Mutex> lock(Anp::GetROOTMutex()); const vector<double> &bvec = xaxis.GetBins(); double *darray = new double[bvec.size()]; for(unsigned int ibin = 0; ibin < bvec.size(); ++ibin) { darray[ibin] = bvec[ibin]; } if(opt.find("profile") != string::npos) { h = new TProfile(name.c_str(), title.c_str(), bvec.size() - 1, darray); } else if(opt.find("float") != string::npos) { h = new TH1F(name.c_str(), title.c_str(), bvec.size() - 1, darray); } else { h = new TH1D(name.c_str(), title.c_str(), bvec.size() - 1, darray); } delete [] darray; } // // Set NULL TDirectory // Anp::SetDir(h, 0); if(xaxis.GetTitle().size() > 0) { h -> GetXaxis() -> SetTitle(xaxis.GetTitle().c_str()); h -> GetXaxis() -> CenterTitle(); } if(yaxis.GetTitle().size() > 0) { h -> GetYaxis() -> SetTitle(yaxis.GetTitle().c_str()); h -> GetYaxis() -> CenterTitle(); } return h; }

TH1 * Anp::HistMan::CreateTH1 (  const std::string &  key, const std::string &  dir = "." )

Definition at line 496 of file PhysicsNtuple/Hist/HistMan.cxx. References Anp::CreateTH1(), DirIter, fDirMap, fMutex, InfoIter, and InfoMap. { Anp::Lock<Mutex> lock(fMutex); // // Extract option // string opt; string dir = dir_; if (dir.find("profile") != string::npos) opt += "profile"; else if(dir.find("float") != string::npos) opt += "float"; if(!opt.empty() && dir.find(opt.c_str()) != string::npos) { dir.erase(dir.find(opt.c_str()), opt.size()); } if(dir.empty()) { dir = "."; } DirIter dit = fDirMap.find(dir); if(dit == fDirMap.end()) { ++fNMiss; return 0; } const InfoMap &imap = dit -> second; InfoIter hit = imap.find(key); if(hit == imap.end()) { ++fNMiss; return 0; } return HistMan::CreateTH1(hit -> second, opt); }

TH1 * Anp::HistMan::CreateTH1 (  int  key, const std::string &  dir = "." )

Definition at line 490 of file PhysicsNtuple/Hist/HistMan.cxx. References Anp::CreateTH1(). Referenced by Anp::MakeTH1(). { return HistMan::CreateTH1(HistMan::Convert(key), dir); }

TH2 * Anp::HistMan::CreateTH2 (  const HistInfo &  info, const std::string &  opt = "float" )  [private]

Definition at line 749 of file PhysicsNtuple/Hist/HistMan.cxx. References Convert(), fKeyWidth, Anp::HistInfo::GetKey(), Anp::AxisInfo::GetMax(), Anp::AxisInfo::GetMin(), Anp::HistInfo::GetName(), Anp::AxisInfo::GetNbins(), Anp::AxisInfo::GetTitle(), Anp::HistInfo::GetTitle(), Anp::HistInfo::GetXaxis(), Anp::HistInfo::GetYaxis(), and Anp::SetDir(). { string strkey = info.GetKey(); string name = info.GetName(); string title = info.GetTitle(); const int intkey = std::atoi(strkey.c_str()); if(intkey != 0) { strkey = Convert(intkey, fKeyWidth); } if(name.empty()) { name = strkey; } if(title.empty()) { title = info.GetName(); } const AxisInfo &xaxis = info.GetXaxis(); const AxisInfo &yaxis = info.GetYaxis(); TH2 *h = 0; if(!name.empty()) { Anp::Lock<Anp::Mutex> lock(Anp::GetROOTMutex()); if(opt.find("double") != string::npos) { h = new TH2D(name.c_str(), title.c_str(), xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax(), yaxis.GetNbins(), yaxis.GetMin(), yaxis.GetMax()); } else { h = new TH2F(name.c_str(), title.c_str(), xaxis.GetNbins(), xaxis.GetMin(), xaxis.GetMax(), yaxis.GetNbins(), yaxis.GetMin(), yaxis.GetMax()); } if(xaxis.GetTitle().size() > 0) { h -> GetXaxis() -> SetTitle(xaxis.GetTitle().c_str()); h -> GetXaxis() -> CenterTitle(); } if(yaxis.GetTitle().size() > 0) { h -> GetYaxis() -> SetTitle(yaxis.GetTitle().c_str()); h -> GetYaxis() -> CenterTitle(); } } // // Set NULL TDirectory // Anp::SetDir(h, 0); return h; }

TH2 * Anp::HistMan::CreateTH2 (  const std::string &  key, const std::string &  dir = "." )

Definition at line 544 of file PhysicsNtuple/Hist/HistMan.cxx. References Anp::Draw::CreateTH2(), DirIter, fDirMap, fMutex, InfoIter, and InfoMap. { Anp::Lock<Mutex> lock(fMutex); // // Extract option // string opt; string dir = dir_; if (dir.find("double") != string::npos) opt += "double"; if(!opt.empty() && dir.find(opt.c_str()) != string::npos) { dir.erase(dir.find(opt.c_str()), opt.size()); } if(dir.empty()) { dir = "."; } DirIter dit = fDirMap.find(dir); if(dit == fDirMap.end()) { ++fNMiss; return 0; } const InfoMap &imap = dit -> second; InfoIter hit = imap.find(key); if(hit == imap.end()) { ++fNMiss; return 0; } return HistMan::CreateTH2(hit -> second, opt); }

TH2 * Anp::HistMan::CreateTH2 (  int  key, const std::string &  dir = "." )

Definition at line 538 of file PhysicsNtuple/Hist/HistMan.cxx. References Convert(), and Anp::Draw::CreateTH2(). Referenced by Anp::MakeTH2(). { return HistMan::CreateTH2(Convert(key), dir); }

const vector< int > Anp::HistMan::GetIntVec (  const std::string &  list  )  const [private]

Definition at line 812 of file PhysicsNtuple/Hist/HistMan.cxx. References UtilString::StringTok(). Referenced by ReadHistogramBlock(). { vector<string> namelist; UtilString::StringTok(namelist, list, ", "); vector<int> ivec; for(vector<string>::const_iterator sit = namelist.begin(); sit != namelist.end(); ++sit) { const int newbase = std::atoi(sit -> c_str()); if(newbase != 0) { ivec.push_back(newbase); } } return ivec; }

const vector< string > Anp::HistMan::GetStringVec (  const std::string &  list  )  const [private]

Definition at line 831 of file PhysicsNtuple/Hist/HistMan.cxx. References UtilString::StringTok(). Referenced by ReadHistogramBlock(). { vector<string> namelist; UtilString::StringTok(namelist, list, ", "); vector<string> nvec; for(vector<string>::const_iterator sit = namelist.begin(); sit != namelist.end(); ++sit) { if(!sit -> empty()) { nvec.push_back(*sit); } } return nvec; }

TH2 * Anp::HistMan::GetTH2 (  const std::string &  xkey, const std::string &  ykey, const std::string &  dir = "." )

Definition at line 592 of file PhysicsNtuple/Hist/HistMan.cxx. References Anp::Draw::CreateTH2(), DirIter, fDirMap, Anp::HistInfo::fKey, fMutex, Anp::HistInfo::fName, Anp::HistInfo::fTitle, Anp::HistInfo::fXaxis, Anp::HistInfo::fYaxis, Anp::HistInfo::GetKey(), Anp::HistInfo::GetXaxis(), Anp::HistInfo::GetYaxis(), InfoIter, InfoMap, and Anp::AxisInfo::Valid(). { Anp::Lock<Mutex> lock(fMutex); DirIter dit = fDirMap.find(dir); if(dit == fDirMap.end()) { ++fNMiss; return 0; } const InfoMap &imap = dit -> second; InfoIter xhit = imap.find(xkey); if(xhit == imap.end()) { ++fNMiss; return 0; } InfoIter yhit = imap.find(ykey); if(yhit == imap.end()) { ++fNMiss; return 0; } const HistInfo &xinfo = xhit -> second; const HistInfo &yinfo = yhit -> second; HistInfo info; info.fKey = "keys_" + xinfo.GetKey() + "_and_" + yinfo.GetKey(); info.fXaxis = xinfo.GetXaxis(); info.fYaxis = yinfo.GetXaxis(); info.fName = info.GetKey(); info.fTitle = info.GetKey(); if(xinfo.GetYaxis().Valid()) { info.fXaxis = xinfo.GetYaxis(); } if(yinfo.GetYaxis().Valid()) { info.fYaxis = yinfo.GetYaxis(); } return HistMan::CreateTH2(info); }

TH2 * Anp::HistMan::GetTH2 (  int  xkey, int  ykey, const std::string &  dir = "." )

Definition at line 586 of file PhysicsNtuple/Hist/HistMan.cxx. References Convert(). { return HistMan::GetTH2(Convert(xkey), HistMan::Convert(ykey), dir); }

Anp::HistMan & Anp::HistMan::Instance (   )  [static]

Definition at line 38 of file PhysicsNtuple/Hist/HistMan.cxx. References fgInstance, and fMutex. Referenced by Anp::MakeTH1(), and Anp::MakeTH2(). { if(!fgInstance) { Anp::Lock<Mutex> lock(fMutex); if(!fgInstance) { fgInstance = new Anp::HistMan(); } } return *fgInstance; }

bool Anp::HistMan::KeyExists (  const std::string &  key, const std::string &  dir = "." )  const

Definition at line 444 of file PhysicsNtuple/Hist/HistMan.cxx. References DirIter, fDirMap, fMutex, and InfoMap. { Anp::Lock<Mutex> lock(fMutex); DirIter dit = fDirMap.find(dir); if(dit == fDirMap.end()) { return false; } const InfoMap &imap = dit -> second; return (imap.find(key) != imap.end()); }

bool Anp::HistMan::KeyExists (  int  key, const std::string &  dir = "." )  const

Definition at line 438 of file PhysicsNtuple/Hist/HistMan.cxx. { return HistMan::KeyExists(HistMan::Convert(key), dir); }

unsigned int Anp::HistMan::NMiss (   )  const

Definition at line 874 of file PhysicsNtuple/Hist/HistMan.cxx. References fMutex. { Anp::Lock<Mutex> lock(fMutex); return fNMiss; }

const HistMan& Anp::HistMan::operator= (  const HistMan &   )  [private]

void Anp::HistMan::Print (  std::ostream &  o = std::cout  )  const

Definition at line 881 of file PhysicsNtuple/Hist/HistMan.cxx. References DirIter, fDirMap, InfoIter, InfoMap, and Anp::AxisInfo::Print(). { o << "HistMan::Print" << endl; for(DirIter dit = fDirMap.begin(); dit != fDirMap.end(); ++dit) { o << "Directory = " << dit -> first << endl; const InfoMap &imap = dit -> second; for(InfoIter it = imap.begin(); it != imap.end(); ++it) { it -> second.Print(o); } } }

bool Anp::HistMan::ReadAxis (  const TiXmlElement &  hist, AxisInfo &  info )  [private]

Definition at line 373 of file PhysicsNtuple/Hist/HistMan.cxx. References TiXmlElement::Attribute(), Anp::AxisInfo::fBins, Anp::AxisInfo::fMax, Anp::AxisInfo::fMin, Anp::AxisInfo::fNbins, Anp::AxisInfo::fTitle, Anp::AxisInfo::fValid, max, min, TiXmlElement::QueryFloatAttribute(), TiXmlElement::QueryIntAttribute(), UtilString::StringTok(), and TIXML_SUCCESS. Referenced by ReadHistogram(). { bool read_bin = false, read_bins = false; int bin = 0; float min = 0.0, max = 0.0; if(axis.QueryIntAttribute("bin", &bin) == TIXML_SUCCESS && axis.QueryFloatAttribute("min", &min) == TIXML_SUCCESS && axis.QueryFloatAttribute("max", &max) == TIXML_SUCCESS) { read_bin = true; } const char *xtitle = axis.Attribute("title"); if(xtitle) { info.fTitle = xtitle; } vector<double> bvec; const char *bins = axis.Attribute("bins"); if(bins) { vector<string> datavec; UtilString::StringTok(datavec, bins, ", "); for(vector<string>::const_iterator sit = datavec.begin(); sit != datavec.end(); ++sit) { if(sit -> size() > 0) { bvec.push_back(std::atof(sit -> c_str())); } } std::sort(bvec.begin(), bvec.end()); //for(unsigned int i = 0; i < bvec.size(); ++i) //{ //cout << "bin = " << i << " edge = " << bvec[i] << endl; //} if(!bvec.empty()) { read_bins = true; } } if(read_bin) { info.fNbins = bin; info.fMin = min; info.fMax = max; info.fValid = true; } else if(read_bins) { info.fBins = bvec; } return true; }

bool Anp::HistMan::ReadFile (  const std::string  file  )

Definition at line 68 of file PhysicsNtuple/Hist/HistMan.cxx. References doc, fMutex, and ReadHistogramBlock(). { // // Parse xml file with histogram definitions. These are top level keys // which must be located within "histograms" object: // keybase - creates multiple keys for one histogram // dirbase - stores keys in "sub-directory" // h1d - histogram definition // Anp::Lock<Mutex> lock(fMutex); TiXmlDocument doc(file.c_str()); const bool loadok = doc.LoadFile(); if(!loadok) { cerr << "HistMan::ReadFile - failed to load: " << file << endl; return false; } for(TiXmlNode *pChild = doc.FirstChild(); pChild != 0; pChild = pChild -> NextSibling()) { if(!pChild) { cerr << "HistMan::ReadFile - null child..." << endl; continue; } if(pChild -> Type() != TiXmlNode::ELEMENT) { continue; } if(pChild -> Value() && strcmp(pChild -> Value(), "histograms")) { continue; } TiXmlElement* pHist = pChild -> ToElement(); if(pHist) { ReadHistogramBlock(*pHist); } else { cerr << "HistMan::ReadFile - failed to cast root element for histograms" << endl; continue; } } return true; }

bool Anp::HistMan::ReadHistogram (  const TiXmlElement &  hist, HistInfo &  info )  [private]

Definition at line 279 of file PhysicsNtuple/Hist/HistMan.cxx. References TiXmlElement::Attribute(), TiXmlNode::FirstChild(), Anp::HistInfo::fKey, Anp::HistInfo::fName, Anp::AxisInfo::fTitle, Anp::HistInfo::fTitle, Anp::HistInfo::fXaxis, Anp::HistInfo::fYaxis, ReadAxis(), and TiXmlNode::Value(). Referenced by ReadHistogramBlock(). { if(hist.Value() && strcmp(hist.Value(), "h1d") != 0) { return false; } AxisInfo xaxis, yaxis; bool axis1dok = false, axis2dok = false; for(const TiXmlNode *pChild = hist.FirstChild(); pChild != 0; pChild = pChild -> NextSibling()) { if(!pChild) { cerr << "Null child..." << endl; continue; } if(pChild -> Type() != TiXmlNode::ELEMENT) { continue; } const TiXmlElement* axis = pChild -> ToElement(); if(!axis) { cerr << "Failed to cast xml node to element..." << endl; continue; } if(!axis -> Value()) { cerr << "Node has no value" << endl; continue; } if(strcmp(axis -> Value(), "axis1d") == 0 && ReadAxis(*axis, xaxis)) { axis1dok = true; } if(strcmp(axis -> Value(), "axis2d") == 0 && ReadAxis(*axis, yaxis)) { axis2dok = true; } } if(!axis1dok) { cerr << "Failed to read 1d axis info" << endl; return false; } else { info.fXaxis = xaxis; } if(axis2dok) { info.fYaxis = yaxis; } const char *key = hist.Attribute("key"); if(key) { info.fKey = string(key); } else { cerr << "HistMan::ReadHistogram() - key string is empty" << endl; return false; } const char *name = hist.Attribute("name"); if(name) { info.fName = name; } const char *title = hist.Attribute("title"); if(title) { info.fTitle = title; } const char *axis = hist.Attribute("axis"); if(axis) { info.fXaxis.fTitle = axis; } return true; }

bool Anp::HistMan::ReadHistogramBlock (  TiXmlElement &  hist  )  [private]

Definition at line 121 of file PhysicsNtuple/Hist/HistMan.cxx. References base, ComputeWidth(), fDirMap, TiXmlNode::FirstChild(), Anp::HistInfo::fKey, fKeyWidth, GetIntVec(), Anp::HistInfo::GetKey(), GetStringVec(), Anp::HistInfo::GetTitle(), InfoMap, and ReadHistogram(). Referenced by ReadFile(). { InfoMap hmap; vector<vector<int> > basevec; vector<string> dirvec; for(TiXmlNode *pChild = block.FirstChild(); pChild != 0; pChild = pChild -> NextSibling()) { if(!pChild) { cerr << "HistMan::ReadFile - null child..." << endl; continue; } if(pChild -> Type() != TiXmlNode::ELEMENT) { continue; } TiXmlElement* hist = pChild -> ToElement(); if(!hist) { cerr << "HistMan::ReadFile - failed to cast xml node to element..." << endl; continue; } if(!hist -> Value()) { cerr << "HistMan::ReadFile - node has no value" << endl; continue; } if(strcmp(hist -> Value(), "keybase") == 0 && hist -> Attribute("keys")) { const vector<int> ivec = GetIntVec(hist -> Attribute("keys")); if(!ivec.empty()) { basevec.push_back(ivec); } continue; } if(strcmp(hist -> Value(), "dirbase") == 0 && hist -> Attribute("dirs")) { const vector<string> svec = GetStringVec(hist -> Attribute("dirs")); if(!svec.empty()) { dirvec.insert(dirvec.end(), svec.begin(), svec.end()); } continue; } HistInfo info; if(strcmp(hist -> Value(), "h1d") != 0 || !ReadHistogram(*hist, info)) { continue; } if(!hmap.insert(InfoMap::value_type(info.GetKey(), info)).second) { cerr << "HistMan::ReadFile - histogram key " << info.GetKey() << " already exists" << endl; } const int intkey = std::atoi(info.GetKey().c_str()); if(intkey != 0) { const unsigned int width = ComputeWidth(intkey); fKeyWidth = max<int>(width, fKeyWidth); if(width != info.GetKey().size()) { cerr << "HistMan::ReadFile - integer width and string size are not equal for key " <<info.GetKey()<<" "<<info.GetTitle()<< endl; } } } vector<int> allvec; for(vector<vector<int> >::const_iterator bit = basevec.begin(); bit != basevec.end(); ++bit) { const vector<int> prevvec = allvec; const vector<int> &currvec = *bit; for(vector<int>::const_iterator cit = currvec.begin(); cit != currvec.end(); ++cit) { if(prevvec.empty()) { allvec.push_back(*cit); } else { for(vector<int>::const_iterator pit = prevvec.begin(); pit != prevvec.end(); ++pit) { allvec.push_back(*cit + *pit); } } } } InfoMap nmap = hmap; for(vector<int>::const_iterator it = allvec.begin(); it != allvec.end(); ++it) { const int base = *it; for(InfoMap::const_iterator hit = hmap.begin(); hit != hmap.end(); ++hit) { const int oldkey = std::atoi((hit -> first).c_str()); if(oldkey == 0) { continue; } const int newkey = base + oldkey; HistInfo info = hit -> second; stringstream strkey; strkey << newkey; info.fKey = strkey.str(); if(!nmap.insert(InfoMap::value_type(info.GetKey(), info)).second) { cerr << "HistMan::ReadFile - histogram key " << info.GetKey() << " already exists" << endl; } const unsigned int width = ComputeWidth(newkey); fKeyWidth = max<int>(width, fKeyWidth); if(width != info.GetKey().size()) { cerr << "HistMan::ReadFile - integer width and string size are not equal" << endl; } } } if(dirvec.empty()) { dirvec.push_back("."); } for(vector<string>::const_iterator dit = dirvec.begin(); dit != dirvec.end(); ++dit) { InfoMap &imap = fDirMap[*dit]; for(InfoMap::const_iterator nit = nmap.begin(); nit != nmap.end(); ++nit) { if(!imap.insert(InfoMap::value_type(*nit)).second) { cerr << "HistMan::ReadFile - failed to add " << nit -> first << " to " << *dit << endl; } } } return true; }

---

Member Data Documentation

DirMap Anp::HistMan::fDirMap [private]

Definition at line 153 of file PhysicsNtuple/Hist/HistMan.h. Referenced by CreateHist1d(), CreateTH1(), CreateTH2(), GetTH2(), KeyExists(), Print(), and ReadHistogramBlock().

Anp::HistMan * Anp::HistMan::fgInstance = 0 [static, private]

do not stream Definition at line 27 of file PhysicsNtuple/Hist/HistMan.cxx. Referenced by Instance().

unsigned short Anp::HistMan::fKeyWidth [private]

Definition at line 151 of file PhysicsNtuple/Hist/HistMan.h. Referenced by CreateTH1(), CreateTH2(), and ReadHistogramBlock().

Anp::Mutex Anp::HistMan::fMutex [static, private]

do not stream Definition at line 28 of file PhysicsNtuple/Hist/HistMan.cxx. Referenced by CreateHist1d(), CreateTH1(), CreateTH2(), GetTH2(), Instance(), KeyExists(), NMiss(), and ReadFile().

unsigned int Anp::HistMan::fNMiss [private]

Definition at line 149 of file PhysicsNtuple/Hist/HistMan.h.

---

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

• PhysicsNtuple/Hist/HistMan.h
• PhysicsNtuple/Hist/HistMan.cxx

---