/* ======================================================================== *\ ! ! * ! * This file is part of MARS, the MAGIC Analysis and Reconstruction ! * Software. It is distributed to you in the hope that it can be a useful ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes. ! * It is distributed WITHOUT ANY WARRANTY. ! * ! * Permission to use, copy, modify and distribute this software and its ! * documentation for any purpose is hereby granted without fee, ! * provided that the above copyright notice appear in all copies and ! * that both that copyright notice and this permission notice appear ! * in supporting documentation. It is provided "as is" without express ! * or implied warranty. ! * ! ! ! Author(s): Thomas Bretz 12/2000 ! ! Copyright: MAGIC Software Development, 2000-2001 ! ! \* ======================================================================== */ ////////////////////////////////////////////////////////////////////////////// // // // MEvtLoop // // // // This class is the core of each event processing. // // First you must set the parameter list to use. The parameter list // // must contain the task list (MTaskList) to use. The name of the task // // list can be specified if you call Eventloop. The standard name is // // "MTaskList". The name you specify must match the name of the MTaskList // // object. // // // // If you call Eventloop first all PreProcess functions - with the // // parameter list as an argument - of the tasks in the task list are // // executed. If one of them returns kFALSE then the execution is stopped. // // If the preprocessing was ok. The Process funtion of the tasks are // // as long as one function returns kSTOP. Only the tasks which are marked // // marked as "All" or with a string which matches the MInputStreamID of // // MTaskList are executed. If one tasks returns kCONTINUE the pending // // tasks in the list are skipped and the execution in continued with // // the first one in the list. // // Afterwards the PostProcess functions are executed. // // // // // // Maybe we can add a TProgressMeter sometimes later to be able to show // // the progress graphically... // // // // // ////////////////////////////////////////////////////////////////////////////// #include "MEvtLoop.h" #include #include #include #include #include "MLog.h" #include "MLogManip.h" #include "MParList.h" #include "MTaskList.h" ClassImp(MEvtLoop); //! //! Maybe we can add a static parameter list to MEvtLoop //! Also we can derive MEvtLoop from MTaskList to have a static tasklist, too //! // -------------------------------------------------------------------------- // // default constructor - emty // MEvtLoop::MEvtLoop() : fParList(NULL), fProgress(NULL) { } // -------------------------------------------------------------------------- // // default destructor - emty // MEvtLoop::~MEvtLoop() { if (TestBit(kIsOwner) && fParList) delete fParList; } // -------------------------------------------------------------------------- // // if you set the Eventloop as owner the destructor of the given parameter // list is calles by the destructor of MEvtLoop, otherwise not. // void MEvtLoop::SetOwner(Bool_t enable) { enable ? SetBit(kIsOwner) : ResetBit(kIsOwner); } // -------------------------------------------------------------------------- // // The proprocessing part of the eventloop. Be careful, this is // for developers or use in special jobs only! // Bool_t MEvtLoop::PreProcess(const char *tlist) { // // check if the needed parameter list is set. // if (!fParList) { *fLog << err << dbginf << "Parlist not initialized." << endl; return kFALSE; } // // check for the existance of the specified task list // the default name is "MTaskList" // fTaskList = (MTaskList*)fParList->FindObject(tlist, "MTaskList"); if (!fTaskList) { *fLog << err << dbginf << "Cannot find tasklist '" << tlist << "' in parameter list." << endl; return kFALSE; } if (fLog != &gLog) fParList->SetLogStream(fLog); // // execute the preprocess of all tasks // connect the different tasks with the right containers in // the parameter list // if (!fTaskList->PreProcess(fParList)) { *fLog << err << "Error detected while PreProcessing" << endl; return kFALSE; } *fLog << endl; return kTRUE; } // -------------------------------------------------------------------------- // // The processing part of the eventloop. Be careful, this is // for developers or use in special jobs only! // void MEvtLoop::Process(Int_t maxcnt) const { // // loop over all events and process all tasks for // each event // *fLog << all <<"Eventloop running ("; if (maxcnt<0) *fLog << "all"; else *fLog << dec << maxcnt; *fLog << " events)..." << flush; if (fProgress && maxcnt>0) fProgress->SetRange(0, maxcnt); Int_t dummy = maxcnt<0 ? 0 : maxcnt; // // start a stopwatch // TStopwatch clock; clock.Start(); // // This is the MAIN EVENTLOOP which processes the data // if maxcnt<0 the number of processed events is counted // else only maxcnt events are processed // if (maxcnt<0) // process first and increment if sucessfull if (fProgress) while (fTaskList->Process()) { fProgress->SetPosition(++dummy); #if ROOT_VERSION_CODE < ROOT_VERSION(3,02,06) gSystem->ProcessEvents(); #else gClient->ProcessEventsFor(fProgress); #endif } else while (fTaskList->Process()) dummy++; else // check for number and break if unsuccessfull if (fProgress) while (dummy-- && fTaskList->Process()) { fProgress->SetPosition(maxcnt - dummy); #if ROOT_VERSION_CODE < ROOT_VERSION(3,02,06) gSystem->ProcessEvents(); #else gClient->ProcessEventsFor(fProgress); #endif } else while (dummy-- && fTaskList->Process()); // // stop stop-watch, print results // clock.Stop(); *fLog << all << "Ready!" << endl << endl; *fLog << dec << endl << "CPU - " << "Time: " << clock.CpuTime() << "s" << " for " << (maxcnt<0?dummy:maxcnt) << " Events" << " --> " << (maxcnt<0?dummy:maxcnt)/clock.CpuTime() << " Events/s" << endl; *fLog << "Real - " << "Time: " << clock.RealTime() << "s" << " for " << (maxcnt<0?dummy:maxcnt) << " Events" << " --> " << (maxcnt<0?dummy:maxcnt)/clock.RealTime() << " Events/s" << endl << endl; } // -------------------------------------------------------------------------- // // The postprocessing part of the eventloop. Be careful, this is // for developers or use in special jobs only! // Bool_t MEvtLoop::PostProcess() const { // // execute the post process of all tasks // return fTaskList->PostProcess(); } // -------------------------------------------------------------------------- // // See class description above. // Bool_t MEvtLoop::Eventloop(Int_t maxcnt, const char *tlist) { Bool_t rc = PreProcess(); // // If all Tasks were PreProcesses successfully start Processing. // if (rc) Process(maxcnt); // // Now postprocess all tasks. Only successfully preprocessed tasks are // postprocessed. If the Postprocessing of one task fail return an error. // if (!PostProcess()) return kFALSE; // // If postprocessing of all preprocessed tasks was sucefully return rc. // This gives an error in case the preprocessing has failed already. // Otherwise the eventloop is considered: successfully. // return rc; }