DiskResidentEventList Class Reference

#include <DiskResidentEventList.h>

Inheritance diagram for DiskResidentEventList:

Public Member Functions
	DiskResidentEventList ()
	DiskResidentEventList (const std::string &fname, const std::string &opt="UPDATE")
	DiskResidentEventList (const std::string &fname, int scales, const std::string &treeName, const std::string &opt="UPDATE")
	DiskResidentEventList (const DalitzEventPattern &pat)
	DiskResidentEventList (const DalitzEventPattern &pat, const std::string &fname, const std::string &opt="RECREATE")
	DiskResidentEventList (const MINT::IMinimalEventList< DalitzEvent > &otherList)
	DiskResidentEventList (const MINT::IMinimalEventList< DalitzEvent > &otherList, const std::string &newFname, const std::string &opt="RECREATE")
virtual	~DiskResidentEventList ()
virtual unsigned int	size () const
virtual bool	empty () const
DalitzEvent	getEvent (unsigned int i) const
virtual bool	Add (const DalitzEvent &evt)
bool	Add (const IDalitzEvent &evt)
bool	Add (const IDalitzEvent *evt)
bool	Add (const MINT::counted_ptr< IDalitzEvent > &evt)
bool	Add (const MINT::IMinimalEventList< DalitzEvent > &otherList)
bool	save ()
bool	Close ()
DalitzHistoSet	histoSet () const
DalitzHistoSet	weightedHistoSet () const
DalitzHistoSet	reWeightedHistoSet (MINT::IReturnRealForEvent< IDalitzEvent > *w) const
DalitzHistoSet	weighedReWeightedHistoSet (MINT::IReturnRealForEvent< IDalitzEvent > *w) const
bool	makePlots (const std::string &filename)
DalitzEvent	operator() (unsigned int i)

Static Public Member Functions
static std::string	generateFname ()

Protected Member Functions
std::string	cName ()
std::string	ntpName ()
bool	openFile ()
bool	makeNtp (const DalitzEvent &evt)
bool	makeNewFile ()
bool	fromFile ()
bool	init ()
Protected Member Functions inherited from MINT::IMinimalEventList< DalitzEvent >
	IMinimalEventList ()

Protected Attributes
std::string	_fname
std::string	_opt
TFile *	_f
TNtupleD *	_ntp
double	_scaleData
bool	_initialised

Static Protected Attributes
static int	__maxBytes = 10000000

Private Attributes
std::string	_cName
std::string	_ntpName

Detailed Description

Definition at line 26 of file DiskResidentEventList.h.

Constructor & Destructor Documentation

DiskResidentEventList() [1/7]

DiskResidentEventList::DiskResidentEventList

(

)

Definition at line 45 of file DiskResidentEventList.cpp.

  : _fname(generateFname())
  , _opt("UPDATE")
  , _f(0)//new TFile(generateFname().c_str(), "UPDATE"))
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  openFile();
  init();
}

DiskResidentEventList() [2/7]

DiskResidentEventList::DiskResidentEventList	(	const std::string &	fname,
		const std::string &	opt = "UPDATE"
	)

Definition at line 79 of file DiskResidentEventList.cpp.

  : _fname(fname)
  , _opt(boost::to_upper_copy(opt))
  , _f(0)//new TFile(fname.c_str(), opt.c_str()))
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  bool dbThis=false;
  openFile();
  cout << "DiskResidentEventList created with file " << fname << endl;
  fromFile();
  if(dbThis)cout << " from file done, now init: " << endl;
  init();
  if(dbThis)cout << " finished creation" << endl;
}

DiskResidentEventList() [3/7]

DiskResidentEventList::DiskResidentEventList	(	const std::string &	fname,
		int	scales,
		const std::string &	treeName,
		const std::string &	opt = "UPDATE"
	)

Definition at line 58 of file DiskResidentEventList.cpp.

  : _fname(fname)
  , _opt(boost::to_upper_copy(opt))
  , _f(0)//new TFile(fname.c_str(), opt.c_str()))
  , _ntp(0)
  , _scaleData(scale)
  , _cName(treeName)
  , _ntpName(treeName+"NTP")
{
  bool dbThis=true;
  openFile();
  cout << "DiskResidentEventList created with file " << fname << endl;
  fromFile();
  if(dbThis)cout << " from file done, now init: " << endl;
  init();
  if(dbThis)cout << " finished creation" << endl;
}

