Changeset 11753 for trunk/FACT++/gui/RawEventsViewer/RawEventsViewer.cc

Timestamp:

08/03/11 14:21:39 (13 years ago)

Author:

lyard

Message:

added display of calibrated data

File:

• trunk/FACT++/gui/RawEventsViewer/RawEventsViewer.cc (modified) (12 diffs)

trunk/FACT++/gui/RawEventsViewer/RawEventsViewer.cc

@@ -514 +514 @@
     float yRange = bboxMax[1] - bboxMin[1];
     glColor3f(1.0,1.0,1.0);
+    float divideMe = (float)(VALUES_SPAN-1);
+    float plusMe = VALUES_SPAN/2;
+    if (drawCalibrationLoaded && calibrationLoaded)
+    {
+        divideMe /=2;
+        plusMe /=2;
+    }
     for (int i=0;i<nRoi-1;i++)
     {
 #ifdef LOAD_RAW
         glVertex2f(bboxMin[0] + xRange*i/(float)nRoi,
-                   bboxMin[1] + yRange*(eventsData[eventNum][selectedPixel][i]+(VALUES_SPAN/2)) /(float)(VALUES_SPAN-1));
+                   bboxMin[1] + yRange*(eventsData[eventNum][selectedPixel][i]+plusMe) /divideMe);
         glVertex2f(bboxMin[0] + xRange*(i+1)/(float)nRoi,
-                   bboxMin[1] + yRange*(eventsData[eventNum][selectedPixel][i+1]+(VALUES_SPAN/2)) /(float)(VALUES_SPAN-1));
+                   bboxMin[1] + yRange*(eventsData[eventNum][selectedPixel][i+1]+plusMe) /divideMe);
 #else
+
         glVertex2f(bboxMin[0] + xRange*i/(float)nRoi,
-                   bboxMin[1] + yRange*(eventData[nRoi*selectedPixel + i]+(VALUES_SPAN/2)) /(float)(VALUES_SPAN-1));
+                   bboxMin[1] + yRange*(eventData[nRoi*selectedPixel + i]+plusMe) /divideMe);
         glVertex2f(bboxMin[0] + xRange*(i+1)/(float)nRoi,
-                   bboxMin[1] + yRange*(eventData[nRoi*selectedPixel + i+1]+(VALUES_SPAN/2)) /(float)(VALUES_SPAN-1));
+                   bboxMin[1] + yRange*(eventData[nRoi*selectedPixel + i+1]+plusMe) /divideMe);
 #endif
     }
@@ -659 +667 @@
     for (int i=0;i<1440;i++)
         updateNeighbors(i);
