#include #include #include #include #include #include "MLog.h" #include "MLogManip.h" #include "MArgs.h" #include "MArray.h" #include "MDirIter.h" #include "MStatusDisplay.h" #include "MSequence.h" #include "MJPedestal.h" #include "MJCalibration.h" #include "MJCalibrateSignal.h" using namespace std; static void StartUpMessage() { gLog << all << endl; // 1 2 3 4 5 // 12345678901234567890123456789012345678901234567890 gLog << "========================================================" << endl; gLog << " Callisto - MARS V" << MARSVER << endl; gLog << " MARS -- CALibrate LIght Signals and Time Offsets" << 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 << " callisto [-c] [-y] [options] sequence.txt" << endl << endl; gLog << " Arguments:" << endl; gLog << " sequence.txt: An ascii file defining a sequence of runs" << 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 kSequence = arg.GetArgumentStr(0); if (gSystem->AccessPathName(kSequence, kFileExists)) { gLog << err << "Sorry, sequence file '" << kSequence << "' doesn't exist." << endl; return -1; } // // Setup sequence and check its validity // MSequence seq(kSequence); if (kPrintSeq) { gLog << all; gLog.Separator(kSequence); seq.Print(); gLog << endl; } if (!seq.IsValid()) { gLog << err << "Sequence read but not valid!" << endl << endl; return -1; } // // Process print options // if (kPrintFiles) { MDirIter Next1, Next2, Next3; seq.SetupPedRuns(Next1, kInpathD); seq.SetupCalRuns(Next2, kInpathD); seq.SetupDatRuns(Next3, kInpathD); gLog << all; gLog.Separator("Pedestal Files"); Next1.Print("all"); gLog << endl; gLog.Separator("Calibration Files"); Next2.Print("all"); gLog << endl; gLog.Separator("Data Files"); Next3.Print("all"); gLog << endl; } // // Check for existance of all files // MDirIter iter; const Int_t n0 = seq.SetupAllRuns(iter, kInpathD); const Int_t n1 = seq.GetNumAllRuns(); if (n0 != n1) { if (kForceExec) gLog << warn << "WARNING"; else gLog << err << "ERROR"; gLog << " - " << n1 << " files in sequence defined, but " << n0 << " found in "; gLog << (kInpathD.IsNull() ? "" : kInpathD.Data()) << endl; if (!kForceExec) return -1; } if (kPrintOnly) return 0; // // Initialize root // MArray::Class()->IgnoreTObjectStreamer(); MParContainer::Class()->IgnoreTObjectStreamer(); TApplication app("Callisto", &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(kSequence); if (kModeC) { // // Calculate pedestal for calibration // MJPedestal job1(Form("MJPedestalC #%d", seq.GetSequence())); job1.SetSequence(seq); job1.SetEnv(kConfig); job1.SetEnvDebug(kDebugEnv); job1.SetDisplay(d);; job1.SetOverwrite(kOverwrite); job1.SetPathData(kInpathD); // job1.SetPathOut(kOutpathC); // not yet needed // job1.SetPathIn(kInpathC); // not yet needed if (!job1.ProcessFile()) { gLog << err << "Calculation of pedestal failed." << endl << endl; return -1; } if (!job1.GetDisplay()) { gLog << warn << "Display closed by user... execution aborted." << endl << endl; return 1; } // // Do calibration // MJCalibration job2(Form("MJCalibration #%d", seq.GetSequence())); job2.SetSequence(seq); job2.SetEnv(kConfig); job2.SetEnvDebug(kDebugEnv); job2.SetDisplay(d);; job2.SetBadPixels(job1.GetBadPixels()); job2.SetOverwrite(kOverwrite); job2.SetPathOut(kOutpathC); job2.SetPathData(kInpathD); // job2.SetPathIn(kInpathC); // not yet needed if (!job2.ProcessFile(job1.GetPedestalCam())) { gLog << err << "Calculation of calibration failed." << endl << endl; return -1; } if (!job2.GetDisplay()) { gLog << warn << "Display closed by user... execution aborted." << endl << endl; return 1; } } if (kModeY) { d->Reset(); // // Calculate starting pedestal for data extraction // MJPedestal job1(Form("MJPedestalY #%d", seq.GetSequence())); job1.SetSequence(seq); job1.SetEnv(kConfig); job1.SetEnvDebug(kDebugEnv); job1.SetDisplay(d);; //job1.SetPathIn(kInpathY); // not yet needed //job1.SetPathOut(kOutpathY); // not yet needed job1.SetPathData(kInpathD); job1.SetOverwrite(kOverwrite); if (!job1.ProcessFile()) { gLog << err << "Calculation of pedestal failed." << endl << endl; return -1; } if (!job1.GetDisplay()) { gLog << warn << "Display closed by user... execution aborted." << endl << endl; return 1; } // // Extract signal and calibrate it // MJCalibrateSignal job2(Form("MJCalibrateSignal #%d", seq.GetSequence())); job2.SetSequence(seq); job2.SetDisplay(d);; job2.SetEnv(kConfig); job2.SetEnvDebug(kDebugEnv); job2.SetOverwrite(kOverwrite); job2.SetPathIn(kInpathY); job2.SetPathOut(kOutpathY); job2.SetPathData(kInpathD); // Where to search for calibration files if (!job2.ProcessFile(job1.GetPedestalCam(), kIsMC)) return -1; if (!job2.GetDisplay()) { gLog << warn << "Display closed by user... execution aborted." << endl << endl; return 1; } } if (kBatch || kQuit) { delete d; return 0; } // 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); return 0; }