DiskResidentEventList() [4/7]

DiskResidentEventList::DiskResidentEventList

(

const DalitzEventPattern &

pat

)

Definition at line 177 of file DiskResidentEventList.cpp.

  : _fname(generateFname())
  , _opt("UPDATE")
  , _f(0)//new TFile(generateFname().c_str(), "RECREATE"))
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  openFile();
  DalitzEvent evt(pat);
  makeNtp(evt); // evt won't be saved, just to use "makeNtupleVarname" (clumsy)
  init();
}

DiskResidentEventList() [5/7]

DiskResidentEventList::DiskResidentEventList	(	const DalitzEventPattern &	pat,
		const std::string &	fname,
		const std::string &	opt = "RECREATE"
	)

Definition at line 192 of file DiskResidentEventList.cpp.

  : _fname(fname)
  , _opt(boost::to_upper_copy(opt))
  , _f(0)//new TFile(fname.c_str(), opt.c_str()))
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  openFile();
  fromFile();
  DalitzEvent evt(pat);
  if(0 == _ntp) makeNtp(evt);
  init();
}

DiskResidentEventList() [6/7]

DiskResidentEventList::DiskResidentEventList

(

const MINT::IMinimalEventList< DalitzEvent > &

otherList

)

Definition at line 99 of file DiskResidentEventList.cpp.

  : _fname(generateFname())
  , _opt("UPDATE")
  , _f(0)//new TFile(generateFname().c_str(), "UPDATE"))
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  openFile();
  Add(otherList);
  init();
}

DiskResidentEventList() [7/7]

DiskResidentEventList::DiskResidentEventList	(	const MINT::IMinimalEventList< DalitzEvent > &	otherList,
		const std::string &	newFname,
		const std::string &	opt = "RECREATE"
	)

Definition at line 125 of file DiskResidentEventList.cpp.

  : _fname(newFname)
  , _opt(boost::to_upper_copy(opt))
  , _f(0)//new TFile(newFname.c_str(), opt.c_str()))
    //  , _counted_ntp(0)
  , _ntp(0)
  , _scaleData(1)
  , _cName("DalitzEventList")
  , _ntpName("DalitzEventList")
{
  openFile();
  cout << " copy with new filename : " << newFname << endl;
  cout << " file pointer: " << _f << endl;
  fromFile();
  Add(otherList);
  init();
}

~DiskResidentEventList()

DiskResidentEventList::~DiskResidentEventList ( )

virtual

Definition at line 211 of file DiskResidentEventList.cpp.

                                             {
  this->Close();
}

Member Function Documentation

Add() [1/5]

bool DiskResidentEventList::Add ( const DalitzEvent &  evt )

virtual

Implements MINT::IMinimalEventList< DalitzEvent >.

Definition at line 322 of file DiskResidentEventList.cpp.

                                                     {
  if(0 == _f) makeNewFile();
  _f->cd();
  if(0 == _ntp) makeNtp(evt);
  unsigned int arraySize = evt.ntupleVarArraySize();
  vector<Double_t> array(arraySize);
  
  bool success = evt.fillNtupleVarArray(array);
  if(! success){
    cout << "ERROR in DiskResidentEventList::Add(const DalitzEvent& evt)"
     << ", call to DalitzEvent::fillNtupleVarArray"
     << " returned failure for event:\n" << evt
     << endl;
  }else{
    //_f = _ntp->GetCurrentFile(); 
    //_f->cd();
    _ntp->Fill(&(array[0]));
    _f = _ntp->GetCurrentFile(); 
    // above lines: to pick up when TNtuple (TTree) changes
    // files in case the current one is too big. All this
    // rubbish is necessary because of root's crazy way of
    // handling this.
  }

  return true;
}

Add() [2/5]

bool DiskResidentEventList::Add

(

const IDalitzEvent &

evt

)

Definition at line 349 of file DiskResidentEventList.cpp.

                                                      {
  DalitzEvent nevt(&evt);
  return Add(nevt);
}

Add() [3/5]

bool DiskResidentEventList::Add

(

const IDalitzEvent *

evt

)

Definition at line 353 of file DiskResidentEventList.cpp.

                                                      {
  DalitzEvent nevt(evt);
  return Add(nevt);
}

