Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Data Structures | File List | Namespace Members | Data Fields | Globals | Related Pages

AnalysisAlgorithm.cpp

00001 #include "AnalysisAlgorithm.h"
00002 #include "AnalysisService.h"
00003 #include "SctData/TestResult.h"
00004 #include "SctData/FitScanResult.h"
00005 #include "SctData/RawScanResult.h"
00006 #include "Sct/UniqueID.h"
00007 #include "SctData/ConfigurationVariable.h"
00008 #include "Sct/SctNames.h"
00009 #include "Sct/LogicErrors.h"
00010 #include "Sct/IOName.h"
00011 #include "Sct/IS/IOManagerIS.h"
00012 #include "Sct/ISProxy/IOManagerISProxy.h"
00013 #include "Sct/ISUtilities.h"
00014 #include "SummaryWriter/SummaryManager.h"
00015 #include "SummaryWriter/SummaryWriter.h"
00016 #include "CalibrationController/IS/TestData.h"
00017 #include "SummaryWriter/TestSummaryIS.h"
00018 #include "ISSummaryWriter/ISSummaryManager.h"
00019 
00020 #include <boost/date_time/posix_time/posix_time_types.hpp>
00021 #include <sstream>
00022 #include <string>
00023 
00024 using namespace std;
00025 using namespace boost;
00026 using namespace Sct;
00027 using namespace SctData;
00028 using namespace Sct::IS;
00029 using namespace boost::posix_time;
00030 
00031 namespace SctAnalysis {
00032 
00033 AnalysisAlgorithm::AnalysisAlgorithm() : prototype(*this), ioTime(0), analysisTime(0), 
00034     overheadTime(0), nCalls(0), nAnalyses(0) {
00035 }
00036 
00037 AnalysisAlgorithm::AnalysisAlgorithm(shared_ptr<const TestData> testData, const string& moduleName, const AnalysisAlgorithm& prototype) throw()
00038     : m_modulename(moduleName), m_test(testData), m_done(false),  prototype(prototype),
00039     ioTime(0), analysisTime(0), overheadTime(0), nCalls(0), nAnalyses(0) {
00040   scanData.resize(getTestData()->nScans);
00041 }
00042 
00043 AnalysisAlgorithm::~AnalysisAlgorithm() throw() {
00044     boost::recursive_mutex::scoped_lock lock(m_access);
00045     //Need to explicitly delete so that lock is in force
00046     m_result.reset();
00047 }
00048 
00049 shared_ptr<const TestData> AnalysisAlgorithm::getTestData() const {
00050   if (!m_test.get()) throw Sct::InvariantViolatedError("Attempt to retrieve null TestData - probably from prototype algorithm",__FILE__,__LINE__);
00051   return m_test;
00052 }
00053 
00054 void AnalysisAlgorithm::checkTestResultExists() {
00055     if (m_result.get()==0){
00056     m_result = createTestResult();
00057     if (m_result.get()==0)
00058         throw IllegalStateError("AnalysisAlgorithm cannot initialize - no test exists", __FILE__, __LINE__);
00059     initializeTestResult();
00060     }
00061 }
00062 
00063 void AnalysisAlgorithm::initializeTestResult() {
00064     m_result->setRunNumber(getTestData()->runNumber);
00065     m_result->setModuleName(m_modulename);
00066     m_result->setTestVariable(*SctData::ConfigurationVariableIOHelper::getFromTypeRep(getTestData()->testVariable));
00067     m_result->getHeader().setStartTimeString(getTestData()->startTime);
00068 
00069     //Now add all the scans:
00070     for (unsigned int i=0; i<getTestData()->testPoints_size; ++i) {
00071         m_result->addScan(getTestData()->startScanNumber+i, getTestData()->testPoints[i]);
00072         //cout << "Analysis algorithm adding scan number ... "<<m_test->startScanNumber+i<<" to test"<<endl;
00073     }
00074 }
00075 
00076 void AnalysisAlgorithm::addFitScanResult(shared_ptr<Sct::IOName> name) {
00077     UniqueID id(name->getUniqueID());
00078     unsigned int index = id.getScanNumber() - getTestData()->startScanNumber;
00079     scanData[index].fitAvailable = true;
00080     scanData[index].fitName = name->getIOName();
00081 }
00082 
00083 void AnalysisAlgorithm::addRawScanResult(shared_ptr<Sct::IOName> name) {
00084     UniqueID id(name->getUniqueID());
00085     unsigned int index = id.getScanNumber() - getTestData()->startScanNumber;
00086     scanData[index].rawAvailable = true;
00087     scanData[index].rawName = name->getIOName();
00088 }
00089 
00090 void AnalysisAlgorithm::setFit(unsigned index, const std::string& fit){
00091     scanData[index].fitName=fit;
00092     scanData[index].fitAvailable=true;
00093 }
00094 
00095 void AnalysisAlgorithm::setRaw(unsigned index, const std::string& raw){
00096     scanData[index].rawName=raw;
00097     scanData[index].fitAvailable=true;
00098 }
00099 
00100 bool AnalysisAlgorithm::hasAllRaws() const {
00101     for (unsigned int i=0; i<scanData.size(); ++i) {
00102     if (!scanData[i].rawAvailable) return false;
00103     }
00104     return true;
00105 }
00106 
00107 bool AnalysisAlgorithm::hasAllFits() const {
00108     for (unsigned int i=0; i<scanData.size(); ++i) {
00109     if (!scanData[i].fitAvailable) return false;
00110     }
00111     return true;
00112 }
00113 
00114 void AnalysisAlgorithm::loadAllRaws() {
00115     for (unsigned int i=0; i<scanData.size(); ++i) {
00116     if (scanData[i].rawAvailable && !scanData[i].raw) 
00117         scanData[i].raw = dynamic_pointer_cast<RawScanResult>(IOManagerIS::instance().read(scanData[i].rawName, 0));
00118     }
00119 }
00120 
00121 
00122 void AnalysisAlgorithm::loadAllFits() {
00123     for (unsigned int i=0; i<scanData.size(); ++i) {
00124     if (scanData[i].fitAvailable && !scanData[i].fit) 
00125         scanData[i].fit = dynamic_pointer_cast<FitScanResult>(IOManagerIS::instance().read(scanData[i].fitName, 0));
00126     }
00127 }
00128 
00129 shared_ptr<RawScanResult> AnalysisAlgorithm::getRaw(unsigned int index) const {
00130     if (!scanData[index].rawAvailable || !scanData[index].raw) throw IllegalStateError("AnalysisAlgorithm::getRaw.  RawScanResult not available", __FILE__, __LINE__);
00131     return scanData[index].raw;
00132 }
00133 
00134 shared_ptr<FitScanResult> AnalysisAlgorithm::getFit(unsigned int index) const {
00135     if (!scanData[index].fitAvailable || !scanData[index].fit) throw IllegalStateError("AnalysisAlgorithm::getFit.  FitScanResult not available", __FILE__, __LINE__);
00136     return scanData[index].fit;
00137 }
00138 
00139 void AnalysisAlgorithm::mergeDefects() {
00140     if (!m_result) return;
00141     for (unsigned int i=0; i<scanData.size(); ++i) {
00142     if (scanData[i].fit) {
00143         m_result->getDefects() += scanData[i].fit->getDefects();
00144     }
00145     }
00146 }
00147 
00148 void AnalysisAlgorithm::finish() {
00149     checkTestResultExists();
00150     m_done=true;
00151     m_result->getHeader().setEndTime(second_clock::universal_time());
00152 
00153     m_result->getPassed() ? addPass() : addFailure();
00154 
00155     std::string options("SCAN_OPTIONS:");
00156     options+=getTestData()->options;
00157     m_result->addComment(options);
00158 
00159     //cout <<"Algorithm done"<<endl;
00160     try {
00161         Sct::IS::IOParamsIS param(AnalysisService::instance().getOutputISServerName());
00162         Sct::ISProxy::IOManagerISProxy::instance().write(*m_result, &param);
00163         //cout <<"TestResult published"<<endl;
00164     } catch (Throwable& e) {
00165         e.sendToMrs(MRS_ERROR);
00166     }
00167 
00168     try{
00169         TestSummaryIS t;
00170         t.dataString = SctData::TestSummary::SummaryManager::instance().write(*m_result);
00171         string name = AnalysisService::instance().getOutputISServerName();
00172         name += ".";
00173         name += "Summary.";
00174         name += (string) m_result->getUniqueID();
00175     Sct::ISUtilities::addOrUpdateOrThrow(name, t, __FILE__, __LINE__, MRS_ERROR);
00176     } catch (Throwable& e) {
00177         e.sendToMrs(MRS_ERROR);
00178     }
00179 
00180     try{
00181 
00182       //  This is a temporary workaround... just for testing purposes. BD 
00183       if ((m_result->getClassName()=="SctData::NPtGainTestResult") || (m_result->getClassName()=="SctData::NoiseOccupancyTestResult")) {
00184     SctData::ISTestSummary::ISSummaryManager::instance().write(*m_result);
00185       }
00186     } catch (Throwable& e) {
00187       e.sendToMrs(MRS_ERROR);
00188     }
00189 
00190 }
00191  
00192  AnalysisAlgorithm::ScanData::ScanData() : rawAvailable(false), fitAvailable(false){
00193  }
00194     
00195 string AnalysisAlgorithm::getStatus() const {
00196   boost::recursive_mutex::scoped_lock lock(m_status_access);
00197     ostringstream oss;
00198     oss << "Status: \n" << nAnalyses << " analyses done in " << analysisTime << "s\n" 
00199     << "I/O: " << ioTime << "s.  Overhead time: " << overheadTime << "s and called " << nCalls << " times\n"
00200     << "Total time: " << (analysisTime+ioTime+overheadTime) << "s average: " << ((analysisTime+ioTime+overheadTime)/nAnalyses) << endl;
00201     oss << "Passed : " << prototype.nPassed << ", Failed : " << prototype.nFailed << endl;
00202     return oss.str();
00203 }
00204 }

Generated on Thu Dec 15 21:14:22 2005 for SCT DAQ/DCS Software - C++ by doxygen 1.3.5