genie::GSimFiles Class Reference

#include <GSimFiles.h>

Collaboration diagram for genie::GSimFiles:

Public Member Functions
	GSimFiles (bool chain=true, const int nmaxmodels=10)
	~GSimFiles (void)
int	NModels (void) const
int	FindModelID (string tag) const
string	ModelTag (int imodel) const
TFile *	XSecFile (int imodel) const
string	XSecFileName (int imodel) const
TChain *	EvtChain (int imodel) const
vector< string > &	EvtFileNames (int imodel) const
const string &	PathToXMLFile (void) const
void	Print (ostream &stream) const
bool	LoadFromFile (string xmlfile)

Private Member Functions
void	Init (const int nmaxmodels)
void	CleanUp (void)

Private Attributes
bool	fDoChain
int	fNModels
vector< string > *	fModelTag
vector< TFile * > *	fXSecFile
vector< string > *	fXSecFileName
vector< TChain * > *	fEvtChain
vector< vector< string > > *	fEvtFileNames
string	fPath2XMLFile

Friends
ostream &	operator<< (ostream &stream, const GSimFiles &gsimf)

Detailed Description

Definition at line 64 of file GSimFiles.h.

Constructor & Destructor Documentation

GSimFiles::GSimFiles	(	bool	chain = true,
		const int	nmaxmodels = 10
	)

Definition at line 34 of file GSimFiles.cxx.

References fDoChain, and Init().

{
  fDoChain = chain;
  this->Init(nmaxmodels);
}

GSimFiles::~GSimFiles

(

void

)

Definition at line 40 of file GSimFiles.cxx.

References CleanUp().

{
  this->CleanUp();
}

Member Function Documentation

void GSimFiles::CleanUp ( void  )

private

Definition at line 273 of file GSimFiles.cxx.

References fPath2XMLFile.

Referenced by ~GSimFiles().

{

   fPath2XMLFile = "";

}

TChain * GSimFiles::EvtChain

(

int

imodel

)

const

Definition at line 76 of file GSimFiles.cxx.

References fEvtChain.

{
  return (*fEvtChain)[imodel];
}

vector< string > & GSimFiles::EvtFileNames

(

int

imodel

)

const

Definition at line 81 of file GSimFiles.cxx.

References fEvtFileNames.

Referenced by Print().

{
  return (*fEvtFileNames)[imodel];
}

int GSimFiles::FindModelID

(

string

tag

)

const

Definition at line 50 of file GSimFiles.cxx.

References fModelTag.

{
  int imodel = 0;
  vector<string>::const_iterator it = fModelTag->begin();
  for( ; it != fModelTag->end(); ++it) {
    if(*it == tag) return imodel;
    imodel++;
  }
  return -1;
}

void GSimFiles::Init ( const int  nmaxmodels )

private

Definition at line 254 of file GSimFiles.cxx.

References fEvtChain, fEvtFileNames, fModelTag, fNModels, fPath2XMLFile, fXSecFile, and fXSecFileName.

Referenced by GSimFiles().

{
  fNModels      = 0;
  fModelTag     = new vector<string>          (nmaxmodels);
  fXSecFile     = new vector<TFile*>          (nmaxmodels);
  fXSecFileName = new vector<string>          (nmaxmodels);
  fEvtChain     = new vector<TChain*>         (nmaxmodels);
  fEvtFileNames = new vector<vector<string> > (nmaxmodels);

  for(int i=0; i<nmaxmodels; i++) {
    (*fModelTag) [i] = "";
    (*fXSecFile) [i] = 0;
    (*fEvtChain) [i] = 0;
  }

  fPath2XMLFile = "";

}

bool GSimFiles::LoadFromFile

(

string

xmlfile

)

Definition at line 91 of file GSimFiles.cxx.

References fDoChain, fEvtChain, fEvtFileNames, fModelTag, fNModels, fPath2XMLFile, fXSecFile, fXSecFileName, genie::gAbortingInErr, LOG, pDEBUG, pERROR, pFATAL, pNOTICE, and genie::utils::str::TrimSpaces().