Add() [4/5]

bool DiskResidentEventList::Add

(

const MINT::counted_ptr< IDalitzEvent > &

evt

)

Definition at line 357 of file DiskResidentEventList.cpp.

                                                                   {
  DalitzEvent nevt(evt.get());
  return Add(nevt);
}

Add() [5/5]

bool DiskResidentEventList::Add

(

const MINT::IMinimalEventList< DalitzEvent > &

otherList

)

Definition at line 361 of file DiskResidentEventList.cpp.

                                                                              {
  if(0 == otherList.size()) return false;
  for(unsigned int i=0; i < otherList.size(); i++){
    Add(otherList.getEvent(i));
  }
  return true;
}

Close()

bool DiskResidentEventList::Close

(

)

Definition at line 397 of file DiskResidentEventList.cpp.

                                 {
  bool dbThis=false;
  bool success = true;

  if(dbThis){
    cout << "DiskResidentEventList::Close() called"
     << " for filename = " << _fname << endl;
  }
  if(0 == _f) return false;
  if(0 == _ntp) return false;
  // Don't understand why sometimes the TFile is already closed by the time the destructor is called,
  // but this avoids crases since _ntp is invalidated.
  if(!_f->IsOpen()) return false ;

  if(_opt.size() != 0 && _opt != "OPEN" && _opt != "READ"){
    success &= save();
    if(dbThis && ! success){
    cout << "DiskResidentEventList::Close() failure in saving when closing" 
        << endl;
    }
    if(dbThis) cout << "saved file with success = " << success << endl;
  }
  //_f->cd();
  //_ntp->Write();
  //_f->Write();
  if(dbThis) {
    cout << "_f " << _f << " _ntp " << _ntp << endl ;
    cout << "_f->IsOpen() " << _f->IsOpen() << endl ;
    cout << "mustcleanup " << _f->TestBit(kMustCleanup) << endl ;
    cout << "_ntp->GetName() " << _ntp->GetName() << endl ;
    cout << "getting current file" << endl;
  }
  _f = _ntp->GetCurrentFile();
  if(dbThis) cout << "_f = " << _f << endl;
  _f->cd();
  if(dbThis) cout << "cd'ed to _f and now about to Write" << endl;
  //  _f->Write();
  if(dbThis) cout << "written,now closing" <<endl;
  _f->Close();
  if(dbThis) cout << "closed." << endl;
  //delete _ntp;
  //delete _f;
  _f   = 0;
  _ntp = 0;

  (void)dbThis;
  return success;
}

cName()

std::string DiskResidentEventList::cName ( )

protected

Definition at line 23 of file DiskResidentEventList.cpp.

                                      {
  return _cName;}

empty()

bool DiskResidentEventList::empty ( ) const

virtual

Definition at line 291 of file DiskResidentEventList.cpp.

                                       {
  return (0 == this->size());
}

fromFile()

bool DiskResidentEventList::fromFile ( )

protected

Definition at line 246 of file DiskResidentEventList.cpp.

                                    {
  bool dbThis=false;
  if(dbThis) cout << "DiskResidentEventList::fromFile() called" << endl;
  if(0 == _f) makeNewFile();
  if(_f->IsZombie()){
    cout << "WARNING in DiskResidentEventList::fromFile():"
     << " current file is zombie - making new one" << endl;
    makeNewFile();
  }
  if(dbThis) cout << " opened file" << endl;
  if(dbThis) _f->ls(); cout << " that's in it" << endl;
  _f->cd();
  if(dbThis){
    cout << " cd'ed to file " << endl;
    cout << " now calling " <<  "_f->Get(" << cName() << ")" << endl;
  }
  _ntp = (TNtupleD*) _f->Get(cName().c_str());
  if(dbThis) cout << " got ntuple : " << _ntp << endl;
  if(0 == _ntp){
    return false;
  }else{
    if(dbThis){
      cout << "Success in DiskResidentEventList::fromFile(): "
       << " found ntuple: " << _ntp->ClassName() << ", " << cName()
       << " in file: ";
      _f->Print();
      cout << endl;
    }
  }
  if(dbThis) cout << " now calling init" << endl;
  init();
  if(dbThis) cout << "DiskResidentEventList::fromFile() returning 'true'" 
          << endl;
  return true;
}

