#include #include #include #include #include #include #include "MLog.h" #include "MLogManip.h" #include "MArgs.h" #include "MArray.h" #include "MDirIter.h" #include "MStatusDisplay.h" #include "MSequences.h" #include "MJCut.h" using namespace std; static void StartUpMessage() { gLog << all << endl; // 1 2 3 4 5 // 12345678901234567890123456789012345678901234567890 gLog << "========================================================" << endl; gLog << " Ganymed - MARS V" << MARSVER << endl; gLog << " MARS -- Gammas Are Now Your Most Exciting Discovery" << endl; gLog << " Compiled on <" << __DATE__ << ">" << endl; gLog << " Using ROOT v" << ROOTVER << endl; gLog << "========================================================" << endl; gLog << endl; } static void Usage() { // 1 2 3 4 5 6 7 8 // 12345678901234567890123456789012345678901234567890123456789012345678901234567890 gLog << all << endl; gLog << "Sorry the usage is:" << endl; gLog << " ganymed [-c] [-y] [options] sequences.txt" << endl << endl; gLog << " Arguments:" << endl; gLog << " sequences.txt: An ascii file defining a collection of sequences" << endl; gLog << " Root Options:" << endl; gLog << " -b Batch mode (no graphical output to screen)" << endl<0) { gLog << warn << "WARNING - Unknown commandline options..." << endl; arg.Print("options"); gLog << endl; } // // check for the right usage of the program // if (arg.GetNumArguments()!=1) { Usage(); return -1; } // // Setup sequence file and check for its existance // const TString kSequences = arg.GetArgumentStr(0); if (gSystem->AccessPathName(kSequences, kFileExists)) { gLog << err << "Sorry, sequences file '" << kSequences << "' doesn't exist." << endl; return -1; } if (kDebugMem) TObject::SetObjectStat(kTRUE); // // Setup sequence and check its validity // MSequences seq(kSequences); if (kNumAnalysis>=0) seq.SetNumAnalysis(kNumAnalysis); if (kPrintSeq || kPrintFiles) { gLog << all; gLog.Separator(kSequences); seq.Print(kPrintFiles?"files":""); gLog << endl; } if (!seq.IsValid()) { gLog << err << "Sequences read but not valid!" << endl << endl; return -1; } // // Initialize root // MArray::Class()->IgnoreTObjectStreamer(); MParContainer::Class()->IgnoreTObjectStreamer(); TApplication app("Ganymed", &argc, argv); if (!gROOT->IsBatch() && !gClient || gROOT->IsBatch() && !kBatch) { gLog << err << "Bombing... maybe your DISPLAY variable is not set correctly!" << endl; return 1; } // // Update frequency by default = 1Hz // MStatusDisplay *d = new MStatusDisplay; // From now on each 'Exit' means: Terminate the application d->SetBit(MStatusDisplay::kExitLoopOnExit); d->SetTitle(kSequences); // // Calculate pedestal for pedestal-calculation and calibration // MJCut job(Form("MJCut #%d", seq.GetNumAnalysis())); job.SetEnv(kConfig); job.SetEnvDebug(kDebugEnv); job.SetDisplay(d);; job.SetOverwrite(kOverwrite); job.SetPathOut(kOutpath); job.SetNameOutFile(kOutfile); job.SetNameSummaryFile(kNameSummary); job.SetNameResultFile(kNameResult); job.EnableWriteOnly(kWriteOnly); if (kWriteSummary) // Don't change flag set in SetNameSummaryFile job.EnableStorageOfSummary(); if (kWriteResult) // Don't change flag set in SetNameSummaryFile job.EnableStorageOfResult(); if (!job.ProcessFile(seq)) { gLog << err << "Calculation of cuts failed." << endl << endl; return -1; } if (!job.GetDisplay()) { gLog << warn << "Display closed by user... execution aborted." << endl << endl; return 1; } if (kBatch || kQuit) delete d; else { // From now on each 'Close' means: Terminate the application d->SetBit(MStatusDisplay::kExitLoopOnClose); // Wait until the user decides to exit the application app.Run(kFALSE); } if (TObject::GetObjectStat()) { TObject::SetObjectStat(kFALSE); gObjectTable->Print(); } return 0; }