{
  LOG("GSimFiles", pNOTICE) << "Loading: " << xmlfile;

  vector<string>  &          model_tag      = *fModelTag;
  vector<TFile*>  &          xsec_file      = *fXSecFile;
  vector<string>  &          xsec_filename  = *fXSecFileName;
  vector<TChain*> &          evt_chain      = *fEvtChain;
  vector< vector<string> > & evt_filenames  = *fEvtFileNames;

  xmlTextReaderPtr reader = xmlNewTextReaderFilename(xmlfile.c_str());
  if(reader == NULL) {
    return false;
  }

  const int kNodeTypeStartElement = 1;
  const int kNodeTypeEndElement   = 15;

  int  imodel            = 0;
  bool is_xsec_file      = false;
  bool is_evt_file       = false;
  bool is_ghep_evt_file  = false;
  bool is_gst_evt_file   = false;
  bool have_ghep_files   = false;
  bool have_gst_files    = false;

  if (reader != NULL) {
     int ret = xmlTextReaderRead(reader);
     while (ret == 1) {
       xmlChar * name  = xmlTextReaderName     (reader);
       xmlChar * value = xmlTextReaderValue    (reader);
       int       type  = xmlTextReaderNodeType (reader);
       int       depth = xmlTextReaderDepth    (reader);

       bool start_element = (type==kNodeTypeStartElement);
       bool end_element   = (type==kNodeTypeEndElement);

       if(depth==0 && start_element) {
         LOG("GSimFiles", pDEBUG) << "Root element = " << name;
         if(xmlStrcmp(name, (const xmlChar *) "genie_simulation_outputs")) {
           LOG("GSimFiles", pERROR)
             << "\nXML doc. has invalid root element! [filename: "
             << xmlfile << "]";
           return false;
         }
       }

       if( (!xmlStrcmp(name, (const xmlChar *) "model")) && start_element) {
         xmlChar * xname = xmlTextReaderGetAttribute(reader,(const xmlChar*)"name");
         string sname    = utils::str::TrimSpaces((const char *)xname);
         model_tag[imodel] = sname;
         LOG("GSimFiles", pNOTICE)
            << "Adding files for model ID: "
            << imodel << " (" << model_tag[imodel] << ")";
         xmlFree(xname);
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "model")) && end_element) {
         LOG("GSimFiles", pNOTICE)
            << "Done adding files for model ID: " << imodel;
         imodel++;
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "xsec_file")) && start_element) {
         is_xsec_file = true;
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "xsec_file")) && end_element) {
         is_xsec_file = false;
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "evt_file")) && start_element) {
         is_evt_file       = true;
         is_ghep_evt_file  = false;
         is_gst_evt_file   = false;
         xmlChar * xfmt = xmlTextReaderGetAttribute(reader,(const xmlChar*)"format");
         string sfmt = utils::str::TrimSpaces((const char *)xfmt);
         if (sfmt.find("gst")  != string::npos)
         {
            is_gst_evt_file  = true;
            if(!have_gst_files) { have_gst_files = true; }
         }
         else
         if (sfmt.find("ghep") != string::npos)
         {
            is_ghep_evt_file = true;
            if(!have_ghep_files) { have_ghep_files = true; }
         }
         if(have_gst_files && have_ghep_files) {
            LOG("GSimFiles", pFATAL)
               << "Oops! You shouldn't mix GHEP and GST event files in GSimFiles";
            LOG("GSimFiles", pFATAL)
               << "Please correct XML file: " << xmlfile;
            gAbortingInErr = true;;
            exit(1);
         }
         xmlFree(xfmt);
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "evt_file")) && end_element) {
         is_evt_file  = false;
       }
       if( (!xmlStrcmp(name, (const xmlChar *) "#text")) && depth==3) {
         string filename = utils::str::TrimSpaces((const char *)value);
         if(is_evt_file) {
           LOG("GSimFiles", pNOTICE)
                << " * Adding event file: " << filename;
           // chain the event trees, if requested
           if(fDoChain) {
             if(!evt_chain[imodel] && is_gst_evt_file) {
                 evt_chain[imodel] = new TChain("gst");
             } else
             if(!evt_chain[imodel] && is_ghep_evt_file) {
                 evt_chain[imodel] = new TChain("gtree");
             }
             if(evt_chain[imodel]) {
                evt_chain[imodel]->Add(filename.c_str());
             }
           }//chain?
           evt_filenames[imodel].push_back(filename);
         }
         if(is_xsec_file) {
           LOG("GSimFiles", pNOTICE)
                << " * Adding cross section file: " << filename;
           xsec_file    [imodel] = new TFile(filename.c_str(), "read");
           xsec_filename[imodel] = filename;
           if(!xsec_file[imodel]) {
             exit(1);
           }
         }
       }

       xmlFree(name);
       xmlFree(value);

       ret = xmlTextReaderRead(reader);

     }//ret==1

     xmlFreeTextReader(reader);

  }//reader!=null

  fNModels = imodel;

  fPath2XMLFile = xmlfile;

  return true;
}