generateFname()

std::string DiskResidentEventList::generateFname ( )

static

Definition at line 29 of file DiskResidentEventList.cpp.

                                              {
  string baseName = "DiskResidentEventList";
  string name = baseName;
  bool exists = true;
  ofstream checkFile;
  int i=0;
  do{
    name = baseName + "_" + anythingToString(i++) + ".root";
    checkFile.open( name.c_str(), ios::in );
    exists = ! checkFile.fail();
    checkFile.close();
  }while(exists);

  return name;
}

getEvent()

DalitzEvent DiskResidentEventList::getEvent ( unsigned int  i ) const

virtual

Implements MINT::IMinimalEventList< DalitzEvent >.

Definition at line 447 of file DiskResidentEventList.cpp.

                                                               {
  if(i > size()){
    cout << "FATAL ERROR in DiskResidentEventList::operator()"
     << " index i=" << i << " out of range " << size()
     << endl;
    throw "index out of range";
  }
  if(0 == _ntp){
    cout << "FATAL ERROR in DiskResidentEventList::operator()"
     << "there is no ntuple "
     << endl;
  }
  _ntp->GetEvent(i);
  return DalitzEvent(_ntp);
}

histoSet()

DalitzHistoSet DiskResidentEventList::histoSet

(

)

const

Definition at line 463 of file DiskResidentEventList.cpp.

                                                    {
  DalitzHistoSet hs;
  for(unsigned int i=0; i< this->size(); i++){
    hs.addEvent(this->getEvent(i));
  }
  return hs;
}

init()

bool DiskResidentEventList::init ( )

protected

Definition at line 282 of file DiskResidentEventList.cpp.

                                {
  if(0 != _ntp) _ntp->SetAutoSave(__maxBytes);
  return (0 != _ntp);
}

makeNewFile()

bool DiskResidentEventList::makeNewFile ( )

protected

Definition at line 295 of file DiskResidentEventList.cpp.

                                       {
  _f = TFile::Open(generateFname().c_str(), "RECREATE");
  // owned by gDirectory, I don't delete this myself.
  return (0 != _f);
}

makeNtp()

bool DiskResidentEventList::makeNtp ( const DalitzEvent &  evt )

protected

Definition at line 301 of file DiskResidentEventList.cpp.

                                                         {
  if(0 == _f){
    cout << "WARNING DiskResidentEventList::makeNtp - no file. "
     << "That shouldn't happen. I'll have one made." << endl;
    makeNewFile();
  }
  _f->cd();
  _ntp = (TNtupleD*) _f->FindObject(ntpName().c_str());
  if(0 == _ntp){
    _ntp = new TNtupleD(cName().c_str(), ntpName().c_str()
            , evt.makeNtupleVarnames().c_str()
            );
    cout << "made new ntuple" << endl;
  }else{
    cout << "found old ntuple" << endl;
  }
  if(0 != _ntp){
    _ntp->SetAutoSave(__maxBytes);
  }
  return 0 != _ntp;
}

makePlots()

bool DiskResidentEventList::makePlots

(

const std::string &

filename

)

Definition at line 502 of file DiskResidentEventList.cpp.

                                                              {
  histoSet().save(filename);
  return true;
}

ntpName()

std::string DiskResidentEventList::ntpName ( )

protected

Definition at line 26 of file DiskResidentEventList.cpp.

                                        {
  return _ntpName;}

openFile()

bool DiskResidentEventList::openFile ( )

protected

Definition at line 215 of file DiskResidentEventList.cpp.

                                    {
  bool dbThis=true;
  if(dbThis){
    cout << "Hello from DiskResidentEventList::openFile()" << endl;
    cout << " calling: TFile(" 
     << _fname << ", " << _opt << ", " << _fname << ")"
     << endl;
  }
  if(0 == _f)_f = TFile::Open(_fname.c_str(), _opt.c_str()
                   , _fname.c_str()
                   );
  
  if(0 == _f) return makeNewFile();
  if(_f->IsZombie()) return makeNewFile();
  if(! _f->IsOpen()) return makeNewFile();
  if(! _f->IsWritable()) {
    //return makeNewFile();
    cout << "Warning in DiskResidentEventList: file "
     << _fname << " not writable!"
     << endl;
  }
  if(dbThis){
    cout << " got _f = " << _f << endl;
    cout << " _f->ls() " << endl;
    _f->ls();
    cout << "?" << endl;
  }

//  if (1 ==1) return makeNewFile();
  return true;
}

