/*************************************************************************** musrview.cpp Author: Andreas Suter e-mail: andreas.suter@psi.ch $Id$ ***************************************************************************/ /*************************************************************************** * Copyright (C) 2007 by Andreas Suter * * andreas.suter@psi.ch * * * * This program is free software; you can redistribute it and/or modify * * it under the terms of the GNU General Public License as published by * * the Free Software Foundation; either version 2 of the License, or * * (at your option) any later version. * * * * This program is distributed in the hope that it will be useful, * * but WITHOUT ANY WARRANTY; without even the implied warranty of * * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * * GNU General Public License for more details. * * * * You should have received a copy of the GNU General Public License * * along with this program; if not, write to the * * Free Software Foundation, Inc., * * 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. * ***************************************************************************/ #include #include #include #include using namespace std; #include #include #include #include "PMusr.h" #include "PStartupHandler.h" #include "PMsrHandler.h" #include "PRunDataHandler.h" #include "PRunListCollection.h" #include "PMusrCanvas.h" //-------------------------------------------------------------------------- /** *

* */ void musrview_syntax() { cout << endl << "usage: musrview | --version | --help"; cout << endl << " : msr/mlog input file"; cout << endl << " 'musrview ' will execute msrfit"; cout << endl << " 'musrview' or 'msrfit --help' will show this help"; cout << endl << " 'musrview --version' will print the msrfit version"; cout << endl << endl; } //-------------------------------------------------------------------------- int main(int argc, char *argv[]) { bool show_syntax = false; int status; bool success = true; switch (argc) { case 1: show_syntax = true; break; case 2: if (strstr(argv[1], "--version")) { cout << endl << "musrview version: $Id$"; cout << endl << endl; return PMUSR_SUCCESS; } else if (strstr(argv[1], "--help")) { show_syntax = true; } else { // check if filename has extension msr or mlog if (!strstr(argv[1], ".msr") && !strstr(argv[1], ".mlog")) { cout << endl << "**ERROR** " << argv[1] << " is not a msr/mlog-file!" << endl; show_syntax = true; } } break; default: show_syntax = true; } if (show_syntax) { musrview_syntax(); return PMUSR_WRONG_STARTUP_SYNTAX; } // get default path (for the moment only linux like) char *pmusrpath; char musrpath[128]; pmusrpath = getenv("MUSRFITPATH"); if (pmusrpath == 0) { // not set, will try default one strcpy(musrpath, "/home/nemu/analysis/bin"); cout << endl << "**WARNING** MUSRFITPATH environment variable not set will try " << musrpath << endl; } else { strncpy(musrpath, pmusrpath, sizeof(musrpath)); } // read startup file char startup_path_name[128]; sprintf(startup_path_name, "%s/musrfit_startup.xml", musrpath); TSAXParser *saxParser = new TSAXParser(); PStartupHandler *startupHandler = new PStartupHandler(); saxParser->ConnectToHandler("PStartupHandler", startupHandler); status = saxParser->ParseFile(startup_path_name); // check for parse errors if (status) { // error cout << endl << "**WARNING** reading/parsing musrfit_startup.xml failed."; cout << endl; /* // clean up if (saxParser) { delete saxParser; saxParser = 0; } if (startupHandler) { delete startupHandler; startupHandler = 0; } */ } startupHandler->CheckLists(); // read msr-file PMsrHandler *msrHandler = new PMsrHandler(argv[1]); status = msrHandler->ReadMsrFile(); if (status != PMUSR_SUCCESS) { switch (status) { case PMUSR_MSR_FILE_NOT_FOUND: cout << endl << "**ERROR** couldn't find '" << argv[1] << "'" << endl << endl; break; case PMUSR_MSR_SYNTAX_ERROR: cout << endl << "**SYNTAX ERROR** in file " << argv[1] << ", full stop here." << endl << endl; default: cout << endl << "**UNKNOWN ERROR** when trying to read the msr-file" << endl << endl; break; } return status; } // make a plot list vector PMsrPlotList *msrPlotList = msrHandler->GetMsrPlotList(); PIntVector plotList; bool runPresent; for (unsigned int i=0; isize(); i++) { for (unsigned int j=0; jat(i).fRuns.size(); j++) { // check that run is not already in the plotList runPresent = false; for (unsigned int k=0; kat(i).fRuns[j].Re() == static_cast(plotList[k])) { runPresent = true; break; } } if (!runPresent) { plotList.push_back(static_cast(msrPlotList->at(i).fRuns[j].Re())); } } } /* cout << endl << ">> plotList:" << endl; for(unsigned int i=0; iGetDataPathList()); else dataHandler = new PRunDataHandler(msrHandler); success = dataHandler->IsAllDataAvailable(); if (!success) { cout << endl << "**ERROR** Couldn't read all data files, will quit ..." << endl; } // generate the necessary histogramms for the view PRunListCollection *runListCollection = 0; if (success) { // feed all the necessary histogramms for the view runListCollection = new PRunListCollection(msrHandler, dataHandler); for (unsigned int i=0; iGetMsrRunList()->size(); i++) { // if run is in plotList add it, otherwise go to the next runPresent = false; for (unsigned int j=0; j(plotList[j]) == i+1) { runPresent = true; break; } } if (runPresent) { //cout << endl << ">> Will add run " << i << endl; success = runListCollection->Add(i, kView); if (!success) { cout << endl << "**ERROR** Couldn't handle run no " << i << " "; cout << (*msrHandler->GetMsrRunList())[i].fRunName.Data(); break; } } } } if (success) { // generate Root application needed for PMusrCanvas TApplication app("App", &argc, argv); vector canvasVector; PMusrCanvas *musrCanvas; bool ok = true; for (unsigned int i=0; iGetMsrPlotList()->size(); i++) { if (startupHandler) musrCanvas = new PMusrCanvas(i, msrHandler->GetMsrTitle()->Data(), 10+i*100, 10+i*100, 800, 600, startupHandler->GetMarkerList(), startupHandler->GetColorList()); else musrCanvas = new PMusrCanvas(i, msrHandler->GetMsrTitle()->Data(), 10+i*100, 10+i*100, 800, 600); if (!musrCanvas->IsValid()) { cout << endl << "**SEVERE ERROR** Couldn't invoke all necessary objects, will quit."; cout << endl; ok = false; break; } // ugly but rootcint cannot handle the spirit-parser framework musrCanvas->SetMsrHandler(msrHandler); musrCanvas->SetRunListCollection(runListCollection); musrCanvas->UpdateParamTheoryPad(); musrCanvas->UpdateDataTheoryPad(); musrCanvas->UpdateInfoPad(); if (!musrCanvas->IsValid()) { // something went wrong ok = false; break; } musrCanvas->Connect("Done(Int_t)", "TApplication", &app, "Terminate(Int_t)"); // keep musrCanvas objects canvasVector.push_back(musrCanvas); } // check that everything is ok if (ok) app.Run(true); // true needed that Run will return after quit so that cleanup works // clean up cout << endl << "clean up canvas vector ..."; char canvasName[32]; for (unsigned int i=0; i> canvasName=" << canvasName; if (gROOT->GetListOfCanvases()->FindObject(canvasName) != 0) { cout << endl << ">> canvasName=" << canvasName << ", found ..."; cout << endl; canvasVector[i]->~PMusrCanvas(); } } canvasVector.empty(); } cout << endl; // clean up plotList.clear(); if (saxParser) { delete saxParser; saxParser = 0; } if (startupHandler) { delete startupHandler; startupHandler = 0; } if (msrHandler) { delete msrHandler; msrHandler = 0; } if (dataHandler) { delete dataHandler; dataHandler = 0; } if (runListCollection) { delete runListCollection; runListCollection = 0; } return 0; }