#include #include #include #include "MParList.h" #include "MTaskList.h" #include "MEvtLoop.h" #include "MRawFileRead.h" #include "MSqlInsertRun.h" #include "MRawFileWrite.h" #include "MReportFileReadCC.h" #include "MWriteRootFile.h" #include "MLog.h" #include "MLogManip.h" #include "MArgs.h" #include "MTime.h" #include "MArray.h" #include "MRawEvtData.h" #include "MRawRunHeader.h" #include "MRawEvtHeader.h" #include "MRawCrateArray.h" #include "MFDataMember.h" using namespace std; ////////////////////////////////////////////////////////////////////////////// // // // This is an easy implementation of the Merging process // // (as compilable prog) // // // // at the moment it reads a binary file ("rawtest.bin") which was written // // in the DAQ raw format. // // // // The data are stored in root container objects (classes derived from // // TObject like MRawRunHeader) // // // // This containers are written to a root file ("rawtest.root") // // // ////////////////////////////////////////////////////////////////////////////// static void StartUpMessage() { gLog << all << endl; // 1 2 3 4 5 // 12345678901234567890123456789012345678901234567890 gLog << "==================================================" << endl; gLog << " MERPP - MARS V" << MARSVER << endl; gLog << " MARS - Merging and Preprocessing Program" << 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 << " merpp [options] inputfile[.rep,[.raw],[.txt]] [outputfile[.root]]" << endl << endl; gLog << " Arguments:" << endl; gLog << " inputfile.raw Magic DAQ binary file." << endl; gLog << " inputfile.rep Magic Central Control report file." << endl; gLog << " inputfile.txt Magic DC currents file." << endl; gLog << " ouputfile.root Merpped root file." << endl << endl; gLog << " Options:" << endl; gLog.Usage(); gLog << " -?, -h, --help This help" << endl << endl; gLog << " File Options:" << endl; gLog << " -c# Compression level #=1..9 [default=2]" << endl; gLog << " -f Force overwrite of an existing file" << endl; gLog << " -u, --update Update an existing file." << endl << endl; gLog << " Raw Data Options:" << endl; gLog << " -ff Force merpp to ignore broken events" << endl; gLog << " --interleave=# Process only each i-th event [default=1]" << endl << endl; // gLog << " --sql=mysql://user:password@url Insert run into database" << endl << endl; gLog << " Report File Options:" << endl; gLog << " --auto-time Take time automatically from MRawRunHeader" << endl; gLog << " (overwrites --start= and/or --stop=)" << endl; gLog << " --start=date/time Start event time" << endl; gLog << " --stop=date/time Stop event time" << endl; gLog << " --run=# Only data corresponding to this run number" << endl; gLog << " (from RUN-REPORT)" << endl; gLog << " --runfile=# Allow only run-control files" << endl; gLog << " (from .rep header)" << endl; gLog << " --sumfile Check for an all night summary file" << endl; gLog << " (from .rep header)" << endl; gLog << " --allfiles Don't check file type " << endl << endl; gLog << " REMARK: - At the moment you can process a .raw _or_ a .rep file, only!" << endl; gLog << " - 'date/time' has the format 'yyyy-mm-dd/hh:mm:ss.mmm'" << endl << endl; } // FIXME: Move to MTime (maybe 'InterpreteCmdline') MTime AnalyseTime(TString str) { Int_t y=0, ms=0, mon=0, d=0, h=0, m=0, s=0; const Int_t n = sscanf(str.Data(), "%d-%d-%d/%d:%d:%d.%d", &y, &mon, &d, &h, &m, &s, &ms); if (n<6 || n>7) { gLog << warn << "'" << str << "' no valid Time... ignored." << endl; return MTime(); } MTime t; t.Set(y, mon, d, h, m, s, ms); return t; } void GetTimeFromFile(const char *fname, MTime &start, MTime &stop) { TFile f(fname, "READ"); TTree *t = (TTree*)f.Get("RunHeaders"); if (t->GetEntries()!=1) { gLog << warn << "WARNING - File " << fname << " contains no or more than one entry in RunHeaders... Times unchanged." << endl; return; } MRawRunHeader *h = 0; t->SetBranchAddress("MRawRunHeader.", &h); t->GetEntry(0); if (!h) { gLog << warn << "WARNING - File " << fname << " did not contain RunHeaders.MRawRunHeader... Times unchanged." << endl; return; } if (!start) start = h->GetRunStart(); if (!stop) stop = h->GetRunEnd(); } int main(const int argc, char **argv) { StartUpMessage(); // // Evaluate arguments // MArgs arg(argc, argv); if (arg.HasOnly("-?") || arg.HasOnly("-h") || arg.HasOnly("--help")) { Usage(); return -1; } gLog.Setup(arg); const Int_t kComprlvl = arg.GetIntAndRemove("-c", 2); const Bool_t kInterleave = arg.GetIntAndRemove("--interleave=", 1); const Bool_t kForce = arg.HasOnlyAndRemove("-f"); const Bool_t kForceProc = arg.HasOnlyAndRemove("-ff"); const Int_t kRunNumber = arg.GetIntAndRemove("--run=", -1); const Bool_t kAutoTime = arg.HasOnlyAndRemove("--auto-time"); Int_t kRunFile = arg.GetIntAndRemove("--runfile=", -1); Bool_t kUpdate = arg.HasOnlyAndRemove("--update") || arg.HasOnlyAndRemove("-u"); MTime kTimeStart; MTime kTimeStop; if (arg.HasOption("--star=")) kTimeStart = AnalyseTime(arg.GetStringAndRemove("--start=")); if (arg.HasOption("--stop=")) kTimeStop = AnalyseTime(arg.GetStringAndRemove("--stop=")); // const TString kSqlDataBase(arg.GetStringAndRemove("--sql=")); if (arg.HasOnlyAndRemove("--sumfile")) kRunFile = 0; if (arg.GetNumOptions()>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 || arg.GetNumArguments()>2) { Usage(); return -1; } // // This is to make argv[i] more readable insidethe code // TString kNamein = arg.GetArgumentStr(0); TString kNameout = arg.GetArgumentStr(1); const Bool_t isreport = kNamein.EndsWith(".rep"); const Bool_t isdc = kNamein.EndsWith(".txt"); const Bool_t israw = !isreport && !isdc; if (!kNamein.EndsWith(".raw") && israw) kNamein += ".raw"; if (kNameout.IsNull()) kNameout = kNamein(0, kNamein.Last('.')); if (!kNameout.EndsWith(".root")) kNameout += ".root"; // if (!kSqlDataBase.IsNull() && !israw) // gLog << warn << "WARNING - Option '--sql=' only valid for raw-files... ignored." << endl; // // Initialize Non-GUI (batch) mode // gROOT->SetBatch(); // // check whether the given files are OK. // if (gSystem->AccessPathName(kNamein, kFileExists)) { gLog << err << "Sorry, the input file '" << kNamein << "' doesn't exist." << endl; return -1; } const Bool_t fileexist = !gSystem->AccessPathName(kNameout, kFileExists); const Bool_t writeperm = !gSystem->AccessPathName(kNameout, kWritePermission); if (fileexist && !writeperm) { gLog << err << "Sorry, you don't have write permission for '" << kNameout << "'." << endl; return -1; } if (fileexist && !kUpdate && !kForce) { gLog << err << "Sorry, file '" << kNameout << "' already existing." << endl; return -1; } if (!fileexist && kUpdate) { gLog << warn << "File '" << kNameout << "' doesn't yet exist." << endl; kUpdate=kFALSE; } // // Evaluate possible start-/stop-time // if (kAutoTime && kUpdate && (isreport || isdc)) GetTimeFromFile(kNameout, kTimeStart, kTimeStop); if (kTimeStart) gLog << inf << "Start Time: " << kTimeStart << endl; if (kTimeStop) gLog << inf << "Stop Time: " << kTimeStop << endl; // // Ignore TObject Streamer (bits, uniqueid) for MArray and MParContainer // MArray::Class()->IgnoreTObjectStreamer(); MParContainer::Class()->IgnoreTObjectStreamer(); // // create a (empty) list of parameters which can be used by the tasks // and an (empty) list of tasks which should be executed // MParList plist; MTaskList tasks; tasks.SetOwner(); plist.AddToList(&tasks); // // ---- The following is only necessary to supress some output ---- /* MRawRunHeader runheader; plist.AddToList(&runheader); MRawEvtHeader evtheader; plist.AddToList(&evtheader); MRawEvtData evtdata; plist.AddToList(&evtdata); MRawCrateArray cratearray; plist.AddToList(&cratearray); MTime evttime; plist.AddToList(&evttime); */ // // create the tasks which should be executed and add them to the list // in the case you don't need parameter containers, all of them can // be created by MRawFileRead::PreProcess // MTask *read = 0; MFilter *filter = 0; MTask *write = 0; const TString option(kUpdate ? "UPDATE" : "RECREATE"); if (isreport || isdc) { MWriteRootFile *w = new MWriteRootFile(kNameout, option, "Magic root-file", kComprlvl); if (isdc) { w->AddContainer("MTimeCurrents", "Currents"); w->AddContainer("MCameraDC", "Currents"); w->AddContainer("MReportCurrents", "Currents"); } else { w->AddContainer("MReportCamera", "Camera"); w->AddContainer("MTimeCamera", "Camera"); w->AddContainer("MCameraAUX", "Camera"); w->AddContainer("MCameraCalibration", "Camera"); w->AddContainer("MCameraCooling", "Camera"); w->AddContainer("MCameraHV", "Camera"); w->AddContainer("MCameraLV", "Camera"); w->AddContainer("MCameraLids", "Camera"); w->AddContainer("MReportTrigger", "Trigger"); w->AddContainer("MTimeTrigger", "Trigger"); w->AddContainer("MTriggerBit", "Trigger"); w->AddContainer("MTriggerIPR", "Trigger"); w->AddContainer("MTriggerCell", "Trigger"); w->AddContainer("MTriggerPrescFact", "Trigger"); w->AddContainer("MTriggerLiveTime", "Trigger"); w->AddContainer("MReportDrive", "Drive"); w->AddContainer("MTimeDrive", "Drive"); w->AddContainer("MReportCC", "CC"); w->AddContainer("MTimeCC", "CC"); w->AddContainer("MReportStarguider", "Starguider"); w->AddContainer("MTimeStarguider", "Starguider"); // w->AddContainer("MReportDAQ", "DAQ"); // w->AddContainer("MTimeDAQ", "DAQ"); } write = w; MReportFileReadCC *r = new MReportFileReadCC(kNamein); r->SetTimeStart(kTimeStart); r->SetTimeStop(kTimeStop); if (isdc) { r->SetHasNoHeader(); r->AddToList("MReportCurrents"); } else { r->SetRunNumber(kRunFile); r->AddToList("MReportCC"); //r->AddToList("MReportDAQ"); r->AddToList("MReportDrive"); r->AddToList("MReportCamera"); r->AddToList("MReportTrigger"); r->AddToList("MReportStarguider"); if (kRunNumber>0) { r->AddToList("MReportRun"); filter = new MFDataMember("MReportRun.fRunNumber", '=', kRunNumber); w->SetFilter(filter); } } read = r; } else { read = new MRawFileRead(kNamein); static_cast(read)->SetInterleave(kInterleave); static_cast(read)->SetForceMode(kForceProc); write = new MRawFileWrite(kNameout, option, "Magic root-file", kComprlvl); } tasks.AddToList(read); if (filter) tasks.AddToList(filter); /* if (israw && !kSqlDataBase.IsNull()) { MSqlInsertRun *ins = new MSqlInsertRun(kSqlDataBase); ins->SetUpdate(); tasks.AddToList(ins); }*/ tasks.AddToList(write); // // create the looping object and tell it about the parameters to use // and the tasks to execute // MEvtLoop magic; magic.SetParList(&plist); // // Start the eventloop which reads the raw file (MRawFileRead) and // write all the information into a root file (MRawFileWrite) // // between reading and writing we can do, transformations, checks, etc. // (I'm think of a task like MRawDataCheck) // if (!magic.Eventloop()) { gLog << err << "ERROR: Merging and preprocessing failed!" << endl; return -1; } tasks.PrintStatistics(); gLog << all << "Merpp finished successfull!" << endl; return 0; }