operator()()

DalitzEvent DiskResidentEventList::operator()

(

unsigned int

i

)

reWeightedHistoSet()

DalitzHistoSet DiskResidentEventList::reWeightedHistoSet

(

MINT::IReturnRealForEvent< IDalitzEvent > *

w

)

const

Definition at line 480 of file DiskResidentEventList.cpp.

                                                                                                  {
  // mainly for diagnostics
  DalitzHistoSet hs;
  if(0 == w) return hs;
  for(unsigned int i=0; i< this->size(); i++){
    DalitzEvent evt( this->getEvent(i));
    hs.addEvent(evt, w->RealVal(evt));
  }
  return hs;
}

save()

bool DiskResidentEventList::save

(

)

Definition at line 375 of file DiskResidentEventList.cpp.

                                {
  bool dbThis=false;
  if(0 == _ntp) return false;
  if(0 == _f) return false;
  bool success=true;
  if(dbThis){
    cout << "DiskResidentEventList::save() called"
       << " for filename = " << _fname << endl;
  }
  //_f->cd();
  // success &=
  //_f = _ntp->GetCurrentFile();
  //_f->cd();
  if(0 == _f || (! _f->IsWritable()) || (! _f->IsOpen())){
    cout << "big problem with f in DiskResidentEventList::save()" << endl;
    throw "Scheibenkartoffel";
  }
  _ntp->AutoSave();
  if(dbThis)cout << " done the saving" << endl;
  return success;
}

size()

unsigned int DiskResidentEventList::size ( ) const

virtual

Implements MINT::IMinimalEventList< DalitzEvent >.

Definition at line 287 of file DiskResidentEventList.cpp.

                                              {
  if(0 == _ntp) return 0;
  return _ntp->GetEntries()/_scaleData;
}

weighedReWeightedHistoSet()

DalitzHistoSet DiskResidentEventList::weighedReWeightedHistoSet

(

MINT::IReturnRealForEvent< IDalitzEvent > *

w

)

const

Definition at line 491 of file DiskResidentEventList.cpp.

                                                                                                         {
  // mainly for diagnostics
  DalitzHistoSet hs;
  if(0 == w) return hs;
  for(unsigned int i=0; i< this->size(); i++){
    DalitzEvent evt(this->getEvent(i));
    hs.addEvent(evt, w->RealVal(evt) * (evt.getWeight()));
  }
  return hs;
}

weightedHistoSet()

DalitzHistoSet DiskResidentEventList::weightedHistoSet

(

)

const

Definition at line 471 of file DiskResidentEventList.cpp.

                                                            {
  // mainly for diagnostics
  DalitzHistoSet hs;
  for(unsigned int i=0; i< this->size(); i++){
    DalitzEvent evt(this->getEvent(i));
    hs.addEvent(evt, evt.getWeight());
  }
  return hs;
}

Member Data Documentation

__maxBytes

int DiskResidentEventList::__maxBytes = 10000000

staticprotected

Definition at line 38 of file DiskResidentEventList.h.

_cName

std::string DiskResidentEventList::_cName

private

Definition at line 111 of file DiskResidentEventList.h.

_f

TFile* DiskResidentEventList::_f

protected

Definition at line 31 of file DiskResidentEventList.h.

_fname

std::string DiskResidentEventList::_fname

protected

Definition at line 28 of file DiskResidentEventList.h.

_initialised

bool DiskResidentEventList::_initialised

protected

Definition at line 36 of file DiskResidentEventList.h.

_ntp

TNtupleD* DiskResidentEventList::_ntp

protected

Definition at line 32 of file DiskResidentEventList.h.

_ntpName

std::string DiskResidentEventList::_ntpName

private

Definition at line 112 of file DiskResidentEventList.h.

_opt

std::string DiskResidentEventList::_opt

protected

Definition at line 29 of file DiskResidentEventList.h.

_scaleData

double DiskResidentEventList::_scaleData

protected

Definition at line 33 of file DiskResidentEventList.h.

---

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

• Mint/DiskResidentEventList.h
• src/Mojito/DalitzEvents/DiskResidentEventList.cpp