+
+    calibrationLoaded = false;
+    drawCalibrationLoaded = false;
+
 }
 /************************************************************
@@ -936 +948 @@
     if (!inputFile->SetPtrAddress("Data", eventData)){
         cout << "Missing column " << "Data" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("EventNum", &eventNum)){
         cout << "Missing column " << "EventNum" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("TriggerType", &triggerType)){
         cout << "Missing column " << "TriggerType" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("SoftTrig", &softTrig)){
         cout << "Missing column " << "SoftTrig" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("PCTime", &pcTime)){
         cout << "Missing column " << "PCTime" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("BoardTime", boardTime)){
         cout << "Missing column " << "BoardTime" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("StartCellData", startPix)){
         cout << "Missing column " << "StartCellData" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     if (!inputFile->SetPtrAddress("StartCellTimeMarker", startTM)){
         cout << "Missing column " << "StartCellTimeMarker" << " Aborting load..." << endl;
-        nRoi = nRows = 0;}
+        nRoi = nRows = 0;return;}
     int backupStep = eventStep;
     rowNum = -1;
@@ -966 +978 @@
     emit signalCurrentPixel(selectedPixel);
 }
-
+void RawDataViewer::openCalibFile(string& file)
+{
+    calibrationLoaded = false;
+    calibInputFile = new fits(file);
+    if (!*calibInputFile)
+    {
+        delete calibInputFile;
+        calibInputFile = NULL;
+        return;
+    }
+
+    if (!calibInputFile->SetPtrAddress("BaselineMean", baseLineMean)){
+        cout << "Missing column " << "BaseLineMean" << " Aborting load..." << endl;
+        return;}
+    if (!calibInputFile->SetPtrAddress("GainMean", gainMean)){
+        cout << "Missing column " << "GainMean" << " Aborting load..." << endl;
+        return;}
+    if (!calibInputFile->SetPtrAddress("TriggerOffsetMean", triggerOffsetMean)){
+        cout << "Missing column " << "TriggerOffsetMean" << " Aborting load..." << endl;
+        return;}
+
+    calibInputFile->GetNextRow();
+
+    delete calibInputFile;
+
+    calibrationLoaded = true;
+    emit newFileLoaded();
+    if (drawCalibrationLoaded)
+        updateGL();
+}
 /************************************************************
  * PLUS EVENT
@@ -1009 +1050 @@
 #endif
 
+    if (drawCalibrationLoaded)
+    {
+        for (int i=0;i<1440;i++)
+            for (int j=0;j<nRoi;j++)
+            {
+                eventData[i*1024+j] /= 2;
+                eventData[i*1024+j] -= (baseLineMean[i*1024+j]+triggerOffsetMean[i*1024+j]);
+                eventData[i*1024+j] /= gainMean[i*1024+j];
+            }
+    }
     updateGL();
     emit signalCurrentEvent(eventNum);
     emit signalCurrentPixel(selectedPixel);
-
 }
 /************************************************************
@@ -1065 +1115 @@
     patchID = 1;
     rescaleWholeCamera = true;
+    currentFile = "none";
+    currentCalibFile = "none";
 }
 void UIConnector::slicesPlusPlus()
@@ -1073 +1125 @@
 {
     viewer->previousSlice();
+}
+void UIConnector::drawCalibratedDataChanged(int state)
+{
+     if (state)
+    {
+        if (viewer->calibrationLoaded)
+        {
+            viewer->drawCalibrationLoaded = true;
+            for (int i=0;i<1440;i++)
+                for (int j=0;j<viewer->nRoi;j++)
+                {
+                    viewer->eventData[i*1024+j] /= 2;
+                    viewer->eventData[i*1024+j] -= (viewer->baseLineMean[i*1024+j]+viewer->triggerOffsetMean[i*1024+j]);
+                    viewer->eventData[i*1024+j] /= viewer->gainMean[i*1024+j];
+                }
+            viewer->updateGL();
+        }
+        else
+        {
+            drawCalibrationCheckBox->setChecked(false);
+        }
+    }
+    else
+    {
+        viewer->drawCalibrationLoaded = false;
+        if (viewer->calibrationLoaded)
+        {
+            for (int i=0;i<1440;i++)
+                for (int j=0;j<viewer->nRoi;j++)
+                {
+                    viewer->eventData[i*1024+j] *= viewer->gainMean[i*1024+j];
+                    viewer->eventData[i*1024+j] += (viewer->baseLineMean[i*1024+j]+viewer->triggerOffsetMean[i*1024+j]);
+                    viewer->eventData[i*1024+j] *= 2;
+               }
+            viewer->updateGL();
+        }
+
+    }
+
 }
 /************************************************************
@@ -1306 +1397 @@
     dialog.exec();
 }
+void UIConnector::loadNewCalibFileClicked()
+{
+    QFileDialog dialog;
+    dialog.setFileMode(QFileDialog::ExistingFile);
+    dialog.open(this, SLOT(calibFileSelected(QString)));
+    dialog.setVisible(true);
+    dialog.exec();
+}
 /************************************************************
  * FILE SELECTED. return of the file open dialog handler
@@ -1315 +1414 @@
         viewer->openFile(currentFile);
 }
+void UIConnector::calibFileSelected(QString file)
+{
+    currentCalibFile = file.toStdString();
+    if (currentCalibFile != "")
+        viewer->openCalibFile(currentCalibFile);
+}
 /************************************************************
  * NEW FILE LOADED. update of the UI after a new file has been loaded
@@ -1322 +1427 @@
     ostringstream str;
     str << "File loaded: " << currentFile << "\n";
+    str << "Calibration file loaded: " << currentCalibFile << "\n";
 //    fileLoadedLabel->setText(QString(str.str().c_str()));
 //    str.str("");
@@ -2123 +2229 @@
     connector.range3 = myUi.colorRange3;
     connector.range4 = myUi.colorRange4;
+    connector.drawCalibrationCheckBox = myUi.calibratedCheckBox;
 
     connector.initHistograms();
@@ -2162 +2269 @@
                      &connector, SLOT(newFileLoaded()));
 
+    QObject::connect(myUi.calibratedCheckBox, SIGNAL(stateChanged(int)),
+                     &connector, SLOT(drawCalibratedDataChanged(int)));
     QObject::connect(myUi.loadNewFileButton, SIGNAL(clicked()),
                      &connector, SLOT(loadNewFileClicked()));
+    QObject::connect(myUi.loadDRSCalibButton, SIGNAL(clicked()),
+                     &connector, SLOT(loadNewCalibFileClicked()));
 
     QObject::connect(myUi.colorRange0, SIGNAL(valueChanged(double)),