string GSimFiles::ModelTag

(

int

imodel

)

const

Definition at line 61 of file GSimFiles.cxx.

References fModelTag.

Referenced by Print().

{
  return (*fModelTag)[imodel];
}

int GSimFiles::NModels

(

void

)

const

Definition at line 45 of file GSimFiles.cxx.

References fNModels.

Referenced by Print().

{
  return fNModels;
}

const string & GSimFiles::PathToXMLFile

(

void

)

const

Definition at line 86 of file GSimFiles.cxx.

References fPath2XMLFile.

{
   return fPath2XMLFile;
}

void GSimFiles::Print

(

ostream &

stream

)

const

Definition at line 236 of file GSimFiles.cxx.

References EvtFileNames(), fPath2XMLFile, ModelTag(), NModels(), XSecFile(), and XSecFileName().

Referenced by genie::operator<<().

{
  stream << endl;
  stream << "loaded from path: " << fPath2XMLFile << endl;
  for(int imodel=0; imodel < this->NModels(); imodel++) {
     stream << "model tag: [" << this->ModelTag(imodel) << "]" << endl;
     if(this->XSecFile(imodel)) {
        stream << "   xsec file  : " << this->XSecFileName(imodel) << endl;
     }
     const vector<string> & filenames = this->EvtFileNames(imodel);
     vector<string>::const_iterator iter = filenames.begin();
     for( ; iter != filenames.end(); ++iter) {
        string filename = *iter;
        stream << "   event file : " << filename << endl;
     }
  }
}

TFile * GSimFiles::XSecFile

(

int

imodel

)

const

Definition at line 66 of file GSimFiles.cxx.

References fXSecFile.

Referenced by Print().

{
  return (*fXSecFile)[imodel];
}

string GSimFiles::XSecFileName

(

int

imodel

)

const

Definition at line 71 of file GSimFiles.cxx.

References fXSecFileName.

Referenced by Print().

{
  return (*fXSecFileName)[imodel];
}

Friends And Related Function Documentation

ostream& operator<< ( ostream &  stream, const GSimFiles &  gsimf  )

friend

Definition at line 27 of file GSimFiles.cxx.

  {
    f.Print(stream);
    return stream;
  }

Member Data Documentation

bool genie::GSimFiles::fDoChain

private

Definition at line 88 of file GSimFiles.h.

Referenced by GSimFiles(), and LoadFromFile().

vector<TChain*>* genie::GSimFiles::fEvtChain

private

Definition at line 93 of file GSimFiles.h.

Referenced by EvtChain(), Init(), and LoadFromFile().

vector<vector<string> >* genie::GSimFiles::fEvtFileNames

private

Definition at line 94 of file GSimFiles.h.

Referenced by EvtFileNames(), Init(), and LoadFromFile().

vector<string>* genie::GSimFiles::fModelTag

private

Definition at line 90 of file GSimFiles.h.

Referenced by FindModelID(), Init(), LoadFromFile(), and ModelTag().

int genie::GSimFiles::fNModels

private

Definition at line 89 of file GSimFiles.h.

Referenced by Init(), LoadFromFile(), and NModels().

string genie::GSimFiles::fPath2XMLFile

private

Definition at line 95 of file GSimFiles.h.

Referenced by CleanUp(), Init(), LoadFromFile(), PathToXMLFile(), and Print().

vector<TFile*>* genie::GSimFiles::fXSecFile

private

Definition at line 91 of file GSimFiles.h.

Referenced by Init(), LoadFromFile(), and XSecFile().

vector<string>* genie::GSimFiles::fXSecFileName

private

Definition at line 92 of file GSimFiles.h.

Referenced by Init(), LoadFromFile(), and XSecFileName().

---

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

• GSimFiles.h
• GSimFiles.cxx