#include #include #include #include #include #include #include "MLog.h" #include "MLogManip.h" #include "MEnv.h" #include "MArgs.h" #include "MArray.h" #include "MDirIter.h" #include "MStatusDisplay.h" #include "MDataSet.h" #include "MJSpectrum.h" using namespace std; static void StartUpMessage() { gLog << all << endl; // 1 2 3 4 5 // 12345678901234567890123456789012345678901234567890123456 gLog << "========================================================" << endl; gLog << " Sponde - MARS V" << MARSVER << endl; gLog << " MARS -- SPectrum ON DEmand" << endl; gLog << " Compiled with ROOT v" << ROOT_RELEASE << " on <" << __DATE__ << ">" << 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 << " sponde [options] inputfile.root mcdataset.txt [outputfile.root]" << endl << endl; gLog << " Arguments:" << endl; gLog << " inputfile.root: Ganymed output [and result] file" << endl; gLog << " mcdataset.txt: Dataset describing the Monte Carlos to be used" << endl; gLog << " For more details see MDataSet." << endl; gLog << " outputfile.root: Optional file in which the result is stored" << endl << endl; gLog << " Root Options:" << endl; gLog << " -b Batch mode (no graphical output to screen)" << endl<" << endl; gLog << " --debug-env[=1] Display untouched resources after program execution" << endl; gLog << " --debug-env=2 Display untouched resources after eventloop setup" << endl; gLog << " --debug-env=3 Debug setting resources from resource file" << endl; gLog << " --debug-mem Debug memory usage" << endl << endl; gLog << endl; gLog << " -q Quit when job is finished" << endl; gLog << " -f Force overwrite of existing files" << endl; gLog << " --print-ds Print Dataset information" << endl; gLog << " --print-files Print Files taken from Sequences ('+' found, '-' missing)" << endl; gLog << " --config=sponde.rc Resource file [default=sponde.rc]" << endl; gLog << endl; gLog << " --version, -V Show startup message with version number" << endl; gLog << " -?, -h, --help This help" << endl << endl; gLog << "Background:" << endl; gLog << " Sponde is a natural satellite of Jupiter. It was discovered by a team" << endl; gLog << " of astronomers from the University of Hawaii led by Scott S. Sheppard," << endl; gLog << " et al in 2001, and given the temporary designation S/2001 J 5. Sponde" << endl; gLog << " is about 2 kilometers in diameter, and orbits Jupiter at an average" << endl; gLog << " of 23,487,000 kilometers. It is also designated as Jupiter XXXVI." << endl; gLog << " It is named after one of the Horae (Hours), which presided over the" << endl; gLog << " seventh hour (libations poured after lunch). The Hours, godesses of" << endl; gLog << " the time of day but also of the seasons, were daughters of Zeus and" << endl; gLog << " Themis." << endl; gLog << " It belongs to the Pasipha‰ group, irregular retrograde moons orbiting" << endl; gLog << " Jupiter at distances ranging between 22.8 and 24.1 Gm, and with" << endl; gLog << " inclinations ranging between 144.5ø and 158.3ø." << endl << endl; } int main(int argc, char **argv) { if (!MARS::CheckRootVer()) return 0xff; // // Evaluate arguments // MArgs arg(argc, argv, kTRUE); gLog.Setup(arg); StartUpMessage(); if (arg.HasOnly("-V") || arg.HasOnly("--version")) return 0; if (arg.HasOnly("-?") || arg.HasOnly("-h") || arg.HasOnly("--help")) { Usage(); return 2; } const TString kConfig = arg.GetStringAndRemove("--config=", "sponde.rc"); const Bool_t kPrintSeq = arg.HasOnlyAndRemove("--print-ds"); const Bool_t kPrintFiles = arg.HasOnlyAndRemove("--print-files"); const Bool_t kDebugMem = arg.HasOnlyAndRemove("--debug-mem"); Int_t kDebugEnv = arg.HasOnlyAndRemove("--debug-env") ? 1 : 0; kDebugEnv = arg.GetIntAndRemove("--debug-env=", kDebugEnv); const Bool_t kQuit = arg.HasOnlyAndRemove("-q"); const Bool_t kBatch = arg.HasOnlyAndRemove("-b"); const Bool_t kOverwrite = arg.HasOnlyAndRemove("-f"); const Bool_t kRefill = arg.HasOnlyAndRemove("--refill"); const Bool_t kSimple = !arg.HasOnlyAndRemove("--accurate"); const Bool_t kRawMc = arg.HasOnlyAndRemove("--raw-mc"); if (arg.GetNumOptions()>0) { gLog << warn << "WARNING - Unknown command line options..." << endl; arg.Print("options"); gLog << endl; return 2; } // // check for the right usage of the program // if (arg.GetNumArguments()<2 || arg.GetNumArguments()>3) { Usage(); return 2; } // // Setup sequence file and check for its existence // const TString kInfile = arg.GetArgumentStr(0); const TString kDataset = arg.GetArgumentStr(1); const TString kOutfile = arg.GetArgumentStr(2); if (gSystem->AccessPathName(kInfile, kFileExists)) { gLog << err << "Sorry, sequences file '" << kInfile << "' doesn't exist." << endl; return 2; } if (gSystem->AccessPathName(kDataset, kFileExists)) { gLog << err << "Sorry, dataset file '" << kDataset << "' doesn't exist." << endl; return 2; } if (kDebugMem) TObject::SetObjectStat(kTRUE); // // Setup sequence and check its validity // MDataSet seq(kDataset); if (kPrintSeq || kPrintFiles) { gLog << all; gLog.Separator(kDataset); seq.Print(kPrintFiles?"files":""); gLog << endl; } if (!seq.IsValid()) { gLog << err << "Sequences read but not valid!" << endl << endl; return 2; } // // Initialize root // MArray::Class()->IgnoreTObjectStreamer(); MParContainer::Class()->IgnoreTObjectStreamer(); TApplication app("sponde", &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(kDataset); // // Calculate pedestal for pedestal-calculation and calibration // MEnv env(kConfig); MJSpectrum job(Form("MJSpectrum #%d", seq.GetNumAnalysis())); job.SetEnv(&env); job.SetEnvDebug(kDebugEnv); job.SetDisplay(d);; job.SetOverwrite(kOverwrite); job.SetPathOut(kOutfile); job.SetPathIn(kInfile); job.EnableRefilling(kRefill); job.EnableSimpleMode(kSimple); job.EnableRawMc(kRawMc); if (!job.Process(seq)) { gLog << err << "Calculation of spectrum failed." << endl << endl; return 2; } if (kDebugEnv>0) env.PrintUntouched(); 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; }