Changeset 14791 for trunk/FACT++/gui

Timestamp:

01/28/13 10:21:39 (12 years ago)

Author:

lyard

Message:

Added Monte-Carlo data support. Disable 3D viewer.

Location:

trunk/FACT++/gui

Files:

• Q3DCameraWidget.cc (modified) (3 diffs)
• RawEventsViewer/RawEventsViewer.cc (modified) (14 diffs)
• RawEventsViewer/RawEventsViewer.h (modified) (2 diffs)

trunk/FACT++/gui/Q3DCameraWidget.cc

@@ -163 +163 @@
     void Q3DCameraWidget::setData(float* ddata)
     {
-        for (int i=0;i<1440;i++)
+        cout << "Info : 3D plotter disabled. requires more work so that less than 300 slices per pixel can be loaded" << endl;
+        cout << "Contact Etienne (etienne.lyard@unige.ch) for more information." << endl;
+        //disabled for now as I am working with 150 slices only
+/*        for (int i=0;i<1440;i++)
             for (int j=0;j<300;j++)
             _data[i*300+j] = (short)(ddata[i*300 + j]);
@@ -169 +172 @@
         if (isVisible())
             updateGL();
-
+*/
     }
     void Q3DCameraWidget::setData(short* ddata)
     {
-        for (int i=0;i<1440;i++)
+        cout << "Info : 3D plotter disabled. requires more work so that less than 300 slices per pixel can be loaded" << endl;
+        cout << "Contact Etienne (etienne.lyard@unige.ch) for more information." << endl;
+/*        for (int i=0;i<1440;i++)
             for (int j=0;j<300;j++)
                 _data[i*300+j] = ddata[i* 300 + j];
@@ -179 +184 @@
         if (isVisible())
             updateGL();
-
+*/
     }
     void Q3DCameraWidget::drawCameraBody()

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

@@ -189 +189 @@
         else
 
-        color = float(eventData[nRoi*hardwareMapping[i] + whichSlice]+(VALUES_SPAN/2))/(float)(VALUES_SPAN-1);
+//        if (_softwareOrdering)
+//            color = float(eventData[nRoi*i + whichSlice] + (VALUES_SPAN/2))/(float)(VALUES_SPAN-1);
+//        else
+            color = float(eventData[nRoi*hardwareMapping[i] + whichSlice]+(VALUES_SPAN/2))/(float)(VALUES_SPAN-1);
         if (logScale)
         {
@@ -291 +294 @@
         }*/
 
-    const int hw = hardwareMapping[selectedPixel];
+//    int mapping = _softwareOrdering ? selectedPixel : hardwareMapping[selectedPixel];
+    int mapping = hardwareMapping[selectedPixel];
+    const int hw = mapping;
     const PixelMapEntry& mapEntry = fPixelMap.index(selectedPixel);
     const int pixelIdInPatch = mapEntry.pixel();
@@ -385 +390 @@
             patchesColor[i][j] = color[j];
     }
+    fZeroArray = NULL;
+
+    _softwareOrdering = false;
 
 }
@@ -402 +410 @@
         delete[] waveLetArray;
     }
