#include #include #include #include #include #include "MArray.h" #include "MLog.h" #include "MLogManip.h" #include "MStatusDisplay.h" #include "MEnv.h" #include "MArgs.h" #include "MJSimulation.h" using namespace std; static void StartUpMessage() { gLog << all << endl; // 1 2 3 4 5 // 12345678901234567890123456789012345678901234567890 gLog << "====================================================" << endl; gLog << " Ceres - MARS V" << MARSVER << endl; gLog << " MARS - Camera Electronics and REflector Simulation" << endl; gLog << " Compiled with ROOT v" << ROOT_RELEASE << " on <" << __DATE__ << ">" << endl; gLog << "====================================================" << endl; gLog << endl; } static void Usage() { gLog << all << endl; gLog << "Sorry the usage is:" << endl; gLog << " ceres [options] inputfile" << endl << endl; gLog << " inputfile MMCS (CORSIKA) binary (cherenkov) file." << 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 and" << endl; gLog << " command line" << endl; gLog << " --debug-mem Debug memory usage" << endl << endl; gLog << " --rc=Name:option Set or overwrite a resource of the resource file." << endl; gLog << " (Note, that this option can be used multiple times)" << endl; gLog << endl; gLog << " Output options:" << endl; gLog << " -q Quit when job is finished" << endl; gLog << " -f Force overwrite of existing files" << endl; gLog << " -ff Force reading of file even if problems occur" << endl; gLog << " --out=path Path to write the all results to [def=local path]" << endl; gLog << " --dev-null Suppress output of files (for test purpose)" << endl; gLog << " --config=ceres.rc Resource file [default=reflector.rc]" << endl; gLog << endl; // gLog << " -f: force reading of runheader" << endl; gLog << " --version, -V Show startup message with version number" << endl; gLog << " -?, -h, --help This help" << endl << endl; gLog << "Background:" << endl; gLog << " Ceres, formal designation 1 Ceres, is the smallest identified dwarf planet in" << endl; gLog << " the Solar System and the only one in the asteroid belt. It was discovered on" << endl; gLog << " January 1, 1801, by Giuseppe Piazzi, and is named after the Roman goddess" << endl; gLog << " Ceres, the goddess of growing plants, the harvest, and of motherly love." << endl; gLog << " With a diameter of about 950km, Ceres is by far the largest and most massive" << endl; gLog << " body in the asteroid belt, and contains a third of the belt's total mass." << endl; gLog << " Recent observations have revealed that it is spherical, unlike the irregular" << endl; gLog << " shapes of smaller bodies with lower gravity. The surface of Ceres is probably" << endl; gLog << " made of a mixture of water ice and various hydrated minerals like carbonates" << endl; gLog << " and clays. Ceres appears to be differentiated into a rocky core and ice mantle." << endl; gLog << " It may harbour an ocean of liquid water underneath its surface, which makes it" << endl; gLog << " a potential target in the search for extraterrestrial life." << endl; gLog << " Ceres' apparent magnitude ranges from 6.7 to 9.3, hence at its brightest is" << endl; gLog << " still too dim to be seen with the naked eye. On Sept. 27, 2007, NASA launched" << endl; gLog << " the Dawn space probe to explore Vesta and Ceres." << endl << endl; gLog << "Example:" << endl; gLog << " ceres -f --out=outpath/ cer000001 cer000002" << endl; gLog << endl; } int main(int argc, char **argv) { if (!MARS::CheckRootVer()) return 0xff; MLog::RedirectErrorHandler(MLog::kColor); // // Evaluate arguments // MArgs arg(argc, argv); 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 Bool_t kBatch = arg.HasOnlyAndRemove("-b"); //const Int_t kCompLvl = arg.GetIntAndRemove("--comp=", 1); const Bool_t kForce = arg.HasOnlyAndRemove("-ff"); const Bool_t kDebugMem = arg.HasOnlyAndRemove("--debug-mem"); const Bool_t kNullOut = arg.HasOnlyAndRemove("--dev-null"); 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 kOverwrite = arg.HasOnlyAndRemove("-f"); const TString kConfig = arg.GetStringAndRemove("--config=", "ceres.rc"); const TString kOutpath = arg.GetStringAndRemove("--out=", "."); // // check for the right usage of the program (number of arguments) if (arg.GetNumArguments()<1/* || arg.GetNumArguments()>2*/) { gLog << warn << "WARNING - Wrong number of arguments..." << endl; Usage(); return 2; } // // Now we access/read the resource file. This will remove all // --rc= from the list of arguments. // MEnv env(kConfig, "ceres.rc"); if (!env.IsValid()) { gLog << err << "ERROR - Reading resource file " << kConfig << "." << endl; return 0xfe; } // And move the resource options from the command line to the MEnv if (!env.TakeEnv(arg, kDebugEnv>2)) return 0xfd; // // check for the right usage of the program (number of options) // if (arg.GetNumOptions()>0) { gLog << warn << "WARNING - Unknown commandline options..." << endl; arg.Print("options"); gLog << endl; return 2; } // // Initialize root // TVector2::Class()->IgnoreTObjectStreamer(); MArray::Class()->IgnoreTObjectStreamer(); MParContainer::Class()->IgnoreTObjectStreamer(); TApplication app("ceres", &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(Form("-- Ceres: %s --", ""/*kNamein.Data()*/)); if (kDebugMem) TObject::SetObjectStat(kTRUE); // // Do star in a block (debug mem) // { MJSimulation job("Ceres"/*Form("MJSimulation #%d", seq.GetSequence())*/); //job.SetSequence(seq); job.SetEnv(&env); job.SetEnvDebug(kDebugEnv); job.SetDisplay(d);; job.SetOverwrite(kOverwrite); job.SetPathOut(kOutpath); job.SetNullOut(kNullOut); job.SetForceMode(kForce); // job.SetPathIn(kInpath); // not yet needed if (!job.Process(arg)) { gLog << err << "Calculation of ceres 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; }