-}
-
+    if (fZeroArray != NULL)
+        delete[] fZeroArray;
+}
+void RawDataViewer::allocateZeroArray()
+{
+    if (fZeroArray == NULL)
+    {
+        fZeroArray = new char[8192];
+    }
+}
 /************************************************************
  * PAINT GL. main drawing function.
@@ -572 +588 @@
         }
     }
+
     nRows =          inputFile->HasKey("NAXIS2") ? inputFile->GetInt("NAXIS2") : 0;
     nRoi =           inputFile->HasKey("NROI") ?  inputFile->GetInt("NROI") : 0;
@@ -584 +601 @@
     nBoards =        inputFile->HasKey("NBOARD") ? inputFile->GetInt("NBOARD") : 0;
     nPixels =        inputFile->HasKey("NPIX") ?  inputFile->GetInt("NPIX") : 0;
-    timeSystem =     inputFile->HasKey("TIMESYS") ? inputFile->GetStr("TIMESYS") : 0;
-    creationDate =   inputFile->HasKey("DATE") ? inputFile->GetStr("DATE") : 0;
+    timeSystem =     inputFile->HasKey("TIMESYS") ? inputFile->GetStr("TIMESYS") : "";
+    creationDate =   inputFile->HasKey("DATE") ? inputFile->GetStr("DATE") : "";
     nightInt =       inputFile->HasKey("NIGHT") ? inputFile->GetInt("NIGHT") : 0;
-    camera =         inputFile->HasKey("CAMERA") ? inputFile->GetStr("CAMERA") : 0;
-    daq =            inputFile->HasKey("DAQ") ? inputFile->GetStr("DAQ") : 0;
+    camera =         inputFile->HasKey("CAMERA") ? inputFile->GetStr("CAMERA") : "";
+    daq =            inputFile->HasKey("DAQ") ? inputFile->GetStr("DAQ") : "";
     adcCount =       inputFile->HasKey("ADCRANGE") ? inputFile->GetFloat("ADCRANGE") : 2000;
+    if (nPixels == 0)
+    {
+        cout << "could not read num pixels from fits header. Assuming 1440 (FACT)." << endl;
+        nPixels = 1440;
+    }
+    if (nRoi == 0 && !inputFile->HasKey("NROI"))
+    {//let's try to figure out the roi from the column's format
+        const fits::Table::Columns& cols = inputFile->GetColumns();
+        if (cols.find("Data") == cols.end())
+        {
+            cout << "ERROR: Column \"Data\" could not be found. abortin load." << endl;
+            return;
+        }
+        const fits::Table::Columns::const_iterator col = cols.find("Data");
+        if (col->second.type != 'I')
+        {
+            cout << "ERROR: Data Column has type " << col->second.type << " while viewer expects I" << endl;
+            return;
+        }
+        if (col->second.num % nPixels != 0)
+        {
+            cout << "ERROR: Num pixels (" << nPixels << ") does not match Data length (" << col->second.num << "). Aborting" << endl;
+            return;
+        }
+        nRoi = col->second.num/nPixels;
+        cout << "Estimate num samples per pixels to be " << nRoi;
+        _softwareOrdering = true;
+    }
+    else
+        _softwareOrdering = false;
 
     if (inputFile->HasKey("OFFSET"))
@@ -598 +645 @@
     nbBad = 0;//inputFile->GetInt("NBEVTBAD");
 
-    eventNum = 0;
+    eventNum = 1;
 
     if (eventData != NULL) {
@@ -612 +659 @@
     {
         inputFile->SetPtrAddress("Data", rawEventData);
-        inputFile->SetPtrAddress("EventNum", &eventNum);
-        inputFile->SetPtrAddress("TriggerType", &triggerType);
-        inputFile->SetPtrAddress("SoftTrig", &softTrig);
-        inputFile->SetPtrAddress("BoardTime", boardTime);
-        inputFile->SetPtrAddress("StartCellData", startPix);
-        inputFile->SetPtrAddress("StartCellTimeMarker", startTM);
-        inputFile->SetPtrAddress("TimeMarker", &rawEventData[1440*nRoi]);
+        if (inputFile->HasColumn("EventNum"))
+            inputFile->SetPtrAddress("EventNum", &eventNum);
+        else
+            cout << "Warning: could not find column \"EventNum\"" << endl;
+        if (inputFile->HasColumn("TriggerType"))
+            inputFile->SetPtrAddress("TriggerType", &triggerType);
+        else
+            cout << "Warning: could not find column \"TriggerType\"" << endl;
+        if (inputFile->HasColumn("SoftTrig"))
+            inputFile->SetPtrAddress("SoftTrig", &softTrig);
+        else
+            cout << "Warning: could not find column \"SoftTrig\"" << endl;
+        if (inputFile->HasColumn("BoardTime"))
+            inputFile->SetPtrAddress("BoardTime", boardTime);
+        else
+            cout << "Warning: could not find column \"BoardTime\"" << endl;
+        if (inputFile->HasColumn("StartCellData"))
+            inputFile->SetPtrAddress("StartCellData", startPix);
+        else
+            cout << "Warning: could not find column \"StartCellData\"" << endl;
+        if (inputFile->HasColumn("StartCellTimeMarker"))
+            inputFile->SetPtrAddress("StartCellTimeMarker", startTM);
+        else
+            cout << "Warning: could not find column \"StartCellTimeMarker\"" << endl;
+        if (inputFile->HasColumn("TimeMarker"))
+            inputFile->SetPtrAddress("TimeMarker", &rawEventData[1440*nRoi]);
+        else
+            cout << "Warning: could not find column \"TimeMarker\"" << endl;
     }
     catch (const runtime_error &e)
@@ -633 +701 @@
     {
         pcTime[0] = pcTime[1] = 0;
-        inputFile->SetPtrAddress("UnixTimeUTC", pcTime);
+        if (inputFile->HasColumn("UnixTimeUTC"))
+            inputFile->SetPtrAddress("UnixTimeUTC", pcTime);
     }
     catch (const runtime_error&)
@@ -639 +708 @@
             try
         {
-            inputFile->SetPtrAddress("PCTime", pcTime);
+            if (inputFile->HasColumn("PCTime"))
+                inputFile->SetPtrAddress("PCTime", pcTime);
+            else
+                cout << "Warning: could not find column \"UnixTimeUTC\" nor \"PCTime\"" << endl;
+
         }
         catch (const runtime_error&)
@@ -824 +897 @@
         return;
     inputFile->GetRow(rowNum);
+    if (_softwareOrdering)
+    {//remap pixels data according to hardware id
+        if (nRoiTM != 0)
+            cout << "Warning: did not expect Time Markers data from Monte-Carlo simulations. These will not be mapped properly." << endl;
+        //first copy the data
+        int16_t* tempData = new int16_t[1440*nRoi];
+        for (int i=0;i<1440*nRoi;i++)
+            tempData[i] = rawEventData[i];
+        //copy back the data and re-map it on the fly
+        for (int i=0;i<1440;i++)
+            for (int j=0;j<nRoi;j++)
+                rawEventData[i*nRoi + j] = tempData[softwareMapping[i]*nRoi + j];
+
+        delete[] tempData;
+    }
 //    cout << "Getting row " << rowNum << endl;
 
@@ -1273 +1361 @@
         return;
     }
-    const int idx =
-        GLWindow->nRoi*GLWindow->hardwareMapping[GLWindow->selectedPixel]
-        + GLWindow->whichSlice;
+//    int mapping = GLWindow->_softwareOrdering ? GLWindow->selectedPixel : GLWindow->hardwareMapping[GLWindow->selectedPixel];
+    int mapping = GLWindow->hardwareMapping[GLWindow->selectedPixel];
+    const int idx = GLWindow->nRoi*mapping + GLWindow->whichSlice;
 
     ostringstream str;
@@ -1874 +1962 @@
         return;
     int softwarePix = pixel;
-    pixel = GLWindow->hardwareMapping[pixel];
+//    if (!GLWindow->_softwareOrdering)
+        pixel = GLWindow->hardwareMapping[pixel];
 
     HwIDBox->setValue(pixel);
@@ -1997 +2086 @@
         return;
 
-    HwIDBox->setValue(GLWindow->hardwareMapping[swid]);
+//    if (GLWindow->_softwareOrdering)
+//        HwIDBox->setValue(swid);
+//    else
+        HwIDBox->setValue(GLWindow->hardwareMapping[swid]);
 }
 

trunk/FACT++/gui/RawEventsViewer/RawEventsViewer.h

@@ -120 +120 @@
     valarray<double> Maxvalues;
     valarray<double> PosOfMaxvalues;
+
+    ///Used to load zero data in case of missing fits columns
+    void allocateZeroArray();
+    char* fZeroArray;
+    //bool telling whether the data is natively ordered in software or hardware id
+    //used to correctly display monte-carlo data.
+    bool _softwareOrdering;
 
 private:
@@ -177 +184 @@
 
 
+
 //    int hardwareMapping[1440];
 //    int softwareMapping[1440];